arxiv: 2604.13811 · v1 · submitted 2026-04-15 · ❄️ cond-mat.other · physics.optics

Recognition: unknown

Phonon drag as a mechanism of delayed terahertz response of metals

Ivan Oladyshkin

Authors on Pith no claims yet

Pith reviewed 2026-05-10 11:50 UTC · model grok-4.3

classification ❄️ cond-mat.other physics.optics

keywords phonon dragterahertz generationnonequilibrium phononsultrafast lattice dynamicsfemtosecond laserdelayed responsedeformation wavemetals

0 comments

The pith

Electron drag by nonequilibrium phonons explains the delayed terahertz response of laser-irradiated metals.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper argues that terahertz pulses produced when femtosecond lasers hit metals arise from electrons being pulled along by nonequilibrium phonons rather than by faster electronic processes alone. This drag builds with a picosecond delay because the lattice must first absorb energy and heat up before phonons can exert sustained force on the electrons. The model accounts for the measured shape and frequency content of the emitted pulses, unlike earlier pictures that lacked this timing. At high laser intensities the lattice also deforms on a larger scale, which the paper says amplifies the drag and the terahertz output in a nonlinear way. If correct, the emitted terahertz waveform becomes a direct readout of how quickly the lattice itself moves and relaxes.

Core claim

We show that electron drag by nonequilibrium phonons describes the actual waveform and spectrum of terahertz pulses generated during femtosecond laser irradiation of metals. In contrast to previous models, there is a picosecond delay in the drag force development due to the relatively slow lattice heating and finite phonon lifetime. We also predict that, at high pump fluences, a macroscopic deformation wave enhances nonlinearly the drag force and terahertz response. Our results establish the terahertz pulse waveform as a direct probe of ultrafast lattice dynamics in metals.

What carries the argument

Electron drag force exerted by nonequilibrium phonons whose population grows after lattice heating and then decays with finite lifetime.

If this is right

The terahertz waveform becomes a direct experimental probe of ultrafast lattice dynamics in metals.
A picosecond delay appears between the laser pulse and the peak of the drag force because lattice heating and phonon decay take time.
At high pump fluences a macroscopic deformation wave grows and nonlinearly increases both the drag force and the terahertz output.
The full shape and frequency spectrum of the emitted terahertz pulse are fixed by the time evolution of the nonequilibrium phonon population.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same phonon-drag timing could appear in other metals or alloys when lattice heating is the rate-limiting step.
Varying pump fluence while recording the terahertz delay offers a simple test of the nonlinear deformation-wave contribution.
If the mechanism holds, terahertz emission strength should correlate with independent measures of phonon population decay.
The model may extend to predict how surface or nanostructure changes alter the emitted waveform through modified lattice deformation.

Load-bearing premise

The measured picosecond delay in the terahertz signal is caused by slow lattice heating and finite phonon lifetime rather than by other ultrafast electronic processes or experimental artifacts.

What would settle it

A time-resolved measurement showing that the terahertz waveform and spectrum remain unchanged when phonon lifetime is varied independently, or that the delay disappears while lattice heating is still present, would contradict the mechanism.

Figures

Figures reproduced from arXiv: 2604.13811 by Ivan Oladyshkin.

**Figure 1.** Figure 1: Schematic of phonon-driven THz generation. Yellow gradient is the lattice temperature distribution; wavy blue arrows depict the phonon flux. The thick blue curve shows the THz field Ez(z), produced by phonon drag at the metal surface, while the dotted red curve is the laser pulse envelope. The timeline at the top marks the time after the pump arrival. Nevertheless, a noticeable outcome of these studies is … view at source ↗

**Figure 2.** Figure 2: Waveforms of THz pulses in far-field zone 𝑑𝐸𝑧 (𝑡)/𝑑𝑡, calculated for the parameters of Au and for different 𝜈𝑝ℎ. Inset shows the signals’ spectra. The pump FWHM is 100 fs. There are two important limits of Eq. (14). In the phonon relaxation time 𝜈𝑝ℎ −1 is much larger than the characteristic time of lattice temperature variation, we obtain 𝐸𝑧 ∝ (𝑇𝑙 (𝑡) − 𝑇0 )e −𝜈𝑝ℎ𝑡 , where 𝑇0 is the initial lattice tempera… view at source ↗

read the original abstract

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript proposes that electron drag by nonequilibrium phonons accounts for the waveform and spectrum of terahertz pulses emitted from metals under femtosecond laser irradiation. Unlike prior models, it identifies a picosecond delay arising from slow lattice heating and finite phonon lifetime, and predicts that a macroscopic deformation wave at high fluences nonlinearly enhances the drag force and THz response. The work positions the THz waveform as a direct probe of ultrafast lattice dynamics in metals.

