Related papers: Markovian Monte Carlo program EvolFMC v.2 for solv…
We develop an Evolutionary Markov Chain Monte Carlo (EMCMC) algorithm for sampling spatial partitions that lie within a large and complex spatial state space. Our algorithm combines the advantages of evolutionary algorithms (EAs) as…
In most sampling algorithms, including Hamiltonian Monte Carlo, transition rates between states correspond to the probability of making a transition in a single time step, and are constrained to be less than or equal to 1. We derive a…
We present a new parallel supercomputer implementation of the Monte-Carlo method for simulating the dynamical evolution of globular star clusters. Our method is based on a modified version of Henon's Monte-Carlo algorithm for solving the…
We propose a multilevel Markov chain Monte Carlo (MCMC) method for the Bayesian inference of random field parameters in PDEs using high-resolution data. Compared to existing multilevel MCMC methods, we additionally consider level-dependent…
Irreversible and rejection-free Monte Carlo methods, recently developed in Physics under the name Event-Chain and known in Statistics as Piecewise Deterministic Monte Carlo (PDMC), have proven to produce clear acceleration over standard…
Uncertainty quantification (UQ) includes the characterization, integration, and propagation of uncertainties that result from stochastic variations and a lack of knowledge or data in the natural world. Monte Carlo (MC) method is a…
We revisit the evolution of generalised parton distributions (GPDs) in momentum space. We formulate the evolution kernels at one-loop in perturbative QCD (pQCD) in a form suitable for numerical implementation and that allows for an accurate…
The solution of DGLAP evolution equation for the twist-3 gluon operators is obtained in the Double Logarithmic Approximation of QCD perturbation theory. The method used for the solution is similar to the reggeon field theory. The…
A new algorithm for implementing the adaptive Monte Carlo method is given. It is used to solve the relativistic Boltzmann equations that describe the time evolution of a nonequilibrium electron-positron pair plasma containing high-energy…
We propose a Markov chain Monte Carlo-based deconvolution method designed to estimate the number of peaks in spectral data, along with the optimal parameters of each radial basis function. Assuming cases where the number of peaks is…
Stochastic PDE eigenvalue problems are useful models for quantifying the uncertainty in several applications from the physical sciences and engineering, e.g., structural vibration analysis, the criticality of a nuclear reactor or photonic…
Markov Chain Monte Carlo algorithms, the method of choice to sample from generic high-dimensional distributions, are rarely used for continuous one-dimensional distributions, for which more effective approaches are usually available (e.g.…
Langevin algorithms are popular Markov chain Monte Carlo (MCMC) methods for large-scale sampling problems that often arise in data science. We propose Monte Carlo algorithms based on the discretizations of $P$-th order Langevin dynamics for…
We investigate two kinds of extensions for the variational Monte Carlo (VMC) method with the Pfaffian in the nuclear shell-model calculations. One is the extension to odd-mass nuclei, for which we find a new Pfaffian expression of the VMC…
Wavefunction correction scheme, which was developed as a variance reduction tool for the pure and fixed-node diffusion Monte Carlo (DMC) computations by Anderson and Freihaut, is applied to the DMC computations of fermions without using the…
We describe a simple algorithm that calculates the distributions of electrons and positrons under multiphoton beamstrahlung at a future linear collider. The evolution equation as given by Chen is solved by a Monte Carlo algorithm. Explicit…
UKQCD's dynamical fermion project uses the Generalised Hybrid Monte-Carlo (GHMC) algorithm to generate QCD gauge configurations for a non-perturbatively O(a) improved Wilson action with two degenerate sea-quark flavours. We describe our…
We present a Monte Carlo algorithm that allows the simultaneous determination of a few extremal eigenpairs of a very large matrix without the need to compute the inner product of two vectors or store all the components of any one vector.…
We analytically solved the QED $\otimes$ QCD coupled DGLAP evolution equations at leading order (LO) quantum electrodynamics (QED) and next to leading order (NLO) quantum chromodynamics (QCD) approximations, using the Laplace transform…
It is shown that the exact dynamics of a composite quantum system can be represented through a pair of product states which evolve according to a Markovian random jump process. This representation is used to design a general Monte Carlo…