Related papers: Loading a relativistic kappa distribution in parti…
Likelihood-free methods, such as approximate Bayesian computation, are powerful tools for practical inference problems with intractable likelihood functions. Markov chain Monte Carlo and sequential Monte Carlo variants of approximate…
The leapfrog integrator is routinely used within the Hamiltonian Monte Carlo method and its variants. We give strong numerical evidence that alternative, easy to implement algorithms yield fewer rejections with a given computational effort.…
Thermal Doppler broadening of spectral profiles for particle populations in the absence or presence of potential fields are described by kappa distributions. The kappa distribution provides a replacement for the Maxwell-Boltzmann…
The kappa-distributed fully ionized plasma with collisional interaction is investigated. The Fokker-Planck equation with Rosenbluth potential is employed to describe such a physical system. The results show that the kappa distribution is…
We discuss non-reversible Markov-chain Monte Carlo algorithms that, for particle systems, rigorously sample the positional Boltzmann distribution and that have faster than physical dynamics. These algorithms all feature a non-thermal…
Rejection sampling is a popular method used to generate numbers that follow some given distribution. We study the use of this method to generate random numbers in the unit interval from increasing probability density functions. We focus on…
A method for generating random $U(1)$ variables with Boltzmann distribution is presented. It is based on the rejection method with transformation of variables. High efficiency is achieved for all range of temparatures or coupling…
In relativistic kinetic theory, which underlies relativistic hydrodynamics, the molecular chaos hypothesis stands at the basis of the equilibrium Maxwell-Juttner probability distribution for the four-momentum $p^{\alpha}$. We investigate…
Modifications to Swisdak [Phys. Plasmas 20, 062110 (2013)]'s rejection sampling algorithm for drawing a Maxwell-J\"{u}ttner distribution in particle simulations are presented. Handy approximations for $e$-folding points and a linear slope…
Acceptance-rejection (AR), Independent Metropolis Hastings (IMH) or importance sampling (IS) Monte Carlo (MC) simulation algorithms all involve computing ratios of probability density functions (pdfs). On the other hand, classifiers…
We investigate the use of the Metropolis-Hastings algorithm to sample posterior distribution in a Bayesian inverse problem, where the likelihood function is random. Concretely, we consider the case where one has full field observations of a…
One of the most demanding calculations is to generate random samples from a specified probability distribution (usually with an unknown normalizing prefactor) in a high-dimensional configuration space. One often has to resort to using a…
We revisit effective scenarios for the origin of heavy tails in stationary velocity distributions. A first analysis combines localization with diffusive acceleration. That gets realized in space plasmas to find the so-called…
In this article, we proposed a new probability distribution named as power Maxwell distribution (PMaD). It is another extension of Maxwell distribution (MaD) which would lead more flexibility to analyze the data with non-monotone failure…
The departure of particle distributions from the Maxwellian is commonly observed in space plasmas. These non-Maxwellian distributions which are typical for plasmas that are not in thermal equilibrium, can be modeled with…
The probability distribution of the order parameter is exploited in order to obtain the criticality of magnetic systems. Monte Carlo simulations have been employed by using single spin flip Metropolis algorithm aided by finite-size scaling…
The system-level dynamics of multivalent biomolecular interactions can be simulated using a rule-based kinetic Monte Carlo method in which a rejection sampling strategy is used to generate reaction events. This method becomes inefficient…
Principal component analysis (PCA) is one of the most popular dimension reduction techniques in statistics and is especially powerful when a multivariate distribution is concentrated near a lower-dimensional subspace. Multivariate extreme…
A Monte Carlo model of radiative transfer in multi-phase dusty media is applied to the situation of stars and clumpy dust in a sphere or a disk. The distribution of escaping and absorbed photons are shown for various filling factors and…
We consider the combined use of resampling and partial rejection control in sequential Monte Carlo methods, also known as particle filters. While the variance reducing properties of rejection control are known, there has not been (to the…