Universally Divergent Grüneisen Ratio

• The universally divergent Grüneisen ratio is a thermodynamic response function that quantifies the critical scaling of entropy near classical and quantum critical points.
• It displays universal divergence forms across various tuning parameters—such as pressure, magnetic field, and strain—indicating its independence from microscopic details.
• Its scaling behavior aids in diagnosing quantum phase transitions and enables practical applications like quantum magnetocaloric cooling and material classification.

The universally divergent Grüneisen ratio is a fundamental thermodynamic response function that encodes the critical scaling of entropy with respect to tuning parameters near both classical and quantum critical points (QCPs). Universality is reflected in both the divergence exponents and the scaling forms, independent of microscopic system details, and extends to various control fields such as pressure, magnetic field, uniaxial strain, and—in anisotropic systems—field orientation (the "rotational Grüneisen ratio"). The emergence of universal Grüneisen divergences provides a direct and quantitative thermodynamic fingerprint of criticality and is central to diagnosing, classifying, and exploiting quantum phase transitions in strongly correlated systems, heavy fermion metals, quantum magnets, and engineered quantum materials.

1. Thermodynamic Definition and Physical Origin

The Grüneisen ratio $\Gamma$ is defined as the ratio of a generalized thermal expansion coefficient $\alpha_g$ to a specific heat $c_g$, both taken at fixed external tuning parameter $g$: $$\Gamma_g = \frac{\alpha_g}{c_g} = -\frac{1}{T} \frac{\left(\partial S/\partial g\right)_T}{\left(\partial S/\partial T\right)_g} = \frac{1}{T} \left(\frac{\partial T}{\partial g}\right)_S.$$ Here, $S$ is the system entropy, and $g$ may refer to pressure, magnetic field, strain, or other external controls. This definition encodes the adiabatic change in temperature under variation of $g$—the caloric effect associated with the corresponding control field. Near a critical point, the accumulation of entropy and softening of energy scales drives both $\alpha_g$ and $c_g$ to display singular behavior, but with distinct exponents, leading to a divergence of $\alpha_g$0 that is determined by the universality class rather than microscopics (Zhou et al., 22 Sep 2025).

2. Universal Critical Scaling and Divergence Shapes

The critical behavior of $\alpha_g$1 is governed by scaling theory. Near a QCP tuned by $\alpha_g$2 and at temperature $\alpha_g$3, the singular part of the free energy density has the scaling form: $\alpha_g$4 where $\alpha_g$5 is the spatial dimension, $\alpha_g$6 the dynamical exponent, and $\alpha_g$7 the correlation length exponent. Direct computation yields for the Grüneisen ratio (Squillante, 2024, Yu et al., 2019, Zhou et al., 22 Sep 2025): $\alpha_g$8 or, equivalently,

$\alpha_g$9

where $c_g$0 and $c_g$1 are universal scaling functions.

The key universal results are:

• Quantum critical regime ($c_g$2, $c_g$3):

$c_g$4

• Quantum disordered/ordered wings ($c_g$5, $c_g$6):

$c_g$7

The exponent $c_g$8 encodes the universality class of the transition; all systems with the same $c_g$9 and $g$0 exhibit the same divergence (Squillante, 2024, Yu et al., 2019, Zhou et al., 22 Sep 2025, Gegenwart, 2016).

3. Types of Grüneisen Ratios and Control Parameters

The Grüneisen ratio formalism generalizes to various tuning fields, each yielding a distinct but related divergent response:

• Volume (pressure-tuned) Grüneisen ratio: $g$1. Sensitive to pressure-driven criticality, with $g$2 the volumetric thermal expansion and $g$3 the specific heat (Gegenwart, 2016, Squillante, 2024).
• Magnetic Grüneisen ratio: $g$4. Diverges at field-tuned QCPs (Gegenwart, 2016, Liu et al., 12 Jan 2026).
• Interaction Grüneisen ratio: In ultracold gases with tunable interaction strength $g$5, $g$6 (Yu et al., 2019).
• Strain/phonon-mode Grüneisen ratio: Probes the softening of a collective mode, e.g., $g$7 and diverges as the soft mode energy vanishes at a structural QCP (Franklin et al., 2020).
• Rotational Grüneisen ratio: When criticality is tuned by field orientation $g$8 in an anisotropic system, $g$9. Displays universal divergence analogous to the standard Grüneisen response and directly probes anisotropic quantum criticality (Yuasa et al., 2024).

