Constant pressure hybrid Monte Carlo simulations in GROMACS
Mario Fernández-Pendás1, Bruno Escribano, Tijana Radivojević
1Basque Center for Applied Mathematics, E-48009, Bilbao, Spain.
Abstract:
Adaptation and implementation of the Generalized Shadow Hybrid Monte Carlo (GSHMC) method for molecular simulation at constant pressure in the NPT ensemble are discussed. The resulting method, termed NPT-GSHMC, combines Andersen barostat with GSHMC to enable molecular simulations in the environment natural for biological applications, namely, at constant pressure and constant temperature. Generalized Hybrid Monte Carlo methods are designed to maintain constant temperature and volume and extending their functionality to preserving pressure is not trivial. The theoretical formulation of NPT-GSHMC was previously introduced. Our main contribution is the implementation of this methodology in the GROMACS molecular simulation package and the evaluation of properties of NPT-GSHMC, such as accuracy, performance, effectiveness for real physical systems in comparison with well-established molecular simulation techniques. Benchmarking tests are presented and the obtained preliminary results are promising. For the first time, the generalized hybrid Monte Carlo simulations at constant pressure are available within the popular open source molecular dynamics software package.
More Related Videos
10:52Multiscale Sampling of a Heterogeneous Water/Metal Catalyst Interface using Density Functional Theory and Force-Field Molecular Dynamics
Published on: April 12, 2019
05:56Exploring Caspase Mutations and Post-Translational Modification by Molecular Modeling Approaches
Published on: October 13, 2022
Related Concept Videos
Constant Pressure Calorimetry
Compartment Models: Single-Compartment Model
Pressure Variation in a Fluid at Rest
When measuring pressure at two different levels within the fluid, the difference in...
Hydrostatic Pressure Force on a Curved Surface
Pressure of Fluids
Homogeneous Equilibria for Gaseous Reactions
For gas-phase reactions, the equilibrium constant may be expressed in terms of either the molar concentrations (Kc) or partial pressures (Kp) of the reactants and products. A relation between these two K values may be simply derived from the ideal gas equation and the definition of molarity. According to the ideal gas equation:
