# Yaroslavl State University Preprint YARU-HE-99/06 astro-ph/9909154 On the Possible Enhancement of the Magnetic Field by Neutrino Reemission Processes in the Mantle of a Supernova

## Abstract

URCA neutrino reemission processes under the conditions in the mantle of a supernova with a strong toroidal magnetic field are investigated. It is shown that parity violation in these processes can be manifested macroscopically as a torque that rapidly spins up the region of the mantle occupied by such a field. Neutrino spin-up of the mantle can strongly affect the mechanism of further generation of the toroidal field, specifically, it can enhance the field in a small neighborhood of the rigid-body-rotating core of the supernova remnant.

PACS numbers: 97.60.Bw, 95.30.Cq

The problem of the shedding of the mantle in an explosion of a type-II supernova is still far from a complete solution . It is known that in several seconds after the collapse of a presupernova an anomalously high neutrino flux with typical luminosities $`L10^{52}`$ ergs/s is emitted from the neutrinosphere, which is approximately of the same size as the remnant core . In principle, such a neutrino flux could initiate a process leading to the shedding of the mantle as a result of the absorption and scattering of neutrinos by nucleons and the $`e^+e^{}`$ plasma of the medium . However, detailed calculations in spherically symmetric collapse models have shown that such processes are too weak for mantle shedding . In a magnetorotational model , mantle shedding is initiated by the outward pressure of a strong toroidal magnetic field generated by the differential rotation of the mantle with the core’s primary poloidal magnetic field ”frozen” into it. Indeed, as calculations show , when a mantle rotates with a millisecond period in a poloidal field $`B_010^{12}10^{13}`$ G a toroidal field $`B10^{15}10^{16}`$ G is generated in a time of order of a second. We note that the model in Ref.  contains a fundamental limitation on the energy of the toroidal field (it cannot exceed the kinetic energy of the ”core + shell” system) and therefore on the maximum field itself.

In the present letter we investigate the possibility that this magnetic field ”frozen in” the mantle is enhanced as the result of elementary neutrino reemission processes occurring in the mantle. We assume the mantle in the vicinity of the neutrinosphere to be a hot ($`T`$ several MeV) and quite dense (though transparent to neutrinos, $`\rho 10^{11}10^{12}`$ g/cm<sup>3</sup>) medium consisting of free nucleons and $`e^+e^{}`$ plasma. Under these conditions the dominant neutrino reemission processes are the URCA processes:

$`p+e^{}n+\nu _e,`$ (1)   
$`n+e^+p+\stackrel{~}{\nu }_e,`$ (2)   
$`n+\nu _ep+e^{},`$ (3)   
$`p+\stackrel{~}{\nu }_en+e^+.`$ (4)   
We note that $`\beta `$ decay is statistically suppressed in such a medium. The basic idea of this letter is as follows. In an external magnetic field, neutrinos are emitted and absorbed asymmetrically with respect to the direction of the magnetic field as a result of the parity violation in the processes (1)–(4. Therefore a macroscopic torque spinning up the mantle can arise in a toroidal field. It is known  that for an equilibrium neutrino distribution function such a neutrino-recoil momentum must be zero. However, the supernova region under consideration is nonequilibrium for neutrinos, so that the torque that arises in it is different from zero. Moreover, as we shall show below, the torque can be large enough to change substantially the distribution of the angular rotational velocities of the mantle in the region filled with a strong magnetic field over the characteristic neutrino emission times. According to the equation governing the generation of a toroidal field , a large change in the gradient of the angular velocities in the region can lead to redistribution of the magnetic field (specifically, enhancement of the field in a small neighborhood of the rigid-body-rotating core).

A quantitative estimate of the effect follows from the expression for the energy-momentum transferred by neutrinos to a unit volume of the mantle per unit time:

$`{\displaystyle \frac{dP_\alpha }{dt}}={\displaystyle \frac{1}{V}}{\displaystyle \underset{i}{}dn_if_i\underset{f}{}dn_f\left(1f_f\right)\frac{\left|S_{if}\right|^2}{𝒯}k_\alpha },`$ (5)   
where $`dn_i`$ and $`dn_f`$ are the number of initial and final states in an element of the phase space, $`f_i`$ and $`f_f`$ are the distribution functions of the initial and final particles, $`k_\alpha `$ is the neutrino momentum, $`|S_{if}|^2/𝒯`$ is the squared $`S`$-matrix element per unit time. It is of interest to calculate the latter under our conditions, since, as far as we know, the URCA processes (1)–(4) have been previously studied for relatively weak  ($`Bm_e^2/e`$) and very strong  ($`Bm_p^2/e`$) fields. Assuming that electrons and positrons in the plasma mainly occupy only the lowest Landau level ($`\mu _e\sqrt{2eB}`$, where $`\mu _e`$ is the chemical potential of electrons), we obtained the following expression for the squared $`S`$-matrix element summed over all proton Landau levels and polarizations of the final particles and averaged over the polarizations of the initial particles:

$`\left|S_{if}\right|^2`$ $`=`$ $`{\displaystyle \frac{G_F^2\mathrm{cos}^2\theta _c\left(2\pi \right)^3𝒯}{2L_yL_zV^2}}{\displaystyle \frac{\mathrm{exp}\left(Q_{}^2/2eB\right)}{4\omega \epsilon }}`$   
$`\times `$ $`{\displaystyle \frac{1}{2}}\left[{\displaystyle \underset{n=0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{\left|M_+\right|^2}{n!}}\left({\displaystyle \frac{Q_{}^2}{2eB}}\right)^n\delta ^{\left(3\right)}+{\displaystyle \underset{n=1}{\overset{\mathrm{}}{}}}{\displaystyle \frac{\left|M_{}\right|^2}{\left(n1\right)!}}\left({\displaystyle \frac{Q_{}^2}{2eB}}\right)^{n1}\delta ^{\left(3\right)}\right],`$   
$`\left|M_\sigma \right|^2`$ $`=`$ $`4\left(1+g_a\sigma \right)^2\left[2\left(up\right)\left(uk\right)\left(pk\right)\left(up\right)\left(k\stackrel{~}{\phi }u\right)\left(uk\right)\left(p\stackrel{~}{\phi }u\right)\right]`$   
$`+`$ $`8g_a^2\left(1+\sigma \right)\left[\left(pk\right)\left(p\stackrel{~}{\phi }k\right)\right],`$   
where $`u_\alpha `$ is a four-velosity of a medium, $`𝐁=(0,0,B)`$, and $`\delta ^{(3)}`$ is the production of the energy delta-function and two delta-functions of the momentum in the direction of the magnetic field and one transverse momentum component which are conserved in the reactions; $`(pk)=\epsilon \omega p_3k_3`$, $`p^\alpha =(\epsilon ,𝐩)`$ and $`k^\alpha =(\omega ,𝐤)`$ are the four-momenta of the electron and neutrino, $`Q_{}^2`$ is the square of the momentum transfer transverse to the to the magnetic field in the reactions (1) and (3) \[with the corresponding substitutions $`pp`$ and $`kk`$ in the crossing reactions (2) and (4)\], $`\stackrel{~}{\phi }_{\alpha \beta }=\stackrel{~}{F}_{\alpha \beta }/B`$ is the dimensionless dual magnetic-field tensor, $`\sigma =\pm 1`$ is the projection of the proton spin on the direction of the magnetic field, $`n`$ is the summation index over the proton Landay levels, $`𝒯V=𝒯L_xL_yL_z`$ is the normalization four-volume, $`g_a`$ is the axial constant of the nucleonic current, $`G_F`$ is the Fermi constant, and $`\theta _c`$ is the Cabibbo angle. We note that in the limit of a strong magnetic field, when the protons occupy only the ground Landay level, expressions (  Yaroslavl State University  Preprint YARU-HE-99/06  astro-ph/9909154  On the Possible Enhancement of the Magnetic Field  by Neutrino Reemission Processes  in the Mantle of a Supernova) and (  Yaroslavl State University  Preprint YARU-HE-99/06  astro-ph/9909154  On the Possible Enhancement of the Magnetic Field  by Neutrino Reemission Processes  in the Mantle of a Supernova) agree with the result obtained previosly in Ref. .

Analysis shows that URCA processes are the fastest reactions in medium considered, and they transfer the medium into a state of $`\beta `$ equilibrium in a time of order $`10^2`$ s. Therefore we employed the condition of $`\beta `$ equilibrium and singled out in the expression for the energy-momentum transfer to the shell (5) the separate contributions from processes involving neutrinos (1), (3) and antineutrinos (2), (4):

$`{\displaystyle \frac{dP_\alpha ^{(\nu ,\stackrel{~}{\nu })}}{dt}}={\displaystyle \frac{d^3k}{\left(2\pi \right)^3}k_\alpha \left[1+\mathrm{exp}\left(\frac{\omega +\mu _{(\nu ,\stackrel{~}{\nu })}}{T}\right)\right]𝒦^{(\nu ,\stackrel{~}{\nu })}\delta f^{(\nu ,\stackrel{~}{\nu })}}.`$ (8)   
Here, $`\delta f^{(\nu ,\stackrel{~}{\nu })}`$ is the deviation of the distribution function from the equilibrium function, $`𝒦^{(\nu ,\stackrel{~}{\nu })}`$ is the (anti)neutrino absorption coefficient, defined as

$`𝒦^{(\nu ,\stackrel{~}{\nu })}={\displaystyle \underset{i}{}dn_if_i\underset{f}{}dn_f\left(1f_f\right)\frac{\left|S_{if}\right|^2}{𝒯}},`$ (9)   
where the integration extends over all states except the neutrino states in the reaction (3) and antineutrino states in the reaction (4), respectively. As follows from Eq. (8), actually, the momentum transferred to the medium is different from zero only if the neutrino distribution function deviates from the equilibrium distribution.

To calculate the absorption coefficient $`𝒦^{(\nu ,\stackrel{~}{\nu })}`$ we assumed that the ultrarelativistic $`e^+e^{}`$ plasma occupies only the ground Landau level, while the protons occupy quite many levels (the dimensionless parameter $`\delta =eB/m_pT1`$). We also used the fact that at the densities under consideration the nucleonic gas is Boltzmannian and nonrelativistic. Then, dropping terms $`\delta `$, we can write expression (9) in the form

$`𝒦^{(\nu ,\stackrel{~}{\nu })}`$ $`=`$ $`{\displaystyle \frac{G_F^2\mathrm{cos}^2\theta _ceBN_{(n,p)}}{2\pi }}\left[\left(1+3g_a^2\right)\left(g_a^21\right)k_3/\omega \right]`$ (10)   
$`\times `$ $`\left[1+\mathrm{exp}\left({\displaystyle \frac{\pm \left(\mu _e\mathrm{}\right)\omega }{T}}\right)\right]^1,`$   
where $`N_n`$, $`N_p`$, and $`m_n`$, $`m_p`$ are the number densities and masses of the neutrons and protons, respectively, $`\mathrm{}=m_nm_p`$, and $`\omega `$ and $`k_3`$ are the neutrino energy and the neutrino momentum in the direction of the magnetic field, respectively.

For further calculations we employed the neutrino distribution function in the model of a spherically symmetric collapse of a supernova in the absence of a magnetic field . This is a quite good approximation when the region occupied by the strong magnetic field is smaller than or of the order of the neutrino mean-free path. By the strong field we mean the field in which $`e^+e^{}`$ plasma occupies only the ground Landau level: $`eB\mu _e^2`$. In the model of Ref.  the region occupied by such a field is no greater than several kilometers in size, and we estimate the neutrino mean-free path in this region as

$`l_\nu 4\mathrm{km}\left({\displaystyle \frac{4.4\times 10^{16}\mathrm{G}}{B}}\right)\left({\displaystyle \frac{5\times 10^{11}\mathrm{g}/\mathrm{cm}^3}{\rho }}\right).`$ (11)   
Therefore the magnetic field cannot strongly alter the neutrino distribution function, and our approximation is quite correct.

As calculations of the components of the energy-momentum (8) transferred to the medium during neutrino reemission showed, the radial force arising is much weaker than the gravitational force and cannot greatly influence the mantle dynamics. However, the force acting in the direction of the magnetic field can change quite rapidly the distribution of the angular velocities in the region occupied by the strong magnetic field. The density of this force can be represented as

$`\mathrm{}_{}^{\left(tot\right)}=\mathrm{}_{}^{\left(\nu \right)}+\mathrm{}_{}^{\left(\stackrel{~}{\nu }\right)}=𝒩\left[\left(3\mu ^2_\nu 1\right)I\left(a\right)e^a+\left(3\mu ^2_{\stackrel{~}{\nu }}1\right)I\left(a\right)\right],`$ (12)   
$`\mu ^2=\left({\displaystyle \mu ^2\omega fd^3k}\right)\left({\displaystyle \omega fd^3k}\right)^1,`$   
where $`\mu `$ is the cosine of the angle between the neutrino momentum and the radial direction, $`a=\mu _e/T`$, and

$`I\left(a\right)={\displaystyle \underset{0}{\overset{\mathrm{}}{}}}{\displaystyle \frac{y^3dy}{\mathrm{e}^{ya}+1}}.`$   
In deriving Eq. (12) we use the one-dimensional factorized neutrino distribution function $`f^{(\nu ,\stackrel{~}{\nu })}=\varphi ^{(\nu ,\stackrel{~}{\nu })}(\omega /T_\nu )\mathrm{\Phi }^{(\nu ,\stackrel{~}{\nu })}(r,\mu )`$ , where $`T_\nu `$ is the neutrino spectral temperature and $`r`$ is the distance from the core center. To estimate the force in the diffusion region we assumed that $`T_\nu T`$ and chose $`\varphi ^{(\nu ,\stackrel{~}{\nu })}(\omega /T_\nu )=\mathrm{exp}(\omega /T_\nu )`$. We determined the dimensional parameter $`𝒩`$ in expression (12) as

$`𝒩`$ $`=`$ $`{\displaystyle \frac{G_F^2\mathrm{cos}^2\theta _c}{\left(2\pi \right)^3}}{\displaystyle \frac{g_a^21}{3}}eBT^4N_N`$   
$``$ $`4.5\times 10^{20}{\displaystyle \frac{\mathrm{dynes}}{\mathrm{cm}^3}}\left({\displaystyle \frac{T}{5\mathrm{MeV}}}\right)^4\left({\displaystyle \frac{B}{4.4\times 10^{16}\mathrm{G}}}\right)\left({\displaystyle \frac{\rho }{5\times 10^{11}\mathrm{g}/\mathrm{cm}^3}}\right),`$   
where $`N_N=N_n+N_p`$ is the total nucleon number density.

The force (12) was estimated numerically in the diffusion region of the supernova atmosphere for typical (excluding the field) values of the macroscopic parameters for this region: $`T=5`$ MeV, $`B=4.4\times 10^{16}`$ G, $`\rho =5\times 10^{11}`$ g/cm<sup>3</sup>. For these values $`a3`$, $`\mu ^2_\nu \mu ^2_{\stackrel{~}{\nu }}0.4`$ (Ref. ), and the force density in the direction of the field can be estimated from Eq. (12) as

$`\mathrm{}_{}^{tot}\mathrm{}_{}^\nu 𝒩.`$ (14)   
We note that the angular acceleration produced by the torque exerted by such a force is large enough to spin up the region of the mantle containing a strong magnetic field to typical angular velocities of a fast pulsar (with the rotational period $`P_010^2`$ s) in a characteristic time of order of a second. In our opinion, this result is of interest in itself and can serve as a basic for a number of applications. However, we shall give a qualitative discussion of only one possible manifestation of this result – the effect of such a fast spin-up of the mantle on the further generation of the toroidal magnetic field. Indeed, if the modification of the gradient of the angular velocities of the mantle is large, the toroidal magnetic field in the mantle at subsequent times will vary according to a law that is different from the linear law . Analysis of the equation governing the generation of a toroidal field with allowance for the force (14), which is linear in this field, leads to the conclusion that its growth in time is much faster (exponential) in the quite small region in which the force acts ($`eB\mu _e^2`$). However, the main source of the magnetic field energy, just as in the case when there is no force, is the kinetic energy of the rigid-body-rotating core. Thus the force (14) can lead to a peculiar rearrangement of the region occupied by the strong field. Specifically, with virtually no change in energy, the magnetic field can become concentrated in a narrower spatial region and can therefore have in this region higher intensities, on average, than in the absence of the spin-up force. We note that the effect under discussion can strongly influence the mantle-shedding process and also the mechanism leading to the formation of an anisotropic $`\gamma `$-ray burst in the explosion of a supernova with a rapidly rotating core . However, in order to perform detailed calculations of the generation of a toroidal field, the enhancement of the field due to the ”neutrino spin-up”, and the effect of this enhancement on the indicated processes, it is necessary to analyze the complete system of MHD equations. That analysis lies far outside the scope of the present work and is a subject of a separate investigation. Qualitative estimates show that the magnetic field with the strength $`B10^{17}`$ G can be generated by the above-described mechanism in a small neighborhood (of order a kilometer) of a rapidly rotating core with the period $`P_05\times 10^3`$ s, and this region decreases with increasing the period.

In summary, we have shown that URCA neutrino reemission processes can produce, in the region of the mantle that is filled with a strong toroidal magnetic field, angular accelerations which are sufficiently large as to greatly influence the mechanism of further generation of the field. Specifically, such rapid redistribution of the angular velocities can enhance the field in a small neighborhood of the rigid-body-rotating core of a remnant.

We are grateful to S. I. Blinnikov for fruitful discussions and for assistance in formulating the problem and to N. V. Mikheev and M. V. Chistyakov for helpful discussions. This work was partially supported by INTAS (Grant No. 96-0659) and the Russian Foundation for Basic Research (Grant No. 98-02-16694).