arxiv: v1 [hep-ph] 2 Oct 2008

Transcription

1 APS/1-QED Meson decay in a corrected P 0 model D. T. da Silva, M. L. L. da Silva, and J. N. de Quadros Instituto de Física, Universidade Federal do Rio Grande do Sul Av. Bento Gonçalves, 9500, Porto Alegre, Rio Grande do Sul, CEP , Brazil D. Hadjimichef Instituto de Física e Matemática, Universidade Federal de Pelotas Campus Universitário, Pelotas, Rio Grande do Sul, CEP , Brazil Dated: September 0, 018) arxiv: v1 hep-ph Oct 008 Extensively applied to both light and heavy meson decay and standing as one of the most successful strong decay models is the P 0 model, in which q q pair production is the dominant mechanism. The pair production can be obtained from the non-relativistic limit of a microscopic interaction Hamiltonian involving Dirac quark fields. The evaluation of the decay amplitude can be performed by a diagrammatic technique for drawing quark lines. In this paper we use an alternative approach which consists in a mapping technique, the Fock-Tani formalism, in order to obtain an effective Hamiltonian starting from same microscopic interaction. An additional effect is manifest in this formalism associated to the extended nature of mesons: bound-state corrections. A corrected P 0 is obtained and applied, as an example, to and a 1 ρ decays. PACS numbers: Tk, 1.9.Jh, 1.5.-k I. INTRODUCTION A great variety of quark-based models are known that describe reasonable success single-hadron properties. A natural question that arises is to what extent a model which gives a good description of hadron properties is, at the same time, able to describe the complex hadron-hadron interaction or by the same principles hadron decay. In particular, the theoretical aspects of strong decay have been challenged by QCD exoticaglueballs and hybrids) a consistent understanding of the mixing schemes for these states is still an open question 1-. The nature of the family of new mesons X,Y,Z 4 is another unsolved puzzle: are they actually new q q mesons, hadronic molecules or something else? In the direction of clarifying these questions is the successful decay model, the P 0 model, which considers only OZIallowed strong-interaction decays. This model was introduced over thirty years ago by Micu 5 and applied to meson decays in the 1970 by LeYaouanc et al. 6. This description is a natural consequence of the constituent quark model scenario of hadronic states. T. Barnes et al have made an extensive survey of meson states in the light of the P 0 model. Two basic parameters of their formulation are the interaction strength) and the wave function s extension parameter). Although they found the optimum values near = 0.5 and = 0.4 GeV, for light 1S and 1P decays, these values lead to overestimates of the widths of higher-l states. In this perspective a modified q q paircreation interaction, = 0.4 was preferred. In the present work, we employ a mapping technique Electronic address: dimiter@ufpel.edu.br; dimihadj@gmail.com in order to obtain an effective interaction for meson decay. A particular mapping technique long used in atomic physics 11, the Fock-Tani formalism FTf), has been adapted, in previous publications 1-16, in order to describe hadron-hadron scattering interactions constituent interchange. Now this technique has been extended in order to include meson decay. We start from the microscopic q q pair-creation interaction, as will be shown, in lower order, the P 0 results are reproduced. An additional and interesting feature appears in higher orders of the formalism: corrections due to the boundstate nature of the mesons and a natural modification in the q q interaction strength. In the Fock-Tani formalism one starts the Fock representation of the system using field operators of elementary constituents which satisfy canonical anti) commutation relations. Composite-particle field operators are linear combinations of the elementary-particle operators and do not generally satisfy canonical anti) commutation relations. Ideal field operators acting on an enlarged Fock space are then introduced in close correspondence the composite ones. Next, a given unitary transformation, which transforms the single composite states into single ideal states, is introduced. Application of the unitary operator on the microscopic Hamiltonian, or on other hermitian operators expressed in terms of the elementary constituent field operators, gives equivalent operators which contain the ideal field operators. The effective Hamiltonian in the new representation has a clear physical interpretation in terms of the processes it describes. Since all field operators in the new representation satisfy canonical anti)commutation relations, the standard methods of quantum field theory can then be readily applied. In this paper we shall extend the FTf to meson decay processes. In the next section we review the basic

