Spin partners of the $Z_b(10610)$ and $Z_b(10650)$ revisited

We study the implications of the heavy-quark spin symmetry for the possible spin partners of the exotic states $Z_b(10610)$ and $Z_b(10650)$ in the spectrum of bottomonium. We formulate and solve numerically the coupled-channel equations for the $Z_b$ states that allow for a dynamical generation of these states as hadronic molecules. The force includes short-range contact terms and the one-pion exchange potential, both treated fully nonperturbatively. The strength of the potential at leading order is fixed completely by the pole positions of the $Z_b$ states such that the mass and the most prominent contributions to the width of the isovector heavy-quark spin partner states $W_{bJ}$ with the quantum numbers $J^{++}$ ($J=0,1,2$) come out as predictions. Since the accuracy of the present experimental data does not allow one to fix the pole positions of the $Z_b$'s reliably enough, we also study the pole trajectories of their spin partner states as functions of the $Z_b$ binding energies. It is shown that, once the heavy-quark spin symmetry is broken by means of the physical $B$ and $B^*$ masses, especially the pion tensor force has a significant impact on the location of the partner states clearly demonstrating the need of a coupled-channel treatment of pion dynamics to understand the spin multiplet pattern of hadronic molecules.


Introduction
The experimental discovery of the charmonium-like state X(3872) by the Belle Collaboration in 2003 [1] inaugurated a new era in hadron spectroscopy. A lot of new candidates for exotic states have been discovered since then in the spectrum of both charmonium and bottomonium -for a review see, for example, Refs. [2,3]. Among those one should especially mention the isovector Z b (10610) and Z b (10650) resonances (for brevity, hereinafter often referred to as Z b and Z b , respectively) [4][5][6]. Their signals are seen in 7 channels, namely Υ(10860) → πZ Since Z b and Z b are charged but decay into final states containing a b and an anti-b quark they must contain at least four quarks and are thus explicitly exotic. It turned out that the decays of the Z b 's into the open-flavour channels almost exhaust their widths, despite the limited phase spaces. This feature is to be regarded as a strong evidence for a molecular nature of the Z b states [7] -for a recent review of the theory of hadronic molecules we refer to Ref. [8].
The corresponding masses and widths extracted from the Breit-Wigner fits for the data and averaged over the above production and decay channels are [9] M Z b = 10607.2 ± 2.0 MeV, Γ Z b = 18.4 ± 2.4 MeV, (1.2) M Z b = 10652.2 ± 1.5 MeV, Γ Z b = 11.5 ± 2.2 MeV, (1.3) to be compared with 10604 MeV and 10649 MeV for the BB * and the B * B * thresholds, respectively. Therefore, both poles are located within a couple of MeV from the respective threshold. The J P C quantum numbers of both states were determined as 1 +− [10] in line with the expectation of the molecular model. Within the molecular picture, the wave functions of the Z b states can be written as [7] 1 + (1 +− ) : where we quote the quantum numbers in the form I G (J P C ) and, for example, 1 − bb denotes the wave function of the bb pair with the total spin 1, and so forth. Such an identification of the Z b 's allows one to explain the dipion transition rates from the Υ(10860) to the vector bottomonia Υ(nS) (n = 1, 2, 3), which appear to be up to two orders of magnitude larger than similar transitions among the lower bottomonia Υ(nS) (n = 1, 2, 3) [11]. In addition, the molecular interpretation of the Z b 's (see Eqs. (1.4)-(1.5)) naturally explains that transitions such as Z ( ) b → πh b (mP ), which are expected to be suppressed by (Λ QCD /m b ) since they involve a change in the heavy-quark spin, happen at a comparable rate to the heavy-quark spin preserving transitions Z ( ) b → πΥ(nS) [4]. Molecular states should in general be located below the most relevant threshold and not above 1 , which seems to be in conflict with the numbers quoted in Eqs. (1.2) and (1.3). However, these numbers need to be interpreted with care. While they show unambiguously that the poles related to the Z b states reside very close to the corresponding open-flavour thresholds, their particular values contain an uncontrollable intrinsic systematic uncertainty since they were determined from sums of the Breit-Wigner functions that ignore the nearby thresholds and which, therefore, appear to be in conflict with analyticity and unitarity. Also, it is argued, for example, in Ref. [14] (see also Ref. [15] for a related discussion in the context of f 0 (980)/a 0 (980)) that the branching fractions extracted from the Breit-Wigner parameterisations do not represent the decay probabilities for near-threshold states. Thus, a more refined data analysis is required. For example, it is demonstrated in Ref. [16] that both Z b 's are compatible with bound state poles in the data for the h b (mP )π channels as soon as the energy dependence of their self-energies is included properly. On the other hand, a combined analysis of the experimental data in all seven channels listed in Eqs. (1.1) consistent with analyticity and unitarity favours both Z b 's as virtual states located within approximately 1 MeV below the respective thresholds [17,18] 2 . In other words, the quality of the existing data does not allow one to draw definite conclusions about the pole locations of these states. Thus, in the present paper, we investigate the fate of the symmetry partners of the Z b states as the binding energies of the latter are 1 Meson-meson dynamics can also lead to above-threshold poles if the interaction is energy-dependent, but such states are not expected to be very narrow [12]. On the other hand, a tetraquark nature of the states is claimed to result in the poles lying slightly above threshold [13].
2 While the formalism of these references is indeed unitary with respect to all seven channels, only the open-flavour and the h b (mP )π channels were included in the fit, since the Υ(nS)π channels call for an additional inclusion of nonresonant production. varied up to zero. As soon as the analysis of Refs. [17,18] is refined to include one-pion exchange interactions, the predictions for the partner states would have become possible directly from the fit for the experimental line shapes. However, a detailed investigation of the line shapes, which also calls for an inclusion of the inelastic channels, is delegated to future studies. In the limit of an infinite b-quark mass, a molecular nature of the Z b states allows one to predict the existence of spin partners with the quantum numbers J P C = J ++ , conventionally denoted in the literature as W bJ , with J = 0, 1, 2 [7,19,20] (see also Ref. [21] for a recent review and for the discussion of the future experimental tasks). In particular, in the notations of Eqs. (1.4) and (1.5), the wave functions of the spin partners take the form Although a detailed microscopic theory is needed to predict individual decay widths of the Z b and W bJ states, the spin symmetry allows one to arrive at various relations between their total widths, for example [19], and [20] (1.12) Additional predictions for the partial decay widths into various hidden-bottom final states can be also found in Refs. [19,20].
To arrive at the above predictions the mass splitting was treated as a large scale (compared to the binding energies) which was integrated out. Hereinafter m and m * denote the B and B * masses, respectively. While this assumption is acceptable for the Z b states themselves 3 , for their spin partner it turns out to be appropriate only within a truncated scheme when only the S-wave interactions are retained. However, it is the central finding of this paper that, as soon as the one-pion exchange (OPE) is included, the effect of the mass difference δ is enhanced significantly via the strong S-D transitions that come with it. This results in binding energies of the spin partner states as large as 20 MeV, which are clearly of the order of the spin symmetry violating parameter δ -see Eq. (1.13). While pion exchanges were already included in studies of the Z b states before [23], to the best of our knowledge its impact on the spin partners has not been investigated so far. Furthermore, since, in Ref. [24], the effect of the OPE was claimed to be somewhat diminished if the one-η exchange (OEE) is added, we also investigate the effect of this addition. It is well known that a nonperturbative inclusion of the OPE interaction in the chiral nuclear effective-field theory allows one to significantly extend the region of applicability of the theory -see, for example, Refs [25,26]. In charmonium-like systems such as the X(3872), pionic effects turn out to be important to predict the dynamical properties such as the light-quark mass dependence of the X-pole [27,28] (for a discussion of the pion mass dependence with perturbative pions we refer to Ref. [29]) and the decay width X(3872) → DDπ [30]. Furthermore, it was shown in Ref. [22] that the inclusion of the OPE interactions had a strong impact on the location of the spin partner states of the X(3872). In particular, the iterations of the OPE interaction to all orders were found to generate strong coupled-channel effects due to the then allowed transitions D * D * ↔ D ( * )D . As a consequence, the binding energy and the width of the spin-2 partner of the X(3872) both appeared to be of the order of several dozens of MeV. Both values were found to be significantly larger than those reported in the literature earlier [31].
In this paper, a calculation for the spin partners of the Z b states is presented in which the pole positions of the Z b 's are treated as input. The underlying coupled-channel model includes short-range contact interactions, OPE, and OEE fully iterated to all orders. As our starting point, we assume both Z b 's to be shallow bound states and investigate the role played by the spin symmetry violation introduced into the system via a nonzero value of the parameter δ defined in Eq. (1.13). In addition, we study the behaviour of the spin partners as the binding energies of the Z b states are reduced (up to zero values) which allows us to make statements about the masses of the spin partners when the Z b 's turn to virtual states. In particular, we find that the tensor spin partner W b2 with the quantum numbers 2 ++ survives in this limit as a bound state and it is expected to produce a visible resonant structure in the BB and BB * line shapes a few MeV below the B * B * threshold.
The key aim of this study is to highlight explicitly the nontrivial impact of the nonperturbative one-pion exchange on the location of the spin partners of the Z b states. For that it appears sufficient to base the study solely on the open-flavour channels introduced below. In particular, we neither include the inelastic channels (πΥ(nS) and πh b (mP )) nor do we aim at a high accuracy description of the data. Clearly, as soon as the analysis of Refs. [17,18] is refined to include one-pion exchange interactions, the predictions for the partner states could be refined further. Such a study, however, is beyond the scope of this paper.
The paper is organised as follows. In Sec. 2, we outline various implications of the heavy-quark spin symmetry (HQSS) in the purely contact theory. Then, in Sec. 3, we discuss the inclusion of the OPE interaction on top of the contact potential introduced in Sec. 2. The results obtained in the theory which incorporates both the contact and the OPE interactions are collected and discussed in Sec. 4. A peculiar alternative solution found in the contact theory which, however, disappears from the theory once the OPE interaction is taken into account is highlighted in Sec. 5. We give an overview of the results obtained in the concluding Sec. 6.

