Related papers: Ab-initio molecular dynamics for high-pressure liq…
We develop a lattice-based Monte Carlo simulation method for charged mixtures capable of treating dielectric heterogeneities. Using this method, we study oil-water mixtures containing an antagonistic salt, with hydrophilic cations and…
We report an ab-initio simulation of liquid Indium in a wide range of pressures and temperatures. We calculate equation of state, thermal expansion and compressibility coefficients. The structure of the system is analyzed by radial…
We report first principle results for nuclear structure and optical responses of high pressure liquid hydrogen along two isotherms in the region of molecular dissociation. We employ Density Functional Theory with the vdW-DF approximation…
In order to find the equilibrium geometries of molecules and solids and to perform ab initio molecular dynamics, it is necessary to calculate the forces on the nuclei. We present a correlated sampling method to efficiently calculate…
Ab initio quantum Monte Carlo (QMC) is a state-of-the-art numerical approach for evaluating accurate expectation values of many-body wavefunctions. However, one of the major drawbacks that still hinders widespread QMC applications is the…
We have developed a Monte Carlo simulation for ion transport in hot background gases, which is an alternative way of solving the corresponding Boltzmann equation that determines the distribution function of ions. We consider the limit of…
Easy physics-inspired approximations of the total and binding energies for the ${\rm H}$ atom and for the molecular ions $${\rm H}_2^{(+)} ({\rm ppe}), {\rm H}_3^{(2+)} ({\rm pppe}), ({\rm HeH})^{++} (\al {\rm p e}), {\rm He}_2^{(3+)} (\al…
An efficient Path Integral Monte Carlo procedure is proposed to simulate the behavior of quantum many-body dissipative systems described within the framework of the influence functional. Thermodynamic observables are obtained by Monte Carlo…
This paper presents an $ab$ $initio$ study of hydrogen dynamics inside nanometer-size voids in $a$-Si within the framework of the density-functional theory for a varying hydrogen load of 10 to 30 H atoms/void at the low and high temperature…
Algorithms to determine transition probabilities in Monte Carlo simulations are tested using a system of classical particles with effective interactions which reproduce Bose-Einstein statistics. The system is appropriate for testing…
In this paper, we consider mathematical modeling and numerical simulation of non-isothermal compressible multi-component diffuse-interface two-phase flows with realistic equations of state. A general model with general reference velocity is…
Proton ordering in water ice is a paradigmatic order-disorder transition in a locally constrained system. The ice rules require exactly two hydrogens close to each oxygen, restricting the disorder to an exponentially large yet strongly…
Molecular crystals play a central role in a wide range of scientific fields, including pharmaceuticals and organic semiconductor devices. However, they are challenging systems to model accurately with computational approaches because of a…
We study a two-dimensional fluid of dipolar hard disks by Monte Carlo simulations in a square with periodic boundary conditions and on the surface of a sphere. The theory of the dielectric constant and the asymptotic behaviour of the…
Monte Carlo methods are used to study the phase transition in ammonium chloride from the orientationally ordered $\delta$ phase to the orientationally disordered $\gamma$ phase. An effective pair potential is used to model the interaction…
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…
Two of the most challenging tasks in molecular simulation consist in capturing the properties of systems with long-range interactions (e.g. electrolyte solutions) as well as systems containing large molecules such as hydrogels. For the…
An ab initio quantum Monte Carlo method is introduced for calculating total rates of muon weak capture in light nuclei with mass number $A \leq 12$. As a first application of the method, we perform a calculation of the rate in $^4$He in a…
A method for computing the thermopower in interacting systems is proposed. This approach, which relies on Monte Carlo simulations, is illustrated first for a diatomic chain of hard-point elastically colliding particles and then in the case…
We examine the fate of fast electrons (with energies E>10 eV) in a thermal gas of primordial composition. To follow their interactions with the background gas, we construct a Monte Carlo model that includes: (1) electron-electron scattering…