2 aspects of the formalism. Section III is dedicated to obtain an effective decay Hamiltonian. In section IV, two light mesons decays examples are calculated and a 1 ρ. The summary and conclusions are followed by appendixes which detail the method employed throughout this work. II. MAPPING OF MESONS This section reviews the formal aspects of the mapping procedure and how it is implemented to quark-antiquark meson states 1. The starting point of the Fock-Tani formalism is the definition of single composite bound states. We write a single-mesonstate in terms of a meson creation operator M α as α = M α 0, 1) 0 is the vacuum state. The meson creation operator M α is written in terms of constituent quark and antiquark creation operators q and q, M α = Φ µν α q µ q ν, ) Φ µν α is the meson wave function and q µ 0 = q ν 0 = 0. The index α identifies the meson quantum numbers of space, spin and isospin. The indices µ and ν denote the spatial, spin, flavor, and color quantum numbers of the constituent quarks. A sum over repeated indices is implied. It is convenient to work orthonormalized amplitudes, α = Φ µν α Φµν = δ α. ) The quark and antiquark operators satisfy canonical anticommutation relations, {q µ,q ν } = { q µ, q ν } = δ µν, {q µ,q ν } = { q µ, q ν } = {q µ, q ν } = {q µ, q ν} = 0. 4) Using these quark anticommutation relations, and the normalization condition of Eq. ), it is easily shown that the meson operators satisfy the following non-canonical commutation relations M α,m = δ α α, M α,m = 0, 5) In addition, α = Φ µν α Φµσ q σ q ν +Φ µν α Φρν q ρ q µ. 6) q µ,m α = Φ µν α q ν, q ν,m α = Φ µν α q µ, q µ,m α = q ν,m α = 0. 7) The presence of the operator α in Eq. 5) is due to the composite nature of the mesons. This term enormously complicates the mathematical description of processes that involve the hadron and quark degrees of freedom. The usual field theoretic techniques used in manybody problems, such as the Green s functions method, Wick s theorem, etc, apply to creation and annihilation operators that satisfy canonical relations. Similarly, the non-vanishing of the commutators q µ,m α and q ν,m α is a manifestation of the lack of kinematic independence of the meson operator from the quark and antiquark operators. Therefore, the meson operators M α and M α are not convenient dynamical variables to be used. A transformation is defined such that a single-meson state α is redescribed by an ideal ) elementary-meson state by α U 1 α = m α 0, 8) m α an ideal meson creation operator. The ideal mesonoperatorsm α andm α satisfy, bydefinition, canonical commutation relations m α,m = δ α, m α,m = 0. 9) The state 0 is the vacuum of both q and m degrees of freedom in the new representation. In addition, in the new representation the quark and antiquark operators q, q, q and q are kinematically independent of the m α and m α q µ,m α = q µ,m α = q µ,m α = q µ,m α = 0. 10) The unitary operator U of the transformation is Ut) = exptf, 11) F is the generator of the transformation and t a parameter which is set to / to implement the mapping. The next step is to obtain the transformed operators in the new representation. The basic operators of the model are expressed in terms of the quark operators. Therefore, in order to obtain the operators in the new representation, one writes qt) = U 1 qu, qt) = U 1 qu. 1) The generator F of the transformation is F = m α M α M αm α 1) M α = n i=0 M i) α, 14) M α, M = δ α + OΦ n+1 ), M α, M = M α, M = 0. 15) It is easy to see from 1) that F = F which ensures that U is unitary. The index i in 14) represents the order of the expansion in powers of the wave function Φ. The M α operator is determined up to a specific order n consistent 15). The examples studied in 1

