Related papers: Path integral Monte Carlo simulation of helium at …
We describe a Monte Carlo procedure which allows sampling of the disjoint configuration spaces associated with crystalline and fluid phases, within a single simulation. The method utilises biased sampling techniques to enhance the…
Thermal properties of graphene monolayers are studied by path-integral molecular dynamics (PIMD) simulations, which take into account the quantization of vibrational modes in the crystalline membrane, and allow one to consider anharmonic…
Variational methods are used to calculate structural and thermodynamical properties of a titrating polyelectrolyte in a discrete representation. The Coulomb interactions are emulated by harmonic repulsive forces, the force constants being…
We introduce a general Monte Carlo scheme for achieving atomistic simulations with monoelectronic Hamiltonians including the thermalization of both nuclear and electronic degrees of freedom. The kinetic Monte Carlo algorithm is used to…
A method for Monte Carlo simulation of 2D spin-polarized electron transport in III-V semiconductor heterojunction FETs is presented. In the simulation, the dynamics of the electrons in coordinate and momentum space is treated…
Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum…
We introduce two novel quantum Monte Carlo methods and employ them to study the superfluid-insulator transition in a two-dimensional system of hard-core bosons. One of the methods is appropriate for zero temperature and is based upon…
This work describes a new 1D hybrid approach for modeling atmospheric pressure discharges featuring complex chemistry. In this approach electrons are described fully kinetically using Particle-In-Cell/Monte-Carlo (PIC/MCC) scheme, whereas…
We develop a series of molecular dynamics computer simulations of liquid water, performed with a polarizable potential model, to calculate the spinodal line and the curve of maximum density inside the metastable supercooled region. After…
Although histogram methods have been extremely effective for analyzing data from Monte Carlo simulations, they do have certain limitations, including the range over which they are valid and the difficulties of combining data from…
In plasma edge simulations, kinetic Monte Carlo (MC) is often used to simulate neutral particles and estimate source terms. For large-sized reactors, like ITER and DEMO, high particle collision rates lead to a substantial computational cost…
The coexistence between the hcp and bcc phases of solid 4He at fixed pressure, and at three different temperatures, is studied by means of the path integral Monte Carlo method. Microscopic results for the main energetic and structure…
Thermodynamic properties of chiral spin liquids are investigated for a variant of the Kitaev model defined on a decorated honeycomb lattice. Using the quantum Monte Carlo simulation, we find that the model exhibits a finite-temperature…
We calculate the hydrogen Hugoniot using ab initio path integral Monte Carlo. We introduce an efficient finite-temperature fixed-node approximation for handling fermions, which includes an optimized mixture of free particle states and…
A new approach to simulating rotational cooling using a direct simulation Monte Carlo (DSMC) method is described and applied to the rotational cooling of ammonia seeded into a helium supersonic jet. The method makes use of ab initio…
A two-dimensional lattice hard-core boson system with a small fraction of bosonic or fermionic impurity particles is studied. The impurities have the same hopping and interactions as the dominant bosons and their effects are solely due to…
We present results of a Monte Carlo study of temperature-programmed desorption in a model system with attractive lateral interactions. It is shown that even for weak interactions there are large shifts of the peak maximum temperatures with…
Using molecular dynamics simulations, we probe a suspended graphene layer by a diamond-like-carbon tip at various temperatures. The force acting on the tip in the sliding direction is measured to be negative at liquid-helium temperature.…
We consider a Mott transition of the Hubbard model in infinite dimensions. The dynamical mean- field theory is employed in combination with a continuous-time quantum Monte Carlo (CTQMC) method for an accurate description at low…
A real-time path integral Monte Carlo approach is developed to study the dynamics in a many-body quantum system until reaching a nonequilibrium stationary state. The approach is based on augmenting an exact reduced equation for the…