Development of Bethe–Salpeter theory for dealing with unstable system

In the framework of relativistic quantum field theory, the solution of homogeneous Bethe–Salpeter equation for two-body bound state can not describe unstable system, so we develop Bethe–Salpeter theory to investigate resonance which is regarded as an unstable two-body system. Based on Bethe–Salpeter wave function, we consider the time evolution of two-body bound state determined by the total Hamiltonian. The total matrix element for arbitrary decay channel is expressed in terms of the Heisenberg picture, and Mandelstam’s approach is generalized to calculate the matrix element between bound states with respect to arbitrary value of the final state energy. Some innovations to Feynman diagram are made so that the key features of dispersion relation can be more clearly exhibited. This new resonance theory in quantum field theory is applied to investigate exotic particle which is considered as an unstable meson–meson molecular state.


I. INTRODUCTION
Many exotic particles have been discovered in experiment and many possible alternative interpretations beyond quark-antiquark state have been proposed in theory [1][2][3][4][5][6].Among these interpretations, homogeneous Bethe-Salpeter (BS) equation is frequently used to investigate the internal structure of exotic particles which are considered as two-body bound states [7][8][9][10][11].In quantum field theory, homogeneous BS equation is a nonperturbative method [12][13][14], which should be only applied to deal with two-body bound state in the strict sense because the solution of homogeneous BS equation can not contain the contribution from decay channels.However, in experiments exotic particles are resonances, so these exotic particles are unstable states which should not be completely treated as stationary two-body bound states and it is more reasonable to regard exotic resonances as unstable two-body systems.More importantly, present field theory seldom takes into account the difference between the internal structure of stable two-body bound state and the one of unstable two-body system, and the theory describing internal structure of resonance has not been constructed in field theory.In this work, two-body bound state is a stable composite particle and resonance is regarded as an unstable composite particle.We develop homogeneous BS theory to describe the internal structure of resonance in the framework of relativistic quantum field theory and illustrate this new theory based on BS equation for exotic meson resonance.
In our previous works [8][9][10], exotic particles were considered as meson-meson bound states.Solving homogeneous BS equations for meson-meson bound states, we obtained masses and BS wave functions.The mass of meson-meson bound state was regarded as mass of exotic meson resonance and the correction for energy level of molecular state due to decay channels has not been considered [8][9][10].As well-known, all decay channels of resonance should contribute to its physical mass.However, the technique about the dynamics of coupled channels based on Schrödinger wave function still remains in nonrelativistic case [15][16][17].In relativistic quantum field theory, the Schrödinger wave function is not a useful quantity to work with, in view of its non-Lorentz-invariant character.So far, the technique about the dynamics of coupled channels based on covariant Bethe-Salpeter wave function has not been established.Therefore, it is necessary and important to seek a development of homogeneous BS theory for dealing with the dynamics of coupled channels in the framework of relativistic quantum field theory.
In this paper, exotic meson resonance is considered as an unstable meson-meson molecular state.Based on BS wave function for meson-meson bound state, we can provide a description for the prepared state and then study the time evolution of meson-meson molecular state determined by the total Hamiltonian.Using dispersion relation, the Heisenberg picture and Mandelstam's approach, we obtain the correction for energy level of resonance and then the physical mass is used to calculate its decay width.An innovative Feynman diagram is introduced, in which the key features of dispersion relation can be exhibited clearly.

