Relativistic shocks in conductive media

Section of Astrophysics, Astronomy and Mechanics, Department of Physics, National and Kapodistrian University of Athens, University Campus, 15784 Zografos, Athens, Greece
e-mail: arloules@phys.uoa.gr; vlahakis@phys.uoa.gr

Received: 24 July 2023
Accepted: 23 October 2023

Abstract

Context. Relativistic shocks are present in all high-energy astrophysical processes involving relativistic plasma outflows interacting with their ambient medium. While they are well understood in the context of relativistic hydrodynamics and ideal magnetohydrodynamics (MHD), there is a limited understanding of the properties related to their propagation in media characterized by finite electrical conductivity.

Aims. This work presents a systematic method for the derivation and solution of the jump conditions for relativistic shocks propagating in MHD media with finite electrical conductivity. This method is applied to the numerical solution of the Riemann problem and the determination of the conditions inside the blastwave that is formed when ultrarelativistic magnetized ejecta interact with the circumburst medium during a gamma-ray burst.

Methods. We derived the covariant relations expressing the jump conditions in a frame-independent manner. The resulting algebraic equations expressing the Rankine-Hugoniot conditions in the propagation medium’s frame were then solved numerically. A variable adiabatic index equation of state was used in order to obtain a realistic description of the post-shock fluid’s thermodynamics. This method was then employed for the solution of the Riemann problem for the case of a forward and a reverse shock, both of which form during the interaction of a gamma-ray burst ejecta with the circumburst medium. This allowed us to determine the kinematics of the resulting blastwave and the dynamical conditions in its interior.

Results. Our solutions clearly depict the impact of the plasma’s electrical conductivity in the properties of the post-shock medium. Two characteristic regimes are identified with respect to the value of a dimensionless parameter that has a linear dependence on the conductivity. For small values of this parameter, the shock affects only the hydrodynamic properties of the propagation medium and leaves its electromagnetic field unaffected. No current layer forms in the shock front; thus, we refer to this as the current-free regime. For large values of this parameter, the ideal MHD regime has been retrieved. We also show that the assumption of a finite electrical conductivity can lead to higher efficiencies in the conversion of the ejecta energy into thermal energy of the blastwave through the reverse shock. The theory developed in this work can be applied to the construction of Riemann solvers for resistive relativistic MHD (RR4MHD).