REVIEW 5 major objections 5 minor 55 references
THz emission from multiple ionized plasma
T0 review · 5 major / 5 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read Including high ionization states in the photocurrent model reproduces THz emission from multiply ionized plasma up to 5e17 W/cm², explaining the staircase energy scaling, spectral shift to 60–70 THz, and multicycle temporal envelope observe
desk verdict Multi-ionization photocurrent model explains high-intensity THz generation and delivers a record 0.2 mJ at 1.1% conversion, but the spectral validation leans on an undetermined absorption length z. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the total THz field obtained from the time derivative of the sum of current densities from each ionization state (Eq. 8), where each current density is built from sequential ADK ionization rates. Its key companions are the plasma-frequency-dependent absorption factor a(ω)=exp(-k_i z) (Eq. 11) and the assumption that the generated spectrum is set only by the time-dependent free-electron density. The ionization rate table (Table 1) supplies the orbital quantum numbers needed to compute ADK rates for each charge state, and the absorption factor converts the bare photocurrent spectrum into the observed spectrum.
What would settle it
Measure the plasma absorption length z and electron density independently (e.g., via interferometry or a separate probe of the plasma column) and check whether the same z reproduces the observed spectra without free adjustment; alternatively, compare spectra acquired at identical electron density but with intentionally different focal geometries, which would change the effective absorption length.
Extended reading notes
Core claim
We show that the photocurrent model, extended to account for all ionization states of the target gas through ADK rates and ground-state depletion, quantitatively reproduces the measured THz energy scaling, spectral shift, and temporal envelope in two-color laser-plasma generation. The 'staircase' rise of THz energy with intensity is caused by successive ionization thresholds; higher ionization states contribute disproportionately because they are created near the pulse peak where the asymmetric push on the electron is maximal. The spectral shift and low-frequency suppression are attributed to plasma reabsorption through the factor a(ω)=exp(-k_i z), with k_i set by the plasma frequency from t
Load-bearing premise
The spectral agreement depends on the absorption factor a(ω)=exp(-k_i z), where the propagation distance z is only stated as 'in the order of few µm' and is not independently measured; if z is tuned to match the data, the spectral comparison does not by itself validate the model.
Editorial extensions
If this is right
- THz energy and spectral shape become predictable from the time-dependent electron density alone, making gas choice, pressure, and laser intensity interchangeable knobs.
- The staircase energy scaling identifies exactly which ionization states contribute at which intensity, allowing optimization of the gas and intensity to reach higher states near the pulse peak.
- The model explains the previously unexplained multicycle temporal envelope of high-intensity THz pulses as a superposition of ionization-state currents that activate at different times within the driving pulse.
- The 1.1% conversion efficiency in krypton with an 800 nm driver suggests that further ionization engineering could push NIR-based sources toward multi-percent efficiency.
- The electron-density clamping due to defocusing sets a practical limit on the low-frequency content of the generated THz pulse, which is directly testable by changing pressure.
- The spectral shift to 60–70 THz is a direct consequence of plasma reabsorption, meaning the THz peak tracks the plasma frequency as sqrt(Ne).
Reading between the lines
- If the electron-density interpretation is correct, the spectral peak position should scale as the square root of the electron density across different gases and pressures; this is a direct test not performed in the paper but implied by Eq. (10).
- The unmeasured propagation distance z in the absorption factor is effectively a free parameter; measuring z independently (e.g., via plasma-length interferometry) would either confirm the model or expose a compensating mechanism.
- The model's prediction that high ionization states contribute most near the pulse peak suggests that pulse shaping or adding a third color to delay the highest ionization states could further boost conversion efficiency beyond the reported 1.1%.
- Applying the same multi-state photocurrent framework to mid-infrared drivers (which already produce higher THz energies at lower intensities) could push conversion efficiencies toward several percent, provided the absorption factor remains valid at those electron densities.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper reports THz generation from two-color laser plasma in noble gases at high driving intensity, reaching 221 µJ with 1.1% conversion efficiency in krypton. The authors extend the photocurrent model to include multiple ionization states up to Kr10+ and Ne8+ using ADK tunneling rates, then multiply the resulting THz spectrum by a plasma-absorption factor exp(-k_i z). They compare the model with their own measurements and claim agreement for the 'staircase' energy scaling, the spectral shift to 60–70 THz, and the multicycle temporal envelope. They attribute all these effects to the time-dependent electron density.
Significance. The experimental results—especially 221 µJ THz pulses with >1% conversion efficiency from an 800 nm driver—are valuable and extend the parameter space of two-color plasma THz sources. Extending the photocurrent model to high ionization stages is a sensible and potentially useful step, and the staircase interpretation linked to sequential ionization thresholds is compelling. The paper would be significant if the spectral validation were independent of a tunable parameter. However, the current manuscript overstates the level of validation: the energy comparison is in arbitrary units and qualitative, and the spectral shift is produced by an absorption factor whose propagation distance z is not independently determined.
major comments (5)
- [Methods, Eq. (11)] The spectral shift and low-frequency suppression are entirely produced by a(ω)=exp(-k_i z). The text states only that z is 'in the order of few µm' and gives no value for any spectrum, no independent measurement or computed estimate, and no test of sensitivity. Because k_i is nonzero only for ω<ω_p, this factor is the sole mechanism suppressing low frequencies. If z is chosen (even implicitly) to match the data, the spectral agreement is a fit rather than an independent validation. Please report z values for each gas/pressure/intensity, justify them from the plasma geometry, and show the spectral variation for a plausible range of z.
- [Fig. 2a and Discussion] The energy comparison in Fig. 2a is in arbitrary units and the text admits a 'small quantitative disagreement' due to propagation and volumetric effects. The 221 µJ and 1.1% conversion efficiency are experimental quantities, not predictions of the model. Yet the abstract and Discussion state that the model 'explains the high conversion efficiency' and shows 'excellent agreement'. Please either provide an absolute calibration of the model energy or reformulate the claims to make clear that the agreement is qualitative and in arbitrary units.
- [Methods, Eq. (8)] Equation (8) as printed is dimensionally inconsistent. It states E_THz(t) ≈ dJ/dt = d/dt Σ_i J_i(t) = d/dt Σ_i ∫ J_i(t') dt'. The integral of the current density is not the current density; the last equality is incorrect. Please correct to E_THz(t) ≈ d/dt Σ_i J_i(t) and ensure all subsequent formulas and figures use the correct expression. This equation is central to the model and the error, if reproduced in the actual manuscript, must be fixed.
- [Methods, ADK ionization rates] The authors use original ADK rates for ionization stages above the barrier-suppression limit, while acknowledging that correction formulas (refs [48],[49]) are not applied to ions. The depletion argument is plausible but unquantified: the staircase positions and relative yields depend on the rate shapes near threshold. Please quantify the sensitivity of the predicted energy and spectra to rate uncertainties, or benchmark against corrected rates at least for the lowest stages where data exist.
- [Fig. 4 and temporal-domain discussion] The simulated temporal profile shows only the first half-cycle: the text explicitly states that plasma oscillations and recombination are not solved, and that the current density remains at a DC value after the pulse. The multicycle character is then inferred from autocorrelation traces of the measured field. This is an indirect comparison. The claim of a 'full description' of the temporal structure is overstated; please state the limitation more prominently and avoid asserting complete temporal agreement.
minor comments (5)
- [Methods, Eq. (9)] Please define all symbols in Eq. (9), including c0 and the branch of the square root used for ω<ω_p. The imaginary part k_i should be written explicitly.
- [Figures 2b, 2c, 3d, 3e] The autocorrelation insets lack time-scale axes and units. Adding a calibrated delay axis would make the multicycle claim easier to evaluate.
- [Experimental spectra] The text notes water absorption between 45–50 THz but the comparison with simulations does not appear to account for it. Please state whether any correction was applied or restrict the comparison to spectral regions unaffected by water lines.
- [Discussion] The phrase 'one and general parameter of THz generation is time-dependent electron density' is too strong: the model still depends on ionization potentials and species-specific rates. Please soften to say that electron density is the dominant controlling parameter in the observed regime.
- [References] Several reference entries are corrupted (e.g., Ref [1], Ref [21], Ref [47]). Please correct the reference list in the final version.
Circularity Check
Spectral shift agreement relies on the undetermined absorption distance z in Eq. (11), making the spectral 'prediction' partly a fit.
-
fitted input called prediction
[Methods, 'Plasma current model of THz generation', Eqs. (9)-(11)]
"To address the plasma frequency propagation cut of the resulting spectra of E_THz(t), we employ the spectrally dependent absorption factor a(ω)=e^{-k_i(ω)z}, where z is the propagation distance in the order of few µm and k_i is the imaginary component of k-vector k(ω)=ω/c0 sqrt(1 - ω_p^2/ω^2). ... The final spectra presented in the paper are obtained applying the frequency dependent absorption S_THz(ω)=a(ω)S'_THz(ω)."
The paper presents the agreement of simulated spectra with measured spectra (Figs. 2b,c; 3d,e) as validation, attributing the low-frequency suppression and 60-70 THz peak shift to plasma reabsorption. That reabsorption is entirely encoded in a(ω)=exp(-k_i z): k_i is nonzero only below ω_p, so the low-frequency cutoff shape is set by ω_p (computed from the model's own electron density) and by z. The paper gives z only as 'in the order of few µm' and does not state a measured or independently estimated value for any gas, pressure, or intensity. Because S_THz(ω)=a(ω)S'_THz(ω), z directly controls the depth and position of the low-frequency suppression. If z is chosen to reproduce the measured spectra, the spectral agreement is not an independent test of the model but a fit of the very data it
full rationale
The central energy-staircase result is not circular: it is computed from ADK rates and the rate-equation photocurrent model (Eqs. 4-8) with no fitted target-data parameters, and the paper itself notes a small quantitative residual due to propagation/volumetric effects. The temporal multicycle structure is also a forward consequence of the time-dependent ionization currents. However, the spectral shift and low-frequency suppression—reported as key successes—are produced by the absorption factor a(ω)=exp(-k_i z) whose propagation distance z is not independently determined. Since the comparison spectra are normalized and the paper does not state the z value(s) used, the spectral agreement is partially a fit. This prevents the spectral validation from being fully self-consistent, but it does not nullify the independent energy-scaling content; hence score 6 rather than 8-10.
Assumptions & free parameters
free parameters (1)
- Propagation length z in plasma absorption factor =
unspecified, 'in the order of few um' (Methods, Eq. 11)
assumptions (4)
- domain assumption ADK ionization rates apply to sequential multi-ionization of noble gases with the l,m assignment in Table 1, even above the barrier-suppression limit.
- domain assumption The THz field is proportional to the time derivative of the total current density, with electron velocity determined only by the laser electric field, neglecting collisions, plasma oscillations, and recombination after the pulse.
- domain assumption The plasma absorption of THz is described by a plane-wave propagation through a homogeneous plasma with refractive index sqrt(1 - omega_p^2/omega^2) over a short distance z.
- domain assumption The two-color driving field is represented as E(t) = E01(t) cos(omega1 t) + E02(t) cos(omega2 t + phi) with 5% SH, phase delay pi/2, and 40 fs pulse duration.
Cite this review
Pith. "Pith review of THz emission from multiple ionized plasma." pith.science (2026). https://pith.science/paper/LZSTDEJL
@misc{pith2026250815462,
author = {Pith},
title = {Pith review of: THz emission from multiple ionized plasma},
year = {2026},
howpublished = {\url{https://pith.science/paper/LZSTDEJL}},
note = {Machine review of arXiv:2508.15462}
}
read the original abstract
Studies employing nonlinear interactions of THz pulses are nowadays a promising scientific research field. To capture these phenomena, THz pulses with energy ranging from hundreds of \muJ to the mJ level are necessary. However, techniques that provide pulses with such energy levels are still not widely established. Upscaling methods of laser-solid interaction is limited by the damage threshold of materials, while the mechanism of THz generation from high intensity laser-gas interactions is not fully understood yet. Here, we establish the photocurrent model of laser-driven plasma THz generation in the high-intensity regime by accounting for high-ionization states of the target gas. Our model shows excellent agreement with experimental observations, provides a clear explanation of phenomena in both spectral and temporal domains, and explains the high conversion efficiency from laser to THz. In the experiments, we achieved a generation of 0.2 mJ THz pulses, driven by a Ti:sapphire laser with a conversion efficiency exceeding 1 %.
Reference graph
Works this paper leans on
-
[48]
ong, X M , et al , „ irical or ula or tatic ield ionization rate o ato and molecules by lasers in the barrier- u re ion regi e,“ J. Phys. B: At. Mol. Opt. Phys. 38 2593, 2005
work page 2005
-
[49]
Zhang, Q , et al , „ irical or ula or over-barrier strong- ield ionization,“ Phys. Rev. A 90, 043410, 2014
work page 2014
-
[1]
onouchi, M , „Cutting-edge terahertz technology,“ Nature Photon 1, 97–105, 2007
work page 2007
-
[2]
Zhang, X , et al , „ tre e terahertz cience,“ Nature Photon 11, 16–18, 2017
work page 2017
-
[3]
alén, P , et al , „Matter ani ulation ith e tre e terahertz light: Progre in the ena ling z technology,“ Physics Reports 836-837, 2019
work page 2019
-
[4]
icario, C , et al , „ u cycle tre e onlinearitie in GaP nduced y an Ultra trong erahertz ield,“ Phys. Rev. Lett. 118, 083901, 2017
work page 2017
-
[5]
hen, Y , et al , „ onlinear Cro -Phase Modulation with Intense Single-Cycle Terahertz Pul e ,“ Phys. Rev. Lett. 99, 043901, 2007
work page 2007
-
[6]
hen, Y , et al , „ lectro-optic time lensing with an intense single-cycle terahertz ul e,“ Phys. Rev. A 81, 053835, 2010
work page 2010
Show all 55 references
-
[7]
Mat unaga, R , et al , „ igg A litude Mode in the BC u erconductor 1-x TixN nduced y erahertz Pul e citation,“ Phys. Rev. Lett. 111, 057002, 2013
2013
-
[8]
Manko ky, R , et al , „ onlinear lattice dyna ic a a a i or enhanced u erconductivity in YBa2Cu O 5,“ Nature volume 516, pages71–73 , 2014
2014
-
[9]
Express 29, 5982-5992, 2021
ort, O , et al,, „ igh-order ara etric generation o coherent XU radiation,“ Opt. Express 29, 5982-5992, 2021
2021
-
[10]
Birulia, A , et al , „Generation o atto econd ul e ith a controlla le carrier-envelope phase via high-order re uency i ing,“ Phys. Rev. A 106, 023514, 2022
2022
-
[11]
Birulia, A , et al , „Macro co ic e ect in generation o atto econd XU ul e via high-order re uency i ing in ga e and la a,“ New J. Phys. 26 023005, 2024
2024
-
[12]
ovác , , et al , „Qua i-Phase-Matching High-Harmonic Radiation Using Chirped THz Pul e ,“ Phys. Rev. Lett. 108, 193903, 2012
2012
-
[13]
óth, G , et al , „ ingle-cycle attosecond pulses by Thomson backscattering of terahertz ul e ,“ J. Opt. Soc. Am. B 35, A103-A109 , 2018
2018
-
[14]
Balogh, , et al , „ ingle atto econd ul e ro terahertz-assisted high-order harmonic generation,“ Phys. Rev. A 84, 023806, 2011
2011
-
[15]
Zhang, D , et al , „ eg ented terahertz electron accelerator and ani ulator AM ,“ Nature Photon 12, 336–342 , 2018
2018
-
[16]
uang, R , et al , „ erahertz-driven, all-o tical electron gun,“ Optica 3, 1209-1212, 2016
2016
-
[17]
anni, A , et al , „ erahertz-driven linear electron acceleration,“ Nat Commun 6, 8486, 2015
2015
-
[18]
óth, G , et al , „ ilted ul e ront u ing techni ue or e icient terahertz ul e generation,“ Light Sci Appl 12, 256, 2023
2023
-
[19]
Express 22, 20155-20163, 2014
ülö , J A , et al , „ icient generation o z ul e ith 4 J energy,“ Opt. Express 22, 20155-20163, 2014
2014
-
[20]
icario, C , et al , „Generation o 9-mJ THz pulses in DSTMS pumped by a Cr:Mg2 iO4 la er,“ Opt. Lett. 39, 6632-6635, 2014
2014
-
[21]
Xie, X , et al , „Coherent Control o z ave Generation in A ient Air,“ Phys. Rev. Lett. 96, 075005, 2006
2006
-
[22]
i , Y , et al , „Coherent control o terahertz u ercontinuu generation in ultra a t laser–ga interaction ,“ Nature Photon 2, 605–609, 2008
2008
-
[23]
Mitro anov, A , et al , „Ultraviolet-to-millimeter-band supercontinua driven by ultrashort mid-in rared la er ul e ,“ Optica 7, 15-19, 2020. 12
2020
-
[24]
uk, D , et al , „Generation o cala le terahertz radiation ro cylindrically ocu ed t o- color la er ul e in air,“ Appl. Phys. Lett. 108, 121106, 2016
2016
-
[25]
oulouklidi , A D , et al , „O ervation o e tre ely e icient terahertz generation ro mid-infrared two-color la er ila ent ,“ Nat Commun 11, 292, 2020
2020
-
[26]
Jang, D , et al , „ icient terahertz and Brunel har onic generation ro air la a via mid-in rared coherent control,“ Optica 6, 1338-1341, 2019
2019
-
[27]
edorov, Y , et al , „ tre e z ield ro t o-color filamentation of midinfrared la er ul e ,“ Phys. Rev. A 97, 063842, 2018
2018
-
[28]
Clerici, M , et al , „ avelength caling o erahertz Generation y Ga onization,“ Phys. Rev. Lett. 110, 253901, 2013
2013
-
[29]
Express 33, 22610-22631, 2025
u, B , et al , „ icient terahertz radiation generation ro air la a ith id- and near-infrared three-color e to econd ul e ,“ Opt. Express 33, 22610-22631, 2025
2025
-
[30]
Photon., 2025
ei z, L , et al , „ ave or -controlled field synthesis of sub-two-cycle pulses at the 100 eak o er level,“ Nat. Photon., 2025
2025
-
[31]
Viotti, A.-L , et al , „Multi-pass cells for post-co re ion o ultra hort la er ul e ,“ Optica 9, 197-216, 2022
2022
-
[32]
oth, , et al , „ YLO la er – the frontier of few-cycle, multi-TW, kHz lasers for ultra a t a lication at e tre e light in ra tructure atto econd light ul e ource,“ J. Phys. Photonics 2 045003, 2020
2020
-
[33]
SPIE PC12577, 2023
Anti enkov, R , et al , „L1 Allegra la er at L Bea line acility a a driver or electron acceleration at 1 k z re etition rate,“ Proc. SPIE PC12577, 2023
2023
-
[34]
Buldt, J , , „Ga -plasma- a ed generation o road and terahertz radiation ith 4 average o er,“ Opt. Lett. 46, 5256-5259, 2021
2021
-
[35]
Express 15, 4577-4584, 2007
i , Y , et al , „ erahertz e i ion ro ultra a t ionizing air in y etry-broken la er ield ,“ Opt. Express 15, 4577-4584, 2007
2007
-
[36]
Zhang, Z , et al , „ ectral tailoring o the terahertz radiation ro air la a e cited y two-color e to econd ul e ,“ Appl. Phys. Lett. 123, 031108, 2023
2023
-
[37]
Zhao, J , et al , „ e oral-spatial manipulation of bi-focal bi-chromatic fields for terahertz radiation ,“ Commun Phys 7, 408, 2024
2024
-
[38]
Ana ta io , L , et al , „ nhanced terahertz ave e i ion ro air-plasma tailored by a ru tly auto ocu ing la er ea ,“ Optica 3, 605-608, 2016
2016
-
[39]
Zhang, Z , et al , „Mani ulation o olarization or road and terahertz ave e itted ro la er la a ila ent ,“ Nature Photon 12, 554–559, 2018
2018
-
[40]
Oh, , et al , „ nten e terahertz generation in t o-color laser filamentation: energy caling ith tera att la er y te ,“ New J. Phys. 15 075002, 2013
2013
-
[41]
edorov, Y , et al , „ z generation y t o-color femtosecond filaments with complex polarization states: four- ave i ing ver u hotocurrent contri ution ,“ Plasma Phys. Control. Fusion 59 014025, 2017
2017
-
[42]
oulouklidi , A D , et al , „ ectral and idth caling la and recon truction o z wave packets generated from two-color la er la a ila ent ,“ Phys. Rev. A 93, 033844, 2016
2016
-
[43]
Plasmas 16, 056706, 2009
i , Y , „Generation o coherent terahertz radiation in ultra a t la er-ga interaction ,“ Phys. Plasmas 16, 056706, 2009
2009
-
[44]
edorov, Yu , et al , „ z generation y t o-color femtosecond filaments with complex polarization states: four- ave i ing ver u hotocurrent contri ution ,“ Plasma Phys. Control. Fusion 59 014025, 2017. 13
2017
-
[45]
Monot, P , et al , „ ocu ing li it o a tera att la er in an underden e la a,“ J. Opt. Soc. Am. B 9, 1579-1584, 1992
1992
-
[46]
Pta hnik, , et al , „ ater va or continuu a or tion in the 2 7 and 25 μ and at decrea ed te erature ,“ Atmos Ocean Opt 29, 211–215, 2016
2016
-
[47]
Tunnelling ionization of complex atoms and of atomic ions in an alternating electromagnetic field.,
Ammosov, M., et al., "Tunnelling ionization of complex atoms and of atomic ions in an alternating electromagnetic field.," Sov. Phys. JETP 64, p. 191–1194, 1986
1986
-
[50]
Re e, , et al , „Pheno enological rate or ula or over-barrier ionization of hydrogen and heliu ato in trong con tant electric ield ,“ arXiv:2502.02697, 2025
2025 arXiv
-
[51]
Bauer, D , et al , „ act ield ionization rate in the arrier-suppression regime from numerical time-dependent Schrödinger-e uation calculation ,“ Phys. Rev. A 59, 569, 1999
1999
-
[52]
o tyukov, Yu , et al , „ ield ionization in hort and e tre ely inten e la er ul e ,“ Phys. Rev. A 98, 043407, 2018
2018
-
[53]
Mironov, A A , et al , „ trong-field ionization in particle-in-cell i ulation ,“ arXiv:2501.11672, 2025
2025 arXiv
-
[54]
uk, D , et al , „Generation o cala le terahertz radiation ro cylindrically ocu ed t o- color la er ul e in air,“ Appl. Phys. Lett. 108, 121106 , 2016
2016
-
[55]
Xie, X , et al , „Coherent Control o z ave Generation in A ient Air,“ Phys. Rev. Lett. 96, 075005 , 2006
2006
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.