# What are the types of Quantum Molecular Dynamics (QMD)?

Matter Modeling Asked on December 4, 2021

In similar spirit to recent questions on Quantum Monte Carlo, ab-initio Molecular Dynamics, types of SCF, and others, I would like to ask:

What are the types of Quantum Molecular Dynamics (QMD)?

As I have recently learned, Quantum Molecular Dynamics is different to AIMD and MD through the fact that in QMD, the nuclei are treated quantum mechanically.

I would like to ask here, if we can summarize the types of QMD in a few paragraphs:

• PIMD (Feynman Path Integral Molecular Dynamics)
• CMD (Centroid Molecular Dynamics)
• RPMD (Ring Polymer Molecular Dynamics) [link to answer]
• *Feel free to let me know of other methods or include them as an answer!

Also, it would be great if we limit one QMD method per answer (and per person).

## RPMD (Ring Polymer Molecular Dynamics)

Introduction (correlation functions and Kubo transforms):

For a time-dependent quantum operator $$hat{A}(t)$$ the auto-correlation function of the operator at time $$t$$ versus at time 0 can be written as:

$$c_{AA}(t) equiv frac{textrm{tr}left( e^{-beta H}hat{A}(0)hat{A}(t) right)}{textrm{tr}left({e^{-beta H}}right)} tag{1}.$$

Applying the Kubo transformation to this, gives the Kubo-transformed auto-correlation function:

$$tilde{c}_{AA}(t) equiv frac{int_0^beta textrm{tr}left( e^{-(beta - lambda) H}hat{A}(0) e^{-lambda H }hat{A}(t) right)textrm{d}lambda}{beta, textrm{tr}left({e^{-beta H}}right)} tag{2}.$$

Remember from high school the chemical reaction rate $$k$$. It can be written in terms of the above two expressions! Reaction rates depend on temperature so we'll write $$k(T)$$:

begin{align} k(t) &=frac{1}{Q_r(T)}int_0^infty c_{ff}(t)textrm{d}t tag{3}\ &=frac{1}{Q_r(T)}int_0^infty tilde{c}_{ff}(t)textrm{d}t. tag{4} end{align}

Review of PIMD (Path Integral MD):

If we have $$N$$ atoms, and we treat each of them as an $$n$$-bead ring of artificial atoms (beads) we can make the approximation:

$${small textrm{tr}left(e^{-beta H}right) approx frac{1}{left( 2pi hbar right)^{3Nn}}int !!!! int cdots int e^{-frac{beta}{n} H_nleft(mathbf{p}_1,ldots,mathbf{p}_{Nn},mathbf{q}_1,ldots,mathbf{q}_{Nn}right)} textrm{d}^3mathbf{p}_1ldots textrm{d}^3mathbf{p}_n textrm{d}^3mathbf{q}_1ldots textrm{d}^3mathbf{q}_ntag{5}, }$$

where $$H_n$$ is the Hamiltonian of the $$N times n$$ beads representing $$N$$ atoms and the $$n$$ beads connected by harmonic spring potentials representing each of the $$N$$ atoms.

