Pith. sign in

REVIEW 4 major objections 6 minor 16 references

High harmonic generation reflecting the sub-cycle evolution of the Mott transition under a mid-infrared electric field

T0 review · 4 major / 6 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read As a mid-infrared field melts the Mott gap in Sr2CuO3, the emitted high-harmonic peaks shift to lower energies and track the carrier density.

desk verdict The experimental redshift is real and worth reporting, but the sub-cycle Mott-transition mechanism is a plausible simulation-based interpretation, not a demonstrated fact. read the letter →

arxiv 2508.00296 v1 pith:HTPUXFS5 submitted 2025-08-01 cond-mat.str-el physics.optics

Ryohei Ikeda , Yuta Murakami , Daiki Sakai , Tatsuya Miyamoto , Toshimitsu Ito , Hiroshi Okamoto This is my paper · ORCID
classification cond-mat.str-elphysics.optics
keywords high-harmonicgenerationMotttransitioninsulatorstronglycorrelatedelectronsdynamicalmean-fieldtheorydoublon-holonrecombinationSr2CuO3mid-infraredspectroscopy
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper aims to show that high-harmonic generation (HHG, the emission of light at integer multiples of the driving-laser frequency) can read out how a strongly correlated material's electronic structure changes from one half-cycle of an intense mid-infrared field to the next, as the field drives it from a Mott insulator (a material whose conduction is blocked by electron-electron repulsion) toward a metal. In the one-dimensional cuprate Sr2CuO3, once the field amplitude exceeds about 6 MV/cm, the 5th through 13th harmonic peaks move down from the odd multiples of the driving frequency, and the shift grows with harmonic order, saturating near 50 meV for the 11th harmonic. Its field dependence tracks the carrier density measured by pump-probe reflectivity, which reaches roughly 0.25 carriers per site at 12 MV/cm. A dynamical mean-field simulation of the driven Hubbard model reproduces the red-shifts and attributes them to a cycle-by-cycle reconstruction of the electronic structure: the Mott gap shrinks, spectral weight moves into the Drude response, and the radiation phase of doublon-holon recombination drifts, so the accumulated emission interferes constructively below the nominal harmonic frequency. If this is right, HH spectroscopy becomes an all-optical, probe-free way to follow ultrafast phase transitions in correlated materials.

What carries the argument

The load-bearing mechanism is the cycle-by-cycle radiation phase of doublon-holon recombination, quantified by a windowed Fourier transform of the current: the phase $\mathrm{Arg}[J(\omega,t)]$ is evaluated in half-cycle windows of width $T_{hp} = \pi/\Omega$, and the accumulation at the formal harmonic frequency $\omega = n\Omega$ is compared with that at the actual peak $\omega = \omega_{\mathrm{peak},n}$. The calculation that carries the argument is time-dependent dynamical mean-field theory (tdDMFT) on the half-filled single-band Hubbard model on a Bethe lattice with $U = 2.6\,\mathrm{eV}$, bandwidth $W = 2.08\,\mathrm{eV}$, and $T = 260\,\mathrm{K}$, solved with the non-crossing approximation; it supplies the sub-cycle single-particle spectra, carrier density, and radiation phase. The interpretive link is a Doppler-like identity: if the doublon-holon pair energy shifts by $\Delta\epsilon$ per half-cycle, recombination radiation gains an extra phase $(\Delta\epsilon/\hbar)(t_r - t_c)$, where $t_r$ and $t_c$ are recombination and creation times; constructive interference then requires the emission frequency to move by $\Delta\omega_r<0$, i.e., a red-shift when the gap closes. This phase mechanism is distinct from carrier-envelope-phase effects because the calculation shows it is insensitive to the carrier-envelope phase.

What would settle it

