Mass of Y(3940) in Bethe-Salpeter equation for quarks

The general form of the Bethe-Salpeter wave functions for the bound states composed of two vector fields of arbitrary spin and definite parity is corrected. Using the revised general formalism, we investigate the observed \emph{Y}(3940) state which is considered as a molecule state consisting of $D^{*0}\bar{D}^{*0}$. Though the attractive potential between $D^{*0}$ and $\bar{D}^{*0}$ including one light meson ($\sigma$, $\pi$, $\omega$, $\rho$) exchange is considered, we find that in our approach the contribution from one-$\pi$ exchange is equal to zero and consider SU(3) symmetry breaking. The obtained mass of \emph{Y}(3940) is consistent with the experimental value.


I. INTRODUCTION
The exotic state Y (3940) was discovered by the Belle collaboration [1] and then confirmed by the BABAR collaboration [2]. The investigation of the structure of Y is of great significance, while the conventional cc charmonium interpretation for this state is disfavored [3]. Then possible alternative interpretations have been proposed, such as hadronic molecule and tetraquark states. Following the experimental results, it is suggested in Ref.
[4] that the Y (3940) is a hadronic molecule state of D * 0D * 0 . However, in previous work [4], the numerical result of the binding energy for the molecule state sensitively depends on the value of typical cutoff in the effective interaction potential between two heavy vector mesons, and this two heavy mesons are considered as pointlike objects. Furthermore, the spin-parity quantum numbers J P of the Y (3940) are not unambiguously determined in experiment, except for C = +. The molecule state hypothesis implies that the quantum numbers of Y (3940) are 0 + or 2 + , but Ref. [4] can not deduce the definite quantum numbers in theory.
Though the general form of the Bethe-Salpeter (BS) wave functions for the bound states consisting of two vector fields of arbitrary spin and definite parity has been given in Ref.
[5], we find that the derivation of this formalism has a serious defect. In this work, the general formalism is firstly corrected. Then we assume that the Y (3940) state is a molecule state composed of two heavy vector mesons D * 0D * 0 and the revised general formalism is applied to investigate this two-body system. To construct the interaction kernel between two heavy mesons derived from one light meson (σ, π, ω, ρ) exchange, we consider that the heavy meson is not a point-like particle but a bound state composed of u-quark and cquark and then investigate the light meson interaction with the u-quark in the heavy meson.
Through the form factor we can obtain the light meson interaction with heavy meson and the potential between two heavy mesons without an extra parameter. [6,7] Obviously, this potential in our approach contains more inspirations of quantum chromodynamics (QCD).
Finally, numerically solving the relativistic Schrödinger-like equation with this potential, we can obtain the mass of the molecule state and then deduce the definite quantum numbers of the Y (3940) system.
In this work, one-π exchange is considered in the interaction kernel between two heavy mesons. When investigating this pseudoscalar meson interaction with the u-quark in the heavy meson, we find that the coupling L I = ig πū γ 5 uπ should have no contribution to this interaction and represent it as the derivative coupling Lagrangian L I = if πū γ µ γ 5 u∂ µ π. In this approach we find that one-π exchange has no contribution to the potential between two heavy vector mesons. Besides, it should be noted that the flavor-SU(3) singlet and octet states of vector mesons mix to form the physical ω and φ mesons, so the exchange-mesons between two heavy mesons should not be the physical mesons but rather the singlet and octet states. Then in the interaction kernel between two heavy mesons SU(3) symmetry breaking should be considered.
This paper has the structure as follows. In Sec. II the revised general form of BS wave functions for the bound states composed of two vector fields with arbitrary spin and definite parity is given. In Sec. III we show the BS wave functions for the molecule states of D * 0D * 0 with J P = 0 + and 2 + . After constructing the interaction kernel between two heavy vector mesons, we obtain the Schrödinger type equations in instantaneous approximation. In Sec.
IV we show how to calculate the form factors of the heavy meson D * . Then the interaction potential and the mass of Y are calculated. Our conclusion is presented in Section V.

