Klemens Process: Anharmonic Decay & Oxide Synthesis

• Klemens Process is a kinetic mechanism describing the anharmonic three-phonon decay of a zone-center optical phonon into two lower-energy phonons, fundamental to Raman and thermal analyses.
• It extends to modeling Umklapp scattering, point-defect and dislocation effects, and anisotropic thermal conductivity via self-energy and relaxation-time formulations.
• A distinct application in oxide whisker synthesis uses a stirring-driven dissolution–recrystallization process to achieve anisotropic fibrous growth under sustained supersaturation.

Searching arXiv for papers related to "Klemens process" and associated phonon/thermal-transport usages. The Klemens process is not a single universally fixed construct in contemporary materials literature. In the most common condensed-matter usage, it denotes the anharmonic three-phonon decay of a zone-center optical phonon into two lower-energy phonons with equal and opposite wavevectors, and it underlies standard Raman linewidth, phonon lineshift, and high-temperature thermal-transport analyses. Closely related usages extend the Klemens framework to Umklapp-limited lattice thermal conductivity, point-defect and dislocation scattering, and anisotropic transport models. In a distinct oxide-whisker literature, the same term is used for a stirring-driven dissolution–recrystallization route that produces fibrous oxides through anisotropic growth under sustained supersaturation (Falkovsky, 2013, D'Souza et al., 2017, Jiang et al., 2016, Jiang et al., 2023, Baranowski et al., 8 Jul 2025, Abakumov et al., 19 Aug 2025).

1. Range of meanings and unifying structure

Across these literatures, the term consistently refers to a rate-limiting microscopic process that is cast in self-energy or relaxation-time form. In phonon spectroscopy and nonequilibrium phonon kinetics, the process is a cubic-anharmonic optical-phonon decay. In lattice-thermal-conductivity theory, it appears as a Klemens or Callaway–Klemens scattering rate for Umklapp or defect scattering. In dislocation theory, it also includes strain-field and core scattering channels. In the fibrous-oxide literature, it denotes a driven kinetic pathway in which mechanical agitation sustains supersaturation and anisotropic recrystallization (Yadav et al., 2021, Jiang et al., 2023, Jiang et al., 2016, Abakumov et al., 19 Aug 2025).

Usage	Core mechanism	Representative systems
Anharmonic phonon decay	Zone-center optical phonon decays into two phonons with opposite momenta	CoS$_2$, rutile oxides, III–V MQWs
Thermal-transport model	Klemens/Callaway relaxation times for Umklapp or defect scattering	h-BN, Half-Heuslers, graphite
Dislocation scattering	Strain-field and core phonon scattering by dislocations	Si nano-films with 60° shuffle-set arrays
Fibrous-oxide synthesis	Stirring-induced ballistic detachment plus anisotropic recrystallization	V$_2$O$_5$ and related oxides

A common structural feature is the decomposition of the observable into additive or effectively separable scattering channels. In the phonon case this is expressed through the phonon self-energy or through Matthiessen-type sums of inverse relaxation times; in the synthesis case it appears as a competition among thermal attachment, thermal detachment, and athermal detachment probabilities. This suggests that the term is best understood as a family of Klemens-type kinetic constructions rather than as a single invariant mechanism.

2. Anharmonic optical-phonon decay

In its canonical form, the Klemens process is the decay of a zone-center optical phonon of frequency $\omega_0$ into two phonons with equal and opposite momenta,

$$(q=0,\omega_0)\rightarrow (q,\omega_q)+(-q,\omega_q), \qquad \omega_q \approx \omega_0/2,$$

with energy and momentum conservation enforced explicitly in the daughter-state phase space (Kumawat et al., 22 May 2026). In polar III–V semiconductors the same process is formulated as intrinsic cubic-anharmonic decay of a zone-center longitudinal optical phonon into two acoustic phonons, with the isotropic linear-dispersion condition giving $|q_1|=|q_2|$ and $\omega_0 \approx 2v|q_1|$ (Baranowski et al., 8 Jul 2025).

For CoS$_2$, the Raman-active phonon at $\omega = 400~\mathrm{cm}^{-1}$ is treated at $k=0$ with acoustic dispersion $_2$0. The imaginary part of the self-energy selects the equal-energy channel $_2$1 with opposite momenta $_2$2 and $_2$3, identified explicitly as the standard Klemens process. The resulting linewidth is

$_2$4

with

$_2$5

while the lineshift is obtained from the real part of the same self-energy and, at low temperature, is expressed in terms of $_2$6 as a negative $_2$7 correction (Falkovsky, 2013). The same paper emphasizes that this three-phonon term accounts for the baseline Raman linewidth and low-temperature lineshift, but not for the full temperature evolution near ferromagnetic ordering.

The III–V multi-quantum-well treatment extends the same decay channel into a coupled nonequilibrium kinetics problem. The LO and acoustic populations satisfy

$_2$8

where $_2$9 is the forward LO$_5$0AA term and $_5$1 the reverse AA$_5$2LO term, with detailed balance preserved explicitly. In this formulation the Klemens channel does not simply drain the LO population: if the generated acoustic phonons do not escape or thermalize rapidly, the reverse channel raises the effective LO lifetime seen by carriers. For the 5.2 nm and 10 nm barrier MQWs, the observed hot-carrier “turn-on” is captured with $_5$3 and $_5$4, while $_5$5 exceeds $_5$6 in the carrier-coupled $_5$7-range (Baranowski et al., 8 Jul 2025).

In Raman practice, the same three-phonon decay is often written in the symmetric Klemens form

$_5$8

which the rutile-oxide study adopts as its conventional baseline before adding an electronic correction for metallic systems (Kumawat et al., 22 May 2026).

3. Klemens–Callaway transport theory

In thermal-transport theory, the Klemens process enters through relaxation times that are inserted into Boltzmann or Callaway conductivity integrals. For single- and multilayer h-BN, the Callaway–Klemens implementation uses the in-plane conductivity

$_5$9

with branch-resolved Umklapp time

$\omega_0$0

This is the Klemens Debye-plus-Grüneisen form used for both 2D and 3D within that framework (D'Souza et al., 2017). The same study shows that the simplified Callaway–Klemens integrals are useful for trends and analytics, while accurate mode-resolved $\omega_0$1 in layered h-BN requires iterative, beyond-RTA BTE treatment with first-principles IFCs.

For ZrCo$\omega_0$2Ir$\omega_0$3Sb, the Klemens–Callaway picture is used to separate point-defect and Umklapp regimes. The point-defect rate is

$\omega_0$4

with

$\omega_0$5

Experimentally, $\omega_0$6 falls from $\omega_0$7 in ZrCoSb to $\omega_0$8 in ZrCo$\omega_0$9Ir$$(q=0,\omega_0)\rightarrow (q,\omega_q)+(-q,\omega_q), \qquad \omega_q \approx \omega_0/2,$$0Sb, while the substituted alloys show $$(q=0,\omega_0)\rightarrow (q,\omega_q)+(-q,\omega_q), \qquad \omega_q \approx \omega_0/2,$$1 below $$(q=0,\omega_0)\rightarrow (q,\omega_q)+(-q,\omega_q), \qquad \omega_q \approx \omega_0/2,$$2 and $$(q=0,\omega_0)\rightarrow (q,\omega_q)+(-q,\omega_q), \qquad \omega_q \approx \omega_0/2,$$3 above $$(q=0,\omega_0)\rightarrow (q,\omega_q)+(-q,\omega_q), \qquad \omega_q \approx \omega_0/2,$$4, identifying point-defect dominance at lower temperature and Umklapp dominance at higher temperature (Yadav et al., 2021).

A further generalization is the anisotropic Klemens model for the thermal-conductivity tensor. There the dispersion is taken as

$$(q=0,\omega_0)\rightarrow (q,\omega_q)+(-q,\omega_q), \qquad \omega_q \approx \omega_0/2,$$5

and the high-temperature Umklapp times along principal axes satisfy

$$(q=0,\omega_0)\rightarrow (q,\omega_q)+(-q,\omega_q), \qquad \omega_q \approx \omega_0/2,$$6

Combined with an ellipsoidal mean-free-path construction, this yields analytical expressions for $$(q=0,\omega_0)\rightarrow (q,\omega_q)+(-q,\omega_q), \qquad \omega_q \approx \omega_0/2,$$7 and $$(q=0,\omega_0)\rightarrow (q,\omega_q)+(-q,\omega_q), \qquad \omega_q \approx \omega_0/2,$$8 that preserve Onsager reciprocity and reproduce the measured $$(q=0,\omega_0)\rightarrow (q,\omega_q)+(-q,\omega_q), \qquad \omega_q \approx \omega_0/2,$$9 scaling of bulk graphite at high temperature, while also describing the thickness dependence of cross-plane transport in graphite films (Jiang et al., 2023).

4. Dislocation and core-scattering extensions

In the silicon nano-film study, the Klemens framework is extended to phonon–dislocation scattering by separating long-range strain-field scattering from core scattering. The strain field outside the core gives $|q_1|=|q_2|$0 rates, while the core, modeled as a line of vacancy-like mass defects, gives a stronger $|q_1|=|q_2|$1 term (Jiang et al., 2016). In the simplified form adopted there,

$|q_1|=|q_2|$2

$|q_1|=|q_2|$3

and

$|q_1|=|q_2|$4

The longitudinal and transverse branches are therefore affected differently through $|q_1|=|q_2|$5 and $|q_1|=|q_2|$6, while the line density $|q_1|=|q_2|$7 sets the overall scale (Jiang et al., 2016).

The central revision relative to the original core picture is quantitative. To reproduce molecular-dynamics thermal conductivities for silicon nano-films with embedded 60° shuffle-set dislocation arrays, the study requires $|q_1|=|q_2|$8 rather than a single-vacancy-line estimate near unity. It also finds that core scattering dominates over strain-field scattering by roughly a factor of $|q_1|=|q_2|$9 for both longitudinal and transverse phonons, with longitudinal modes most strongly degraded. The molecular-dynamics analysis further infers an effective scattering diameter of about $\omega_0 \approx 2v|q_1|$0 around each core, which produces a density threshold: the conductivity reduction becomes distinct only when $\omega_0 \approx 2v|q_1|$1 (Jiang et al., 2016).

These results make the Klemens process in defected semiconductors strongly local rather than spatially homogeneous. For $\omega_0 \approx 2v|q_1|$2, the model gives bulk silicon $\omega_0 \approx 2v|q_1|$3 and nano-film $\omega_0 \approx 2v|q_1|$4 at $\omega_0 \approx 2v|q_1|$5. Extrapolation to dislocation spacing $\omega_0 \approx 2v|q_1|$6 with $\omega_0 \approx 2v|q_1|$7 yields $\omega_0 \approx 2v|q_1|$8, about $\omega_0 \approx 2v|q_1|$9 of the perfect-film value and $_2$0 of bulk silicon (Jiang et al., 2016).

5. Beyond the standard Klemens picture: magnetic and electronic self-energies

Several recent studies use the Klemens term as a necessary baseline but not a sufficient explanation. In CoS$_2$1, the full linewidth and lineshift are written additively as

$_2$2

For the $_2$3 Raman phonon, the two-magnon decay is kinematically forbidden because $_2$4 exceeds $_2$5, so $_2$6. The paper concludes that Klemens anharmonicity governs the intrinsic low-temperature linewidth and lineshift, whereas interband electron transitions activated by the temperature-dependent Weiss field are required for quantitative agreement near and above $_2$7 (Falkovsky, 2013).

The rutile-oxide study reaches an analogous conclusion for metallic RuO$_2$8 and IrO$_2$9, but in a non-magnetic context. Conventional three-phonon fits reproduce the qualitative temperature trends, yet direct fits return anharmonic parameters inconsistent with the expectation $\omega = 400~\mathrm{cm}^{-1}$0 for modes above $\omega = 400~\mathrm{cm}^{-1}$1 when the cutoff is $\omega = 400~\mathrm{cm}^{-1}$2. The proposed resolution is a modified Klemens form,

$\omega = 400~\mathrm{cm}^{-1}$3

$\omega = 400~\mathrm{cm}^{-1}$4

where the $\omega = 400~\mathrm{cm}^{-1}$5 term encodes a non-adiabatic electronic contribution to the real part of the phonon self-energy (Kumawat et al., 22 May 2026).

Experimentally, metallic RuO$\omega = 400~\mathrm{cm}^{-1}$6 and IrO$\omega = 400~\mathrm{cm}^{-1}$7 show $\omega = 400~\mathrm{cm}^{-1}$8–$\omega = 400~\mathrm{cm}^{-1}$9 on cooling from $k=0$0 to $k=0$1, whereas insulating TiO$k=0$2 and SnO$k=0$3 show only $k=0$4–$k=0$5. By contrast, linewidth changes $k=0$6–$k=0$7 do not separate cleanly into metallic and insulating classes. The study interprets this asymmetry as a consequence of the weak and phase-space-limited temperature dependence of $k=0$8 in good metals, so that the extra electronic effect appears primarily in the frequency shift rather than in the linewidth (Kumawat et al., 22 May 2026).

A broader implication is that the phrase “Klemens process” in current spectroscopy often refers to the anharmonic reference model against which additional self-energy channels are diagnosed. The absence of Fano asymmetry, for example, does not preclude finite electron–phonon coupling if the real part of the electronic self-energy still renormalizes $k=0$9, as argued explicitly for metallic rutiles (Kumawat et al., 22 May 2026).

6. Distinct synthesis usage: fibrous-oxide growth under stirring

In the whisker and fiber community, the term “Klemens process” has a distinct meaning. It refers to producing oxide whiskers or fibers by mechanically agitating oxide powders in water so that dissolution–recrystallization under sustained supersaturation yields highly anisotropic crystals. The V$_2$00O$_2$01 Monte Carlo study states that it does not use the name explicitly, but models the reported fabrication route by adding a stirring-controlled, athermal ballistic detachment channel to a TLK crystal-growth scheme (Abakumov et al., 19 Aug 2025).

The thermodynamic driving force is written as

$_2$02

with $_2$03 and, for the anisotropic ballistic channel,

$_2$04

capped to the interval $_2$05. For orthorhombic V$_2$06O$_2$07, the reported anisotropic surface energies are $_2$08, $_2$09, and $_2$10, leading to terrace energies $_2$11 with $_2$12, $_2$13, and $_2$14 at $_2$15 (Abakumov et al., 19 Aug 2025).

Because $_2$16 is smallest, growth on the $_2$17 terrace is fastest, so the model predicts elongation along the $_2$18 direction. Under stirring, ballistic detachments feed the solution and maintain supersaturation, while favorable ledges and kinks on $_2$19 capture re-attaching quasiatoms. In ensemble simulations the number of clusters decreases with time, but the mean fiber length grows and the total surface energy increases. The study labels this antiripening and reports that the mean fiber-axis length grows roughly as $_2$20 rather than the $_2$21 associated with LSW ripening (Abakumov et al., 19 Aug 2025).

This usage is conceptually separate from phonon decay and thermal transport. The shared term therefore reflects historical naming rather than a shared microscopic object. A plausible implication is that encyclopedia treatment of the Klemens process must distinguish the phonon-theory meaning from the whisker-growth meaning to avoid a false unification.

7. Assumptions, limitations, and domain of validity

The standard phonon-decay form relies on highly specific approximations: a Raman phonon at $_2$22, symmetric daughter states with opposite momenta, and Debye-like or linearized acoustic dispersions. In CoS$_2$23, the low-temperature lineshift formula is derived under simplified dispersion and density-of-states assumptions, and the paper explicitly warns that Menéndez and Cardona had noted discrepancies “mainly from the poor description of the phonon dispersion curves”; it therefore states that the low-$_2$24 fitting formula should be used with caution. The same work also notes that fourth-order anharmonicity may be required at high temperatures $_2$25, although no explicit four-phonon term is included there (Falkovsky, 2013).

Transport implementations inherit parallel limitations. The h-BN Callaway–Klemens treatment uses isotropic or simplified branch dispersions and, in its analytical form, retains only resistive Umklapp scattering; the authors therefore rely on an iterative beyond-RTA BTE to recover strong ZA contributions, length dependence, and hydrodynamic corrections not captured by the Umklapp-only integrals (D'Souza et al., 2017). The anisotropic graphite model addresses one major formal limitation—violation of Onsager reciprocity by gradient-direction-dependent relaxation times—by insisting that $_2$26 remain a scalar function of $_2$27, but it still assumes a high-temperature Umklapp regime, anisotropic Debye dispersion, and diffuse black boundaries for films (Jiang et al., 2023).

Defect-based extensions are likewise parameter sensitive. In dislocation-scattering theory, $_2$28, $_2$29, $_2$30, and the assumed $_2$31-power laws strongly affect predicted conductivities, while the $_2$32 effective diameter in silicon is inferred from molecular dynamics rather than encoded directly in the analytical relaxation times (Jiang et al., 2016). In the V$_2$33O$_2$34 synthesis model, the ballistic probability is phenomenological, hydrodynamics are unresolved, and the mapping between stirring speed and $_2$35 is not established quantitatively (Abakumov et al., 19 Aug 2025).

Taken together, these constraints define the status of the Klemens process in present research. It is a compact and productive phenomenology for anharmonic decay, transport-limiting scattering, and, in a separate materials-synthesis context, driven anisotropic recrystallization. Its broad utility derives from the simplicity of self-energy and relaxation-time formulations, but quantitative fidelity depends on when additional channels—four-phonon terms, non-adiabatic electron–phonon corrections, magnon processes, explicit boundary scattering, or chemically specific surface kinetics—must be restored.