Record the emitted harmonic field with sub-cycle temporal resolution (for example, by electro-optic sampling) under the same pump conditions and compare the radiation phase across consecutive half-cycles at $\omega = 9\Omega$ and at the measured peak $\omega_{\mathrm{peak},9}$. The dynamical mean-field picture predicts that at 12 MV/cm the phase at the nominal harmonic wanders cycle by cycle while the phase at the peak stays nearly constant; if both phases stay locked to $9\Omega$, or wander identically, the phase-drift mechanism is ruled out. A complementary test is to run the identical HHG measurement on a rigid-band semiconductor with a comparable gap and field, where the proposed mechanism predicts no red-shift.

Watch

Extended reading notes

Core claim

The central experimental discovery is that once the mid-infrared field exceeds about 6 MV/cm, the harmonic peaks of orders 5 to 13 in Sr2CuO3 shift away from the odd multiples $n\hbar\Omega$: the red-shift grows with $n$ and saturates at roughly 50 meV for the 11th harmonic, and its dependence on the field matches the measured carrier density $n_c$, about 0.25 doublons plus holons per site at 12 MV/cm. The paper's interpretation is that the shift is the direct fingerprint of the Mott transition taking place inside the pulse: with each half-cycle the single-particle spectrum is rebuilt, the Mott gap shrinks, spectral weight transfers to the Drude component, and the antiferromagnetic correlation is reduced, so the phase of the radiation emitted by recombining doublons and holons changes cycle by cycle. In the time-dependent dynamical mean-field calculation, that phase drift makes successive sub-cycle emissions interfere constructively at a frequency lower than $n\Omega$, producing the red-shifted HH peaks. The third harmonic, which arises from in-gap and intraband processes, stays at $3\hbar\Omega$, and the fifth is intermediate.

Load-bearing premise

The argument stands on the assumption that the theoretical model of a single-band Hubbard material on a Bethe lattice, solved with the non-crossing approximation, really behaves like the copper-oxide chains in Sr2CuO3; the cycle-by-cycle phase evolution that explains the red-shift is calculated, not directly measured.

Editorial extensions

If this is right

  • HHG can act as an all-optical, probe-free monitor of the degree of metallisation in a driven Mott insulator: the harmonic red-shift is predicted to track the carrier density cycle by cycle.
  • Rigid-band materials should not show this red-shift under the same excitation conditions, so the effect gives a fingerprint that separates correlated gap dynamics from static band-structure response.
  • The mechanism is insensitive to the carrier-envelope phase and does not require long scattering times, unlike previously reported harmonic shifts in gases, semiconductors, and topological insulators.
  • The harmonic orders separate by physical origin: the third harmonic stays unshifted, the fifth is intermediate, and the seventh and higher harmonics red-shift, so future HHG measurements on correlated materials should analyse each order separately.
  • If the shift indeed measures the carrier density, HH spectroscopy could estimate the density of photo-doped doublon-holon pairs without a separate probe pulse.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Read as an extension, the red-shift magnitude should encode the per-cycle gap-closing rate: measuring $\omega_{\mathrm{peak},n}$ as a function of harmonic order could provide an all-optical estimate of how fast the Mott gap collapses without any probe pulse.
  • A two-pulse version of the experiment (weak pre-pump to create carriers, followed by a weaker mid-infrared field) would test the mechanism by producing a red-shift at lower peak fields than the single-pulse threshold, separating transition-induced phase drift from field-driven nonlinearities.
  • Because the mechanism is set by carrier generation and electronic-structure reconstruction rather than by band rigidity, the same sub-cycle phase-drift picture should apply to other light-induced ordered phases in correlated materials, whose HH fingerprints remain to be worked out.
  • Direct sub-cycle sampling of the emitted field (electro-optic detection of the harmonic transient) would test whether the phase at $n\Omega$ really wanders cycle-by-cycle while the phase at $\omega_{\mathrm{peak},n}$ stays stable; the time-integrated spectra reported here cannot distinguish phase drift from intensity suppression.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. This paper reports an experimental and theoretical study of high-harmonic generation in the one-dimensional cuprate Sr2CuO3 under intense mid-infrared excitation. The experiments show that above a field amplitude of about 6 MV/cm the peak energies of the 5th through 13th harmonics shift downward from odd multiples of the MIR photon energy, with the shift increasing with harmonic order and saturating near 50 meV for the 11th harmonic, while pump-probe reflectivity measurements indicate a threshold-like growth of carrier density in the same field range. Time-dependent dynamical mean-field theory on the single-band Hubbard model with a Bethe lattice and non-crossing approximation is used to argue that the redshifts reflect a cycle-by-cycle reconstruction of the electronic structure—gap collapse, spectral-weight transfer, and doublon-holon phase modification—during the field-induced Mott transition. The paper further argues that this mechanism is CEP-insensitive and negligible in rigid-band systems.

