Temperature-Dependent Vibrational Spectra

• Temperature-dependent vibrational spectra are sets of vibrational mode frequencies, intensities, broadenings, and coupling effects that vary with temperature due to thermal expansion and anharmonicity.
• Experimental observations reveal systematic red/blue shifts, linewidth changes, and mode coupling signatures across Raman, IR, and EELS techniques, enabling extraction of key parameters.
• Computational methods such as QHA, AIMD, and self-consistent phonon theory provide practical approaches to model and interpret the complex temperature effects in vibrational spectra.

Temperature-dependent vibrational spectra are the set of vibrational mode frequencies, intensities, broadenings, and couplings in molecules, clusters, or extended solids as functions of temperature. The temperature dependence arises from several mechanisms: thermal expansion (volume or bond-length increase), occupancy of vibrational excited states, anharmonicity in the interatomic potential, and inter-mode (vibron–vibron or phonon–phonon) coupling. These effects manifest in systematic shifts (red or blue), mode broadenings, intensity changes, and mode coupling signatures across Raman, infrared (IR), inelastic neutron/X-ray scattering, and ultrafast spectroscopies.

1. Origin and Theoretical Frameworks

The basic theoretical description starts from the harmonic approximation, in which vibrational modes (phonons in extended systems, vibrons in clusters/molecules) are derived as eigenmodes of the mass-weighted dynamical matrix. This framework defines normal modes with well-defined zero-temperature frequencies and displacement patterns. Temperature dependence enters at several levels:

• Thermal expansion: The change in equilibrium geometry (cell parameters or bond lengths) with temperature modifies the interatomic force constants and thus the mode frequencies.
• Anharmonicity: Higher-order terms in the potential energy surface introduce shifts, broadenings, and multiphonon processes not captured by the harmonic model.
• Mode coupling: Different vibrational modes exchange energy via anharmonic terms, leading to finite lifetimes and complex energy flow within the spectrum.
• Statistical occupation: At finite T, vibrational states are thermally populated according to the Boltzmann distribution, enabling hot-band transitions and altering lineshapes.

The Grüneisen parameter $\gamma_i$ provides a linearized description of frequency shifts due to bond length change for mode $i$,

$$\gamma_i = -\frac{V}{\omega_i}\left( \frac{\partial \omega_i}{\partial V} \right)_T,$$

leading to a first-order frequency shift under isotropic expansion:

$$\Delta\omega_i(T) \simeq -\gamma_i \, \omega_i^0 \frac{\Delta L(T)}{L_0}$$

where $\omega_i^0$ is the zero-temperature frequency and $\Delta L(T)/L_0$ is the fractional bond expansion (Han et al., 2013).

2. Temperature-induced Frequency Shifts and Broadening

The temperature dependence of vibrational frequencies is both mode- and material-specific, governed by the competition between thermal expansion and intrinsic anharmonicity.

Representative examples:

