Pith. sign in

REVIEW 4 major objections 5 minor 56 references

Adding first-principles electron-phonon collision integrals to real-time TDDFT yields dissipative dynamics that reproduce measured lifetimes in silicon and WS2.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

T0 review · deepseek-v4-flash

2026-08-01 18:31 UTC pith:RSG2FEDM

load-bearing objection A credible and genuinely implemented e-ph collision-integral extension of rt-TDDFT, but the benchmark agreement is not a clean test because the Markovian limit is stretched on the very timescales used for validation. the 4 major comments →

arxiv 2607.17265 v2 pith:RSG2FEDM submitted 2026-07-19 cond-mat.mtrl-sci

First-principles electron-phonon scattering in real-time TDDFT

classification cond-mat.mtrl-sci
keywords electron-phonon scatteringreal-time TDDFTBorn-Markov approximationcollision integraldensity matrixcarrier relaxationintervalley scatteringtr-ARPES
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The paper sets out to remove a standing limitation of real-time time-dependent density functional theory: its unitary evolution cannot describe irreversible carrier relaxation, thermalization, or decoherence. It proposes to treat the electronic system as an open quantum system weakly coupled to a thermal phonon bath, deriving a collision integral for the reduced one-body density matrix from first-principles electron-phonon matrix elements within the Born-Markov approximation. The resulting equations keep the coherent real-space Kohn-Sham propagation intact while adding phonon-mediated transitions that redistribute carriers in energy and crystal momentum and damp coherences. If the framework works as claimed, it gives a practical first-principles route to simulate carrier cooling, intervalley scattering, and time-resolved spectroscopic signals in realistic materials, with validation shown for bulk silicon and monolayer WS2.

Core claim

The central claim is that equations (8) and (9) — coherent real-time Kohn-Sham propagation supplemented by a Born-Markov electron-phonon collision integral built from first-principles transition rates — correctly capture irreversible relaxation and decoherence in solids. The collision integral's diagonal part reduces to a Boltzmann gain-loss equation with Pauli blocking, while its off-diagonal part damps interband coherences using the same microscopic scattering channels. The paper validates this by matching simulated Si hot-carrier energy lifetimes (40 fs at 2.6 eV, 106 fs at 1.6 eV) to experimental values (30 fs and 120 fs), and simulated WS2 intervalley scattering (17 fs under linear pola

What carries the argument

