QM/MM simulation of liquid water with an adaptive quantum region

Noam Bernstein, Csilla Várnai, Iván Solt, Steven A. Winfield, Mike C. Payne, I. Simon, M. Fuxreiter, Gábor Csányi

Research output: Contribution to journalArticle

60 Citations (Scopus)

Abstract

The simulation of complex chemical systems often requires a multi-level description, in which a region of special interest is treated using a computationally expensive quantum mechanical (QM) model while its environment is described by a faster, simpler molecular mechanical (MM) model. Furthermore, studying dynamic effects in solvated systems or bio-molecules requires a variable definition of the two regions, so that atoms or molecules can be dynamically re-assigned between the QM and MM descriptions during the course of the simulation. Such reassignments pose a problem for traditional QM/MM schemes by exacerbating the errors that stem from switching the model at the boundary. Here we show that stable, long adaptive simulations can be carried out using density functional theory with the BLYP exchange-correlation functional for the QM model and a flexible TIP3P force field for the MM model without requiring adjustments of either. Using a primary benchmark system of pure water, we investigate the convergence of the liquid structure with the size of the QM region, and demonstrate that by using a sufficiently large QM region (with radius 6 ) it is possible to obtain radial and angular distributions that, in the QM region, match the results of fully quantum mechanical calculations with periodic boundary conditions, and, after a smooth transition, also agree with fully MM calculations in the MM region. The key ingredient is the accurate evaluation of forces in the QM subsystem which we achieve by including an extended buffer region in the QM calculations. We also show that our buffered-force QM/MM scheme is transferable by simulating the solvated Cl - ion.

Original languageEnglish
Pages (from-to)646-656
Number of pages11
JournalPhysical Chemistry Chemical Physics
Volume14
Issue number2
DOIs
Publication statusPublished - Jan 14 2012

Fingerprint

Water
Liquids
liquids
water
simulation
Molecules
Angular distribution
stems
radial distribution
ingredients
field theory (physics)
Density functional theory
molecules
Buffers
angular distribution
buffers
adjusting
Boundary conditions
Ions
boundary conditions

ASJC Scopus subject areas

  • Physical and Theoretical Chemistry
  • Physics and Astronomy(all)

Cite this

Bernstein, N., Várnai, C., Solt, I., Winfield, S. A., Payne, M. C., Simon, I., ... Csányi, G. (2012). QM/MM simulation of liquid water with an adaptive quantum region. Physical Chemistry Chemical Physics, 14(2), 646-656. https://doi.org/10.1039/c1cp22600b

QM/MM simulation of liquid water with an adaptive quantum region. / Bernstein, Noam; Várnai, Csilla; Solt, Iván; Winfield, Steven A.; Payne, Mike C.; Simon, I.; Fuxreiter, M.; Csányi, Gábor.

In: Physical Chemistry Chemical Physics, Vol. 14, No. 2, 14.01.2012, p. 646-656.

Research output: Contribution to journalArticle

Bernstein, N, Várnai, C, Solt, I, Winfield, SA, Payne, MC, Simon, I, Fuxreiter, M & Csányi, G 2012, 'QM/MM simulation of liquid water with an adaptive quantum region', Physical Chemistry Chemical Physics, vol. 14, no. 2, pp. 646-656. https://doi.org/10.1039/c1cp22600b
Bernstein, Noam ; Várnai, Csilla ; Solt, Iván ; Winfield, Steven A. ; Payne, Mike C. ; Simon, I. ; Fuxreiter, M. ; Csányi, Gábor. / QM/MM simulation of liquid water with an adaptive quantum region. In: Physical Chemistry Chemical Physics. 2012 ; Vol. 14, No. 2. pp. 646-656.
@article{c5549a1c5c33470baf6bcac028cae27c,
title = "QM/MM simulation of liquid water with an adaptive quantum region",
abstract = "The simulation of complex chemical systems often requires a multi-level description, in which a region of special interest is treated using a computationally expensive quantum mechanical (QM) model while its environment is described by a faster, simpler molecular mechanical (MM) model. Furthermore, studying dynamic effects in solvated systems or bio-molecules requires a variable definition of the two regions, so that atoms or molecules can be dynamically re-assigned between the QM and MM descriptions during the course of the simulation. Such reassignments pose a problem for traditional QM/MM schemes by exacerbating the errors that stem from switching the model at the boundary. Here we show that stable, long adaptive simulations can be carried out using density functional theory with the BLYP exchange-correlation functional for the QM model and a flexible TIP3P force field for the MM model without requiring adjustments of either. Using a primary benchmark system of pure water, we investigate the convergence of the liquid structure with the size of the QM region, and demonstrate that by using a sufficiently large QM region (with radius 6 ) it is possible to obtain radial and angular distributions that, in the QM region, match the results of fully quantum mechanical calculations with periodic boundary conditions, and, after a smooth transition, also agree with fully MM calculations in the MM region. The key ingredient is the accurate evaluation of forces in the QM subsystem which we achieve by including an extended buffer region in the QM calculations. We also show that our buffered-force QM/MM scheme is transferable by simulating the solvated Cl - ion.",
author = "Noam Bernstein and Csilla V{\'a}rnai and Iv{\'a}n Solt and Winfield, {Steven A.} and Payne, {Mike C.} and I. Simon and M. Fuxreiter and G{\'a}bor Cs{\'a}nyi",
year = "2012",
month = "1",
day = "14",
doi = "10.1039/c1cp22600b",
language = "English",
volume = "14",
pages = "646--656",
journal = "Physical Chemistry Chemical Physics",
issn = "1463-9076",
publisher = "Royal Society of Chemistry",
number = "2",

}