Contact theory
We start from the purely contact theory. The basis of states -see Ref.
[32] -can be adapted directly from the c-sector to the b-sector, Here, the individual partial waves are labelled as 2S+1 L J with S, L, and J denoting the total spin, the angular momentum, and the total momentum of the two-meson system, respectively. The C-parity eigenstates are defined as 4 and comply with the convention for the C-parity transformationĈM =M.
In the basis of Eq. (2.1) and for a given set of the J P C quantum numbers the leadingorder EFT potentials V (JP C) LO which respect the heavy-quark spin symmetry read [32][33][34] where {C 1a , C 1b } and {C 1 , C 1 } are two alternative sets of low-energy constants (contact terms).
With the help of the appropriate unitary transformations, the potential matrices (2.3) and (2.4) can be diagonalised to take the form diag(C 1 , C 1 ). Therefore, in the strict heavyquark limit (δ = 0) and in the purely contact theory, if the potential C ( ) 1 is strong enough to bind the system in one channel, it is inevitably strong enough to produce bound state(s) in the sibling channel(s). Furthermore, since the binding energies in different channels are governed by the same combinations of the contact terms, the corresponding molecular states appear to be degenerate in mass. This observation allowed the authors of Refs. [22,35] to predict, in the strict HQSS limit, the existence of multiple degenerate isosinglet states in the c-sector. Obviously, this conclusion trivially translates into the b-sector and holds for the isovector states, too. In particular, one finds where the superscript (0) indicates the exact HQSS limit which implies the relation between the spin symmetry violating parameter δ, the typical binding energy E B in the Z b and W bJ states, and the heavy-meson mass. Therefore, the input provided by the masses of the two Z b states is sufficient to fix the parameters of the contact interaction at leading order. In Refs. [7,19,20,36] a different relation between the relevant scales of the problem was assumed, while the contact interaction at this order was still treated as spin symmetric -see Eqs. (2.3)-(2.6). Then, the degeneracy (2.7) is lifted and instead one arrives at some relations between the binding momenta of the partners, which, in particular, yield the relations (1.11)-(1.12) for the widths. It should be noted that, to this order, the states Z b and W b0 cannot yet decay into the (open) channels BB * and BB, respectively, and the widths (1.11)-(1.12) originate entirely from inelastic transitions. As was shown in Ref. [22], already at order O(δ), relations (2.10) acquire corrections linear with the cutoff introduced to regularise the divergent integrals in the Lippmann-Schwinger equations. Then, to absorb this cutoff dependence, the renormalisation in principle calls for the presence of HQSS violating contact interactions, in addition to the HQSS preserving potentials (2.3)-(2.6). Since the HQSS violating counter terms are not known, as an alternative approach, one may vary the cutoff in some reasonable range to estimate the uncertainty in the position of the partner states due to the omission of these interactions. In practice, the cutoff dependence appears to be very marginal (at least in the pionless theory) which can be, in part, justified by the large b-quark mass (much larger than the characteristic cutoff varied in the range from the pion mass and up to the chiral symmetry breaking scale ∼ 1 GeV) and by the additional suppression of such terms by a factor γ/m (at least for the uncoupled partial waves) [22]. Note that the transitions B * B * → BB * and B * B * → BB do not vanish already in the pionless theory, such that the coupled-channel dynamics for the states Z b and W b0 produce an additional shift of the binding momenta and of the imaginary part of the order O(γ 2 / √ mδ) [22]. Moreover, as it will be shown in this paper, once the OPE interaction is included, the state W b2 can also participate in the coupled-channel transitions, the latter having a significant impact on the location of this state.

