Grüneisen parameter

The Grüneisen parameter, γ, named after Eduard Grüneisen, describes the effect that changing the volume of a crystal lattice has on its vibrational properties, and, as a consequence, the effect that changing temperature has on the size or dynamics of the lattice. The term is usually reserved to describe the single thermodynamic property γ, which is a weighted average of the many separate parameters γ_i entering the original Grüneisen's formulation in terms of the phonon nonlinearities.^[1]

Thermodynamic definitions

Because of the equivalences between many properties and derivatives within thermodynamics (e.g. see Maxwell Relations), there are many formulations of the Grüneisen parameter which are equally valid, leading to numerous distinct yet correct interpretations of its meaning.

Some formulations for the Grüneisen parameter include:

$\gamma = V \left(\frac{dP}{dE}\right)_V = \frac{\alpha K_T}{C_V \rho} = \frac{\alpha K_S}{C_P \rho} = \frac{\alpha v_s^2}{C_P}$

where V is volume, $C_P$ and $C_V$ are the principal (i.e. per-mass) heat capacities at constant pressure and volume, E is energy, α is the volume thermal expansion coefficient, $K_S$ and $K_T$ are the adiabatic and isothermal bulk moduli, $v_s$ is the speed of sound in the medium, and ρ is density. The Grüneisen parameter is dimensionless.

Grüneisen constant for perfect crystals with pair interactions

The expression for the Grüneisen constant of a perfect crystal with pair interactions in $d$-dimensional space has the form:^[2]

