No Arabic abstract
We study the thermalization, injection, and acceleration of ions with different mass/charge ratios, $A/Z$, in non-relativistic collisionless shocks via hybrid (kinetic ions-fluid electrons) simulations. In general, ions thermalize to a post-shock temperature proportional to $A$. When diffusive shock acceleration is efficient, ions develop a non-thermal tail whose extent scales with $Z$ and whose normalization is enhanced as $(A/Z)^2$, so that incompletely-ionized heavy ions are preferentially accelerated. We discuss how these findings can explain observed heavy-ion enhancements in Galactic cosmic rays.
It is widely believed that the maximum energy of synchrotron photons when electrons are accelerated in shocks via the Fermi process is about 50 MeV (in plasma comoving frame). We show that under certain conditions, which are expected to be realized in relativistic shocks of gamma-ray bursts, synchrotron photons of energy much larger than 50 MeV (comoving frame) can be produced. The requirement is that magnetic field should decay downstream of the shock front on a length scale that is small compared with the distance traveled by the highest energy electrons before they lose half their energy; photons of energy much larger than 50 MeV are produced close to the shock front whereas the highest Lorentz factor that electrons can attain is controlled by the much weaker field that occupies most of the volume of the shocked plasma.
We perform first-principles relativistic particle-in-cell simulations of aligned pulsar magnetosphere. We allow free escape of particles from the surface of a neutron star and continuously populate the magnetosphere with neutral pair plasma to imitate pair production. As pair plasma supply increases, we observe the transition from a charge-separated electrosphere solution with trapped plasma and no spin-down to a solution close to the ideal force-free magnetosphere with electromagnetically-dominated pulsar wind. We calculate the magnetospheric structure, current distribution and spin-down power of the neutron star. We also discuss particle acceleration in the equatorial current sheet.
We present first-principles relativistic particle-in-cell simulations of the oblique pulsar magnetosphere with pair formation. The magnetosphere starts to form with particles extracted from the surface of the neutron star. These particles are accelerated by surface electric fields and emit photons capable of producing electron-positron pairs. We inject secondary pairs at locations of primary energetic particles, whose energy exceeds the threshold for pair formation. We find solutions that are close to the ideal force-free magnetosphere, with the Y-point and current sheet. Solutions with obliquities $lt 40^{circ}$ do not show pair production in the open field line region, because the local current density along magnetic field is below the Goldreich-Julian value. The bulk outflow in these solutions is charge separated, and pair formation happens in the current sheet and return current layer only. Solutions with higher inclinations show pair production in the open field line region, with high multiplicity of the bulk flow and the size of pair-producing region increasing with inclination. We observe the spin-down of the star to be comparable to MHD model predictions. The magnetic dissipation in the current sheet ranges between 20% for the aligned rotator and 3% for the orthogonal rotator. Our results suggest that for low obliquity neutron stars with suppressed pair formation at the light cylinder, the presence of phenomena related to pair activity in the bulk of the polar region, e.g., radio emission, may crucially depend on the physics beyond our simplified model, such as the effects of curved space-time or multipolar surface fields.
We present 3D general relativistic magnetohydrodynamic (GRMHD) simulations of the accretion flow surrounding Sagittarius A* that are initialized using larger-scale MHD simulations of the $sim$ 30 Wolf--Rayet (WR) stellar winds in the Galactic center. The properties of the resulting accretion flow on horizon scales are set not by ad hoc initial conditions but by the observationally constrained properties of the WR winds with limited free parameters. For this initial study we assume a non-spinning black hole. Our simulations naturally produce a $sim 10^{-8} M_odot$ yr$^{-1}$ accretion rate, consistent with previous phenomenological estimates. We find that a magnetically arrested flow is formed by the continuous accretion of coherent magnetic field being fed from large radii. Near the event horizon, the magnetic field is so strong that it tilts the gas with respect to the initial angular momentum and concentrates the originally quasi-spherical flow to a narrow disk-like structure. We also present 230 GHz images calculated from our simulations where the inclination angle and physical accretion rate are not free parameters but are determined by the properties of the WR stellar winds. The image morphology is highly time variable. Linear polarization on horizon scales is coherent with weak internal Faraday rotation.
The Sun is an active star that can launch large eruptions of magnetised plasma into the heliosphere, called coronal mass ejections (CMEs). These ejections can drive shocks that accelerate particles to high energies, often resulting in radio emission at low frequencies (<200 MHz). To date, the relationship between the expansion of CMEs, shocks and particle acceleration is not well understood, partly due to the lack of radio imaging at low frequencies during the onset of shock-producing CMEs. Here, we report multi-instrument radio, white-light and ultraviolet imaging of the second largest flare in Solar Cycle 24 (2008-present) and its associated fast CME (3038+/-288 km/s). We identify the location of a multitude of radio shock signatures, called herringbones, and find evidence for shock accelerated electron beams at multiple locations along the expanding CME. These observations support theories of non-uniform, rippled shock fronts driven by an expanding CME in the solar corona.