One-pion exchange interaction
We start from the lowest-order nonrelativistic interaction Lagrangian [37,38], and employ the heavy-quark flavour symmetry to equate the dimensionless coupling constant g b to the similar constant g c which parametrises the D * Dπ vertex. It can be determined directly from the observable D * → Dπ decay width via where q is the centre-of-mass momentum in the final state. The numerical value extracted from data [9] reads and it agrees within 10% with the recent lattice QCD determination of the B * Bπ coupling constant [39].
Depending on the channel, the OPE potential contains the pion propagator and two vertices of the type B * → B ( * ) π following directly from Lagrangian (3.1), , and * stand for the polarisation vectors of the initial and final B * mesons, q is the pion momentum, τ a is the isospin Pauli matrix, and f π = 92.2 MeV is the pion decay constant -see also Refs. [22,40]. Then, the OPE potentials connecting the heavy-meson B ( * ) B ( * ) pairs in the initial and in the final state read 5 where p (p ) denotes the centre-of-mass momentum of the initial (final) heavy-meson pair and the pion momentum is q = p + p . For the isovector states, the isospin factor which appears from the operator τ 1 · τ c 2 is Note that we use a convention where the integration weight is absorbed into the potential. Accordingly, the usual factor (2π) 3 does not appear in the integral equation (3.12) given below.
where τ c stands for the generator of the antifundamental representation, τ c = τ 2 τ T τ 2 = −τ . Note that in contrast to what is given above, in the isoscalar channel, studied, for example, in Ref. [22], the isospin factor is +3. Furthermore, the expressions for the propagators in Eq. (3.5) read where E π = q 2 + m 2 π is the pion energy and √ s defines the total energy of the system which is conveniently represented as √ s = m + m * + E, with E counted with respect to the BB * threshold. The time-reversed transition potentials are trivially obtained from Eqs. (3.5) by interchanging p and p and using that D B * Bπ (p, p ) = D BB * π (p , p).
The three-body effects are incorporated in the above OPE potentials via the heavymeson recoil corrections and via the energy-dependent terms in the time-ordered propagators from Eqs. (3.7)-(3.9). The leading-order static OPE potentials are then easily obtained from the full results by neglecting the three-body terms, that is by the replacement (3.10) We checked numerically that, for the mass spectra which are the focus of this work, the three-body terms are basically irrelevant. Accordingly, one could as well have used the static pion propagators. The reason is that the B * -B mass difference δ -see Eq. (1.13) -is about three times smaller than the pion mass and, therefore, there are no three-body cuts located in the vicinity of B ( * ) B ( * ) thresholds. This situation is, however, very different in the c-sector: the DDπ threshold lies about 7 MeV below the DD * threshold and, hence, three-body effects play a prominent role for understanding the dynamics of the X(3872) state treated as a DD * molecule [27,28,30]. The OEE can be calculated straightforwardly from the expressions given above if one makes the following replacements: (i) the isospin factor of Eq. (3.6) by +1; (ii) the pion mass by the η mass, and (iii) the pion coupling constant g c by the η coupling constant g c / √ 3 which can be found, for example, in Ref. [41]. As soon as the OPE potential is included, it enables transitions to D and even G waves which, therefore, have to be included in the extended set of basis states, As before, the C-parity of the states is indicated explicitly in parentheses whenever necessary.
The coupled-channel integral equations for the scattering amplitude take the form (3.12) where i, j, and n label the basis vectors, as they appear in Eq. (3.11), and are the propagator and the reduced mass of the heavy meson-antimeson pair in the given channel evaluated in its centre of mass. The potential V (JP C) ij (p, p ) includes the contact term as well as the meson exchange interactions in the channel with the given quantum numbers J P C and in the given partial wave. The Lippmann-Schwinger-type equations derived are regularised with the sharp cutoff, where the value of the latter is chosen to be of the order of a natural hard scale in the problem -see Ref. [22]. All necessary details of the coupled-channel formalism, of the partial wave projection, and so on can be found in Refs. [22,42]. However, it is important to pinpoint a crucial difference between the potentials used in Refs. [22] and the ones from this work: while the OPE potential from Ref. [22] describes the interaction in the isosinglet channels, here it operates in the isovector systems. Thus, the trivial modification one needs to make in the equations is to change the isospin coefficient 3, which multiplies the OPE potential in the isosinglet channels, for the isospin coefficient -1, which corresponds to the isovector channels -see Eq. (3.6). In addition, it has to be noticed that the sign of the OPE potential also depends on the Cparity of the system. Then, as a net result, the central (S-wave) part of the OPE potential is attractive in the C-odd Z b 's channel while it is repulsive in the C-even channels where the W bJ states are predicted. It is important to note, however, that the net effect of the OPE can be attractive also in those channels due to the nondiagonal transition potentials -in particular, due to the S-D transitions induced by the tensor force.
The inclusion of the OPE interaction modifies the potentials (2.3)-(2.6). Nevertheless, as demonstrated in Ref. [22], in the strict heavy-quark limit, if all coupled channels and all relevant partial waves are taken into account, for a given set of quantum numbers J P C , one can find the unitary operator that block diagonalises the potential into the form (for the sake of transparency, the size of the blocks is quoted explicitly in parentheses), This conclusion does not depend on the isospin of the system and it is, therefore, equally valid for both isosinglets and isovectors. Since the contact interactions contribute to both matrix A (via C 1 = C 1a + C 1b ) and matrix B (via C 1 = C 1a − 3C 1b ) and since at least one of these two matrices enters all four potentials for all quantum numbers one expects that the existence of the Z b and Z b as bound states in the 1 +− channel entails the existence of the bound states in the other three channels, too, with the degeneracy pattern given by Eq. (2.7) as in case of the purely contact theory. As explained in the previous section, this pattern is expected to be lifted as soon as the physical masses of the B and B * mesons are used in the calculations. Meanwhile, the HQSS constraints are still expected to work quite well for the contact potentials.