The central object is the reduced single-particle density matrix rho(1) of the Kohn-Sham system, evolved by d rho(1)/dt = -i[H_KS, rho(1)] + I_e-ph[rho(1)]. The collision integral I_e-ph is the load-bearing piece: its diagonal entries are first-principles electron-phonon transition rates W(m,k)<-(n,k') that transfer carriers between bands and momenta with Pauli blocking, and its off-diagonal entries decay coherences with a rate gamma_n(k,t) built from the same W. In the implementation, coherent dynamics run on a real-space grid while the collision integral acts in a time-dependent adiabatic active space; the active-space density matrix is mapped back to natural orbitals with a phase- and rot

Load-bearing premise

The load-bearing premise is that the randomizing effect of lattice vibrations can be treated as instantaneous and memory-free, with the phonons staying in thermal equilibrium at 300 K, on the few-tens-of-femtoseconds timescales used to validate the method — a premise the paper itself flags as strained for the shortest momentum lifetimes.

What would settle it

A clean falsifier: measure the momentum-resolved population decay in photoexcited silicon under identical pump conditions but at a different lattice temperature; if the extracted sub-20 fs momentum lifetimes change substantially with bath temperature in a way the static-equilibrium-rate prediction cannot reproduce, the Markovian thermal-bath assumption is violated. Alternatively, an exact non-Markovian memory-kernel calculation on the same system that yields energy lifetimes more than a factor of two from the 40 fs and 106 fs values would undercut the validation.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • Conventional rt-TDDFT can be extended to open systems without abandoning real-space propagation, so relaxation and decoherence become first-principles outputs rather than phenomenological add-ons.
  • Because population relaxation and coherence decay arise from the same microscopic transition rates, pump-induced dephasing, damping of coherent oscillations, and carrier cooling can be described consistently in one simulation.
  • The validated energy lifetimes in silicon and intervalley lifetime in WS2 imply the method can predict hot-carrier cooling timescales relevant to photovoltaics and valleytronics.
  • Combined with a photoemission simulator, the framework turns internal carrier-population dynamics into computed time- and angle-resolved photoemission intensities, letting experiments and theory be compared directly.
  • The active-space construction keeps extra cost limited, making dissipation feasible for realistic materials and extendable to other reservoirs such as photons, magnons, or substrates.

Where Pith is reading between the lines

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

  • If the Born-Markov limit is genuinely pushed to its edge at about 10 fs, the cleanest next test is to extend the collision integral with a memory kernel and see whether the shortest momentum lifetimes change; until then, the sub-20 fs numbers should be read as trends rather than precise rates.
  • The same framework could be applied to solid-state high-harmonic generation, where decoherence competes with coherent driving; one testable prediction is that including electron-phonon scattering suppresses high harmonic yields at long pulse durations in a computable way.
  • The reported separation of fast momentum redistribution and slower energy relaxation in silicon suggests a two-stage hot-carrier picture that could be probed by transient absorption or two-photon photoemission with femtosecond resolution.
  • Because only equilibrium phonon rates enter, the method currently cannot describe hot-phonon feedback; a straightforward extension would couple the phonon occupations self-consistently, which should slow relaxation under strong excitation — a prediction testable against the pump-intensity dependence of valley lifetimes.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

4 major / 5 minor

Summary. The paper develops a dissipative real-time TDDFT framework in which the unitary Kohn-Sham propagation in real space is supplemented by a Born-Markov electron-phonon collision integral acting on a reduced single-particle density matrix in an adiabatic active space. The collision integral is built from first-principles electron-phonon matrix elements, coarse-grained onto a computationally tractable k-grid, and applied to the diagonal (Boltzmann-like) and off-diagonal (dephasing) sectors. The method is benchmarked on hot-carrier relaxation in bulk Si and intervalley scattering in monolayer WS2, and it is combined with a tSURFF-based tr-ARPES simulation to connect the microscopic dynamics to a photoemission observable. The reported quantitative agreements are Si energy lifetimes of 40/106 fs versus experimental 30/120 fs, and a WS2 intervalley lifetime of 17 fs versus an experimental exciton depolarization time of 16 fs.

Significance. If the central claims hold, the framework would be a valuable practical tool: it extends rt-TDDFT beyond unitary dynamics while retaining real-space nonlinear and spectroscopic capabilities, and it derives the dissipative kernel from first-principles electron-phonon interactions rather than from phenomenological parameters. The paper ships open code and input files, which is a notable strength for reproducibility. The method's main value lies in providing a single framework for population relaxation, coherence decay, and time-resolved observables. However, the validation is not yet as clean as the text suggests: the benchmark lifetimes are comparable to the phonon bath correlation time, the WS2 comparison is against an excitonic observable rather than free-carrier scattering, and the strong-field applicability is asserted rather than quantitatively bounded.

major comments (4)
  1. [Ultrafast carrier relaxation dynamics in silicon (Fig. 2d) and Ultrafast Valley Dynamics in WS2 (Fig. 3d)] The paper itself states that the extracted Si momentum lifetimes (8 and 14 fs) approach the Markovian validity limit set by the maximum Si phonon energy (~65 meV, ~10 fs bath correlation time). The WS2 intervalley lifetime of 17 fs is similarly only marginally longer than typical optical-phonon periods. Since the central validation rests on quantitative agreement with experiment for these short lifetimes, the agreement is not a clean test of the Born-Markov approximation. Please provide a quantitative estimate of the error incurred by the Markov and equilibrium-bath assumptions, e.g., by comparing the full memory integral in Eq. (S13) with the Markov result for the same rates, or by showing convergence of the extracted lifetimes with respect to the memory-time cutoff τ_E. Without such a check, the claim that the numerical agreement verifies the framework is overstated.
  2. [Ultrafast Valley Dynamics in WS2, comparison with Ref. [36]] The experiment quoted for the 16 fs lifetime probes exciton valley depolarization, while the simulation describes single-particle free-carrier intervalley scattering and explicitly omits long-range electron-hole interactions. The paper acknowledges this simplification but still presents the 17 fs versus 16 fs agreement as validation. This comparison is not apples-to-apples: the agreement could be coincidental or arise from cancellation of errors. Please either provide a theoretical argument for why the single-particle intervalley rate should match the exciton depolarization time, or reframe the WS2 result as a qualitative/semi-quantitative benchmark and temper the validation claim accordingly.
  3. [Eq. (9) and SI Sec. S.1] The main text states that Eq. (9) is the single-particle reduction of an underlying CPTP open-system evolution. The SI derives the general nonlinear quantum kinetic equation (S17) and then applies the secular approximation to obtain Eq. (S20) and the diagonal Boltzmann form. It is not demonstrated that the resulting diagonal + dephasing structure, with the state-dependent γ_n of Eq. (10), preserves complete positivity. The citation to Ref. [26] is not a proof. Since complete positivity is a central formal claim, either provide a direct proof for Eq. (9), state that positivity is only approximate and must be monitored, or explain how Ref. [26] covers this specific nonlinear secular form.
  4. [Coupling active-space dissipation to real-space propagation; Discussion] The transition rates W are computed from static equilibrium electron-phonon matrix elements, while the collision integral is applied in the instantaneous adiabatic basis of the driven system. The paper asserts that the adiabatic states remain continuously connected to the reference band states, but this is not quantified. For strong-field driving, the adiabatic Houston states are shifted in crystal momentum by the vector potential, so the e-ph matrix elements in the active space should in principle be evaluated at the instantaneous shifted momenta. The present benchmarks use moderate intensities, but the abstract and Discussion advertise the method for strong-field phenomena. Please specify the maximum field strength / vector-potential shift for which the static-rate approximation is valid, or describe how the rates are updated during the propagation.
minor comments (5)
  1. [Methods heading] The heading 'A veraged Transition Rate' contains a typo; it should be 'Averaged Transition Rate'.
  2. [Eq. (16) and Methods 'Basis Construction'] The residual states |d_μ(t)⟩ are constructed from the finite set of TDKS orbitals |ψ_μ(t)⟩, so Eq. (17) is a decomposition of the subspace spanned by the TDKS orbitals, not of the full single-particle Hilbert space. The statement 'we strictly preserve the Hilbert space completeness' is too strong and should be rephrased to refer to the completeness of the propagated orbital manifold.
  3. [Fig. 3 and text] The text reports a CP-light lifetime of 20 fs and an LP-light lifetime of 17 fs, but the fits are not shown in Fig. 3(c,d). Please indicate the fit windows and, ideally, the fit uncertainties, since the extracted lifetimes are central to the validation.
  4. [Eq. (14)] The convolution model for the tr-ARPES intensity includes parameters A1, A2, and τ0, but no details are given about how these are constrained or how sensitive τ0 is to the choice of pump/probe envelopes. A brief sensitivity statement would be useful.
  5. [Author contributions / Competing interests] The manuscript contains placeholders 'State individual contributions' and 'Declare competing interests or state none.' These should be completed before publication.

Circularity Check

0 steps flagged

No significant circularity: e-ph rates are first-principles inputs and the simulated lifetimes are read out, not fitted to the experimental targets.

full rationale

The derivation chain is self-contained. The collision integral of Eqs. (8)-(9) is built from first-principles electron-phonon matrix elements g (Eq. 7), and the transition rates W entering the collision integral are computed by Fermi's Golden Rule and coarse-graining in Eqs. (18)-(19) using EPW; they are not fitted parameters. The reported Si lifetimes (40 fs, 106 fs, 45 fs) and WS2 intervalley lifetimes (20 fs, 17 fs) are obtained by exponential fits to the simulated population dynamics and then compared with independent experimental values (30 fs, 120 fs, 16 fs). No equation or parameter in the method is defined in terms of those experimental lifetimes, so the validation is not circular by construction. The paper's own caveat that the extracted Si momentum lifetimes (8 fs and 14 fs) approach the Markovian validity limit is a limitation of the Born-Markov approximation in that regime, not a circular reduction of the prediction to its inputs. Citations to Octopus, tSURFF, and the customized EPW interface are implementation/code references, not load-bearing self-citations of an unverified result. No uniqueness theorem or ansatz is imported from the authors' prior work to force the central conclusion. Therefore no specific circular step can be exhibited.

Axiom & Free-Parameter Ledger

3 free parameters · 8 axioms · 0 invented entities

The framework is a standard open-quantum-system construction layered on DFT inputs; the main unstated or lightly tested elements are the active-space and coarse-graining choices, which are not converged, and the Markovian phonon-bath assumption.

free parameters (3)
  • Active-space manifold P (selected adiabatic bands and k points) = material-specific; not exhaustively quantified (subset of bands around the excitation window)
    The collision integral is evaluated only in the P-P, P-Q, Q-P blocks (Eq. 17, Methods); size and composition are user choices that determine which scattering channels can occur, but no convergence scan is reported.
  • Coarse-graining patch assignment C(k) = fine-to-coarse mapping: WS2 180x180x1 -> 30x30x1; Si 40^3 -> 20^3
    Effective rates W in Eq. (19) are algebraic averages over fine-grid patches; patch geometry and coarse-grid density affect energy/momentum conservation and rates, with no sensitivity analysis provided.
  • Exponential fit windows and tr-ARPES convolution parameters (A1, A2, tau0) = tau0 values: 40, 106, 45, 8, 14, 20, 17, 13 fs
    These are read-out parameters fitted to the simulated dynamics/spectra; they do not enter the equations of motion, but the central 'quantitative agreement' claim is defined in terms of them, and no uncertainties are given.
axioms (8)
  • domain assumption Born (weak-coupling) approximation: rho(t) approximately equals rho_S(t) tensor rho_E(t)
    Invoked in SI before Eq. (S6) to trace out the phonon bath; standard but limits the method to weak e-ph coupling.
  • domain assumption Markov approximation with infinite memory time and energy-conserving delta functions
    SI after Eq. (S14): density matrix and phonon occupations are taken constant over the bath memory time, and the upper integration limit is extended to infinity; the paper acknowledges the Si momentum lifetimes approach the validity boundary.
  • domain assumption Thermal phonon bath in equilibrium; phonon occupations locked to Bose-Einstein distribution (no hot-phonon feedback)
    SI Eq. (S14) and Discussion: n_q is locked to its initial value; the method excludes hot-phonon dynamics and feedback from the electronic system onto the bath.
  • domain assumption Secular approximation: non-secular coherence couplings average to zero
    SI Eq. (S20); required to reduce the full non-linear kinetic equation to the simple coherence-decay form in Eq. (9).
  • domain assumption Mean-field factorization of four-operator correlators
    SI Eq. (S10); closes the BBGKY hierarchy but discards genuine two-body electron-phonon correlations beyond mean field.
  • domain assumption Static equilibrium e-ph matrix elements and band states remain valid for rates during laser-driven dynamics
    Methods and Discussion: rates W are computed with EPW from equilibrium states and applied in the instantaneous adiabatic basis; assumes continuity and no strong-field modification of the scattering kernel.
  • domain assumption LDA exchange-correlation and HGH pseudopotentials give accurate band structure and e-ph matrix elements for Si and WS2
    Computational Details; LDA is known to underestimate band gaps (e.g., WS2 gap), which can shift scattering phase space.
  • domain assumption CPTP property of the reduced collision integral is cited to ref [26] rather than proven here
    Main text after Eq. (9) and SI: 'Both Eqs (S17,S21) preserve CPTP ... [26]'; the present paper does not provide an independent proof.

pith-pipeline@v1.3.0-alltime-deepseek · 17446 in / 19331 out tokens · 167192 ms · 2026-08-01T18:31:38.591168+00:00 · methodology

0 comments
read the original abstract

Real-time time-dependent density functional theory provides a first-principles description of coherent electron dynamics in laser-driven solids, but its unitary formulation cannot capture the irreversible scattering, relaxation, and decoherence processes that drive excited carriers toward equilibrium. Here, we develop a dissipative rt-TDDFT framework in which first-principles electron-phonon interactions enter the evolution of the reduced one-body density matrix through self-energy-derived collision integrals within the Born-Markov approximation. The approach retains the quantum-coherent real-time propagation of the electronic system while introducing phonon-mediated transitions that redistribute carriers in energy and crystal momentum, thereby incorporating the microscopic momentum-transfer processes responsible for relaxation in real materials. The resulting framework provides a practical first-principles route to simulate relaxation, decoherence, and time-resolved spectroscopic signatures in realistic crystalline materials.

Figures

Figures reproduced from arXiv: 2607.17265 by Alexander Buccheri, Hannes H\"ubener, Marti L\"uders, Shunsuke A. Sato, Subhojit Pal, Umberto De Giovannini, Zhengwei Nie.

Figure 1
Figure 1. Figure 1: Real-space/active-space propagation cycle. [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Ultrafast carrier relaxation dynamics in bulk silicon. a [PITH_FULL_IMAGE:figures/full_fig_p008_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Ultrafast valley scattering dynamics in monolayer WS [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: ARPES simulations of valley relaxation dynamics. a [PITH_FULL_IMAGE:figures/full_fig_p012_4.png] view at source ↗

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

56 extracted references · 2 canonical work pages

  1. [1]

    & Kohn, W

    Gross, E. & Kohn, W. Time-dependent density-functional theory. InAdvances in quantum chemistry, vol. 21, 255–291 (Elsevier, 1990)

  2. [2]

    Marques, M. A. & Gross, E. K. Time-dependent density functional theory.Annu. Rev. Phys. Chem.55, 427–455 (2004)

  3. [3]

    Schultze, M.et al.Attosecond band-gap dynamics in silicon.Science346, 1348 – 1352 (2014)

  4. [4]

    Physical Review Letters(2014)

    Wachter, G.et al.Ab-initio simulation of optical-field induced currents in dielectrics. Physical Review Letters(2014)

  5. [5]

    D., Kärtner, F

    Tancogne-Dejean, N., Mücke, O. D., Kärtner, F. X. & Rubio, A. Ellipticity dependence of high-harmonic generation in solids originating from coupled intraband and interband dynam- ics.Nature Communications8, 745 (2017). URL https://www.nature.com/articles/ s41467-017-00764-5

  6. [6]

    D., Kärtner, F

    Tancogne-Dejean, N., Mücke, O. D., Kärtner, F. X. & Rubio, A. Impact of the electronic band structure in high-harmonic generation spectra of solids.Physical Review Letters118, 087403 (2017). URLhttps://journals.aps.org/prl/abstract/10.1103/PhysRevLett. 118.087403

  7. [7]

    D., Hübener, H

    Giovannini, U. D., Hübener, H. & Rubio, A. Monitoring electron-photon dressing in WSe2. Nano Letters16, 7993–7998 (2016)

  8. [8]

    A., Giovannini, U

    Hübener, H., Sentef, M. A., Giovannini, U. D., Kemper, A. F. & Rubio, A. Creating stable floquet–weyl semimetals by laser-driving of 3d dirac materials.Nature Communications8, 13940 (2017). URLhttp://www.nature.com/articles/ncomms13940

  9. [9]

    Neufeld, O.et al.Time- and angle-resolved photoelectron spectroscopy of strong-field light-dressed solids: Prevalence of the adiabatic band picture.Physical Review Research4, 033101 (2022)

  10. [10]

    Fan, B.et al.Floquet optical selection rules in black phosphorus.Science Advances11, eadw2744 (2025). 5

  11. [11]

    Choi, D.et al.Observation of floquet–bloch states in monolayer graphene.Nature Physics 1–6 (2025)

  12. [12]

    Science353, 916 – 919 (2016)

    Lucchini, M.et al.Attosecond dynamical franz-keldysh effect in polycrystalline diamond. Science353, 916 – 919 (2016)

  13. [13]

    Lucchini, M.et al.Unravelling the intertwined atomic and bulk nature of localised excitons by attosecond spectroscopy.Nature Communications12, 1021 (2021)

  14. [14]

    Neb, S.et al.Local fields reveal atomic-scale nonadiabatic carrier-phonon dynamics.Science 391, 75–78 (2026)

  15. [15]

    Electron-phonon interactions from first principles.Reviews of Modern Physics 89, 015003 (2017)

    Giustino, F. Electron-phonon interactions from first principles.Reviews of Modern Physics 89, 015003 (2017)

  16. [16]

    Bernardi, M., Vigil-Fowler, D., Lischner, J., Neaton, J. B. & Louie, S. G. Ab initio study of hot carriers in the first picosecond after sunlight absorption in silicon.Physical review letters112, 257402 (2014)

  17. [17]

    & Burgdörfer, J

    Floss, I., Lemell, C., Yabana, K. & Burgdörfer, J. Incorporating decoherence into solid-state time-dependent density functional theory.Physical Review B99, 224301 (2019)

  18. [18]

    & Peralta, J

    Oz, A., Nitzan, A., Hod, O. & Peralta, J. E. Electron dynamics in open quantum systems: The driven liouville-von neumann methodology within time-dependent density functional theory.Journal of Chemical Theory and Computation19, 7496–7504 (2023)

  19. [19]

    Nature communications11, 2780 (2020)

    Xu, J.et al.Spin-phonon relaxation from a universal ab initio density-matrix approach. Nature communications11, 2780 (2020)

  20. [20]

    Marques, M. A. L., Castro, A., Bertsch, G. F. & Rubio, A. octopus: a first-principles tool for excited electron–ion dynamics.Computer Physics Communications151, 60 – 78 (2003)

  21. [21]

    Castro, A.et al.octopus: a tool for the application of time-dependent density functional theory.phys. stat. sol.(b)243, 2465 – 2488 (2006)

  22. [22]

    Andrade, X.et al.Time-dependent density-functional theory in massively parallel computer architectures: the octopus project.Journal of Physics: Condensed Matter24, 233202 (2012)

  23. [23]

    Andrade, X.et al.Real-space grids and the octopus code as tools for the development of new simulation approaches for electronic systems.Physical Chemistry Chemical Physics17, 31371 – 31396 (2015)

  24. [24]

    Tancogne-Dejean, N.et al.Octopus, a computational framework for exploring light-driven phenomena and quantum dynamics in extended and finite systems.The Journal of chemical physics152(2020)

  25. [25]

    D., Hübener, H

    Giovannini, U. D., Hübener, H. & Rubio, A. A first-principles time-dependent density functional theory framework for spin and time-resolved angular-resolved photoelectron spectroscopy in periodic systems.Journal of Chemical Theory and Computation13, 265 – 273 (2017). URLhttp://pubs.acs.org/doi/abs/10.1021/acs.jctc.6b00897

  26. [26]

    & Ping, Y

    Simoni, J., Riva, G. & Ping, Y. First-principles open quantum dynamics for solids based on density-matrix formalism.The Journal of Chemical Physics163(2025). 6

  27. [27]

    & Queisser, H

    Shockley, W. & Queisser, H. Detailed balance limit of efficiency of p–n junction solar cells. InRenewable energy, Vol2_35–Vol2_54 (Routledge, 2018)

  28. [28]

    Ross, R. T. & Nozik, A. J. Efficiency of hot-carrier solar energy converters.Journal of Applied Physics53, 3813–3818 (1982)

  29. [29]

    Marini, A. Competition between the electronic and phonon–mediated scattering channels in the out–of–equilibrium carrier dynamics of semiconductors: an ab-initio approach.Journal of Physics: Conference Series427, 012003 (2013)

  30. [30]

    Fischetti, M. V. & Laux, S. E. Monte carlo analysis of electron transport in small semicon- ductor devices including band-structure and space-charge effects.Physical Review B38, 9721 (1988)

  31. [31]

    & Zhao, J

    Wang, Z., Zheng, Z., Zheng, Q. & Zhao, J. Real-time ab initio investigation on hot electron relaxation dynamics in silicon.The Journal of Physical Chemistry Letters15, 3907–3913 (2024)

  32. [32]

    & Vast, N

    Tanimura, H., Kanasaki, J., Tanimura, K., Sjakste, J. & Vast, N. Ultrafast relaxation dynamics of highly excited hot electrons in silicon.Physical Review B100, 035201 (2019)

  33. [33]

    & Zhao, H

    Ceballos, F. & Zhao, H. Ultrafast laser spectroscopy of two-dimensional materials beyond graphene.Advanced Functional Materials27, 1604509 (2017)

  34. [34]

    Bertoni, R.et al.Generation and evolution of spin-, valley-, and layer-polarized excited carriers in inversion-symmetric wse 2.Physical review letters117, 277201 (2016)

  35. [35]

    Kumar, A.et al.Spin/valley coupled dynamics of electrons and holes at the mos2–mose2 interface.Nano Letters21, 7123–7130 (2021)

  36. [36]

    Wallauer, R.et al.Momentum-resolved observation of exciton formation dynamics in monolayer ws2.Nano letters21, 5867–5873 (2021)

  37. [37]

    V.et al.Sub-100 fs formation of dark excitons in monolayer ws2.Nano Letters24, 14663–14670 (2024)

    Kolesnichenko, P. V.et al.Sub-100 fs formation of dark excitons in monolayer ws2.Nano Letters24, 14663–14670 (2024)

  38. [38]

    Zhu, X.et al.A holistic view of the dynamics of long-lived valley polarized dark excitonic states in monolayer ws2.Nature Communications16, 6385 (2025)

  39. [39]

    & Yao, W

    Xiao, D., Liu, G.-B., Feng, W., Xu, X. & Yao, W. Coupled spin and valley physics in monolayers of mos 2 and other group-vi dichalcogenides.Physical review letters108, 196802 (2012)

  40. [40]

    M.et al.Optical generation of excitonic valley coherence in monolayer wse2

    Jones, A. M.et al.Optical generation of excitonic valley coherence in monolayer wse2. Nature nanotechnology8, 634–638 (2013)

  41. [41]

    Timmer, D.et al.Ultrafast coherent exciton couplings and many-body interactions in monolayer ws2.Nano Letters24, 8117–8125 (2024)

  42. [42]

    & Rubio, A

    De Giovannini, U., Hübener, H. & Rubio, A. A first-principles time-dependent density functional theory framework for spin and time-resolved angular-resolved photoelectron spectroscopy in periodic systems.Journal of chemical theory and computation13, 265–273 (2017). 7

  43. [43]

    De Giovannini, U., Larsen, A. H. & Rubio, A. Modeling electron dynamics coupled to continuum states in finite volumes with absorbing boundaries.The European Physical Journal B88, 56 (2015)

  44. [44]

    H., Larsen, D

    Van Stokkum, I. H., Larsen, D. S. & Van Grondelle, R. Global and target analysis of time-resolved spectra.Biochimica et Biophysica Acta (BBA)-Bioenergetics1657, 82–104 (2004)

  45. [45]

    Acceleration of electrons in a crystal lattice.Physical Review57, 184 (1940)

    Houston, W. Acceleration of electrons in a crystal lattice.Physical Review57, 184 (1940)

  46. [46]

    Sato, S. A. & Yabana, K. Efficient basis expansion for describing linear and nonlinear electron dynamics in crystalline solids.Physical Review B89, 224305 (2014). URL http://journals.aps.org/prb/abstract/10.1103/PhysRevB.89.224305

  47. [47]

    The averaged transition rates for the real-time dissipative electronic dyanmics for octopus code from epw.https://gitlab.com/supal/epw-octopus-interface(2026)

  48. [48]

    Optimally smooth norm-conserving pseudopotentials.Physical Review B32, 8412 (1985)

    Vanderbilt, D. Optimally smooth norm-conserving pseudopotentials.Physical Review B32, 8412 (1985)

  49. [49]

    Perdew, J. P. & Zunger, A. Self-interaction correction to density-functional approximations for many-electron systems.Physical review B23, 5048 (1981)

  50. [50]

    & Jellinek, F

    Schutte, W., De Boer, J. & Jellinek, F. Crystal structures of tungsten disulfide and diselenide. Journal of Solid State Chemistry70, 207–209 (1987)

  51. [51]

    Monkhorst, H. J. & Pack, J. D. Special points for brillouin-zone integrations.Physical review B13, 5188 (1976)

  52. [52]

    Journal of physics: Condensed matter29, 465901 (2017)

    Giannozzi, P.et al.Advanced capabilities for materials modelling with quantum espresso. Journal of physics: Condensed matter29, 465901 (2017)

  53. [53]

    & Van Leeuwen, R.Nonequilibrium many-body theory of quantum systems: a modern introduction(Cambridge University Press, 2013)

    Stefanucci, G. & Van Leeuwen, R.Nonequilibrium many-body theory of quantum systems: a modern introduction(Cambridge University Press, 2013)

  54. [54]

    Haug, H., Koch, S. W. & Keldysh, L. V. Quantum theory of the optical and electronic properties of semiconductors (1994)

  55. [55]

    & Jauho, A.-P.Quantum kinetics in transport and optics of semiconductors (Springer, 2008)

    Haug, H. & Jauho, A.-P.Quantum kinetics in transport and optics of semiconductors (Springer, 2008)

  56. [56]

    Korolev, V.et al.Unveiling the role of electron-phonon scattering in dephasing high-order harmonics in solids.arXiv preprint arXiv:2401.12929(2024). 8