The strong attractive interaction of the φ meson and the proton has recently been reported by the ALICE Collaboration. The corresponding scattering length f0 is given as ${\rm{Re}}({f}_{0})=0.85\pm 0.34({\rm{stat}})\pm 0.14({\rm{syst}})$ and ${\rm{Im}}({f}_{0})=0.16\pm 0.10({\rm{stat}})\pm 0.09({\rm{syst}})$ fm. The fact that the real part is significant in contrast to the imaginary part indicates a dominant role of elastic scattering, whereas the inelastic process is less important. In this work, such scattering processes are inspected based on a unitary coupled-channel approach inspired by the Bethe–Salpeter equation. The φp scattering length is calculated based on this approach, and it is found that the experimental value of the φp scattering length can be obtained only if the attractive interaction of the φ meson and the proton is taken into account. A significant outcome of such attractive interaction is a two-pole structure in the φp scattering amplitude. One of the poles, located at (1969 − i283) MeV might correspond to N(1895)1/2− or N(1875)3/2− as listed in the review of the Particle Data Group (PDG). The other one, located at 1949 − i3 MeV should be a φN bound state, which has no counterpart in the PDG data.
Bao-Xi Sun, Ying-Ying Fan, Qin-Qin Cao. The φp bound state in the unitary coupled-channel approximation[J]. Communications in Theoretical Physics, 2023, 75(5): 055301. DOI: 10.1088/1572-9494/acc31d
1. Introduction
The strong interaction of the φ meson and the proton at low energies, due to its non-perturbative attribute, has led to rich structures of matter as reflected by an abundance of hadronic states. Understanding the nature of those hadrons has posed a long-lasting challenge to physicists. Among the attempts to tackle this challenge, chiral perturbation theory provides a systematic approach that effectively accounts for the interaction among hadrons. By asserting the unitarity condition, the chiral interaction between hadrons results in scattering amplitudes possessing poles in unphysical Riemann sheets in the complex energy plane. Such poles can be connected with heavier hadronic states as listed in the review of the Particle Data Group [1]. This strategy has had great success in interpreting a host of resonances as dynamically generated from the chiral interaction between the baryon and meson ground states, i.e. the spin 1/2 baryon octet and the pseudoscalar Goldstone-boson octet under the SU(3) flavor symmetry [2–9]. Noticeably, the chiral interaction between the spin 3/2 decuplet and the Goldstone bosons is also found to be important to understanding some higher resonances [10].
Besides the above hadronic degrees of freedom, previous works have shown that vector-meson degrees of freedom may be of crucial significance in understanding the hadron spectra [11–18]. Remarkably, a recent result of the heavy-ion collision experiment at ALICE indicates an attractive strong interaction potential between the φ meson and the proton [19]. But on the other hand, it is troublesome to accommodate those vector-meson degrees of freedom in an effective Lagrangian due to their intermediate masses and their unique decay patterns. Different models have been proposed for this enterprise [20]. Among them, the hidden-gauge symmetry model provides plausible accounts of several important phenomenological laws involving the vector mesons [20–24]. Based on this approach, the interaction of vector mesons with baryon octet was studied in detail in reference [15]. This work suggested a vanishing elastic potential of the φ meson with the proton. This conclusion is, however, in contradiction with the recent ALICE result. From the experiment, the strong interaction between the φ meson and the proton is shown to be attractive as indicated by the corresponding scattering length as given by the Lednicky–Lyuboshits fit. The real and imaginary parts of the scattering length read, respectively,
Apparently, the imaginary part of the scattering length is negligible given the error budget. This fact indicates that the elastic scattering is playing a dominant role in the interaction between φ and the proton, whereas the inelastic scattering turns out to be unimportant. Interestingly, such an attractive potential may be attributed to the K+K−p three-body correlation as suggested by reference [25]. But for the time being, the dynamical aspects of such a vector-meson-proton system are still far from being comprehended and urgently require more studies.
The present work is dedicated to the further study of the dynamics of the φ-proton interacting system using an effective Lagrangian approach. Both the vector-meson nonet and baryon ground-state octet under flavor SU(3) are incorporated in the effective Lagrangian. Based on the Lagrangian, we scrutinize the coupled-channel effects of such a system in more depth. Especially, we will explore the impacts of an attractive φ-proton potential on the couple-channel dynamics.
In the following section, we first discuss the construction of the effective Lagrangian. The vector-meson degrees of freedom will be introduced in accordance with the hidden-gauge formalism. Next, we will discuss the coupled-channel approach with a restoration of the two-body unitarity. The latter is implemented by a bubble-loop resummation inspired by the Bethe–Salpeter equation. Such a resummation involves a potential that can be obtained directly from the effective Lagrangian. However, given our current leading-order approach, the effective Lagrangian will not lead to a direct φ-proton potential, which should be significant and attractive as suggested by ALICE. We will incorporate such effect by proposing a Yukawa-type potential in addition to the one directly obtained from the Lagrangian. The corresponding discussion is referred to in section 4. This work will then end with an application of our approach to various two-body coupled-channel systems involving the vector-meson nonet and the baryon octet. Special attention will be paid to a possible resonance dynamically generated in the φN scattering channel.
2. Hidden-gauge formalism
The interaction of the vector meson with the octet of baryons can be involved according to the hidden-gauge formalism. The Lagrangian of vector mesons can be written as
with $g=\tfrac{{M}_{V}}{2f}$ and the pion decay constant f = 93 MeV [14, 15].
The Lagrangian of vector mesons in equation (2) supplies an interacting vertex of three vector mesons in the interaction of the vector meson with the octet of baryons, which comes from
With the interactions in equations (4) and (6), the t-channel interaction of the vector meson and the baryon can be constructed, where a vector meson is exchanged between them. Since the scattering process is studied on the energy region near the φ-proton threshold, the momentum of the exchanged vector meson can be neglected, and thus only the mass term is left in the propagator. Therefore, the potential of the vector meson with the octet of baryons takes the form of
where ${k}^{0}(\vec{\epsilon })$ and ${k}^{{\prime} }{}^{0}({\vec{\epsilon }}^{{\prime} })$ are energies (polarization vectors) of the incoming and outgoing vector mesons, respectively, and the Cij coefficient values can be obtained according to the Lagrangian in equations (4) and (6). Furthermore, due to the mixing of the octet and the singlet of vector mesons in the SU(3) space, the physical states of the ω and φ mesons can be written as
Since the singlet state ω1 does not couple to other vector mesons, and only the octet state ω8 gives a contribution to the amplitude, the Cij coefficients of the ω and φ mesons must be multiplied by factors $\sqrt{\tfrac{1}{3}}$ and $-\sqrt{\tfrac{2}{3}}$, respectively.
In the channel of isospin I = 1/2, the isospin states for ρN, ωN, φN, K*Λ, K*Σ can be written as
respectively. We have used the phase convention ρ+ = − ∣1, 1 > and Σ+ = − ∣1, 1 > for the isospin states in equations (10) and (14), which is consistent with the structure of the Vμ and B matrices. It is apparent that the interaction between isospin states with isospin orientation Iz = − 1/2 would generate the same resonances, so only the isospin states with Iz = 1/2 are discussed in this section.
The coefficients Cij for the fixed strangeness S = 0 and isospin I = 1/2 are listed in table 1.
Table 1. Coefficients CijIS for the sector I = 1/2, S = 0.
ρN
ωN
φN
K*Λ
K*Σ
ρN
2
0
0
$\tfrac{3}{2}$
$-\tfrac{1}{2}$
ωN
0
0
$-\tfrac{3}{2}\tfrac{1}{\sqrt{3}}$
$-\tfrac{3}{2}\tfrac{1}{\sqrt{3}}$
φN
0
$-\tfrac{3}{2}\left(-\sqrt{\tfrac{2}{3}}\right)$
$-\tfrac{3}{2}\left(-\sqrt{\tfrac{2}{3}}\right)$
K*Λ
0
0
K*Σ
2
3. Unitary coupled-channel approach
A full scattering amplitude can be expressed in an integral representation according to the Bethe–Salpeter equation. By applying the on-shell condition to the potential V involved, the equation can be reduced to a resummation of bubble-loop series and we end up with the scattering amplitude of φp → φp, reading as
with G being the bubble-loop function of the intermediate stable vector meson and baryon. According to this approach, coupled-channel effects are well incorporated. In the dimensional regularization scheme, the loop function G takes the form of
where the loop function follows the form in reference [6], and μ = 630 MeV and ai(μ) = −2.0 in the calculation, and the two-body unitarity is therefore preserved in equation (15).
Since the vector mesons, particularly the ρ and the K*, are rather broad, a proper account for their widths is required when performing the loop integral. We follow reference [15] and convolute the vector-meson-baryon loop function with the mass distribution of such a meson, i.e. by replacing the G function appearing in equation (16) by $\tilde{G}$
where m1 and Γ1 are the mass and width of the vector meson in the loop. We only take into account the widths of the ρ and the K*. In the case of the ω or φ, one or both of the kernels of these integrals will reduce to a delta function $\delta ({\tilde{m}}^{2}-{M}^{2})$. The ${\tilde{{\rm{\Gamma }}}}_{i}$ function is energy dependent and has the form of
m1 = mπ, m2 = mK and $m={m}_{{K}^{* }}$ for the K*, where λ is the Källen function, $\lambda (x,y,z)={\left(x-y-z\right)}^{2}-4{yz}$, and Γ0 is the nominal width of the ρ or the K*.
at the pole position $\sqrt{{s}_{0}}$ in the complex energy plane approximately. Therefore the coupling to the ith channel gi can be calculated with the residue of the scattering amplitude Tii at the pole position. After gi is obtained, the couplings to the other channels gj are easy to be evaluated with the residues of inelastic scattering amplitudes Tij.
When the Bethe–Salpeter equation is solved in the unitary coupled-channel approximation, the vertices of the baryon octet and vector mesons are used, and the outlines in Feynman diagrams are cut off, which implies the polarization vector of mesons is eliminated in the calculation. More detailed discussions on this question can be found in reference [26]. In the present work, we adopt the original formula in references [13–15], where $\vec{\epsilon }{\vec{\epsilon }}^{{\prime} }=-1$ in equation (8) is assumed. Actually, this scheme is consistent with the formula discussed in reference [26]. As for the spin of the nucleon, the interaction would vanish if the spin orientations of the initial and final nucleons are different from each other. Therefore, the vertices used in the Bethe–Salpeter equation are independent of the spin of the system. At this point, it is different from the case in the chiral SU(3) quark model [27].
4. Yukawa-type potential of the φ meson with the proton
The Yukawa-type potential of the φ meson with the proton takes the form of
with A = 0.021 ± 0.009 ± 0.006(syst) and α = 65.9 ± 38.0(stat) ± 17.5(syst) MeV[19]. Therefore, the φ-proton potential in the momentum space can be obtained with a Fourier transformation of equation (23), which is written as
with $\tfrac{{g}_{\phi N}^{2}}{4\pi }=A$. Apparently, the potential in equation (24) is a nonrelativistic form supposing that the three-momentum of the proton(phi-meson) is far lower than the proton(phi-meson) mass. In this case, the zero component of the momentum transfer in the phi-proton t-channel interaction tends to zero in the nonrelativistic approximation. More detailed iterations can be found in the appendix of reference [28]. It should be noticed that the phi–proton coupling in this work ${g}_{\phi N}^{2}/4\pi =A$ is different from that in reference [19], where ${g}_{\phi N}^{2}=A$.
When the Bethe–Salpeter equation is solved, the Yukawa-type potential in equation (24) is multiplied by a mass of φ meson. and the three-momentum transfer $\vec{q}$ vanishes at the φ- proton threshold. Therefore,
which can be used as the potential involved in the Bethe–Salpeter resummation.
In the calculation of the present work, we set A = 0.021 and α = 65.9 MeV, and the corresponding potentials of the vector meson and baryon octet are depicted in figure 1. There are no direct interactions of ωN, φN, and K*Λ when the hidden-gauge symmetry is taken into account. However, the attractive φN potential in equation (25) is a constant and about half of the ρN potential, and it is no doubt that it would play a critical role in the sector of strangeness S = 0 and isospin I = 1/2.
Figure 1. The potential of the vector meson and baryon octet as a function of the total energy of the system ECM for elastic processes with strangeness S = 0 and isospin I = 1/2. The φN, ρN, and K*Σ cases are represented by solid, dashed, and dotted lines, respectively, while the ωN and K*Λ potentials are zero.
In reference [19], it is assumed that the elastic process of the φ meson and the proton plays a dominant role, while the inelastic process is not important.
5. φN bound state
With the coefficients listed in table 1, the scattering amplitude T of the φN → φN reaction is evaluated by solving the Bethe–Salpeter equation in the unitary coupled-channel approximation, and then the φN scattering length can be obtained according to the φN elastic scattering amplitude at the threshold, i.e.
With μ = 630 MeV and al(μ) = −2.0, the real part of the φN scattering length is about −0.15 fm, and the imaginary part is 0.11 fm approximately. Apparently, the φN scattering length takes a small value and can not be consistent with the value announced by the ALICE Collaboration recently, as given in equation (1). Even if the subtraction constants are adjusted, it is still far from the experimental value.
In order to obtain the φN scattering length consistent with the experimental value supplied by the ALICE Collaboration, we include the attractive φN potential as given by equation (25) when performing the Bethe–Salpeter resummation. By adjusting the subtraction constant values, the φN scattering length is obtained, and it is consistent with the experimental value in the range of uncertainty. The subtraction constant values and the corresponding φN scattering length are listed in table 2.
Table 2. Subtraction constants and the real and imaginary parts of the corresponding φN scattering length obtained in the unitary coupled channel in units of megaelectronvolts, and the regularization scale μ = 630 MeV.
aρN
aωN
aφN
${a}_{{K}^{* }{\rm{\Lambda }}}$
${a}_{{K}^{* }{\rm{\Lambda }}}$
f0
−2.0
−2.0
−2.4
−1.9
−1.8
0.86 + i0.19
With the subtraction constants listed in table 2, the scattering amplitudes ∣Tii∣2 of the vector meson and the baryon octet are calculated in the unitary coupled-channel approximation. In figure 2, the scattering amplitude with strangeness S = 0 and isospin I = 1/2 is displayed. In this case, a peak around 1950 MeV is visible in the φN, K*Λ, and K*Σ channels. Since it is below the φN threshold, it could be associated with a φN bound state. Actually, there is another pole of ∣Tii∣2 detected around 1970 MeV in the K*Λ and K*Ξ channels. However, it is far from the real axis, and therefore leads to no visible effects as shown in figure 2. Since it is higher than the φN threshold, it can be regarded as a φN resonance.
Figure 2. The amplitude squared ∣T∣2 as a function of the total energy of the system ECM for elastic processes with strangeness S = 0 and isospin I = 1/2. The φN, K*Λ, and K*Σ cases are represented by solid, dashed, and dotted lines, respectively, while the amplitudes squared of ρN and ωN are small and the corresponding curves are invisible in this figure.
The pole positions and corresponding couplings of these two states to the different channels, obtained from the residues at the poles are shown in table 3. Since the inelastic scattering amplitudes are far lower than the elastic ones, all couplings are calculated according to elastic processes. It can be seen that the state at 1949 − i3 MeV couples strongly to the φN channel, and the state at 1969 − i283 MeV mainly couples to the K*Λ and K*Σ channels.
Table 3. Pole positions in units of megaelectronvolts and coupling constants to various channels of resonances with strangeness S = 0 and isospin I = 1/2.
Pole positions
ρN
ωN
φN
K*Λ
K*Σ
1949 − i3
0.0 + i0.0
0.1 + i0.0
2.1 + i0.1
1.6 + i0.3
1.8 + i0.0
1969 − i283
0.1 + i0.1
0.0 + i0.2
0.1 − i0.1
0.3 − i0.4
0.2 − i0.0
The Yukawa-type potential of the phi-proton system is only valid when the total energy of the system is close to the phi-proton threshold. At this point, the three-momentum transfer in the phi-proton potential is not large and can be neglected in the calculation, and a phi-nucleon bound state is found when the scattering amplitude is analyzed. Since the bound state is about 10 MeV lower than the phi-nucleon threshold, it is reasonable to neglect the three-momentum transfer in the phi-proton potential in equation (24) when the Bethe–Salpeter equation is solved.
Actually, we also studied the case in which the three-momentum transfer $\vec{q}$ is not zero. By solving the Bethe–Salpeter equation, it is found that two poles appear at 1949 − i3 MeV and 1970 − i283 MeV on the complex energy plane, respectively. It is apparent that these results are almost the same as those with $\vec{q}=0$, where the three-momentum transfer is neglected. Moreover, the same value of the phi-nucleon scattering length is obtained with the subtraction constants listed in table 2. The coupling constants to different channels are calculated again, and it shows that their values change slightly. However, these two states still couple strongly to the φN, K*Λ, and K*Σ channels.
In the S-wave approximation, two poles of the scattering amplitude with strangeness S = 0 and isospin I = 1/2 are found in the unitary coupled-channel approximation, and they are degenerate in JP = 1/2−, 3/2−. The resonance state at 1969 − i283 MeV might correspond to either the baryon N(1895)1/2− or N(1875)3/2− particles, while the pole at 1949 − i3 MeV has no counterpart in the PDG data since its decay width is too small. It might be a φN bound state caused by the strong attractive interaction between the φ meson and the proton. The pole position and corresponding PDG data are listed in table 4.
Table 4. Pole position and corresponding PDG data with strangeness S = 0 and isospin I = 1/2, all in units of megaelectronvolts.
Pole position
name
JP
status
mass
width
1949 − i3
—
—
—
—
—
1969 − i283
N(1895)
1/2−
****
1870–1920
80–200
N(1875)
3/2−
***
1850–1920
120–250
The pole position in the complex energy plane is closely related to subtraction constants when the Bethe–Salpeter equation is solved. With the subtraction constant a = −2 for variant channels, a pole at 1979 − i56 MeV is generated in the sector of strangeness S = 0 and isospin I = 1/2, which is similar to the result in reference [15]. However, in this work, the subtraction constants are determined in order to obtain the experimental phi-nucleon scattering length, as listed in table 2. Therefore, a pole at 1969 − i283 MeV is generated, which is higher than the threshold of the phi-meson and the nucleon, and can be regarded as a phi-N resonance state. While the particles N(1895)1/2− and N(1875)3/2− is under the phi-N threshold, might be bound states of the phi-meson and the nucleon.
The discrepancy in the calculation results from the experimental data might come from the SU(3) hidden-gauge symmetry, by which the vector mesons are included in the interaction Lagrangian. Anyway, the masses of vector mesons are different from each other, so the SU(3) hidden-gauge symmetry is only reliable approximately. Moreover, perhaps some high-order corrections on the kernel should be taken into account when the Bethe–Salpeter equation is solved in the unitary coupled-channel approximation. These topics would be investigated continuously in our future work.
The φ-nucleon bound state is also studied in the chiral SU(3) quark model [27], where the exchange of scalar meson σ is dominant in the attractive potential of the phi-meson and the nucleon. In the phi-nucleon potential in equation (24), the three-momentum transfer is negligible in the calculation, and thus only the ratio ${g}_{\phi N}^{2}/{\alpha }^{2}$ is critical in the calculation, which equals 6.0 × 10−5 MeV−2 approximately. In the chiral SU(3) quark model, the attractive potential via a scalar meson can be obtained when gφN and α in equation (24) are replaced with the coupling constant gch and the mass of the scalar meson σ in reference [27], respectively. In the region of the phi-nucleon threshold, the three-momentum transfer can be neglected and the value of ${g}_{{ch}}^{2}/{m}_{\sigma }^{2}\approx 2.4\times {10}^{-5}$ MeV−2 is obtained in the extended chiral SU(3) quark model of type II, which has the same order of magnitude as that of ${g}_{\phi N}^{2}/{\alpha }^{2}$ in this work. Moreover, the bound state is about 10 MeV lower than the phi-nucleon threshold, and it is consistent with the binding energy of the phi-nucleon bound state obtained in the chiral SU(3) quark model [27].
6. Summary
In the present work, we investigate the implications of a strong attractive potential of the φ meson and proton as recently reported by the ALICE Collaboration. An effective description of such a potential is obtained by introducing a Yukawa-type formalism on top of the vector-meson-nonet and baryon-octet interaction given by the effective Lagrangian. We deliberate on the impacts of such a potential to the coupled-channel effects of vector-meson and baryon scattering systems. A unitary approach has been employed to implement the coupled-channel effects. This approach involves a series resummation reduced from Bethe–Salpeter equation, and we arrive at a φN scattering length consistent with the experiment. Based on such potential, we investigate the coupled-channel effect, and special attention has been paid to the coupled-channel effects of the vector-meson and baryon scattering system with strangeness S = 0 and isospin I = 1/2. Two poles have been found in the complex energy plane of the corresponding scattering amplitude. One pole locates at $\sqrt{s}=(1969-i283)$ MeV, and might correspond to either the baryon N(1895)1/2− or N(1875)3/2−. In addition, we predict a possible φN bound state, which has not been found in the PDG. This state leads to another pole structure at $\sqrt{s}=(1949-i3)$ MeV in the scattering amplitude.
Acknowledgments
Bao-Xi Sun would like to thank Xiao-Yu Guo for reading and revising the whole manuscript. Moreover, we appreciate the constructive and helpful advice from the referees of Communications in Theoretical Physics.
DanilkinI VLutzM F MLeupoldSTerschlusenC2013 Photon-fusion reactions from the chiral Lagrangian with dynamical light vector mesons Eur. Phys. J. C73 2358
Del GrandeRŠerkšnytėLFabbiettiLSartiV MMihaylovD2022 A method to remove lower order contributions in multi-particle femtoscopic correlation functions Eur. Phys. J. C82 244