---
title: Temperature Dependent Effective Potential
url: https://www.emergentmind.com/topics/temperature-dependent-effective-potential-tdep
type: topic
---

# Temperature Dependent Effective Potential

The temperature dependent effective potential (TDEP) framework is a first-principles approach for capturing thermodynamic, vibrational, and transport properties of solids with strong anharmonicity by extracting effective force constants directly from finite-temperature atomistic trajectories. Rather than relying on the static, zero-temperature potential energy surface, TDEP replaces the bare Taylor expansion with a temperature-renormalized model Hamiltonian, whose parameters are fitted via statistical inference on molecular dynamics (MD) or quantum path-integral MD data. This variational, data-driven pipeline enables a self-consistent, quantitatively accurate description of phonons, thermodynamics, and elastic and transport coefficients across the full temperature range, including dynamically stabilized and quantum-fluctuation-dominated regimes. The method has been generalized to extract higher-order interatomic force constants, extended for quantum nuclei using path-integral sampling, and incorporated into multiple implementations for materials ranging from elemental solids to complex oxides and paramagnetic compounds [1303.1145, 1308.5436, 1604.08855, 2211.14474, 2409.09551].

## 1. Mathematical Formalism and Model Hamiltonian

TDEP begins from a general expansion of the Born–Oppenheimer or Born–von Kármán potential energy $U(\{ \mathbf{u} \})$ in atomic displacements about temperature-dependent equilibrium positions:
\[
U(\{ \mathbf{u} \}) = U_0 + \sum_i \Phi_i u_i + \frac{1}{2}\sum_{ij}\Phi_{ij}^{(2)}u_iu_j + \frac{1}{3!}\sum_{ijk}\Phi_{ijk}^{(3)}u_iu_ju_k + \dots
\]
where $u_i$ is the displacement of atom $i$, and $\Phi^{(n)}$ are the $n$th-order interatomic force constants (IFCs), with $U_0$ an energy reference [1303.1145, 2409.09551].

In TDEP, these IFCs are replaced by *effective*, temperature-dependent parameters $\Phi^{(n)}(T)$, defining an optimized model Hamiltonian
\[
\hat H_{TDEP}(T) = U_0 + \sum_i \frac{p_i^2}{2m_i} + \frac{1}{2}\sum_{ij}\Phi_{ij}^{(2)}(T)u_iu_j + \frac{1}{3!}\sum_{ijk}\Phi_{ijk}^{(3)}(T)u_iu_ju_k + \dots
\]
which is fit to force-displacement data sampled from MD or quantum MD at temperature $T$. By construction, $\Phi^{(2)}(T)$ incorporates anharmonic and quantum fluctuations renormalized at that temperature [1303.1145, 2211.14474, 2409.09551].

## 2. Fitting Procedure and Statistical Extraction from Simulation

To extract the $\Phi^{(n)}(T)$ effective force constants, TDEP employs a linear least-squares minimization of the difference between the true forces (from ab initio MD or PIMD) and the model forces, for a set of $N$ uncorrelated snapshots $\{ \mathbf{u}(t), \mathbf{F}^{MD}(t)\}$ [1303.1145, 1308.5436, 2409.09551]:
\[
\min_{\Phi^{(2)},\Phi^{(3)},\dots} \frac{1}{N} \sum_{t=1}^{N} \left\| \mathbf{F}^{MD}(t) - \mathbf{F}^{model}(t; \Phi^{(2)},\Phi^{(3)},\dots) \right\|^2
\]
where the model force for atom $i$ is
\[
F_i^{model}(t) = -\sum_j \Phi_{ij}^{(2)} u_j(t) - \frac{1}{2} \sum_{jk}\Phi_{ijk}^{(3)}u_j(t)u_k(t) - \cdots
\]
The covariance structure and force constants are reduced by imposing crystal translation, point group, and permutation symmetries, as well as acoustic sum rules. For higher-order IFCs, the least-squares system is enlarged to include cubic and (optionally) quartic terms [1308.5436, 2305.10918]. In quantum solids or at low temperature, path-integral MD sampling is essential to capture nuclear quantum effects without artifacts [2211.14474].

