Related papers: Hard-sphere melting and crystallization with event…
Microscopic processes on surfaces such as adsorption, desorption, diffusion and reaction of interacting particles can be simulated using kinetic Monte Carlo (kMC) algorithms. Even though kMC methods are accurate, they are computationally…
We carry out Monte Carlo simulations to investigate the effect of molecular shape on liquid-crystal order. In our approach, each model mesogen consists of several soft spheres bonded rigidly together. The arrangement of the spheres may be…
A simple and efficient algorithm of the molecular-dynamics simulation of the hard disk system based on the Event-Driven method is developed. From the analysis of algorithm, the complexity is O(log N) per 1 event, and the constant…
The hard-sphere model is one of the most extensively studied models in statistical physics. It describes the continuous distribution of spherical particles, governed by hard-core interactions. An important quantity of this model is the…
This chapter is devoted to the computation of equilibrium (thermodynamic) properties of quantum systems. In particular, we will be interested in the situation where the interaction between particles is so strong that it cannot be treated as…
Event-driven molecular dynamics is a valuable tool in condensed and soft matter physics when particles can be modeled as hard objects or more generally if their interaction potential can be modeled in a stepwise fashion. Hard spheres model…
Diffusion Monte Carlo (DMC) simulations for fermions are becoming the standard to provide high quality reference data in systems that are too large to be investigated via quantum chemical approaches. DMC with the fixed-node approximation…
The approach to the ergodic limit in Monte Carlo simulations is studied using both analytic and numerical methods. With the help of a stochastic model, a metric is defined that enables the examination of a simulation in both the ergodic and…
Self-learning hybrid Monte Carlo (SLHMC) is a first-principles simulation that allows for exact ensemble generation on potential energy surfaces based on density functional theory. The statistical sampling can be accelerated with the…
The crystal nucleation from liquid in most cases is too rare to be accessed within the limited timescales of the conventional molecular dynamics (MD) simulation. Here, we developed a "persistent embryo" method to facilitate crystal…
A kinetic Monte Carlo (KMC) method is used to study the structural properties and dynamics of a supercooled binary Lennard-Jones liquid around the glass transition temperature. This technique permits us to explore the potential energy…
We prove bounds on the variance of a function $f$ under the empirical measure of the samples obtained by the Sequential Monte Carlo (SMC) algorithm, with time complexity depending on local rather than global Markov chain mixing dynamics.…
We perform molecular dynamics (MD) simulations of the crystallization process in binary Lennard-Jones systems during heating and cooling to investigate atomic-scale crystallization kinetics in glass-forming materials. For the cooling…
We propose a Multi-Cell Monte Carlo algorithm, or (MC)^2, for predicting stable phases in chemically complex crystalline systems. Free atomic transfer among cells is achieved via the application of the lever rule, where an assigned molar…
We consider a theoretical model for a binary mixture of colloidal particles and spherical emulsion droplets. The hard sphere colloids interact via additional short-ranged attraction and long-ranged repulsion. The droplet-colloid interaction…
We use specialized Monte Carlo simulation methods and moment free energy calculations to provide conclusive evidence that dense polydisperse spheres at equilibrium demix into coexisting fcc phases, with more phases appearing as the spread…
A common way to simulate the transport and spread of pollutants in the atmosphere is via stochastic Lagrangian dispersion models. Mathematically, these models describe turbulent transport processes with stochastic differential equations…
This paper presents a seamless algorithm for the application of the multilevel Monte Carlo (MLMC) method to the ensemble transform particle filter (ETPF). The algorithm uses a combination of optimal coupling transformations between coarse…
From basic considerations of the Lie group that preserves a target probability measure, we derive the Barker, Metropolis, and ensemble Markov chain Monte Carlo (MCMC) algorithms, as well as variants of waste-recycling Metropolis-Hastings…
We examine, using molecular dynamics simulation, the structure and thermodynamics of the (100) and (111) disordered face-centered cubic (FCC) crystal/melt interfaces for a binary hard-sphere system. This study is an extension of our…