TY - JOUR

T1 - QM/MM simulation of liquid water with an adaptive quantum region

AU - Bernstein, Noam

AU - Várnai, Csilla

AU - Solt, Iván

AU - Winfield, Steven A.

AU - Payne, Mike C.

AU - Simon, I.

AU - Fuxreiter, M.

AU - Csányi, Gábor

PY - 2012/1/14

Y1 - 2012/1/14

N2 - The simulation of complex chemical systems often requires a multi-level description, in which a region of special interest is treated using a computationally expensive quantum mechanical (QM) model while its environment is described by a faster, simpler molecular mechanical (MM) model. Furthermore, studying dynamic effects in solvated systems or bio-molecules requires a variable definition of the two regions, so that atoms or molecules can be dynamically re-assigned between the QM and MM descriptions during the course of the simulation. Such reassignments pose a problem for traditional QM/MM schemes by exacerbating the errors that stem from switching the model at the boundary. Here we show that stable, long adaptive simulations can be carried out using density functional theory with the BLYP exchange-correlation functional for the QM model and a flexible TIP3P force field for the MM model without requiring adjustments of either. Using a primary benchmark system of pure water, we investigate the convergence of the liquid structure with the size of the QM region, and demonstrate that by using a sufficiently large QM region (with radius 6 ) it is possible to obtain radial and angular distributions that, in the QM region, match the results of fully quantum mechanical calculations with periodic boundary conditions, and, after a smooth transition, also agree with fully MM calculations in the MM region. The key ingredient is the accurate evaluation of forces in the QM subsystem which we achieve by including an extended buffer region in the QM calculations. We also show that our buffered-force QM/MM scheme is transferable by simulating the solvated Cl - ion.

AB - The simulation of complex chemical systems often requires a multi-level description, in which a region of special interest is treated using a computationally expensive quantum mechanical (QM) model while its environment is described by a faster, simpler molecular mechanical (MM) model. Furthermore, studying dynamic effects in solvated systems or bio-molecules requires a variable definition of the two regions, so that atoms or molecules can be dynamically re-assigned between the QM and MM descriptions during the course of the simulation. Such reassignments pose a problem for traditional QM/MM schemes by exacerbating the errors that stem from switching the model at the boundary. Here we show that stable, long adaptive simulations can be carried out using density functional theory with the BLYP exchange-correlation functional for the QM model and a flexible TIP3P force field for the MM model without requiring adjustments of either. Using a primary benchmark system of pure water, we investigate the convergence of the liquid structure with the size of the QM region, and demonstrate that by using a sufficiently large QM region (with radius 6 ) it is possible to obtain radial and angular distributions that, in the QM region, match the results of fully quantum mechanical calculations with periodic boundary conditions, and, after a smooth transition, also agree with fully MM calculations in the MM region. The key ingredient is the accurate evaluation of forces in the QM subsystem which we achieve by including an extended buffer region in the QM calculations. We also show that our buffered-force QM/MM scheme is transferable by simulating the solvated Cl - ion.

UR - http://www.scopus.com/inward/record.url?scp=83455203339&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=83455203339&partnerID=8YFLogxK

U2 - 10.1039/c1cp22600b

DO - 10.1039/c1cp22600b

M3 - Article

C2 - 22089416

AN - SCOPUS:83455203339

VL - 14

SP - 646

EP - 656

JO - Physical Chemistry Chemical Physics

JF - Physical Chemistry Chemical Physics

SN - 1463-9076

IS - 2

ER -