Variational, self-consistent algorithms (sometimes with re-weighting for all prior samples) are used to ensure convergence and statistical robustness [2110.04771]. Symmetry-imposed variants (e.g., SIFC-TDEP for paramagnetic or disordered systems) enforce ideal lattice point-group averaging for force constants, increasing efficiency and statistical convergence [1604.08855].

## 3. Phonons, Elastic Constants, and Free Energy from Effective Force Constants

The fitted effective harmonic force constants $\Phi_{ij}^{(2)}(T)$ define the temperature-dependent dynamical matrix:
\[
D_{ij}(\mathbf{q};T) = \frac{1}{\sqrt{m_i m_j}} \sum_{\mathbf{R}} \Phi^{(2)}_{ij}(\mathbf{R};T) e^{i\mathbf{q}\cdot\mathbf{R}}
\]
Diagonalization yields temperature-renormalized phonon frequencies $\omega_{\mathbf{q},\nu}(T)$ and eigenvectors. The Helmholtz free energy is then constructed from harmonic and anharmonic (residual) terms:
\[
F_{TDEP}(T, V) = U_0(T,V) + F_{vib}[g(\omega;T)]
\]
\[
F_{vib} = \int_0^\infty g(\omega)\left[ k_BT\ln(1-e^{-\hbar\omega/k_BT}) + \frac{1}{2}\hbar\omega \right] d\omega
\]
with $g(\omega)$ the phonon density of states. Optionally, further anharmonic corrections are computed via direct averaging or perturbative expansions in higher-order IFCs [1303.1145, 2110.04771].

The elastic constants at finite temperature are extracted using the long-wavelength limit (Born–Huang relations) of the dynamical matrix or via direct sums over IFCs:
\[
C_{ijkl}(T) = -\frac{1}{2V} \sum_{n} \Phi_{0n}^{ij}(T) r_n^k r_n^l
\]
Voigt–Reuss–Hill averaging provides polycrystalline moduli [1604.08855, 2409.09551].

## 4. Incorporation of Anharmonicity, Quantum Effects, and Sampling Schemes

TDEP, by design, captures the leading effect of strong anharmonicity via the thermal renormalization of force constants. Extensions have been developed to explicitly fit third- and fourth-order IFCs, enabling direct calculation of Grüneisen parameters, three- and four-phonon scattering rates, and phonon linewidths/lifetimes relevant for transport [1308.5436, 2305.10918, 2409.09551]. The quantum-classical boundary is addressed by full temperature-dependent potential (FTDP) methods, which use ab initio path-integral MD sampling to extract $\Phi^{(n)}(T)$ unbiased with respect to the quantum distribution, ensuring physicality below the Debye temperature and for light-element compounds [2211.14474, 2409.09551].

The workflow for quantum TDEP (FTDP) is as follows:
1. Generate PIMD trajectories with a machine-learned or ab initio potential.
2. For each snapshot and “bead”, record atomic displacements and forces (excluding spring terms).
3. Fit the force constants as in the classical case, discarding bead–bead spring contributions [2211.14474, 2409.09551].
4. Check convergence with supercell and bead number; statistical averaging and error bars are recommended.

## 5. Applications and Benchmark Results

TDEP and its variants have enabled quantitatively accurate predictions for multiple classes of materials:
- **Anharmonic metals and dynamically unstable phases**: Dynamical stabilization of bcc Zr, which is harmonically unstable at $T=0$ but stabilized by TDEP-renormalized force constants at $T > 1000$ K, producing real phonon frequencies, correct free energy surfaces, and finite-$T$ EOS [1303.1145].
- **Paramagnetic and disordered magnets**: Assessment of elastic constants in high-$T$ paramagnetic CrN via SIFC-TDEP + DLM-MD, revealing a 14% reduction in Young’s modulus from 300 to 1200 K and increasing isotropy with $T$ [1604.08855].
- **Quantum and ultra-anhamonic solids**: In metallic hydrogen, FTDP enabled detection of colossal nuclear quantum effects and the breakdown of the quasi-phonon picture [2211.14474].
- **Thermal transport**: Calculation of lattice thermal conductivity in ceramics (e.g., CeO$_2$, MgO) requires temperature-renormalized phonons and explicit inclusion of both three- and four-phonon processes, with TDEP dramatically improving agreement with experiment [2305.10918, 2409.09551].
- **Critical and soft-mode systems**: Modified TDEP formalism remedies unphysical limitations (e.g., absence of diverging susceptibility at second-order transitions), enabling correct dynamical and static response functions near ferroelectric or structural instabilities [2407.03090].

