Related papers: A Formulation of the Ring Polymer Molecular Dynami…
Spatial multiscale methods have established themselves as useful tools for extending the length scales accessible by conventional statics (i.e., zero temperature molecular dynamics). Recently, extensions of these methods, such as the…
It was recently shown that path integral Monte Carlo can be used to directly compute partition functions of Hamiltonians with vibronic coupling [J. Chem. Phys. 148, 194110 (2018)]. While the importance sampling Monte Carlo integration…
This work combines for the first time ab initio molecular dynamics (AIMD) within the Born-Oppenheimer approximation, with a global natural orbital functional (GNOF), an approximate functional of the one-particle reduced density matrix. The…
This paper presents a multi-temporal formulation for simulating elastoplastic solids under cyclic loading. We leverage the proper generalized decomposition (PGD) to decompose the displacements into multiple time scales, separating the…
Replica exchange molecular dynamics (REMD) becomes more efficient as the frequency of swap between the temperatures is increased. Recently in [Plattner et al, J. Chem. Phys. 135, 134111 (2011)] a method was proposed to implement infinite…
A new approach for integration of motion in many-body systems of interacting polyatomic molecules is proposed. It is based on splitting time propagation of pseudo-variables in a modified phase space, while the real translational and…
Peridynamics (PD) is a non-local continuum formulation. The original version of PD was restricted to bond-based interactions. Bond-based PD is geometrically exact and its kinematics are similar to classical continuum mechanics (CCM).…
Methods to efficiently determine the relative stability of polymorphs of organic crystals are highly desired in crystal structure predictions (CSPs). Current methodologies include use of static lattice phonons, quasi-harmonic approximation…
Path Integral Molecular Dynamics (PIMD) is a well established simulation technique to compute exact equilibrium properties for a quantum system using classical trajectories in an extended phase space. Standard PIMD simulations are…
We present an analysis of the Feynman path centroid density that provides new insight into the correspondence between the path integral and the Schr\"odinger formulations of statistical mechanics. The path centroid density is a central…
Matsubara dynamics has recently emerged as the most general form of a quantum-Boltzmann-conserving classical dynamics theory for the calculation of single-time correlation functions. Here, we present a generalization of Matsubara dynamics…
We present a new formulation of ab initio molecular dynamics which exploits the efficiency of plane waves in adaptive curvilinear coordinates, and thus provides an accurate treatment of first-row elements. The method is used to perform a…
Path integral quantum Monte Carlo (PIMC) is a method for estimating thermal equilibrium properties of stoquastic quantum spin systems by sampling from a classical Gibbs distribution using Markov chain Monte Carlo. The PIMC method has been…
We investigate the anomalous dynamics of unentangled polymer melts. The proposed equation of motion formally relates the anomalous center-of-mass diffusion, as observed in computer simulations and experiments, to the nature of the effective…
In this paper we explore the use of block coordinate descent (BCD) to optimize the centroidal momentum dynamics for dynamically consistent multi-contact behaviors. The centroidal dynamics have recently received a large amount of attention…
A Molecular Dynamics (MD) study of static and dynamic properties of molten and glassy germanium dioxide is presented. The interactions between the atoms are modelled by the classical pair potential proposed by Oeffner and Elliott (OE)…
We have implemented the Centroid Molecular Dynamics scheme (CMD) into the Grand Canonical-like version of the Adaptive Resolution Simulation Molecular Dynamics (GC-AdResS) method. We have tested the implementation on two different systems,…
An approach is described to calculate Generalized Parton Distributions (GPDs) in Constituent Quark Models (CQM). The GPDs are obtained from wave functions to be evaluated in a given CQM. The general relations linking the twist-two GPDs to…
ProtoMD is a toolkit that facilitates the development of algorithms for multiscale molecular dynamics (MD) simulations. It is designed for multiscale methods which capture the dynamic transfer of information across multiple spatial scales,…
The Matrix Product method (MPM) has been used in the past to generate variational ansatzs of the ground state (GS) of spin chains and ladders. In this paper we apply the MPM to study the GS of conjugated polymers in the valence bond basis,…