Premium
Efficiency of nested Markov chain Monte Carlo for polarizable potentials and perturbed Hamiltonians
Author(s) -
Calvo Florent
Publication year - 2010
Publication title -
international journal of quantum chemistry
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 0.484
H-Index - 105
eISSN - 1097-461X
pISSN - 0020-7608
DOI - 10.1002/qua.22527
Subject(s) - markov chain monte carlo , statistical physics , monte carlo method , hamiltonian (control theory) , quantum monte carlo , hybrid monte carlo , markov chain , quantum , free energy perturbation , physics , parallel tempering , mathematics , quantum mechanics , mathematical optimization , molecular dynamics , statistics
Abstract Nested Markov chain Monte Carlo is a rigorous way to enhance sampling of a given energy landscape using an auxiliary, approximate potential energy surface. Its practical efficiency mainly depends on how cheap and how different are the auxiliary potential with respect to the reference system. In this article, a combined efficiency index is proposed and assessed for two important families of energy surfaces. As illustrated for water clusters, many‐body polarizable potentials can be approximated by simplifying the polarization contribution and keeping only the two‐body terms. In small systems, neglecting polarization entirely is also acceptable. When the reference potential energy is obtained from diagonalization of a quantum mechanical Hamiltonian, a first‐order perturbation scheme can be used to estimate the energy difference occuring on a Monte Carlo move. Our results indicate that this perturbation approximation performs well provided that the number of steps between successive diagonalization is adjusted beforehand. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem 110:2342–2346, 2010