Results and discussion
In this section, we present and discuss the results of our calculations with special emphasis on the role of the pion dynamics for the location of the poles of the W bJ 's -the spin partners of the Z b states. Furthermore, motivated by the fact that in the isovector channel the relative importance of the η-exchange is larger than in the isoscalar channel, we also include it explicitly in our calculations.
As our starting point, we assume both Z b and Z b to be shallow bound states and treat their binding energies as input parameters. This allows us to fix completely the HQSSconstrained leading-order contact potential -see Eqs. (2.3)-(2.6). For definiteness, as an example, we fix the binding energies of the Z b 's to be in line with Ref. [16]. For both Z b 's the energy is counted relative to the respective reference threshold explicitly stated in Eq. (4.1) in parentheses. Ideally, one would need to extract pole locations and residues for the spin partner states of the Z b 's directly from the calculated T -matrices. However, this requires an analytic continuation of the amplitudes into the complex plain that goes beyond the scope of the present paper. Therefore, to extract the binding energies and the widths of the spin partners, W bJ 's, we mimic the experimental procedure and calculate the production rates in the elastic channels for the spin partner states. In what follows, similarly to Eq. (4.1), the binding energies of the partner states will be always defined relative to their reference thresholds, namely, Threshold: Spin partner: (4.2)

