Related papers: Diffusive Shock Acceleration Simulations: Comparis…
The acceleration of energetic particle transport in high amplitude magnetosonic and Alfvenic turbulence is considered using the method of Monte Carlo particle simulations, involving integration of particle equations of motion. We derive the…
In this paper, a diffusion-aggregation equation with delta potential is introduced. Based on the global existence and uniform estimates of solutions to the diffusion-aggregation equation, we also provide the rigorous derivation from a…
A straightforward analytical scheme is proposed for computing the long-time, asymptotic mean velocity and dispersivity (effective diffusivity) of a particle undergoing a discrete biased random walk on a periodic lattice amongst an array of…
This topical review describes the methodology of continuum variational and diffusion quantum Monte Carlo calculations. These stochastic methods are based on many-body wave functions and are capable of achieving very high accuracy. The…
Collisionless shocks are efficient particle accelerators. At Earth, ions with energies exceeding 100 keV are seen upstream of the bow shock when the magnetic geometry is quasi-parallel, and large-scale supernova remnant shocks can…
We use kinetic hybrid simulations (kinetic ions - fluid electrons) to characterize the fraction of ions that are accelerated to non-thermal energies at non-relativistic collisionless shocks. We investigate the properties of the shock…
For problems in astrophysics, planetary science and beyond, numerical simulations are often limited to simulating fewer particles than in the real system. To model collisions, the simulated particles (aka superparticles) need to be inflated…
We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust,…
We present an efficient Monte Carlo method to simulate reaction-diffusion processes with spatially varying particle annihilation or transformation rates as it occurs for instance in the context of motor-driven intracellular transport. Like…
We present various techniques that make orbit-following Monte Carlo simulations faster and more reliable when assessing collisionless fast particle losses due to magnetic field perturbations. These techniques are based on identifying…
We use kinetic simulations of diffusive shock acceleration (DSA) to study the time-dependent evolution of plane, quasi-parallel, cosmic-ray (CR) modified shocks. Thermal leakage injection of low energy CRs and finite Alfv\'en wave…
An efficient Path Integral Monte Carlo procedure is proposed to simulate the behavior of quantum many-body dissipative systems described within the framework of the influence functional. Thermodynamic observables are obtained by Monte Carlo…
We present simulations of magnetized astrophysical shocks taking into account the interplay between the thermal plasma of the shock and supra-thermal particles. Such interaction is depicted by combining a grid-based magneto-hydrodynamics…
A shock wave propagating perpendicularly to an ambient magnetic field accelerates particles considerably faster than in the parallel propagation regime. However, the perpendicular acceleration stops after the shock overruns a circular…
Computational chemistry allows researchers to experiment in sillico: by running a computer simulations of a biological or chemical processes of interest. Molecular dynamics with molecular mechanics model of interactions simulates N-body…
Finding the injection threshold for diffusive shock acceleration (DSA) of electrons in collisionless shocks has been a longstanding unsolved problem. Using first-principles kinetic simulations, we identify the conditions for electron…
Advanced algorithms are necessary to obtain faster-than-real-time dynamic simulations in a number of different physical problems that are characterized by widely disparate time scales. Recent advanced dynamic Monte Carlo algorithms that…
The spectral index $s$ of particles diffusively accelerated in a relativistic shock depends on the unknown angular diffusion function $\mathcal{D}$, which itself depends on the particle distribution function $f$ if acceleration is…
Cosmic Ray transport in curved background magnetic fields is investigated using numerical Monte-Carlo simulation techniques. Special emphasis is laid on the Solar system, where the curvature of the magnetic field can be described in terms…
We study the physics of electron acceleration at collisionless shocks that move through a plasma containing large-scale magnetic fluctuations. We numerically integrate the trajectories of a large number of electrons, which are treated as…