All these ratios are constrained by universal identities in integrable models (e.g., $$\Gamma_g = \frac{\alpha_g}{c_g} = -\frac{1}{T} \frac{\left(\partial S/\partial g\right)_T}{\left(\partial S/\partial T\right)_g} = \frac{1}{T} \left(\frac{\partial T}{\partial g}\right)_S.$$0 in $$\Gamma_g = \frac{\alpha_g}{c_g} = -\frac{1}{T} \frac{\left(\partial S/\partial g\right)_T}{\left(\partial S/\partial T\right)_g} = \frac{1}{T} \left(\frac{\partial T}{\partial g}\right)_S.$$1 dimensions) (Yu et al., 2019).

4. Experimental Realizations and Signatures

Experiments confirm universal Grüneisen divergences in a diverse set of systems:

Material/System	Tuning Variable	Universality Class	Divergence	Reference
CeRhSn, CeIrSn	Field orientation	Anisotropic heavy fermion	$$\Gamma_g = \frac{\alpha_g}{c_g} = -\frac{1}{T} \frac{\left(\partial S/\partial g\right)_T}{\left(\partial S/\partial T\right)_g} = \frac{1}{T} \left(\frac{\partial T}{\partial g}\right)_S.$$2, critical angle $$\Gamma_g = \frac{\alpha_g}{c_g} = -\frac{1}{T} \frac{\left(\partial S/\partial g\right)_T}{\left(\partial S/\partial T\right)_g} = \frac{1}{T} \left(\frac{\partial T}{\partial g}\right)_S.$$3, $$\Gamma_g = \frac{\alpha_g}{c_g} = -\frac{1}{T} \frac{\left(\partial S/\partial g\right)_T}{\left(\partial S/\partial T\right)_g} = \frac{1}{T} \left(\frac{\partial T}{\partial g}\right)_S.$$4	(Yuasa et al., 2024)
Nd$$\Gamma_g = \frac{\alpha_g}{c_g} = -\frac{1}{T} \frac{\left(\partial S/\partial g\right)_T}{\left(\partial S/\partial T\right)_g} = \frac{1}{T} \left(\frac{\partial T}{\partial g}\right)_S.$$5BWO$$\Gamma_g = \frac{\alpha_g}{c_g} = -\frac{1}{T} \frac{\left(\partial S/\partial g\right)_T}{\left(\partial S/\partial T\right)_g} = \frac{1}{T} \left(\frac{\partial T}{\partial g}\right)_S.$$6	Magnetic field	3D Ising (CEP)	$$\Gamma_g = \frac{\alpha_g}{c_g} = -\frac{1}{T} \frac{\left(\partial S/\partial g\right)_T}{\left(\partial S/\partial T\right)_g} = \frac{1}{T} \left(\frac{\partial T}{\partial g}\right)_S.$$7 with $$\Gamma_g = \frac{\alpha_g}{c_g} = -\frac{1}{T} \frac{\left(\partial S/\partial g\right)_T}{\left(\partial S/\partial T\right)_g} = \frac{1}{T} \left(\frac{\partial T}{\partial g}\right)_S.$$8	(Liu et al., 12 Jan 2026)
SrTiO$$\Gamma_g = \frac{\alpha_g}{c_g} = -\frac{1}{T} \frac{\left(\partial S/\partial g\right)_T}{\left(\partial S/\partial T\right)_g} = \frac{1}{T} \left(\frac{\partial T}{\partial g}\right)_S.$$9 (strained)	Uniaxial strain	Soft-phonon QCP	$S$0	(Franklin et al., 2020)
1D/2D quantum magnets	Field, pressure	Ising, Potts, Heisenberg, O(3)	Systematic scaling function collapse, exponents $S$1 determined numerically	(Zhou et al., 22 Sep 2025)
Molecular Fabre salts, V$S$2O$S$3, $S$4He, ultracold atoms	$S$5, $S$6, $S$7	Classical, Mott/BEC, quantum gases	Universally $S$8	(Squillante, 2024, Souza et al., 2016, Yu et al., 2019)

