Neutrino propagation in an electron background with an inhomogeneous magnetic field

We study the electromagnetic coupling of a neutrino that propagates in a two-stream electron background medium. Specifically, we calculate the electromagnetic vertex function for a medium that consists of a normal electron background plus another electron stream background that is moving with a velocity four-vector vμ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$v^\mu $$\end{document} relative to the normal background. The results can be used as the basis for studying the neutrino electromagnetic properties and various processes in such a medium. As an application, we calculate the neutrino dispersion relation in the presence of an external magnetic field (B\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf {B}$$\end{document}), focused in the case in which B is inhomogeneous, keeping only the terms of the lowest order in 1/mW2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$1/m^2_W$$\end{document} and linear in the B and its gradient. We show that the dispersion relation contains additional anisotropic terms involving the derivatives of B\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf {B}$$\end{document}, such as the gradient of k^·(v×B)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{k}}\cdot (\mathbf {v}\times \mathbf {B})$$\end{document}, which involve the stream background velocity, and a term of the form k^·(∇×B)\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{k}}\cdot (\nabla \times \mathbf {B})$$\end{document} that can be present in the absence of the stream background, in addition to a term of the form k^·v\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{k}}\cdot \mathbf {v}$$\end{document} and the well known term k^·B\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$${\hat{k}}\cdot \mathbf {B}$$\end{document} that arises in the constant B\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathbf {B}$$\end{document} case. The derivative-dependent terms are even under a CP transformation. As a result, in contrast to the latter two just mentioned, they depend on the sum of the particle and antiparticle densities and therefore can be non-zero in a CP-symmetric medium in which the particle and antiparticle densities are equal.


