REVIEW 5 major objections 7 minor 33 references
Role of phonon coupling in driving photo-excited Mott insulators towards a transient superconducting steady state
T0 review · 5 major / 7 minor · reviewed 2026-08-11 · deepseek-v4-flash
Pith's one-line read A laser-pumped Mott insulator coupled to a cold phonon bath settles into a long-lived prethermalized state with persistent eta-pairing superconducting order, and the nonequilibrium steady-state formalism reproduces its spectral functions.
desk verdict Plausible and well-motivated, but the central long-time eta-SC plateau rests on an unchecked memory-truncation cutoff; needs convergence evidence before I'd trust it. 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 machinery is memory-truncated Kadanoff-Baym nonequilibrium DMFT on the Bethe lattice, with entropy-cooling pulses to prepare a photodoped state and a Holstein phonon bath that acts as an energy sink. The central object is the eta order parameter $\eta_x = \frac{1}{2}(c^\dagger_\uparrow c^\dagger_\downarrow + \mathrm{H.c.})$ with alternating sign between sublattices, measured in the Nambu Green's function; its plateau signals the prethermal superconducting state. The NESS counterpart replaces the laser and phonon history by permanent weak coupling to cold fermion baths at $\pm U/2$ plus the same phonon bath, so that the two-time Green's function depends only on the time difference. The load-bearing connection is that both routes give the same spectral functions when the effective temperature and doublon density match.
What would settle it
Recompute the $U=10$, $d\approx0.4$ case with progressively larger truncation windows $t_c$ (for example, doubling it twice) and check that the $\eta_x$ plateau and the extracted $\beta\approx16$ do not shift, and repeat the run with an exact impurity solver to confirm the spectral functions still match the NESS result.
Extended reading notes
Core claim
The central claim of the paper is that in the single-orbital repulsive Hubbard model on the Bethe lattice, after a chirped laser pulse creates roughly $d\approx 0.4$ doublons, the eta-pairing order parameter $\eta_x$ (a staggered superconducting order in which electron pairs have alternating phase between sublattices) grows during the pulse and then decays if the system is isolated; with a Holstein phonon coupling ($\omega_0=0.4$, $g=0.2$) and $U \gtrsim 10$, the order parameter instead rises to $\eta_x \approx 0.2$–$0.3$ and stays there until the end of the simulation at $t=2000$ (in units of inverse hopping). The system reaches a quasi-steady state with an almost constant spectral function and an effective inverse temperature $\beta \approx 16$ extracted identically from the normal and anomalous components. Nonequilibrium steady-state (NESS) DMFT with cold fermion baths, tuned to the same doping and temperature, reproduces the total and occupied spectral functions from the real-time simulation, with closer agreement at lower photodoping $d\approx 0.14$ than at $d\approx 0.4$. This establishes that the NESS construction is a valid description of the long-lived prethermalized eta-paired state in large-gap photodoped Mott insulators.
Load-bearing premise
The long-time results assume that all memory of the initial state dies out within the self-energy truncation window $t_c$, and that the non-crossing impurity solver is accurate; if correlations persist beyond $t_c$, the apparent superconducting plateau could be a numerical artifact.
Editorial extensions
If this is right
- For $U \gtrsim 10$ and photodoping $d \sim 0.4$, the $\eta_x$ order parameter saturates near $0.2$–$0.3$ and persists beyond $t=2000$ inverse hoppings, a timescale of several picoseconds in real materials.
- The NESS formalism reproduces the real-time spectral functions, so steady-state calculations can be used in place of expensive real-time simulations to hunt for hidden orders.
- Phonon coupling plays two roles: it cools the photo-carriers in large-gap systems, but in small-gap systems ($U \lesssim 8$) it accelerates doublon-holon recombination and prevents $\eta$ order from forming.
- Without phonon coupling, the photodoped system slowly heats toward a negative-temperature state; with phonon coupling, that heating is suppressed on the simulated timescales, leaving a well-defined effective temperature.
Reading between the lines
- A natural follow-up is a systematic convergence study in the memory-truncation window $t_c$, which would confirm that the $\eta_x$ plateau is a physical prethermal state rather than a truncation artefact.
- The same steady-state-plus-phonon protocol could be applied to multi-orbital or frustrated Hubbard models to search for other photo-induced orders, such as spin-triplet or chiral pairing, by independently controlling doping and effective temperature.
- An experimental test could look for a long-lived plateau in the transient optical or terahertz response of a large-gap Mott material with strong electron-phonon coupling; the predicted lifetime of thousands of inverse hoppings corresponds to picoseconds.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The authors study the long-time dynamics of photo-doped Mott insulators using nonequilibrium DMFT with memory-truncated Kadanoff-Baym equations and the NCA impurity solver, after preparing photodoped states via entropy cooling with chirped pulses. They show that without phonon coupling the η-pairing order decays due to heating, while coupling to a cold phonon bath can cool the carriers and, for U≳10 with doping d∼0.4, stabilize ηx at a plateau that persists beyond t=2000 (in units of inverse hopping). They further argue that the resulting quasi-steady state is well reproduced by a steady-state (NESS) DMFT calculation with suitably adjusted fermion-bath parameters.
Significance. If the central claim holds, the paper establishes a concrete mechanism — dissipative phonon cooling — by which light-induced η-paired superconducting order can be made long-lived in large-gap Mott insulators, and it demonstrates that the cheaper NESS formalism can capture the relevant quasi-steady state. The paper's strengths are its systematic parameter scans in U and photodoping, the two-time spectral-function analysis including the anomalous component, and the use of a standard, documented numerical framework. The main numerical evidence is nevertheless missing a set of convergence and validation checks that are needed before the long-time plateau can be regarded as established.
major comments (5)
- [§II and §III B] The long-time plateau in Fig. 5(b) is computed within the memory-truncation scheme described in Sec. II, but no value of the cutoff tc is reported anywhere and no convergence test with respect to tc is shown. This is not a merely numerical detail: in Eq. (7) the phonon contribution contains the free-boson propagator D(t,t′), which oscillates without decaying, so the decay of the memory kernel relies entirely on the decay of G(t,t′); in the η-SC state the anomalous component of G is precisely the long-lived quantity whose plateau is the central result. If tc is shorter than the decay time of that anomalous component, the truncation can artificially force the appearance of a steady state. Please report the value of tc used in Figs. 5(b) and 7 and provide a convergence study of ηx(t) and the spectral functions as tc is increased.
- [§II, §III A, and §III D] The manuscript states in Sec. II that 'In all our calculations we take ω0 = 0.4', but Sec. III A and Fig. 2 describe the phonon coupling as having frequency ω0 = 0.2; Sec. III D again uses ω0 = 0.4. Because g^2/ω0 changes from 0.1 to 0.2 between these values and the phonon cooling rate enters the central stabilization argument, this inconsistency must be resolved by specifying the value used in each figure or correcting the typo.
- [§III D] The comparison in Fig. 7 is presented as evidence that the NESS approach 'well describes' the real-time quasi-steady state, but the fermion-bath coupling Γ and temperature Tb are adjusted so that the NESS state matches the effective temperature and doublon density extracted from the real-time simulation. The agreement is therefore a consistency check rather than an independent confirmation. To make the NESS claim meaningful, the authors should report the adjusted bath parameters, test the sensitivity of the spectra to these parameters, and ideally predict at least one observable (e.g. the anomalous spectral function or the order parameter) without fitting.
- [§II, Eq. (5)] The order parameter is measured with a seed field Pseed = 0.001, and Sec. II asserts that in a symmetry-broken phase the order parameter becomes independent of the seed. No such independence is demonstrated. A seed-field dependence study (e.g. Pseed = 0.0001 and 0.01) is required to rule out that the plateau in Fig. 5(b) is a driven response rather than spontaneous η-SC order.
- [§II (impurity solver)] The noncrossing approximation (NCA) is used for both real-time and NESS calculations, but no accuracy assessment is provided for the parameter regime U ≈ 10, d ≈ 0.4, g = 0.2. NCA is an uncontrolled approximation, and the quantitative values of ηx and of the decay rates in Figs. 3 and 5 may depend on it. Please add a comparison with an alternative solver (e.g. IPT or a numerically exact method for one representative case) or at least a discussion of the expected NCA accuracy in this regime.
minor comments (7)
- [§III C] The section title 'Effect fo η-SC order' should read 'Effect of η-SC order'.
- [§II] The text 'coordination numberz' should read 'coordination number z'.
- [§III A] The text says 'Ωfin > Ωin' but the variable introduced is Ωini; please fix the typo.
- [§II, Eq. (3)] The next-nearest-neighbor hopping vNN_h is introduced through the self-consistency equation, but no corresponding term is added to the Hamiltonian in Eq. (1); please clarify whether vNN_h is part of the model Hamiltonian or only an effective hopping introduced in the DMFT self-consistency.
- [Fig. 2(d)] The arrow indicating the fully thermalized kinetic energy should be defined in the caption (e.g. the value from a thermal equilibrium calculation at β = 30).
- [Fig. 7 caption] The abbreviation 'KB contour' is used without definition; spell out 'Kadanoff-Baym contour' at first use.
- [§III A] The word 'consitent' should be 'consistent'.
Circularity Check
Partial circularity only in the NESS comparison: bath parameters are tuned to reproduce the real-time effective temperature and doublon density before the spectra are compared; the central real-time eta-SC plateau is otherwise independent.
-
fitted input called prediction
[Sec. III D, Fig. 7; cf. Sec. II definition of effective temperature]
"By adjusting the coupling to the Fermion bath Γ and the temperature of the Fermion baths Tb, we are able to create photodoped steady states with an effective temperature very close to that obtained in the quasi-steady state of the real-time simulations. Figure 7 compares these spectral functions... We find a close correspondence between the NESS steady state result and the real-time simulation result for both A(ω) and A<(ω), as shown in Fig. 7."
The real-time effective temperature is itself defined by fitting a Fermi distribution to the ratio A<(ω)/A(ω) inside the Hubbard bands. The NESS calculation is then run with Γ and Tb adjusted until it has an effective temperature and doublon density close to those of the real-time state. Because the comparison in Fig. 7 includes A<(ω), the occupied spectral weight is matched by construction through the chosen bath temperature and coupling, rather than independently predicted. The agreement is therefore a consistency check for the total spectral function A(ω) and for the emergent η-SC order, but part of the claimed 'close correspondence' is enforced by the parameter adjustment.
full rationale
The central physical claim—that coupling to a cold phonon bath produces a long-lived quasi-steady η-paired state in large-gap Mott insulators—comes from real-time memory-truncated DMFT simulations (Figs. 2–5) and is not fitted to any target result. No parameter is adjusted to force the plateau in ηx(t); the order grows from a small seed and then decays or stabilizes according to U and photodoping. The NESS comparison in Sec. III D is the only partially circular element: its bath parameters are explicitly adjusted to reproduce the real-time effective temperature and doublon density, so the close match in the occupied spectra is partly by construction. This is a consistency check with independent content in the total spectral line shape and in the ordered state, so the score is moderate rather than high. Self-citations to entropy cooling, NESS, and memory-truncation methods are methodological, and reference [25] provides an independent numerically exact benchmark, so they are not load-bearing circularity. Two non-circular but important gaps remain: the memory-truncation cutoff tc is never reported or convergence-tested (Sec. II states only that 'the self energy is assumed to decay within a cutoff window tc'), and the phonon frequency is given as ω0 = 0.4 in Sec. II and Sec. III D but as ω0 = 0.2 in Sec. III A. These gaps weaken confidence in the long-time plateau and reproducibility, but they are not circularity, and they are weighed here mainly as caveats on the central numerical result.
Assumptions & free parameters
free parameters (6)
- next_nearest_neighbor_hopping_ratio =
v_NN/v_h = 0.25
- phonon_frequency =
omega_0 = 0.4 (stated as 0.2 in Sec. III A)
- electron_phonon_coupling =
g = 0.2, with g^2/omega_0 < 1
- seed_field =
P_seed = 0.001
- NESS_bath_parameters =
Gamma and T_b adjusted to match real-time beta_eff and d
- chirped_pulse_parameters =
not specified (A0, envelope, tramp, Omega_fin per curve)
assumptions (5)
- domain assumption DMFT self-consistency on the infinite-coordination Bethe lattice
- domain assumption NCA impurity solver accuracy
- ad hoc to paper Memory truncation validity
- domain assumption Cold phonon bath with fixed temperature
- domain assumption Migdal-type electron-phonon self-energy
Cite this review
Pith. "Pith review of Role of phonon coupling in driving photo-excited Mott insulators towards a transient superconducting steady state." pith.science (2026). https://pith.science/paper/XEX4EV3M
@misc{pith2026241219205,
author = {Pith},
title = {Pith review of: Role of phonon coupling in driving photo-excited Mott insulators towards a transient superconducting steady state},
year = {2026},
howpublished = {\url{https://pith.science/paper/XEX4EV3M}},
note = {Machine review of arXiv:2412.19205}
}
abstract
Understanding light-induced hidden orders is relevant for nonequilibrium materials control and future ultrafast technologies. Hidden superconducting order, in particular, has been a focus of recent experimental and theoretical efforts. In this study, we investigate the stability of light-induced $\eta$ pairing. Using a memory truncated implementation of nonequilibrium dynamical mean field theory (DMFT) and entropy cooling techniques, we study the long-time dynamics of the photoinduced superconducting state. In the presence of coupling to a cold phonon bath, the photodoped system reaches a quasi-steady state, which can be sustained over a long period of time in large-gap Mott insulators. We show that this long-lived prethermalized state is well described by the nonequilibrium steady state implementation of DMFT.
Figures
Figures from the paper (4 more)
Reference graph
Works this paper leans on
-
[1]
Cavalleri, Photo-induced superconductivity, Contemp
A. Cavalleri, Photo-induced superconductivity, Contemp. Phys. 59, 31 (2018)
work page 2018
- [2]
-
[3]
T. Suzuki, T. Someya, T. Hashimoto, S. Michimae, M. Watan- abe, M. Fujisawa, T. Kanai, N. Ishii, J. Itatani, S. Kasahara, Y . Matsuda, T. Shibauchi, K. Okazaki and S. Shin, Photoinduced possible superconducting state with long-lived disproportionate band filling in FeSe. Commun. Phys. 2, 115 (2019)
work page 2019
- [4]
-
[5]
M. Mitrano, A. Cantaluppi, D. Nicoletti, S. Kaiser, A. Perucchi, S. Lupi, P. Di Pietro, D. Pontiroli, M. Ricco, S. R. Clark, et al., Possible light-induced superconductivity in K 3C60 at high temperature, Nature 530, 461 (2016)
work page 2016
- [6]
-
[7]
Y . Murakami, D. Golez, M. Eckstein, P.Werner, Photo-induced nonequilibrium states in Mott insulators, arXiv:2310.05201
-
[8]
N. Strohmaier, D. Greif, R. J ¨ordens, L. Tarruell, H. Moritz, T. Esslinger, R. Sensarma, D. Pekker, E. Altman and E. Demler, Observation of Elastic Doublon Decay in the Fermi-Hubbard Model, Phys. Rev. Lett. 104, 080401 (2010)
work page 2010
Show all 33 references
-
[9]
Sensarma, D
R. Sensarma, D. Pekker, E. Altman, E. Demler, N. Strohmaier, D. Greif, R. J ¨ordens, L. Tarruell, H. Moritz and T. Esslinger, Lifetime of double occupancies in the Fermi-Hubbard model, Phys. Rev. B 82, 224302 (2010)
2010
-
[10]
Eckstein and P
M. Eckstein and P. Werner, Thermalization of a pump-excited Mott insulator, Phys. Rev. B 84, 035122 (2011)
2011
-
[11]
Lenar ˇciˇc and P
Z. Lenar ˇciˇc and P. Prelov ˇsek, Ultrafast charge recombination in a photoexcited mott-hubbard insulator. Phys. Rev. Lett. 111, 016401 (2013)
2013
-
[12]
Lenarcic and P
Z. Lenarcic and P. Prelovsec, Charge recombination in undoped cuprates, Phys. Rev. B 90, 235136 (2014)
2014
-
[13]
Murakami, S
Y . Murakami, S. Takayoshi, T. Kaneko, A. M. L ¨auchli, and P. Werner, Spin, Charge, and η-Spin Separation in One- Dimensional Photodoped Mott Insulators, Phys. Rev. Lett.130, 106501 (2023)
2023
-
[14]
Murakami, S
Y . Murakami, S. Takayoshi, T. Kaneko, Z. Sun, D. Gole ˚A ¨Y, A. J. Millis, and P. Werner, Exploring nonequilibrium phases of photo-doped Mott insulators with generalized Gibbs ensembles, Comm. Phys. 5, 23 (2022)
2022
-
[15]
J. Li, D. Golez, P. Werner, and M. Eckstein, η-paired supercon- ducting hidden phase in photodoped Mott insulators, Phys. Rev. B 102, 165136 (2020)
2020
-
[16]
S. Ray, Y . Murakami, and P. Werner, Nonthermal superconduc- tivity in photodoped multiorbital Hubbard systems, Phys. Rev. B 108, 174515 (2023)
2023
-
[17]
Ray and P
S. Ray and P. Werner, Photoinduced ferromagnetic and super- conducting orders in multiorbital Hubbard models, Phys. Rev. B 110, L041109 (2024)
2024
-
[18]
Ejima, T
S. Ejima, T. Kaneko, F. Lange, S. Yunoki, and H. Fehske, Pho- toinduced η-pairing at finite temperatures, Phys. Rev. Research 2, 032008(R) (2020)
2020
-
[19]
Kaneko, S
T. Kaneko, S. Yunoki, and A. J. Millis, Charge stiffness and long-range correlation in the optically induced η-pairing state of the one-dimensional Hubbard model, Phys. Rev. Research2, 032027(R) (2020)
2020
-
[20]
J. Li, M. M ¨uller, A. Kim, A. L ¨aeuchli, and P. Werner, Twisted chiral superconductivity in photodoped frustrated Mott insula- tors, Phys. Rev. B 107, 205115 (2023)
2023
-
[21]
Sch ¨uler, D
M. Sch ¨uler, D. Golelez, Y . Murakami, N. Bittner, A. Herrmann, H. U. R. Strand, P. Werner, and M. Eckstein, NESSi: The Non- Equilibrium Systems Simulation package, Comput. Phys. Com- mun. 257, 107484 (2020)
2020
-
[22]
Werner, M
P. Werner, M. Eckstein, M. M ¨uller, and G. Refael, Light- induced evaporative cooling of holes in the Hubbard model, Nat. Comm. 10, 5556 (2019)
2019
-
[23]
Werner, J
P. Werner, J. Li, D. Golez, and M. Eckstein, Entropy-cooled nonequilibrium states of the Hubbard model, Phys. Rev. B100, 155130 (2019). 9
2019
-
[24]
Li and M
J. Li and M. Eckstein, Nonequilibrium steady-state theory of photodoped Mott insulators, Phys. Rev. B 103, 045133 (2021)
2021
-
[25]
K ¨unzel, A
F. K ¨unzel, A. Erpenbeck, D. Werner, E. Arrigoni, E. Gull, G. Cohen, and M. Eckstein, Numerically Exact Simulation of Photodoped Mott Insulators, Phys. Rev. Lett. 132, 176501 (2024)
2024
-
[26]
Sch ¨uler, M
M. Sch ¨uler, M. Eckstein, and P. Werner, Truncating the mem- ory time in nonequilibrium dynamical mean field theory calcu- lations, Phys. Rev. B 97, 245129 (2018)
2018
-
[27]
Stahl, N
C. Stahl, N. Dasari, J. Li, A. Picano, P. Werner, and M. Eck- stein, Memory truncated Kadanoff-Baym equations, Phys. Rev. B 105, 115146 (2022)
2022
-
[28]
Georges, G
A. Georges, G. Kotliar, W. Krauth, and M. J. Rozenberg, Dy- namical mean-field theory of strongly correlated fermion sys- tems and the limit of infinite dimensions, Rev. Mod. Phys. 68, 13 (1996)
1996
-
[29]
Eckstein and P
M. Eckstein and P. Werner, Nonequilibrium dynamical mean- field calculations based on the noncrossing approximation and its generalizations, Phys. Rev. B 82, 115115 (2010)
2010
-
[30]
Eckstein and P
M. Eckstein and P. Werner, Photoinduced States in a Mott Insu- lator, Phys. Rev. Lett. 110, 126401 (2013)
2013
-
[31]
Thoenniss, M
J. Thoenniss, M. Sonner, A. Lerose, and D. A. Abanin, Effi- cient method for quantum impurity problems out of equilib- rium, Phys. Rev. B 107, L201115 (2023)
2023
-
[32]
Eckstein, Solving quantum impurity models in the non- equilibrium steady state with tensor trains, arXiv:2410.19707 (2024)
M. Eckstein, Solving quantum impurity models in the non- equilibrium steady state with tensor trains, arXiv:2410.19707 (2024)
2024 arXiv
-
[33]
Kim and P
A. Kim and P. Werner, Strong coupling impurity solver based on quantics tensor cross interpolation, arXiv:2411.19026 (2024)
2024 arXiv
Reviewed August 11, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.