Relativistic Expansion of Magnetic Loops at the Self-similar Stage II: Magnetized outflows interacting with the ambient plasma

TL;DR

This paper derives and validates self-similar solutions for relativistically expanding magnetic loops interacting with ambient plasma, including shock discontinuities, and explores their stability and behavior under different ambient density profiles through simulations.

Contribution

It extends previous models by incorporating shock discontinuities and performs simulations to verify stability and analyze effects of ambient density variations.

Findings

Solutions are stable over simulation time.

Shock Lorentz factor increases with radius, proportional to r^{0.25}.

Shock Lorentz factor varies with time as t^{( extdelta-3)/2} depending on ambient density profile.

Abstract

We obtained self-similar solutions of relativistically expanding magnetic loops by assuming axisymmetry and a purely radial flow. The stellar rotation and the magnetic fields in the ambient plasma are neglected. We include the Newtonian gravity of the central star. These solutions are extended from those in our previous work (Takahashi, Asano, & Matsumoto 2009) by taking into account discontinuities such as the contact discontinuity and the shock. The global plasma flow consists of three regions, the outflowing region, the post shocked region, and the ambient plasma. They are divided by two discontinuities. The solutions are characterized by the radial velocity, which plays a role of the self-similar parameter in our solutions. The shock Lorentz factor gradually increases with radius. It can be approximately represented by the power of radius with the power law index of 0.25. We also carried out magnetohydrodynamic simulations of the evolution of magnetic loops to study the stability and the generality of our analytical solutions. We used the analytical solutions as the initial condition and the inner boundary conditions. We confirmed that our solutions are stable over the simulation time and that numerical results nicely recover the analytical solutions. We then carried out numerical simulations to study the generality of our solutions by changing the power law index \delta of the ambient plasma density \rho_0 \propto r^{-\delta}. We alter the power law index \delta from 3.5 in the analytical solutions. The analytical solutions are used as the initial conditions inside the shock in all simulations. We observed that the shock Lorentz factor increases with time when \delta is larger than 3, while it decreases with time when \delta is smaller than 3. The shock Lorentz factor is proportional to t^{(\delta-3)/2}. These results are consistent with the analytical studies by Shapiro (1979).