I. INTRODUCTION AND SUMMARY
The subject matter of this paper is the propagation of neutrinos in a medium in the presence of an external electromagnetic field.There are various problems of interest associated with this subject that have been well studied in previous works.In most previous studies the interest has been on a medium consisting of a thermal background of various particle species, which can be taken to be at rest, in the presence of an external magnetic field in the same frame, which is assumed to be homogeneous.In those studies typically the focus is on the dispersion relation of a neutrino that propagates in such environments.The assumptions underlying the previous works do not allow us to consider situations in which the thermal backgrounds of the different particle species move with some relative velocity relative to each other, and/or situations in which the external field is not homogeneous.
There are several reasons why considering the more general situations just mentioned above are of interest.For example, the propagation of photons in two-stream plasma systems is a well studied subject in the context of plasma physics, particularly with regard to the so-called two-stream instabilities [1].In recent works similar studies have been carried out for magnetized two-stream plasma systems [2].In these works the focus is typically the dispersion relation of the photon when it propagates in the environment that is being considered.
From a fundamental and conceptual point of view the problem we want to consider is the counterpart for neutrinos.The problem of the propagation of neutrinos in magnetized media is relevant in several physical contexts, such as pulsars [3], supernovas [4] and gamma-ray bursts [5], where the magnetic fields are believed to have important implications.In those contexts, the effects of stream backgrounds and/or inhomogeneous fields can be of practical interest.
In a recent work [6] we initiated the study of the propagation of neutrinos in medium along these lines, calculating the self-energy and dispersion relation of a neutrino that propagates in a magnetized two-stream background medium.Specifically, we considered a medium composed of an electron background, which can be taken to be at rest, and a second electron background that moves with a velocity four-vector v µ relative to the first.We refer to them as the normal and stream backgrounds, respectively.In addition we assumed that, in the rest frame of the normal background, there is a magnetic field ( B) that is homogeneous.The calculation was based on the local limit of the weak interactions, and therefore restricted to the leading O(1/m 2 W ) terms, and on the application of the Schwinger propagator method, adapted to the two-stream background, but keeping only up to the linear terms in B.
The main results obtained in Ref. [6] are summarized as follows.For a neutrino that propagates in a two-electron background and a constant magnetic field, as described above, the dispersion relation acquires an anisotropic contribution of the form k • v, in addition to the well known term k • B [7,8].As discussed and explained in Ref. [6], a term of the form k • ( v × B) does not appear in the constant B case.The physical reason behind this result is that such a term is odd under time-reversal and there is no source of time-reversal breaking effects in the context of our calculations.However it was noted that terms of similar form, but involving the derivative of the electromagnetic field, are even under time-reversal and therefore could be present in the case that the electromagnetic field is not homogeneous.
As a continuation of the above work, here we study the electromagnetic vertex function of a neutrino that propagates in a two-stream electron background.This complements and extends our previous work in at least two ways.On one hand, the knowledge of the vertex function allows us to determine the neutrino electromagnetic properties and to calculate the rate for various processes involving neutrinos in such media, in analogy with the study of electromagnetic neutrino processes in ordinary media [9].On the other hand, which we pursue here, by considering the effective neutrino interaction with an external electromagnetic field, it is used to determine the self-energy and dispersion relation of a neutrino that propagates in a two-stream medium with an inhomogeneous magnetic field.This complements our previous calculation of the dispersion relation based on the Schwinger propagator method, which is restricted to a uniform magnetic field.The dispersion relation obtained can be used as the basis for studying the effects of inhomogeneous fields on the neutrino oscillations in several environments such as pulsars, supernovas, and gamma-ray bursts that have been considered in the literature cited, as well as several related application contexts where the inhomogeneity of the magnetic fields may have a prominent role such as transition radiation induced by a magnetic field [10], neutrino-induced plasma instabilities in supernova [11] and the effects of asymmetric neutrino propagation in proto-neutrons stars [12].
It is appropriate to mention here that the calculation of the neutrino electromagnetic vertex function in the two-stream electron background is based on the local limit of the weak interactions, i.e., it is limited to the O(1/m 2 W ) contributions.Moreover, in the application to the calculation of neutrino self-energy and dispersion relation mentioned above we retain only the terms that are at most linear in the derivatives of the field.
The results of the calculations confirm that in the case of an inhomogeneous field, the dispersion relation acquires additional terms, such as the gradient of k • ( v × B) as already mentioned above, that involve the derivatives of the magnetic field.Moreover, these additional terms are even under the CP transformation and as a result they are proportional to the sum of the particle and antiparticle densities.This is in contrast with the k • v and k • B terms, which are even under CP and depend on the difference of the particle and antiparticle densities.In situations where the medium is CP -symmetric and the particle and antiparticle densities are equal, the O(1/m 2 W ) the contributions from the k • v and k • B terms vanish, and the contributions from the the terms involving the derivatives of the magnetic field could gain more importance, depending on the degree of inhomogeneity of the magnetic field.
It is worth mentioning that, in order to calculate properly the stream contribution to the vertex function, and more specifically the vertex function's zero photon momentum limit (which is related to the neutrino index of refraction), the screening effects of the electron background must be taken into account.The technical reason is that the electric form factors (those that couple to the electric components of the electromagnetic field ∼ k•( v× B)) diverge in the zero photon-momentum limit, unless the screening effects are taken into account.In the case in which there is only one background ( v = 0), or the magnetic field is uniform, only the magnetic form factor couplings enter in the effective interaction with the electromagnetic field, for which the screening corrections are not relevant.An important ingredient of the present work is the proper inclusion of the background screening effects that are present in the kind of medium that we are envisaging, in the calculation of the neutrino index of refraction.
In Section II A we summarize some of the notations and conventions that are used throughout.The 1-loop formulas for the vertex function are given in Section II B, generalizing the formulas given in Ref. [7], adapted to the present notation and context.As already mentioned, they are based on the local limit of the weak interactions, i.e., they are restricted to the O(1/m 2 W ) contributions.The vertex function is expressed in terms of a set of form factors that are given as integrals over the distribution functions of the background electrons.Since the calculation of the self-energy and dispersion relation in a non-homogeneous external field involves evaluating the vertex function in the static limit appropriately, specially in the context of the two-stream system, in Section II C we define precisely this limiting operation.There we also summarize the static limit value of the integrals involved in such formulas, which are relevant in the calculation of the self-energy and dispersion relation.Some of the calculational details regarding those formulas are provided in Appendix B. The actual calculation of the self-energy in the presence of an external inhomogeneous field is carried out in Section III, retaining the terms that are at most linear in the derivatives of the field, and paying attention to the treatment and incorporation of the screening effects of the electron background.There we first enumerate the possible terms that may appear in the B-dependent part of the self-energy under the specified conditions, and write down its generic form in terms of a set of structure tensors with corresponding coefficients to be determined.The coefficients are then identified by considering the transition amplitude in the presence of an external field, using the results of the one-loop expression for the neutrino electromagnetic vertex function.The need to include the screening effects for properly determining the self-energy in the two-stream background case is explained there in more detail.The corresponding dispersion relations are obtained and discussed in Section IV, focusing on some of the features that illustrate the salient implications of the results for the self-energy.In Section V we review our work and summarize the main results.

