TY - JOUR

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

AU - Chen, Bin

AU - Siepmann, J. Ilja

PY - 1999/12

Y1 - 1999/12

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

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

KW - Adiabatic sampling

KW - Monte Carlo algorithms

KW - Polarizable force fields

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

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

U2 - 10.1007/s002140050519

DO - 10.1007/s002140050519

M3 - Article

AN - SCOPUS:0033261015

SN - 1432-881X

VL - 103

SP - 87

EP - 104

JO - Theoretical Chemistry Accounts

JF - Theoretical Chemistry Accounts

IS - 2

ER -