Related papers: A guiding center implementation for relativistic p…
We study a diffusion approximation for a model of stochastic motion of a particle in one spatial dimension. The velocity of the particle is constant but the direction of the motion undergoes random changes with a Poisson clock. Moreover,…
The relativistic spinning particle model, proposed in [3,4], is analyzed in a Hamiltonian framework. The spin is simulated by extending the configuration space by introducing a light-like four vector degree of freedom. The model is heavily…
Despite their discovery about half a century ago, the Gamma-ray burst (GRB) prompt emission mechanism is still not well understood. Theoretical modeling of the prompt emission has advanced considerably due to new computational tools and…
Context. Relativistic jets in active galactic nuclei are known for their exceptional energy output, and imaging the synthetic synchrotron emission of numerical jet simulations is essential for a comparison with observed jet polarization…
It is generally accepted that the dynamics of relativistic particles in the lab frame can be described by taking into account the relativistic dependence of the particles momenta on the velocity, with no reference to Lorentz…
The design and optimization of realistic architectures for fault-tolerant quantum computation requires error models that are both reliable and amenable to large-scale classical simulation. Perhaps the simplest and most practical…
At the core of some of the most important problems in plasma physics -- from controlled nuclear fusion to the acceleration of cosmic rays -- is the challenge to describe nonlinear, multi-scale plasma dynamics. The development of reduced…
Force-free electrodynamics describes the electromagnetic field of the magnetically dominated plasma found near pulsars and active black holes, but gives no information about the underlying particles that ultimately produce the observable…
Lattice spin models are useful for studying critical phenomena and allow the extraction of equilibrium and dynamical properties. Simulations of such systems are usually based on Monte Carlo (MC) techniques, and the main difficulty is often…
We consider Monte Carlo algorithms for the simulation of charged lattice gases with purely local dynamics. We study the mobility of particles as a function of temperature and show that the poor mobility of particles at low temperatures is…
A new adaptive scheme to be used in Particle-In-Cell codes for carrying out gyroaveraging operations with matrices is presented. This new scheme uses an intermediate velocity grid whose resolution is adapted to the local thermal Larmor…
Deterministically solving charged particle transport problems at a sufficient spatial and angular resolution is often prohibitively expensive, especially due to their highly forward peaked scattering. We propose a model order reduction…
We present a new general relativistic hydrodynamics code specifically designed to study magneto-rotational, relativistic, stellar core collapse. The code is an extension of an existing (and thoroughly tested) hydrodynamics code, which has…
For an axisymmetric tokamak plasma, Hamiltonian theory predicts that the orbits of charged particles must stay on invariant tori of conserved energy in the moving frame of reference of a wave that propagates along the torus with a fixed…
In this paper we present a new approach to the solution to a generalized version of Hughes' models for pedestrian movements based on a follow-the-leader many particle approximation. In particular, we provide a rigorous global existence…
Astrophysical observations suggest that magnetic reconnection in relativistic plasmas plays an important role in the acceleration of energetic particles. Modeling this accurately requires numerical schemes capable of addressing large scales…
The electromagnetic field is typically measured by the charged particle motion observation. Generally in the experiments, position, velocity and other physical parameters concerning relativistic particle beams, are estimated evaluating the…
We present the Fortran code SIMLA, which is designed for the study of charged particle dynamics in laser and other background fields. This can be done classically via the Landau-Lifshitz equation, or alternatively, via the simulation of…
We present a new implicit asymptotic preserving time integration scheme for charged-particle orbit computation in arbitrary electromagnetic fields. The scheme is built on the Crank-Nicolson integrator and continues to recover full-orbit…
We use the Los Alamos VPIC code to investigate particle acceleration in relativistic, unmagnetized, collisionless electron-ion plasmas. We run our simulations both with a realistic proton-to-electron mass ratio m_p/m_e = 1836, as well as…