3 requiredthedeterminationof below M 0) α = M α ; M1) α = 0 M ) α M i) α uptoorderasshown = 1 αm ; M) α = 1 M T α M, 16) T α = M α,. In the zero-order approximation, overlap among mesons is neglected and terms of the same power in the bound-state wave function Φ α Φ α) are collected. In order to have a consistent power counting scheme, the implicit Φ α Φ α ) entering via Eq. ) are not counted. The consequence of this is that the equations for m α and M α are manifestly symmetric, dm α t) dt d M α t) dt = m α t),f = M α t), = M α t),f = m α t), 17) and their solutions involve only trigonometric functions of t, m α t) = M α sint+m α cost, M α t) = M α cost m α sint. 18) The equations of motion for the quark operators q and q can be obtained by making use of Eq. 7) in a similar way, dq µ t) dt = q µ t),f ; d q µ t) dt = q µ t),f.19) In the zero-order approximation, the effects of the meson structure are neglected resulting In first order one has q 0) µ t) = q µ, q 0) ν t) = q ν, m 0) α t) = m α cost+m α sint, M 0) α t) = M α cost m α sint. 0) q 1) µ t) = Φ µν1 α q ν 1 m α sint+m α 1 cost), q ν 1) t) = Φµ1ν α q µ 1 m α sint+m α 1 cost), m 1) α t) = 0, M1) α t) = 0. 1) The second and third order solutions to 19) were calculated in reference1 and appear again, for completeness, in appendix A, together the higher order operators required in our calculation. Once a microscopic interaction Hamiltonian H is defined, at the quark level, a new transformed Hamiltonian can be obtained. This effective interaction we shall call the Fock-Tani Hamiltonian and is evaluated by the application of the unitary operator U on the microscopic Hamiltonian H FT = U 1 HU. The transformed Hamiltonian H FT describes all possible processes involving mesons and quarks. The general structure of H FT is of the following form H FT = H q +H m +H mq, ) the first term involves only quark operators, the secondoneinvolvesonly idealmeson operators,and H mq involves quark and meson operators. In H FT there are higher order terms that provide bound-state corrections also called orthogonality corrections) to the lower order ones. The basic quantity for these corrections is the bound-state kernel ρτ; λν) defined as ρτ;λν) = Φ ρτ α Φ λν α. ) To discuss the physical meaning of the bound-state corrections and how they modify the fundamental quark interaction we shall present an example, in a toy model similar to the model studied in 1, the basic arguments are outlined. In this example, the starting point is a two-body microscopic quark-antiquark Hamiltonian of the form H q = T µ)q µ q µ +T ν) q ν q ν +V q q µν;σρ)q µ q ν q ρq σ + 1 V qqµν;σρ)q µ q ν q ρq σ + 1 V q qµν;σρ) q µ q ν q ρ q σ. 4) The transformation H FT = U 1 H q U is implemented again by transforming each quark and antiquark operator in Eq. 4), a similar structure to Eq. ) is obtained. In free space, the wave function Φ of Eq. ) satisfy the following equation Hµν;σρ)Φ σρ α = ǫ α Φ µν α, 5) Hµν; σρ) is the Hamiltonian matrix Hµν;σρ) = δ µσ δ νρ Tσ)+Tρ) +V q q µν;σρ), 6) ǫ α is the total energy of the meson. There is no sum over repeated indices inside square brackets. The effective quark Hamiltonian H q has an identical structure to the microscopic quark Hamiltonian, Eq. 4), except that the term corresponding to the quark-antiquark interaction is modified as follows V q q = V q q H H + H, 7) V q q V q q µν;σρ) and the contraction H Hµν;τξ) τξ;σρ). An important property of the bound-state kernel is µν;σρ)φ σρ α = Φ µν α, 8) which follows from the wave function s orthonormalization, Eq. ). In the case that Φ is a solution of Eq. 5), the new quark-antiquark interaction term becomes V q q µν;σρ) = V q q µν;σρ) α ǫ α Φ µν α Φσρ α. 9)

4 4 The spectrum of the modified quark Hamiltonian, H q, is positive semi-definite and hence has no bound-states11. This result is exactly the same as in Weinberg s quasiparticle method 17, the bound-states are redescribed by ideal particles. The new V q q is a weaker potential, modified in such a way that no new bound-states are formed. In the quark-meson sector of Eq. ) in H mq appears a term related to spontaneous meson break-up H m q q = Vµν;α)q µ q ν m α 0) Vµν;α) = Hµν;σρ)Φ σρ α µν;σρ)hσρ;τλ)φ τλ α. 1) Again, in the case that Φ is a solution of Eq. 5), a straightforward calculation demonstrates that H m q q = 0. When there is no external interaction, this result is a direct consequence of the bound-state s stability against spontaneous break-up. This term can be interesting in studies related to dense hadronic mediums. For these systems the wave function is, in general, not a solution of Eq. 5) and the strength of the potential Vµν;α) is now only decreased 1. In the ideal meson sector H m many similar approaches to FTf 1 have obtained the meson-meson scattering interaction in the Born approximation: Resonating Group Method RGM) 18, Quark Born Diagram Formalism QBDF) 19, H mm = T mm +V mm, ) T mm is the kinetic term and V mm is the mesonmeson interaction potential constituent interchange. This potential is given by V mm = V dir mm +Vexc mm +Vint mm, ) Vmm dir is the direct potentialno quark interchange), Vmm exc int the quark exchange term and Vmm the intraexchange term. As shown in Ref. 1 and 1, if one extends the FT calculation to higher orders a new mesonmeson Hamiltonian is obtained H mm = H mm +δh mm 4) δh mm is the bound-state correction Hamiltonian. If the wavefunction Φ is chosen to be an eigenstate ofthe microscopic quark Hamiltonian, then the intra-exchange term Vmm int is cancelled V int mm +δh mm = 0. 5) In summary, these examples reveal an important and common feature of bound-state corrections: they weaken the quark-antiquark potential. In the next section we shall follow the same procedure for a quark pair creation interaction, which is fundamental for the description of meson decay. Similar to the toy model, the resulting interaction that describes meson decay, will contain a Born order contribution and a bound-state correction. III. THE P 0 DECAY MODEL IN THE FOCK-TANI FORMALISM Inthe paperofe.s.ackleh, T.BarnesandE.S.Swanson 7 a formulation of the P 0 model is presented. It regards the decay of an initial state meson in the presence of a q q pair created from the vacuum. The pair production is obtained from the non-relativistic limit of the interaction Hamiltonian H I involving Dirac quark fields H I = m q d x ψ x)ψ x), 6) is the pair production strength. For a q q meson A to decay to mesons B + C we must have q q) A q q) B + q q) C. To determine the decay rate a matrix element of 6) is evaluated BC H I A = δ P A P B P C )h fi. 7) Theevaluationofh fi isperformedbydiagrammatictechnique for drawing quark lines. The h fi decay amplitude is combined relativistic phase space, resulting in the differential decay rate dγ A BC dω = P E B E C M A h fi 8) which after integration in the solid angle Ω a usual choice for the meson momenta is made: PA = 0 P = P B = P C ). In our approach, the starting point for the Fock-Tani h fi is also the microscopic Hamiltonian H I in 6). The momentum expansion of the quark fields, color and flavor are not represented explicitly, is ψ x) = s d k ) /u s k)q s k) +v s k) q s k)e i k x. 9) In the product ψx)ψx) we shall retain only the q q term, which yields from Eq. 6) a Hamiltonian in a compact form, H I = V µν q µ q ν 40) sum integration) is again implied over repeated indexes. In the compact notation, the quark and antiquark momentum, spin, flavor and color are written as µ = p µ,s µ,f µ,c µ ); ν = p ν,s ν,f ν,c ν ), while the pair creation potential V µν is given by V µν m q δ p µ + p ν )ū sµf µc µ p µ )v sνf νc ν p ν ). 41) It should be noted that since Eq. 6) is meant to be taken in the nonrelativistic limit, Eq. 41) should be as well. In the meson decay calculations, of the next section, this limit is considered.

