Related papers: Test-Particle Diffusive Shock Acceleration in Rela…
We present results of test-particle simulations on both the first- and the second-order Fermi acceleration for relativistic parallel shock waves. Our studies suggest that the role of the second-order mechanism in the turbulent downstream of…
We present an extensive Monte Carlo study on light transport in optically thin slabs, addressing both axial and transverse propagation. We completely characterize the so-called ballistic-to-diffusive transition, notably in terms of the…
We consider a phenomenological model for the thermal leakage injection in the diffusive shock acceleration (DSA) process, in which suprathermal protons and electrons near the shock transition zone are assumed to have the so-called…
Energetic particles are ubiquitous in space and astrophysical plasmas, and interplanetary shocks are widely regarded as one of the main particle accelerators in the heliosphere. Indeed, in-situ measurements typically show that energetic…
A major attraction of diffusive shock acceleration is the prediction of power-law spectra for energetic particle distributions. However, this property is not fundamental to the theory. We demonstrate that for planar shocks with an oblique…
The aim is to assess the combined effect of diffusion and dispersion on shocks in the moderate dispersion regime. For a diffusive dispersive approximation of the equations of one-dimensional elasticity (or p-system), we study convergence of…
We present a novel way of performing kinetic Monte Carlo simulations which does not require an {\it a priori} list of diffusion processes and their associated energetics and reaction rates. Rather, at any time during the simulation,…
Diffusion-Limited Aggregation (DLA), the canonical model for non-equilibrium fractal growth, emerges from the simple rule of irreversible attachment by random walkers. Despite four decades of study, a unified computational framework…
We propose a multilevel Monte Carlo method for a particle-based asymptotic-preserving scheme for kinetic equations. Kinetic equations model transport and collision of particles in a position-velocity phase-space. With a diffusive scaling,…
One of the main features of astrophysical shocks is their ability to accelerate particles to extremely high energies. The leading acceleration mechanism, the diffusive shock acceleration is reviewed. It is demonstrated that its efficiency…
The hybrid Monte Carlo (HMC) algorithm is a ubiquitous method in computational physics with applications ranging from condensed matter to lattice QCD and beyond. However, HMC simulations often suffer from long autocorrelation times,…
We have developed a cosmic ray (CR) shock code in one dimensional spherical geometry with which the particle distribution, the gas flow and their nonlinear interaction can be followed numerically in a frame comoving with an expanding shock.…
Non-relativistic collisionless shock waves are widespread in space and astrophysical plasmas and are known as efficient particle accelerators. However, our understanding of collisionless shocks, including their structure and the mechanisms…
Energetic positrons have been observed in the interstellar medium, and high-energy positrons with relativistic energies up to approximately 1 TeV have been detected in Galactic cosmic rays. We conducted a study on the acceleration of…
Computer simulation plays a central role in modern day materials science. The utility of a given computational approach depends largely on the balance it provides between accuracy and computational cost. Molecular crystals are a class of…
Kinetic approaches are routinely employed to simulate the dynamics of systems that are too rarified to be described by the Navier-Stokes equations. However, generally they are far too computationally expensive to be applied for systems that…
Radio relics in galaxy clusters are giant diffuse synchrotron sources powered in cluster outskirts by merger shocks. Although the relic-shock connection has been consolidated in recent years by a number of observations, the details of the…
We developed a Monte Carlo simulation method to calculate incoherent Thomson scattering spectra in high temperature plasmas. The basic idea is to treat the entire scattering process as the superposition of individual photon-electron…
Radio relics in galaxy clusters are thought to be associated with powerful shock waves that accelerate particles via diffusive shock acceleration (DSA). Among the particles accelerated by DSA, relativistic protons should outnumber electrons…
Turbulent, relativistic nonthermal plasmas are ubiquitous in high-energy astrophysical systems, as inferred from broadband nonthermal emission spectra. The underlying turbulent nonthermal particle acceleration (NTPA) processes have…