Significance. The main strength is the experimental connection between two independent observables: the HH peak positions and the pump-probe carrier density, both of which display a threshold near 6 MV/cm. The DMFT calculation is not fitted to the measured shifts; the theoretical redshift appears as an output of stated parameters, and the sub-cycle phase analysis is a concrete, falsifiable prediction of the model. If the mechanism survives the fidelity checks noted below, the work would open a useful direction: HH spectroscopy as a probe of sub-cycle electronic-structure dynamics during nonequilibrium phase transitions in correlated materials, beyond the static band-structure information extracted in weakly correlated solids. The experimental analysis is careful, with a plausible assignment of the 3rd harmonic to in-gap/current processes and higher harmonics to the doublon-holon three-step recombination, and the manuscript is clearly written.

major comments (4)
  1. [Methods (DMFT calculations); Fig. 5; Supplementary S5.2] The central mechanistic claim—that the experimental HH redshifts originate from sub-cycle electronic-structure reconstructions—is carried by the tdDMFT simulation on a single-band Hubbard model on a Bethe lattice with the NCA impurity solver (Methods, U=2.6 eV, W=2.08 eV, T=260 K). The experimental spectra are time-integrated, and the cycle-by-cycle phase ΔArg[J(ω,t)] shown in Figs. 5j–l and S9 is a simulation output, not a measured observable. For this attribution to be load-bearing, the fidelity of the model needs independent support: Sr2CuO3 is a charge-transfer insulator whose lowest excitation is an odd-parity exciton (Fig. 2a), and the mapping of the O-2p band to the lower Hubbard band is a heuristic; moreover, NCA is uncontrolled at U/W≈1.25. Please add benchmarks of NCA against a controlled method for the equilibrium spectral function and transient doublon density at these parameters, a sensitivity study in U/W and T, and a comparison of the simulated HH spectra for a noninteracting (rigid-band) reference. Without these checks, the sub-cycle mechanism is plausible but not established.
  2. [Fig. 4h; Methods (HH spectrum measurements)] The experimental redshift is the anchor of the paper, but the peak positions ω_peak,n are presented without error bars or an explicit fitting protocol in Fig. 4h. Please report the statistical uncertainty of each Gaussian centroid, the number of independent spectra per field value, and the sensitivity to the fitting window; also quantify systematic contributions from the 22 meV MIR spectral width, spectrometer calibration, and the wavelength-dependent correction factor. The maximum reported shift (~50 meV at the 11th harmonic) is about twice the MIR spectral width, so this is not a trivial check: the harmonic-order dependence and the saturation near 12 MV/cm should be shown to be significant against these uncertainties.
  3. [Supplementary S4; Eq. (1); Fig. 3e] The conversion from the transient reflectivity change to the absolute carrier density n_c uses two Lorentzians for the exciton and continuum and the relation −ΔI/I = 2n_c (Supplementary S4). The threshold behavior near 6 MV/cm is likely robust, but the claim that the HH redshift 'tracks' n_c quantitatively (Figs. 4g,h) depends on this conversion, which carries an unknown systematic error because the doublon–holon continuum is represented by a Lorentz oscillator and the high-energy extrapolation of the spectral-weight sum is not specified. Please state the uncertainty in n_c from the Lorentzian fits and from the 2n_c assumption, and show whether the apparent offset between the n_c threshold and the HH-shift onset is within that uncertainty.
  4. [Discussion; Abstract] The paper states that the observed redshift mechanism is 'negligible in rigid-band systems,' but no control experiment or noninteracting-band simulation with the same pulse parameters is presented. Because known HH shifts in solids have been attributed to CEP and nonadiabatic effects (refs. 8, 50), and because the present experiment does not measure the sub-cycle phase directly, the distinction from rigid-band behavior should be demonstrated rather than assumed. A measurement on a conventional semiconductor under identical conditions, or a calculation of HH shifts in a noninteracting band model, would test whether gap collapse is uniquely needed to produce the observed harmonic-order-dependent redshift.
