Related papers: Multi-level Monte Carlo path integral molecular dy…
Correlated fermions are of high interest in condensed matter (Fermi liquids, Wigner molecules), cold atomic gases and dense plasmas. Here we propose a novel approach to path integral Monte Carlo (PIMC) simulations of strongly degenerate…
We present a new approach to path integral Monte Carlo (PIMC) simulations based on the worm algorithm, originally developed for lattice models and extended here to continuous-space many-body systems. The scheme allows for efficient…
Path integral Monte Carlo with Green's function analysis allows the sampling of quantum mechanical properties of molecules at finite temperature. While a high-precision computation of the energy of the Born-Oppenheimer surface from path…
In a typical finite temperature quantum Monte Carlo (QMC) simulation, estimators for simple static observables such as specific heat and magnetization are known. With a great deal of system-specific manual labor, one can sometimes also…
We present a method to compute real-time path integrals numerically, by Monte-Carlo sampling on near-Lefschetz thimbles. We present a collection of tools based on the Lefschetz thimble methods, which together provide an alternative to…
We present a path-integral Monte Carlo estimator for calculating the dipole polarizability of interacting Coulomb plasma in the long-wavelength limit, i.e., the optical region. We present comprehensive details and method validation studies…
We report a novel multi-scale simulation methodology to quantitatively predict the thermodynamic behaviour of polymer mixtures, that exhibit phases with broken orientational symmetry. Our system consists of a binary mixture of oligomers and…
In this work, we introduce a simple modification of the Monte Carlo algorithm, which we call step Monte Carlo (sMC). The sMC approach allows to simulate processes far from equilibrium and obtain information about the dynamic properties of…
The nested sampling algorithm has been shown to be a general method for calculating the pressure-temperature-composition phase diagrams of materials. While the previous implementation used single-particle Monte Carlo moves, these are…
Being motivated by the surge of fermionic quantum Monte Carlo simulations at finite temperature, we present a detailed analysis of the permutation-cycle properties of path integral Monte Carlo (PIMC) simulations of degenerate electrons.…
Precise knowledge of the static density response function (SDRF) of the uniform electron gas (UEG) serves as key input for numerous applications, most importantly for density functional theory beyond generalized gradient approximations.…
Convergence with respect to imaginary-time discretization is an essential part of any path-integral-based calculation. However, an unfortunate property of existing non-preconditioned numerical integration schemes for path-integral molecular…
In a recent publication [S. Groth \textit{et al.}, PRB (2016)], we have shown that the combination of two novel complementary quantum Monte Carlo approaches, namely configuration path integral Monte Carlo (CPIMC) [T. Schoof \textit{et al.},…
In this article, we study the application of Multi-Level Monte Carlo (MLMC) approaches to numerical random homogenization. Our objective is to compute the expectation of some functionals of the homogenized coefficients, or of the…
We extend the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), a semiclassical method for computing real-time correlation functions, to electronically nonadiabatic systems using the Meyer-Miller-Stock-Thoss (MMST) Hamiltonian…
The advances in materials and biological sciences have necessitated the use of molecular simulations to study polymers. The Markov chain Monte Carlo simulations enable the sampling of relevant microstates of polymeric systems by traversing…
In this article we consider computing expectations w.r.t.~probability laws associated to a certain class of stochastic systems. In order to achieve such a task, one must not only resort to numerical approximation of the expectation, but…
In this paper we present a Path Integral Monte Carlo (PIMC) simulation of the orthorhombic phase of crystalline polyethylene, using an explicit atom force field with unconstrained bond lengths and angles. This work represents a quantum…
The basic problem in equilibrium statistical mechanics is to compute phase space average, in which Monte Carlo method plays a very important role. We begin with a review of nonlocal algorithms for Markov chain Monte Carlo simulation in…
We present a novel hybrid computational method to simulate accurately dendritic solidification in the low undercooling limit where the dendrite tip radius is one or more orders of magnitude smaller than the characteristic spatial scale of…