General structure of neutral $\rho$ meson self energy and its spectral properties in hot and dense magnetized medium

The one loop self energy of the neutral rho meson is obtained for the effective $\rho\pi\pi$ and $\rho NN$ interaction at finite temperature and density in presence of a constant background magnetic field of arbitrary strength. In our approach, the eB-dependent vacuum part of the self energy is extracted by means of dimensional regularization where the ultraviolet divergences corresponding to the pure vacuum self energy manifest as the pole singularities of gamma as well as Hurwitz zeta functions. This improved regularization procedure consistently reproduces the expected results in the vanishing magnetic field limit and can be used quite generally in other self energy calculations dealing with arbitrary magnetic field strength. In presence of the external magnetic field, the general Lorentz structure for the in-medium vector boson self energy is derived which can also be implemented in case of the gauge osons such as photons and gluons. It has been shown that with vanishing perpendicular momentum of the external particle, essentially two form factors are sufficient to describe the self energy completely. Consequently, two distinct modes are observed in the study of the effective mass, dispersion relations and the spectral function of $\rho^0$ where one of the modes possesses two fold degeneracy. For large baryonic chemical potential, it is observed that the critical magnetic field required to block the $\rho^0\rightarrow\pi^+\pi^-$ decay channel increases significantly with temperature. However, in case of smaller values reaching down to vanishing chemical potential, the critical field follows the opposite trend. The requirement of smaller critical field makes the stability of $\rho^0$ against the $\pi^+\pi^-$ decay more likely to be observed in heavy-ion collision experiments in LHC.


I. INTRODUCTION
In non-central heavy-ion collisions at LHC, the relative motion of the ions themselves can generate strong decaying magnetic pulse of the order eB ∼15m 2 π (B ∼ 5 × 10 15 Tesla) [1]. While some of the studies support rapid decrease in the magnitude [2,3], an adiabatic decay is expected [4][5][6] due to the high conductivity of the produced medium. In spite of the ambiguities, the intensity of the produced magnetic field being much larger than the typical QCD scale, the possibility of magnetic modifications of different properties of the produced extreme state of matter can not be refuted completely. In general, high intensity magnetic fields can play a significant role in many astrophysical and cosmological phenomena [7][8][9][10]. Moreover, the magnetic influence on the properties of magnetars adds to the motivation of studying high density matter in presence of extreme magnetic fields [11][12][13][14][15][16][17][18].
The study of ρ meson properties like the effective mass and dispersion relations are important in the context of magnetic field induced vacuum superconductivity [19][20][21][22][23][24][25][26]. Using NJL model in presence of magnetic background, Liu. et.al. have shown that the charged rho condensation in vacuum occurs at critical magnetic field eB c ∼ 0.2GeV 2 [25]. Generalization of the study to finite temperature and density shows that the condensation survives even in presence of finite temperature and density [27]. At vanishing chemical potential, the corresponding critical magnetic field is observed to lie in the range 0.2 -0.6 GeV 2 for temperatures in between 0.2-0.5 GeV. However, the neutral ρ meson in vacuum, having no trivial Landau shifts in the energy eigenvalue, shows a slow decrease in the effective mass [28] in weak magnetic field region. Thus, if neutral rho condensation is possible, extremely large magnetic field values will be required to observe the condensation. It should be mentioned here that it has been shown using NJL model that the effective mass of ρ 0 meson in fact increases at higher values of magnetic fields showing no possibility of condensation [25]. In this scenario, ρ 0 → π + π − decay may serve as an important probe to observe the influence of the magnetic field. As argued in Ref. [19], even if point like ρ 0 meson is considered without any influence by magnetic field, there exists a critical value of the external magnetic field for which the ρ 0 to π + π − decay stops due to the trivial enhancement of the charged pion mass. Later the magnetic modification arising from the loop corrections are taken into account at weak [28,29] as well as at strong field limits [30] at zero temperature. An immediate generalization of the previous works will be to incorporate the medium effects of the ρ 0 meson which may reflect in the modification of the decay rate and the required critical magnetic field. It should be noted here that apart from being important in the study of dense hadronic matter at extreme conditions usually expected to be present within compact stars, the incorporation of the medium effects is also essential for the proper estimation of pion production in non-central heavy ion collisions.
In this work we focus on the temperature and density modifications of neutral ρ meson properties in presence of a static homogeneous magnetic background. The one loop self energy of ρ meson is calculated for the effective ρππ and ρN N interaction with magnetically modified pion and nucleon propagators corresponding to general field strength. After decomposing the self energy in terms of the form factors, the decay width for ρ 0 → π + π − channel is obtained. It should be mentioned here that the spectral properties of rho meson in presence of finite temperature and magnetic field have been studied in our earlier work [31]. However, unlike the previous case, dimensional regularization technique is used here to extract the ultraviolet divergence as pole singularities of gamma and Hurwitz zeta functions [32]. Also, instead of considering only the spin averaged thermal self energy contribution, the general Lorentz structure has been addressed in detail. Apart from the technical differences, the density dependence arising from the charged nucleon loop serves as the most important extension of the previous study. Its importance can be understood as follows. It is well known that the general expression of decay width is related to the imaginary part of the self energy. Now, as far as the ρ 0 → π + π − decay is concerned, the invariant mass regime of interest does not allow the nucleon loop to directly contribute to the imaginary part as the unitary cut threshold of N N loop begins at much higher value. However, it should be noted that in the rest frame of the decaying particle, the decay width depends on its effective mass. The contribution from the nucleon loop incorporates significant modification in the effective mass of ρ 0 which in turn influences the decay. As we shall see, the critical field required to stabilize the neutral ρ against the π + π − decay has a non-trivial dependence on the baryonic chemical potential.
The article is organized as follows. In Sec. II the vacuum self energy of ρ is discussed followed by evaluation of the inmedium ρ self-energy at zero magnetic field in Sec. III. Next in Sec. IV, the in-medium self energy at non-zero external magnetic field is presented. Sec. V is devoted to the discussion of the general Lorentz structure of the in-medium self energy function in presence of a constant background magnetic field. After addressing the Lorentz structure of the interacting ρ propagator in Sec.VI, the analytic structure of the self energy is discussed in Sec. VII. Sec. VIII contains the numerical results. Finally we summarize and conclude in Sec. IX. Some of the relevant calculational details are provided in the Appendix.
The evaluations of the momentum integrals in Eqs. (3) and (4) (Π µν N ) pure-vac (q) = (q 2 g µν − q µ q ν ) where ∆ π and ∆ N are defined in Eqs. (B6) and (B7). As can be seen from the above equations, the vacuum self energy is divergent and scale dependent which renormalizes the bare ρ 0 mass to its physical mass after adding proper vacuum counter terms in the Lagrangian. The particular Lorentz structure in the above equations renders the self energy transverse to the ρ 0 momentum i.e. q µ Π µν vac = 0.