Dependence on the HQSS breaking scale
To highlight the impact of the OPE and OEE on the pole locations of the spin partner states, in Fig. 1, we show the evolution of the extracted binding energies as functions of δ -the HQSS breaking mass splitting between the B and B * mesons -for four different scenarios: • Scenario A: purely contact potential; • Scenario B: contact potential plus the central (S-wave) OPE interaction; • Scenario C: contact potential plus the full OPE interaction; • Scenario D: contact potential plus the full OPE and OEE interactions altogether.
It follows from Eqs. (2.7) and (3.14) that, in the strict HQSS limit (δ = 0), the Z b 's and the W bJ 's populate two families of states and that this conclusion holds both with and without OPE and OEE interactions included. Given the input from Eq. (4.1), this explains why, for all four scenarios above, the 2 ++ , 1 ++ , and one of the 0 ++ states have the same binding energy E B = 5 MeV while the other 0 ++ spin partner has the binding energy E B = 1 MeV in this limit. However, for δ > 0, HQSS is lifted and the exact degeneracy (2.7) is lost -see the discussion in Secs. 2 and 3 above. Then, as δ grows from zero to the physical value of 45 MeV, the spin partners from the first family (W b0 , W b1 , and W b2 ) tend to become more bound while the state W b0 becomes unbound fast and its pole moves to the second Riemann sheet. 6 Note that while varying δ we keep adjusting the parameters of the contact potential such that the pole locations of Z b and Z b stay fixed at the values given in Eqs. (4.1).
From Fig. 1 one can conclude that the central (S-wave) part of the OPE has no influence on the position of the poles 7 while its tensor part has a significant impact, especially on the 2 ++ state W b2 . Indeed, due to the additional attraction which stems from the OPE-related tensor forces, the binding energy of the W b2 grows by almost a factor of two as compared to the results of the pionless theory (from 13 to 23 MeV). Furthermore, by comparing the results in the strict HQSS limit (δ = 0) with those for the physical mass splitting (δ = 45 MeV), one can conclude that spin symmetry is violated quite strongly and that a substantial amount of this violation for the 2 ++ state stems from the tensor forces.
Finally, it can also be seen in Fig. 1 that the η-exchange does not play a prominent role for the systems under study, for it only slightly diminishes the OPE. This can be seen in all four plots in Fig. 1.