5 5 Applying the Fock-Tani transformation to H I one obtains the effective Hamiltonian H FT = U 1 H I U. 4) The physical quantities in the FTf appear in a second quantization notation. The effective decay amplitude will be a product of the ideal meson operators the following structure in the ideal meson sector: m m m. To obtain this product corresponds to expand, in powers of the wave function, up to third order. A Hamiltonian that describes this decayprocess, which we shall callh m, can be extracted from the mapping 4) by the following products H m = V µν q ) µ q 0) ν +V µν q 1) µ q ) ν. 4) After the substitution of Eqs. 0), 1), A1) and A) into4) results in the effective meson decay Hamiltonian H m = f µν α,,)v µν m α m m 44) f µν α,,) = Φ µτ α Φ ρν Φρτ. 45) In the ideal meson space the new initial and final states involve only ideal meson operators A = m 0 and BC = m αm 0. The P 0 amplitude is obtained in the FTf by an expression equivalent to Eq. 7),! BC H I A = 0 m α m H m m 0 = f µν α,,)v µν +f µν,α,)v µν. 46) The term f µν,α,) of 46) is shown in Fig. 1a), the term f µν α,,) corresponds to the same diagram α. In the FTf perspective a new aspect is introduced to meson decay: bound-state corrections. The lowest order correction is one that involves only one bound-state kernel µν; σρ). This implies that the Hamiltonian representing this correction must be of fifth order in the power expansion of the wave function. We shall call this new Hamiltonian, the same basic operatorial structure m α m m, of δh m. The only combinationsq i) q j) thatresultsinafifthorderhamiltonian are δh m = q µ ) q ν ) +q µ 1) q ν 4) +q µ 5) q ν 0) V µν. 47) Details of this calculation is found in the appendix B. The bound-state corrected P 0 Hamiltonian, which shall be called the C P 0 Hamiltonian, is H CP0 = H m +δh m = Φ ρξ α Φ λτ Φ σ V CP0 m αm m 48) V CP0 is a condensed notation for V CP0 = δ µλ δ ξν δ ρ δ στ 1 δ σξ δ λ ρτ;µν) δ σξ δ µλ ρτ;ν) δ ξν δ λ ρτ;µσ) V µν. 49) IV. LIGHT MESON DECAY EXAMPLES The light meson sector is an interesting test ground the effects of the bound-state correction can be compared to the usual P 0 model. In particular, as examples, two specific decay processes will be studied: and a 1 ρ. The wave function and details of the matrix elements are found in the appendix C. The general decay amplitude can be written as h CP0 fi = 1/4 MCP0 1/ fi. 50) For the first decay process,, results in a decay amplitude given by C / M CP0 fi = C 01 Y 00 Ω x )+C 1 Y 0 Ω x ), 51) C 1 x 11/ x = P/ and ) e 1 x) = exp x e x 1 x) / 9/ e 1x) 17/ 7 7/ 9 e x) ; e x) = exp 1 x e x) 5) ) 9x.5) 8 The decay rate has a straightforward evaluation, by substituting 51), 5) in 50) and then in 8) obtaining Γ b1 = x E E C ) 01 M +C 1. 54) b1 The second decay process, a 1 ρ, is similar to the former one and results in the following amplitude C 01 9/ 5/ M CP0 fi = C 01 Y 00 Ω x )+C 1 Y 0 Ω x ), 55) C 1 x x e 1 x) 11/ / e 1x) 7 5 9/ 7 e x) 7/ and by a similar procedure one obtains 1 x e x) 56) Γ a1 ρ = x E ρe C01 +C ) 1. 57) M a1