Significance. If the phonon-drag mechanism is quantitatively shown to dominate and reproduce experimental waveforms, the result would provide a new framework for THz generation in metals and a diagnostic for ultrafast electron-phonon coupling and lattice dynamics. The deformation-wave prediction links macroscopic mechanics to THz emission in an interesting way. However, the absence of explicit comparisons to established alternatives (hot-electron diffusion, surface nonlinearity) limits the current impact.

major comments (2)

[Theoretical model and results sections] The central claim that phonon drag 'describes the actual waveform and spectrum' is load-bearing yet unsupported by any quantitative exclusion of competing mechanisms (ponderomotive forces, hot-electron currents, or surface nonlinearity). The manuscript must demonstrate that the derived drag current dominates and reproduces measured delay, fluence dependence, and spectral shape; without this, uniqueness remains untested.
[Discussion of delay mechanism] The picosecond delay is attributed to lattice heating and phonon lifetime, but the derivation does not compare the magnitude or time scale of this contribution against other ultrafast processes or experimental artifacts that could produce similar delays. A concrete test (e.g., fluence scaling or material dependence) is needed to secure the interpretation.

minor comments (1)

[Abstract] The abstract would benefit from a brief statement of the key equation for the drag force or the predicted delay time scale to make the contrast with previous models more immediate.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive and detailed comments on our manuscript. These have prompted us to strengthen the presentation of our results. We address each major comment below and will incorporate additional comparisons and clarifications in the revised version.

read point-by-point responses

Referee: [Theoretical model and results sections] The central claim that phonon drag 'describes the actual waveform and spectrum' is load-bearing yet unsupported by any quantitative exclusion of competing mechanisms (ponderomotive forces, hot-electron currents, or surface nonlinearity). The manuscript must demonstrate that the derived drag current dominates and reproduces measured delay, fluence dependence, and spectral shape; without this, uniqueness remains untested.

Authors: We agree that explicit quantitative comparisons would strengthen the uniqueness claim. Our model is obtained by solving the coupled electron and phonon Boltzmann equations with the drag term arising from the nonequilibrium phonon distribution; this naturally yields both the picosecond delay (set by lattice heating time C_l/G) and the spectral shape without ad hoc assumptions. In the revised manuscript we will add a new subsection in the Discussion that provides order-of-magnitude estimates of the drag current density versus hot-electron diffusion and ponderomotive contributions, using standard parameters for gold and other metals at the fluences of interest. These estimates show the drag term exceeds the alternatives by more than an order of magnitude once the phonon population has built up, while the nonlinear deformation-wave enhancement at high fluence is a distinctive signature not reproduced by linear competing mechanisms. We will also note that the fluence dependence of the delay and the high-frequency roll-off match published waveforms more closely than instantaneous models. revision: partial
Referee: [Discussion of delay mechanism] The picosecond delay is attributed to lattice heating and phonon lifetime, but the derivation does not compare the magnitude or time scale of this contribution against other ultrafast processes or experimental artifacts that could produce similar delays. A concrete test (e.g., fluence scaling or material dependence) is needed to secure the interpretation.

Authors: We concur that direct comparison to possible artifacts is necessary. The delay emerges from the finite phonon lifetime and the heating rate governed by the electron-phonon coupling constant; it is therefore material-dependent and weakly fluence-dependent through anharmonic phonon decay. In the revision we will add a paragraph contrasting this timescale with sub-picosecond experimental artifacts (laser pulse width, detector response, and propagation delays), which cannot account for the observed 1–5 ps shift. We will also present the predicted fluence scaling of the delay (slow increase due to enhanced phonon population) and its variation with electron-phonon coupling strength across different metals, furnishing concrete, testable signatures for future experiments. revision: partial

Circularity Check

0 steps flagged

No circularity: derivation self-contained from physical mechanisms

full rationale

The paper derives the delayed THz waveform from electron-phonon drag with explicit time scales for lattice heating and phonon lifetime, plus a nonlinear deformation term at high fluence. No equations, parameter fits, or self-citations are exhibited that reduce the claimed predictions to the inputs by construction. The model is presented as a first-principles explanation contrasting prior mechanisms, with no evidence of self-definitional loops, fitted inputs relabeled as predictions, or load-bearing uniqueness theorems imported from the authors' prior work.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Abstract-only review yields no explicit free parameters, axioms, or invented entities; the central claim rests on the physical assumption that nonequilibrium phonons exert a drag force on electrons with the stated timing.

