English

An ab initio path integral Monte Carlo simulation method for molecules and clusters: application to Li_4 and Li_5^+

Computational Physics 2009-10-31 v1 Statistical Mechanics Atomic and Molecular Clusters Chemical Physics

Abstract

A novel method for simulating the statistical mechanics of molecular systems in which both nuclear and electronic degrees of freedom are treated quantum mechanically is presented. The scheme combines a path integral description of the nuclear variables with a first-principles adiabatic description of the electronic structure. The electronic problem is solved for the ground state within a density functional approach, with the electronic orbitals expanded in a localized (Gaussian) basis set. The discretized path integral is computed by a Metropolis Monte Carlo sampling technique on the normal modes of the isomorphic ring-polymer. An effective short-time action correct to order τ4\tau^4 is used. The validity and performance of the method are tested in two small Lithium clusters, namely Li4_4 and Li5+_5^+. Structural and electronic properties computed within this fully quantum-mechanical scheme are presented and compared to those obtained within the classical nuclei approximation. Quantum delocalization effects are significant but tunneling turns out to be irrelevant at low temperatures.

Keywords

Cite

@article{arxiv.physics/9802020,
  title  = {An ab initio path integral Monte Carlo simulation method for molecules and clusters: application to Li_4 and Li_5^+},
  author = {Ruben O. Weht and Jorge Kohanoff and Dario A. Estrin and Charusita Chakravarty},
  journal= {arXiv preprint arXiv:physics/9802020},
  year   = {2009}
}

Comments

11 text pages, 7 figures, to be published in J. Chem. Phys

R2 v1 2026-07-22T19:17:29.995Z