III. ρ 0 SELF ENERGY IN THE MEDIUM
In order to calculate the ρ 0 self energy at finite temperature and density, we employ the real time formalism of finite temperature field theory where all the two point correlation functions such as the propagator and the self energy become 2 × 2 matrices in the thermal space [34,35]. However, they can be put in a diagonal form where the diagonal elements can be obtained from any one component (say the 11-component) of the said 2 × 2 matrix. The 11-components of real time thermal pion and nucleon propagators are where and f ± (x) are respectively the Bose-Einstein and Fermi-Dirac distribution functions corresponding to pions and nucleons: Here, Θ(x) is the unit step function, u µ is the medium four-velocity; T and µ B are respectively the temperature and baryon chemical potential of the medium. In the local rest frame (LRF) of the medium, u µ LRF ≡ (1, 0). For the evaluation of the 11-component of the thermal self energy matrix, the vacuum pion and nucleon propagators in Eqs. (3) and (4) are replaced by the respective 11-components of the thermal propagators given in Eqs. (11) and (12) as [35] (Π µν π ) 11 (q) = i The analytic thermal self energy function of ρ 0 denoted by a bar ReΠ µν (q 0 , q) = ReΠ µν π (q 0 , q)+ReΠ µν N (q 0 , q) is related to the above quantities by the relations [35] ReΠ µν π,N (q 0 , q) = ReΠ µν π,N 11 ImΠ µν π,N (q 0 , q) = sign q 0 tanh where, sign (x) = Θ(x) − Θ(−x). Rewriting Eqs. (11) and (12) as and substituting into Eqs. (14) and (15) and performing the dk 0 integration (using the Dirac delta functions) followed by using Eqs. (16) and (17) we get the real parts as, and the imaginary parts as ImΠ µν π (q 0 , q) = −sign q 0 tanh where, P denotes the Cauchy Principal value integration, ω k = m 2 π + k 2 , Ω k = m 2 N + k 2 and N N (q, k) is defined in Eq. (B3).

IV. ρ 0 SELF ENERGY IN THE MAGNETIZED MEDIUM
In presence of the external magnetic field B = Bẑ, the propagations of the charged pion and proton are modified. One of the possible ways to incorporate the effect of external magnetic field is the Schwinger proper time formalism in which the 11-components of charged pion and proton propagators respectively become [36,37] where, ∆ B (k, m π ) and S B (k, m N ) denote the momentum space vacuum (zero temperature) Schwinger proper time propagators for charged pion and proton respectively [36]: In the above equations, e = |e| is the charge of the proton; the four-vector k is decomposed into k = (k + k ⊥ ) where k µ = g µν k ν and k µ ⊥ = g µν ⊥ k ν corresponding to the decomposition of the metric tensor g µν = (g µν + g µν ⊥ ) with g µν = diag(1, 0, 0, −1) and g µν ⊥ = diag(0, −1, −1, 0). The above decomposition can be done in a Lorentz covariant way by introducing another four-vector where G µν = 1 2 ǫ µναβ F αβ is the dual of the electromagnetic field tensor F µν . In the local rest frame of the medium, b µ LRF ≡ (0, 0, 0, 1), which is the direction of the external magnetic field. Using b µ , we can write It is important to note that, the coordinate space Schwinger propagator contains a gauge dependent translationally non-invariant phase factor. However, for the one-loop graphs containing equally charged particle in the loop, the phase factor gets canceled and the momentum space propagator is sufficient for the calculation of the self energy. The proper time integral in Eqs. (26) and (27) can be performed in order to express the propagators as a sum over discrete Landau levels as where, with α k = −k 2 ⊥ /eB. Here, L a l (z) denotes the generalized Laguerre polynomial with L a −1 (z) = 0 and L l (z) = L 0 l (z). We now rewrite Eqs. (24) and (25) using Eqs. (30) and (31) as where we have defined the Landau level dependent "dimensionally reduced effective masses" (as a consequence of dimensional reduction) of pion and proton as m l = m 2 π + (2l + 1)eB and M l = m 2 N + 2leB .
We now replace the 11-component of the charged pion and proton propagators in Eqs. (14) and (15) as D 11 → D 11 B , S 11 p → S 11 B i.e by the respective magnetized ones given in Eqs. (33) and (34) and then perform the dk 0 integrations (using the Dirac delta functions). Following Eqs. (16) and (17) we get the thermal self energy functions under external magnetic field which we will denote by a double bar to distinguish them from the thermal self energy functions in the absence of magnetic field. Their explicit expressions are given by where, The first terms on the RHS of Eqs. (37) and (39) are the contributions from the neutron-neutron loop which are not affected by the external magnetic field. The last terms on the RHS of Eqs. (36) and (37) are the contributions from ππ and proton-proton loop which depend on the external magnetic field but independent of temperature. Their explicit forms are given by It is important to note that, the above quantities respectively contain the divergent pure vacuum contributions (Π µν π ) pure-vac (q) and 1 2 (Π µν N ) pure-vac (q) in a nontrivial way (as the above equations seem to appear non-perturbative in eB). In contrast, for the case of weak magnetic field expansion of the Schwinger propagator, the pure vacuum contribution to the self energy trivially decouples from the magnetic field dependent terms. Since we are working with the full propagator including all the Landau levels, we have to properly regularize the above expressions in order to extract the pure vacuum contributions from these quantities. We use dimensional regularization in which the ultraviolet divergence appear as the pole of Gamma and Hurwitz zeta function the details of which are provided in the Appendices C and D. Here, we take the transverse momentum of ρ 0 to be zero i.e. q ⊥ = 0 which makes substantial simplifications of the analytic calculations. The final result can be read off from Eqs. (C8) and (D8) as where, the scale dependent divergent pure-vacuum parts are completely decoupled as the first term on the RHS of the above equation; the scale independent and finite "eB-dependent vacuum contribution" to the real part of the self energy functions are Eqs. (46) and (47) imply that the vacuum counter terms are sufficient to renormalize the theory and thus the external magnetic field does not create additional divergences. For q ⊥ = 0, the d 2 k ⊥ integrals in Eqs. (36)-(39) can be analytically performed (see Appendix E) and we finally get, ImΠ µν π (q 0 , q z ) = −sign q 0 tanh where,Ñ µν π,nl (q , k ) andÑ µν p,nl (q , k ) can be read off from Eq. (E5) and (E8). The presence of Kronecker delta functions in the expressions ofÑ µν π,nl (q , k ) andÑ µν p,nl (q , k ) has eliminated one of the double sums or in other words, the sum over index l now runs from (n − 1) to (n + 1).

