Monte Carlo algorithms for simulating systems with adiabatic separation of electronic and nuclear degrees of freedom

Research output: Contribution to journalArticlepeer-review

29 Scopus citations

Abstract

Two Monte Carlo algorithms for the adiabatic sampling of nuclear and electronic degrees of freedom are introduced. In these algorithms the electronic degrees of freedom are subject to a secondary low-temperature thermostat in close analogy to the extended Lagrangian formalism used in molecular dynamics simulations. Numerical tests are carried out for two model systems of coupled harmonic oscillators, and for two more realistic systems: the water dimer and bulk liquid water. A statistical-mechanical discussion of the partition function for systems with adiabatic separation of electronic and nuclear degrees of freedom, but with finite electronic temperature, is presented. The theoretical analysis shows that the algorithms satisfy the adiabatic limit using suitable choices of the electronic temperature, Telec, the number of electronic moves, Relec, and the maximum step sizes used for displacements of nuclear coordinates. For quadratic coupling of the nuclear and electronic degrees of freedom, the electronic phase-space volume is independent of the nuclear coordinates. In this case, the sampling of the nuclear coordinate phase-space recovers the correct Born-Oppenheimer limit independent of Telec, but each electronic degree of freedom contributes an offset of 0.5kBT′elec (with T′-1elec = T-1 + T-1elec) to the average total energy. For nonquadratic coupling, satisfactory sampling of the nuclear coordinate phase-space requires a low Telec to limit the ratio of the electronic phase-space volumes at Telec and T′elec to be close to unity.

Original languageEnglish (US)
Pages (from-to)87-104
Number of pages18
JournalTheoretical Chemistry Accounts
Volume103
Issue number2
DOIs
StatePublished - Dec 1999

Keywords

  • Adiabatic sampling
  • Monte Carlo algorithms
  • Polarizable force fields

Fingerprint Dive into the research topics of 'Monte Carlo algorithms for simulating systems with adiabatic separation of electronic and nuclear degrees of freedom'. Together they form a unique fingerprint.

Cite this