We present a fully Hamiltonian and computationally efficient scheme to include the electrostatic effects due to the classical environment in a Car-Parrinello mixed quantum Mechanics/molecular mechanics (QM/MM) method. The polarization due to the MM atoms close to the quantum system is described by a Coulombic potential modified at short range. We show that the functional form of this potential has to be chosen carefully in order to obtain the correct interaction properties and to prevent an unphysical escape of the electronic density to the MM atoms (the so-called spill-out effect). The interaction between the QM system and the more distant MM atoms is modeled by a Hamiltonian term explicitly coupling the multipole moments of the quantum charge distribution with the classical point charges. Our approach remedies some of the well known deficiencies of current electrostatic coupling schemes in QM/MM methods, allowing molecular dynamics simulations of mixed systems within a fully consistent and energy conserving approach.

A Hamiltonian electrostatic coupling scheme for hybrid Car-Parrinello molecular dynamics simulations / Laio, A.; Vandevondele, J.; Rothlisberger, U.. - In: THE JOURNAL OF CHEMICAL PHYSICS. - ISSN 0021-9606. - 116:16(2002), pp. 6941-6947. [10.1063/1.1462041]

A Hamiltonian electrostatic coupling scheme for hybrid Car-Parrinello molecular dynamics simulations

Laio, A.;
2002-01-01

Abstract

We present a fully Hamiltonian and computationally efficient scheme to include the electrostatic effects due to the classical environment in a Car-Parrinello mixed quantum Mechanics/molecular mechanics (QM/MM) method. The polarization due to the MM atoms close to the quantum system is described by a Coulombic potential modified at short range. We show that the functional form of this potential has to be chosen carefully in order to obtain the correct interaction properties and to prevent an unphysical escape of the electronic density to the MM atoms (the so-called spill-out effect). The interaction between the QM system and the more distant MM atoms is modeled by a Hamiltonian term explicitly coupling the multipole moments of the quantum charge distribution with the classical point charges. Our approach remedies some of the well known deficiencies of current electrostatic coupling schemes in QM/MM methods, allowing molecular dynamics simulations of mixed systems within a fully consistent and energy conserving approach.
2002
116
16
6941
6947
Laio, A.; Vandevondele, J.; Rothlisberger, U.
File in questo prodotto:
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/20.500.11767/12853
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 574
  • ???jsp.display-item.citation.isi??? 538
social impact