Galvano-rotational effect induced by electroweak interactions in pulsars

Galvano-rotational effect induced by electroweak interactions in pulsars


We study electroweakly interacting particles in rotating matter. The existence of the electric current along the axis of the matter rotation is predicted in this system. This new galvano-rotational effect is caused by the parity violating interaction between massless charged particles in the rotating matter. We start with the exact solution of the Dirac equation for a fermion involved in the electroweak interaction in the rotating frame. This equation includes the noninertial effects. Then, using the obtained solution, we derive the induced electric current which turns out to flow along the rotation axis. We study the possibility of the appearance of the galvano-rotational effect in dense matter of compact astrophysical objects. The particular example of neutron and hypothetical quark stars is discussed. It is shown that, using this effect, one can expect the generation of toroidal magnetic fields comparable with poloidal ones in old millisecond pulsars. We also briefly discuss the generation of the magnetic helicity in these stars. Finally we analyze the possibility to apply the galvano-rotational effect for the description of the asymmetric neutrino emission from a neutron star to explain pulsars kicks.

a,b,c]Maxim Dvornikov

Prepared for submission to JCAP

Galvano-rotational effect induced by electroweak interactions in pulsars

  • Institute of Physics, University of São Paulo, CP 66318, CEP 05314-970 São Paulo, SP, Brazil

  • Pushkov Institute of Terrestrial Magnetism, Ionosphere and Radiowave Propagation (IZMIRAN), 142190 Troitsk, Moscow, Russia

  • Physics Faculty, National Research Tomsk State University, 36 Lenin Ave., 634050 Tomsk, Russia


Keywords: neutron stars, magnetic fields, gravity

ArXiv ePrint: 1503.00608



1 Introduction

The importance of noninertial effects for various areas in modern physics cannot be underestimated. Some of the examples of these effects are mentioned in ref. [1]. One of the most common manifestations of noninertial effects is the description of physical processes in a rotating frame. One can expect the appearance of additional interesting phenomena if, besides the matter rotation, there is a parity violating interaction in the system. In the present work we will show that, in this situation, an electric current flowing along the rotation axis can be induced. We shall call this phenomenon as the new galvano-rotational effect (GRE).

Previously the generation of an electric current due to nontrivial topological effects was studied mainly in connection to the chiral magnetic effect (CME) [2]. In that case a nonzero current can be induced in the system of massless particles embedded in an external magnetic field [3], provided there is an imbalance between left and right particles. Recently, in ref. [4], we showed that the electric current can be generated even at zero chiral imbalance if charged particles are involved in the parity violating interaction.

We will apply the new GRE in astrophysical media to generate a toroidal magnetic field (TMF) inside a compact star. Although stellar TMFs cannot be observed directly, they are an internal ingredient of various astrophysical objects. For example, the most plausible explanation of the solar cycle is the oscillation between poloidal and toroidal components of the solar magnetic field [5]. Moreover, purely poloidal or toroidal stellar magnetic fields were shown in refs. [6, 7] to be unstable. Thus TMF is inherent to a magnetized star. In this work we show how TMF can be generated in an old millisecond pulsar basing on GRE.

Since macroscopic fluxes of electroweakly interacting particles are produced by GRE, we can try to apply this effect to explain linear velocities of millisecond pulsars. It is known from astronomical observations [8] that pulsars possess great linear velocities. Nevertheless physical processes underlaying pulsar kicks are still unclear. It might be reasonable to use GRE to account for pulsars kicks due to, e.g., anisotropic neutrino emission since linear velocities of pulsars were reported in ref. [9] to be correlated with their angular velocities.