minor comments (6)
  1. [Results, first paragraph] 'd band' should be 'd_{x^2-y^2} band'.
  2. [Methods, DMFT calculations] The definition of the pump vector potential uses F_Gauss but does not distinguish the vector-potential amplitude from the reported field amplitude |E_MAX| used in Fig. 5; define E0 and its relation to |E_MIR|.
  3. [Fig. 5j–l and Supplementary S5.2] The reference time for ΔArg is given as t_i=−73.2 fs in the main text and t_ref=−1.59 fs in the supplement; harmonize the notation.
  4. [References and Methods] Reference 2 is formatted as a sentence ('For a recent review...') rather than a standard citation title, and 'Nessi' and 'NESSi' are used inconsistently in the Methods.
  5. [Fig. 3e caption] The caption mentions light-blue and purple broken lines, but the line styles in the figure are not clearly legible in the provided version; please ensure all line styles are identified.
  6. [Methods, HH spectrum measurements] The correction procedure with the standard lamp is described, but the uncertainty of the correction factor near the harmonic energies is not quoted; please comment on its magnitude.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the HH redshifts are an independent DMFT output, not fitted to the measured shifts.

full rationale

The paper's central claim—that HH peak redshifts track carrier density and reflect sub-cycle Mott-gap reconstructions—rests on two independent legs: (i) experimental HH spectra and pump-probe reflectivity measurements of Sr2CuO3 (Figs. 4g,h) and (ii) tdDMFT simulations of the single-band Hubbard model on the Bethe lattice with stated parameters U=2.6 eV, W=2.08 eV, T=260 K and a Gaussian MIR pulse (Methods). The simulated HH peak shifts (Fig. 4j) and carrier density (Fig. 4i) are outputs of the time evolution, not parameters fitted to the experimental redshift; there is no evidence in the text that U, W, T, or the pulse shape were adjusted to reproduce the ~50 meV shift. The sub-cycle radiation-phase analysis (Figs. 5j-l) is a diagnostic applied to the same simulated current, so it illustrates the proposed mechanism rather than providing an independent proof, but this is a limitation of internal consistency, not circularity. Self-citations, such as refs. 28 and 29 for the three-step doublon-holon HHG mechanism, are used as background context and are not load-bearing for the redshift claim, which is carried by the paper's own independent simulation and experimental observables. The possible mismatch between the Bethe-lattice single-band NCA model and the actual 1D charge-transfer cuprate is a correctness/fidelity concern, not a circularity in the derivation chain.

Assumptions & free parameters 5 free parameters · 6 assumptions · 0 invented entities

The central claim rests on mapping the experiment to a single-band Hubbard model solved by tdDMFT, with material-matched parameters (U, W, T), an approximate impurity solver, and a carrier-density conversion rule (-Delta I / I = 2 n_c). The theoretical redshift is not fitted to the target data, so circularity burden is low, but the parameter choices and approximations carry most of the model dependence. No new particles, forces, dimensions, or conserved quantities are introduced.