A summary table of selected applications is provided:

| System               | Extension/Variant    | Remarkable Feature/Result                                                              |
|----------------------|---------------------|----------------------------------------------------------------------------------------|
| bcc Zr               | TDEP, classical     | Dynamical stabilization at high T, force localization with T [1303.1145]               |
| CrN paramagnet       | SIFC-TDEP + DLM-MD  | Single/polycrystal $C_{ij}(T)$, enhanced isotropy, ~14% $E$ reduction [1604.08855]     |
| Metallic H           | FTDP (quantum)      | Colossal NQE, phonon nonconvergence, lattice quantum fluctuations [2211.14474]         |
| Si, FeSi             | TDEP + cubic        | Mode Grüneisen parameters, $T$-dependent lifetimes [1308.5436]                         |
| CeO$_2$, MgO         | TDEP + quartic      | Four-phonon rates, $T$-renormalized $\kappa(T)$ matching experiment [2305.10918]       |
| PbTe, CsSnI$_3$      | Modified TDEP       | Critical-mode stabilization, correct spectral and static response [2407.03090]         |

## 6. Advantages, Limitations, and Recent Developments

TDEP provides direct inclusion of anharmonic effects, captures thermodynamic stabilization of otherwise unstable phases, and is extensible to arbitrary order with symmetry-imposed reduction in parameter space [1303.1145, 1604.08855, 2409.09551]. The method is robust with respect to statistical noise (by leveraging symmetry and block-averaging), highly parallelizable, and compatible with machine-learned force fields [2409.09551].

Limitations include: (i) standard TDEP models are quadratic (or, with extensions, cubic/quartic) and so only implicitly capture higher-order correlations unless higher IFCs are explicitly fit [1308.5436, 2305.10918]; (ii) classical TDEP is invalid below the Debye temperature without quantum sampling; (iii) conventional TDEP does not produce divergent susceptibilities at critical points—a new correction (subtracting the static bubble from the self-energy) is required to restore physicality in the dynamical spectra [2407.03090].

Recent work incorporates PGIMD for quantum sampling, machine-learning potentials for efficiency in large systems, self-consistent ensemble reweighting for improved convergence, and hybrid approaches to join TDEP with self-consistent phonon or renormalization group frameworks [2110.04771, 2211.14474, 2409.09551, 2407.03090].

## 7. Extensions and Variants: SIFC-TDEP, FTDP, and Beyond

- **SIFC-TDEP**: Enforces full point-group symmetry of the underlying lattice for efficient force-constant reduction and faster convergence, especially in large/disordered or magnetic supercells [1604.08855].
- **FTDP (Full Temperature-Dependent Potential)**: Generalizes TDEP by fitting to path-integral ab initio trajectories, ensuring validity in the quantum limit and for systems with significant nuclear quantum effects [2211.14474].
- **Step-by-step higher-order fitting**: Successively fits cubic and quartic terms to force residuals for improved accuracy in scattering and thermal transport properties [1308.5436, 2305.10918].
- **Modified self-energies for criticality**: Corrects the static limit of bubble diagrams to ensure divergence of susceptibility at phase transitions and proper frequency dependences of dynamical response [2407.03090].
- **Integration with machine-learned potentials and compressive sensing**: Enables tractable extraction of effective force constants from long, large-scale quantum MD trajectories [2409.09551].

TDEP thus remains a leading approach for first-principles investigation of strongly anharmonic, thermally activated, and quantum-fluctuating solids across the material spectrum. Its implementation, validation, and further development continue to be active research areas [1303.1145, 1308.5436, 1604.08855, 2110.04771, 2211.14474, 2305.10918, 2407.03090, 2409.09551].

Source: https://www.emergentmind.com/topics/temperature-dependent-effective-potential-tdep