Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userNonequilibrium path-ensemble averages for symmetric protocols
158
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
According to the nonequilibrium work relations, path-ensembles generated by irreversible processes in which a system is driven out of equilibrium according to a predetermined protocol may be used to compute equilibrium free energy differences and expectation values. Estimation has previously been improved by considering data collected from the reverse process, which starts in equilibrium in the final thermodynamic state of the forward process and is driven according to the time-reversed protocol. Here, we develop a theoretically rigorous statistical estimator for nonequilibrium path-ensemble averages specialized for symmetric protocols, in which forward and reverse processes are identical. The estimator is tested with a number of model systems: a symmetric 1D potential, an asymmetric 1D potential, the unfolding of deca-alanine, separating a host-guest system, and translocating a potassium ion through a gramicidin A ion channel. When reconstructing free energies using data from symmetric protocols, the new estimator outperforms existing rigorous unidirectional and bidirectional estimators, converging more quickly and resulting in smaller error. However, in most cases, using the bidirectional estimator with data from a forward and reverse pair of asymmetric protocols outperforms the corresponding symmetric protocol and estimator with the same amount of simulation time. Hence, the new estimator is only recommended when the bidirectional estimator is not feasible or is expected to perform poorly. The symmetric estimator has similar performance to a unidirectional protocol of half the length and twice the number of trajectories.
rate research
Read More
The sensitivity of molecular dynamics on changes in the potential energy function plays an important role in understanding the dynamics and function of complex molecules.We present a method to obtain path ensemble averages of a perturbed dynamics from a set of paths generated by a reference dynamics. It is based on the concept of path probability measure and the Girsanov theorem, a result from stochastic analysis to estimate a change of measure of a path ensemble. Since Markov state models (MSM) of the molecular dynamics can be formulated as a combined phase-space and path ensemble average, the method can be extended toreweight MSMs by combining it with a reweighting of the Boltzmann distribution. We demonstrate how to efficiently implement the Girsanov reweighting in a molecular dynamics simulation program by calculating parts of the reweighting factor on the fly during the simulation, and we benchmark the method on test systems ranging from a two-dimensional diffusion process to an artificial many-body system and alanine dipeptide and valine dipeptide in implicit and explicit water. The method can be used to study the sensitivity of molecular dynamics on external perturbations as well as to reweight trajectories generated by enhanced sampling schemes to the original dynamics.
We study stochastic processes in which the trajectories are constrained so that the process realises a large deviation of the unconstrained process. In particular we consider stochastic bridges and the question of inequivalence of path ensembles between the microcanonical ensemble, in which the end points of the trajectory are constrained, and the canonical or s ensemble in which a bias or tilt is introduced into the process. We show how ensemble inequivalence can be manifested by the phenomenon of temporal condensation in which the large deviation is realised in a vanishing fraction of the duration (for long durations). For diffusion processes we find that condensation happens whenever the process is subject to a confining potential, such as for the Ornstein-Uhlenbeck process, but not in the borderline case of dry friction in which there is partial ensemble equivalence. We also discuss continuous-space, discrete-time random walks for which in the case of a heavy tailed step-size distribution it is known that the large deviation may be achieved in a single step of the walk. Finally we consider possible effects of several constraints on the process and in particular give an alternative explanation of the interaction-driven condensation in terms of constrained Brownian excursions.
We present generalized-ensemble algorithms for isobaric-isothermal molecular simulations. In addition to the multibaric-multithermal algorithm and replica-exchange method for the isobaric-isothermal ensemble, which have already been proposed, we propose a simulated tempering method for this ensemble. We performed molecular dynamics simulations with these algorithms for an alanine dipeptide system in explicit water molecules to test the effectiveness of the algorithms. We found that these generalized-ensemble algorithms are all useful for conformational sampling of biomolecular systems in the isobaric-isothermal ensemble.
The emph{ab initio} path integral Monte Carlo (PIMC) approach is one of the most successful methods in quantum many-body theory. A particular strength of this method is its straightforward access to imaginary-time correlation functions (ITCF). For example, the well-known density-density ITCF $F(mathbf{q},tau)$ allows one to estimate the linear response of a given system for all wave vectors $mathbf{q}$ from a single simulation of the unperturbed system. Moreover, it constitutes the basis for the reconstruction of the dynamic structure factor $S(mathbf{q},omega)$ -- a key quantity in state-of-the-art scattering experiments. In this work, we present analogous relations between the nonlinear density response in quadratic and cubic order of the perturbation strength and generalized ITCFs measuring correlations between up to four imaginary-time arguments. As a practical demonstration of our new approach, we carry out simulations of the warm dense electron gas and find excellent agreement with previous PIMC results that had been obtained with substantially larger computational effort. In addition, we give a relation between a cubic ITCF and the triple dynamic structure factor $S(mathbf{q}_1,omega_1;mathbf{q}_2,omega_2)$, which evokes the enticing possibility to study dynamic three-body effects on an emph{ab initio} level.
Recent work by Teifel and Mahler [Eur. Phys. J. B 75, 275 (2010)] raises legitimate concerns regarding the validity of quantum nonequilibrium work relations in processes involving moving hard walls. We study this issue in the context of the rapidly expanding one-dimensional quantum piston. Utilizing exact solutions of the time-dependent Schru007fodinger equation, we find that the evolution of the wave function can be decomposed into static and dynamic components, which have simple semiclassical interpretations in terms of particle-piston collisions. We show that nonequilibrium work relations remains valid at any finite piston speed, provided both components are included, and we study explicitly the work distribution for this model system.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