This paper is organized in the following way. First, in section 2, we briefly describe the Standard Model interaction between leptons and quarks in flat space-time. Then, in section 3, we derive the Dirac equation for a fermion which interacts electroweakly with a rotating background matter. For this purpose we write down this Dirac equation in the corotating frame, using the method of an effective curved space-time. The exact solution of the Dirac equation for an ultrarelativistic fermion, accounting for the noninertial effects, is obtained in section 3. In section 4, we establish GRE, which, in this situation, consists in the appearance of the electric current along the rotation axis. We calculate this current in section 4 using the exact solution of the Dirac equation obtained in section 3. In section 5, we apply GRE for the generation of TMF and the magnetic helicity in compact rapidly rotating stars. Finally, in section 6, we try to use GRE to produce anisotropic neutrino emission from pulsars to explain their great linear velocities.

In section 7, we summarize our results and compare them with the findings of other authors. Some details of the derivation of the electric current in rotating matter are provided in appendix A.

2 Electroweak interaction of fermions in flat space-time

In this section we shall briefly remind the description of the electroweak interaction between leptons and quarks in flat space-time.

Let us consider a medium consisting of electrons, electron neutrinos as well as and quarks. We shall assume that quarks are both in confined states, forming nucleons, and hypothetical free particles. The effective Lagrangian for the electroweak interaction in this system in the Fermi approximation has form [10],


where is the Fermi constant, is the charged current, and is the neutral current.

In the considered system of elementary particles, the charged current has the form,


where are the wave functions of and quarks, are the wave functions of and an electron, , are the Dirac matrices, , and is the element of the Cabbibo-Kobayashi-Maskawa matrix. The neutral current can be expressed as


where and


Here and is the Weinberg angle.

In this section we shall consider the case when background fermions are at rest and unpolarized. While averaging the currents in eqs. (2.2) and (2.3) over the Fermi-Dirac distributions , in this case we get that only , where is the invariant number density of these fermions. The quantities and are vanishing since they are proportional to the macroscopic velocity and the polarization.

After averaging over the ensemble of background particles, we can rewrite eq. (2.1) in the form,


where is the wave function of a test fermion which undergoes a scattering off background particles and the effective potentials are given in table 1 for any scattering channels.

When we consider the electron and scattering off nucleons, and quarks are confined inside neutrons or protons. In the case of electron-nucleons interaction, only the neutral current contributes to eq. (2.5). If we study the scattering quarks off quarks, and vice versa, as well as the scattering off background fermions, both the charged and the neutral currents give the contributions to eq. (2.5). To derive the expression for for , , and interactions on the basis of eq. (2.2), we use the Fierz transformation. We also note that we consider the interaction with electroneutral matter where .

Test         particle Background particles
electron nucleons
nucleons & electrons 0
quark quarks
quark quarks
Table 1: The values of the effective potentials in eq. (2.5) for various channels of the scattering of a test fermion off background particles. Here is the electron density, are the densities of neutrons and protons, and are the densities of and quarks.

3 Dirac equation for a fermion interacting with rotating matter

In this section we shall find the exact solution of the Dirac equation for a fermion interacting with a rotating matter by means of the electroweak forces. The obtained solution includes noninertial effects.

Let us discuss the interaction of a fermion with matter rotating with the constant angular velocity . In this case we cannot directly apply the results of section 2 taking , where is the fermions velocity, while averaging over the ensemble of background particles. In the situation, when matter moves with an acceleration, one should account for possible noninertial effects.

Nevertheless we can still choose a noninertial reference frame where matter is at rest. Assuming that matter is unpolarized, we get that only in this reference frame. Thus, formally we can use the effective potentials derived in section 2. For the first time this approach was put forward in ref. [11], where the neutrino interaction with a rotating matter was studied.

It is known that the description of a particle in an accelerated frame is analogous to the motion of this particle in the curved space-time or the interaction with an effective gravitational field. For example, when we study the motion in the rotating frame, the interval takes the form [12],


where is the metric tensor of the effective gravitational field. Here we use the cylindrical coordinates .

Using eq. (2.5), we get that the Dirac equation for a test fermion, with the mass , involved in the parity violating interaction and moving in a curved space-time, has the form (see also ref. [11]),


