Related papers: Loading a relativistic kappa distribution in parti…
For random variables produced through the inverse transform method, approximate random variables are introduced, which are produced by approximations to a distribution's inverse cumulative distribution function. These approximations are…
We consider the problem of learning a target probability distribution over a set of $N$ binary variables from the knowledge of the expectation values (with this target distribution) of $M$ observables, drawn uniformly at random. The space…
The various distribution functions can encompass the diverse characteristics of the magnetospheric plasma of surrounding neutron stars in both hot and cold environments; however, the Maxwell-J\"uttner distribution is so far widely used to…
Semi-linear response theory determines the absorption coefficient of a driven system using a resistor network calculation: Each unperturbed energy level of a particle in a vibrating trap, or of an electron in a mesoscopic ring, is regarded…
The non-Maxwellian $\kappa$-distributions have been detected in the solar transition region and flares. These distributions are characterized by a high-energy tail and a near-Maxwellian core and are known to have significant impact on the…
When it comes to active particles, even an ideal-gas model in a harmonic potential poses a mathematical challenge. An exception is a run-and-tumble model (RTP) in one-dimension for which a stationary distribution is known exactly. The case…
Monte Carlo simulations of systems with a complex action are known to be extremely difficult. A new approach to this problem based on a factorization property of distribution functions of observables has been proposed recently. The method…
Monte Carlo simulations are an essential tool in particle physics data analysis. Events are typically generated alongside weights that redistribute the cross section of the simulated process across the phase space. These weights can be…
Metastability is a common obstacle to performing long molecular dynamics simulations. Many numerical methods have been proposed to overcome it. One method is parallel replica dynamics, which relies on the rapid convergence of the underlying…
The quantum computer offers significant advantages in simulating physical systems, particularly those with exponentially large state spaces, such as quantum systems. Stochastic reaction-diffusion systems, characterized by their stochastic…
We find the laws for the spreading of the spatial widths (parallel and transverse to the direction of average motion) of the relativistic position probability density for a massive, spinless particle. We find that when the momentum width of…
Particle-mesh methods, such as the particle-in-cell (PIC) method, cannot retain exact pairwise interaction at sub-cell scales. For dense nonneutral relativistic electron bunches, this makes it difficult to accurately capture the…
Approximate Bayesian Computation (ABC) enables statistical inference in simulator-based models whose likelihoods are difficult to calculate but easy to simulate from. ABC constructs a kernel-type approximation to the posterior distribution…
Probabilistic principal component analysis (PCA) and its Bayesian variant (BPCA) are widely used for dimension reduction in machine learning and statistics. The main advantage of probabilistic PCA over the traditional formulation is…
Particle-based methods include a variety of techniques, such as Markov Chain Monte Carlo (MCMC) and Sequential Monte Carlo (SMC), for approximating a probabilistic target distribution with a set of weighted particles. In this paper, we…
We investigate the random walk process in relativistic flow. In the relativistic flow, photon propagation is concentrated in the directions of the flow velocity due to relativistic beaming effect. We show that, in the pure scattering case,…
Delayed-acceptance Markov chain Monte Carlo (DA-MCMC) samples from a probability distribution via a two-stages version of the Metropolis-Hastings algorithm, by combining the target distribution with a "surrogate" (i.e. an approximate and…
Starting from the model of continuous time random walk, we focus our interest on random walks in which the probability distributions of the waiting times and jumps have fat tails characterized by power laws with exponent between 0 and 1 for…
Accurate radiative transfer coefficients (emissivities, absorptivities, and rotativities) are needed for modeling radiation from relativistically hot, magnetized plasmas such as those found in Event Horizon Telescope sources. Here we…
Polarized models of relativistically hot astrophysical plasmas require transport coefficients as input: synchrotron absorption and emission coefficients in each of the four Stokes parameters, as well as three Faraday rotation coefficients.…