II. TIME EVOLUTION OF BS WAVE FUNCTION USED TO STUDY ENERGY LEVEL AND DECAY WIDTH OF RESONANCE
Let us begin with the interaction Lagrangian for the coupling of light quark fields to light meson fields as in effective theory at low energy QCD.According to the effective theory at low energy QCD, non-vanishing vacuum condensate causes the spontaneous breaking of chiral symmetry, which leads to the appearance of Goldstone bosons [18].At low energy QCD, the effective interaction Lagrangian can be regarded as Lagrangian for the interaction of light mesons with quarks [9] where Q = (ū, d, s), Q′ = (ū, d), g represents the corresponding meson-quark coupling constant, P and V are the octet pseudoscalar and nonet vector meson matrices, respectively.
From this Lagrangian, we have investigated the light meson interaction with quarks in heavy mesons and obtained the interaction of heavy meson with light meson through the heavy meson form factor [8,19].
Using path integrals, one can obtain a homogeneous integral equation for arbitrary bound state composed of two mesons.In symbolic notation homogeneous BS equation may be written as where χ represents BS wave function, the kernel V is the sum of all irreducible graphs, S (1) and S (2) represent meson propagators, respectively.Solving this BS equation, one can obtain the mass M 0 and BS wave function χ P (x 1 , x 2 ) for this meson-meson bound state with momentum P = (P, i P 2 + M 2 0 ).We emphasize that the kernel V in homogeneous BS equation expressed as Eq. ( 2) is defined in two-body channel so V is not complete interaction.The kernel in homogeneous BS equation (2) plays a central role for making two-body system to be a stable bound state, and the solution of homogeneous BS equation (2) should only describe bound state.
Since resonance decays spontaneously into other particles, we can suppose that at the times t 1 = 0 and t 2 = 0 this unstable state has been prepared to decay.This prepared state (ps) can be described by the ground-state BS wave function which has the form where E(p) = p 2 + m 2 and η 1 +η 2 = 1.Then the time evolution of this system determined by the total Hamiltonian H has the explicit form where G(ǫ) = (ǫ − H) −1 is the Green's function and the contour C 2 runs from ic r + ∞ to ic r − ∞ in energy-plane.The positive constant c r is sufficiently large that no singularity of (ǫ − H) −1 lies above C 2 .The Green's function can be represented by scattering matrix [20] G aa (ǫ) = (χ ps a , G(ǫ) where χ ps a represents (2π) −3/2 [2E(P )] −1/2 χ P (x 1 − x 2 ) in Eq. ( 3).The proof of Eq. ( 5) has been given by Ref. [20].This work will give T aa (ǫ) in the framework of relativistic quantum field theory.In field theory the operator T (ǫ) is just the scattering matrix with energy ǫ, and T aa (ǫ) is the T -matrix element between two bound states, which should be defined as Because of the analyticity of T aa (ǫ), we define where ǫ approaches the real axis from above, D and I are the real and imaginary parts, respectively.When there is only one decay channel, we can use the unitarity of T aa (ǫ) to obtain [20] 2I where P b = (P b , iE b ) is the total energy-momentum vector of all particles in the final state and the T -matrix element T ba (ǫ) is defined as b out|a in = −i(2π) 4 δ (3) (P b − P)δ(E b − ǫ)T ba (ǫ).The delta-function in Eq. ( 7) means that the energy ǫ in scattering matrix is equal to the total energy E b of the final state, and b represents summing over momenta and spins of all particles in the final state.For E b = ǫ, we also denote the total energy of the final state by ǫ and I(ǫ) becomes a function of the final state energy.Using dispersion relation for the function T aa (ǫ), we obtain The symbol P means that this integral is a principal value integral and the variable of integration is the total energy ǫ ′ of the final state.In order to obtain the real part D(ǫ), we should calculate the function I(ǫ ′ ) of value of the final state energy ǫ ′ , which is an arbitrary real number over the real interval ǫ M < ǫ ′ < ∞.As usual the momentum of initial bound state a is set as P = (0, 0, 0, iM 0 ) in the rest frame and ǫ M denotes the sum of all particle masses in the final state.
Let us suppose that there are several decay channels and the final state b may contain n composite particles and n ′ elementary particles in decay channel c ′ .From the Heisenberg field operator corresponding to elementary particle in Lagrangian, one can obtain the infield and out-field by means of the Yang-Feldman equations.Asymptotically an elementary particle is represented by the in-field and out-field satisfying the corresponding free field equations of motion.Then the in-state and out-state consisting of elementary particle can be generated and the ordinary S-matrix element has been well defined.For the composite particle composed by elementary particles, it is clearly impossible to introduce a Heisenberg field operator corresponding to this composite particle.However, asymptotically composite particle can also be represented by the in-field and out-field satisfying the corresponding free field equations of motion [21], and then the in-state and out-state consisting of composite particle can also be generated.As far as the asymptotic condition is concerned, there is no clear distinction between elementary and composite particles in the final state b.
This work investigates the time evolution of unstable composite particle and the initial bound state a must be different from the final composite particle, so we have to consider the internal structure of composite particle in the final state b.From Eq. ( 7), we have where ..Q n are the momenta of final elementary and composite particles, respectively; P ǫ ′ = (0, 0, 0, iǫ ′ ), T (c ′ ;b)a (ǫ ′ ) is the T -matrix element with respect to ǫ ′ , spins represents summing over spins of all particles in the final state, and c ′ represents summing over all open and closed channels.In Eq. ( 9) the energy in scattering matrix is equal to the total energy ǫ ′ of the final state b, which is an arbitrary real number over the real interval The mass M 0 and BS amplitude of initial bound state a have been specified and the value of the initial state energy in the rest frame is a specified value M 0 .From Eq. ( 9), we have I(ǫ ′ ) > 0 for ǫ ′ > ǫ M and I(ǫ ′ ) = 0 for ǫ ′ ǫ M , which is the reason that the integration in dispersion relation (8) ranges from ǫ M to +∞.Therefore, in order to obtain the real part D(ǫ) due to all open and closed channels, we have to calculate with respect to value of the final state energy ǫ ′ , which is an arbitrary real number over the In our theory T (c ′ ;b)a (ǫ ′ ) must involve bound state, so this matrix element can not be calculated as an ordinary S-matrix element.Because it is impossible to introduce a Heisenberg field operator corresponding to composite particle, we define the in-field and out-field to represent composite particle.For example, we suppose that the initial bound state is a scalar composite particle composed of two elementary particles.Asymptotically this composite particle can be represented by in-field φ M 0 in (X) satisfying the corresponding free field equation of motion [21] ( where M 0 is the mass of bound state, X is the centre of mass coordinate and X = η 1 x 1 +η 2 x 2 .
Then we can define the creation operator a M 0 † P,in in φ M 0 in , which generates the in-state |P in = a M 0 † P,in |0 .Similarly, n composite particles in the final state can be generated by n out-fields, respectively, and the final state including n composite particles is out-state Q 1 ...Q n out|.
Therefore, the initial bound state a is |P in and the final state b including n composite particles and n ′ elementary particles is Using the Heisenberg picture, we can obtain the total matrix element between a final state Q Here these Heisenberg fields φ represent the boson fields appearing in Lagrangian.The functions f are solutions to the corresponding free field equations of motion, and these S ′ represent free boson propagators.In Eq. ( 11), we have used reduction formulae to deal with elementary particles and of great interest is the matrix element of a time-order product of Heisenberg field operators between bound states.Mandelstam's approach is a technique based on BS wave function for evaluating the general matrix element between bound states.Applying Mandelstam's approach, one can express the general matrix element between bound states in terms of BS wave functions and a two-particle irreducible Green's function, and the proof has been given by Ref. [21].In this work, Mandelstam's approach is generalized to evaluate the bound state matrix element with respect to arbitrary value of the final state energy ǫ where T(y 1 ...y 2n ; z 1 ...z i ...z n ′ ; x 1 , x 2 ) is the two-particle irreducible Green's function, χ and χ are BS wave functions for the final and initial bound states, respectively.The function T can, in principle, be evaluated by means of perturbation theory.It is necessary to emphasize that the general matrix element ( 12) is calculated with respect to ǫ ′ , and the energy in T is equal to the final state energy ǫ ′ , which is an arbitrary real number over the real interval , we obtain T (c ′ ;b)a (ǫ ′ ).To illustrate this, we imagine that the initial bound state (MS) is composed of two heavy vector mesons (V M and V M ) and the final state contains a heavy meson (HM) and a light meson (LM).If a bound state with spin j is created by two massive vector fields, its BS wave function can be defined as χ j P (λτ ) (x 1 , x 2 ) = 0|T A λ (x 1 )A † τ (x 2 )|P, j and we have given the general form for this BS wave function χ j λτ (P, p) in the momentum representation, where p is the relative momentum of two vector fields [8,10].This BS wave function should satisfy the equation where V θθ ′ ,κ ′ κ is the interaction kernel, p ′ 1 and p ′ 2 are the momenta carried by two vector fields, ∆ F λθ (p ′ 1 ) and ∆ F κτ (p ′ 2 ) are the propagators for the spin 1 fields.Owing to the effective interaction Lagrangian at low energy QCD (1), we have to consider that the heavy meson is a bound state composed of a quark and an antiquark and investigate the interaction of light meson with quarks in heavy meson.Through the heavy meson form factor describing the heavy meson structure, we have obtained the interaction kernel between two heavy vector mesons (V M and V M ) derived from one light meson (σ, ω, ρ, φ) exchange in Refs.[8,22].
In our previous works [8,19], BS equation ( 13) has been solved and the mass M 0 and BS wave function χ j λτ (P, p) for the bound state composed of two heavy vector mesons have been obtained.In this paper, we are interested only in mass correction for molecular state and do not repeat the procedure for solving BS equation.
Taking into account the internal structure of heavy mesons (V M, V M and HM) and retaining the lowest order value of T [9], we can obtain the T -matrix element with respect to arbitrary value of the final state energy ǫ ′ in the momentum representation where p 1 , p 3 , p 4 , p 2 are the momenta of four quarks; p ′ 1 and p ′ 2 are the momenta of two heavy vector mesons; q, k and k ′ are the relative momenta between quark and antiquark in heavy mesons, respectively; ε(p) is the polarization vector of vector meson with momentum p, Γ H (K, k) represents BS amplitude of heavy meson, S F (p) is the quark propagator and its superscript is a flavor label, shown as Figure 2. In our approach, the initial bound state is considered as a four-quark state [9], so the generalized BS amplitude of initial bound state should be Γ V λ (p ′ 1 , k)χ j λτ (P, p)Γ V τ (p ′ 2 , k ′ ), which has been specified.In Figure 2(a), the energy in T is equal to the energy of final state, and then the quark momenta in left-hand side of crosses depend on the final state energy and the momenta in right-hand side depend on the initial state energy, i.e., p 1 − p 2 − p 3 + p 4 = Q + Q ′ = P ǫ ′ and p ′ 1 − p ′ 2 = P .In the rest frame, we have P = (0, 0, 0, iM 0 ), P ǫ ′ = (0, 0, 0, iǫ ′ ) and ǫ M < ǫ ′ < ∞.When ǫ ′ = M 0 , the crosses in Figure 2 In general, the decay width of resonance is very small compared with its energy level, i.e., (2π) 3 I(M 0 ) ≪ M 0 .This situation is ordinarily interpreted as implying that both (2π) 3 |D(ǫ)| and (2π) 3 I(ǫ) are also very small quantities, as compared to M 0 .Finally, we can expect that G aa (ǫ) has a pole on the second Riemann sheet from Eq. ( 5) where ∆M = (2π) 3 D(M 0 ) is the correction for energy level of resonance and FIG.2: Matrix element between bound states in the momentum representation.The momenta in the final state satisfy Q + Q ′ = P ǫ ′ and the momentum of the initial state is P .The solid lines denote quark propagators.In the rest frame, we have P = (0, 0, 0, iM 0 ), P ǫ ′ = (0, 0, 0, iǫ ′ ) and  11) and obtain the decay width Γ for physical resonance.
Up to now, a theoretical approach to investigate resonance in the framework of relativistic quantum field theory has been established, which is a new way based on BS theory.
In this paper, we only explore exotic meson resonance which is considered as an unstable meson-meson molecular state.The extension of our approach to more general resonances is straightforward, while the interaction Lagrangian may be modified.
In experiments two strong decay modes of χ c0 (3915) have been observed: J/ψω and D + D − .
Here, we assume that the isoscalar χ c0 (3915) is a mixed state of two unstable molecular states ; and decay channel J/ψω can be exhibited by these two Feynman diagrams.Here, we consider that the light vector meson ω is an elementary particle and the heavy vector meson and ω, respectively.Applying Eq. ( 9), we obtain the function where ) is the bound state matrix element with respect to ǫ ′ .These heavy mesons D * and D * are considered as quark-antiquark bound states, and T (c ′ 1 ;b)a (ǫ ′ ) has been given by Eq. ( 14), where flavor labels C = D and A = B represent c-quark and light quark, respectively.
The meson-quark coupling constant g ′ 0 becomes g ω and g 2 ω = 2.42/2 was obtained within QCD sum rules approach [26].BS amplitudes of heavy vector mesons J/ψ and D * have , where ω J/ψ =0.826GeV [9,27] and ω D * =1.50GeV [14].These momenta in Figure 2(a) become , Q+Q ′ = P ǫ ′ = (0, 0, 0, iǫ ′ ) and P = (0, 0, 0, iM 0 ).Using Eq. ( 14), we can calculate the T -matrix element T (c ′ 1 ;b)a (ǫ ′ ) with respect to arbitrary energy ǫ ′ for channel c ′ 1 .From Eq. ( 16), we obtain the function I 1 (ǫ ′ ) for channel c ′ 1 and dispersion relation (8) becomes where ) and the function I 1 (ǫ ′ ) for channel c ′ 1 are calculated over the real interval ǫ c ′ 1 ,M < ǫ ′ < ∞, and we obtain the mass correction ∆M 1 = (2π) 3 D 1 (M 0 ) due to channel c ′ 1 .For decay channel D + D − , we consider that these heavy pseudoscalar mesons D + and D − are quark-antiquark bound states.Q 1 and Q 2 represent the momenta of final particles, 9), we obtain the function I 2 (ǫ ′ ) where T (c ′ 2 ;b)a (ǫ ′ ) represents the bound state matrix element with ǫ ′ .Considering the lowest order term of T, we obtain T (c ′ 2 ;b)a (ǫ ′ ) represented graphically by Figure 3, where 1 − p ′ 2 = P , and the crosses mean that the momenta of quark propagators and the momentum w of the exchanged light meson depend on Q 1 and Q 2 .The T -matrix element with respect to ǫ ′ for channel c ′ 2 becomes where q, q ′ , k and k ′ are the relative momenta between quark and antiquark in heavy mesons, respectively; the meson-quark coupling constants g were obtained within QCD sum rules approach [8,22,26], O dl (p, q) represents the meson-quark vertex, ∆ F (w) is the light meson propagator, the superscript of quark propagator S F (p) is flavor label, and l = u, d represents P ǫ ′ , the momentum of the initial state is P , and we have P = (0, 0, 0, iM 0 ), P ǫ ′ = (0, 0, 0, iǫ ′ ) and , where ω D =1.50GeV [14].The meson-quark vertex O dl (p, q) is unit matrix for one-σ exchange; and it becomes γ µ for one light vector meson exchange.From Eqs. ( 18), ( 8) and (15), we obtain the mass where ǫ c ′ 2 ,M = M D + + M D − .The T -matrix element T (c ′ 2 ;b)a (ǫ ′ ) and the function I 2 (ǫ ′ ) for channel c ′ 2 are calculated over the real interval ǫ c ′ 2 ,M < ǫ ′ < ∞.Following the same procedure as for channels c ′ 1 and c ′ 2 , we can calculate the mass corrections ∆M 3 and ∆M 4 due to the channel c ′ 3 and closed channel c ′ 4 , respectively.Considering the isospin conservation, we employ the constituent quark masses m u = m d = 0.33GeV, the heavy quark mass m c = 1.55GeV [6] and the meson masses [28].By doing the numerical calculation, we obtain the mass corrections ∆M i (i = 1, 2, 3, 4) due to three open decay channels J/ψω, D + D − , D 0 D0 and one closed channel D * D * , respectively.Subsequently, the mass M for physical resonance χ c0 (3915) can be applied to calculate its decay width.Replacing M 0 by M in Eq. ( 14) and setting ǫ ′ = M, we calculate the matrix element T (c ′ 1 ;b)a (M) and obtain that the width for physical decay model χ c0 (3915) → J/ψω is Γ 1 = 2(2π) 3 I 1 (M).Replacing M 0 by M in Eq. ( 19) and setting ǫ ′ = M, we calculate the matrix element T (c ′ 2 ;b)a (M) and obtain that the width for physical decay model χ For the isospin conservation, it is easy to obtain the width Γ 3 for physical decay model χ c0 (3915) → D 0 D0 .Our numerical results are presented in Table I, and the mass M and full width Γ are in good agreement with experimental data.Furthermore, the calculated D + D − width Γ 2 is very small compared with the calculated J/ψω width Γ 1 , and then we can explain why the decay model χ c0 (3915) → D + D − had not been observed in experiments for a long time.Therefore, this work provides a further verification for the molecular hypothesis of χ c0 (3915) and predicts the exact values of these strong decay widths Γ 1 (χ c0 (3915) → J/ψω), Γ 2 (χ c0 (3915) → D + D − ) and Γ 3 (χ c0 (3915) → D 0 D0 ).
In the actual calculation, we require the meson-quark coupling constants g and the parameters ω H in BS amplitudes of heavy mesons to calculate the mass and decay width of physical resonance.The meson-quark coupling constants can be determined by QCD sum rules approach [26], and these parameters in BS amplitudes of heavy mesons are fixed by providing fits to observables [13,14,27].Our approach also involves the constituent quark masses m u , m d , and the heavy quark mass m c .According to the spontaneous breaking of chiral symmetry, the light quarks (u, d, s) obtain their constituent masses because the vacuum condensate is not equal to zero, and the heavy quark mass m c is irrelevant to vacuum condensate.Normally, the value slightly greater than a third of nucleon mass is employed as the constituent mass of light quark.The value of heavy quark mass m c can be determined by the experimental mass of charmonium system J/ψ.Therefore, there is not an adjustable parameter in our approach.Of course, the values of these parameters, including g, ω H , m u , m d and m c , are values in respective ranges.Simultaneously varying these parameters in respective ranges, we find that the of numerical results are at most 5%.
Despite the large uncertainty of meson mass M σ , it has been found that the uncertainties of numerical results from meson mass M σ are also very small in our previous works [8][9][10]19] and Refs.[11,29].Thus in our approach the calculated mass and decay width are uniquely determined.In this paper we emphatically illuminate the physical meaning of new resonance theory in quantum field theory, and the details in computational process will be shown in our future article.