V. LORENTZ STRUCTURE OF THE VECTOR BOSON SELF ENERGY IN MAGNETIZED MEDIUM
In this section, we will derive the tensorial decomposition of the massive vector boson self energy. We note that, the self energy Π µν (q) being a second rank tensor, has sixteen components which will mix among themselves with the change of frame. It is useful to use linearly independent basis tensors (constructed with the available vectors and tensors) to express Π µν (q) so that the form factors (corresponding to each basis) remain Lorentz invariant. This will also enable one to solve the Dyson-Schwinger equation in order to obtain the complete interacting vector boson propagator. In order to proceed, we first note that the vector boson self energy satisfies the following constrain Π µν (q) = Π νµ (q) and q µ Π µν (q) = 0 .
Let us first consider the pure vacuum case i.e. for zero temperature and zero external magnetic field. In this case, the only available vector is the momentum q µ along with the metric tensor g µν so that Π µν (q) is a linear combination of q µ q ν and g µν i.e Π µν (q) = (α 1 g µν + α 2 q µ q ν ). Imposing the constrains of Eq. (54), we get α 1 + α 2 q 2 = 0 which makes the only possible Lorentz structure of the self energy as where the Lorentz invariant form factor α 1 = α 1 (q 2 ) = 1 3 Π µ µ . Note that, with q µ and g µν , the only possible Lorentz scalar that can be formed by contracting with Π µν (q) is the quantity g µν Π µν = Π µ µ implying the existence of only one form factor.
We now consider the case with finite temperature but zero magnetic field. In this case we have an additional four vector u µ (medium four-velocity) along with q µ and g µν . This makes Π µν to be a linear combination of g µν , q µ q ν , u µ u ν , q µ u ν and q ν u µ i.e. Π µν (q) = (α 1 g µν + α 2 q µ q ν + α 3 u µ u ν + α 4 q µ u ν + α 5 q ν q µ ) However, imposing the constrains in Eq. (54), we find the following relationship among the coefficients which makes only two of the coefficients independent. Choosing α 1 and α 2 as independent, we get, where the Lorentz invariant form factors α 1 = α 1 (q 2 , q · u) and α 2 = α 2 (q 2 , q · u) can be obtained by contracting both side of the above equations with g µν and u µ u ν so that the form factors will become functions of the Lorentz scalars g µν Π µν = Π µ µ and u µ u ν Π µν . Note that, with q µ , u µ and g µν , only two possible Lorentz scalars that can be formed by contracting with Π µν (q) are the quantities Π µ µ and u µ u ν Π µν implying the existence of only two form factors. Unlike the pure vacuum case given in Eq. (55), here the decomposition of Π µν in Eq. (60) is not unique. As already mentioned, it is useful to construct linearly independent (and mutually orthogonal) basis tensors (note that the basis tensors within square brackets in Eq. (60) are not mutually orthogonal). One such choice of orthogonal tensor basis could be which is constructed from u µ by subtracting out its projection along q µ . It is easy to check that P µν 1 and P µν 2 satisfy all the properties of projection tensors i.e.
Therefore, Π µν can be written as where the form factors are Care should be taken when considering the special case like q = 0 [35]. To see this, let us consider q i = | q|n i so that the spatial components of the projectors at q = 0 become (in the LRF) This implies that the spatial components of self energy at vanishing three momentum depend on the direction of q even at | q| = 0. This ambiguity is eliminated by setting additional constraint on the form factors as Π 1 (q 0 , q = 0) = Π 2 (q 0 , q = 0). Following the same strategy, we now construct suitable orthogonal tensor basis for the vector bososn self energy at finite temperature under external magnetic field. In this case we have an additional four vector b µ (corresponding to the magnetic field direction) along with q µ , u µ and g µν . This makes the symmetric Π µν to be a linear combination of seven tensors as However, imposing the constrains in Eq. (54), we find the following relationship among the coefficients which makes only (7-3=4) four of the coefficients independent. The Lorentz invariant form factors α i = α i (q 2 , q·u, q·b) with i = 1, 2, ..., 7 can be obtained by contracting both side of the above equations separately with g µν , u µ u ν , b µ b ν and u µ b ν so that the form factors will become functions of the Lorentz scalars Π µ µ , u µ u ν Π µν , b µ b ν Π µν and u µ b ν Π µν . Note that, with q µ , u µ , b µ and g µν , only four possible Lorentz scalars that can be formed by contracting with Π µν (q) are the quantities Π µ µ , u µ u ν Π µν , b µ b ν Π µν and u µ b ν Π µν implying the existence of only four form factors. Like the finite temperature case, here the the decomposition of Π µν is also not unique. One convenient choice of tensor basis could be The basis tensors in Eqs. (72)-(75) satisfy the following relations: Using the basis given in Eqs. (72)-(75), the self energy at finite temperature under external magnetic can be written as where the form factors are obtained as Analogous to the case of only finite temperature, care should be taken while considering the special case q ⊥ = 0. To see this, let us consider q i ⊥ = | q ⊥ |n i with i = 1, 2 so that the following components of self energy at vanishing q ⊥ become (in the LRF) which depend on the direction of q ⊥ even at q ⊥ = 0. This ambiguity is eliminated by setting additional constraints on the form factors as

VI. THE INTERACTING ρ MESON PROPAGATOR AND ITS LORENTZ STRUCTURE
Let us first consider the zero temperature and zero magnetic field case for which the complete interacting ρ propagator D µν is obtained by solving the Dyson-Schwinger equation is the free vacuum Feynman propagator and Π µν is the one-loop self energy of ρ meson which has the Lorentz structure given in Eq. (55) as with the form factor Π = 1 3 Π µ µ . In order to solve Eq. (93), we rewrite it as Substituting Π µν from Eq. (95) in the above equation, we get the inverse of the complete propagator which can be inverted using the relation D µα (D αν ) −1 = g µ ν to obtain the complete propagator as We now consider the case of finite temperature and zero magnetic field. As already mentioned in Sec. III, in RTF of finite temperature field theory all the two point correlation functions become 2 × 2 matrices in thermal space. In this case the Dyson-Schwinger equation also becomes a matrix equation [35] Each term of the above equation can be diagonalized in terms of the respective analytic functions (denoted by a bar) so that the above equation becomes an algebric one where ∆ µν = ∆ µν . The above equation can be rewritten as In this case, the Lorentz structure of the thermal self energy function is given in Eq. (64) as where the projection tensors and form factors are respectively defined in Eqs. (61) and (65). Substituting the above equation in Eq. (100), we get the inverse of the complete propagator. In order to obtain the complete propagator, we write and use the relation D µα D αν −1 = g µ ν to extract A 1 , A 2 and ξ. The final form of the complete interacting thermal propagator is obtained as Finally we consider the case with both finite temperature and external magnetic field. In this case we need to solve the Dyson-Schwinger equation where a double bar is used to denote thermal self energy function and complete propagator under external magnetic field as discussed in Sec. IV. In this case, the Lorentz structure of the thermal self energy function is given in Eq. (85) as where the basis tensors and form factors are given in Eqs. (72)-(75) and (86)-(89). Substituting the above equation in Eq. (104), we get the inverse of the complete propagator. In order to obtain the complete propagator, we write and use the relation D µα D αν −1 = g µ ν to extract the coefficients as (111)

