Toward the correction of effective electrostatic forces in explicit-solvent molecular dynamics simulations: restraints on solvent-generated electrostatic potential and solvent polarization

Verfasser / Beitragende:
[Maria Reif, Chris Oostenbrink]
Ort, Verlag, Jahr:
2015
Enthalten in:
Theoretical Chemistry Accounts, 134/2(2015-02-01), 1-19
Format:
Artikel (online)
ID: 605488355
LEADER caa a22 4500
001 605488355
003 CHVBK
005 20210128100455.0
007 cr unu---uuuuu
008 210128e20150201xx s 000 0 eng
024 7 0 |a 10.1007/s00214-014-1600-8  |2 doi 
035 |a (NATIONALLICENCE)springer-10.1007/s00214-014-1600-8 
245 0 0 |a Toward the correction of effective electrostatic forces in explicit-solvent molecular dynamics simulations: restraints on solvent-generated electrostatic potential and solvent polarization  |h [Elektronische Daten]  |c [Maria Reif, Chris Oostenbrink] 
520 3 |a Despite considerable advances in computing power, atomistic simulations under nonperiodic boundary conditions, with Coulombic electrostatic interactions and in systems large enough to reduce finite-size associated errors in thermodynamic quantities to within the thermal energy, are still not affordable. As a result, periodic boundary conditions, systems of microscopic size and effective electrostatic interaction functions are frequently resorted to. Ensuing artifacts in thermodynamic quantities are nowadays routinely corrected a posteriori, but the underlying configurational sampling still descends from spurious forces. The present study addresses this problem through the introduction of on-the-fly corrections to the physical forces during an atomistic molecular dynamics simulation. Two different approaches are suggested, where the force corrections are derived from special potential energy terms. In the first approach, the solvent-generated electrostatic potential sampled at a given atom site is restrained to a target value involving corrections for electrostatic artifacts. In the second approach, the long-range regime of the solvent polarization around a given atom site is restrained to the Born polarization, i.e., the solvent polarization corresponding to the ideal situation of a macroscopic system under nonperiodic boundary conditions and governed by Coulombic electrostatic interactions. The restraints are applied to the explicit-water simulation of a hydrated sodium ion, and the effect of the restraints on the structural and energetic properties of the solvent is illustrated. Furthermore, by means of the calculation of the charging free energy of a hydrated sodium ion, it is shown how the electrostatic potential restraint translates into the on-the-fly consideration of the corresponding free-energy correction terms. It is discussed how the restraints can be generalized to situations involving several solute particles. Although the present study considers a very simple system only, it is an important step toward the on-the-fly elimination of finite-size and approximate-electrostatic artifacts during atomistic molecular dynamics simulations. 
540 |a The Author(s), 2015 
690 7 |a Computer simulation  |2 nationallicence 
690 7 |a Molecular dynamics  |2 nationallicence 
690 7 |a Electrostatic artifacts  |2 nationallicence 
690 7 |a Ion solvation  |2 nationallicence 
690 7 |a Solvent polarization  |2 nationallicence 
700 1 |a Reif  |D Maria  |u Institute for Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Muthgasse 18, 1190, Vienna, Austria  |4 aut 
700 1 |a Oostenbrink  |D Chris  |u Institute for Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Muthgasse 18, 1190, Vienna, Austria  |4 aut 
773 0 |t Theoretical Chemistry Accounts  |d Springer Berlin Heidelberg  |g 134/2(2015-02-01), 1-19  |x 1432-881X  |q 134:2<1  |1 2015  |2 134  |o 214 
856 4 0 |u https://doi.org/10.1007/s00214-014-1600-8  |q text/html  |z Onlinezugriff via DOI 
898 |a BK010053  |b XK010053  |c XK010000 
900 7 |a Metadata rights reserved  |b Springer special CC-BY-NC licence  |2 nationallicence 
908 |D 1  |a research-article  |2 jats 
949 |B NATIONALLICENCE  |F NATIONALLICENCE  |b NL-springer 
950 |B NATIONALLICENCE  |P 856  |E 40  |u https://doi.org/10.1007/s00214-014-1600-8  |q text/html  |z Onlinezugriff via DOI 
950 |B NATIONALLICENCE  |P 700  |E 1-  |a Reif  |D Maria  |u Institute for Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Muthgasse 18, 1190, Vienna, Austria  |4 aut 
950 |B NATIONALLICENCE  |P 700  |E 1-  |a Oostenbrink  |D Chris  |u Institute for Molecular Modeling and Simulation, University of Natural Resources and Life Sciences, Vienna, Muthgasse 18, 1190, Vienna, Austria  |4 aut 
950 |B NATIONALLICENCE  |P 773  |E 0-  |t Theoretical Chemistry Accounts  |d Springer Berlin Heidelberg  |g 134/2(2015-02-01), 1-19  |x 1432-881X  |q 134:2<1  |1 2015  |2 134  |o 214