where are the coordinate dependent Dirac matrices, is the covariant derivative, is the spin connection, , is the covariant antisymmetric tensor in curved space-time, , and are the effective potentials. Note that, since we choose a corotating frame, then and , where are given in table 1. Analogous Dirac equation was also discussed in ref. [13].

One can check that, using the following vierbein vectors:


the metric tensor in eq. (3.1) can be diagonalized, , where is the metric in a locally Minkowskian frame.

Let us introduce the constant Dirac matrices in a locally Minkowskian frame by . As shown in ref. [11], does not depend on coordinates. The spin connection in the Dirac eq. (3.2) has the form [14],


where are the generators of the Lorentz transformations in a locally Minkowskian frame and the semicolon stays for the covariant derivative. The explicit calculation on the basis of eq. (3.4) shows that the nonzero components of the connection one-form are


Using eqs. (3.4) and (3.5) we get that .

Using the definition of , the Dirac eq. (3.2) takes the form,


where are the vector and axial parts of the effective potentials. Note that the operator in eq. (3) is analogous to that recently studied in ref. [15].

Since eq. (3) does not explicitly depend on , , and , we shall look for its solution in the form,


where is the wave function depending on the radial coordinate. The values of in eq. (3.7) were found in ref. [16] to be .

It is convenient to rewrite eq. (3) as


where , is the effective electric charge, ,


is the potential of the effective electromagnetic field, and .

Let us look for the solution of eq. (3.8) in the form, , where is the new spinor. The equation for reads


where .

The solution of eq. (3.10) can be found for ultrarelativistic particles. In the limit , we can represent in eq. (3.10), where is a scalar function and is a constant spinor satisfying and , with and , since both and now commute with the operator of eq. (3.10).

Let us first study left particles, , corresponding to . The case of right particles with can be studied analogously. Using the new variable we can write the equation for as




In eq. (3.11) we choose , where .

Assuming that at , we get that and , for , as well as and , for . Here , , and is the Laguerre function111The definition of the Laguerre function is given, e.g., in ref. [11].. The energy spectrum can be found if we take that in eq. (3). We can present is the form,


One can see that the neutrino energy depends on the sign of . Note that one should understand in eq. (3.13) in the perturbative sense.

The total radial wave function can be found in the explicit form if we choose the spinors as


In eq. (3.14) we assume that the Dirac matrices are in the chiral representation [17],


where are the Pauli matrices.

Using eq. (3.14), we get the radial wave functions corresponding to different spin projections,






In eq. (3.17) the upper signs stay for and the lower ones for .

The properly normalized two component spinors in eq. (3.16) have the form,


for , and


for . The signs in eqs. (3.19) and (3.20) are correlated with the sign in eq. (3.13).

It should be noted that the spinors and in eqs. (3.19) and (3.20) are linearly dependent as it should be for ultrarelativistic particles. Thus, one can use any independent pair of and . As usual, we shall attribute with the upper sign to a particle degree of freedom and with the lower sign to antiparticles. This choice of independent spinors is convenient for . If , one can better use for and for as independent degrees of freedom.

Using the normalization condition for the total wave function,


which includes the dependence on and , we get the coefficients in eqs. (3.19) and (3.20) as


which is valid for any sign of .

On the basis of eqs. (3.19) and (3.20) one can notice that, at , is correlated with the particle helicity. Using eq. (3.13) with as well as eqs. (3.19) and (3.20), one can find the possible values of at . For the convenience, they are listed in table 2. Note that at , .

Values of Allowed value of
Left particles
Left antiparticles
Right particles
Right antiparticles
Table 2: Allowed values of in the ground state for different values of the effective potentials for left and right particles and antiparticles.

It should be noted that in table 2 is a formal quantum number. The physical value of for antiparticles is opposite to that shown in table 2: . Otherwise the electric charge would not be conserved.