VII. ANALYTIC STRUCTURE OF THE SELF ENERGY
In this work, we have considered the transverse momentum of the rho meson to be zero i.e. q ⊥ = 0. As shown in Eq. (92), for the special case q ⊥ = 0, the additional constraints to be imposed on the form factors are Using the above constraints, we get from Eqs. (86)-(89) which imply that we need to calculate only the two quantities quantities Π µ µ and u µ u ν Π µν = Π 00 . These are obtained from Eqs. (50)-(53) by contracting them with g µν and u µ u ν . This essentially means replacing N µν for all the loops with N µ µ or N 00 , an explicit list for which has been provided in Appendix F. Let us now discuss the analytic structure of the self energy functions. We first consider the zero magnetic field case. The imaginary part of the self energy function for ππ and NN loops as given in Eqs. (22) and (23) each contains four Dirac delta functions. These delta functions represent energy-momentum conservation and they are non vanishing in certain kinematic domain. They are termed as the Unitary-I, Unitary-II, Landau-II and Landau-I cuts as they appear in those equations. The kinematic regions for the Unitary-I and Unitary-II cuts are given by [35] q 2 + 4m 2 L < q 0 < ∞ and −∞ < q 0 < − q 2 + 4m 2 L whereas the same for the two Landau cuts are |q 0 | < | q| where m L is the mass of the loop particle i.e. m L = m π or m N . These cuts correspond to different physical processes such as decay or scattering. For example, Unitary cuts correspond to the decay of ρ 0 into a π + π − or NN pair and the Landau cuts correspond to the scattering of a ρ 0 with a pion or nucleon producing the same in the final state along with their time reversed processes. If we restrict ourselves to the physical timelike kinematic regions defined in terms of q 0 > 0 and q 2 > 0, then only the Unitary-I cut contributes. It is important to note that, a non-trivial Landau cut appears in the physical timelike region only if the loop particles have different masses and lie in the kinematic domain | q| < q 0 < q 2 + ∆m 2 where ∆m is the mass difference of the loop particles.
Let us now consider the case of both finite temperature and non zero external magnetic field. In this case the imaginary parts of the self energy as given in Eqs. (52) and (53) also contain four Dirac delta functions corresponding to the Unitary and Landau cuts. It is important to note that the arguments of the delta functions contain only the longitudinal dynamics (because of dimensional reduction) which implies that the analytic structure of the self energy functions will only depend on the longitudinal momentum of ρ. On the other hand, the transverse dynamics has appeared as Landau level dependent "dimensionally reduced effective mass" to the loop particles as given in Eq. (35). Therefore, even if the loop particles have the same masses, a non-trivial Landau cut may appear in the physical timelike kinematic domain if the two loop particles reside in different Landau levels. Physically, this means that ρ 0 can get absorbed in a scattering with a pion or a proton in a lower Landau level producing another pion or proton in a higher Landau level as the final state. A detailed discussions on the analytic structure in presence of external magnetic field can be found in Refs. [31,38]. The Unitary-I and Unitary-II terms for the ππ loop are non-vanishing in the kinematic domains q 2 z + 4(m 2 π + eB) < q 0 < ∞ and −∞ < q 0 < − q 2 z + 4(m 2 π + eB) whereas the kinematic domain for both the Landau cuts is The corresponding kinematic domains for the NN loop are q 2 z + 4m 2 N < q 0 < ∞ and −∞ < q 0 < − q 2 z + 4m 2 N for the Unitary-I and Unitary-II cuts respectively and for the Landau cuts. Note that, the threshold of the Landau cuts appears when the "dimensionally reduced effective mass" difference between the loop particles is the maximum. As can be seen from Eqs. (52) and (53), for a particular value of the index n, the sum over the index l runs only for three values (n − 1), n and (n + 1) which implies that, the Landau level difference between the loop particles can be at most one. Thus the maximum difference in their "dimensionally reduced effective mass" appears when one of them is at the lowest Landau level and the other one is at the first Landau level which in turn defines the Landau cut threshold in Eqs. (116) and (117). We now simplify the expressions of the imaginary parts given in Eqs. (22), (23), (52) and (53) by evaluating one of the integrals using the Dirac delta functions. For the imaginary parts at zero magnetic field, we evaluate the d(cos θ) integrals and get (after imposing the kinematic restrictions discussed above), where, with λ(x, y, z) = x 2 + y 2 + z 2 − 2xy − 2yz − 2zx being the Källén function. For the imaginary parts at finite magnetic field, we evaluate the dk z integrals in Eqs. (52) and (53) using the Dirac delta functions. The imaginary part due to ππ loop simplifies to where, with, The corresponding expression of the imaginary part due to NN loop reads where, with, The first term on the RHS of Eq. (137) is the contribution from the neutron-neutron loop (which is not affected by the external magnetic field ) whose simplified form is given in Eq. (118).