II. REVISED GENERAL FORM OF THE BS WAVE FUNCTIONS
If a bound state of spin j and parity η P ′ is composed of two vector fields with masses M 1 and M 2 , respectively, its BS wave function is a 4 × 4 matrix which can be written as where and η 1,2 are two positive quantities such that η 1,2 = M 1,2 /(M 1 + M 2 ). Then one has the BS wave function in the momentum representation where P ′ is the momentum of the bound state, p ′ is the relative momentum of two vector fields and we have P ′ = p ′ 1 + p ′ 2 , p ′ = η 2 p ′ 1 − η 1 p ′ 2 , p ′ 1 and p ′ 2 are the momenta of two vector fields, respectively. This can be shown as Fig. 1.
Then we will reduce the general form without any assumption and approximation.
Using the subsidiary condition ∂ µ A µ (x) = 0 for the massive vector field and the equaltime commutation relation, we get The proof has been given by Ref.
[5]. The BS wave function in Eq.
(1) has this relation and in momentum representation To satisfy Eqs. (8) and (11), the BS wave function of the bound state becomes where φ i and ψ i are the linear combinations of f i (i = 1, · · · , 16) in Eq. (6), respectively.
This derivation makes use of the fact that P ′ , p ′ and η µ 1 ···µ j are linearly independent.
Now, under space reflection with x ′′ = (−x ′ , x ′ 0 ) and P ′′ = (−P ′ , P ′ 0 ). We obtain the properties of the BS wave function under space reflection from Eqs. (1) and (14): and in the momentum representation From Eqs. (4), (6), (12), (13) and (16), it is easy to derive that, for η P ′ = (−1) j , and, for η P ′ = (−1) j+1 ,  17) and (18), we conclude that the BS wave function of a bound state composed of two massive vector fields with spin j = 0 and parity (−1) j has only four independent components and that of a bound state with J P = 0 + has only two independent components.
From Eqs. (19) and (20), we conclude that the BS wave function of a bound state with spin j and parity (−1) j+1 has only five independent components except for j = 0 and one for j = 0. Up to now, all the above analyses are model independent. In the next section we will apply the general formalism to investigate molecule states composed of two vector mesons.

