Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userParticle transport in hybrid PIC shock simulations: A comparison of diagnostics
119
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
Recent in-situ and remote observations suggest that the transport regime associated with shock accelerated particles may be anomalous {i.e., the Mean Square Displacement (MSD) of such particles scales non-linearly with time}. We use self-consistent, hybrid PIC plasma simulations to simulate a quasi-parallel shock with parameters compatible with heliospheric shocks, and gain insights about the particle transport in such a system. For suprathermal particles interacting with the shock we compute the MSD separately in the upstream and downstream regions. Tracking suprathermal particles for sufficiently long times up and/or downstream of the shock poses problems in particle plasma simulations, such as statistically poor particle ensembles and trajectory fragments of variable length in time. Therefore, we introduce the use of time-averaged mean square displacement (TAMSD), which is based on single particle trajectories, as an additional technique to address the transport regime for the upstream and downstream regions. MSD and TAMSD are in agreement for the upstream energetic particle population, and both give a strong indication of superdiffusive transport, consistent with interplanetary shock observations. MSD and TAMSD are also in reasonable agreement downstream, where indications of anomalous transport are also found. TAMSD shows evidence of heterogeneity in the diffusion properties of the downstream particle population, ranging from subdiffusive behaviour of particles trapped in the strong magnetic field fluctuations generated at the shock, to superdiffusive behaviour of particles transmitted and moving away from the shock.
rate research
Read More
In this paper, by performing a two-dimensional particle-in-cell simulation, we investigate magnetic reconnection in the downstream of a quasi-perpendicular shock. The shock is nonstationary, and experiences a cyclic reformation. At the beginning of reformation process, the shock front is relatively flat, and part of upstream ions are reflected by the shock front. The reflected ions move upward in the action of Lorentz force, which leads to the upward bending of magnetic field lines at the foot of the shock front, and then a current sheet is formed due to the squeezing of the bending magnetic field lines. The formed current sheet is brought toward the shock front by the solar wind, and the shock front becomes irregular after interacting with the current sheet. Both the current sheet brought by the solar wind and the current sheet associated with the shock front are then fragmented into many small filamentary current sheets. Electron-scale magnetic reconnection may occur in several of these filamentary current sheets when they are convected into the downstream, and magnetic islands are generated. A strong reconnection electric field and energy dissipation are also generated around the X line, and high-speed electron outflow is also formed.
A shock tube problem is solved numerically by using one-dimensional full particle-in-cell simulations under the condition that a relatively tenuous and weakly magnetized plasma is continuously pushed by a relatively dense and strongly magnetized plasma having supersonic relative velocity. A forward and a reverse shock and a contact discontinuity are self-consistently reproduced. The spatial width of the contact discontinuity increases as the angle between the discontinuity normal and ambient magnetic field decreases. The inner structure of the discontinuity shows different profiles between magnetic field and plasma density, or pressure, which is caused by a non-MHD effect of the local plasma. The region between the two shocks is turbulent. The fluctuations in the relatively dense plasma are compressible and propagating away from the contact discontinuity, although the fluctuations in the relatively tenuous plasma contain both compressible and incompressible components. The source of the compressible fluctuations in the relatively dense plasma is in the relatively tenuous plasma. Only compressible fast mode fluctuations generated in the relatively tenuous plasma are transmitted through the contact discontinuity and propagate in the relatively dense plasma. These fast mode fluctuations are steepened when passing the contact discontinuity. This wave steepening and probably other effects may cause the broadening of the wave spectrum in the very local interstellar medium plasma. The results are discussed in the context of the heliospheric boundary region or heliopause.
Cosmic rays (CRs) have critical impacts in the multiphase interstellar medium (ISM), driving dynamical motions in low-density plasma and modifying the ionization state, temperature, and chemical composition of higher-density atomic and molecular gas. We present a study of CR propagation between the ionized ISM and a neutral cloud. Using one-dimensional magnetohydrodynamic particle-in-cell simulations which include ion-neutral drag to damp Alfv$acute{text{e}}$n waves in the cloud, we self-consistently evolve the kinetic physics of CRs and fluid dynamics of the multiphase gas. By introducing the cloud in our periodic domain, our simulations break translational symmetry and allow the emergence of spatial structure in the CR distribution function. A negative spatial gradient forms across the fully-ionized ISM region while a positive gradient forms across the neutral cloud. We connect our results with CR hydrodynamics formulations by computing the wave-particle scattering rates as predicted by quasilinear, fluid, and Fokker-Planck theory. For momenta where the mean free path is short relative to the box size, we find excellent agreement among all scattering rates. By exploring different cloud sizes and ion-neutral collision rates, we show that our results are robust. Our work provides a first-principles verification of CR hydrodynamics when particles stream down their pressure gradient, and opens a pathway toward comprehensive calibrations of transport coefficients from self-generated Alfv$acute{text{e}}$n wave scattering with CRs.
Non-thermal pickup ions (PUIs) are created in the solar wind (SW) by charge-exchange between SW ions (SWIs) and slow interstellar neutral atoms. It has long been theorized, but not directly observed, that PUIs should be preferentially heated at quasi-perpendicular shocks compared to thermal SWIs. We present in situ observations of interstellar hydrogen (H+) PUIs at an interplanetary shock by the New Horizons Solar Wind Around Pluto (SWAP) instrument at ~34 au from the Sun. At this shock, H+ PUIs are only a few percent of the total proton density but contain most of the internal particle pressure. A gradual reduction in SW flow speed and simultaneous heating of H+ SWIs is observed ahead of the shock, suggesting an upstream energetic particle pressure gradient. H+ SWIs lose ~85% of their energy flux across the shock and H+ PUIs are preferentially heated. Moreover, a PUI tail is observed downstream of the shock, such that the energy flux of all H+ PUIs is approximately six times that of H+ SWIs. We find that H+ PUIs, including their suprathermal tail, contain almost half of the total downstream energy flux in the shock frame.
The Earths magnetotail is characterized by stretched magnetic field lines. Energetic particles are effectively scattered due to the field-line curvature, which then leads to isotropization of energetic particle distributions and particle precipitation to the Earths atmosphere. Measurements of these precipitation at low-altitude spacecraft are thus often used to remotely probe the magnetotail current sheet configuration. This configuration may include spatially localized maxima of the curvature radius at equator (due to localized humps of the equatorial magnetic field magnitude) that reduce the energetic particle scattering and precipitation. Therefore, the measured precipitation patterns are related to the spatial distribution of the equatorial curvature radius that is determined by the magnetotail current sheet configuration. In this study, we show that, contrary to previous thoughts, the magnetic field line configuration with the localized curvature radius maximum can actually enhance the scattering and subsequent precipitation. The spatially localized magnetic field dipolarization (magnetic field humps) can significantly curve magnetic field lines far from the equator and create off-equatorial minima in the curvature radius. Scattering of energetic particles in these off-equatorial regions alters the scattering (and precipitation) patterns, which has not been studied yet. We discuss our results in the context of remote-sensing the magnetotail current sheet configuration with low-altitude spacecraft measurements.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
