Rho meson decay in presence of magnetic field

We find a general expression for the one-loop self-energy function of neutral $\rho$-meson due to $\pi^+\pi^-$ intermediate state in a background magnetic field, valid for arbitrary magnitudes of the field. The pion propagator used in this expression is given by Schwinger, which depends on a proper-time parameter. Restricting to weak fields, we calculate the decay rate $\Gamma(\rho^0 \rightarrow \pi^+ +\pi^-)$, which changes negligibly from the vacuum value.


I. INTRODUCTION
Though some steller objects (like neutron star) were long known to possess magnetic fields [1][2][3], the realization, that such fields are created in noncentral collisions of heavy ions [4][5][6][7], has initiated looking for effects of this background magnetic field on various observables [8][9][10][11]. Thus its effect on dilepton production [12][13][14] and on resonances created in the hadron phase [15][16][17][18][19] are investigated in detail. A more involved effect of this background field, called the chiral magnetic effect, demonstrates the topological nature of the QCD vacuum [4,20,21]. Apart from these effects in heavy ion collisions, the magnetic field can enhance the symmetry breaking of a theory, e.g. it increases the magnitude of the quark condensate, which breaks the flavor symmetry of QCD [22][23][24].
Here we investigate the effect of an external magnetic field in the decay of ρ-meson in the dominant channel ρ 0 → π + + π − [15], which may affect the estimate of pion production in noncentral heavy ion collisions. This decay rate may be obtained from the imaginary part of the self energy graph of ρ-meson with two pion intermediate state (Fig III.4). The effect of the external field can be included in the decay process by taking the modified pion propagation in this field.
Such a modified (scalar and spinor) propagator in coordinate space has been derived long ago by Schwinger [25], to all orders in the external electromagnetic field, as an integral over proper time. Working in quantum electrodynamics, he used it to find corrections to Maxwell Lagrangian. But for the electron self energy function, he wrote the usual form, namely an integral over intermediate momentum k µ with the electron propagator that depends on the kinematical momentum Π µ , containing the electromagnetic potential A µ . Then the shift of the origin in k space, necessary in carrying out the k integration, cannot be made, owing to the noncommutativity of the components of Π µ . He circumvented this difficulty by an ingenius ξ-device 1 and evaluated the self energy function analytically for weak and strong magnetic fields.
In this work we write the pion self energy in coordinate space with pion propagators as given by Schwinger. There is no difficulty here as it contains no operators. The resulting expression is Fourier transformed to go over to momentum space. Having obtained the ρ 0 self energy for a general external field, we find its decay rate in a weak magnetic field.
In Section II we outline Schwinger's derivation of scalar propagator in coordinate space.
In Section III we find the ρ meson self energy, first in a general background field and then 1 Ref [26], vol II, p. 224; vol III, p. 145 specialize it to magnetic field. In Section IV we calculate the decay rate to order quadratic in the magnetic field. Finally a general discussion of our method is given in the last Section V. An Appendix evaluates the relevant integrals.

II. SCALAR PROPAGATOR
The Lagrangian for charged pions of mass m interacting with an external electromagnetic giving the equation of motion for the pion field as The pion propagator is defined as where T represents the time ordering and |0 is the vacuum for the quantum fields. The propagator satisfies We now review the steps arising in Schwinger's derivation of the exact propagator. If we introduce states labeled by space-time coordinates, G(x, x ′ ) may be written as the matrix element of an operator G when we can express Eq.(II.4) as an operator equation where Π µ = p µ − eA µ , p µ = i∂ µ and they satisfy the commutation relations, The operator equation (II.6) has the formal solution (II.9) This notation emphasizes that U(s) may be regarded as the operator describing the dynamics of a particle governed by the Hamiltonian 'H' in the proper time parameter 's'. The spacetme coordinate x µ = (t, x) of the particle depends on this parameter.
In the Heisenberg representation the operator x µ and Π µ have the 'time' dependence and the base ket and bra evolve as Then the construction of G(x ′ , x) reduces to the evaluation of which is the transformation function from a state in which the position operator x µ (s = 0) has value x ′ µ , to a state in which x µ (s) has the value x ′′ µ . The equations of motion for the operators following from Eq.(II.10) are The transformation function itself can be found by solving the differential equation satisfied by it, i d ds x ′′ ; s|x ′ ; 0 = x ′′ ; s|H(x(s), Π(s))|x ′ ; 0 , (II.14) and along with a similar one for Π µ (0), with boundary conditions Using (II.18) and (II. 19) and the antisymmetry of the field tensor (F µν = −F νµ ) we get To evaluate the matrix element on the right of Eq.(II.14), we need to order the operator x(s) to the left of x(0), which would require the commutator Then we get with tr indicating the trace over 4 × 4 matrices. Eq.(II.14) can now be solved as The s-independent function C can be found by solving Eq.(II.15) has vanishing curl, it can be written as where the integration in the phase factor Φ runs on a straight line between x ′ and x ′′ and the F term vanishes. The constant C is given by the first boundary condition of (II.16) as This integral is evaluated in the Appendix, to give C = −i/(4π) 2 .
We finally get the transformation function as which in turn gives the propagator In our work below we shall encounter the product of propagators, derivatives acting on them. Without the derivatives, the phase factor in G would cancel out mutually. In presence of derivatives we can still get rid of the phase factors, if we make a gauge choice in the potential, replacing A µ with when Eq.(II.4) is satisfied by G without the phase factor Φ 4 . In the following we choose this gauge to write G(x ′′ , x ′ ) without this phase.