Right particles can be treated in the same way as left ones. That is why we just present only the final results. The expression for the energy has the same structure as eq. (3.13) with the replacement . In eq. (3.16) one has , where


for , and


for . The argument of the Laguerre functions is now. The new normalization constant in eq. (3.23) and (3.24) reads


The signs in eqs. (3.23)-(3.25) are correlated with the signs in the expression for the energy.

As in the case of left particles, for right fermions we have that at . The allowed values of at are shown in table 2.

4 Induced electric current along the rotation axis

In this section we show that there is a nonzero induced electric current flowing along the rotation axis in the system of electroweakly interacting particles. In our calculation we shall use the exact solution of the Dirac equation obtained in section 3.

In section 3 we already mentioned that there is a correlation between and the helicity at . Thus one expects that there can be macroscopic fluxes of particles in the rotating matter. Let us first examine this issue for left fermions. We shall calculate the mean hydrodynamic currents of particles and antiparticles with respect to the coordinates in the rotating frame. These currents have the form,


where is the Fermi-Dirac distribution for fermions, with the lower sign staying for antifermions, is the reciprocal temperature of the fermion gas, and is the chemical potential of left particles. The spinor in eq. (4.1) corresponds to the exact solution of the Dirac equation in eqs. (3.19) and (3.20). Note that, for the first time, this method for the calculation of the current was proposed in ref. [3].

We will be interested in the expression for linear in . That is why we use instead of the total particle energy, cf. eq. (3.13), in the distribution function. The contribution of the noninertial part of the energy , see eq. (3.13), to the currents is computed in appendix A. We should study only since it is this component of the current that is linear in .

Using the orthogonality of Laguerre functions,


as well as eqs. (3.3), (3.7), (3.15), (3.16), (3.19), (3.20), and (3.22) we get that


where the upper sign stays for particles and the lower one for antiparticles.

On the basis of eq. (4.3) and table 2 we find that only the lowest level contributes to the current. Finally we get that


Analogously to eq. (4.4) we can obtain the following expression:


which is valid for antifermions.

The contribution to the hydrodynamic current from right femions is analogous to eqs. (4.4) and (4.5). It is


for particles, and


for antiparticles. Here are the distributions of right particles and antiparticles which can be obtained from by replacing , where is the chemical potential of right fermions.

Now, using eqs. (4.4)-(4.7), we can obtain the expression for the third component of the electric current as , where is the electric charge of the fermion including the sign. For example, for an electron, for an quark, and for a quark. Here is the absolute value of the elementary electric charge. In the expression for , we use the convention that the direction of the electric current coincides with the motion of the positive electric charge. Finally we get for the electric current,


where we restore vector notations. We remind that in eq. (4.8) we keep only the terms linear in and . Some details of the derivation of eq. (4.8) from eqs. (4.4)-(4.7) are provided in appendix A.

We can attribute the existence of the induced electric current in rotating matter, where the parity violating interaction is present, to the new GRE; cf. section 1.

5 Generation of TMF and magnetic helicity in a pulsar

In this section we apply GRE for the calculation of TMF inside a pulsar. We also briefly consider the generation of the magnetic helicity in a compact rotating star.

If we consider a rapidly rotating compact astrophysical object, like a neutron star (NS) or even a hypothetical quark star (QS), then the mechanism described in section 4 will induce the electric current along the rotation axis of such a star. We shall suppose that this current forms a closed circuit connected somewhere at the stellar surface. Thus, using the Maxwell equation , we get that this current should induce a TMF, , where is the stellar radius.

It should be noted that a compact star typically has a poloidal magnetic field , which is measured in astronomical observations. For instance, the radiation of a pulsar can be explained by the emission of electromagnetic waves by the rotating magnetic dipole associated with , provided there is a nonzero angle between and . However, as shown in ref. [18], using general arguments for the magnetohydrodynamic equilibrium of an axisymmetric NS, a purely poloidal magnetic field configuration turns out to be unstable. Thus an internal nonzero TMF should exist in a compact star.