• In Si$_{10}$H$_{16}$ clusters, TA-like Si–Si modes blue-shift with temperature ($\gamma<0$), while optical Si–Si and all Si–H modes red-shift ($\gamma>0$) (Han et al., 2013). Mode-specific Grüneisen parameters can be extracted from the shifts and bond-length expansions; for example, $\gamma_\text{TA} \approx -1.8$, $\gamma_\text{optical} \approx 1.2$.
• In h-BN, momentum-resolved electron energy loss spectroscopy (EELS) reveals softening (negative $\Delta\omega$) of TO and LO modes with increasing T, accelerating at high T due to strong cubic/quartic anharmonicity; DFT-based self-consistent phonon (SCP) theory with explicit thermal-expansion and multi-phonon corrections quantitatively reproduces these shifts (O'Hara et al., 2023).

In molecular crystals such as aspirin and paracetamol, low-frequency collective modes ($<300$ cm$^{-1}$) are most sensitive to lattice expansion and anharmonicity, shifting by up to 10–20 cm$^{-1}$ over 0 K–300 K, while high-frequency intra-molecular modes are comparatively rigid (1901.10587).

The temperature-dependence of linewidths arises primarily from anharmonic phonon-phonon scattering. The linewidth $\Gamma(T)$ typically follows

$$\Gamma(T) = \Gamma_0 + B\left[1 + 2n\left(\hbar\omega_0/2k_BT\right)\right]$$

with Bose population factor $n(x)$ and anharmonic coefficient $B$ (Mal et al., 2018). Broadening can be extremely strong in dynamically disordered lattices, e.g., methylammonium lead iodide shows cage-mode linewidths increasing from 10–20 cm$^{-1}$ (low $T$) to $>$80 cm$^{-1}$ (300 K), reflecting rapid phonon decay and motional disorder (Leguy et al., 2016).

3. Mechanisms: Quasi-harmonic, Anharmonic, and Coupled Mode Models

Quasi-harmonic approximation (QHA): Expands the Helmholtz free energy in cell parameters, allowing extraction of temperature-dependent lattice geometry and thus “quasi-harmonic” frequencies. QHA captures most frequency shifts due to expansion, accurate in the low-$T$ regime and for modes weakly affected by anharmonicity (1901.10587).

Anharmonic corrections: Enter via explicit treatment in ab initio molecular dynamics (AIMD), time-correlation formalism, self-consistent phonon theory, or perturbative phonon self-energies:

$$\omega(T) = \omega_0 + \Delta\omega_{\text{anh}}(T).$$

Three- and four-phonon interactions contribute to both frequency shifts and lifetime shortening, with quartic terms often dominating for optical modes at high T (O'Hara et al., 2023).

Mode coupling extraction: Projection of AIMD trajectories onto normal modes permits calculation of mode occupation autocorrelation functions $C_\nu(t)$; their Fourier transforms reveal Rabi-like oscillations corresponding to dominant vibron–vibron or phonon–phonon couplings ($g_{\nu j}$), which are themselves weakly temperature-dependent (increase $\sim$5% from 50 K to 150 K for Si$_{10}$H$_{16}$, with maximum $g_{\nu j} \simeq 2.5$ THz) (Han et al., 2013).

Moment-based methods: Classical spectral moments computed from static ensemble averages yield effective harmonic frequencies and line widths without explicit dynamics, enabling efficient, non-dynamical computation of temperature-dependent spectra (Ljungberg et al., 2012).

4. Experimental Observation in Diverse Material Classes

Cluster/molecular:

• Si$_{10}$H$_{16}$ clusters: Blue shift (TA), red shift (optical/Si-H), mode-dependent broadening, ultrafast mode coupling (Han et al., 2013).
• Organic excimer and charge-transfer states: Temperature-dependent emission lineshapes described by Franck–Condon displacement models with T-weighted Boltzmann factors; hot-band structure and linewidth evolution enable extraction of vibrational and binding parameters across 10–350 K (Hammer et al., 2022).
• Double-well anharmonic modes (e.g., thiourea lattices): Temperature-dependent “soft mode” terahertz signatures captured by explicit quantum solution of vibrational double-well Hamiltonians, revealing broadening and blue shift of low-frequency IR peaks with increasing T (Mitoli et al., 2024).

Extended solids:

• h-BN: Direct tracking of phonon softening and multi-phonon (Umklapp) processes up to 1300 K via momentum-resolved vibrational EELS; non-linear T-dependence and necessity to include fourth-order force constants (O'Hara et al., 2023).
• FeGa$_3$: Infrared and Raman spectra with exceptionally weak T-dependence (shifts $<$1%, broadening $<$2 cm$^{-1}$ from 4–300 K), characteristic of weakly anharmonic, weakly correlated semiconductors (Martin et al., 2023).
• Hybrid perovskites: Dramatic linewidth broadening and phonon mode softening across symmetry-breaking phase transitions, with dynamic disorder and mode–mode coupling reflected in temperature-dependent Raman/THz spectra (Leguy et al., 2016, Brivio et al., 2015).

Liquids:

• Water: Raman and IR spectra from 273–373 K mapped as linear combinations of two T-independent "basis spectra", with band shifts and intensities linearly correlated to the instantaneous local tetrahedral order parameter $q$ (Morawietz et al., 2017, Vuilleumier et al., 2023). Stretching modes blue-shift with T at +0.2–0.6 cm$^{-1}$ K$^{-1}$; hydrogen bond loss per K scales all peak shifts and intensities.

5. Computational Approaches and Algorithms

Methodology	Key Features	Material Applicability
Quasi-harmonic approximation	Lattice-expansion-induced frequency shifts	Crystals with moderate anharmonicity
AIMD + Mode Projection	Individual mode energy tracking; autocorrelation-based coupling extraction	Clusters, molecular crystals
Classical Moment Expansion	Effective frequencies, widths from static averages	Generic classical many-body systems
Self-consistent phonon theory	Explicit anharmonic self-energies (cubic, quartic)	High-T solids, strong multi-phonon effects
Franck–Condon progression (T>0)	Boltzmann-weighted quantum or semiclassical emission/absorption spectra	Organic molecular aggregates, excimers
Thermofield Dynamics & Split-Operator	Numerically exact finite-T vibronic spectra by mapping mixed states to doubled-pure-states	Molecules and clusters (low-dimensional)
Herman–Kluk Thermofield Approximations	Semiclassical trajectory ensemble captures hot-band emergence and anharmonic T-dependence	General multidimensional anharmonic systems

The choice of method depends on system size, degree of anharmonicity, desired spectral features (peak positions, widths, line shapes), and computational cost. For strongly anharmonic or dynamically disordered materials, explicit molecular dynamics or thermofield-based quantum dynamics are typically required (Zhang et al., 2023, Kröninger et al., 29 Jul 2025).

6. Spectroscopic Consequences and Applications

Red/blue shifts and broadening: Analysis of temperature dependent spectral data allows extraction of Grüneisen parameters, anharmonic constants, mode-specific couplings, and thermodynamic quantities (e.g., vibrational free energy, heat capacity, entropy, thermal expansion) (Han et al., 2013, O'Hara et al., 2023, 1901.10587, Zhuang et al., 2021).

Mode lifetimes and energy transfer: Shortening of lifetimes with T and increasing mode-coupling strengths result in ultrafast energy transfer processes, directly impacting phonon-limited electronic, thermal, and optical properties in nanostructures, hybrid perovskites, and soft solids (Leguy et al., 2016, O'Hara et al., 2023).

Structure–spectra correlation: In water and complex liquids, the mapping of spectral shifts and intensities to structural order parameters (e.g., tetrahedrality $q$) provides a basis for non-invasive structure determination and monitoring of local environments in aqueous/biological systems (Morawietz et al., 2017).

Coherent vs. spontaneous Raman thermometry: In spontaneous Raman, the Stokes/anti-Stokes intensity ratio is directly sensitive to temperature via Bose–Einstein factors. For coherent stimulated Raman (SRS, CARS), this sensitivity is suppressed for harmonic modes due to destructive pathway interference, but restored for anharmonic ("hot-band") transitions, which can be quantitatively analyzed to extract temperature without calibration of polarizability derivatives (Batignani et al., 2024).

Discrimination of vibronic vs. electronic coherence: The temperature-dependence of relative oscillation phases in two-dimensional vibrational spectra enables clear distinction of vibrational (vibronic) from electronic coherence, with vibronic beat phases showing strong T-dependence due to vibrational Boltzmann factors (Perlík et al., 2013).

7. Implications, Limitations, and Outlook

Temperature-dependent vibrational spectra encode vital information on bonding, structure, and dynamics across scales and materials. Accurate modeling and interpretation require inclusion of thermal expansion, explicit anharmonicities, mode–mode coupling, and quantum/thermal population statistics. Conventional quasi-harmonic and harmonic treatments are insufficient for systems with strong soft-mode or double-well physics, high-T, or dynamic disorder.

Emerging computational frameworks—moment-based static approaches, semiclassical thermofield methods, and explicit high-dimensional quantum dynamics—continue to extend the accessible range of systems and temperatures. Cross-validation against high-resolution, momentum-resolved experiments is essential for quantifying the role of multiphonon/Umklapp processes in transport and phase transitions (O'Hara et al., 2023, Zhang et al., 2023, Kröninger et al., 29 Jul 2025).

In sum, temperature-dependent vibrational spectroscopy, underpinned by advanced theoretical models and computational methods, provides a powerful window into the underlying atomic-scale mechanisms controlling functional properties of materials, clusters, and complexes across scientific disciplines.