A. Notation and conventions
We borrow some of the notation from Ref. [6], which we briefly summarize here for convenience as follows.We use the symbols e and e ′ to refer to the electrons in the normal and stream backgrounds, respectively, while the symbol f stands for either e or e ′ .Denoting by u µ f the velocity four-vector of each background, the convention stated above means that the velocity four vector of the normal background is set to where, as usual, while for the stream with (2.4) The relevant diagrams for the calculation of the electron background contributions to the neutrino electromagnetic vertex are shown in Fig. 1.For the calculation we will need the following neutral current couplings where We denote the momentum vectors of the incoming and outgoing neutrino by k, k ′ , respectively, and denotes the momentum vector of the incoming photon [13].
The form factors of each background contribution are functions of the scalar variables Physically, Ω f represents the energy of the photon in the rest frame of the normal background, while Q f is the magnitude of the 3-momentum vector in the same frame, which we denote by Q f .

B. 1-loop formulas
As already mentioned, the relevant diagrams for the calculation of the electron background contributions to the neutrino electromagnetic vertex function are shown in Fig. 1.
FIG. 1.The diagrams that contribute to the neutrino electromagnetic vertex in an electron background to the lowest order for a given neutrino flavor ν ℓ (ℓ = e, µ, τ ).Diagram (a) contributes only to the ν e vertex function, while Diagram (b) contributes for the three neutrino flavors.
We denote by Γ (W,Z) f µ the contribution from diagrams (a) and (b), respectively, and write the total vertex function as (2.9) where We now rely on the results obtained in Ref. [7], adapted for our present purposes.The results of the one-loop calculation of Γ is that Γ f µ can be written in the form where the tensors T f µν do not contain any gamma matrices and have the decomposition with T eT,L,P and T e ′ T,L,P being scalar functions of Ω e , Q e and Ω e ′ , Q e ′ , respectively.In writing the last term in Eq. (2.12) we have taken into account the fact that the definition of q here [Eq.(2.7)] is the opposite to that in Ref. [7], as already mentioned in Section II A. The basis tensors R f µν , Q f µν , P f µν that appear in Eq. (2.12) are defined by where gµν = g µν − q µ q ν q 2 , (2.14) and ũfµ = gµν u ν f . (2.15) The tensors R f µν , Q f µν , P f µν satisfy the relations while the contractions of anyone of them with the others vanish.
The functions T f T,L,P that appear in Eq. (2.12) are given by where with In these formulas, m is the electron mass, and f f, f are the electron and positron thermal distribution functions From Eq. (2.17), X e (ν µ,τ ) −Y e (ν µ,τ ) (2.23)

