Radiation mediated shocks in gamma-ray bursts: Pair creation

Relativistic sub-photospheric shocks are a possible mechanism for producing prompt gamma-ray burst (GRB) emission. Such shocks are mediated by scattering of radiation. We introduce a time-dependent, special relativistic code which dynamically couples Monte Carlo radiative transfer to the flow hydrodynamics. The code also self-consistently follows electron–positron pair production in photon–photon collisions. We use the code to simulate shocks with properties relevant to GRBs. We focus on plane-parallel solutions, which are accurate deep below the photosphere. The shock generates a power-law photon spectrum through the first-order Fermi mechanism, extending upward from the typical upstream photon energy. Strong (high Mach number) shocks produce rising νF ν spectra. We observe that in non-relativistic shocks the spectrum extends to , where v is the speed difference between the upstream and downstream. In relativistic shocks the spectrum extends to energies where its slope softens due to Klein–Nishina effects. Shocks with Lorentz factors γ > 1.5 are prolific producers of electron–positron pairs, yielding hundreds of pairs per proton. The main effect of pairs is to reduce the shock width by a factor of . Most pairs annihilate far downstream of the shock, and the radiation spectrum relaxes to a Wien distribution, reaching equilibrium with the plasma at a temperature determined by the shock jump conditions and the photon number per proton. We discuss the implications of our results for observations of radiation generated by sub-photospheric shocks.

Note:

14 pages, 8 figures, submitted to ApJ

gamma-ray burst: general
plasmas
radiation mechanisms: non-thermal
radiative transfer
scattering

References(40)

Figures(9)

[1]

Non-relativistic RMSs energize photons via the first-order Fermi process, which produces power-law photon spectra inside the shock transition region extending up to max = Emax/mec2 ~ v2/c2, where v is the upstream speed relative to the downstream. At energy Emax the photon energy gain from the Fermi process is balanced by the energy loss due to electron scattering recoil

[2]

GRB jets contain large numbers of photons, and the main RMS photon source is upstream advection. The average downstream photon energy {\bar{\epsilon }}_{d} is therefore known from the shock jump conditions and conservation of photon number. The power-law spectrum extends from the average upstream photon energy up to max, and the power-law index has the value that gives the predicted {\bar{\epsilon }}_{d}

[3]

Photon energy gain by scattering on the thermal electrons (and positrons) inside the RMS can be neglected as long as nγ/np mp/me, which is easily satisfied for GRBs

[4]

Mildly relativistic shocks with (v/c)γ 1 have max 1 and produce pairs via photon-photon annihilation

[5]

The power-law spectrum inside the RMS starts to curve downward at ≈ 10-1 due to Klein-Nishina effects. Photon-photon annihilation affects the spectrum strongly at 1

[6]

Pair production and annihilation is balanced inside the shock. This implies λγγ ≈ λsc, where λsc is the photon scattering mean free path and λγγ is the mean free path to photon-photon absorption at energies 1

[7]

GRB RMSs can generate huge amounts of pairs due to their large photon-to-proton ratios

[7]

we found Z± 200 for a shock with (v/c)γ = 3 and nγ/np = 2 × 105

[8]

The spatial shock width is modified by a factor {Z}_{\pm }^{-1}

[9]

The pairs annihilate at an optical depth of τ ~ ln Z± downstream of the RMS

[10]

The photons reach kinetic equillibrium (i.e., a Wien spectrum) at an optical depth \tau \sim 1/3{\bar{\epsilon }}_{d} downstream of the RMS

Crossref: