27.12.2012 Views

Proc. Neutrino Astrophysics - MPP Theory Group

Proc. Neutrino Astrophysics - MPP Theory Group

Proc. Neutrino Astrophysics - MPP Theory Group

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

So far the damping rate corresponds to the zero field case, but the generalization to a<br />

background magnetic field is straightforward after identifying the different factors in Eq. (3).<br />

The δ-functions come from the imaginary part of the usual electron and neutrino propagators<br />

(or rather from the thermal part of the real-time propagators). We now have to use the<br />

electron propagator in an external field instead, and here it is convenient to use Schwinger’s<br />

formulation<br />

i(�P + me)<br />

P 2 − m 2 e + iǫ →<br />

� ∞<br />

ds<br />

0<br />

eieBsσz<br />

cos(eBs) exp<br />

�<br />

is(P 2 tan(eBs)<br />

� +<br />

�<br />

�P� + e−ieBsσz<br />

cos(eBs) �P⊥<br />

�<br />

+ me<br />

eBs P 2 ⊥ − m2 e<br />

�<br />

+ iǫ)<br />

81<br />

, (5)<br />

where a·b � = a0bb −azbz and a·b⊥ = −axbx −ayby. This expression is quite difficult to use in<br />

general but the complicated parts, cos(eBs) and tan(eBs)/eBs, have no linear dependence<br />

on eB and can therefore be approximated by 1 for weak fields. The s-integral can then be<br />

performed and after taking the imaginary part we end up with the replacement rule<br />

(�P + me)δ(P 2 − m 2 e ) → (�P � + me)eBσzsign(P 2 − m 2 e )δ((P2 − m 2 e )2 − (eB) 2 ) . (6)<br />

To obtain this rule I have kept the full eB dependence in the exponentials in Eq. (5) since<br />

they are necessary for the IR convergence of the energy integrals. In the final result, after<br />

performing all integrals, the damping rate has a very clear linear dependence on eB for<br />

weak fields.<br />

To compute the damping from absorption and creation of neutrinos in URCA processes<br />

one can use the same formalism but with neutron and proton propagators in the loop. In the<br />

present case absorption is dominted by ν + n → p + e and creation by the inverse process.<br />

The nucleons are also essentially nonrelativistic which simplifies the matrix elements. On the<br />

other hand, since the general expression comes out for free and the trace is easily evaluated<br />

with a symbolic algebraic program, we might as well keep the full relativistic expression and<br />

all possible channels when we perform the numerical integrals.<br />

Results<br />

The physical conditions in a supernova vary considerably during the short time of the explosion.<br />

I have as an illustration evaluated the damping asymmetry for two more or less typical<br />

conditions. The first case has a high electron chemical potential (µe = 200 MeV) and neutrino<br />

energy q0 = 100 MeV, and the other one has µe = 50 MeV and q0 = 30 MeV. In both cases<br />

I used T = 10 MeV and Ye = 0.3. Local thermal equilibrium was assumed and the proton<br />

chemical potential was determined from charge neutrality. These two cases correspond to<br />

the densities 2×10 14 g/cm 3 and 4×10 12 g/cm 3 . For simplicity the incoming neutrino has its<br />

momentum parallel to the magnetic field (the linear term in eB is of course antisymmetric<br />

in qz). Writing the damping factor as Γ = Γ (0) + eB Γ (1) , I found the relative asymmetry for<br />

URCA and ν-e processes to be (B14 is the field strength in units of 10 14 Gauss)<br />

�<br />

�<br />

�eB<br />

Γ<br />

�<br />

�<br />

(1)<br />

URCA<br />

Γ (0)<br />

�<br />

�<br />

�<br />

�<br />

�<br />

URCA<br />

≃ 3.6×10−6 B14 ,<br />

�<br />

�<br />

�<br />

�<br />

eB Γ<br />

�<br />

�<br />

(1)<br />

ν−e<br />

Γ (0)<br />

�<br />

�<br />

�<br />

�<br />

�<br />

ν−e � ≃ 9.8×10−7 B14 , (7)

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!