Related papers: Self-learning hybrid Monte Carlo method for isothe…
Metropolis Monte Carlo simulation is a powerful tool for studying the equilibrium properties of matter. In complex condensed-phase systems, however, it is difficult to design Monte Carlo moves with high acceptance probabilities that also…
Hamiltonian Monte Carlo is a prominent Markov Chain Monte Carlo algorithm, which employs symplectic integrators to sample from high dimensional target distributions in many applications, such as statistical mechanics, Bayesian statistics…
The self-healing diffusion Monte Carlo algorithm (SHDMC) [Phys. Rev. B {\bf 79}, 195117 (2009), {\it ibid.} {\bf 80}, 125110 (2009)] is shown to be an accurate and robust method for calculating the ground state of atoms and molecules. By…
We report the development of a pattern-recognition scheme that takes into account both fcc and hcp adsorption sites in performing self-learning kinetic Monte Carlo (SLKMC-II) simulations on the fcc(111) surface. In this scheme, the local…
The recently-introduced self-learning Monte Carlo method is a general-purpose numerical method that speeds up Monte Carlo simulations by training an effective model to propose uncorrelated configurations in the Markov chain. We implement…
The Kinetic Monte Carlo (KMC) method has become an important tool for examination of phenomena like surface diffusion and thin film growth because of its ability to carry out simulations for time scales that are relevant to experiments. But…
The shell-model Monte Carlo (SMMC) technique transforms the traditional nuclear shell-model problem into a path-integral over auxiliary fields. We describe below the method and its applications to four physics issues: calculations of sdpf-…
We consider parallel asynchronous Markov Chain Monte Carlo (MCMC) sampling for problems where we can leverage (stochastic) gradients to define continuous dynamics which explore the target distribution. We outline a solution strategy for…
Complex soft matter systems can be efficiently studied with the help of adaptive resolution simulation methods, concurrently employing two levels of resolution in different regions of the simulation domain. The non-matching properties of…
We show how accurate benchmark values of the surface formation energy of crystalline lithium hydride can be computed by the complementary techniques of quantum Monte Carlo (QMC) and wavefunction-based molecular quantum chemistry. To…
We describe a new Hybrid Monte Carlo (HMC) algorithm for dynamical overlap fermions, which improves the rate of topological index changes by adding an additional (intensive) term to the action for the molecular dynamics part of the…
Lattice-switch Monte Carlo (LSMC) is a method for evaluating the free energy between two given solid phases. LSMC is a general method, being applicable to a wide range of problems and interatomic potentials. Furthermore it is extremely…
Predicting the thermodynamic properties of mixtures is crucial for process design and optimization in chemical engineering. Machine learning (ML) methods are gaining increasing attention in this field, but experimental data for training are…
Warm dense matter (WDM) is an active field of research, with applications ranging from astrophysics to inertial confinement fusion. Ionization degree and continuum lowering are important quantities to understand how materials behave under…
In the last few years we have been developing a Monte Carlo simulation method to cope with systems of many electrons and ions in the Born-Oppenheimer (BO) approximation, the Coupled Electron-Ion Monte Carlo Method (CEIMC). Electronic…
The hybrid Monte Carlo (HMC) algorithm is applied for the Bayesian inference of the stochastic volatility (SV) model. We use the HMC algorithm for the Markov chain Monte Carlo updates of volatility variables of the SV model. First we…
We present a novel way of performing kinetic Monte Carlo simulations which does not require an {\it a priori} list of diffusion processes and their associated energetics and reaction rates. Rather, at any time during the simulation,…
The hybrid Monte Carlo algorithm (HMCA) is applied for Bayesian parameter estimation of the realized stochastic volatility (RSV) model. Using the 2nd order minimum norm integrator (2MNI) for the molecular dynamics (MD) simulation in the…
In online clustering problems, there is often a large amount of uncertainty over possible cluster assignments that cannot be resolved until more data are observed. This difficulty is compounded when clusters follow complex distributions, as…
An overview of recent work on Monte Carlo simulations of a granular binary mixture is presented. The results are obtained numerically solving the Enskog equation for inelastic hard-spheres by means of an extension of the well-known direct…