6 6 FIG. 1: Diagrams representing the C P 0 model. Diagram a) corresponds to the P 0 amplitude. Diagrams b), c) and d) are the bound state corrections. The complete h fi amplitude includes the diagrams above plus diagrams α. In the former equations, e x) = 0, recovers the original P 0 results. In addition to the decay widths Γ, and a 1 mesons have D/S ratios, which give a sensitive test for decay models. By definition, these quantities are obtained from the ratios of C 1 and C 01 coefficients, in equations 5) and 56). D S D S a1 ρ b1 = = x 1/ e 1 x) 1/ 5/ 5e 7 7/ x) 1 9 x e 1 x) / / 1 x e x) x / e 1 x) 9/ 1/ e 7 7/ x) 1 9 x e 1 x) / / 1 x e x). 58) The meson masses assumed in the numerical calculation were M = 18 MeV; M ρ = 775 MeV; M a1 = 10 MeV; M b1 = 19 MeV; M = 78 MeV 0. The choice of SHO wave functions allow exact evaluations of the decay amplitudes even in the corrected model. A first new aspect that appears is the presence of a new dependence in the exponential of the corrected term. This implies in a different range for the boundstate correction due to the fact that e x)/e 1 x) 0 as x. The correction introduces the bound-state kernel, Eq. ), to the calculation of the decay processes. A general sum over the meson index α is present and as stated before, this index represents the quantum numbers of space, spin and isospin. The OZI-allowed decays represent, flavor conserved continuous anti)quark lines. A direct consequence of this fact is the possibility to sum over a larger set of mesons in the α index. In our calculation the sum was restricted only to the final state mesons. In the b decay, there are two boundstate kernel contributions one associated to meson and the other to +. Similarly, the a + 1 ρ+ 0 decay has two bound-state kernel contributions one associated to ρ + meson and the other to the 0. In this example, the parameters were chosen in order to give a closer fit to the experimental data. In the decay, width and partial waves are known accuracy. The P 0 model s optimum fit for the data Γ and D/S ratio) is achieved = and = 0.97 GeV. In the C P 0 model a similar fit is obtained = 0.59 and = 0.96 GeV. These parameters are used in the two models to describe the a 1 decay. The results for Γ as a function of appear in figure and specific values are presented in table I. In figure, the D/S ratios for the two models are plotted. V. SUMMARY AND CONCLUSIONS In this paper we have presented an alternative approach for meson decay which consists in a mapping

7 7 TABLE I: Decay rates P 0 = e = 0.97 GeV ) and C P 0 = 0.59 e = 0.96 GeV ) Γ MeV) D/S Decay Exp 0 P 0 C P 0 Exp 0 P 0 C P ) a 1 ρ 50 to ) a 1 ρ Exp > Γ GeV) P 0 C P 0 D/S a 1 Exp a 1 -> ρ P0 C P GeV) GeV) FIG. : Decay rates for and a 1 ρ decays, for P 0 = 0.506) and C P 0 = 0.59) models. FIG. : D/S ratios in and a 1 ρ decays, for P 0 and C P 0 models. Acknowledgments technique, known as the Fock-Tani formalism, long used in atomic physics. This formalism has been applied to hadron-hadron scattering interactions constituent interchange. The challenge, resided in extending the approach to include meson decay. After demonstrating that in lower order the result obtained was equivalent to the P 0 model, an additional feature pointed out was the appearance of bound-state corrections in the effective decay Hamiltonian. These corrections present a natural modification in the q q interaction strength. As an example, we studied two decay processes and a 1 ρ. The D/S ratios, in Fig. ), show that a common range of values for mesons is obtained. In a new calculation the inclusion of other decay processes it might require different values 1. The corrected model presents an interesting feature that for these two mesons the decay width differs slightly when compared the P 0, but D/S ratios are improved. The examples studied here are encouraging but a more extensive survey of the light meson sector would be a necessary next step. The inclusion of the full meson octet, in the evaluation of the bound-state correction, may provide a fine tuning for the model. The authors would like to thank H. Stöcker, J. Aichelin and W. Greiner for important and enlightening discussions. This research was supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico CNPq), Universidade Federal do Rio Grande do Sul UFRGS) and Universidade Federal de Pelotas UFPel). APPENDIX A: SECOND AND THIRD ORDER OPERATORS The second order operators q µ ) t) = 1 Φ µν1 α Φ µν1 M α q µ q ν ) t) = 1 Φ µ1ν α Φ µ1ν M α q ν, A1) M α = m αm sintcost m αm sin t M α M cost sin t) M α m sint sintcost). A)

