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 -