free parameters (5)
  • U (onsite Coulomb repulsion) = 2.6 eV
    Set in the DMFT simulation to be consistent with Sr2CuO3; not computed from first principles.
  • W (bandwidth, related to hopping) = 2.08 eV
    Chosen to match the Sr2CuO3 band structure; directly controls gap size and dynamics.
  • T (temperature) = 260 K
    Selected so that the model Mott gap is about 1.5 eV and antiferromagnetic order develops, consistent with the material.
  • E_th (tunneling threshold field) = 13.4 MV/cm
    Fitted to the carrier-density data using Gamma proportional to exp(-pi E_th / |E|) in Eq. (1); yields the doublon-holon correlation length xi ~ 6.57 Angstrom.
  • Lorentz oscillator parameters (LO1, LO2, eps_inf) = Table S2 values
    Fitted to the linear reflectivity and epsilon_2 spectra to convert -Delta R / R into carrier density n_c, under the additional assumption -Delta I / I = 2 n_c.
assumptions (6)
  • domain assumption The single-band Hubbard model with Peierls coupling describes the relevant charge dynamics of Sr2CuO3 despite its charge-transfer character.
    Main text Fig. 1d and Methods: the O-2p band is mapped onto the lower Hubbard band, a standard but approximate treatment for CuO chains.
  • domain assumption DMFT on a Bethe lattice with the non-crossing approximation gives a faithful description of the dynamics; lattice geometry is insensitive.
    Methods paragraph: 'we used the Hubbard model on a Bethe lattice to reduce computational cost... Previous studies have shown that DMFT dynamics are generally insensitive to lattice geometry.' NCA is an approximate impurity solver.
  • domain assumption Carrier density is related to reflectivity reduction by -Delta I / I = 2 n_c.
    Supplementary S4: the paper cites rigorous calculations giving a slightly steeper slope but uses the simple relation to convert -Delta R/R into n_c.
  • domain assumption Harmonics with n >= 7 arise from the doublon-holon three-step mechanism, while the 3rd harmonic arises from in-gap and intraband processes.
    Assigned from prior theoretical work on HHG in Mott insulators (refs 28, 29) and used to interpret the different phase behavior of n = 3 versus n >= 7.
  • standard math Spatial inversion symmetry restricts the harmonic spectrum to odd orders.
    Used to explain the observation of only odd-order harmonics in Sr2CuO3.
  • standard math Kramers-Kronig transformation of reflectivity yields valid optical conductivity spectra.
    Used in Methods to derive sigma and epsilon_2 from the polarized reflectivity spectrum.

how reviews work

0 comments
Cite this review

Pith. "Pith review of High harmonic generation reflecting the sub-cycle evolution of the Mott transition under a mid-infrared electric field." pith.science (2026). https://pith.science/paper/HTPUXFS5

@misc{pith2026250800296,
  author       = {Pith},
  title        = {Pith review of: High harmonic generation reflecting the sub-cycle evolution of the Mott transition under a mid-infrared electric field},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HTPUXFS5}},
  note         = {Machine review of arXiv:2508.00296}
}
read the original abstract

Solids in an intense laser field show high-harmonic generation (HHG), which can provide information on carrier dynamics and band structures in weakly correlated systems. In strongly correlated systems, a laser field can induce a transition between the various electronic phases formed by the entanglement of charge, spin, and orbital degrees of freedom via carrier generation. The HHG accompanying this process should contain information on the nonequilibrium electronic-state dynamics along the oscillating field - an aspect that remains unresolved to date. Here, we show that an intense mid-infrared (MIR) pulse induces a Mott insulator-metal transition in a one-dimensional cuprate, Sr2CuO3, the evolution of which is reflected by the spectral features of HHs. When the electric-field amplitude exceeds 6 MV/cm, carriers are efficiently generated and each harmonic frequency decreases from odd multiples of the MIR frequency. Dynamical mean-field theory indicates that these redshifts originate from a series of electronic-structure reconstructions in each electric-field cycle during the melting of the Mott-insulator state, which modifies the radiation phase from carrier recombination cycle-by-cycle. This phenomenon is negligible in rigid-band systems. This experimental-theoretical study confirms that HH spectroscopy research can potentially unravel the sub-cycle dynamics of nonequilibrium phase transitions in correlated materials.

