---
title: Second-generation Car–Parrinello Methods
url: https://www.emergentmind.com/topics/second-generation-car-parrinello-methods
type: topic
---

# Second-generation Car–Parrinello Methods

Searching arXiv for the provided topic and cited papers.
Second-generation Car–Parrinello methods are ab-initio molecular dynamics formalisms devised to occupy the intermediate ground between conventional Born–Oppenheimer molecular dynamics and original Car–Parrinello molecular dynamics. In the narrow CP2K sense, the term denotes a predictor–corrector electronic propagation scheme that keeps the electronic state very close to the instantaneous Born–Oppenheimer ground state while requiring only one preconditioned gradient calculation per AIMD step, thereby providing a Car–Parrinello-like route to Born–Oppenheimer-quality dynamics without full self-consistent minimization at every step [2205.14741]. In broader usage, the label also encompasses Car–Parrinello-inspired extended-Lagrangian Born–Oppenheimer frameworks that replace direct SCF-driven BOMD by reversible auxiliary electronic dynamics on a shadow Born–Oppenheimer surface [1705.10845]. Across these variants, the common objective is to retain large BOMD-like timesteps, eliminate or sharply reduce iterative electronic optimization, and avoid the fictitious-mass limitations of first-generation CPMD [1303.2067] [1201.5945].

## 1. Conceptual origin and historical placement

The central motivation is the classical AIMD tradeoff between electronic relaxation accuracy and time-to-solution. In conventional BOMD, the electronic structure is minimized to the instantaneous ground state at every nuclear step, typically through an SCF procedure. The reviews emphasize that this makes BOMD computationally expensive, especially when long trajectories and large cells are required for slowly converging observables such as diffusion, viscosity, and hydrogen-bond kinetics in liquids [1303.2067] [1201.5945].

Original CPMD solves a different problem and creates a new one. By introducing fictitious dynamics for the electronic variables, it avoids full electronic minimization at each step, but accuracy then depends on a fictitious electron mass parameter and on adiabatic separation between ionic and electronic motion. The resulting timestep must be smaller than in BOMD, and metallic or small-gap systems are particularly difficult because the admissible timestep and BO proximity depend on the electronic gap [1201.5945] [1705.10845].

Second-generation methods were therefore designed to combine the best features of both approaches: large timesteps as in BOMD, no expensive SCF minimization per step as in CPMD, and no fictitious electron mass parameter in the predictor–corrector branch [1303.2067]. The 2017 extended-Lagrangian formulation explicitly presents itself as a successor to CPMD, “in the spirit of Car–Parrinello molecular dynamics,” but with a different Lagrangian, higher-order accuracy, and no orthonormalization or idempotency constraints on the auxiliary extended variables [1705.10845]. This terminological breadth explains why “second-generation Car–Parrinello” may refer either to the Kühné-style predictor–corrector CPMD used in CP2K or, more broadly, to Car–Parrinello-inspired extended-Lagrangian BOMD.

## 2. Predictor–corrector second-generation CPMD

The canonical predictor–corrector formulation exploits the smooth time evolution of the occupied electronic subspace. In the review treatment, the one-electron density operator is written as
$$
\rho = \sum_i |\psi_i\rangle \langle \psi_i |,
$$
and the next-step occupied subspace is extrapolated from the previous $K$ steps rather than recomputed from scratch [1303.2067]. The predictor is
$$
|\psi_{i}^{p}(t_{n}) \rangle = \sum_{m=1}^{K} \rho(t_{n-m}) |\psi_{i}(t_{n-1}) \rangle \alpha_{m},
$$
with
$$
\alpha_{m} = (-1)^{m+1} m \frac{\binom{2K}{K-m}}{\binom{2K-2}{K-1}},
$$
followed by the corrector
$$
|\psi_{i}(t_{n}) \rangle = (1 - \omega) |\psi_{i}^{p}(t_{n}) \rangle + \omega | \delta \psi_{i}^{p}(t_{n}) \rangle,
$$
where
$$
\omega = \frac{K}{2K-1}, \qquad K \ge 2.
$$
The correction is a single preconditioned electronic gradient step, commonly implemented through the orbital transformation method [1303.2067].