8 8 The third order operators are q µ ) t) = 1 Φ ρσ α Φ µσ Φρσ1 q σ 1 M α 1 Φ ρσ α Φ µν1 α Φ ρσ1 q ν 1 q σ 1 q σ M 1 Φ ρσ α Φ µν1 α Φ ρ1σ q ν 1 q ρ 1 q ρ M δf µν α,,) = 1 8 Φ ρσ α Φ µτ ρτ;λν)φ λσ. B6) The q µ 5) q ν 0) combination implies in a fifth order generator to obtain the complete set of equations of motion 17) and 19) q ν ) t) = 1 Φ ρσ α Φ ρν Φρ1σ q ρ 1 M α + 1 Φ ρσ α Φ µ1ν α Φ ρσ1 q µ 1 q σ 1 q σ M M 5) α = M Z α M M M W αδ M δm B7) + 1 Φ ρσ α Φµ1ν α Φρ1σ q µ 1 q ρ 1 q ρ M, A) M α = m α m m sin t+m α M m sint sin t ) +M α m M sint sintcost sin t ) + M α m m +m α M m ) cost+cos t ) +m α m M cost+cos t+sin t ) +M α M M cost cos t sin t ) +m αm M sint sintcost sin t ) M = M 1 cost)+m sint. APPENDIX B: THE δh m HAMILTONIAN A4) TheδH m HamiltonianisevaluatedfromEq. 47). The q µ ) q ν ) combination can be obtained from A1) and A) δh 1 = q ) µ q ) ν V µν = δf µν 1 α,,)v µν m α m m B1) δf µν 1 α,,) = 1 4 Φ ρσ α Φ µτ ρτ;λν)φ λσ. B) The q µ 1) q ν 4) combination has an important feature: a contribution from a higher order operator. A new generator M α has to be evaluated, the inclusion of the following fourth order term M 4) α = 8 α M 1 8 M α, δ M δ M 1 4 M M α,t δ M M δ. B) The only relevant term in the q 4) ν for meson decay is q ν 4) t) 1 8 Φ ση α σν;ρτ)φ ρη M 0) α t) q τ M 0) t). The resulting contribution is then B4) δh = q 1) µ q 4) ν V µν = δf µν α,,)v µν m αm m B5) Z α = 8 T αδ δ 5 8 δt αδ 1 4 T δ δα W αδ = M α, Q δ, T αδ Q αδ = 1 α, δ M α,t δ. B8) The only relevant terms in the q 5) µ for meson decay are q 5) µ t) The resulting contribution is 1 ρτ;µ)φ ρσ α Φ λτ Φ λσ 1 4 ρτ;µσ)φ ρ α Φ λτ Φ λσ 8 ρτ;λ)φ µτ α Φ ρσ Φ λσ M 0) α t)m 0) t) q M 0) t). B9) δh = δf µν α,,)v µν m α m m B10) δf µν α,,) = 1 Φ ρσ α Φ λτ ρτ;µν)φ λσ The complete δh m Hamiltonian is 1 4 Φ στ α Φ ρν ρτ;µλ)φσλ 8 Φ ρσ α Φ µτ ρτ;λν)φ λσ.b11) δh m = δh 1 +δh +δh = δf µν α,,)v µν m α m m B1) δf µν α,,) = δf µν 1 +δf µν +δf µν 1 = Φ ρσ α Φ λτ ρτ;µν)φ λσ 1 4 Φ ρσ α Φ µτ ρτ;λν)φ λσ. 1 4 Φ στ α Φ ρν ρτ;µλ)φσλ. B1)