Dependence on the one-boson-exchange strength parameter
To further illustrate the role played by the dynamics governed by the pion and η-meson exchanges, in Fig. 2, we plot the variation of the binding energies for the spin partners W bJ with the coupling constant g b varied from 0 (the results of the purely contact theory are naturally recovered in this limit) to its physical value quoted in Eq. (3.3). In this calculation, the masses of the B and B * mesons are fixed to their physical values [9]. One can draw several conclusions from  Figure 1. Evolution of the binding energies of the Z b 's spin partners calculated with and without π-and η-meson exchanges as functions of the mass splitting δ between the B * and B mesons. The contact terms are re-fitted for each value of δ to provide the given binding energies of the Z b and Z b states used as input -see Eq. (4.1). The binding energies of the W bJ states are defined relative to their reference thresholds quoted in parentheses. The red dotted curves correspond to the pionless (purely contact) theory -Scenario A; the blue dashed curves are obtained for the central (S-wave) part of the OPE included -Scenario B; the blue dashed-dotted lines represent the results for the full OPE, including tensor forces -Scenario C; the black solid curves show the results of the full calculation with both OPE and OEE included on top of the contact interactions -Scenario D. The physical limit corresponds to the right edge of the plots. The results are obtained with the sharp cutoff Λ = 1 GeV in the integral equations (3.12). The uncertainty caused by the residual Λ-dependence of the equations can be estimated using the results presented in Table 1. hand, starting from g b 0.3, the dashed-dotted line starts to deviate from the dashed line indicating that the role of the tensor forces increases fast thus providing a substantial shift in the binding energy in the physical limit for the g b -this effect is best seen for the tensor partner W b2 .

