Search

An oblique-rotating pulsar generates a relativistic striped wind in a pulsar wind nebula (PWN). The termination shock of the PWN compresses the Poynting-flux-dominated flow and drives magnetic reconnection. By carrying out particle-in-cell (PIC) simulations of the termination shock of the PWN, we study the shock structure as well as the energy conversion processes and particle acceleration mechanisms. With the recent advances in the numerical methods, we extend the simulations to the ultra-relativistic regime with bulk Lorentz factor up to \gamma_{0}=10^{6}. Magnetic reconnection at the termination shock is highly efficient at converting magnetic energy to particle kinetic energy and accelerating particles to high energies. We find that the resulting energy spectra crucially depend on \lambda/d_{e}. When \lambda/d_{e} is large (\lambda\gtrsim40d_{e}) , the downstream particle spectra form a power-law distribution in the magnetically dominated relativistic wind regime with upstream magnetization parameter \sigma_{0}=10. By analyzing particle trajectories and statistical quantities relevant to particle energization, we find that Fermi-type mechanism dominates the particle acceleration and power-law formation. We find that the results for particle acceleration are scalable as \gamma_{0} and \sigma_{0} increase to large values. The maximum energy for electrons and positrons can reach hundreds of TeV if the wind has a bulk Lorentz factor \gamma_{0}\approx10^{6} and magnetization parameter \sigma_{0}=10, which can explain the recent observations of high-energy gamma-rays from pulsar wind nebulae (PWNe)., Comment: 17 pages, 10 figures, submitted to The Astrophysical Journal

A rotating pulsar creates a surrounding pulsar wind nebula (PWN) by steadily releasing an energetic wind into the interior of the expanding shockwave of supernova remnant or interstellar medium. At the termination shock of a PWN, the Poynting-flux-dominated relativistic striped wind is compressed. Magnetic reconnection is driven by the compression and converts magnetic energy into particle kinetic energy and accelerating particles to high energies. We carrying out particle-in-cell (PIC) simulations to study the shock structure as well as the energy conversion and particle acceleration mechanism. By analyzing particle trajectories, we find that many particles are accelerated by Fermi-type mechanism. The maximum energy for electrons and positrons can reach hundreds of TeV., Comment: 6 pages, 2 figures. 49th International Symposium on Multiparticle Dynamics (ISMD 2019)

Three-dimensional FLASH radiation-magnetohydrodynamics (radiation-MHD) modeling is carried out to study the hydrodynamics and magnetic fields in the shock-shear derived platform. Simulations indicate that fields of tens of Tesla can be generated via Biermann battery effect due to vortices and mix in the counter-propagating shock-induced shear layer. Synthetic proton radiography simulations using MPRAD and synthetic X-ray image simulations using SPECT3D are carried out to predict the observable features in the diagnostics. Quantifying the effects of magnetic fields in inertial confinement fusion (ICF) and high-energy-density (HED) plasmas represents frontier research that has far-reaching implications in basic and applied sciences., Comment: 14 pages, 7 figures

Proton radiography is used in various high-energy-density (HED) plasma experiments. In this paper, we describe a Monte Carlo and ray-tracing simulation tool called MPRAD that can be used for modeling the deflection of proton beams in arbitrary three dimensional electromagnetic fields, as well as the diffusion of the proton beams by Coulomb scattering and stopping power. The Coulomb scattering and stopping power models in cold matter and fully ionized plasma are combined using interpolation. We discuss the application of MPRAD in a few setups relevant to HED plasma experiments where the plasma density can play a role in diffusing the proton beams and affecting the prediction and interpretation of the proton images. It is shown how the diffusion due to plasma density can affect the resolution and dynamical range of the proton radiography., Comment: 11 pages, 4 figures

The WT scheme, a piecewise polynomial force interpolation scheme with time-step dependency, is proposed in this paper for relativistic particle-in-cell (PIC) simulations. The WT scheme removes the lowest order numerical Cherenkov instability (NCI) growth rate for arbitrary time steps allowed by the Courant condition. While NCI from higher order resonances is still present, the numerical tests show that for smaller time steps, the numerical instability grows much slower than using the optimal time step found in previous studies. The WT scheme is efficient for improving the quality and flexibility of relativistic PIC simulations., Comment: 11 pages, single column, one figure

Inertial confinement fusion fuel suffers increased X-ray radiation losses when carbon from the capsule ablator mixes into the hot-spot. Here we present one and two-dimensional ion Vlasov-Fokker-Planck simulations that resolve hot-spot self heating in the presence of a localized spike of carbon mix, totalling 1.9% of the hot-spot mass. The mix region cools and contracts over tens of picoseconds, increasing its alpha particle stopping power and radiative losses. This makes a localized mix region more severe than an equal amount of uniformly distributed mix. There is also a purely kinetic effect that reduces fusion reactivity by several percent, since faster ions in the tail of the distribution are absorbed by the mix region. Radiative cooling and contraction of the spike induces fluid motion, causing neutron spectrum broadening. This artificially increases the inferred experimental ion temperatures and gives line of sight variations.

A Maxwell solver derived from finite element method with \mathcal{O}(N) computing cost is developed to improve the numerical dispersion properties in relativistic particle-in-cell (PIC) simulations. The correction of the dispersion relation of the electromagnetic wave is achieved using the neighboring cells via an iteration scheme without decomposing into Fourier modes. The local nature of the communication is ideally suited to massively parallel computer architectures. This Maxwell solver constrains the Numerical Cherenkov instability (NCI) for the ultra-relativistic drifting pair plasma in x direction to large wave vectors for two dimensional grid. The growth rate of NCI is suppressed by using the low pass filtering., Comment: 16 pages, 4 figures, submitted to Computer Physics Communications