Measurement protocols typically involve either monitoring temperature change under adiabatic parameter sweeps (magnetocaloric, barocaloric, mechanocaloric effects) or direct extraction of the thermal expansion and specific heat under near-critical conditions. Data collapse analyses on universal scaling functions confirm the theoretical predictions for systems across dimensions and field types (Zhou et al., 22 Sep 2025, Yuasa et al., 2024, Squillante, 2024).

5. Breakdown, Limitations, and Generalizations

The universal divergence assumes well-defined second-order transitions, sufficiently weak coupling to secondary degrees of freedom (e.g., modest magnetoelastic coupling), and the validity of scaling/hyperscaling. Several factors may evade or regularize the divergence:

• Finite size and disorder: Rounds the divergence as the correlation length saturates.
• First-order or preempted transitions: No true criticality, possible cutoff of divergence.
• Background contributions and crossovers: Non-singular heat capacity or thermal expansion terms can mask or reduce the observed exponent over limited windows (Squillante, 2024, Gegenwart, 2016).
• Non-extensive entropy and $S$9-generalizations: In the presence of long-range correlations that invalidate extensive Boltzmann-Gibbs entropy, generalizing to non-additive $g$0-entropy leads to non-diverging Grüneisen ratios at the unique $g$1 for which entropy regains extensivity. The apparent divergence in conventional statistics is then an artifact of improper entropy assignment (Soares et al., 2024). However, under standard (additive) entropy, the observed divergence remains robust.

6. Extensions: Quantum Information, Entanglement, and Rotational Criticality

At strictly zero temperature, the classical thermodynamic Grüneisen ratio is undefined, but an exact quantum analogue may be formulated in terms of ground-state energy derivatives, and, via the Hellmann–Feynman relation, derivatives of entanglement measures such as the von Neumann entropy. If the ground-state energy $g$2 is nonlinear in the tuning parameter at criticality, the second derivative diverges, and the quantum Grüneisen ratio $g$3 inherits the universal divergence: $g$4 This divergence is strongly linked to breakdowns of the Hellmann–Feynman theorem at QCPs and encodes universal sensitivity of ground-state entanglement to external control (Squillante et al., 2023).

For highly anisotropic systems, orientational tuning by field rotation (rotational Grüneisen ratio $g$5) generates an entire line of QCPs, with divergence controlled by the easy-axis field component. Data collapse in heavy fermion materials such as CeRhSn and CeIrSn confirms the scaling predictions for $g$6 and its role as a universal probe of strong Ising anisotropy (Yuasa et al., 2024).

7. Impact and Applications

The universality of the Grüneisen divergence supplies a unifying, experimentally accessible handle on quantum criticality across a broad spectrum of physical platforms:

• Quantum magnetocaloric refrigeration: Exploits the enhanced $g$7 near critical fields for efficient sub-Kelvin cooling (Liu et al., 12 Jan 2026, Zhou et al., 22 Sep 2025).
• Quantum materials diagnostics: Directly identifies which collective mode drives criticality (e.g., the transverse soft-mode in SrTiO$g$8 (Franklin et al., 2020)).
• Universality class extraction: Enables determination of $g$9 through scaling of $g$0 versus $g$1 and control parameters (Squillante, 2024, Gegenwart, 2016).
• Criticality in ultracold gases and unconventional systems: Bethe-Ansatz exact scaling in 1D gases (Yu et al., 2019); criticality in quantum entanglement landscapes (Squillante et al., 2023).
• Development of scaling maps: Systematic mapping of quantum critical regimes and crossovers via Grüneisen amplitude and collapse enables robust material comparisons and confirmation of universal thermodynamics (Zhou et al., 22 Sep 2025).

The universally divergent Grüneisen ratio is thus both a quantitative diagnostic of critical phenomena and a tool for exploring and utilizing criticality-driven phenomena in quantum materials.