Dependence on the input for the Z b 's binding energies
Finally, we investigate the dependence of the binding energies of the spin partners W bJ 's on the input used for binding energies of the Z b 's. To simplify the presentation of the results, we choose coinciding binding energies of the latter and vary them in a sufficiently wide range compatible with the values found in the literature. The results for the binding energies of the spin partners W bJ 's are given in Fig. 3. To guide the eye, the constraint  Figure 2. The binding energies of the Z b 's spin partners as functions of the coupling constant g b for the physical mass splitting δ between B * and B mesons. The binding energies of the Z b and Z b states are used as input -see Eq. (4.1). The notation of curves is the same as in Fig. 1. The results are obtained with the sharp cutoff Λ = 1 GeV in the integral equations (3.12). Figure 3. Dependence of the binding energies of the W bJ 's on the input used for the binding energies of the Z b 's for the purely contact interaction (the first plot) and for the full theory, including the OPE and OEE (the second plot). To guide the eye, the (equal) binding energies of the Z b 's are shown as the grey dashed lines. The results are obtained with the sharp cutoff Λ = 1 GeV in the integral equations (3.12).
→ 0 which provides a smooth matching between the case of the Z b 's as bound states and as virtual levels. The results for the widths of the partners are shown in Fig. 4. It is instructive to note that, in agreement with Eq. (2.7), in the strict HQSS Figure 4. The decay width of the W b2 (2 ++ ) state from coupled-channel transitions B * B * → BB and B * B * → BB * evaluated in the full theory, including the OPE and OEE.
, all six partner states are strictly degenerate. Therefore, deviations of the curves for the spin partner states W bJ 's from the grey dashed line in Fig. 3 illustrate the importance of the HQSS breaking corrections related to the nonvanishing B * -B mass splitting. In particular, had the strict HQSS limit been realised in nature all six spin partner states would have become unbound in unison for Actually, as it is seen from the left plot in Fig. 3, a similar pattern is observed in the purely contact theory beyond the strict HQSS limit. On the other hand, in the full theory with pions and η-mesons (see the right plot in Fig. 3), the 1 ++ and 2 ++ partners survive as bound states in the limit The width of the 2 ++ state is expected at the level of a few MeV -see Fig. 4 -so that one should be able to resolve it from the B * B * threshold. Meanwhile, the poles for both 0 ++ partners move to the unphysical Riemann sheets.
We note that the inclusion of the OPE interaction leads to selfconsistent results only if all relevant partial waves as well as particle channels which are coupled via the pionexchange potential are taken into account. Conversely, the partial neglect of the coupledchannel dynamics not only significantly affects the predictions for the partner states but also leads to the regulator(cutoff)-dependent results already in the strict HQSS limit. This is in full agreement with the results in the c-quark sector reported in Ref. [22].
In agreement with the discussion in Sec. 2, the results of the full coupled-channel problem beyond the HQSS limit may also exhibit some cutoff dependence which is however quite mild and which can be absorbed into a redefinition of higher-order HQSS violating contact interactions. To illustrate that the effect of the cutoff variation on our results is minor, in Table 1 we estimate the uncertainty in the binding energies and widths of the partner states when the cutoff is varied in the range from 800 to 1500 MeV.
Note that, in the vicinity of the peak, the energy dependence of the line shape for the 2 ++ partner state has a clear Breit-Wigner form, from which the parameters quoted in Table 1 were extracted. On the contrary, since the Z b state resides very close to the B * B * threshold, the Breit-Wigner distribution in the resonance region is strongly distorted by threshold effects. Then, to arrive at the quantity that can be compared to the width from -----6.2 ± 1.1 Table 1. The binding energies and the widths of the spin partners for the two sets of the Z b and Z b binding energies used as input parameters. The uncertainty in the results is due to the variation of the cutoff in the Lippmann-Schwinger equations from 800 to 1500 MeV. experiment one would need to convolute the energy distribution from the coupled-channel amplitudes with the resolution function, integrate it over the energy bins and then analyse with the standard Breit-Wigner techniques. Since this goes beyond the scope of the current work, we do not quote the width of the Z b in Table 1.

Comment on the pionless theory
In the previous section, we demonstrated that nonperturbative pion exchange has a significant impact on the pole locations of the spin partner states for the Z b and Z b . In this section, we show that it plays an additional important role to identify correctly the realistic solution from the pair of solutions present in the purely contact theory as long as input from the 1 +− channel is used to fix the parameters. Indeed, in the potential (2.4) written in terms of the low-energy constants C 1 and C 1 , an interchange of the latter only changes the sign of the off-diagonal elements which has no effect on the observables as long as Eq. (2.4) is the only contribution to the potential. Therefore, in the purely contact theory there are two solutions that both lead to the same masses of the Z b states. Accordingly, there is an alternative solution to the one given in Eq. (2.7) for the spin partner states in the strict HQSS limit, namely The two solutions presented in Eqs. (2.7) and (5.1) and their evolution for finite δ's are different in a nontrivial way. This is demonstrated in Fig. 5 where, for definiteness, we used Eq. (4.1) as input for the Z b 's. 8 However, as soon as the pion exchange is added, the strength of the off-diagonal elements changes and the system gets sensitive to the interchange C 1 ↔ C 1 . It turns out that, in the presence of the OPE with its physical strength (3.3), the Z b 's binding energies can be fixed to their physical values only for the first solution which in the strict HQSS limit corresponds to Eq. (2.7) and which is discussed in the previous section. Thus, the OPE naturally removes an ambiguity from the predictions for the spin partner states.  Figure 5.
Binding energies of the W b2 (2 ++ ) and W b0 (0 ++ ) states versus the mass splitting δ for two different solutions for the contact terms in the pionless theory: the red dotted lines correspond to the first solution discussed in Sec. 4 while the blue solid lines represent the results of the alternative scenario with the interchanged contact terms (see the text for the further details).