C. Static limit
As we have mentioned, the results for the electromagnetic vertex function will be used as the starting point to determine the self-energy and dispersion relation in an external field, which involves evaluating the vertex function in static limit.It is appropriate to state precisely what we mean by this limit, specially in the context of our calculation that includes the effects of the stream background and possibly a non-homogeneous external field.
Let us look first at the case considered in Ref. [7], namely the normal electron background contribution to the neutrino index of refraction in the presence of an external constant B field, that is a field that is constant in time and homogeneous in space.This requires the evaluation of the normal background contribution to the vertex function in the zero momentum limit, which operationally is implemented by first setting maintaining Q e fixed, and then taking the limit We indicate this two-step process by the notation The idealization involved here is that the νν transition amplitude is being calculated over a region that is microscopically large, but macroscopically sufficiently small such that the external field is constant over the region.In this situation, the terms in the νν transition amplitude that contain factors of q multiplying the external field do not contribute.
In the present work we consider the possibility that the external field is not necessarily homogeneous.This is taken into account by interpreting each factor of q µ multiplying the external field in the νν amplitude as a derivative acting on the external field.
In addition, in the presence of the stream, the limit Q e → 0 is complicated by the fact that the stream contributions to the neutrino electromagnetic vertex function depend on the variables where Ω e ′ represents the energy of the photon in the rest frame of the stream background while Q e ′ is the magnitude of the 3-momentum vector in that frame.
For Ω e = 0, they are given by (2.29) Therefore, there is a separate dependence on u e ′ • Q e , and not just on the magnitude Q e , and as a consequence the process of taking the zero momentum limit Q e → 0 is not unique.
For our purposes (calculating the self-energy in the presence of an external field), we take the zero momentum limit in this case in the following manner.First, after setting Ω e = 0, make an expansion of the stream contribution to the vertex function in powers of u e ′ • Q e , and then in harmony with the treatment of the terms with q µ specified above in Eq. (2.27), interpret each such factor as a derivative acting on the external field.Since the remaining factors then depend only on Q e after making this replacement, the Q e → 0 limit can be taken subsequently in an unambiguous way.In particular, the stream contribution form factors, which are functions of Ω e ′ and Q e ′ , are evaluated in this limit according to a rule analogous to Eq. (2.26), that is In this work we retain the terms that are at most linear in the derivatives acting on the external magnetic field after making the identifications made in Eqs.(2.27) and (2.30).In the idealized situation that the external field is strictly homogeneous all such terms vanish.
For easy reference we quote here the following formulas that are given in Eqs.(2.27-2.29)and Eq.(3.2) of Ref. [7], where The derivation of the above formulas is sketched in Appendix A, and in Appendix B we derive the formula for the static limit value of where which will be relevant in the discussion in Section III.
The integrals defined in Eqs.(2.33) and (2.35) can be performed straightforwardly once the distribution functions are specified.For guidance and reference purposes we give below the results of their evaluation in the particular case that the distribution functions can be taken to be those of the classical ideal gas.Using the fact that in that the classical (independently of whether the gas is relativistic or not), it follows simply that (2.37) In the case of C 0 and A ′ 0 the results in the relativistic and non-relativistic cases are different.
In the non-relativistic limit (β f m ≫ 1) while in the extremely relativistic limit (β f m ≪ 1) The derivation of Eqs.(2.38) and (2.39) is sketched in Appendix C.

III. NEUTRINO SELF-ENERGY IN A STATIC MAGNETIC FIELD A. General form
The chirality of the neutrino interactions implies that the background contribution to the neutrino self-energy, Σ eff , is of the form and we will decompose Σ in the form where Σ (B) stands for part that depends on B and Σ (m) for the B-independent part.The neutrino dispersion relation is obtained by solving the equation As is well known, in the two-stream electron background Σ (m) is of the form where, to order 1/m 2 W , The coefficients a f are O(1/m 4 W ) and therefore we will discard them.The issue that we address now is the enumeration of the possible terms that can appear in Σ (B) for the two-stream background with the external electromagnetic field.The situation we consider is that in the rest frame of the normal background there is a magnetic field B = B b, and in that frame we define We can then write the corresponding EM tensor in the form and its dual, as usual, is given by In the enumeration of the possible terms that may appear in the result of the 1-loop calculation of Σ (B) , we must remember the following working conditions: 1. restrict ourselves to the terms that are most linear in the derivatives of the field; 2. omit the terms that depend on the neutrino momentum k since they would be of O(1/m 4 W ) that we are not considering; 3. in the 1-loop calculation each background contributes independently, so that terms involving the products of vectors u µ f corresponding to different backgrounds do not appear.
The following is then the list of the terms that can appear: (a) Terms with no derivatives of the electromagnetic field: F µν u f ν γ µ (b) Terms with one derivative of the electromagnetic field: Schematic diagrams for the effective neutrino electromagnetic vertex taking into account the polarization effects of the background electrons as expressed in Eq. (3.13).
The coefficients defined here will be determined by calculating the νν transition amplitude in the presence of an external electromagnetic field, using the 1-loop formulas for the vertex function given in Section II.

