Path integral formulation: Difference between revisions
Carl McBride (talk | contribs) (Further tidy.) |
Carl McBride (talk | contribs) m (Removed duplicated reference.) |
||
Line 16: | Line 16: | ||
:<math>\Phi_P (x_1...x_P;\beta)= \frac{mP}{2\beta^2 \hbar^2} \sum_{t=1}^P (x_t - x_{t+1})^2 + \frac{1}{P} \sum_{t=1}^P V(x_t)</math>. | :<math>\Phi_P (x_1...x_P;\beta)= \frac{mP}{2\beta^2 \hbar^2} \sum_{t=1}^P (x_t - x_{t+1})^2 + \frac{1}{P} \sum_{t=1}^P V(x_t)</math>. | ||
where <math>P</math> is the Trotter number. In the Trotter limit, where <math>P \rightarrow \infty</math> these equations become exact. In the case where <math>P=1</math> these equations revert to a classical simulation. It has long been recognised that there is an isomorphism between this discretised quantum mechanical description, and the classical [[statistical mechanics]] of polyatomic fluids, in particular flexible ring molecules (Ref. | where <math>P</math> is the Trotter number. In the Trotter limit, where <math>P \rightarrow \infty</math> these equations become exact. In the case where <math>P=1</math> these equations revert to a classical simulation. It has long been recognised that there is an isomorphism between this discretised quantum mechanical description, and the classical [[statistical mechanics]] of polyatomic fluids, in particular flexible ring molecules (Ref. 3), due to the periodic boundary conditions in imaginary time. | ||
==Rotational degrees of freedom== | ==Rotational degrees of freedom== | ||
In the case of systems having (<math>d</math>) rotational [[degree of freedom | degrees of freedom]] the [[Hamiltonian]] can be written in the form (Ref. | In the case of systems having (<math>d</math>) rotational [[degree of freedom | degrees of freedom]] the [[Hamiltonian]] can be written in the form (Ref. 8 Eq. 2.1): | ||
:<math>\hat{H} = \hat{T}^{\mathrm {translational}} + \hat{T}^{\mathrm {rotational}}+ \hat{V}</math> | :<math>\hat{H} = \hat{T}^{\mathrm {translational}} + \hat{T}^{\mathrm {rotational}}+ \hat{V}</math> | ||
where the rotational part of the kinetic energy operator is given by (Ref. | where the rotational part of the kinetic energy operator is given by (Ref. 8 Eq. 2.2) | ||
:<math>T^{\mathrm {rotational}} = \sum_{i=1}^{d^{\mathrm {rotational}}} \frac{\hat{L}_i^2}{2\Theta_{ii}}</math> | :<math>T^{\mathrm {rotational}} = \sum_{i=1}^{d^{\mathrm {rotational}}} \frac{\hat{L}_i^2}{2\Theta_{ii}}</math> | ||
where <math>\hat{L}_i</math> are the components of the angular momentum operator, and <math>\Theta_{ii}</math> are the moments of inertia. | where <math>\hat{L}_i</math> are the components of the angular momentum operator, and <math>\Theta_{ii}</math> are the moments of inertia. | ||
==Techniques== | ==Techniques== | ||
====Path integral Monte Carlo==== | ====Path integral Monte Carlo==== | ||
Line 65: | Line 62: | ||
#[http://dx.doi.org/10.1080/014423597230190 Charusita Chakravarty "Path integral simulations of atomic and molecular systems", International Reviews in Physical Chemistry '''16''' pp. 421-444 (1997)] | #[http://dx.doi.org/10.1080/014423597230190 Charusita Chakravarty "Path integral simulations of atomic and molecular systems", International Reviews in Physical Chemistry '''16''' pp. 421-444 (1997)] | ||
# M. J. Gillan "The path-integral simulation of quantum systems" in "Computer Modelling of Fluids Polymers and Solids" eds. C. R. A. Catlow, S. C. Parker and M. P. Allen, NATO ASI Series C '''293''' pp. 155-188 (1990) ISBN 978-0-7923-0549-1 | # M. J. Gillan "The path-integral simulation of quantum systems" in "Computer Modelling of Fluids Polymers and Solids" eds. C. R. A. Catlow, S. C. Parker and M. P. Allen, NATO ASI Series C '''293''' pp. 155-188 (1990) ISBN 978-0-7923-0549-1 | ||
#[http://dx.doi.org/10. | #[http://dx.doi.org/10.1088/0953-8984/11/11/003 Dominik Marx and Martin H Müser "Path integral simulations of rotors: theory and applications", Journal of Physics: Condensed Matter '''11''' pp. R117-R155 (1999)] | ||
[[Category: Monte Carlo]] | [[Category: Monte Carlo]] | ||
[[category: Quantum mechanics]] | [[category: Quantum mechanics]] |
Revision as of 15:25, 22 October 2008
The Path integral formulation is an elegant method by which quantum mechanical contributions can be incorporated within a classical simulation using Feynman path integrals (Refs. 1-7)
Principles
In the path integral formulation the canonical partition function (in one dimension) is written as (Ref. 4 Eq. 1)
where is the Euclidian action, given by (Ref. 4 Eq. 2)
where is the path in time and is the Hamiltonian. This leads to (Ref. 4 Eq. 3)
where the Euclidean time is discretised in units of
and ((Ref. 4 Eq. 4)
- .
where is the Trotter number. In the Trotter limit, where these equations become exact. In the case where these equations revert to a classical simulation. It has long been recognised that there is an isomorphism between this discretised quantum mechanical description, and the classical statistical mechanics of polyatomic fluids, in particular flexible ring molecules (Ref. 3), due to the periodic boundary conditions in imaginary time.
Rotational degrees of freedom
In the case of systems having () rotational degrees of freedom the Hamiltonian can be written in the form (Ref. 8 Eq. 2.1):
where the rotational part of the kinetic energy operator is given by (Ref. 8 Eq. 2.2)
where are the components of the angular momentum operator, and are the moments of inertia.
Techniques
Path integral Monte Carlo
Path integral Monte Carlo (PIMC)
Path integral molecular dynamics
Path integral molecular dynamics (PIMC)
Centroid molecular dynamics
Centroid molecular dynamics (CMD)
- Jianshu Cao and Gregory A. Voth "The formulation of quantum statistical mechanics based on the Feynman path centroid density. II. Dynamical properties", Journal of Chemical Physics 100 pp. 5106- (1994)
- Seogjoo Jang and Gregory A. Voth "A derivation of centroid molecular dynamics and other approximate time evolution methods for path integral centroid variables", Journal of Chemical Physics 111 pp. 2371- (1999)
Ring polymer molecular dynamics
Ring polymer molecular dynamics (RPMD)
- Ian R. Craig and David E. Manolopoulos "Quantum statistics and classical mechanics: Real time correlation functions from ring polymer molecular dynamics", Journal of Chemical Physics 121 pp. 3368- (2004)
- Bastiaan J. Braams and David E. Manolopoulos "On the short-time limit of ring polymer molecular dynamics", Journal of Chemical Physics 125 124105 (2006)
- Thomas E. Markland and David E. Manolopoulos "An efficient ring polymer contraction scheme for imaginary time path integral simulations", Journal of Chemical Physics 129 024105 (2008)
Grand canonical Monte Carlo
A path integral version of the Widom test-particle method for grand canonical Monte Carlo simulations:
Applications
Phase transitions, quantum dynamics, centroids etc.
- J. R. Melrose and K. Singer "An investigation of supercooled Lennard-Jones argon by quantum mechanical and classical Monte Carlo simulation", Molecular Physics 66 1203-1214 (1989)
- Jianshu Cao and Gregory A. Voth "The formulation of quantum statistical mechanics based on the Feynman path centroid density. I. Equilibrium properties", Journal of Chemical Physics 100 pp. 5093-5105 (1994)
- Jianshu Cao and Gregory A. Voth "Semiclassical approximations to quantum dynamical time correlation functions", Journal of Chemical Physics 104 pp. 273-285 (1996)
- Rafael Ramírez and Telesforo López-Ciudad "The Schrödinger formulation of the Feynman path centroid density", Journal of Chemical Physics 111 pp. 3339-3348 (1999)
- C. Chakravarty and R. M. Lynden-Bell "Landau free energy curves for melting of quantum solids", Journal of Chemical Physics 113 pp. 9239-9247 (2000)
External links
- Density matrices and path integrals computer code on SMAC-wiki.
References
- R. P. Feynman and A. R. Hibbs "Path-integrals and Quantum Mechanics", McGraw-Hill, New York (1965) ISBN 0-07-020650-3
- R. P. Feynman "Statistical Mechanics", Benjamin, Reading, Massachusetts, (1972) ISBN 0805325085
- David Chandler and Peter G. Wolynes "Exploiting the isomorphism between quantum theory and classical statistical mechanics of polyatomic fluids", Journal of Chemical Physics 74 pp. 4078-4095 (1981)
- B. J. Berne and D. Thirumalai "On the Simulation of Quantum Systems: Path Integral Methods", Annual Review of Physical Chemistry 37 pp. 401-424 (1986)
- D. M. Ceperley "Path integrals in the theory of condensed helium", Reviews of Modern Physics 67 279 - 355 (1995)
- Charusita Chakravarty "Path integral simulations of atomic and molecular systems", International Reviews in Physical Chemistry 16 pp. 421-444 (1997)
- M. J. Gillan "The path-integral simulation of quantum systems" in "Computer Modelling of Fluids Polymers and Solids" eds. C. R. A. Catlow, S. C. Parker and M. P. Allen, NATO ASI Series C 293 pp. 155-188 (1990) ISBN 978-0-7923-0549-1
- Dominik Marx and Martin H Müser "Path integral simulations of rotors: theory and applications", Journal of Physics: Condensed Matter 11 pp. R117-R155 (1999)