Let us first consider the generation of TMF in NS composed of degenerate electrons and nucleons, like neutrons and protons. In this case and quarks are confined inside nucleons. The typical electron density in NS is , which corresponds to the electron fraction . It gives the chemical potential of electrons . Therefore electrons are ultrarelativistic whereas neutrons and protons are nonrelativistic. Note that the nonzero electron mass was shown in ref. [19] to slightly contribute to the induced electric current. Thus we can assume that electrons are approximately massless in NS and the results of sections 3 and 4 are valid.

Since the chiral symmetry is unbroken, we can consider left and right chiral projections as independent degrees of freedom. For simplicity we shall take that left and right electrons are in equilibrium with . The situations, when the chiral imbalance is important in NS, are studied, e.g., in refs. [4, 19]. Since in NS, for . Eventually, taking that and as well as using eq. (4.8), we get that can be generated in a rotating NS.

The obtained value of is comparable with in weakly magnetized old millisecond pulsars [20]. It should be noted that the stability of magnetic fields in NS can be reached if  [18], where is the energy of TMF and is the total magnetic energy. Our estimate for satisfies this criterion.

Let us discuss the creation of TMF in a hypothetical QS. Although QSs have not been observed yet, their properties are actively studied theoretically [21]. Various models of QS predict that it consists of free and quarks with some admixture of quarks. After the analysis of various equations of state of QS matter, the strangeness fraction was found in ref. [22] not to exceed . Thus we can approximately omit the contribution of quarks in the calculation of the induced current.

We shall suppose that inside QS we have and , where is the nuclear density. At such high densities the chiral symmetry can be unbroken in QS [23], allowing one to consider the independent chiral projections of the wave functions of and quarks. Therefore we can again use the results of sections 3 and 4. As in case of NS, we can also assume that left and right quarks are in equilibrium, , just for simplicity. Here and , where .

Finally, using, in eq. (4.8), the adopted values of densities and chemical potentials of quarks, the values of in table 1 as well as for and , we get that can exist inside a rotating QS. The obtained value of TMF is slightly greater than for NS. Note that the derived strength of TMF is also in agreement with stability criterion obtained in ref. [18] for old weakly magnetized millisecond pulsars [20].

It should be noted that, in our estimate of TMF in QS, we account for only and contributions to the electric current. Using the analogy of the rotating electroweak matter with the presence of an effective magnetic field [24] (see also eqs. (3.8) and (3.9) in section 3) and the results of ref. [25], we get that and interactions do not contribute to the current in eq. (4.8).

We should mention that, in the generation of TMF in a compact star, we discuss a thermally relaxed stage in the evolution of this astrophysical object. This approximation is valid since we consider old millisecond pulsars with ages  [20]. It means that we discard any possible effects related to turbulence which should be treated on the basis of the Navier-Stokes equation. In our analysis we also do not consider a differential rotation either.

The creation of TMF in a compact star is closely related to the problem of the generation of the magnetic helicity defined as


where is the 3D vector potential. If there is configuration of magnetic fields in a star consisting of toroidal and poloidal fields, then in eq. (5.1) has the form [26], , where and are the fluxes of toroidal and poloidal fields and is the linkage number. The magnetic helicity is a conserved quantity in a perfectly conducting medium. It is this fact which provides the stability of a poloidal field in a compact star. Note that another mechanism for the generation of the magnetic helicity in a nonrotating NS, based on the electron-nucleon electroweak interaction, was recently proposed in refs. [4, 27].

6 Pulsar kicks due to the asymmetric neutrino emission

In this section we shall use GRE for the description of the asymmetry in the neutrino emission from NS. We shall also consider the applicability of our results to explain linear velocities of pulsars.

