Related papers: Loading a relativistic kappa distribution in parti…
Propagation of charged cosmic-rays in the Galaxy depends on the transport parameters, whose number can be large depending on the propagation model under scrutiny. A standard approach for determining these parameters is a manual scan,…
We analyze the line radiative transfer in protoplanetary disks using several approximate methods and a well-tested Accelerated Monte Carlo code. A low-mass flaring disk model with uniform as well as stratified molecular abundances is…
The time of flight distribution for a cloud of cold atoms falling freely under gravity is considered. We generalise the probability current density approach to calculate the quantum arrival time distribution for the mixed state describing…
We propose a new method, probabilistic divide-and-conquer, for improving the success probability in rejection sampling. For the example of integer partitions, there is an ideal recursive scheme which improves the rejection cost from…
For a wide class of applications of the Monte Carlo method, we describe a general sampling methodology that is guaranteed to converge to a specified equilibrium distribution function. The method is distinct from that of Metropolis in that…
The efficiency of Hamiltonian Monte Carlo (HMC) can suffer when sampling a distribution with a wide range of length scales, because the small step sizes needed for stability in high-curvature regions are inefficient elsewhere. To address…
Principal component analysis (PCA) is perhaps the most widely used method for data dimensionality reduction. A key question in PCA is deciding how many factors to retain. This manuscript describes a new approach to automatically selecting…
Models of stochastic processes are widely used in almost all fields of science. Theory validation, parameter estimation, and prediction all require model calibration and statistical inference using data. However, data are almost always…
A recent proposal to connect the loop quantization with the spin foam model for cosmology via the path integral is hereby adapted to the case of mechanical systems within the framework of the so called polymer quantum mechanics. The…
Relativistic magnetic reconnection has been proposed as an important nonthermal particle acceleration (NTPA) mechanism that generates power-law spectra and high-energy emissions. Power-law particle spectra are in general characterized by…
A stationary distribution function that describes the entire processes of propagation of relativistic particles, including the transition between the ballistic and diffusion regimes, is obtained. The spacial component of the constructed…
Synchrotron emission and absorption determine the observational appearance of many astronomical systems. In this paper, we describe a numerical scheme for calculating synchrotron emissivities and absorptivities in all four Stokes parameters…
Importance sampling is a Monte Carlo method which designs estimators of expectations under a target distribution using weighted samples from a proposal distribution. When the target distribution is complex, such as multimodal distributions…
Computing reaction rates in biomolecular systems is a common goal of molecular dynamics simulations. The reactions considered often involve conformational changes in the molecule, either changes in the structure of a protein or the relative…
Stochastic reaction-diffusion models are employed to represent many complex physical, biological, societal, and ecological systems. The macroscopic reaction rates describing the large-scale kinetics in such systems are effective,…
We study a generalization of the standard trapping problem of random walk theory in which particles move subdiffusively on a one-dimensional lattice. We consider the cases in which the lattice is filled with a one-sided and a two-sided…
In Bayesian inverse problems, the posterior distribution is used to quantify uncertainty about the reconstructed solution. In practice, Markov chain Monte Carlo algorithms often are used to draw samples from the posterior distribution.…
We consider a particle performing a stochastic motion on a one-dimensional lattice with jump widths distributed according to a power-law with exponent $\mu + 1$. Assuming that the walker moves in the presence of a distribution $a(x)$ of…
We propose a new algorithm for sampling the $N$-body density $|\Psi({\bf R})|^2/\int_{\mathbb{R}^{3N}} |\Psi|^2$ in the Variational Monte Carlo (VMC) framework. This algorithm is based upon a modified Ricci-Ciccotti discretization of the…
The Quasiparticle Random Phase Approximation (QRPA) is used in evaluation of the total muon capture ratesfor the final nuclei participating in double-beta decay. Several variants of the method are used, depending on the size of the single…