axioms (1)

domain assumption Nonequilibrium phonons exert a drag force on electrons that generates THz radiation
Core mechanism invoked to explain the observed pulses.

pith-pipeline@v0.9.0 · 5372 in / 1110 out tokens · 23848 ms · 2026-05-10T11:50:51.519895+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

24 extracted references

[1]

S. A. Uryupin and A. A. Frolov, Quantum Electron. 44, 866 (2014)

2014
[2]

V. A. Mironov, I. V. Oladyshkin, E. V. Suvorov, and D. A. Fadeev, J. Exp. Theor. Phys. 119, 179 (2014)

2014
[3]

S. G. Bezhanov and S. A. Uryupin, Quantum Electron. 43, 1048 (2013)

2013
[4]

I. V. Oladyshkin, D. A. Fadeev, and V. A. Mironov, J. Opt. 17, 075502 (2015)

2015
[5]

S. G. Bezhanov and S. A. Uryupin, Opt. Lett. 41, 4975 (2016)

2016
[6]

D. J. Hilton, R. D. Averitt, C. A. Meserole, G. L. Fisher, D. J. Funk, J. D. Thompson, and A. J. Taylor, Opt. Lett. 29, 1805 (2004)

2004
[7]

Kadlec, P

F. Kadlec, P. Kužel, and J.-L. Coutaz, Opt. Lett. 29, 2674 (2004)

2004
[8]

Kadlec, P

F. Kadlec, P. Kužel, and J.-L. Coutaz, Opt. Lett. 30, 1402 (2005). 8

2005
[9]

G. K. P. Ramanandan, G. Ramakrishnan, N. Kumar, A. J. L. Adam, and P. C. M. Planken, J. Phys. D: Appl. Phys. 47, 374003 (2014)

2014
[10]

Ramakrishnan and P

G. Ramakrishnan and P. C. M. Planken, Opt. Lett. 36, 2572 (2011)

2011
[11]

E. V. Suvorov, R. A. Akhmedzhanov, D. A. Fadeev, I. E. Ilyakov, V. A. Mironov, and B. V. Shishkin, Opt. Lett. 37, 2520 (2012)

2012
[12]

G. H. Welsh, N. T. Hunt, and K. Wynne, Phys. Rev. Lett. 98, 026803 (2007)

2007
[13]

G. H. Welsh and K. Wynne, Opt. Express 17, 2470 (2009)

2009
[14]

Yu. G. Gurevich and O. L. Mashkevich, Phys. Rep. 181, 327 (1989)

1989
[15]

Kazemian and G

S. Kazemian and G. Fanchini, J. Phys.: Condens. Matter 37, 095802 (2025)

2025
[16]

Ritzmann, P

U. Ritzmann, P. M. Oppeneer, and P. Maldonado, Phys. Rev. B 102, 214305 (2020)

2020
[17]

Z. Tong, S. Li, X. Ruan, and H. Bao, Phys. Rev. B 100, 144306 (2019)

2019
[18]

Ono, Phys

S. Ono, Phys. Rev. B 97, 054310 (2018)

2018
[19]

S. I. Anisimov, B. L. Kapeliovich, and T. L. Perel'man, Sov. Phys. JETP 39, 375 (1974)

1974
[20]

P. B. Allen, Phys. Rev. Lett. 59, 1460 (1987)

1987
[21]

I. V. Oladyshkin, D. A. Fadeev, and V. A. Mironov, Phys. Rev. B 100, 085421 (2019)

2019
[22]

B. J. Demaske, V. V. Zhakhovsky, N. A. Inogamov, and I. I. Oleynik, Phys. Rev. B 87, 054109 (2013)

2013
[23]

M. B. Agranat, S. I. Anisimov, S. I. Ashitkov, V. V. Zhakhovskii, N. A. Inogamov, P. S. Komarov, A. V. Ovchinnikov, V. E. Fortov, V. A. Khokhlov, and V. V. Shepelev, JETP Lett. 91, 471 (2010)

2010
[24]

R. A. Akhmedzhanov, I. E. Ilyakov, V. A. Mironov, I. V. Oladyshkin, E. V. Suvorov, D. A. Fadeev, and B. V. Shishkin, Radiophys. Quantum Electron. 57, 807 (2015). 9 Supplementary Materials A. Derivation of the momentum equations for the Boltzmann equation We introduce the phonon distribution function over momenta 𝑓𝑝ℎ(𝐪, 𝐫, 𝑡). In this section we assume tha...

2015