III. THE EXTENDED BETHE-SALPETER EQUATION
Assuming that the Y (3940) is a S-wave molecule state consisting of two heavy vector mesons D * 0 andD * 0 , one can have J P = 0 + or 2 + for this system.
[4] From Eqs. (17) and (18), we can obtain the BS wave function describing this bound state, for J P = 0 + , or, for J P = 2 + , The BS wave function of this bound state satisfies the Bethe-Salpeter equation where V αθ,βκ is the interaction kernel, the propagators for the spin 1 fields ∆ F λα (p ′ 1 ) = (δ λα + 2 +M 2 2 −iǫ and the bound state momentum is set as P ′ = (0, 0, 0, iM) in the rest frame. In Ref.
[5], we have considered that the effective interaction between two heavy mesons is derived from one light meson (σ, ω, ρ) exchange and obtained that the molecule state D * 0D * 0 lies above the threshold. In this work, one-π exchange is also considered and one-ω exchange is reconsidered, shown as in Fig. 2. Now, we construct the kernel between two heavy vector mesons from one-π exchange.
The charmed meson D * is composed of a heavy quark c and a light antiquarkū. Owing to the large mass of c-quark, the Lagrangian representing the interaction of π-meson triplet with quarks should be or where g π and f π are the π-meson-quark coupling constants, g π = 340 97 [8,9]. Because the contribution of Fig. 2 is from the term iūγ 5 π 0 u or iūγ µ γ 5 ∂ µ π 0 u, we set f π = gπ 2mu and m u is u-quark mass. Then the effective quark current is J − = iūγ 5 u or J + µ = iūγ µ γ 5 u and the S-matrix element between two heavy mesons is or where V |J − |V and V |J + µ |V represent the vertices of the pseudoscalar meson interaction with the heavy vector meson, respectively. The matrix elements of these quark currents can be expressed as , we obtain the interaction kernel from one-π exchange Then we reconsider one-ω exchange between two heavy vector mesons. In hadronic physics, the physical ω and φ mesons are linear combinations of the SU(3) octet V 8 = (uū + dd − 2ss)/ √ 6 and singlet V 1 = (uū + dd + ss)/ √ 3 as g φ = −g 8 cosθ + g 1 sinθ, g ω = g 8 sinθ + g 1 cosθ, where g ω and g φ are the corresponding meson-quark coupling constants, g 2 ω = 2.42 and g 2 φ = 13.0 [11]. Since the SU(3) is broken, the masses of the singlet V 1 and octet V 8 states are approximatively identified with the two physical masses of ω and φ mesons, respectively. The interaction kernel derived from one scalar meson exchange and one vector meson exchange has been given in Ref. [5], and the kernel from one light meson (π, σ, ρ, V 1 and where k = (k, 0).
Firstly, we assume that the Y (3940) is a molecule state with J P = 0 + . Substituting its , which are small for the masses of the heavy mesons are large. It is difficult to strictly solve the BS equation, in this paper we use a simple approach to solve it as follow. Ignoring the cross terms, one can obtain two individual equations: And solving this two equations respectively, one can obtain two series of eigenvalues and eigenfunctions. Because the cross terms are small, we can take the ground state BS wave function to be a linear combination of two eigenstates F 10 λτ and F 20 λτ corresponding to lowest energy in Eqs. (36) and (37). Then in the basis provided by F 10 , the BS wave function χ 0 + λτ is considered as Substituting Eq. (38) into the BS equation (23) and comparing the tensor structures in both sides, we obtain a eigenvalue equation where the eigenvalues are different from the eigenvalues in (36) and (37). From this equation, we can obtain the eigenvalues and eigenfunctions which contain the contribution from the cross terms.
Comparing the terms p ′ 2λ p ′ 1τ in the left and right sides of Eq. (36), we obtain where V 1 (p ′ , q ′ ; P ′ ) contains all coefficients of the term p ′ 2λ p ′ 1τ in the right side of Eq. (36). In this paper, we set k = (k, 0). Then the fourth components of momenta of two heavy mesons have no change: . To simplify the potential, we replace the heavy meson energies The potential depends on the three-vector momentum V (p ′ , q ′ ; P ′ ) ⇒ V (p ′ , q ′ , M). Integrating both sides of Eq. (40) over p ′ 0 and multiplying by (M + ω 1 + ω 2 )(M 2 − (ω 1 − ω 2 ) 2 ), we obtain the equation of the Schrödinger type and the potential between D * 0 andD * 0 up to the second order of the p ′ /M H expansion , ω 1 = p ′2 + M 2 1 and ω 2 = p ′2 + M 2 2 . And comparing the terms p ′ 1λ p ′ 1τ in both sides of Eq. (37), we obtain and the potential between D * 0 andD * 0 up to the second order of the p ′ /M H expansion In instantaneous approximation the eigenfunctions in Eqs. (36) and (37) can be calculated and the eigenvalue equation (39) becomes where the matrix elements and b 2 10 (M)/(2µ R ) and b 2 20 (M)/(2µ R ) are the eigenvalues corresponding to lowest energy in Eqs. (41) and (44), respectively; Ψ 0 + 10 and Ψ 0 + 20 are the corresponding eigenfunctions. In Eqs. (42), (45) and (47) the contribution from one-π exchange to the potential between two heavy vector mesons has vanished, but we still give the heavy meson form factors h (p) (k 2 ) and h (p) (k 2 ) in Sec. IV. Then applying the method as above, we can investigate the alternative J P = 2 + assignment for the Y state.

IV. FORM FACTORS OF HEAVY VECTOR MESONS
To calculate these heavy vector meson form factors h(k 2 ) describing the heavy meson structure, we have to know the wave function of heavy vector meson in instantaneous approximation. The heavy vector meson D * 0 is regarded as a resonance and its three-vector and n (i) are three four-vectors defined by In Eq. (50) the first term represents the light pseudoscalar meson interaction with the uantiquark in D * 0 , while the second term represents its interaction with the u-quark inD * 0 .
Substituting the vertex function Γ (1) given by Eq. (50) into the matrix element (49) and then comparing the resulting expressions with the form factor decompositions (28) and (29), Then the quark current becomes J + µ = iūγ µ γ 5 u and the matrix element between heavy meson states has the form and the corresponding vertex function is The form factors corresponding to one light pseudoscalar meson exchange are obtained where Ψ V is the wave function of heavy vector meson expressed as Eq. (48). The function Fig. 4.
Finally, we obtain the potentials between two heavy vector mesons for J P = 0 where F 1 (k 2 ) and F 2 (k 2 ) represent the form factors corresponding to one light scalar and vector meson exchange, respectively, which have been given in our previous work [5].
The lower energy should be the ground state mass of this molecule state which lies below the threshold and the calculated mass is consistent with the experimental data, while in the experiment the mass of Y (3940) is 3.943GeV [1]. This result is different from Ref.
[5], this is because in this paper we consider the SU(3) symmetry. At the same time, the molecule state D * 0D * 0 with J P = 2 + lies above the threshold. Therefore, we can have a conclusion that if the Y (3940) state is a molecule state composed of D * 0D * 0 , its quantum numbers should be J P = 0 + . Besides, we can predict that two bottom mesons B * 0B * 0 can be bound by one light meson exchange as a molecule state with J P = 0 + .

V. CONCLUSION
In this work, we assume that the exotic state Y (3940) is a D * 0D * 0 molecule state and calculate its mass. The revised general formalism of the BS wave functions for the bound states composed of two vector fields is applied to investigate this issue. Considering one-π exchange and SU(3) symmetry breaking, we construct the interaction potential between D * 0 andD * 0 through the heavy vector meson form factors describing the heavy meson structure.
The calculated mass of the molecule state with J P = 0 + is consistent with the mass of the Y state in experiment, so we conclude that if the Y (3940) state is a molecule state composed of D * 0D * 0 , its quantum numbers should be J P = 0 + . Adopting this method, we investigate the molecule states composed of two bottom mesons.    [5] X. Chen, R. Liu, R. Shi, and X. Lü, Phys. Rev. D 87, 065013 (2013).  [7] X. Chen, R. Liu, R. Shi, Y. Yin, Z. Shi, A. Yang, and X. Lü, Commun. Theor. Phys. 57, 833 (2012).