9 9 APPENDIX C: DECAY, WAVE FUNCTION AND MATRIX ELEMENTS. L q q = 1 Φ 1m p) = φ p) Y 1m Ω p ) C8) We will use the decay b + 1 +ẑ) +ẑ)+ to illustrate the nature of our formalism and, simply quote the other case in the text. φ p) = 5 1 p exp p 8. C9) 1. The wave function The general meson wave function can be written as Φ µν α = χ s1s S α f f1f f α C c1c Φ P α p 1 p nl, C1) a direct product of the spin χ s1s S α the indexes s 1 and s are the quark antiquark) spin projections s = 1 and s = ); the index S α denotes the meson spin; flavor f f1f f α ; color C c1c P and space Φ α p 1 p nl components. In all our calculations the color component will be given by C c1c = 1 δ c1c. C) The spatial part is defined as harmonic oscillator wave functions Φ P α p 1 p nl = δ P α p 1 p ) Φ nl p 1, p ) C) Φ nl p i, p j ) is given by Φ nl p i, p j ) = 1 )l N nl p i p j l exp p i p j ) 8 L l+1 pi p j ) n 4 Y lm Ω pi p j ) C4) p ij) theinternalmomentum, thesphericalharmonic Y lm, a scale parameter, N nl the normalizationconstant dependent on the radial and orbital quantum numbers N nl = n!) Γn+l+/) 1. C5) The Laguerre polynomials L l+1 n p) are defined as L l+1 n p) = n k=0 ) k Γn+l+/) n k)! k!γk +l+/) p k. C6) In this paper two kinds of light non-strange mesons will be studied: 1. L q q = 0 ϕ p) Φ 00 p) = 1 exp p /4 / 8 C7) Returning to our example the pion, has J = 0 and J = 1. We choose the +ẑ) direction for this calculation, so the spin wave functions become = 1 ) = = 1 ) or in the χ notation χ 11 = 1, χ 1 = χ 1 = χ = 0 C10) χ 11, = χ, = 0 ; χ 1, = χ 1, = 1. C11) The flavor component f fµfν f α follows the same logic as the spin part b + 1 = + = u d = 1 uū + d d ) C1) The b + 1 and the + mesons have the same flavor contribution f 1 = 1;f 11 = f 1 b + 1,+ b + 1,+ For, one has = b + f 1,+ b + 1 f 11 = f = 1 ; f 1 = f1 = 0.. The spin matrix elements,+ = 0. C1) C14) In the evaluation of a decay amplitude, the following spin matrix element is necessary χ s σ P ) χ c s C15) ) ) ) ) 1 0 χ 1 = ; χ 0 = ; χ 1 c 0 1 = ; χ 1 c 1 =. 0 C16) By direct calculation one can show χ 1 σ P ) χ c 1 = P x ip y χ 1 σ P ) χ c = P z ) χ χ σ P σ P χ c 1 = P z ) χ c = P x +ip y ). C17)

10 10. Matrix elements: b Decay The transition considered is of the form m m α + m, the initial state is A = m 0 and the final state is given by BC = m α m 0. The matrix element of the uncorrected part results in d 1 and d are defined as BC H m A = d 1 d, d 1 Φ ρν α Φ µλ Φ ρλ V µν d Φ µλ α Φ ρν Φρλ V µν. C18) C19) Equations C19) can be decomposed according to the sector of the wave function they correspond: flavor, color, spin-space: d 1 = d f 1 dc 1d s e 1 d = d f dc ds e. C0) The matrix elements of the bound-state correction refer to diagrams b), c) and d) of figure 1). The boundstate kernel s definition as Φ ρτ α Φ λν α implies in an additional element, due to the contraction in the α index, a sum over species requirement 11. A question that naturally arises is, which states to include in this sum? We shall adopt in our calculation a restrictive choice: include in the sum only the particles that are present in the final state. For the b + 1 decay, ρτ;λν) will have two contributions: and +. Similarly, the a + 1 decay shall be corrected by the final state mesons ρ + and 0. Duetotheparityassignmentofthespatialpart, theintegration of diagram 1b) is zero. Spatial symmetry also implies that the matrix elements of diagrams 1c) and 1d) are equal. This simplifies our calculation, reducing the problem to the evaluation of diagram 1c) only. The bound-state correction bsc) matrix element reduces to evaluate the following expression BC δh m A = d bsc 1 d bsc C1) d bsc 1 = 1 4 d 1 +d 1 ) d bsc = 1 4 d +d ) C) d 1j = Φ ρσ α Φ µτ j ρτ;λν)φ λσ V µν d f 1j dc 1j ds e 1j d j = Φ µτ α Φ ρσ j ρτ;λν)φ λσ V µν d f j dc j ds e j. In C) j refers to mesons and +. C) Flavor: Color: 4. b Decay uncorrected) d f 1 = d f = f fρfν f fµf λ f fρf λ δ fµf ν = 1. C4) d c 1 = d c = 1 δcρcν δ cµc λ δ cρc λ δ cµcν = 1. Spin-space: The spin matrix element is d s 1 = ds = χsρsν χ sµs λ C5) χ sρs λ Vs s e µs ν = 1 V 11 s e p µ, p ν ) C6) V s e 11 p µ, p ν ) = δ p µ + p ν )χ 1 σ p µ p ν )χ C 1. C7) Using C17) to evaluate C7) and after integrating momentum conservation deltas one arrives in ) d1 s e = d K K x ik y )ϕ P K φ P K ) ) )Y 11 Ω P ˆ K ϕ P K. C8) Introducing the spatial wave function and integrating ) d1 s e 7/ = ){1 9 5/ 1/4 1/ x Y 00 Ω x ) + } 5 x Y 0 Ω x ) e 1 x). C9) d s e is obtained from d1 s e by P P. The decay amplitude results ) 4 h fi = ){1 9 5/ 1/4 1/ x Y 00 Ω x ) + } 5 x Y 0 Ω x ) e 1 x). C0) 5. b Decay bound-state corrected) The quantities between... in the following expressions are related to the bound-state kernel. Flavor: d f 1 = f fρfσ f fµfτ f f ρf τ f f λf µ f f λ f σ = 1 d f 1 = f fρfσ f fµfτ f f ρf τ f f λf µ f f λ f σ = 1 d f = ffµfτ f fρfσ f f ρf τ f f λf µ f f λ f σ = 1 d f = ffµfτ f fρfσ f f ρf τ f f λf µ f f λ f σ = 0. C1)