$\Gamma_0 = -\frac{1}{2d}\frac{\Pi'''(a)a^2 + (d-1)\left[\Pi''(a)a - \Pi'(a)\right]}{\Pi''(a)a + (d-1)\Pi'(a)},$

where $\Pi$ is the interatomic potential, $a$ is the equilibrium distance, $d$ is the space dimensionality. Relations between the Grüneisen constant and parameters of Lennard-Jones, Morse, and Mie potentials are presented in the table below.

Lattice	Dimensionality	Lennard-Jones potential	Mie Potential	Morse potential
Chain	$d=1$	$10\frac{1}{2}$	$\frac{m+n+3}{2}$	$\frac{3\alpha a}{2}$
Triangual lattice	$d=2$	$5$	$\frac{m+n+2}{4}$	$\frac{3\alpha a - 1}{4}$
FCC, BCC	$d=3$	$\frac{19}{6}$	$\frac{n+m+1}{6}$	$\frac{3\alpha a-2}{6}$
"Hyperlattice"	$d=\infty$	$-\frac{1}{2}$	$-\frac{1}{2}$	$-\frac{1}{2}$
General formula	$d$	$\frac{11}{d}-\frac{1}{2}$	$\frac{m+n+4}{2d}-\frac{1}{2}$	$\frac{3\alpha a + 1}{2d}-\frac{1}{2}$

The expression for the Grüneisen constant of a 1D chain with Mie potential exactly coincides with the results of MacDonald and Roy.^[3] Using the relation between Grüneisen parameter and interatomic potential one can derive the simple necessary and sufficient condition for Negative Thermal Expansion in perfect crystals with pair interactions

$\Pi'''(a)a > -(d-1)\Pi''(a).$

Microscopic definition via the phonon frequencies

The physical meaning of the parameter can also be extended by combining thermodynamics with a reasonable microphysics model for the vibrating atoms within a crystal. When the restoring force acting on an atom displaced from its equilibrium position is linear in the atom's displacement, the frequencies ω_i of individual phonons do not depend on the volume of the crystal or on the presence of other phonons, and the thermal expansion (and thus γ) is zero. When the restoring force is non-linear in the displacement, the phonon frequencies ω_i change with the volume $V$. The Grüneisen parameter of an individual vibrational mode $i$ can then be defined as (the negative of) the logarithmic derivative of the corresponding frequency $\omega_i$:

$\gamma_i= - \frac{V}{\omega_i} \frac{\partial \omega_i}{\partial V}.$

Relationship between microscopic and thermodynamic models

Using the quasi-harmonic approximation for atomic vibrations, the macroscopic Grüneisen parameter (γ) can be related to the description of how the vibration frequencies (phonons) within a crystal are altered with changing volume (i.e. γ_i's). For example, one can show that

$\gamma = \frac{\alpha K_T}{C_V \rho}$

if one defines $\gamma$ as the weighted average

$\gamma = \frac{\sum_i \gamma_i c_{V,i} }{ \sum_i c_{V,i} },$

where $c_{V,i}$'s are the partial vibrational mode contributions to the heat capacity, such that $C_{V} = \frac{1}{\rho V} \sum_i c_{V,i} .$

Proof

To prove this relation, it is easiest to introduce the heat capacity per particle $\tilde{C}_V = \sum_i c_{V,i}$; so one can write

$\frac{\sum_i \gamma_i c_{V,i}}{\tilde{C}_V} = \frac{\alpha K_T}{C_V \rho} = \frac{\alpha V K_T}{\tilde{C}_V}$.

This way, it suffices to prove

$\sum_i \gamma_i c_{V,i} = \alpha V K_T$.

Left-hand side (def):

$\sum_i \gamma_i c_{V,i} = \sum_i \left[- \frac{V}{\omega_i} \frac{\partial \omega_i}{ \partial V} \right] \left[ k_B \left(\frac{\hbar \omega_i}{k_B T}\right)^2 \frac{\exp\left( \frac{\hbar \omega_i}{k_B T} \right)}{\left[\exp\left(\frac{\hbar \omega_i}{k_B T}\right) - 1\right]^2} \right]$

Right-hand side (def):

$\alpha V K_T = \left[ \frac{1}{V} \left(\frac{\partial V}{ \partial T}\right)_P \right] V \left[-V \left(\frac{\partial P}{\partial V}\right)_T\right] = - V \left( \frac{\partial V}{\partial T} \right)_P \left(\frac{\partial P}{\partial V}\right)_T$

Furthermore (Maxwell relations):

$\left( \frac{\partial V}{\partial T} \right)_P = \frac{\partial}{\partial T} \left(\frac{\partial G}{\partial P}\right)_T = \frac{\partial}{\partial P} \left(\frac{\partial G}{\partial T}\right)_P = - \left( \frac{\partial S}{\partial P} \right)_T$

Thus

$\alpha V K_T = V \left( \frac{\partial S}{\partial P} \right)_T \left(\frac{\partial P}{\partial V}\right)_T = V \left( \frac{\partial S}{\partial V} \right)_T$

This derivative is straightforward to determine in the quasi-harmonic approximation, as only the ω_i are V-dependent.

$\frac{\partial S}{\partial V} = \frac{\partial }{\partial V} \left\{ - \sum_i k_B \ln\left[ 1 - \exp\left( -\frac{\hbar\omega_i (V)}{k_BT} \right) \right] + \sum_i \frac{1}{T} \frac{\hbar\omega_i (V)}{\exp\left(\frac{\hbar\omega_i (V)}{k_BT}\right) - 1} \right\}$

$V \frac{\partial S}{\partial V} = - \sum_i \frac{V}{\omega_i} \frac{\partial \omega_i}{\partial V} \;\; k_B \left(\frac{\hbar \omega_i}{k_BT}\right)^2 \frac{\exp\left( \frac{\hbar \omega_i}{k_B T} \right)}{\left[\exp\left(\frac{\hbar \omega_i}{k_B T}\right) - 1\right]^2} = \sum_i \gamma_i c_{V,i}$

This yields

$\gamma = \dfrac{\sum_i \gamma_i c_{V,i}}{\sum_i c_{V,i}} = \dfrac{\alpha V K_T}{\tilde{C}_V}.$

See also

• Debye model
• Negative thermal expansion
• Mie–Gruneisen equation of state

External links

• Definition from Eric Weisstein's World of Physics

References

1. Grüneisen, E. (1912), "Theorie des festen Zustandes einatomiger Elemente", Annalen der Physik 344: 257–306, doi:10.1002/andp.19123441202
2. Krivtsov, A.M.; Kuzkin, V.A. (2011), "Derivation of Equations of State for Ideal Crystals of Simple Structure", Mechanics of Solids 46 (3): 387–399, doi:10.3103/S002565441103006X
3. MacDonald, D. K. C.; Roy, S.K. (1955), "Vibrational Anharmonicity and Lattice Thermal Properties. II", Phys. Rev. 97: 673–676, Bibcode:1955PhRv...97..673M, doi:10.1103/PhysRev.97.673