It is well established that some pulsars have great linear velocities up to  [8]. There are various models for the explanation of this phenomenon based on, e.g., the asymmetric electromagnetic radiation [28] and the asymmetric explosion leading to the anisotropic neutrino emission [29]. We also mention ref. [30], where the asymmetry in neutrino oscillations in matter and an external magnetic field was used to account for pulsar kicks. The idea that anisotropically emitted electrons, owing to CME, pass the momentum to NS, was discussed in ref. [19]. Nevertheless the origin of peculiar velocities of pulsars is still unclear.

In ref. [9] it was established that linear velocities of pulsars are correlated with their angular velocities. Thus we may try to apply GRE, which predicts particle fluxes along the rotation axis, to explain pulsar kicks. However, unlike ref. [19], we shall examine the possibility of asymmetric neutrino emission since it is not very clear how charged particles can escape NS.

Using eqs. (4.4) and (4.5), we get the total hydrodynamic current of neutrinos along the rotation axis as


where and is the neutrino temperature. To derive eq. (6.1) we take into account that for (see table 1) and . Note that both and are equal to zero for ultrarelativistic neutrinos. That is why we account for the contribution of only left neutrinos in eq. (6.1).

The more complicated structure of the hydrodynamic current in eq. (6.1) compared to that of the electric current in eq. (4.8) is owing to the fact that has the same direction as , whereas is directed oppositely to . The function is plotted in figure 1. It is interesting to mention that . Thus, there is a nonzero neutrino flux even at .

Figure 1: The function in eq. (6.1) versus .

The simulations carried out in ref. [31] show that there is a significant nonzero neutrino asymmetry , owing to direct Urca processes, which lasts up to after the onset of the supernova collapse. When , the typical neutrino energy and chemical potential of is  [32]. At other neutrino species start to be emitted resulting in the diminishing of . When only pairs can be emitted in modified Urca processes [33] leading to . For simplicity we shall assume that and  [33] at this stage of the NS evolution.

The total momentum carried away by neutrinos during the time interval is , where is area of the equatorial cross section of NS and and is the NS radius. Using eq. (6.1), one gets that the greatest is achieved at . As a result, NS will get a recoil velocity , where is the NS mass. We shall take , where is the solar mass, to get the upper bound for . Taking that , we get that . The obtained value is sure to be beyond the possibility of modern astronomical observations. Therefore, despite there is a nonzero anisotropy in the neutrino emission in NS owing to GRE, this effect is unlikely to result in any testable phenomena.

7 Discussion

In conclusion we note that in the present work we have studied the evolution of particles, involved in the parity violating electroweak interaction, in the rotating matter. In section 3, we have obtained the new exact solution of the Dirac equation for a test ultrarelativistic particle, which account for the noninertial effects. Then, in section 4, we have computed the induced electric current along the rotation axis on the basis of the exact solution of the Dirac equation. In section 5, we have applied our results for the generation of TMF and the magnetic helicity in compact rotating stars. Finally, in section 6, we have considered the production of the anisotropy in the neutrino emission from NS, owing to GRE, and examined the applicability of this effect for the explanation of peculiar velocities of pulsars.

Several new results have been obtained in this work. Firstly, we mention that the vierbein vectors in eq. (3.3) have never been used previously in the Dirac eq. (3.2), which accounts for the electroweak interaction with background matter in curved space-time. Another veirbein was recently used in ref. [11]. However, the choice of the vierbein in the present work is likely to be more appropriate for ultrarelativistic particles in a rotating frame. In particular, here we have obtained the correct form of the “centrifugal” energy, or the energy of the rotation—angular momentum coupling, ; cf. eq. (3.13). The obtained expression for coincides with the result of ref. [34] derived on the basis of the general analysis. The form of obtained in ref. [11], where another vierbein was used, is slightly different. Therefore, the vierbein adopted in ref. [11] is likely to be more appropriate for the description of nonrelativistic particles in a rotating frame; cf. ref. [35].

