Quasi-Harmonic Approximation¶
The quasi-harmonic approximation computes crystal free energy from harmonic phonons whose frequencies depend on a global static constraint, commonly volume.
Core Idea¶
The quasi-harmonic approximation (QHA) treats lattice vibrations as independent harmonic phonons at each chosen value of a global static constraint \(X\), while allowing their frequencies \(\omega_{\mathbf q\nu}(X)\) to change with that constraint. Volume is the usual choice, but \(X\) may also describe anisotropic strain, an external field, or internal cell distortion. The free energy combines static lattice energy, zero-point vibration, and thermal occupation of these constraint-indexed harmonic modes. In the common volume case, minimizing the appropriate free energy at stated temperature and pressure yields equilibrium volume and thermal expansion; comparing phase free energies can yield a phase ranking. The implicit anharmonic response comes through changing \(X\), not explicit phonon–phonon interactions at fixed \(X\).[1]
Structural Signature¶
Sig role-phrases: constraint-indexed lattice; harmonic phonons at each constraint; vibrational free energy; thermodynamic comparison; omitted intrinsic anharmonicity.
- Constraint-indexed lattice: a family of static crystal configurations supplies \(X\), usually volume \(V\) but not necessarily only volume.
- Harmonic phonons at each constraint: the force-constant calculation is repeated as \(X\) changes; its mode frequencies are \(\omega_{\mathbf q\nu}(X)\).
- Vibrational free energy: \(F(X,T)=U_0(X)+\tfrac12\sum_{\mathbf q\nu}\hbar\omega_{\mathbf q\nu}(X)+k_BT\sum_{\mathbf q\nu}\log[1-e^{-\hbar\omega_{\mathbf q\nu}(X)/(k_BT)}]\). The three terms are static, zero-point, and thermal, respectively.
- Thermodynamic comparison: the relevant potential is minimized over allowed \(X\); for \(X=V\) at fixed \(T,p\), this is \(F(V,T)+pV\), or analogous phase potentials are compared.
- Omitted intrinsic anharmonicity: explicit mode interactions at fixed \(X\) do not enter this harmonic-spectrum ansatz.[1]
What It Is Not¶
QHA is not a calculation at one fixed constraint: without a family of \(X\)-dependent spectra it cannot supply the constraint response that defines the method. A single-volume harmonic calculation therefore cannot by itself predict QHA thermal expansion, but absence of a volume grid does not exclude a QHA study parameterized by strain or another global constraint. QHA is not a fully anharmonic treatment of phonon lifetimes or strongly renormalized modes. It also does not ensure that every predicted phase boundary is accurate; small free-energy differences depend on electronic-structure and vibrational approximations.[1][2]
Scope of Application¶
First-principles work on silver calculated thermal properties from phonon dispersions within QHA, using volume dependence to estimate expansion. A separate study of MgSiO3 perovskite and post-perovskite combined static energies with quasiharmonic vibrational free energies to compare phase stability. These are both \(X=V\) applications; they do not exhaust the more general static-constraint formulation. They differ in whether the downstream target is one material's expansion or competition between phases.[3][2][1]
Clarity¶
The word “harmonic” describes each fixed-\(X\) calculation; “quasi” signals that changing \(X\) changes the harmonic force constants and phonon frequencies. Temperature changes phonon occupations and hence vibrational free energy, which can favor a different \(X\). Frequencies are recalculated from each constrained static configuration; they are not arbitrary temperature-dependent interaction shifts applied to one fixed structure.[1]
Manages Complexity¶
QHA replaces a fully interacting finite-temperature lattice problem with tractable harmonic calculations across selected values of \(X\). The resulting \(F(X,T)\) compresses many modes into a thermodynamic surface. A volume-based implementation still needs enough volume points, stable phonon spectra, and a reliable static-energy method; interpolation and numerical differences can matter when phases are close.[1][2]
Abstract Reasoning¶
Choose the physically relevant global static constraint, compute \(U_0(X)\) and a harmonic phonon spectrum at each sampled \(X\), then add the zero-point and thermal terms to build \(F(X,T)\). Minimize the potential appropriate to the controlled conditions. For \(X=V\), expansion follows from how the vibrational term shifts the minimum of \(F+pV\), not from motion in one fixed harmonic well. If a sampled mode becomes unstable or fixed-\(X\) interactions matter, the harmonic free-energy expression loses reliability there; a more explicit anharmonic treatment or stability analysis is needed before trusting that region.[1]
Knowledge Transfer¶
The same volume-parameterized calculation transfers from elemental silver to mineral phases because both supply static energies and phonon spectra across volumes. The general QHA structure also permits another global static \(X\), provided the harmonic spectrum is recomputed across its values and used in the free-energy comparison. Material-specific frequencies, pressure range, phase candidates, and the work term conjugate to a different constraint do not transfer automatically. A broad principle of controlled approximation appears elsewhere, but the named QHA mechanism requires lattice phonons and constraint-dependent vibrational thermodynamics.
Examples¶
Silver thermal expansion¶
Xie and colleagues calculated silver phonon dispersions and used QHA to obtain thermal properties. For this \(X=V\) application, the static input is the calculated silver lattice energy at sampled volumes; the vibrational input is the harmonic dispersion recomputed at those volumes. The zero-point and thermal mode terms turn those inputs into \(F(V,T)\), whose pressure-adjusted minimum gives an equilibrium volume and thus expansion as temperature changes.[3][1]
Mapped back: silver volumes supply \(X\) and \(U_0(V)\); harmonic dispersions supply \(\omega_{\mathbf q\nu}(V)\); their zero-point/thermal free energy shifts the \(F(V,T)+pV\) minimum; explicit phonon scattering at fixed volume is outside this calculation.
MgSiO3 phase stability¶
The perovskite/post-perovskite study combined static energies for competing MgSiO3 structures with QHA vibrational free energies. For each phase, volume is the chosen \(X\): a constraint-indexed harmonic spectrum supplies zero-point and thermal terms, and pressure-adjusted free-energy surfaces can be compared at stated \(T,p\). The downstream result is a phase-stability boundary rather than expansion of one phase, subject to small-energy and anharmonicity limits.[2][1]
Mapped back: each phase has its own \(U_0(V)\) and \(\omega_{\mathbf q\nu}(V)\); zero-point and thermal terms give competing free energies; phase ranking follows from the appropriate pressure-dependent potential, not from a frequency shift assigned at one fixed volume.
Structural Tensions¶
T1: Computational economy versus intrinsic anharmonic fidelity. Keeping modes harmonic at each \(X\) makes a broad temperature/constraint scan tractable and retains zero-point effects, but can miss strong fixed-\(X\) interactions and fail near an unstable soft mode. Adding explicit anharmonic dynamics can test that regime, but costs far more sampling and may obscure simple mode-by-mode attribution. If all relevant modes remain stable and comparison with measured or higher-level properties is good, QHA may be an adequate bounded model; if a mode is unstable, its harmonic free-energy contribution is not a trustworthy shortcut. Diagnostic: Are the modes stable across the selected \(X\) range, and is the omitted fixed-\(X\) coupling likely to change the target quantity?
T2: Decisive phase ranking versus uncertainty discipline. Taking the lowest computed free energy yields a usable phase boundary, but a tiny phase gap can reverse with static-energy, phonon, or interpolation errors. Demanding convergence and cross-method checks delays a single boundary estimate and costs computation, yet prevents false precision where phases are near degenerate. A robust gap supports ranking under the stated method; an unresolved gap calls for an uncertainty band or additional calculations, not an unqualified boundary. Diagnostic: Is the phase free-energy difference larger than plausible numerical and model uncertainty at the conditions of interest?
Structural–Framed Character¶
QHA is strongly structural within crystal thermodynamics. Its steps and neglected interactions are explicit; evaluative weight enters when judging whether its regime is adequate. The method is a scientific convention, not an institutional category, and depends on human choices of static constraint, sampling, and computational model. Its phonon and lattice vocabulary travels across materials and from volume to other global constraints, but applying “quasi-harmonic” to any mildly nonlinear system would require a separate mathematical definition. Recognition rests on constraint-indexed harmonic spectra and free energy, not an analogy to generic approximation. Its character: a formal physical approximation with constraint- and material-dependent validity.
Structural Core vs. Domain Accent¶
The skeletal relation is solving a parameterized family of locally simple models and letting variation of the parameter supply a thermodynamic response. QHA's domain-bound mechanism is harmonic crystal phonons recomputed at each global static constraint, then assembled into a temperature-dependent free-energy surface. The named entry fails the prime bar because its defining tests require phonons, a physically specified static constraint, and equilibrium thermodynamics; replacing them with arbitrary simple models leaves an analogy, not QHA. Whether the parameterized-family skeleton deserves a future prime is a separate question; no live strict parent is asserted here.
Instantiates / Related Primes¶
No strict parent edge is asserted in the unparented root. Live Debye Model fixes a different spectral simplification, Thermal Expansion is a possible result, and the live prime Approximation requires controlled error that QHA does not necessarily quantify. A harmonic-lattice-dynamics or constraint-indexed phonon-model intermediate remains an open taxonomy question.
Neighborhood in Abstraction Space¶
Quasi-Harmonic Approximation sits in a sparse region of the domain-specific corpus (82nd percentile for distinctiveness): few abstractions share its structure, so a faithful description tends to retrieve it precisely.
Family — Unclustered & Miscellaneous (2551 abstractions)
Nearest neighbors
- Davydov Soliton — 0.83
- Crystal momentum — 0.82
- Crystal Lattice — 0.82
- Electron backscatter diffraction — 0.81
- Rigidity Theory (Physics) — 0.81
Computed from structural-signature embeddings · 2026-10-08
Not to Be Confused With¶
Harmonic approximation: fixed-volume quadratic vibration model. Anharmonic phonon theory: includes interactions and temperature-dependent renormalization beyond ordinary QHA. Debye Model: an idealized phonon spectrum that may be combined with thermodynamics but does not define QHA. Thermal expansion: a physical response, not the approximation used to predict it.[1]
References¶
[1] Stefano Baroni, Paolo Giannozzi, and Eyvaz Isaev, “Thermal properties of materials from ab-initio quasi-harmonic phonons”, original investigator review, §2 Eq. (4) and discussion following it. registry ↩a ↩b ↩c ↩d ↩e ↩f ↩g ↩h ↩i ↩j
[2] “Equations of state and stability of MgSiO3 perovskite and post-perovskite phases from quantum Monte Carlo simulations”, original QHA phase-comparison study. registry ↩a ↩b ↩c ↩d
[3] Jianjun Xie et al., “First-principles calculation of the thermal properties of silver”, original study. registry ↩a ↩b