This approximation becomes exact when $$nrightarrow infty$$, which would mean we have an infinite-dimensional integral (known as a Feynman integral or "path integral", in this case actually a double-Feynman-integral or Feynman double-integral since there's two entirely different "path" sets over which a Feynman integral is being done).

The "RPMD" approximation:

It might sound bizarre, because really PIMD is MD on "ring polymers" (ring polymer just being another name for the set of beads representing each atom), so PIMD could be called RPMD. However when people use the term "RPMD" they are referring to this approximation:

begin{align}{tiny !!!!!!!!tilde{c}_{AA}(t) approx frac{int !!! int cdots int e^{-frac{beta}{n} H_nleft(mathbf{p}_1(t),ldots,mathbf{p}_{Nn}(t),mathbf{q}_1(t),ldots,mathbf{q}_{Nn}(t)right)}hat{A}_nleft( mathbf{q}_1(0),ldots,mathbf{q}_{Nn}(0)right) hat{A}_n left( mathbf{q}_1(t),ldots,mathbf{q}_{Nn}(t) right) textrm{d}^3mathbf{p}_1(0)ldots textrm{d}^3mathbf{p}_n(0) textrm{d}^3mathbf{q}_1(0)ldots textrm{d}^3mathbf{q}_n(0)}{int !!! int cdots int e^{-frac{beta}{n} H_nleft(mathbf{p}_1,ldots,mathbf{p}_{Nn},mathbf{q}_1,ldots,mathbf{q}_{Nn}right)} textrm{d}^3mathbf{p}_1ldots textrm{d}^3mathbf{p}_n textrm{d}^3mathbf{q}_1ldots textrm{d}^3mathbf{q}_n }} tag{6}, end{align}

where $$hat{A}_nleft( mathbf{q}_1(t),ldots,mathbf{q}_{Nn}(t)right)$$ involves for each of the $$N$$ atoms, an average over all $$n$$ of its beads at time $$t$$:

$${small hat{A}_nleft( mathbf{q}_1(t),ldots,mathbf{q}_{Nn}(t)right) equiv frac{1}{n}sum_{j=1}^n hat{A}left( mathbf{q}_{j}(t),mathbf{q}_{n+j}(t),mathbf{q}_{2n+j}(t),ldots ,mathbf{q}_{(N-1)n+j}(t) right).tag{7} }$$

Unlike PIMD which calculates a static property exactly in the limit of $$nrightarrow infty$$, RPMD approximates a function of time (the auto-correlation function) even in the limit as $$nrightarrow infty$$. However the approximation is exact (in the limit where $$nrightarrow infty$$) in some limits:

• infinitely high temperature,
• $$t approx 0$$, (short-time limit),
• harmonic limit (where the $$N$$ atoms interact via harmonic potentials) if $$A$$ is linear,
• $$hat{A}=1$$ (the identity operator).

Pros:

• It is relatively cheap to calculate (compared to exact real-time quantum dynamics). It costs only a bit more than doing classical MD several ($$n$$) times,
• it allows one to incorporate the effects of tunneling and zero-point-energy effects in an MD calculation.

Cons:

• "coherence" between paths in the Feynman integral is not taken into account,
• the formulation above does not give attention to how non-adiabatic effects would be treated. Several proposals have bee made for how to treat non-adiabatic effects, but no single proposal has stood out as a "gold standard" as far as I know.
• it's at least $$n$$ times more expensive than doing classical MD, so it's applicability is limited to cases where quantum mechanical effects of nuclei are important enough for it to be worth the extra cost.

Answered by Nike Dattani on December 4, 2021

## Related Questions

### What are typical RPA capabilities of plane-wave codes?

1  Asked on August 19, 2021 by michael-f-herbst

### Quadrature over three Euler Angles for orientation averaging

1  Asked on August 19, 2021 by emil-zak

### What are some open-source all-electron DFT alternatives to Wien2K?

4  Asked on August 19, 2021

### Real space projection vs reciprocal space projection in DFT calculations

1  Asked on August 19, 2021

### How to calculate the Fock matrix in the molecular orbital basis PySCF?

1  Asked on August 19, 2021 by wychh

### Bond Order: When and how is it used today?

2  Asked on August 19, 2021

### Is it possible to calculate/estimate the value of the J parameter to be used in the Heisenberg/Ising Hamiltonians?

1  Asked on August 19, 2021

### What are some methods for modeling bulk phase infrared spectra?

2  Asked on August 19, 2021

### What is the closest thing we have to “the” universal density functional?

4  Asked on August 19, 2021

### Molden installation on macOS?

1  Asked on August 19, 2021 by dmitry-eremin

### What are the types of bond orders?

4  Asked on August 19, 2021

### What is band inversion and how to recognize it in band structure?

1  Asked on August 19, 2021

### What are the types of charge analysis?

10  Asked on August 19, 2021 by dgkang

### What are the types of DFT?

5  Asked on August 19, 2021

### Supercomputers around the world

14  Asked on August 19, 2021

### What are the types of pseudopotentials?

2  Asked on August 19, 2021

### What are the types of ab initio Molecular Dynamics?

3  Asked on August 19, 2021 by etienne-palos

### Are there any MD packages that do proper free energy sampling in the NPT ensemble?

1  Asked on August 19, 2021 by b-kelly

### VASP vs MedeA-VASP

1  Asked on August 19, 2021

### Derivation of Slater-Koster equations

2  Asked on August 19, 2021

### Ask a Question

Get help from others!