11 11 Color: After integration one finds d c 1 = d c 1 = d c = d c = 1 9 δcρcσ δ cµcτ δ cρcτ δ c λc ν δ c λc σ δ cµcν = 1. C) Spin-space: The spin matrix element is d1 s = χ sρsσ χ sµsτ χ sρsτ = 1 V s e 11 p µ, p ν ) d1 s = χ sρsσ χ sµsτ χ sρsτ = 1 4 V s e 11 p µ, p ν ) d s = χsµsτ χ sρsσ = 0 d s = χ sµsτ χ sρsσ χ sρsτ χ sρsτ χ s λs ν χ s λs ν χ s λs ν χ s λs ν χ s λs σ Vs s e µs ν χ s λs σ Vs s e µs ν χ s λs σ Vs s e µs ν χ s λs σ Vs s e µs ν = 1 4 V s e 11 p µ, p ν ). C) Due to symmetries in the spatial part the following relations are true d s e 1 = d s e = 1 ds e 1 ; d1 s e is given by d1 s e = d Kd qk x ik y ) ϕ q P ) ϕ K P ) ϕ q + K P ) ϕ q + K ds e = 0, C4) ) φ q)y 11 Ω q ). C5) d1 s e 11/ = 7 5/ + 7/ 1 ) 1/4 1/ 1 10 x Y 0 Ω x ) ){ } x Y 00 Ω x ) e x). The decay amplitude for the bound-state correction h bsc fi = ) 4 7 5/ 1/4 1/ } + 7/ 1 x Y 0 Ω x ) e x). The total amplitude will be h CP0 fi = h fi +h bsc fi = 1/4 which are expressions 50) and 51). C6) ) {1 81 x Y 00 Ω x ) MCP0 1/ fi C7), C8) 1 F. E. Close and A. Kirk, Eur. Phys. J. C 1, ). U. Thoma, Eur. Phys. J. A 18, 15 00). E. Klempt, A. Zaitsev, Phys. Rept. 454, 1 007). 4 T. Barnes, Int. J. Mod. Phys. A 1, ). 5 L. Micu, Nucl. Phys. B 10, ). 6 A. Leyaouanc, L. Oliver, O. Pène and J. Raynal, Phys. Rev. D 8, 197); ibid. D 9, ); D 11, ); D 11, ); Phys. Rev. Lett. B 71, ). 7 E. S. Ackleh, T. Barnes and E. S. Swanson, Phys. Rev. D 54, ). 8 T. Barnes, F. E. Close, P. R. Page and E. S. Swanson, Phys. Rev. D 55, ). 9 T. Barnes, N. Black and P.R. Page, Phys. Rev. D 68, ). 10 T. Barnes, S. Godfrey and E. S. Swanson, Phys. Rev. D 7, ). 11 M.D. Girardeau, Phys. Rev. Lett. 7, ), ibid. J. Math. Phys. 16, ); Phys Rev. A 6, ). 1 D. Hadjimichef, G. Krein, S. Szpigel and J. S. da Veiga, Ann. of Phys. 68, ); ibid. Phys. Lett. B 67, ). 1 S. Szpigel, Interação Méson-Méson no Formalismo Fock- Tani. PhD thesis, Instituto de Física, Universidade de São Paulo, São Paulo, D. Hadjimichef, J. Haidenbauer, G. Krein, Phys. Rev. C 6, ); ibid. C 66, ). 15 D.T. da Silva, D. Hadjimichef, J. Phys. G 0, ). 16 M. L. L. Silva, D. Hadjimichef, C. A. Z. Vasconcellos, B. E. J. Bodmann, J. Phys. G, ). 17 S. Weinberg, Phys. Rev. 10, ); ibid. 11, 440

12 1 196); M. Scadron and S. Weinberg, Phys. Rev. 1, B ); M. Scadron, S. Weinberg and J. Wright, Phys. Rev. 15 B0 1964). 18 M. Oka and K. Yazaki, Prog. Theor. Phys ); ibid ). 19 E. S. Swanson, Ann. of Phys. 0, 7 199); T. Barnes and E. S. Swanson, Phys. Rev. D 46, ); T. Barnes, S. Capstick, M. D. Kovarik and E. S. Swanson, Phys. Rev. C 48, ). 0 Particle Data Group, W.-M. Yao et al., J. Phys. G, 1 006). 1 M. Strohmeier-Prešiček, T. Gutsche, R. Vinh Mau, and A. Faessler, Phys. Rev. D 60, ).