Figures

Figures reproduced from arXiv: 2508.00296 by the authors.

Figure 2
Figure 2. The laser source and generation method for [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

16 extracted references · 16 canonical work pages

  1. [1]

    Kishida, H. et al. Gigantic optical nonlinearity in one-dimensional Mott–Hubbard insulators. Nature 405, 929–932 (2000)

    work page 2000

  2. [2]

    Ono, M. et al. Linear and nonlinear optical properties of one-dimensional Mott insulators consisting of Ni-halogen chain and CuO-chain compounds. Phys. Rev. B 70, 085101 (2004)

    work page 2004

  3. [3]

    Direct Observation of Excitons and a Continuum of One-Dimensional Mott Insulators: A Reflection-Type Third- Harmonic-Generation Studyof Ni-Halogen Chain Compounds

    Ono, M., Kishida, H., and Okamoto, H. Direct Observation of Excitons and a Continuum of One-Dimensional Mott Insulators: A Reflection-Type Third- Harmonic-Generation Studyof Ni-Halogen Chain Compounds. Phys. Rev. Lett. 95, 087401 (2005)

    work page 2005

  4. [4]

    Miyamoto, T. et al. Biexciton in one-dimensional Mott insulators. Commun. Phys 2, 131 (2019)

    work page 2019

  5. [5]

    Yamakawa, T. et al. Dynamical aspects of excitonic Floquet states generated by a phase-locked mid-infrared pulse in a one-dimensional Mott insulator. New J. Phys. 25, 093044 (2023)

    work page 2023

  6. [6]

    & Kitamura, S

    Oka, T. & Kitamura, S. Floquet Engineering of Quantum Materials. Annu. Rev. Condens. Matter Phys. 10, 387 (2019)

    work page 2019

  7. [7]

    Wang, Y . H. et al. Observation of Floquet-Bloch States on the Surface of a Topological Insulator. Science 342, 453-457 (2013)

    work page 2013

  8. [8]

    Merboldt, M. et al. Observation of Floquet states in graphene. Nat. Phys. (2025)

    work page 2025

Show all 16 references
  1. [9]

    et al., Observation of excitonic Floquet states in a one-dimensional organic Mott insulator using mid-infrared pump near-infrared probe reflection spectroscopy

    Ikeda, R. et al., Observation of excitonic Floquet states in a one-dimensional organic Mott insulator using mid-infrared pump near-infrared probe reflection spectroscopy. Preprint at https://arxiv.org/abs/2407.17759

    arXiv

  2. [10]

    Zhou, S. et al. Pseudospin-selective Floquet band engineering in black phosphorus., Nature 614, 75-80 (2023)

    2023

  3. [11]

    Shan, J-Y . et al. Giant modulation of optical nonlinearity by Floquet engineering. Nature 600, 235–239 (2021)

    2021

  4. [12]

    Okamoto, H. et al. Photoinduced transition from Mott insulator to metal in the undoped cuprates Nd2CuO4 and La2CuO4. Phys. Rev. B 83, 125102 (2011)

    2011

  5. [13]

    Okamoto, H. et al. Photoinduced Metallic State Mediated by Spin-Charge Separation in a One-Dimensional Organic Mott Insulator. Phys. Rev. Lett. 98, 037401 (2007)

    2007

  6. [14]

    & Andrzej, M

    Eskes, H. & Andrzej, M. O. Two Hubbard Bands: Weight Transfer in Optical and One-Particle Spectra. Phys. Rev. Lett. 73, 1279 (1994)

    1994

  7. [15]

    Schubert, O. et al. Sub-cycle control of terahertz high-harmonic generation by dynamical Bloch oscillations. Nat. Photon. 8, 119–123 (2014). 18

    2014

  8. [16]

    Schmid, C. P. et al. Tunable non-integer high-harmonic generation in a topological insulator. Nature 593, 385–390 (2021)

    2021

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.