B. νν transition amplitude in an external electromagnetic field
We are now set to consider the νν transition amplitude in the presence of an external electromagnetic field.The external electromagnetic potential is represented by and the corresponding field by where The diagram for the process is shown in Fig. 2. As shown schematically in that figure, it includes the photon polarization tensor in order to take into account the screening effects of the background electrons.The off-shell ν − ν scattering amplitude in the presence the external electromagnetic potential is then given by where Γ µ is the total neutrino electromagnetic vertex function given in Eq. (2.9) and [14] with ∆ µν ef f (q) being the photon propagator in the medium.The latter quantity is determined by solving where π µν is the two-background contribution to the photon polarization tensor.Denoting by π eµν and π e ′ µν the contributions due to the normal and stream backgrounds, respectively, then Each term in the previous relation can be decomposed according to where the photon self-energy functions π f T,L are given by with A f , B f being the integrals defined in Eq. (2.19) [15].
Let us look at the case considered in Ref. [7], an external magnetic field and only the normal background.In this case and Eq. (3.20) follows from the fact that for a pure magnetic field ũ • a = 0, which can be seen in various ways.For example, with the usual convention in which A µ = (0, A) with ∇• A = 0, it follows that u • a = 0 and q • a = 0, and therefore ũ • a = 0.More generally, ũ • a = i q 2 q µ u ν f µν , while u ν f µν = 0 when f µν corresponds to a magnetic field.Thus, remembering that π eT → 0 (and T eT → 0) in the static limit, Eqs.(3.19) and (3.20) imply that D λµ a µ → a λ in the static limit, and therefore the screening corrections are not relevant in that case.
With the stream contributions the situation is different.The stream electrons, in their own rest frame, "see" an electric field in addition to the magnetic field, and the screening corrections are relevant in that case.The present situation is complicated by the fact that in the presence of the two backgrounds the inversion required in Eq. (3.15) is not as simple in the general case as the one leading to the one-background result given in Eq. (3.19).
We overcome this difficulty here as follows.In the limit u e ′ → 0, each of the tensors R e ′ , Q e ′ , P e ′ coincides with its corresponding counterpart R e , Q e , P e .It is straightforward to show that, in general, the corresponding primed and unprimed tensors differ by terms with analogous relations for Q e ′ µν and P e ′ µν .Since, as stated in Section I, in this work we will retain only up to the linear terms u e ′ • Q e in the calculation of the self-energy, we then can then write where We reiterate that Eq. (3.22) is valid assuming that we are dropping the terms proportional to ( u e ′ • Q e ) 2 and higher powers, which translate to terms containing second and higher derivatives of the external field, by making the identification shown in Eq. (2.30).Using Eq.
(2.12) and the multiplication rules of the tensors R, Q, P , then where with TfT,P = q 2 T f T,P q 2 − π T , Using Eqs.(3.13) and (3.24), the νν amplitude is then given by An equivalent expression for the functions Γ f µ , which is more convenient for the purpose of the interpretation of the form factors and for taking the static limit, is [7] Γ where, using the formulas given in Eq. (2.13) for the tensors R f µν , Q f µν , P f µν , It follows from Eq. (3.28) that where f µν is defined in Eq. (3.12) and the M f µν are given by and from Eq. (3.27) Eq. (3.32) is a useful starting point to determine the contribution to the neutrino self-energy in a static field, including the case of an inhomogeneous field, which we consider next.

C. Self-energy
We now consider the νν transition amplitude for the case of a static magnetic field F µν .
As stated in Section II C, we make the idealization that we are calculating it over a region that is microscopically large but macroscopically sufficiently small such that the external field and its space derivatives are constant over that region.In addition we retain only the terms that are most linear in the derivatives.Operationally this means that in Eq. (3.32) we can take while neglecting the terms with higher powers of q, which would translate to terms with higher order derivatives of the external field.Then, Eq. (3.32) becomes where is obtained from Eq. (3.31) by keeping the terms that are at most linear in ( u e ′ • Q e ) and/or q and following the procedure outlined in Section II C: make the identification stated in Eqs.(2.27) and (2.30), and then take the q → 0 limit in the remaining terms as indicated in Eqs.(2.26) and (2.31).The B-dependent part of the self-energy, Σ (B) , which is identified by writing is then given by where where the coefficients t f T,L,P are defined by X e (ν µ,τ ) X e (ν µ,τ ) where A 0 f , C 0 f , A ′ 0 f are the integrals defined in Eqs.(2.33) and (2.35).Eqs.(3.38) and (3.40) summarize the result of our calculation of the contribution to the neutrino self-energy due to its interaction with an external electromagnetic field that is not necessarily homogeneous.The result given in Eq. (3.38) corresponds to the general form given in Eq. (3.9), with the coefficients given specifically by (3.41)

