Related papers: Loading a relativistic kappa distribution in parti…
Generating random variates from high-dimensional distributions is often done approximately using Markov chain Monte Carlo. In certain cases, perfect simulation algorithms exist that allow one to draw exactly from the stationary…
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…
Rejection Sampling is a fundamental Monte-Carlo method. It is used to sample from distributions admitting a probability density function which can be evaluated exactly at any given point, albeit at a high computational cost. However,…
Recently-proposed particle MCMC methods provide a flexible way of performing Bayesian inference for parameters governing stochastic kinetic models defined as Markov (jump) processes (MJPs). Each iteration of the scheme requires an estimate…
We provide an algorithm for generation of momenta (or energies) of relativistic particles according to the relativistic Bose-Einstein or Fermi-Dirac distributions. The algorithm uses rejection method with effectively selected comparison…
A number of distributions that arise in statistical applications can be expressed in the form of a weighted density: the product of a base density and a nonnegative weight function. Generating variates from such a distribution may be…
Particle velocity distribution functions (VDF) in space plasmas often show non Maxwellian suprathermal tails decreasing as a power law of the velocity. Such distributions are well fitted by the so-called Kappa distribution. The presence of…
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…
In this paper, a mathematical formulation of the probabilistic available transfer capability (PATC) problem is proposed to incorporate uncertainties from the large-scale renewable energy generation (e.g., wind farms and solar PV power…
Light propagation in turbid media is driven by the equation of radiative transfer. We give a formal probabilistic representation of its solution in the framework of biological tissues and we implement algorithms based on Monte Carlo methods…
An extension of Maxwell's original prescription for an ideal gas is adopted to derive a broad class of Kappa-type velocity distributions, encompassing both fat and short-tailed forms. Within this general framework, a physically consistent…
Random sequential adsorption is an irreversible surface deposition of extended objects. In systems with continuous degrees of freedom coverage follows a power law, theta(t) = theta_J - c t^{-alpha}, where the exponent alpha depends on the…
We analyze here in some detail, the derivation of the Particle and Monte Carlo methods of plasma simulation, such as Particle in Cell (PIC), Monte Carlo (MC) and Particle in Cell / Monte Carlo (PIC/MC) from formal manipulation of transport…
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…
The generally accepted representation of $\kappa$-distributions in space plasma physics allows for two different alternatives, namely assuming either the temperature or the thermal velocity to be $\kappa$-independent. The present paper aims…
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 find the probability density function $\mathcal{P}(V_{\texttt{r}})$ of the relativistic relative velocity for two colliding particles in a non-degenerate relativistic gas. The distribution reduces to Maxwell distribution for the relative…
Rejection sampling is a technique for sampling from difficult distributions. However, its use is limited due to a high rejection rate. Common adaptive rejection sampling methods either work only for very specific distributions or without…
An algorithm is demonstrated that performs first-principles tracking of relativistic charged-particles. A covariant approach is used which relies on retarded vector potentials for trajectory integration instead of performing electromagnetic…
We develop Monte Carlo methods for sampling random states and corresponding bit strings in qubit systems. To this end, we derive exact probability density functions that yield the Porter-Thomas distribution in the limit of large systems. We…