pith. sign in

arxiv: 2605.01388 · v1 · pith:OUVW5RB7new · submitted 2026-05-02 · ⚛️ physics.optics

Tracking doublon-holon dynamics in high-harmonic generation from Mott insulators

Tao-Yuan Du , Hui-Ru Li , Bo Li , Ruifeng Lu This is my paper

Pith reviewed 2026-05-09 18:22 UTC · model grok-4.3

classification ⚛️ physics.optics
keywords high-harmonic generationMott insulatorsdoublon-holon dynamicsHubbard modelstrong-field physicsintraband responseinterband transitions
0
0 comments X

The pith

Mott insulators exhibit a filling-dependent crossover in high-harmonic generation from intraband Bloch response at low density to interband doublon-holon processes near half filling.

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

The paper examines high-harmonic generation in a laser-driven one-dimensional Hubbard chain modeling Mott insulators. By tracking doublon population through projections onto equilibrium Hubbard bands, it separates intraband spin-wave-like currents from interband doublon-holon creation. A clear crossover appears with electron filling: dilute systems show Bloch-like intraband motion, intermediate fillings mix both channels, and near half filling the response becomes interband-dominated, producing a plateau and cutoff set by the Mott gap. Increasing interaction strength U enlarges this gap and localizes carriers, strongly suppressing interband contributions. The analysis links equilibrium spin-charge separation to nonequilibrium strong-field emission and shows how dephasing affects coherence and net doublon accumulation.

Core claim

Projecting the driven dynamics onto equilibrium Hubbard bands reveals that doublon population and its time evolution diagnose two distinct excitation channels: intraband (spin-wave-like) motion at low filling and interband doublon-holon pair creation near half filling. This produces a filling-dependent crossover in the HHG spectrum, with interband processes generating the characteristic plateau and cutoff once the system approaches half filling, while larger U suppresses interband weight through gap enlargement and correlation-induced localization.

What carries the argument

Doublon population dynamics obtained by projecting the time-dependent state onto equilibrium Hubbard bands, used as a diagnostic to decompose intraband (spin-wave-like) and interband (doublon-holon creation) current channels.

If this is right

  • At low filling the HHG response follows Bloch-like intraband oscillations without a clear cutoff.
  • Near half filling interband doublon-holon creation dominates and produces a plateau whose width tracks the Mott gap.
  • Increasing U suppresses interband contributions by enlarging the gap and localizing carriers.
  • Intra- and interband currents flow in opposite directions below the Mott gap, while dephasing selectively damps interband coherence and enhances net doublon accumulation.
  • Time-frequency analysis of the emitted light reveals distinct quantum trajectories that change with filling and appear as filling-dependent emission below the Mott gap.

Where Pith is reading between the lines

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

  • The same doublon diagnostic could be applied to two-dimensional or three-dimensional Hubbard models to test whether the crossover survives beyond one dimension.
  • Doping level may serve as an experimental knob to switch HHG between intraband and interband regimes in real correlated materials.
  • The opposing current flows below the gap suggest a possible route to control net charge transfer or photocurrent direction via laser parameters.

Load-bearing premise

The one-dimensional Hubbard chain treated with exact diagonalization and time-dependent density-matrix propagation captures the essential physics of real three-dimensional Mott insulators driven by strong lasers.

What would settle it

Experimental HHG spectra from doped Mott insulators that show no filling-dependent shift in plateau position or no suppression of high harmonics with increasing interaction strength.

Figures

Figures reproduced from arXiv: 2605.01388 by Bo Li, Hui-Ru Li, Ruifeng Lu, Tao-Yuan Du.

Figure 1
Figure 1. Figure 1: (c) shows the doublon number associated with each eigenstate. The distinct doublon distributions under the two filling conditions imply that the corresponding light-driven doublon dynamics lead to different harmonic radiation. In the following, we analyze the filling-dependent doublon dy￾namics and their relation to the nonlinear current and HHG spectra. Figures 2(a-c) then show the filling-dependent cur￾r… view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2 view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3 view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4 view at source ↗
read the original abstract