Secondly, we have predicted the new GRE. This effect consists in the appearance of the electric current in the rotating matter composed of massless particles involved in the parity violating electroweak interaction. This electric current flows along the rotation axis. The new GRE is analogous to CME, known in QED, which consists in the generation of the electric current of massless charged particles along the external magnetic field [2, 3]. It should be noted that, for the first time, the analogy between the motion in a rotating electroweak matter and in an external magnetic field was mentioned in ref. [24].

Note that the appearance of the electric tension in a rotating conductor owing to the noninertial effects was also discussed in ref. [36]. The electric tension, predicted in ref. [36], is induced mainly by the Coriolis force acting on charged particles in a rotating conductor. If is chosen along the -axis, this tension is found in ref. [36] to be along the azimuthal direction. In our case, the electric current is owing to both the matter rotation and the presence of the parity violating electroweak interaction. We predict that the induced electric current flows along the rotation axis.

We have used the solution of a Dirac equation in the rotating frame for the calculation of the induced electric current. It means that this current flows inside the rotating matter since the quantum states of charged particles are measured by a corotating observer. For example, if one used the wave functions obtained in ref. [37], although they look similar to those in eqs. (3.19) and (3.20), we would get the electric current with respect the a nonrotating observer, which cannot be applied for the generation of the internal TMF.

We have used the calculated electric current to generate TMF and the magnetic helicity inside a rotating compact star, like NS or QS. The strength of TMF generated turned out to be moderate, , for both NS and QS. However, such TMF is comparable with a poloidal field in weakly magnetized old millisecond pulsars [20]. Note that the obtained strength of TMF is in agreement with a criterion for the magnetic field stability derived in ref. [18]. It should be noted that our model for the generation of TMF does not require the existence of a significant chiral imbalance between left and right charged particles. Such an imbalance is essential if CME is used to create TMF; cf. ref. [19].

Finally, in section 6, we analyzed the possibility to apply GRE to explain linear velocities of pulsars by the asymmetric neutrino emission from a rotating NS. We have estimated the total neutrino flux along the rotation axis as well as the recoil velocity of NS. It turned out that a pulsar kick caused by GRE is outside the observationally tested region.

At the end of this section we should make a comment on the influence of nonzero masses of charged particles on the generation of an electric current in a rotating star. It was mentioned in ref. [19] that the value of the current, induced by CME, slightly diminishes if a small, compared to the energy, but nonzero electron mass is accounted for. A more detailed analysis of the influence on the current, induced in frames of CME, from nonzero masses of flavored fermions, in case of strong interactions, was made in ref. [38]. It was shown in ref. [38] that CME disappears only in the great masses limit. The effect of the nonzero electron mass on the generation of strong magnetic fields in magnetars was also studied in ref. [39].

We should mention that, in case of GRE, a nonzero current in eq. (4.8) exists even at zero chiral imbalance: . Therefore, even if an initial is washed out owing to spin-flip processes taking place at a nonzero mass, as predicted in refs. [4, 27, 39], the current will be nonzero due to . For example, in section 5, we assumed that to simplify the estimates. In the model for the generation of strong magnetic fields in magnetars, elaborated in refs. [4, 27], appears inside NS in the course of the magnetic fields evolution. In this case, our estimates for obtained in section 5 will slightly change since .

GRE at and is likely to exist in all particular cases we considered in the present work since we discussed ultrarelativistic particles having , where is the fermion mass and is the typical fermion energy. Nevertheless, a more detailed quantum field theory analysis of this fact, like in ref. [38], is required.


I am thankful to V.G. Bagrov and V.B. Semikoz for useful discussions, to A.I. Studenikin for communications, to FAPESP (Brazil) for the Grant No. 2011/50309-2, to the Tomsk State University Competitiveness Improvement Program and to RFBR (research project No. 15-02-00293) for partial support.

Appendix A Details of the electric current calculation

In this appendix we derive the electric current in eq. (4.8) and discuss the approximations made.

Let us first consider left fermions. Using eq. (4.1), one gets that the contribution of particles and antiparticles to the electric current along the rotation axis is