Related papers: A simple procedure for generating a Kappa distribu…
A procedure for loading particle velocities from a relativistic kappa distribution in particle-in-cell (PIC) and Monte Carlo simulations is presented. It is based on the rejection method and the beta prime distribution. The rejection part…
A random number generator for the Kappa velocity distribution in particle simulations is proposed. Approximating the cumulative distribution function with the q-exponential function, an inverse transform procedure is constructed. The…
Simulating from a gamma distribution with small shape parameter is a challenging problem. Towards an efficient method, we obtain a limiting distribution for a suitably normalized gamma distribution when the shape parameter tends to zero.…
The Gamma distribution is well-known and widely used in many signal processing and communications applications. In this letter, a simple and extremely efficient accept/reject algorithm is introduced for the generation of independent random…
Kappa-distributed velocities in plasmas are common in a wide variety of settings, from low-density to high-density plasmas. To date, they have been found mainly in space plasmas, but are recently being considered also in the modelling of…
Kappa distributions are widely used in space plasma physics to model velocity distribution functions with heavy tails. Parameter estimation in these distributions is, however, complicated by the fact that the kappa distribution does not…
Numerical procedures for generating non-Maxwellian velocity distributions in particle simulations are presented. First, Monte Carlo methods for an $(r,q)$ distribution that generalizes flattop and kappa distributions are discussed. Then,…
Algorithms for generating random numbers that follow a gamma distribution with shape parameter less than unity are proposed. Acceptance-rejection algorithms are developed, based on the generalized exponential distribution. The squeeze…
In this paper we present a method to generate independent samples for a general random variable, either continuous or discrete. The algorithm is an extension of the acceptance-rejection method, and it is particularly useful for kinetic…
We demonstrate that a nonthermal distribution of particles described by a kappa distribution can be accurately approximated by a weighted sum of Maxwell-Boltzmann distributions. We apply this method to modeling collision processes in…
The kappa distribution of velocities appears routinely in the study of collisionless plasmas present in Earth's magnetosphere, the solar wind among other contexts where particles are unable to reach thermal equilibrium. Originally justified…
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…
In the literature different so-called $\kappa$-distribution functions are discussed to fit and model the velocity (or energy) distributions of solar wind species, pickup ions or magnetospheric particles. Here we introduce a generalized…
Particle acceleration in collisionless plasma systems is a central question in astroplasma and astroparticle physics. The structure of the acceleration regions, electron-ion energy equilibration, preacceleration of particles at shocks to…
In particle-in-cell simulations, excessive or even unfeasible computational demands can be caused by the growth of the number of particles in the course of prolific ionization or cascaded pair production due to the effects of quantum…
From the perspective of non-equilibrium statistical mechanics, modeling the velocity distribution of particles in non-equilibrium, steady-state plasmas presents a significant challenge. Under this context, a family of kappa distributions…
The use of explicit particle-in-cell (PIC) method for relativistic plasma simulations is restricted by numerical heating and instabilities that may significantly constrain the choice of time and space steps. To partially eliminate these…
We present the BACCO project, a simulation framework specially designed to provide highly-accurate predictions for the distribution of mass, galaxies, and gas as a function of cosmological parameters. In this paper, we describe our main…
Generation of pseudorandom numbers from different probability distributions has been studied extensively in the Monte Carlo simulation literature. Two standard generation techniques are the acceptance-rejection and inverse transformation…
An integrated simulation approach fully based upon particle-in-cell (PIC) model is proposed, which involves both fast particle generation via laser solid-density plasma interaction and transport and energy deposition of the particles in…