IV. DISPERSION RELATIONS
For the purpose of determining the dispersion relation, it is convenient to express the total self-energy, Eq. (3.2), in the form with The formula for the V µ f follows from Eqs. (3.4) and (3.9), which we summarize in the form where is non-zero only if the field is inhomogeneous.In writing Eq. ( 4.4) we have dropped the terms that vanish according to the results we have obtained in Eq. (3.41).We can express In the expression for V which in the rest frame of the matter background have components where The equation for the propagating modes, Eq. (3.3), implies that the dispersion relations are given by Thus to the lowest order in 1/m 2 W considered in this work, which among other things implies that V µ does not depend on k, and the solutions are k 0 = ω ± ( k), where (4.10) In Eq. (4.10) k denotes the unit vector along the direction of propagation.The neutrino and antineutrino dispersion relations, which are identified in the usual way as, are then given by where According to the decomposition in Eq. ( 4.3), we can write and from Eq. (4.5), For the two-stream electron background specifically, using Eqs.(2.2)-(2.4)and Eq.(3.41), and where In this case, Together with Eq. (4.12), Eqs.(4.14) and (4.16) provide a general and concise expression for the neutrino and antineutrino dispersion relations in the kind of situation that we are envisaging, and Eqs.(4.17), (4.18) and (4.20) for the two-stream electron background we are specifically considering.It is useful to consider some special situations that illustrate particular features of the general results.

A. Homogeneous external field
It is simple to verify that the results obtained in Ref. [6] are reproduced as a special case.
Thus, if the external field is homogeneous, which is the result obtained in Ref. [6].In particular, in the absence of the stream background, which leads by Eq. (4.12) to the result obtained in Ref. [7] for the dispersion relation in a magnetized electron background.The angular asymmetry of the dispersion relation in this case has been subject of much interest in the literature in connection to the problem of pulsar kick and related issues.The terms in Eq. (4.17) due to the stream background (in the same case of a homogeneous field) give an additional asymmetric contribution that depends on the direction of propagation relative to the stream background velocity.

B. Inhomogeneous magnetic field
In the case of a non-homogeneous field, the additional terms given by δ e ′ can be present.As an example, let us consider the case in which ∇ × B = 0.In this case, only the term involving h e ′ 3 in Eq. (4.18) survives and therefore, from Eq. (4.20), where n e ′ has been defined in Eq. (4.19).In particular, in addition to the angular dependence involving the direction of propagation relative to B and v, there is a dependence involving a third vector v × n e ′ .On the other hand, if the conditions are such that t e ′ L terms in Eq.
(4.20) are negligible, then C. Discussion One point that stands out in these results is the following.To the order 1/m 2 W , the coefficients b e,e ′ , c e,e ′ are proportional to the electron-positron asymmetry in the corresponding backgrounds.Therefore in a CP -symmetric medium b e,e ′ and c e,e ′ vanish to the order 1/m 2 W , and in such cases the order 1/m 4 W contributions must be included.In contrast, to the order 1/m 2 W , the t e ′ L,T terms in δ e ′ are proportional to the sum of the electron-positron number densities (in the stream background), and whence need not be zero even in a CP -symmetric medium.Thus, in a CP -symmetric medium (e.g., the Early Universe) the dominant contribution to the neutrino index of refraction could be due to the terms contained in δ where and similarly for b e ′ .Eq. (4.25) holds for the case the electron gas can be adequately described by a classical distribution in the relativistic limit (T e ≫ m).The question we ask is under what conditions the term t e ′ L k • ( v × n e ′ ) could be as important, or more important, than the b (4)  e term in the dispersion relation.Under the example idealized conditions that we have assumed for the purpose of this discussion (classical relativistic electron backgrounds), using Eqs.(2.37) and (3.40) the question translates to see for what parameters it is possible where we have put E ν ∼ | k|.Taking the quantity in parenthesis to be O (1) and N e ∼ T 3 e , the condition would be that As this particular case illustrates, the contributions to the dispersion relation due to the inhomogeneity of the magnetic field do not seem to be significant if the background is electron-positron asymmetric.
However, more generally, it is not inconceivable that those contributions may be relevant under the appropriate environmental conditions including an electron-positron symmetric background.While we have not considered a specific application, the example estimates above show that the gradient-dependent contributions could be comparable to the standard terms in such environments.Since they give rise to distinctive kinematic features in the dispersion relation (e.g., angular asymmetries) the possible need to include their effects in some specific application contexts should be kept in mind.

