We use large hybrid simulations to study ion acceleration and generation of magnetic turbulence due to the streaming of particles that are self-consistently accelerated at non-relativistic shocks. ...When acceleration is efficient, we find that the upstream magnetic field is significantly amplified. The total amplification factor is larger than 10 for shocks with Alfvenic Mach number M = 100, and scales with the square root of M. The spectral energy density of excited magnetic turbulence is determined by the energy distribution of accelerated particles, and for moderately strong shocks (M lap 30) agrees well with the prediction of resonant streaming instability, in the framework of quasilinear theory of diffusive shock acceleration. For M lap 30, instead, Bell's non-resonant hybrid (NRH) instability is predicted and found to grow faster than resonant instability. NRH modes are excited far upstream by escaping particles, and initially grow without disrupting the current, their typical wavelengths being much shorter than the current ions' gyroradii. Then, in the nonlinear stage, most unstable modes migrate to larger and larger wavelengths, eventually becoming resonant in wavelength with the driving ions, which start diffuse. Ahead of strong shocks we distinguish two regions, separated by the free-escape boundary: the far upstream, where field amplification is provided by the current of escaping ions via NRH instability, and the shock precursor, where energetic particles are effectively magnetized, and field amplification is provided by the current in diffusing ions. The presented scalings of magnetic field amplification enable the inclusion of self-consistent microphysics into phenomenological models of ion acceleration at non-relativistic shocks.
Abstract
We use fully kinetic particle-in-cell simulations with unprecedentedly large transverse box sizes to study particle acceleration in weakly magnetized mildly relativistic shocks travelling at ...a velocity ≈ 0.75c and a Mach number of 15. We examine both subluminal (quasi-parallel) and superluminal (quasi-perpendicular) magnetic field orientations. We find that quasi-parallel shocks are mediated by a filamentary non-resonant (Bell) instability driven by returning ions, producing magnetic fluctuations on scales comparable to the ion gyroradius. In quasi-parallel shocks, both electrons and ions are accelerated into non-thermal power laws whose maximum energy grows linearly with time. The upstream heating of electrons is small, and the two species enter the shock front in rough thermal equilibrium. The shock’s structure is complex; the current of returning non-thermal ions evacuates cavities in the upstream that form filaments of amplified magnetic fields once advected downstream. At late times, 10 per cent of the shock’s energy goes into non-thermal protons and ${\gtrsim }10{{\ \rm per\ cent}}$ into magnetic fields. We find that properly capturing the magnetic turbulence driven by the non-thermal ions is important for properly measuring the energy fraction of non-thermal electrons, εe. We find εe ∼ 5 × 10−4 for quasi-parallel shocks with v = 0.75c, slightly larger than what was measured in simulations of non-relativistic shocks. In quasi-perpendicular shocks, no non-thermal power-law develops in ions or electrons. The ion acceleration efficiency in quasi-parallel shocks suggests that astrophysical objects that could host mildly relativistic quasi-parallel shocks – for example, the jets of active galactic nuclei or microquasars – may be important sources of cosmic rays and their secondaries, such as gamma-rays and neutrinos.
We used unprecedentedly large two-dimensional and three-dimensional hybrid (kinetic ions-fluid electrons) simulations of non-relativistic collisionless strong shocks in order to investigate the ...effects of self-consistently accelerated ions on the overall shock dynamics. The current driven by suprathermal particles streaming ahead of the shock excites modes transverse to the background magnetic field. The Lorentz force induced by these self-amplified fields tends to excavate tubular, underdense, magnetic-field-depleted cavities that are advected with the fluid and perturb the shock surface, triggering downstream turbulent motions. These motions further amplify the magnetic field, up to factors of 50-100 in knot-like structures. Once downstream, the cavities tend to be filled by hot plasma plumes that compress and stretch the magnetic fields in elongated filaments; this effect is particularly evident if the shock propagates parallel to the background field. Highly magnetized knots and filaments may provide explanations for the rapid X-ray variability observed in RX J 1713.7-3946 and for the regular pattern of X-ray bright stripes detected in Tycho's supernova remnant.
Abstract Magnetized plasma columns and extended magnetic structures with both footpoints anchored to a surface layer are an important building block of astrophysical dissipation models. Current loops ...shining in X-rays during the growth of plasma instabilities are observed in the corona of the Sun and are expected to exist in highly magnetized neutron star magnetospheres and accretion disk coronae. For varying twist and system sizes, we investigate the stability of line-tied force-free flux tubes and the dissipation of twist energy during instabilities using linear analysis and time-dependent force-free electrodynamics simulations. Kink modes ( m = 1) and efficient magnetic energy dissipation develop for plasma safety factors q ≲ 1, where q is the inverse of the number of magnetic field line windings per column length. Higher-order fluting modes ( m > 1) can distort equilibrium flux tubes for q > 1 but induce significantly less dissipation. In our analysis, the characteristic pitch μ ˜ 0 of flux-tube field lines determines the growth rate ( ∝ μ ˜ 0 3 ) and minimum wavelength of the kink instability ( ∝ μ ˜ 0 − 1 ). We use these scalings to determine a minimum flux tube length for the growth of the kink instability for any given μ ˜ 0 . By drawing analogies to idealized magnetar magnetospheres with varying regimes of boundary shearing rates, we discuss the expected impact of the pitch-dependent growth rates for magnetospheric dissipation in magnetar conditions.
Abstract
One scenario for the generation of fast radio bursts (FRBs) is magnetic reconnection in a current sheet of the magnetar wind. Compressed by a strong magnetic pulse induced by a magnetar ...flare, the current sheet fragments into a self-similar chain of magnetic islands. Time-dependent plasma currents at their interfaces produce coherent radiation during their hierarchical coalescence. We investigate this scenario using 2D radiative relativistic particle-in-cell simulations to compute the efficiency of the coherent emission and to obtain frequency scalings. Consistent with expectations, a fraction of the reconnected magnetic field energy,
f
∼ 0.002, is converted to packets of high-frequency fast magnetosonic waves, which can escape from the magnetar wind as radio emission. In agreement with analytical estimates, we find that magnetic pulses of 10
47
erg s
−1
can trigger relatively narrowband GHz emission with luminosities of approximately 10
42
erg s
−1
, sufficient to explain bright extragalactic FRBs. The mechanism provides a natural explanation for a downward frequency drift of burst signals, as well as the ∼100 ns substructure recently detected in
FRB 20200120E
.
Collisionless shocks can be produced as a result of strong magnetic fields in a plasma flow, and therefore are common in many astrophysical systems. The Weibel instability is one candidate mechanism ...for the generation of sufficiently strong fields to create a collisionless shock. Despite their crucial role in astrophysical systems, observation of the magnetic fields produced by Weibel instabilities in experiments has been challenging. Using a proton probe to directly image electromagnetic fields, we present evidence of Weibel-generated magnetic fields that grow in opposing, initially unmagnetized plasma flows from laser-driven laboratory experiments. Three-dimensional particle-in-cell simulations reveal that the instability efficiently extracts energy from the plasma flows, and that the self-generated magnetic energy reaches a few percent of the total energy in the system. This result demonstrates an experimental platform suitable for the investigation of a wide range of astrophysical phenomena, including collisionless shock formation in supernova remnants, large-scale magnetic field amplification, and the radiation signature from gamma-ray bursts.
We explore the physics of shock evolution and particle acceleration in non-relativistic collisionless shocks using hybrid simulations. We analyze a wide range of physical parameters relevant to the ...acceleration of cosmic rays (CRs) in astrophysical shock scenarios. We show that there are fundamental differences between high and low Mach number shocks in terms of the electromagnetic turbulence generated in the pre-shock zone; dominant modes are resonant with the streaming CRs in the low Mach number regime, while both resonant and non-resonant modes are present for high Mach numbers. Energetic power-law tails for ions in the downstream plasma account for up to 15% of the incoming upstream flow energy, distributed over ~5% of the particles in a power law with slope --2 ? 0.2 in energy. Quasi-parallel shocks with Delta *c <= 45? are good ion accelerators, while power laws are greatly suppressed for quasi-perpendicular shocks, Delta *c > 45?. The efficiency of conversion of flow energy into the energy of accelerated particles peaks at Delta *c = 15?-30? and MA = 6, and decreases for higher Mach numbers, down to ~2% for MA = 31. Accelerated particles are produced by diffusive shock acceleration (DSA) and by shock drift acceleration (SDA) mechanisms, with the SDA contribution to the overall energy gain increasing with magnetic inclination. We also present a direct comparison between hybrid and fully kinetic particle-in-cell results at early times. In supernova remnant (SNR) shocks, particle acceleration will be significant for low Mach number quasi-parallel flows (MA < 30, Delta *c < 45). This finding underscores the need for an effective magnetic amplification mechanism in SNR shocks.
Revised Pulsar Spin-down Contopoulos, I; Spitkovsky, A
The Astrophysical journal,
06/2006, Letnik:
643, Številka:
2
Journal Article
Recenzirano
Odprti dostop
We address the issue of electromagnetic pulsar spin-down by combining our experience from the two limiting idealized cases that have been studied thoroughly in the past: that of an aligned rotator, ...in which ideal MHD conditions apply, and that of a misaligned rotator in vacuum. We construct a spin-down formula that takes into account the misalignment of the magnetic and rotation axes, and the magnetospheric particle acceleration gaps. We show that near the death line, aligned rotators spin down more slowly than orthogonal ones. In order to test this approach, we use a simple Monte Carlo method to simulate the evolution of pulsars and find a good fit to the observed pulsar distribution in the P- super(P) diagram without invoking magnetic field decay. Our model may also account for individual pulsars spinning down with braking index n < 3, by allowing the corotating part of the magnetosphere to end inside the light cylinder. We discuss the role of magnetic reconnection in determining the pulsar braking index. We show, however, that n 6 3 remains a good approximation for the pulsar population as a whole. Moreover, we predict that pulsars near the death line have braking index values n > 3, and that the older pulsar population has preferentially smaller magnetic inclination angles. We discuss possible signatures of such alignment in the existing pulsar data.
Abstract
The ability of collisionless shocks to efficiently accelerate nonthermal electrons via diffusive shock acceleration (DSA) is thought to require an injection mechanism capable of ...preaccelerating electrons to high enough energy where they can start crossing the shock front potential. We propose, and show via fully kinetic plasma simulations, that in high-Mach-number shocks electrons can be effectively injected by scattering in kinetic-scale magnetic turbulence produced near the shock transition by the ion Weibel, or current filamentation, instability. We describe this process as a modified DSA mechanism where initially thermal electrons experience the flow velocity gradient in the shock transition and are accelerated via a first-order Fermi process as they scatter back and forth. The electron energization rate, diffusion coefficient, and acceleration time obtained in the model are consistent with particle-in-cell simulations and with the results of recent laboratory experiments where nonthermal electron acceleration was observed. This injection model represents a natural extension of DSA and could account for electron injection in high-Mach-number astrophysical shocks, such as those associated with young supernova remnants and accretion shocks in galaxy clusters.