Cosmic-ray Magnetohydrodynamics


Background and Motivation

Cosmic rays (CRs) are extremely energetic particles that propagate throughout the universe, consisting of protons, electrons, heavy nuclei, etc. They can be accelerated to energies beyond \(10^{20}~{\rm eV}\) with possible sources including relativistic jets and powerful astrophysical shocks. Much effort has been made starting from the early twentieth century to present day in understanding their origin, properties, and role in various environments in the universe, from the interstellar medium (ISM), to galaxies, and to galaxy clusters. The study of CRs has aided in advances in a variety of fields, including particle physics, cosmology, high-energy astrophysics, and plasma astrophysics. Open questions still remain to this day, such as the origin of the ultra-high-energy CRs and the nature of particle acceleration. Continual advances, however, in numerical methods, high-performance computing architectures, observational technology, and laboratory-astrophysics techniques show promise for further unraveling the mysteries behind CRs.

The two main deviations from a pure power-law are the knee, where the energy spectrum steepens slightly near \(10^{15}~{\rm eV}\), and at the ankle, where the spectrum slightly flattens near \(10^{18}~{\rm eV}\). These deviations in the power-law may indicate a transition of different populations of CRs, where it is theorized that CRs below the knee originate in the Galaxy and CRs above the ankle are extra-galatic in origin (Aharonian 2004, Beatty and Westerhoff 2009).

One of the widely accepted sources of CRs is supernova remnants (SNRs) through diffusive shock acceleration (DSA), where particles are accelerated from thermal energies to non-thermal energies at the shock front. DSA has been mostly studied in the context of CR transport in collisionless shocks, where injected particles into the shock excite flow instabilities as DSA progresses. A complete understanding of how the mechanism manifests in supernovae has yet to be achieved. Surveys by the Chandra X-Ray Observatory and Fermi Gamma-Ray Telescope have provided a wealth of data for validating theoretical models and studying supernova processes at thermal and non-thermal energies.

Supernovae are stellar explosions that generate strong shocks. This work does not focus on the cause of such explosions, but rather on the conditions after the explosion that may be suitable for CR production. Supernovae explosions eventually become supernova remnants (SNRs). Not only are supernovae of interest to the astrophysics community, they are also of interest to the high energy-density physics (HEDP) community. Technological advancements have made it possible to control and study astrophysical processes in the laboratory.

Supernova-remnant theory

Supernovae can occur either by the thermonuclear disruption of a white dwarf or by the core-collapse of a massive star. The first case is referred to as a type Ia supernova, where it is assumed that the explosion occurs when a massive white dwarf accretes matter from a companion and the mass exceeds the Chandrasekhar mass limit. Here, the explosion mechanism is reasonable well understood. The second case is referred to as a type II supernova, where the process is reasonably well understood with the exception of the late stages of the explosion mechanism. Several models have been proposed for this last stage. For example, the proper motion of a neutron star can be investigated because it depends on the kick a neutron star gains during the explosion when they are formed. Therefore, the proper motion strongly depends on the final explosion mechanism. The focus of this work will be the remnants of such explosions. SNRs are generally grouped by their morphology into shell-type SNRs, plerions, composite, and thermal-composite SNRs. A shell-like structure is created when the shock wave propagates through the ISM and heats the plasma. Hence, an SNR that shows such a structure is named a shell-type SNR and this group includes most known remnants. A plerion is a highly energetic wind of relativistic electrons and positrons, accelerated to ultra-relativistic energies by a rapidly rotating neutron star (pulsar), also referred to as a pulsar-wind nebulae (Kargaltsev et al. 2013). Remnants with a pulsar-wind nebula surrounded