III. ρ SELF-ENERGY IN EXTERNAL FIELD
We now express the self energy graph of Fig

A. General field
We take a phenomenological Lagrangian forρππ interaction, The coupling g can be found from the experimental decay width Γ(ρ 0 → π + π − ) = 149 MeV in vacuum to give g = 6.0. We now work out the complete ρ-propagator to order g 2 in the interaction representation. We take the ρ field as free 5 , but the pion field lives in the background electromagnetic field. Contracting the fields for graph of Fig. 1 where Σ is the self-energy tensor involving the two pion propagators. The propagators are distinguished only by their proper times s 1 and s 2 , over which they are integrated. Carrying out the derivatives contained in L int on these propagators, we get [27,28] Here L 1 = L(s 1 ), L 2 = L(s 2 ) and similarly for R 1 , R 2 .
Having obtained equation (III.2) in configuration space, we go over to momentum space by taking Fourier transforms. Letting K µν (x) to denote any of the D µν (x), D ′ µν (x) and Σ µν (x), their Fourier transforms are defined as and Eq. (III.2) becomes The vacuum and the complete propagator of ρ meson are given by where D ′ (q 2 ) is the function we want to find.
A simplification results on noting that the ρ field is coupled to a conserved pion current in the interaction lagrangian given by Eq. (III.1). As a result, contracting q µ and q ν of the ρ propagators with Σ µν in the second term of Eq.(III.5) yield zero. We are then left with the metric tensor in ρ propagator. Contracting further the indices σ and λ we get Here Σ(q) is the Fourier transform of Σ µν (x) after contracting the indices, Clearly the expression for Σ is divergent, to which we have to add renormalization counterterms. Beside cancelling the divergent pieces, we shall choose the finite pieces in the counterterms, such that the total self energy Σ tot satisfies Then m ρ will remain the physical ρ-meson mass and g, the renormalized ρππ coupling. We shall come back to this renormalization in Section IV below.
Including the sum of all reducible graphs in Eq.(III.8) we get the Dyson-Schwinger equation for the ρ propagator as .

B. Pure magnetic field
We derived above the expressions for the pion propagator and the consequent ρ-meson self energy in a general background external field F µν . We now specialize this field to magnetic field B in the z direction, i.e.  , 19) which is evaluated in the Appendix yielding Note that all the parameters α, β, γ and δ are positive. The where Λ ′ is related to Λ by the change of variable. To write Λ ′ , we define a new set of variables from Eq.(III.17), .
(III. 22) In terms of these variables, we have Another form of Σ, which will be useful for the discussion below may be obtained by scaling t with eB, that is, we sett = eBt, when Σ becomes . (III.26) IV. ρ-MESON DECAY As already stated, the self energy function in Eq.(III.21) is valid for momenta, for which the ρ-meson cannot decay into a pion pair. To calculate this decay rate we therefore need to continue Eq.(III.21) beyond such momenta [29]. This process of analytic continuation is immediate, if we can evaluate the t integral analytically. However the exact expression (III.24) for Σ containing various hyperbolic functions, makes it difficult to do so. The procedure here is to consider separately weak (eB < m 2 ) and strong (eB > m 2 ) fields 7 .
As strong fields are considered extensively in the literature [15][16][17][18], we take up the case of weak fields, which is realized, in particular, in the hadronic phase of noncentral heavy ion collisions. In this case the exponential in Eq.(III.25) shows that only correspondingly small values of t can contribute. We can then expand the different functions in powers of t. To get the leading effect, we need to keep only the first two terms in their expansions. After some algebra, we get the self energy as Here a numerical factor of 4 4 has been put into the coupling constant factor. Also for convenience, the variable u is replaced by 1 2 (1 + v).
As the terms in D do not depend on the magnetic field and we want to calculate the change in ρ meson decay width in this magnetic field, it is clear that our calculation will not involve D. However we want to show the nature of those terms. To this end, consider the first term in D, behaving as t −1 . Integrating partially w.r.t v, it gives g 2π where the divergence at t = 0 is isolated in the first term. The second term in D behaving as t −2 can also be put in a similar form after integrating twice partially w.r.t v. These are the local divergent terms, which we expected in Section III for the general field case and are analogous to those appearing in loop integrals over intermediate momenta in conventional field theory. As a consistency check in our calculation, let us note that though individual terms in Λ ′ given by Eq.(III.23) do contain (eB) 2 dependent divergent terms, they cancel out in the complete expression for Λ ′ , showing that divergences originate only from the vacuum piece of self energy, as expected.
Going back to the eB dependent self energy given by the F terms in Eq.(IV.2), we rewrite it as We are now in a position to carry out the analytic continuation mentioned at the beginning of the section. If we hold the q 2 variable in the region q 2 < 4m 2 , the t integration is well defined and can be integrated trivially to give .
It can now be continued for q 2 < 4m 2 in the q 2 plane with a cut along the real axis for 4m 2 < q 2 < ∞. To display this cut structure explicitly, we write Σ(q 2 ) as a dispersion integral by changing the integration variable v to q ′2 given by v = 1 − 4m 2 /q ′2 , getting It's imaginary part is given by the discontinuity across the cut Taking m 2 /m 2 ρ = 1/30, it becomes (IV.10) For eB < m 2 and q 2 ⊥ < m 2 ρ , it gives Γ eB < 0.6 MeV. The smallness of Γ eB may be explained by the fact that while the (small) pion mass is the scale entering in the self energy loop, it is evaluated at a (large) external momentum of ρ meson mass. Also note that there is no pion mass in the denominator of Eq.(IV.8). It is protected by chiral symmetry (m → 0), according to which physical quantities must be finite in this limit.
In passing, we note that the effect of temperature on the decay width of ρ-meson has been discussed extensively in the literature [30]. Here we have investigated the effect of weak magnetic fields on the same quantity and found it to be negligible with respect to the thermal effects.

