Related papers: Hybrid simulation of energetic particles interacti…
In partially ionised plasmas, the magnetic field can become decoupled from the neutral gas and diffuse through it in a process known as ambipolar diffusion. Although ambipolar diffusion has been implemented in several grid codes, we here…
We present particle-in-cell simulations of one dimensional relativistic electromagnetic shocks in a uniform magnetic field, for a range of magnetic field strengths, plasma temperatures and numerical initial conditions. We show that the…
Our focus is on simulating the dynamics of non-interacting particles including the effects of an external potential, which, under certain assumptions, can be formally described by the Dean-Kawasaki equation. The Dean-Kawasaki equation can…
We investigate a hybrid numerical algorithm aimed at the large-scale cosmological N-body simulation for the on-going and the future high precious sky surveys. It makes use of a truncated Fast Multiple Method (FMM) for short-range gravity,…
We discuss the fundamentals of the implicit moment method for Particle In Cell (PIC) simulation as presently implemented in the CELESTE3D code. We present the method in its fully electromagnetic and fully kinetic version. The application of…
Collisionless, magnetized turbulence offers a promising framework for the generation of non-thermal high-energy particles in various astrophysical sites. Yet, the detailed mechanism that governs particle acceleration has remained subject to…
Numerical studies of shock waves in large scale systems via kinetic simulations with millions of particles are too computationally demanding to be processed in serial. In this work we focus on optimizing the parallel performance of a…
We present a super-polynomial improvement in the precision scaling of quantum simulations for coupled classical-quantum systems in this paper. Such systems are found, for example, in molecular dynamics simulations within the…
This paper describes a new code for simulating astrophysical plasmas that solves a hybrid model composed of gyrokinetic ions (GKI) and an isothermal electron fluid (ITEF) [A. Schekochihin et al., Astrophys. J. Suppl. \textbf{182}, 310…
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…
The implicit 2D3V particle-in-cell (PIC) code developed to study the interaction of ultrashort pulse lasers with matter [G. M. Petrov and J. Davis, Computer Phys. Comm. 179, 868 (2008); Phys. Plasmas 18, 073102 (2011)] has been parallelized…
A massively parallel simulation code, called \textit{dHybrid}, has been developed to perform global scale studies of space plasma interactions. This code is based on an explicit hybrid model; the numerical stability and parallel scalability…
We present a concurrent Monte Carlo (MC) - molecular dynamics (MD) approach to modeling of matter response to excitation of its electronic system. The two methods are combined on-the-fly at each time step in one code, TREKIS-4. The MC model…
The role of MHD turbulence in the cosmic ray acceleration process in a volume with a reconnecting magnetic field is studied by means of Monte Carlo simulations. We performed modelling of proton acceleration with the 3D analytic model of…
Quantum computers have the potential to simulate chemical systems beyond the capability of classical computers. Recent developments in hybrid quantum-classical approaches enable the determinations of the ground or low energy states of…
We present preliminary performance results of gPLUTO, the new GPU-optimized implementation of the PLUTO code for computational plasma astrophysics. Like its predecessor, gPLUTO employs a finite-volume formulation to numerically solve the…
This paper investigates hybrid kinetic-MHD models, where a hot plasma (governed by a kinetic theory) interacts with a fluid bulk (governed by MHD). Different nonlinear coupling schemes are reviewed, including the pressure-coupling scheme…
Multi-particle collision dynamics (MPCD) is a mesoscopic simulation method to simulate fluid particle-like flows. MPCD has been widely used to simulate various problems in condensed matter. In this study, hydrostatic behavior of gas in the…
Particle accelerator modeling is an important field of research and development, essential to investigating, designing and operating some of the most complex scientific devices ever built. Kinetic simulations of relativistic, charged…
We report on the development of a test-particle based kinetic Monte Carlo code for large systems and its application to simulate matter in the continuum regime. Our code combines advantages of the Direct Simulation Monte Carlo and the…