Constant pressure hybrid Monte Carlo simulations in GROMACS

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...

ver descrição completa

Detalhes bibliográficos
Autores: Fernández-Pendás, M., Escribano, B., Radivojevic, T., Akhmatskaya, E.
Formato: artículo
Estado:Versión aceptada para publicación
Fecha de publicación:2014
País:España
Recursos:Basque Center for Applied Mathematics (BCAM)
Repositorio:BIRD. BCAM's Institutional Repository Data
OAI Identifier:oai:bird.bcamath.org:20.500.11824/167
Acesso em linha:http://hdl.handle.net/20.500.11824/167
Access Level:acceso abierto
Palavra-chave:accuracy
Article
barostat
constant pressure
generalized shadow hybrid monte carl simulation
molecular dynamics
Monte Carlo method
pressure
temperature
controlled study
Generalized Shadow Hybrid Monte Carlo
measurement accuracy
methodology
oscillation
priority journal
quality control
simulation
theoretical study
villin
Descrição
Resumo: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.