High-harmonic generation (HHG) in strongly correlated Mott insulators is investigated using exact diagonalization and time-dependent density-matrix propagation of a laser-driven one-dimensional Hubbard chain. By projecting onto equilibrium Hubbard bands, we use the doublon population and its dynamics as a diagnostic to analyze intraband (spin-wave-like) and interband (doublon-holon creation) excitation channels. A filling-dependent crossover emerges: Bloch-like intraband response at dilute filling, mixed dynamics at intermediate filling, and interband-dominated HHG with plateau and cutoff near half filling. In the considered parameter range, increasing interaction strength $U$ strongly suppresses interband contributions through the enlarged Mott gap and correlation-induced localization. Intra- and interband current decomposition reveals opposing flows below the Mott gap (\Delta_{\mathrm{Mott}}) and selective dephasing suppression of interband coherence, enhancing net doublon accumulation. Time-frequency analysis uncovers the filling-dependent features of quantum trajectories, manifesting in distinct below-\Delta_{\mathrm{Mott}} emission. This doublon-based analysis provides a transparent link between equilibrium spin-charge separation and nonequilibrium strong-field response, and clarifies how dephasing modifies interband coherence and doublon accumulation.

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 / 3 minor

Summary. The manuscript investigates high-harmonic generation (HHG) in Mott insulators by simulating a laser-driven one-dimensional Hubbard chain using exact diagonalization and time-dependent density-matrix propagation. It introduces a doublon population diagnostic projected onto equilibrium Hubbard bands to distinguish intraband (spin-wave-like) and interband (doublon-holon) excitation channels. The central claims are a filling-dependent crossover in HHG response—from Bloch-like intraband at dilute fillings, through mixed dynamics, to interband-dominated with plateau and cutoff near half filling—and strong suppression of interband contributions with increasing interaction strength U due to the enlarged Mott gap and correlation-induced localization. Additional analyses include intra/interband current decomposition, dephasing effects, and time-frequency analysis linking to quantum trajectories and below-gap emission.

Significance. If the results hold, this work provides a valuable microscopic understanding of how electron correlations and filling control HHG in strongly correlated systems, offering a bridge between equilibrium spin-charge separation and nonequilibrium strong-field phenomena. The doublon-based analysis could inform experimental interpretations of HHG spectra in materials like transition metal oxides. The suppression with U highlights the role of correlations in controlling emission channels.

major comments (2)
  1. [Methods] The central claims of a filling-dependent crossover and U-induced suppression rest entirely on simulations of the 1D Hubbard chain (Methods section). The manuscript invokes 1D spin-charge separation to connect equilibrium bands to nonequilibrium doublon dynamics, but provides no quantitative discussion or test of how these features change in 2D or 3D lattices where spin-charge separation is absent; this is load-bearing for the generality asserted in the abstract.
  2. [Results] In the Results section, convergence checks with respect to chain length (number of sites) and time-step size for the time-dependent density-matrix propagation are not reported for the key observables (doublon population, intra/interband currents, and HHG spectra), particularly near half filling where finite-size effects and strong correlations are pronounced.
minor comments (3)
  1. [Abstract] The abstract refers to 'the considered parameter range' without listing explicit values of U, filling, laser frequency, or intensity; these should appear in a table or the first paragraph of the Model section for reproducibility.
  2. [Throughout] Notation for the Mott gap (denoted Δ_Mott) and related quantities is introduced without a dedicated definition paragraph; a consistent first-use definition would improve clarity.
  3. [Discussion] The time-frequency analysis of quantum trajectories is presented without reference to prior HHG trajectory literature in correlated systems; adding 1-2 citations would place the below-Δ_Mott emission features in context.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their careful reading and constructive comments on our manuscript. We address each major comment point by point below and have revised the manuscript where appropriate to improve clarity and rigor.

read point-by-point responses

  1. Referee: [Methods] The central claims of a filling-dependent crossover and U-induced suppression rest entirely on simulations of the 1D Hubbard chain (Methods section). The manuscript invokes 1D spin-charge separation to connect equilibrium bands to nonequilibrium doublon dynamics, but provides no quantitative discussion or test of how these features change in 2D or 3D lattices where spin-charge separation is absent; this is load-bearing for the generality asserted in the abstract.

    Authors: The manuscript explicitly studies the one-dimensional Hubbard model, as stated in the abstract and Methods, where spin-charge separation is exact and enables a direct mapping from equilibrium bands to doublon-holon dynamics. We do not provide quantitative tests in 2D or 3D because exact diagonalization does not scale to those dimensions. To address the concern, we have revised the abstract to stress the one-dimensional setting and added a paragraph in the Discussion clarifying that while the Mott-gap suppression mechanism is expected to be relevant more broadly, quantitative details may differ in higher dimensions due to the absence of spin-charge separation. Extensions would require other methods such as DMFT. revision: partial

  2. Referee: [Results] In the Results section, convergence checks with respect to chain length (number of sites) and time-step size for the time-dependent density-matrix propagation are not reported for the key observables (doublon population, intra/interband currents, and HHG spectra), particularly near half filling where finite-size effects and strong correlations are pronounced.

    Authors: We agree that explicit convergence checks strengthen the results. In the revised manuscript we have added a new subsection (and corresponding Supplementary Material) reporting convergence tests for the doublon population, intra- and interband currents, and HHG spectra with respect to chain length (N = 8–16 sites) and time-step size. These confirm that the filling-dependent crossover and U-suppression features remain qualitatively robust near half filling, although quantitative finite-size corrections are larger in the strongly correlated regime. revision: yes

standing simulated objections not resolved
  • Quantitative tests or detailed discussion of how the reported crossover and suppression change in 2D or 3D lattices, as these lie outside the reach of the exact-diagonalization methods used in the present work.

Circularity Check

0 steps flagged

No significant circularity in simulation-based derivation

full rationale

The paper computes HHG spectra and doublon-holon dynamics directly via exact diagonalization plus time-dependent density-matrix propagation on the 1D Hubbard chain. The filling-dependent crossover, intra/interband decomposition, and U-suppression of interband channels are numerical outputs of these simulations, not algebraic reductions or fitted parameters renamed as predictions. The doublon diagnostic is constructed from equilibrium bands but applied to driven trajectories without self-definition or tautological closure. No load-bearing self-citations, uniqueness theorems, or ansatzes imported from prior author work appear in the provided derivation chain. The analysis remains self-contained as forward simulation results.

Axiom & Free-Parameter Ledger

2 free parameters · 2 axioms · 0 invented entities

Central claim rests on the standard 1D Hubbard Hamiltonian, projection onto equilibrium bands, and the validity of doublon counting as a proxy for excitation channels; no new entities postulated.

free parameters (2)
  • interaction strength U
    Varied to show suppression of interband contributions; fitted or chosen to enlarge Mott gap.
  • electron filling
    Scanned across dilute to half-filling to reveal crossover behavior.
axioms (2)
  • domain assumption The 1D Hubbard model captures essential Mott physics under laser driving.
    Invoked throughout the simulation setup and interpretation of doublon dynamics.
  • domain assumption Projection onto equilibrium Hubbard bands remains valid during strong-field driving.
    Used to define doublon population diagnostic.

pith-pipeline@v0.9.0 · 5518 in / 1471 out tokens · 48979 ms · 2026-05-09T18:22:29.840120+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

43 extracted references · 43 canonical work pages

  1. [1]

    Ferray, A

    M. Ferray, A. L’Huillier, X. F. Li, L. A. Lompré, G. Mainfray, and C. Manus, Multiple-harmonic conver- sion of 1064 nm radiation in rare gases, J. Phys. B: At., Mol. Opt. Phys. 21, L31 (1988)

    work page 1988

  2. [2]

    P. B. Corkum, Plasma perspective on strong field multi- photon ionization, Phys. Rev. Lett. 71, 1994 (1993)

    work page 1994

  3. [3]

    Goulielmakis and T

    E. Goulielmakis and T. Brabec, High harmonic genera- tion in condensed matter, Nat. Photon. 16, 411 (2022)

    work page 2022

  4. [4]

    Heide, Y

    C. Heide, Y. Kobayashi, S. R. U. Haque, and S. Ghimire, Ultrafast high-harmonic spectroscopy of solids, Nat. Phys. 20, 1546–1557 (2024)

    work page 2024

  5. [5]

    Silva, I

    R. Silva, I. V. Blinov, A. N. Rubtsov, O. Smirnova, and M. Ivanov, High-harmonic spectroscopy of ultra- fast many-body dynamics in strongly correlated systems, Nat. Photon. 12, 266 (2018)

    work page 2018

  6. [6]

    A. J. Uzan-Narovlansky, Á. Jiménez-Galán, G. Oren- stein, R. E. Silva, T. Arusi-Parpar, S. Shames, B. D. Bruner, B. Yan, O. Smirnova, M. Ivanov, et al. , Obser- vation of light-driven band structure via multiband high- harmonic spectroscopy, Nat. Photon. 16, 428 (2022)

    work page 2022

  7. [7]

    Ghimire and D

    S. Ghimire and D. A. Reis, High-harmonic generation from solids, Nat. Phys. 15, 10 (2019)

    work page 2019

  8. [8]

    Zhang, Z

    J. Zhang, Z. Wang, F. Lengers, D. Wigger, D. E. Reiter, T. Kuhn, H. J. Wörner, and T. T. Luu, High-harmonic spectroscopy probes lattice dynamics, Nat. Photon. 18, 792 (2024)

    work page 2024

  9. [9]

    Vampa, C

    G. Vampa, C. R. McDonald, G. Orlando, D. D. Klug, P. B. Corkum, and T. Brabec, Theoretical analysis of high-harmonic generation in solids, Phys. Rev. Lett. 113, 073901 (2014)

    work page 2014

  10. [11]

    S. Imai, A. Ono, and S. Ishihara, High harmonic gener- ation in a correlated electron system, Phys. Rev. Lett. 124, 157404 (2020)

    work page 2020

  11. [12]

    Murakami, T

    Y. Murakami, T. Hansen, S. Takayoshi, L. B. Madsen, and P. Werner, Many-body effects on high-harmonic gen- eration in hubbard ladders, Phys. Rev. Lett. 134, 096504 (2025)

    work page 2025

  12. [13]

    C. Shao, H. Lu, X. Zhang, C. Yu, T. Tohyama, and R. Lu, High-harmonic generation approaching the quantum crit- ical point of strongly correlated systems, Phys. Rev. Lett. 128, 047401 (2022)

    work page 2022

  13. [14]

    C. S. Lange, T. Hansen, and L. B. Madsen, Electron- correlation-induced nonclassicality of light from high- order harmonic generation, Phys. Rev. A 109, 033110 (2024)

    work page 2024

  14. [15]

    Pizzi, A

    A. Pizzi, A. Gorlach, N. Rivera, A. Nunnenkamp, and I. Kaminer, Light emission from strongly driven many- body systems, Nat. Phys. 19, 551 (2023)

    work page 2023

  15. [16]

    Udono, K

    M. Udono, K. Sugimoto, T. Kaneko, and Y. Ohta, Exci- tonic effects on high-harmonic generation in mott insu- lators, Phys. Rev. B 105, L241108 (2022)

    work page 2022

  16. [17]

    Masur, D

    J. Masur, D. I. Bondar, and G. McCaul, Optical dis- tinguishability of mott insulators in the time versus fre- quency domain, Phys. Rev. A 106, 013110 (2022)

    work page 2022

  17. [18]

    C. Ma, C. Huang, Y. You, H. Liu, Z. Ding, M. Ding, C. Liu, J. Zhang, G. Li, Z. Sun, S. Wu, C. Ma, E. Wang, H. Hong, and K. Liu, Phase matching of high harmonic generation in twisted van der waals crystals, Phys. Rev. Lett. 135, 166902 (2025)

    work page 2025

  18. [19]

    Gabriele, C

    F. Gabriele, C. Ortix, M. Cuoco, and F. Forte, Hall- marks of spin textures for high-harmonic generation in two-dimensional materials, Phys. Rev. B 112, 155421 (2025)

    work page 2025

  19. [20]

    Murakami, M

    Y. Murakami, M. Eckstein, and P. Werner, High- harmonic generation in mott insulators, Phys. Rev. Lett. 121, 057405 (2018)

    work page 2018

  20. [21]

    Lysne, Y

    M. Lysne, Y. Murakami, and P. Werner, Signatures of bosonic excitations in high-harmonic spectra of mott in- sulators, Phys. Rev. B 101, 195139 (2020)

    work page 2020

  21. [22]

    Murakami, S

    Y. Murakami, S. Takayoshi, A. Koga, and P. Werner, High-harmonic generation in one-dimensional mott insu- lators, Phys. Rev. B 103, 035110 (2021)

    work page 2021

  22. [23]

    Murakami, D

    Y. Murakami, D. Golež, M. Eckstein, and P. Werner, Photoinduced nonequilibrium states in mott insulators, Rev. Mod. Phys. 97, 035001 (2025)

    work page 2025

  23. [24]

    Ikeda, Y

    R. Ikeda, Y. Murakami, D. Sakai, T. Miyamoto, T. Ito, and H. Okamoto, High harmonic generation reflecting the sub-cycle evolution of the mott transition under a 7 mid-infrared electric field (2025), arXiv:2508.00296

    work page arXiv 2025

  24. [25]

    Oka, Nonlinear doublon production in a mott insu- lator: Landau-dykhne method applied to an integrable model, Phys

    T. Oka, Nonlinear doublon production in a mott insu- lator: Landau-dykhne method applied to an integrable model, Phys. Rev. B 86 (2012)

    work page 2012

  25. [26]

    Queisser, K

    F. Queisser, K. Krutitsky, P. Navez, and R. Schützhold, Doublon–holon pair creation in mott–hubbard systems in analogy to qed, New J. Phys. 27, 023030 (2025)

    work page 2025

  26. [28]

    V. L. Chang, L. Yue, M. B. Gaarde, Y.-h. Chan, and D. Y. Qiu, Many-body enhancement of high-harmonic generation in monolayer MoS 2, Nat. Commun. 15, 6228 (2024)

    work page 2024

  27. [29]

    Uchida, G

    K. Uchida, G. Mattoni, S. Yonezawa, F. Nakamura, Y. Maeno, and K. Tanaka, High-order harmonic gen- eration and its unconventional scaling law in the mott- insulating ca2ruo4, Phys. Rev. Lett. 128, 127401 (2022)

    work page 2022

  28. [30]

    R. E. F. Silva, I. V. Blinov, A. N. Rubtsov, O. Smirnova, and M. Ivanov, High-harmonic spectroscopy of ultra- fast many-body dynamics in strongly correlated systems, Nat. Photon. 12, 266–270 (2018)

    work page 2018

  29. [31]

    Huang and T.-Y

    K. Huang and T.-Y. Du, Quasiparticle picture for high-harmonic generation in correlated electron systems, Phys. Rev. B 108, 155125 (2023)

    work page 2023

  30. [32]

    Ma, X.-B

    C. Ma, X.-B. Bian, and T.-Y. Du, Role of symme- try breaking in high-order harmonic generation from su-schrieffer-heeger systems, Phys. Rev. B 106, 125117 (2022)

    work page 2022

  31. [33]

    Murakami, K

    Y. Murakami, K. Uchida, A. Koga, K. Tanaka, and P. Werner, Anomalous temperature dependence of high- harmonic generation in mott insulators, Phys. Rev. Lett. 129, 157401 (2022)

    work page 2022

  32. [34]

    M. S. Mrudul, N. Tancogne-Dejean, A. Rubio, and G. Dixit, High-harmonic generation from spin-polarised defects in solids, npj Comput. Mater. 6, 10 (2020)

    work page 2020

  33. [35]

    F. H. L. Essler, H. Frahm, F. Göhmann, A. Klüm- per, and V. E. Korepin, The one-dimensional Hubbard model (Cambridge University Press, Cambridge, Eng- land, 2005)

    work page 2005

  34. [36]

    G. D. Mahan, Many-particle physics , 3rd ed. (Springer, New York, 2000)

    work page 2000

  35. [37]

    Hansen and L

    T. Hansen and L. B. Madsen, Doping effects in high- harmonic generation from correlated systems, Phys. Rev. B 106, 235142 (2022)

    work page 2022

  36. [38]

    Hansen, L

    T. Hansen, L. B. Madsen, and Y. Murakami, Charge creation via quantum tunneling in one-dimensional mott insulators: A numerical study of the extended hubbard model, Phys. Rev. B 112, 035144 (2025)

    work page 2025

  37. [39]

    Vampa, C

    G. Vampa, C. R. McDonald, G. Orlando, P. B. Corkum, and T. Brabec, Semiclassical analysis of high harmonic generation in bulk crystals, Phys. Rev. B 91, 064302 (2015)

    work page 2015

  38. [40]

    Wang, Y.-X

    G. Wang, Y.-X. Chen, M. Li, and T.-Y. Du, High-order harmonic spectroscopy for probing intraband and inter- band dynamics, Phys. Rev. A 109, 063511 (2024)

    work page 2024

  39. [41]

    Du, Control of high-order harmonic emission in solids via the tailored intraband current, Phys

    T.-Y. Du, Control of high-order harmonic emission in solids via the tailored intraband current, Phys. Rev. A 104, 063110 (2021)

    work page 2021

  40. [42]

    Zhang, Z

    X. Zhang, Z. Deng, L. Zhao, P. Wang, C. Wang, L. Zhou, G. Wang, T.-Y. Du, and M. Li, Anomalous in-plane anisotropic harmonic generation in few-layer pdse2, ACS Nano 19, 33557 (2025)

    work page 2025

  41. [43]

    Hogger, A

    W. Hogger, A. Riedel, D. Roy, A. Knothe, C. Gorini, J.-D. Urbina, and K. Richter, Coherent suppression and dephasing-induced reentrance of high harmonics in Dirac materials, arXiv 2503, 11432 (2025) , arXiv:2503.11432 [cond-mat.mes-hall]

    work page arXiv 2025

  42. [44]

    Orthodoxou, A

    C. Orthodoxou, A. Zaïr, and G. H. Booth, High har- monic generation in two-dimensional Mott insulators, npj Quantum Mater. 6, 76 (2021)

    work page 2021

  43. [45]

    M. Wu, S. Ghimire, D. A. Reis, K. J. Schafer, and M. B. Gaarde, High-harmonic generation from bloch electrons in solids, Phys. Rev. A 91, 043839 (2015)

    work page 2015