VIII. NUMERICAL RESULTS
We begin this section by presenting the real and imaginary parts of the in-medium self energy functions of ρ 0 . As can be seen from Eqs. (89)-(115), we have only two non-zero form factors for the self energy which are Π α and Π β for q ⊥ = 0. Let us first consider the zero magnetic field case for which the imaginary and real parts of Π α and Π β are depicted in Figs. 2 and 3 respectively. In Fig. 2(a), ImΠ α and ImΠ β due to ππ loop are plotted as a function of invariant mass ( q 2 ) of ρ 0 for vacuum as well as for medium (T = 160 MeV and µ B = 400 MeV) with q z = 250 MeV. It is to be understood that in the case of vacuum the two form factors are equal. In this case, the only contribution comes from the Unitary-I cut which starts at 2m π in the invariant mass axis. With the increase in temperature , the degeneracy between the form factor get lifted as well as they are enhanced with respect to the vacuum. This is due to the enhancement of the thermal factor in Eq. (120) which increases the available phase space with the increase in temperature. The corresponding results for the NN loop is shown in Fig. 2(b) for which the threshold of the Unitary-I cut is 2m N . In this case, with the increase in temperature and density, the imaginary part decreases slightly with respect to the vacuum which can be understood from Eq. We now turn on the external magnetic field. For the check of consistency of the calculation at non-zero magnetic field, it is essential that eB → 0 limit of non-zero magnetic field results reproduces the eB = 0 one. In order to take the eB → 0 limit numerically, we have considered upto 500 Landau levels for a convergent result. We have shown the imaginary part of the self energy as a function of invariant mass of ρ 0 with longitudinal momentum q z = 250 MeV at temperature T = 130 MeV and at baryon chemical potential µ B = 300 MeV for the two cases: eB = 0 and eB → 0 in Fig. 4 separately for the ππ and NN loops. Fig. 4(a) shows ImΠ α for the ππ loop in which the eB → 0 graph has a series of spikes infinitesimally separated from each other all over the whole invariant mass region whereas the eB = 0 graph is finite and well behaved. Interestingly, the eB → 0 graph does not miss the eB = 0 curve which implies that when average is done, the eB = 0 line will be exactly reproduced. The appearance of these spikes are due to the "threshold singularities" [31,38,39] at each Landau level as can be understood from Eq. (132) where the Källén function goes to zero at each threshold of the Unitary and Landau cuts defined in terms of the unit step functions therein, which is a consequence of the dimensional reduction. In order to extract physical and finite results out of these spikes, we have used Ehrenfest's coarse-graining (CG) [38,40,41]. In this method, the whole invariant mass region has been discretized in small bins followed by bin averages. In other words, the self energy at a given q 2 is approximated by its average over the neighbourhood around that point. This in turn smears out the spike like structures. As can be seen in the figure, after CG, ImΠ α exactly matches with the analytic eB = 0 graph. The corresponding comparison of eB → 0 and eB = 0 result for ImΠ β due to ππ loop is shown in Fig. 4(b). In this case, eB → 0 graph is finite and free from the threshold singularities and it matches exactly with the eB = 0 graph. The absence of the threshold singularities in this case is due to an overall factor of Källén functions coming fromÑ 00 π,nl in Eq. (133) which cancels the Källén functions in the denominator of Eq. (132). Thus the ImΠ β due to the ππ loop does not require to be coarse grained. The corresponding results for the NN loop is depicted in Figs. 4(c) and 4(d). In this case, both the ImΠ α and ImΠ β suffer threshold singularities as there is no overall Källén functions coming fromÑ µν p,nl . So both the form factors have to be coarse grained after which they exactly reproduce the eB = 0 graphs.
We now turn our attention to the real part of the self energy at non-zero magnetic field and show how a numerical limit of eB → 0 agrees with the eB = 0 results. This has been shown in Fig. 5 where the real part of the form factors is shown as a function of ρ 0 invariant mass with longitudinal momentum q z = 250 MeV at temperature T = 130 MeV and at baryon chemical potential µ B = 300 MeV for the two cases eB → 0 and eB = 0. The contributions from the ππ and NN loops are shown separately. Fig. 5(a) depicts ReΠ α whereas Fig. 5(b) shows ReΠ β . As can be seen from the figure, the eB → 0 graphs exactly reproduce the eB = 0 for the case of NN loop. Whereas, for the ππ loop, eB → 0 is slightly deviated from the eB = 0 graph but with an excellent qualitative agreement in their behaviour with respect to the variation of invariant mass of ρ 0 . This small disagreement between the eB → 0 and eB = 0 graph is due to the inaccuracy in the numerical principal value integration of Eqs. (20) and (50) for which the two particle bound state threshold q 2 > 2m π = 280 MeV is less than the ρ 0 mass pole m ρ = 0.770 (in contrast, for the NN loop, the two particle bound state threshold is at q 2 > 2m N = 1.878 GeV much higher than the range of the plot).
Having checked the consistency of the non-zero magnetic field calculations, we now proceed to present the imaginary part of the self energy for nonzero values of the magnetic field. In Fig. 6, the variation of ImΠ α is shown as a function of ρ 0 invariant mass with longitudinal momentum q z = 250 MeV at temperature T = 130 MeV and at baryon chemical potential µ B = 300 MeV. We have plotted the self energy upto q 2 = 1.5 GeV for which the Unitary cut of the NN loop does not contribute. Fig. 6(a) depicts ImΠ α at magnetic field eB = 0.05 GeV 2 in which the spikes get separated from each other by a finite value and it oscillates about the eB = 0 graph. This is more clearly visible in the CG points which are used to obtain a coarse-grained interpolated (CGI) graph. Fig. 6(b) shows the CGI imaginary part at two different values of the magnetic field (eB = 0.05 and 0.10 GeV 2 respectively); both of them are found to oscillate about the eB = 0 graph. Moreover, with the increase in magnetic field, the oscillation frequency decreases with an increase in the oscillation amplitude. This behavior of the imaginary part with increasing magnetic field is consistent with Fig. 4, where for the eB → 0 case, the oscillation frequency becomes infinite and amplitude becomes zero, thus reproducing the eB = 0 graph. Also with the increase in magnetic field, the threshold of the unitary cut moves towards the higher invariant mass value as discussed in Sec. VII. This has been shown clearly in the inset plot.  The corresponding results for the ImΠ β due to ππ loop as a function of ρ 0 invariant mass with longitudinal momentum q z = 250 MeV at temperature T = 130 MeV and at baryon chemical potential µ B = 300 MeV are shown in Fig. 7 for the two different values of the magnetic field eB = 0.10 and 0.20 GeV 2 . Analogous to ImΠ α , ImΠ β also oscillates about eB = 0 curve, but in this case the oscillation frequency is much smaller as compared to ImΠ α . The threshold of the Unitary cut moves towards higher invariant mass with the increase in magnetic field as clearly depicted in the inset plot.
As discussed in Sec. VII, a non trivial Landau cut contribution in presence of external magnetic field may appear even if the loop particles have the same mass. In this case, we have observed Landau cut contribution only in ImΠ α , whereas the Landau cut does not appear in ImΠ β . This can be understood from the expressions of trace and 00 component ofÑ µν π,nl andÑ µν p,nl as given in Appendix F. It can be noticed that, for both the ππ and proton-proton loops, the expression for the trace (i.eÑ  )). This implies that, for ImΠ α , the loop particles can be in different Landau levels whereas for ImΠ β the loop particles will always stay in the same Landau levels. Thus, as discussed in Se. VII, the non-trivial Landau cuts will appear only in ImΠ α and not in ImΠ β . The contribution of the CGI Landau cuts to ImΠ α as a function of ρ 0 invariant mass with longitudinal momentum q z = 250 MeV is shown in Fig. 8. It is to be noted that, the Landau cuts also contain the threshold singularities and thus have to be coarse grained. Fig. 8(a) shows the variation of ImΠ α at a temperature T = 130 MeV and at baryon chemical potential µ B = 300 MeV for three different values of the magnetic field (eB = 0.05, 0.07 and 0.10 GeV 2 respectively), whereas Fig. 8(b) shows the corresponding variation at magnetic field (eB = 0.10 GeV 2 ) for two different values of temperature (T = 100 and 130 MeV respectively). The contributions due to ππ loop and proton-proton loops are shown separately and in Fig. 8(b); the contribution due to proton-proton loop is shown for two different values of baryon chemical potential (µ B = 200 and 300 MeV). As can be seen from the figures, the threshold of the Landau cuts due to ππ loop is different (greater) than that of proton-proton loop which can be understood from the discussions of Sec. VII. The threshold for ππ loop is q 2 < m 2 π + eB − m 2 π + 3eB , whereas the same for proton-proton loop is q 2 < m N − m 2 N + 2eB . The shift of the Landau cut threshold towards the higher invariant mass values with the increase in magnetic field can be clearly seen in Fig. 8(a). It is observed the the magnitude of the Landau cut contribution due to proton-proton loop is much less than that of ππ loop at lower values of the magnetic field and they become comparable to each other only at eB 0.10 GeV 2 . In Fig. 8(a), we observe that with the increase in temperature and density, the Landau cut contribution increases without changing its threshold in the invariant mass axis. We now turn our attention to the real part of the self energy at finite temperature under external magnetic field. In Fig. 9, we have shown the thermal contribution to the real part of the self energy as a function of invariant mass with ρ 0 longitudinal momentum q z = 250 MeV at temperature T = 130 MeV and at baryon chemical potential µ B = 300 MeV for two different values of the magnetic field (eB = 0.05 and 0.10 GeV 2 respectively). The contributions from the ππ and NN loops are summed up in this figure. We notice that, with the increase in magnetic field, the thermal contribution to the real part of the self energy oscillates about the eB = 0 curve. The oscillation frequency and the oscillation amplitude respectively decreases and increases with the magnetic field.
Next in Fig. 10, the "eB-dependent vacuum" contribution to the real part of the self energy is shown as a function of ρ 0 invariant mass with longitudinal momentum q z = 250 MeV for two different values of magnetic field (eB = 0.10 and 0.20 GeV 2 respectively). Figs. 10(a) and 10(b) show the contributions from Π α and Π β respectively. The contributions due to ππ and proton-proton loops are shown separately. First of all, we note that at eB = 0, these term will vanish. With the increase of the magnetic field, the eB-dependent vacuum term also increases and the contribution of Π β is more than Π α .
Having obtained the real and imaginary parts of the self energy, we now proceed to evaluate the in-medium spectral functions of ρ 0 under external magnetic field. We have from Eq. (106), the complete ρ 0 propagator as where the coefficients are given in Esq. (107)-(111) and the basis tensors are provided in Eqs. (72)-(75). Since we will be considering the special case q ⊥ = 0 for which Π α = Π γ and Π δ = 0 as given in Eq. (112), the coefficients in the above equation become so that the complete in-medium interacting propagator is given by It is clear from the above equation, that there will be three modes for the propagation of ρ 0 meson in magnetized medium for vanishing transverse momentum of ρ 0 . Of the three modes, two are found to be degenerate (the first and third term in the RHS of above equation) leaving two distinct modes for the propagation of ρ 0 which we denote as Mode-A and Mode-B respectively. We now define the spectral function S ρ of ρ 0 for the two distinct modes as the the imaginary part of the complete propagator which is obtained from Eq. (148) as and In Fig. 11, the spectral function for the two modes at zero magnetic field is shown as a function of ρ 0 invariant mass with ρ 0 longitudinal momentum q z = 250 MeV at baryon chemical potential µ B = 300 MeV for three different values of temperature (T = 100, 130 and 160 MeV respectively). The vacuum spectral function (which is same for the two modes) is also shown for comparison. We find that, the spectral functions have a nice Breit-Wigner shape around the ρ 0 mass pole with a width O(150 MeV) corresponding to the decay of ρ 0 → π + π − . With the increase in temperature, the width of the spectral function increases and the peak decreases. Physically, it corresponds to the enhancement of the decay process in the medium implying that the ρ 0 become more unstable at a high temperature. It is important to note that, for the invariant mass region shown in the plot, the imaginary part of the self energy that enters in the calculation of spectral function is completely due to the Unitary-I cut of ππ loop. On the other hand, the real part of the self energy that enters in the spectral function calculation has contributions from both the ππ and NN loops. We now turn on the external magnetic field and show the spectral function of ρ 0 as a function of its invariant mass for the two modes in Fig. 12. The range of the invariant mass axis is taken as 0.5-1.2 GeV which is dominated by the Unitary cut contributions from the ππ loop. In Fig. 12(a), the spectral function with ρ 0 longitudinal momentum q z = 250 MeV at temperature T = 130 MeV and at baryon chemical potential µ B = 300 MeV is shown for three different values of the magnetic field (eB = 0.10, 0.15 and 0.20 GeV 2 respectively). It is observed that, with the increase in the magnetic field, the two modes get well separated from each other and the threshold of the spectral function moves towards higher values of invariant mass corresponding to the magnetic field dependent Unitary cut threshold of the imaginary part of the self energy. At sufficiently high values of the magnetic field, the spectral function misses the ρ 0 mass pole (770 MeV) so that it looses its Breit-Wigner shape which may be termed as ρ 0 "melting" in presence of magnetic field. The critical value of the magnetic field for a given temperature and baryon chemical potential for which the ρ 0 will melt is discussed later.
In Fig. 12(b), the spectral function with ρ 0 longitudinal momentum q z = 250 MeV at a magnetic field eB = 0.10 GeV 2 and at a baryon chemical potential µ B = 300 MeV is shown for three different values of temperature (T = 100, 130 and 160 MeV respectively). In this case, the threshold of the spectral function remains fixed and for both the modes, the spectral function becomes shorter and wider with the increase in temperature with a marginal shift of its peak. The shift of the peak is due to the modification in the real part of the self energy with the change in temperature. Fig. 12(c) depicts the spectral function with ρ 0 longitudinal momentum q z = 250 MeV at a magnetic field eB = 0.10 GeV 2 and at a temperature T = 160 MeV for three different values of the baryon chemical potential (µ B = 200, 300 and 400 MeV respectively). Analogous to the previous case, the threshold of the spectral function remains fixed for both the modes. Since the baryon chemical potential only affects the real part of the self energy in the given kinematic region, the peak of the spectral function changes its position (keeping the width almost same) with the change in baryon chemical potential. It can be noticed, that in contrast to Fig. 12(b), the peak position of the spectral function is more sensitive to µ B as compared to the temperature which is due to the dominant contribution coming from NN loop.
In Fig. 12(d), the spectral function at a magnetic field eB = 0.10 GeV 2 and at a temperature T = 130 MeV with baryon chemical potential µ B = 300 MeV is shown for two different values of ρ 0 longitudinal momentum (q z = 0 and 500 MeV). In this case, the threshold of the spectral function remains same and the height of the spectral function increases with the increase of the longitudinal momentum.
We have already mentioned that, a non-trivial Landau cut in the physical kinematic region would appear in presence of the external magnetic field. In our case, the non-zero contribution to the Landau cut comes only from the form factor ImΠ α which is reflected in the the spectral function of Mode-(A). In Fig. 13, the spectral function as a function of ρ 0 invariant mass with ρ 0 longitudinal momentum q z = 250 MeV is shown in the low invariant mass region which is dominated by the Landau cut contribution. It can be observed that the magnitude of the spectral function in this region is much lower as compared to the Unitary cut regions. Fig. 13(a) shows the spectral function at temperature T = 130 MeV and at baryon chemical potential µ B = 300 MeV for three different values of magnetic field (eB = 0.10, 0.15 and 0.20 GeV 2 respectively). As can be seen in the graph, the threshold of the Landau cut moves towards the higher values of invariant mass with the increase in magnetic field as a consequence of similar behaviour of the Landau cut contribution to the imaginary part as shown in Fig. 8. Also the height of the spectral function is enhanced with the increase in eB. Fig. 13(b) shows the corresponding plots of spectral function at magnetic field eB = 0.10 GeV 2 for four different combinations of temperature and baryon chemical potential ((T = 100 MeV, µ B = 300 MeV), (T = 130 MeV, µ B = 300 MeV), (T = 160 MeV, µ B = 300 MeV) and (T = 160 MeV, µ B = 400 MeV) respectively). As can be seen in the graph, the height of the spectral function increases with the increase in temperature and density owing to an enhancement of the corresponding scattering processes in presence of external magnetic field. We now proceed to obtain the effective mass and dispersion relation of the ρ 0 in a magnetized medium. They follow from the pole of the complete ρ 0 propagator given in Eq. (148) which are obtained by solving the following transcendental equations ω 2 − q 2 z − m 2 ρ + ReΠ α (q 0 = ω, q z , eB, T, µ B ) = 0 (151) sponding to ρ 0 propagation in the magnetized medium. The effective mass m * ρ of ρ 0 is obtained from the dispersion relation by setting q z = 0 i.e. m * ρ (eB, T, µ B ) = ω(q z = 0, eB, T, µ B ). Fig. 14(a) depicts the variation of m * ρ /m ρ as a function of magnetic field at a temperature T = 130 MeV and at a baryon chemical potential µ B = 300 MeV. The effective mass for the two modes starts from the same value arround eB = 0 and with the increase in magnetic field, they get separated. For both the modes, the effective ρ 0 mass decreases with the increase in the magnetic field which is due to the strong positive contribution coming from the dominating eB-dependent vacuum part. The effect of magnetic field is found to be more in Mode-(B) as compared to Mode-(A). At a magnetic field value eB = 0.20 GeV 2 , the effective ρ 0 mass in Mode-(A) decreases by about 2% whereas for the Mode-(B) it decreases by about 10%. Fig. 14(b) depicts the corresponding variation of effective mass with temperature at a magnetic field eB = 0.10 GeV 2 and at a baryon chemical potential µ B = 300 MeV. We find that, for both the modes effective mass of ρ 0 get enhanced by a small amount with the increase in temperature. Even at T = 160 MeV the change in effective mass is less than 2%. In Fig. 14(c), the variation of effective ρ 0 mass is shown as a function of baryon chemical potential at a magnetic field eB = 0.10 GeV 2 and at a temperature T = 130 MeV. In this case also, we observe an enhancement of the effective mass for both the modes with the increase in baryon density. Though the effect of µ B on effective mass is more at a higher value of µ B the change in the effective mass remains less than 2% even at µ B = 500 MeV.
Next, we present the dispersion curves of ρ 0 propagation in magnetized medium for both the modes in Fig. 15. We have plotted the energy ω of the ρ 0 scaled with the inverse of the vacuum rho mass m ρ = 770 MeV as a function of the longitudinal momentum of ρ 0 . Fig. 15(a) depicts the dispersion curves at temperature T = 130 MeV and at baryon chemical potential µ B = 300 MeV for two different values of magnetic field (eB = 0.10 and 0.20 GeV 2 respectively). Fig. 15(b) shows the same at a magnetic field eB = 0.10 GeV 2 and baryon chemical potential µ B = 300 MeV for two different temperatures (T = 100 and 160 MeV respectively). Finally, Fig. 15(c) shows the corresponding graphs at a magnetic field eB = 0.10 GeV 2 and at a temperature T = 130 MeV for two different values of baryon chemical potential (µ B = 200 and 400 MeV respectively). In all the cases, the dispersion curves are well separated from each other at lower transverse momentum. With the increase in q z , the loop correction becomes subleading with respect to the kinetic energy of ρ 0 and thus, it approaches to a light-like dispersion.
Finally we calculate the decay width of ρ 0 for the decay into charged pions which is defined for the two modes as In Fig. 16, the variation of the decay width Γ of ρ 0 scaled with inverse of its vacuum width (Γ 0 = 156 MeV) for the two modes is shown as a function of magnetic field. Note that the vacuum decay width is obtained from the imaginary part of the vacuum self energy as Results are presented for two different combinations of temperature and baryon chemical potential ((T = 130 MeV, µ B = 300 MeV) and (T = 160 MeV, µ B = 400 MeV) respectively). Because of the presence of threshold singularity in ImΠ α , Γ (A) also suffers from the presence of threshold singularity for which it needs to be coarse grained. However, ImΠ β and hence Γ (B) is finite and free from the singularities. As can be seen from the figure, the ratio Γ/Γ 0 starts from a value greater than unity near eB = 0 which is due to the enhancement of the decay width over its vacuum value due to the effect of finite temperature and density. Also for a particular value of magnetic field, larger decay width is observed at higher temperature and density. Near eB = 0, the two modes have almost the same decay widths which begin to differ from each other with the increase in the magnetic field. An oscillatory behaviour of the decay width can be clearly seen throughout the magnetic field range. One should also notice that, for both the modes, the oscillation amplitude increases whereas oscillation frequency decreases with eB. Finally at a critical value of the magnetic field, the decay width becomes zero. This is because of fact that, the eB-dependent Unitary cut threshold for the ππ loop has to satisfy m * ρ (eB) > 2 m 2 π + eB (156) for a kinematically favorable decay of ρ 0 → π + π − . But, with the increase in magnetic field, the RHS of the above equation increases, whereas m * ρ in the LHS decreases so that at some critical value of magnetic field, the above inequality is violated and the decay width becomes zero. Physically it means, that ρ 0 becomes stable against the decay into π + π − pair. This critical value of the field may be considered as the critical value of the magnetic field required for the "melting" of the spectral function of ρ 0 .
In order to calculate the critical value of the magnetic field eB c for a given temperature T and baryon chemical potential µ B , we need to solve the transcendental equation The green dash-dotted curve in Fig. 14(a) corresponds to m * ρ /m ρ = 2 m 2 π + eB so that, the intersection of this curve with the m * ρ = m * ρ (eB) represents the solution of the above equation. In Fig. 17, we show the variation of the critical magnetic field eB c for the two decay modes. Fig. 17(a) depicts eB c as a function of temperature for two different values of baryon chemical potential (µ B = 50 and 200 MeV) whereas Fig. 17(b) shows the corresponding variation with baryon chemical potential at two different values of temperature (T = 100 and 160 MeV). Although, with fixed temperature, the variation with respect to µ B shows monotonically increasing trend, both the plots suggests non-monotonic variations of the critical magnetic field with respect to the temperature. More specifically, there exists a maximum value of chemical potential (see Fig. 17(b)) below which the critical field decreases with the temperature there by requiring relatively weaker magnetic field to completely stop the particular decay channel. However, for even larger values of µ B , a significant increase with temperature can be observed for both of the decay modes.