Conclusions
In this paper, we study the spin partners of the isovector bottomonium-like states Z b (10610) and Z b (10650) in the molecular model. For definiteness, we assume both Z b states to be shallow bound states and treat their binding energies as input parameters of the model. This allows us to fix two low-energy constants describing the leading-order short-range interactions in the 1 +− channel. Then HQSS is used to construct the contact potentials for the spin partner states of the Z b and Z b with the quantum numbers J ++ (J = 0, 1, 2), conventionally denoted as W bJ (J = 0, 1, 2). Motivated by the important role played by the nonperturbative pion dynamics in nuclear chiral EFTs -especially in the deuteron channel which demonstrates certain similarities to the Z b states (see, for example, Refs. [25,26]) -and in the c-quark sector [22], we performed a fully nonperturbative calculation with the one-pion exchange interaction added to the short-range terms. In addition, we include one-η exchange, also iterated to all orders for the sake of completeness, to see if it indeed diminishes the effect of the OPE as claimed in the literature [24].
Our results revealing the influence of the OPE and OEE interactions on the properties of the bottomonium-like systems under study can be summarised as follows.
• In contrast to the charmonium-like systems -see, for example, Refs. [22,30] three-body effects in the bottomonium sector play basically no role and, in analogy with the N N problem, they can be treated perturbatively. The reason why the static approximation for the OPE works well in the bottomonium-like systems is twofold: first, the B * -B mass splitting is significantly smaller than the pion mass; second, the masses of the B mesons are large enough to make the typical three-body momenta (proportional to m π /m × soft scale -see, for example, Refs. [43,44] for the details) sufficiently small compared to the other relevant scales of the problem such as the binding momenta and the pion mass.
• D waves (particularly the S-D transitions) play an important role and cannot be disregarded. On the contrary, if only the central S-wave part of the OPE interaction is retained, the problem is practically indistinguishable from the purely contact one after the parameters of the system are re-adjusted to place the Z b poles, used as input, to the prescribed locations. Since the leading non-trivial effect from the pions stems from the S-D tensor forces, the perturbative inclusion of the OPE is not sufficient.
• In agreement with the results of Ref. [22], the inclusion of the OPE interaction leads to selfconsistent results only if all relevant partial waves as well as particle channels which are coupled via the pion-exchange potential are taken into account. Due to the larger b-quark mass compared to the c-quark mass, consequences of uncontrolled omissions of partial waves are not that severe in the bottomonium sector. However, any partial neglect of this kind results in a strong cutoff dependence which reduces noticeably the predictive power of the approach.
• The most important effect of the coupled-channel dynamics reveals itself in the 2 ++ channel where, due to the B * B * ↔ B ( * )B transitions driven by the S-D tensor forces from the OPE, the corresponding partner state W b2 acquires a width at the level of a few MeV and it remains bound even in the limit of vanishing binding energies of the Z b states.
• One-η exchange, iterated to all orders, plays a minor role for the problem at hand slightly diminishing the effect of the OPE.
We, therefore, confirm the results of our previous studies of the OPE interaction which demonstrate its sizeable effect on the properties of the near-threshold resonances -see, for example, Refs. [22,30]. On the other hand, our findings suggest that the claim made in Ref. [45] that the influence of the static S-wave OPE on the line shapes of the Z b states can be as large as 30% is not correct. On the contrary, one expects that also for the line shapes the effect of the central S-wave OPE interaction can be largely absorbed into re-definition of the parameters of the model. However, we delegate a detailed investigation of the line shapes, which calls for an inclusion of the inelastic channels, to a later work.
The predicted masses of the spin partners of the Z b states are quoted in Table 1. In the same table we give the width of the W b2 state due to the coupled-channel transitions to the open-flavour channels BB and BB * . From this table one can conclude that three sibling states with the quantum numbers 0 ++ , 1 ++ , and 2 ++ are expected to exist near the BB, BB * , and B * B * thresholds, respectively, and to produce visible peak-like structures in the line shapes. A similar 0 ++ state near the B * B * threshold has the pole located on a remote Riemann sheet. It remains to be seen whether or not this state can be observed experimentally.
The probably most important finding of this work is that we predict the existence of an isovector 2 ++ tensor state lying a few MeV below the B * B * threshold. Unlike its scalar partner states, this tensor state should reveal itself as a bound-state peak in the B * B * line shape even in the limit of vanishing binding energies of the Z b states. The width of this structure is expected at the level of a few MeV so in principle one should be able to resolve it in the experiment.
It should be stressed, however, that in this work we did not include any inelastic channels. The W b2 width just reported -see Table 1 -stems from the possible decays of the W b2 to the BB ( * ) final states in D waves. Therefore, the total width of this state should be somewhat bigger due to its decays to, for example, Υ(nS)ρ, as proposed in Ref. [21], or to χ b1 π and χ b2 π. However, since inelastic channels play only a minor role for the Z b 's, it is also natural to expect their contribution to be small enough for the W b2 which must make this state resolvable from the threshold and, therefore, detectable in the line shape.