V. CONCLUSIONS AND OUTLOOK
In this work we have calculated the electromagnetic vertex function of a neutrino that propagates in a medium consisting of a normal electron background plus another electron stream background which is moving with a velocity relative to the normal background.The results obtained were used to determine the neutrino self-energy and dispersion relation in such a medium in the presence of an external magnetic field ( B), paying special attention to the case in which B is inhomogeneous, keeping the terms that are linear in B and its spatial derivatives.
The 1-loop formulas for the vertex function were given in Section II.The calculation is based on the local limit of the weak interactions, that is, it is restricted to the order 1/m 2 W terms only.The formulas generalize those given in Ref. [7], adapted to the present context.
The vertex function is expressed in terms of a set of form factors that are given as integrals over the distribution functions of the background electrons.We also summarized the static limit value of the integrals involved in such formulas, which are subsequently used in the calculation of the neutrino self-energy and dispersion relation in the two-stream medium in the presence of a non-homogeneous external field.
In Section III we used the results for the vertex function in the two-stream medium to determine neutrino self-energy in the presence of a static external magnetic field.In contrast to the previous calculations of the neutrino index of refraction in magnetized media, we took into account and emphasized the case in which the field is not homogeneous.There The corresponding dispersion relations were obtained and discussed in Section IV, focusing on some of the features that depend on the presence of the stream electron background and/or the inhomogeneity of the B field.In particular, we showed that the terms that depend on the field derivatives, in contrast to those that depend on B itself, are proportional to the sum of the electron and positron densities, and therefore are non-zero to order 1/m 2 W in a CP -symmetric medium in which the particle and antiparticle densities are equal.Thus, in a CP -symmetric medium the dominant contribution to the neutrino index of refraction could be due to the terms that depend on the derivatives of B, which are of order 1/m 2 W , in contrast with the constant B terms which to that order vanish and are of order 1/m 4 W in that case.
From a more general point of view, the present work is a step in our effort to study problems related to the propagation of neutrinos in a medium that consists of various particle backgrounds that may be streaming with different velocities.The results of our first step in this direction were presented in Ref. [6], in which we considered the propagation of a neutrino in a magnetized two-stream electron background medium.There we considered the self-energy and dispersion relation of neutrino that propagates in a two-stream electron medium, that is a medium composed of an electron background taken to be at rest (to which we refer as the normal background), and a second electron background that moves with some velocity relative to the first.In addition we assumed that, in the rest frame of the normal background, there is a magnetic field that is homogeneous.Here we have extended that work by considering the neutrino electromagnetic vertex function in the two-stream electron medium.As already mentioned in the Introduction, the knowledge of the vertex function allows us to determine the neutrino electromagnetic properties and to calculate the rate for various processes involving neutrinos in such media, but we do not consider these applications here.Alternatively, by considering the effective neutrino interaction with an external electromagnetic field, we have used the results for the vertex function to determine the self-energy and dispersion relation of a neutrino that propagates in a two-stream medium with an inhomogeneous magnetic field.In particular this extends the previous works on neutrino propagation in magnetized media which are restricted to the case of a homogeneous magnetic field.There is an extensive literature related to the effects of an external magnetic field in the propagation of neutrinos in dense media in a variety of astrophysical and cosmological contexts.The results of this work provide a firm setting for exploring the effects of the combined presence of stream backgrounds and inhomogeneities of external fields along the same lines.
where we have defined and Since we are interested in eventually taking the limit Q f → 0, we expand where in the second line we have used the fact that N depends on E (and not explicitly on P ) so that Substituting Eq. (A7) in Eq. (A4) we obtain Writing m 2 = E 2 − P 2 and using Eq.(A8) once more, the term in curly brackets above can be rewritten, and since the second term on the right-hand side integrates to zero we thus finally obtain the result quoted in Eq. (2.32).The results for B f and C f quoted in Eq. (2.32) can be obtained in similar fashion.In particular, since A f and B f are equal at Appendix B: Evaluation of the integral A ′ f in the static limit Here we calculate the static limit value of A ′ f including the O(Q 2 f ) term, which in turn determines T f T (0, Q f ) and in the end t f T (0, Q f ) as defined in Eq. (3.40).We wish to state explicitly here that the momentum thermal distribution functions are assumed to be isotropic (in each background's rest frame), so that the symmetric integration replacements such as that in Eq. (B14) below are valid.Thus the result given in Eq. (B20) [which is quoted in Eq. (2.34)] is subject to this restriction.From Eq. (2.19), and making the change of variable indicated in Eq. (A3), = − 1 2 where we have defined and we have used the definitions given in Eqs.(A5) and (A6) once more.In this case we expand N ∓ up to the cubic terms in Q f , and obtain Ã′ = d 3 P (2π) 3 where Using Eq. (A8), and therefore, from Eq. (B2), where have used the fact that 3N + P • ∇N = ∇ • ( P N) , (B9) which integrates to zero and therefore does not contribute to A ′ .
For the evaluation of I (2) , we first write all the derivatives with respect to P i in terms of derivatives with respect to E using Eq.(A8), which we express in the form ∇ i P N = P i N ′ , ∇ i P ∇ j P N = δ ij N ′ + P i P j N ′′ , ∇ i P ∇ j P ∇ k P N = (P i δ jk + P j δ ik + P k δ ij )N ′′ + P i P j P k N ′′′ , (B10) where we define and so on.Using Eq. (B10), the first term in the integral I (2) in Eq. (B6) is reduced using while for the second term we use 1 The second line in Eq. (B13) indicates the replacement that can be made in the integral over P , which allows us to replace in the integrand.Thus I (2) becomes This integral can be simplified further by using the following identities where the arrow symbol means that the relations hold under the integral sign after dropping the terms that integrate to zero.Eq. (B16) can be shown by contracting the second formula in Eq. (B10) with δ ij and the third one with P i δjk, ∇ 2 N = 3N ′ + P 2 N ′′ , ( P • ∇)∇ 2 N = 5P 2 N ′′ + P 2 P 2 N ′′′ , (B17) and noticing that the left-hand side in each of them can be writen as a total divergence and therefore integrate to zero.Therefore, using the same notation, so that, and from Eq. (B8) we arrive at the final result quoted in Eq. (2.34), Appendix C: Evaluation of A ′ 0 f , C 0 f in the classical limit Let us consider A ′ 0 , defined in Eq. (2.35).Performing the angular integration, and using followed by an integration by parts, Let us consider the integral over f f (E), where in the classical limit that we are considering.Making the change of variable P = m cosh ξ , (C4) The second step follows by expanding cosh ξ up to the quadratic term in ξ and performing the Gaussian integration, and in the thrird step we have used the explicit relationship between the total number density n f and the chemical potential (e α ), which is different in the nonrelativistic (β f m ≫ 1) and the extremely-relativistic (β f m ≪ 1) limits.The integral over f f can be expressed in analogous fashion.
Thus, using the result given in Eq. (C5) (and the corresponding result for f f ) in Eq. (C2) we arrive at the formulas for A ′ 0 quoted in Eqs.(2.38) and (2.39).The formulas for C 0 quoted in there follow in similar fashion.
.21) where 1/β e,e ′ and µ e,e ′ are the temperature and chemical potential of the normal and stream background electrons, respectively.The corresponding formulas for the functions T (W ) f T,L,P corresponding to diagram (a) in Fig. 1 are obtained from Eq. (2.18) by making the replace- (c) Terms similar to those in (a) and (b), with F µν → Fµν Thus, under these conditions the most general form of Σ (B) is have used Eq.(3.8), and for V (i) f µ we have introduced the vectors are of order 1/m 2 W .To quantify somewhat these statements, recall that the O(1/m 4 W ) contribution to b e is[16] we explained in some detail the need to include the screening effects of the background electrons in the calculation of the self-energy in the two-stream medium case.The results for the B-dependent part of the self-energy are summarized in Eqs.(3.38) and (3.40).
We can ask the same question for the term t e ′ T ( k • m).Using Eqs.(2.39) and (3.40), the condition that this term be of the same order as b(4)On the other hand, if the medium is electron-positron asymmetric, the effects of the inhomogeneous terms seem to be unimportant compared to the standard terms.As an specific example, let us compare the t e ′ T ( k • m) against the term c e B, which is also a source of an asymmetry in the dispersion relation.The condition that it be of the same order as , and defined ∆N e = n e − n ē.If we assume, for example, that ∆N e ∼ N e ′ and T e ∼ T e ′ , the condition is