V. DISCUSSION
In earlier calculations of hadron properties in a magnetic field [15][16][17][18], the majority of works consider strong fields, taking the contribution of the leading Landau level for the system. A result to note at this point is that for strong enough fields the main decay channel, such as ρ 0 → π + + π − that we are considering, may become closed. It is due to the generation of an effective pion massm 2 = m 2 + eB, causing the phase space for the process to shrink as the magnetic field becomes stronger [15].
In the present work we investigate the decay by setting up a general framework, valid for both weak and strong magnetic fields. It is obtained by writing the ππ loop in the correction to the ρ propagator in configuration space, with pion propagator as given by Schwinger [25]. When Fourier transformed, it gives the ρ meson self energy (III.21) as an integral over proper times, which is defined for momenta below the two-pion threshold.
If we now restrict the general representation Eq.(III.21) to weak fields (eB < m), the exponential factor in it (or quivalently Eq.(III.25)) shows the leading contribution to arise from the neighbourhood of proper time t = 0, when we can expand the hyperbolic functions in powers of t. Still remaining below the two-pion threshold, we can integrate the resulting terms to get a series in powers of (eB) 2 . These terms can be simply continued beyond the threshold and the imaginary part of the self energy giving the decay width can be determined. In this work we retain only the (eB) 2 terms, though calculation of higher order terms is also straightforward. As we show at the end of Section IV, the change in the decay width from the vacuum value turns out to be negligibly small.
So far we only discussed the effect of weak magnetic fields. But as already emphasized, strong field effects can also be obtained from the same general formula Eq.(III.21). For eB > m 2 , the exponential in this formula shows that large values of t would also contribute.
It is thus simple to keep the leading term in different hyperbolic functions. Collecting the exponentials in Eq.(III.21), we get giving the effective pion mass, as mentioned above.
There are at least two other methods of calculating the decay rate. One is Schwinger's ξ-device mentioned in Section I and the other is the Ritus method of eigenfunction expansion [31]. It will be interesting to get comparable values from these methods.

VI. ACKNOWLEDGEMENT
The work of AB is supported by the Depatment of Atomic Energy (DAE), India.
Here we evaluate the integrals in Eq. (II.28) and (III. 19), paying attention to the phases appearing in the manipulations. First consider For J we put x 2 1 /4s = u to get To avoid oscillations in the integrand, we take the contour of Λ(s, u, q , q ⊥ ) = e −im 2 s d 4 x e iq·x −αx 2 ⊥ + βx 2 + 16i(γ + δ) exp iγx 2 ⊥ − iδx 2 . Here the basic integrals are in terms of which the Fourier transform may be written as The basic integrals are generally of the same form as J and K, if we complete the squares in the exponents. Thus on substituting x ⊥ → x ⊥ + q ⊥ /2γ and using polar coordinates. Taking a contour in the first quadrant, we get L 1 = e iπ/2 π γ e −iq 2 ⊥ /4γ . (A.7) In the same way we can evaluate the integrals L 2 and L 3 by taking contours respectively in the fourth and first quadrants, L 2 = e −iπ/4 √ π √ δ e iq 2 0 /4δ , (A.8) Putting these values of integrals L i , i = 1, 2, 3, in (A.6) we get Λ as given by Eq.(III.20) in the text.