Related papers: Hard-sphere melting and crystallization with event…
A Monte Carlo method to sample the classical configurational canonical ensemble is introduced. In contrast to the Metropolis algorithm, where trial moves can be rejected, in this approach collisions take place. The implementation is…
We introduce a new method to simulate the physics of rare events. The method, an extension of the Temperature Accelerated Molecular Dynamics, comes in use when the collective variables introduced to characterize the rare events are either…
We review the application of lattice QCD techniques, most notably the Hybrid Monte-Carlo (HMC) simulations, to first-principle study of tight-binding models of crystalline solids with strong inter-electron interactions. After providing a…
Thermalization of heavy quarks in the quark-gluon plasma (QGP) is one of the most promising phenomena for understanding the strong interaction. The energy loss and momentum broadening at low momentum can be well described by a stochastic…
Fixed-node diffusion Monte Carlo (FNDMC) is a stochastic quantum many-body method that has a great potential in electronic structure theory. We examine how FNDMC satisfies exact constraints, linearity and derivative discontinuity of total…
A kinetic Monte Carlo (KMC) method is presented to simulate the electrodeposition of a metal on a single crystal surface of the same metal under galvanostatic conditions. This method utilizes the multi-body embedded-atom method (EAM)…
We investigate Monte Carlo simulation strategies for determining the effective ("depletion") potential between a pair of hard spheres immersed in a dense sea of much smaller hard spheres. Two routes to the depletion potential are…
Monte Carlo (MC) simulations are powerful computational tools for investigating thermodynamic behavior and validating analytical approaches in complex physical systems. Here we present ETHER (Efficient Tool for THermodynamics Exploration…
We construct asymptotic arguments for the relative efficiency of rejection-free Monte Carlo (MC) methods compared to the standard MC method. We find that the efficiency is proportional to $\exp{({const} \beta)}$ in the Ising, $\sqrt{\beta}$…
We present a lattice Monte Carlo algorithm based on the one originally proposed by Maggs and Rossetto for simulating electrostatic interactions in inhomogeneous dielectric media. The original algorithm is known to produce attractive…
Crystallization processes at the mesoscopic scale, where faceted, dendritic growth, and multigrain formation can be observed, are of particular interest within materials science and metallurgy. These processes are highly nonlinear,…
We propose a modification, based on the RESTART (repetitive simulation trials after reaching thresholds) and DPR (dynamics probability redistribution) rare event simulation algorithms, of the standard diffusion Monte Carlo (DMC) algorithm.…
Competing phases or interactions in complex many-particle systems can result in free energy barriers that strongly suppress thermal equilibration. Here we discuss how extended ensemble Monte Carlo simulations can be used to study the…
We use replica exchange Monte-Carlo simulations to measure the equilibrium equation of state of the disordered fluid state for a binary hard sphere mixture up to very large densities where standard Monte-Carlo simulations do not easily…
A method for making realistic estimates of the density of levels in even-even nuclei is presented making use of the Monte Carlo shell model (MCSM). The procedure follows three basic steps: (1) computation of the thermal energy with the…
This work addresses uncertainty quantification of electromagnetic devices determined by the eddy current problem. The multilevel Monte Carlo (MLMC) method is used for the treatment of uncertain parameters while the devices are discretized…
Accurately and efficiently estimating system performance under uncertainty is paramount in power system planning and operation. Monte Carlo simulation is often used for this purpose, but convergence may be slow, especially when detailed…
We investigated the nematic to smectic transition undergone by parallel hard spherocylinders in the framework provided by the residual multi-particle entropy (RMPE) formalism. The RMPE is defined as the sum of all contributions to the…
The investigation of freezing transitions of single polymers is computationally demanding, since surface effects dominate the nucleation process. In recent studies we have systematically shown that the freezing properties of flexible,…
We introduce the energy-stepping Monte Carlo (ESMC) method, a Markov chain Monte Carlo (MCMC) algorithm based on the conventional dynamical interpretation of the proposal stage but employing an energy-stepping integrator. The…