The same idea can be expressed in CP2K’s matrix language. If $\mathbf C$ is the molecular-orbital coefficient matrix, then
$$
\mathbf P = \mathbf C \mathbf C^T,
$$
and, for a non-orthogonal basis, the fermionic constraint is
$$
\mathbf P = \mathbf P \mathbf S \mathbf P.
$$
The key observation is that $\mathbf P\mathbf S$ evolves more smoothly than the raw orbitals $\mathbf C$, so the ASPC predictor advances the occupied subspace using previous density-kernel information rather than direct orbital extrapolation [1201.5945]. In the CP2K perspective, this is the defining operational reduction: electrons are advanced dynamically and corrected once, with only one preconditioned electronic gradient per ionic step, while remaining very close to the instantaneous Kohn–Sham ground state [2205.14741].

Time-reversibility is approximate but high order. The reviews state that the electron propagation is time reversible up to $\mathcal{O}(\Delta t^{2K-2})$; for $K=7$, this becomes $\mathcal{O}(\Delta t^{12})$ [1303.2067]. In an earlier CP2K implementation discussed for liquid silicon, silica, and water, $K=4$ was used, corresponding to reversibility up to $O(\Delta t^6)$ [1201.5945].

## 3. Energy, forces, and constraint handling

Second-generation CPMD does not evaluate the exact BO energy at each step. Instead, it uses an approximation to the Harris–Foulkes functional built from the predicted density:
$$
\begin{aligned}
E^{PC}[\rho^{p}] &= \sum_{i} \langle \psi_{i} | \hat{H}[\rho^{p}] |\psi_{i}\rangle - \frac{1}{2} \int d\mathbf{r} \int d\mathbf{r}' \frac{\rho^{p}(\mathbf{r}) \rho^{p}(\mathbf{r}')}{|\mathbf{r} - \mathbf{r}'|} \\
&\quad - \int d\mathbf{r} \, v_{XC}[\rho^{p}] \rho^{p} + E_{XC}[\rho^{p}] + E_{II}.
\end{aligned}
$$
Here the Hamiltonian is built from the predicted density $\rho^p$, while the orbitals entering the band-structure term are the corrected orbitals [1303.2067].

The forces are analytic gradients of this predictor–corrector functional rather than exact Hellmann–Feynman forces on a fully converged BO surface. Because the corrected density $\rho$ differs from the predicted density $\rho^p$,
$$
\Delta \rho \equiv \rho - \rho^p \neq 0,
$$
the Hellmann–Feynman and Pulay contributions must be supplemented by an extra term,
$$
-\int d\mathbf r \left\{ \left[ \left( \frac{\partial V_{\mathrm{XC}}[\rho^p]}{\partial \rho^p} \right)\Delta\rho + V_{\mathrm{H}}[\Delta\rho] \right] \left(\nabla_{\mathbf R_I}\rho^p\right) \right\},
$$
which is specific to the use of a predicted, rather than fully self-consistent, density [1201.5945].

Constraint handling is moved out of fictitious orbital dynamics and into the minimization geometry. The CP2K implementation uses the orbital transformation method, parameterizing corrected orbitals as
$$
\mathbf{C}(\mathbf X) = \mathbf C^p(t_n)\cos(\mathbf U) + \mathbf X \mathbf U^{-1}\sin(\mathbf U),
$$
with
$$
\mathbf U = (\mathbf X^T \mathbf S \mathbf X)^{1/2},
$$
and tangent-space constraint
$$
\mathbf X^T \mathbf S \mathbf C^p(t_n)=0.
$$
If the predicted orbitals are orthonormal, this parameterization yields an idempotent density matrix for any $\mathbf X$ satisfying the tangent-space constraint. The reviews note, however, that ASPC only approximately preserves idempotency, so explicit purification iterations may occasionally be necessary [1201.5945].

## 4. Dissipation, modified Langevin sampling, and dynamical fidelity

A defining caveat of the predictor–corrector branch is that it is not exact BO propagation. The reviews state explicitly that the scheme is slightly dissipative over long times, probably because the electron propagator is not symplectic [1303.2067] [1201.5945]. The standard Langevin equation is written as
$$
M_{I} \ddot{\mathbf R}_{I} = \mathbf F_{I}^{BO} - \gamma M_{I}\dot{\mathbf R}_{I} + \mathbf{\Xi}_{I},
$$
with fluctuation–dissipation relation
$$
\left< \mathbf{\Xi}_{I}(0) \mathbf{\Xi}_{I}(t) \right> = 2 \gamma k_{B} T M_{I} \delta(t).
$$
The intrinsic dissipation of second-generation CPMD is modeled by
$$
\mathbf F_{I}^{PC} = \mathbf F_{I}^{BO} - \gamma_{D}M_{I}\dot{\mathbf R}_{I},
$$
which yields the practical modified Langevin-like equation
$$
M_{I} \ddot{\mathbf R}_{I} = \mathbf F_{I}^{PC} + \mathbf{\Xi}_{I}.
$$
This construction is used to recover exact canonical sampling despite the dissipative character of the deterministic predictor–corrector propagation [1303.2067].

The intrinsic damping coefficient $\gamma_D$ is not known a priori. The liquid-water review states that it can be determined in a preliminary run, inspired by Krajewski and Parrinello, using a Berendsen-like algorithm until equipartition holds,
$$
\left< \frac{1}{2} M_{I} \dot{\mathbf R}_{I}^{2} \right> = \frac{3}{2} k_{B} T.
$$
For the 128-water simulations reported there, the value used was
$$
\gamma_{D} = 8.65 \times 10^{-5}\ \mathrm{fs}^{-1}.
$$
The later practical guide generalizes this viewpoint in CP2K syntax through an intrinsic `NOISY_GAMMA`, an optional overlay `GAMMA`, and species-wise temperature checks as a stringent sampling diagnostic [1303.2067] [2601.12191].

The principal misconception addressed by these papers is that second-generation CPMD is simply BOMD with a better initial guess. It is not exact BO minimization at every step, and its energy conservation is not that of ideal microcanonical BO dynamics. Conversely, it is also not original CPMD with a modified fictitious mass, because the defining electronic update is predictor–corrector propagation with one gradient correction and no fictitious mass parameter [1303.2067] [1201.5945].

## 5. Performance, implementation in CP2K, and representative applications

Within CP2K, second-generation Car–Parrinello is implemented in Quickstep, the electronic-structure engine based on the Gaussian and plane wave method and its all-electron GAPW extension. The 2022 CP2K perspective presents SGCP as an improved coupled electron–ion dynamics scheme and as a “unique selling point” of CP2K because it combines efficient GPW/GAPW electronic-structure evaluation, analytic forces, and propagation that avoids repeated SCF minimization [2205.14741].

The performance claims are qualitative but substantial. The 2022 perspective states that CP2K can “routinely conduct nanosecond long DFT-based AIMD simulations with thousands of atoms” because SGCP requires only one preconditioned gradient calculation per AIMD step [2205.14741]. The 2012 review states that the superior efficiency is between one and two orders of magnitude, depending on the system, and that the method makes it possible to simulate medium-sized systems up to a few thousand atoms for as long as a couple of nanoseconds [1201.5945]. The 2013 water review attributes efficiency gains of one to two orders of magnitude to the absence of SCF cycles and diagonalization, together with large BOMD-like timesteps [1303.2067].

Water is the most extensively documented validation case. The review reports direct comparison with BOMD for 128-water simulations, with radial distribution functions, velocity autocorrelation functions, and vibrational densities of states in excellent agreement with BOMD references; the preconditioned mean gradient deviation from the BO surface was $10^{-5}$ a.u., only slightly larger than typical fully converged BOMD values [1303.2067]. In that review, simulation cells with up to 128 molecules and trajectories totaling more than 1 ns of AIMD were used to estimate self-diffusion, shear viscosity, and hydrogen-bond lifetime and relaxation [1303.2067].

The practical 2026 guide turns these ideas into a CP2K workflow. It recommends short BOMD pre-equilibration, explicit storage of density-matrix history, restart from that history, and short NVE tuning runs to minimize energy dissipation. In CP2K, the predictor is activated by `EXTRAPOLATION ASPC`, the number of electronic correction steps is controlled through `MAX_SCF_HISTORY`, and the practical corrector strength is set by OT `STEPSIZE` rather than the idealized $\omega$ of the 2007 theory [2601.12191]. The guide states that values between 1 and 2 usually suffice for `MAX_SCF_HISTORY`, and that the final optimum for `EXTRAPOLATION_ORDER` often lies between 0 and 3 [2601.12191]. Its worked example uses 32 water molecules at $T=300\ \mathrm{K}$ and $P=1\ \mathrm{bar}$ in a cubic cell with $a=9.8528\ \mathrm{\AA}$, PBE, GTH pseudopotentials, DZVP-GTH basis sets, `CUTOFF 240`, timestep $0.5$ fs, `EXTRAPOLATION_ORDER 3`, and `STEPSIZE 0.15` [2601.12191].

## 6. Extended-Lagrangian and related variants

A broader branch of second-generation Car–Parrinello methodology is extended-Lagrangian Born–Oppenheimer molecular dynamics. In this formulation, the exact universal functional is linearized around an auxiliary density $n(\mathbf r)$, producing a shadow Born–Oppenheimer surface that is cheap to evaluate and can be integrated exactly within a reversible extended dynamics [1705.10845]. The auxiliary electronic field is not the physical wavefunction manifold; it is an extended variable whose dynamics is governed by a kernel-preconditioned harmonic oscillator centered on the variational minimum of the linearized functional.

In the adiabatic limit, the equations of motion become
$$
M_I {\ddot R}_I = \left.-\frac{\partial U^{(1)}_{\rm BO}[n]({\bf R})}{\partial R_I}\right|_{n},
$$
and
$$
{\ddot n}({\bf r}) = -\omega^2 \int K({\bf r},{\bf r'})\left(\rho_{\rm min}^{(1)}[n]({\bf r'})-n({\bf r'})\right)d{\bf r'}.
$$
The kernel is the inverse Jacobian of the residual map, written in coarse-grained notation as
$$
K = \left[\frac{\delta \rho[n]}{\delta n} - I\right]^{-1},
$$
and serves both as the inner-product metric of the harmonic extension and as a nonlocal preconditioner [1705.10845].

This branch is explicitly positioned as higher-order than CPMD. The paper states that CPMD has electronic error $\sim \delta t$ and potential-energy-surface error $\sim \delta t^2$, whereas XL-BOMD has electronic error $\sim \delta t^2$ and shadow-surface error $\sim \delta t^4$ [1705.10845]. It also claims timesteps of the same order as in direct BOMD, no requirement of iterative nonlinear electronic ground-state optimization prior to force evaluations, and no systematic drift in the total energy because the approximate forces derive from a genuine shadow Hamiltonian or free-energy functional [1705.10845]. At finite electronic temperature, the surface becomes
$$
\Omega^{(1)}_{\rm KS}[n]({\bf R}) = U^{(1)}_{\rm KS}[n]({\bf R}) - T_e S[{\bf f}],
$$
with
$$
S[{\bf f}] = -k_B \sum_i \left(f_i \ln f_i + (1-f_i)\ln(1-f_i)\right),
$$
so fractional occupations become intrinsic to stability and conservative finite-temperature forces [1705.10845].

A related but distinct line of work addresses one of first-generation CPMD’s core weaknesses—the hand-tuned fictitious mass—without adopting predictor–corrector propagation or a shadow BO surface. The adaptive-mass algorithm for CP and Ehrenfest dynamics introduces a time-dependent artificial mass chosen from local gap and velocity information, motivated by Landau–Zener transition probabilities and supported by an Ehrenfest gap-based asymptotic estimate [1409.4992]. This work belongs to the broader agenda of improving Car–Parrinello-like dynamics, but it is not a predictor–corrector second-generation CPMD formulation in the CP2K sense, nor an XL-BOMD reformulation [1409.4992].

Taken together, these strands define second-generation Car–Parrinello methods as a family of SCF-avoiding, Born–Oppenheimer-oriented AIMD schemes. In the narrower and most widely cited CP2K form, the family is characterized by density-matrix-guided ASPC prediction, one OT-based preconditioned correction step, approximate Harris–Foulkes energy and force evaluation, and a modified Langevin treatment of residual dissipation [2205.14741] [1303.2067] [1201.5945]. In the broader extended-Lagrangian interpretation, it denotes a higher-order, kernel-preconditioned successor to CPMD that replaces direct BOMD’s incomplete-SCF workflow by exact integration on a nearby shadow BO or free-energy surface [1705.10845].

Source: https://www.emergentmind.com/topics/second-generation-car-parrinello-methods