IV. CONCLUSION
We recognize that resonance can not be completely treated as a stationary bound state and provide a reasonable and feasible scheme to describe unstable system in the framework of relativistic quantum field theory.Based on BS wave function, we provide a description of the prepared state and investigate the time evolution of two-body bound state as determined by the total Hamiltonian.According to dispersion relation, the total matrix elements for all decay channels should be calculated with respect to arbitrary value of the final state energy, and these matrix elements are expressed in terms of the Heisenberg picture.Mandelstam's approach is generalized to calculate the matrix element between bound states with arbitrary value of the final state energy, which is exhibited in extended Feynman diagram.Finally, the mass and decay width for physical resonance are obtained.In this paper, we illustrate this new resonance theory in quantum field theory by reference to the example of exotic meson which is considered as an unstable meson-meson molecular state, and obviously our work can be extended to more general resonances and creates a new paradigm for investigating hadron resonances.

FIG. 1 :
FIG.1: General matrix element between two bound states Q out|T φ(z 1 )...φ(z i )...φ(z n ′ )|P in with respect to ǫ ′ .The filled blob represents the two-particle irreducible Green's function, and the unfilled ellipses represent BS amplitudes.The crosses mean that the final state energy extends from ǫ M to +∞ while the initial state energy is specified, which is different from the traditional Feynman diagram.
(a) disappear and Figure 2(a) becomes the traditional Feynman diagram, shown as Figure 2(b).From Figure 2(b), we have calculated the matrix element T (c ′ ;b)a (M 0 ) and the decay width Γ(M 0 ) with the mass of meson-meson bound state [9].
the final state energy extends from ǫ M to +∞ while the initial state energy is specified, and the crosses mean that the momenta of quark propagators depend on the final state energy; in diagram (b) the crosses disappear when ǫ ′ = M 0 .(2π) 3 D(M 0 ) is the physical mass for resonance.The mass M 0 of two-body bound state is obtained by solving homogeneous BS equation (2), which should not be the mass of physical resonance.Γ(M 0 ) with mass M 0 also should not be the width of physical resonance, which should depend on its physical mass M. Using the mass M 0 and BS wave function of bound state, we can calculate the imaginary part I(M 0 ) and obtain the correction ∆M = (2π) 3 D(M 0 ) due to decay channels.Replacing M 0 in the momentum of initial bound state by M and setting ǫ ′ = M, we can calculate the matrix element T (c ′ ;b)a (M) from Eq. (

D * 0
D * 0 and D * + D * − with spin-parity quantum numbers 0 + .Firstly, we consider the mixed state of two bound states D * 0 D * 0 and D * + D * − , which can be denoted as 1/ √ 2|D * 0 D * 0 + 1/ √ 2|D * + D * − .In Refs.[8, 10, 19], we have obtained the mass M 0 and BS wave function χ 0 + λτ (P, p) for this mixed state of two bound states D * 0 D * 0 and D * + D * − without an adjustable parameter.In this paper, our attention is focused on the mass correction due to all decay channels and the decay width of physical resonance.Let us list all open and closed channels.The narrow state χ c0 (3915) was discovered in 2005 [23] and for a long time a series of experiments only observed one strong decay mode of χ c0 (3915): J/ψω denoted as c ′ 1 .In 2020 LHCb Collaboration observed another decay channel D + D − [25] denoted as c ′ 2 .Though the neutral channel D 0 D0 still has not been observed, this neutral channel should exist for the isospin conservation, which is denoted as c ′ 3 .Because the total energy ǫ ′ of the final state extends from ǫ M to +∞, we obtain one closed channel D * D * derived from the interaction Lagrangian (1), denoted as c ′ 4 .Since bound state lies below the threshold, i.e., M 0 < M D * + M D * , the closed channel c ′ 4 can not occur inside the physical world.Then we can apply Eqs.(11) and (12) to evaluate the T -matrix element T (c ′ ;b)a (ǫ ′ ) for arbitrary decay channel.In Figure 2, V M and V M become D * and D * , respectively; HM becomes J/ψ with momentum Q

FIG. 3 :
FIG. 3: Matrix element for decay channel D + D − .The momenta in the final state satisfy Q 1 +Q 2 = w represents the momentum of the exchanged light meson.The crosses mean that the momenta of quark propagators and the momentum w of the exchanged light meson depend on the final state energy.the u, d-antiquark in heavy vector meson D * 0 or D * + , respectively.BS amplitude of heavy pseudoscalar meson D + has the form Γ

TABLE I :
Mass M and width Γ for physical resonance χ c0 (3915).M 0 is the mass of mixed state of two bound states D * 0 D * 0 and D * + D * − , ∆M i is the calculated correction due to ith decay channel.