IX. SUMMARY AND CONCLUSIONS
In this work, the spectral properties of the neutral rho meson is studied at finite temperature and density in a constant external magnetic field using the real time formalism of finite temperature field theory. The effective ρππ and ρN N interactions are considered for the evaluation of the one loop self energy of ρ 0 . Accordingly, the magnetically modified in-medium propagators for pions and protons are used which contain infinite sum over the Landau levels implying no constraint on the strength of the external magnetic field. From the self-energy, the eB-dependent vacuum part is extracted by means of dimensional regularization in which the ultraviolet divergence corresponding to the pure vacuum self energy is isolated as the pole of gamma and Hurwitz zeta functions. It is shown that the external magnetic field does not create additional divergences so that the vacuum counter terms required in absence of the background field remain sufficient to renormalize the theory at non zero magnetic field.
The general Lorentz structure for the in-medium massive vector boson self energy in presence of external magnetic field has been constructed with four linearly independent basis tensors out of which three form a mutually orthogonal set. Thus, the the extraction of the form factors from the self energy becomes considerably simple. Moreover, it is shown that with vanishing perpendicular momentum of the external particle, one can arrive at new set of constraint relations among the form factors which essentially leave only two form factors to be determined from the self energy. As a consistency check, the numerical B → 0 limit of the real as well as imaginary parts of the form factors are shown to reproduce the zero field results. Solving the Dyson-Schwinger equation with the one loop self energy, the complete interacting ρ 0 propagator is obtained. Consequently, two distinct modes are observed in the study of the effective mass, dispersion relations and the spectral function of ρ 0 where one of the modes ( Mode-A) possesses two fold degeneracy. It is known [31,38] that non trivial Landau cuts appear in presence of external magnetic field along with finite temperature even if the loop particles are of equal mass which is completely a magnetic field effect. However, in contrast to Mode-A, the non-trivial Landau cut is found to be absent in case of Mode-B. Also, sharper decrease in the effective mass is observed for the later which essentially stems from the dominant eB-dependent vacuum contribution in the real part of the corresponding form factor.
Finally, the decay width for ρ 0 → π + π − channel is obtained for the two distinct modes and is found to become zero at certain critical values of magnetic field depending upon the temperature and baryon chemical potential. The corresponding variation of the critical field with these external parameters shows increasing trend for large baryonic chemical potential. However, it is observed that, both the distinct modes possess a maximum value of µ B below which the temperature dependence gets reversed. Especially, at a given temperature (say T = 160 MeV) , eB c attains the lowest values (123 MeV 2 for Mode-A and 116 MeV 2 for Mode-B) in case of zero chemical potential. In Ref. [27], charged rho meson condensation has been studied at finite temperature and density. For charged rho mesons, the critical field for which the vector meson mass vanishes is observed to lie in the range of 0.2-0.6 GeV 2 at zero density with temperature in the range 0.2-0.5GeV. However, in case of ρ 0 , the absence of the trivial Landau shift in the energy eigenvalue results in much slower decrease in the effective mass. As a consequence, unrealistically high magnetic field values are required to observe neutral rho condensation in presence of temperature and medium (see Fig.14). In this scenario, the suppression in the ρ 0 → π + π − channel can serve as an important observable. It should be mentioned here that magnetic modification of rho meson properties studied in this work deals with effective hadronic interactions. Thus, the observable modification can only occur if the initial burst of magnetic field survives up to hadronization retaining an appreciable field strength. The requirements can be ensured only after a comprehensive knowledge of the space time dependence of the produced magnetic field is obtained which definitely demands a significant amount of research work in future.
In order to evaluate the momentum integrals in Eqs. (3) and (4), they are rewritten as where, N µν N (q, k) contains the trace over Dirac matrices: Applying standard Feynman paramerization, the denominators of Eqs. (B1) and (B2) are combined to get, where, and the space-time dimension has been changed from 4 to d in order to work with the dimensional regularization so that the additional scale parameters Λ π and Λ N of dimension GeV 2 have been introduced to keep the overall dimension of the self energy same. It is now straight forward to perform the momentum integrals of the above equations after a momentum shift k → (k − xq) using the identities provided in Appendix A, so that, the vacuum self energies becomes (Π µν π ) pure-vac (q) = (q 2 g µν − q µ q ν ) where ε = (2 − d/2). Expanding the above equations about ε = 0, we get (Π µν π ) pure-vac (q) = (q 2 g µν − q µ q ν ) −g 2 ρππ q 2 32π 2 1 0 dx∆ π 1 ε − γ E + 1 − ln ∆ π 4πΛ π ε→0 (B10) (Π µν N ) pure-vac (q) = (q 2 g µν − q µ q ν ) where, γ E is the Euler-Mascheroni constant.