Full Quantum and Mixed Quantum--Classical Dynamics of Hot Exciton Cooling in Semiconductor Nanocrystals

Bokang Hou; Eran Rabani; Hans-Dieter Meyer; Johan E. Runeson; Michael Thoss; Salvatore Gatto; Samuel L. Rudge

arxiv: 2605.27708 · v1 · pith:SIHG6BTUnew · submitted 2026-05-26 · ⚛️ physics.chem-ph · cond-mat.mes-hall· cond-mat.mtrl-sci

Full Quantum and Mixed Quantum--Classical Dynamics of Hot Exciton Cooling in Semiconductor Nanocrystals

Bokang Hou , Johan E. Runeson , Samuel L. Rudge , Salvatore Gatto , Hans-Dieter Meyer , Michael Thoss , Eran Rabani This is my paper

Pith reviewed 2026-06-29 14:26 UTC · model grok-4.3

classification ⚛️ physics.chem-ph cond-mat.mes-hallcond-mat.mtrl-sci

keywords hot exciton coolingsemiconductor nanocrystalsCdSeexciton-phonon couplingsurface hoppingquantum master equationmixed quantum-classical dynamicsdiabatic state mixing

0 comments

The pith

Mapping approach to surface hopping matches full quantum dynamics for hot exciton cooling in CdSe nanocrystals.

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

The paper benchmarks perturbative quantum master equations and mixed quantum-classical methods against exact quantum dynamics for hot-exciton relaxation in atomistically parameterized CdSe and CdSe/CdS nanocrystal models. It establishes that bare CdSe shows an ultrafast initial decay from diabatic state mixing driven by low-frequency phonon fluctuations, followed by slower cooling, while core-shell systems are dominated by the slower process. The quantum master equation reproduces the fast decay but often fails on slower relaxation in the diabatic picture. The mapping approach to surface hopping yields the closest match to both the benchmark trajectories and the correct equilibrium populations.

Core claim

Using fully quantum mechanical reference calculations on CdSe core and CdSe/CdS core-shell models, the work shows that the ultrafast component of exciton cooling arises from rapid diabatic state mixing induced by thermal fluctuations of low-frequency phonons rather than from nuclear-assisted energy relaxation. The perturbative quantum master equation captures this initial fast decay but can fail for the subsequent slower relaxation when formulated in the diabatic representation. Among the tested approximate methods, the mapping approach to surface hopping produces the most consistent agreement with both the reference dynamics and the equilibrium state populations.

What carries the argument

The mapping approach to surface hopping (MASH), which approximates the quantum exciton-phonon dynamics through classical trajectories with stochastic surface hops.

Load-bearing premise

The atomistically parameterized exciton-phonon Hamiltonians for the CdSe and core-shell models are accurate enough to separate diabatic mixing from energy relaxation and to serve as a reliable benchmark.

What would settle it

A higher-accuracy full-quantum calculation or direct experimental measurement on the identical CdSe model Hamiltonian that produced decay rates or final populations differing from MASH predictions would falsify the claim of consistent agreement.

Figures

Figures reproduced from arXiv: 2605.27708 by Bokang Hou, Eran Rabani, Hans-Dieter Meyer, Johan E. Runeson, Michael Thoss, Salvatore Gatto, Samuel L. Rudge.

**Figure 2.** Figure 2: Real-time diabatic population dynamics P (d) n (t) for CdSe core and CdSe/CdS core–shell NCs for the highest n = 9 and lowest n = 1 excitonic states. Thermofield ML-MCTDH method is used as the benchmark. Three classes of approximate methods are compared: perturbative QME, mean-field (Ehrenfest, spin-LSC), and surface hopping (MASH and FSSH) approaches. The black dotted line is the diabatic equilibrium pop… view at source ↗

**Figure 3.** Figure 3: State-resolved diabatic equilibrium populations for (a) CdSe core and (b) CdSe/CdS core–shell NCs at [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗

**Figure 4.** Figure 4: Results for a reduced two-level exciton–phonon model. (a) Population dynamics from ML-MCTDH, MASH, QME [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗

**Figure 5.** Figure 5: Basis-independent excitonic energy relaxation in (a) CdSe core and (b) CdSe/CdS core–shell NCs. The plotted [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗

read the original abstract

Hot-exciton relaxation in semiconductor nanocrystals (NCs) is often described using perturbative theories, but their accuracy is difficult to assess for realistic exciton--phonon Hamiltonians. Here, we benchmark the perturbative quantum master equation (QME) and several mixed quantum--classical (MQC) methods against fully quantum mechanical dynamics. Using atomistically parameterized models for CdSe core and CdSe/CdS core--shell NCs, we find that bare CdSe exhibits an ultrafast initial decay followed by slower cooling, whereas the core--shell system is dominated by the slower component. Analysis of reduced models shows that the ultrafast component arises from rapid diabatic state mixing driven by thermal fluctuations of low-frequency phonons, rather than from nuclear-assisted energy relaxation. The QME captures the initial fast decay but can fail for the slower relaxation in the diabatic representation, while the mapping approach to surface hopping (MASH) gives the most consistent agreement with both benchmark dynamics and equilibrium populations. These results establish a benchmark for exciton-cooling dynamics in NCs and clarify the physical regimes in which widely used approximate methods are reliable.

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.

Desk Editor's Note private letter to a colleague

Benchmarks show MASH tracks full quantum dynamics best on these NC models, with ultrafast decay traced to diabatic mixing rather than relaxation, but the Hamiltonians have no external validation.

read the letter

The paper's core contribution is a side-by-side comparison of full quantum dynamics against QME and several MQC methods on atomistically parameterized CdSe and CdSe/CdS exciton-phonon models. It reports that bare CdSe shows an initial fast decay from thermal phonon-driven diabatic mixing, while the core-shell case is slower overall, and that MASH reproduces both the time traces and the long-time populations more reliably than the other approximations.

The reduced-model analysis that isolates mixing from energy relaxation is a clear plus; it gives a concrete mechanistic picture instead of just reporting numbers. The full quantum results serve as an internal benchmark, which is the right way to test the approximate methods.

The main limitation is that the Hamiltonians themselves are not checked against experiment or higher-level calculations. The abstract and description give no comparison to measured absorption spectra, relaxation times, or phonon spectra, so the separation between mixing and relaxation, and therefore the ranking of the methods, stays tied to these particular fitted models. If the low-frequency couplings or diabatic gaps are systematically off, the claimed regimes of validity shift.

This is a methods paper aimed at people who already work on exciton dynamics in nanocrystals or who need to choose between QME and surface-hopping variants for similar systems. It is not a broad advance in theory or a new framework, but the concrete benchmark numbers and the mixing insight are useful within the subfield.

It is worth sending to peer review. The comparisons are grounded enough to merit referee attention even if the model validation needs strengthening.

Referee Report

2 major / 1 minor

Summary. The manuscript benchmarks the perturbative quantum master equation (QME) and several mixed quantum-classical (MQC) methods, including the mapping approach to surface hopping (MASH), against fully quantum mechanical dynamics for hot-exciton cooling in CdSe core and CdSe/CdS core-shell nanocrystals. Using atomistically parameterized exciton-phonon Hamiltonians, it reports an ultrafast initial decay in bare CdSe driven by diabatic state mixing from low-frequency phonons (rather than energy relaxation), slower cooling in the core-shell system, limitations of QME in the diabatic representation, and superior consistency of MASH with both the benchmark dynamics and equilibrium populations.

Significance. If the Hamiltonians are faithful, the work supplies a useful benchmark set for approximate methods on realistic NC models and clarifies the physical regimes separating diabatic mixing from phonon-assisted relaxation. The direct comparison to full quantum dynamics on atomistically derived Hamiltonians is a methodological strength.

major comments (2)

[Abstract] Abstract and model-construction paragraphs: the central claims (ultrafast component arises from diabatic mixing, MASH is most consistent) rest on the atomistically parameterized exciton-phonon Hamiltonians faithfully separating diabatic gaps from low-frequency phonon couplings. No external validation (computed vs. measured relaxation timescales, absorption spectra, or phonon densities of states) is supplied, rendering the full-quantum reference and method ranking model-dependent rather than general.
[Reduced models] Reduced-model analysis section: the attribution of the ultrafast decay specifically to thermal fluctuations of low-frequency phonons (rather than nuclear-assisted relaxation) requires quantitative demonstration that the full quantum trajectories on the chosen Hamiltonians reproduce independent observables; absent such checks, the mechanism identification cannot be cleanly separated from possible Hamiltonian inaccuracies.

minor comments (1)

The abstract would be clearer if it stated the quantitative metrics (e.g., population RMSE, equilibrium deviation) used to rank method agreement.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their positive assessment of the significance of our work and for the constructive major comments. We address each point below and propose revisions to clarify the scope and limitations of our study.

read point-by-point responses

Referee: [Abstract] Abstract and model-construction paragraphs: the central claims (ultrafast component arises from diabatic mixing, MASH is most consistent) rest on the atomistically parameterized exciton-phonon Hamiltonians faithfully separating diabatic gaps from low-frequency phonon couplings. No external validation (computed vs. measured relaxation timescales, absorption spectra, or phonon densities of states) is supplied, rendering the full-quantum reference and method ranking model-dependent rather than general.

Authors: We agree that our conclusions are tied to the specific Hamiltonians used and that the absence of direct comparison to experimental data makes the results model-dependent. The manuscript's focus is on providing a benchmark comparison of methods using the same atomistically parameterized models for CdSe and CdSe/CdS NCs, as detailed in the methods. The parameterization follows established procedures from prior literature on these systems. To address this, we will revise the abstract and relevant paragraphs to explicitly note that the findings pertain to these models and to emphasize the model dependence. We cannot provide new external validation within the scope of this study. revision: partial
Referee: [Reduced models] Reduced-model analysis section: the attribution of the ultrafast decay specifically to thermal fluctuations of low-frequency phonons (rather than nuclear-assisted relaxation) requires quantitative demonstration that the full quantum trajectories on the chosen Hamiltonians reproduce independent observables; absent such checks, the mechanism identification cannot be cleanly separated from possible Hamiltonian inaccuracies.

Authors: The reduced models are constructed by selectively including or excluding phonon modes from the full Hamiltonian to isolate the role of low-frequency modes in driving diabatic mixing. Since the full quantum dynamics are solved exactly on the Hamiltonian, the observed ultrafast decay is a direct consequence of the dynamics on that model. We recognize that without matching to independent observables, the physical interpretation remains conditional on the Hamiltonian's accuracy. We will add clarifying text in the reduced-model section to state this explicitly and to note that the mechanism attribution is within the context of the parameterized model. revision: partial

standing simulated objections not resolved

External validation of the atomistically parameterized Hamiltonians against experimental data such as relaxation timescales or spectra.

Circularity Check

0 steps flagged

No circularity detected; derivation self-contained against independent benchmarks

full rationale

The paper computes full quantum dynamics on atomistically parameterized exciton-phonon Hamiltonians for CdSe and CdSe/CdS models, then benchmarks QME and MQC methods (including MASH) against those dynamics and equilibrium populations. No quoted step reduces a claimed prediction to a fitted parameter by construction, invokes a self-citation as the sole justification for a uniqueness theorem, or renames an input as an output. The full quantum reference is generated on the same Hamiltonians but is not equivalent to the approximate methods' results; the analysis of diabatic mixing versus relaxation is performed via reduced models without self-referential definitions. This is the normal case of an independent benchmark comparison.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

Only abstract available; ledger entries inferred from stated modeling choices.

axioms (1)

domain assumption Atomistically parameterized exciton-phonon Hamiltonians for CdSe and CdSe/CdS NCs accurately represent real physical systems.
All dynamics and benchmarking rest on these models being faithful.

pith-pipeline@v0.9.1-grok · 5766 in / 1164 out tokens · 37776 ms · 2026-06-29T14:26:29.384511+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

1 extracted references · 1 canonical work pages

[1]

1X. Zhou, D. Ward, J. Martin, F. Van Swol, J. Cruz-Campa, and D. Zubia, Phys. Rev. B88, 085309 (2013). 2S. Plimpton, J. Comput. Phys.117, 1 (1995). 3E. Rabani, B. Hetenyi, B. J. Berne, and L. E. Brus, J. Chem. Phys.110, 5355 (1999). 4D. Jasrasaria, D. Weinberg, J. P. Philbin, and E. Rabani, J. Chem. Phys.157, 020901 (2022). 5S. Toledo and E. Rabani, J. Co...

work page doi:10.1063/5.0106546 2013

[1] [1]

1X. Zhou, D. Ward, J. Martin, F. Van Swol, J. Cruz-Campa, and D. Zubia, Phys. Rev. B88, 085309 (2013). 2S. Plimpton, J. Comput. Phys.117, 1 (1995). 3E. Rabani, B. Hetenyi, B. J. Berne, and L. E. Brus, J. Chem. Phys.110, 5355 (1999). 4D. Jasrasaria, D. Weinberg, J. P. Philbin, and E. Rabani, J. Chem. Phys.157, 020901 (2022). 5S. Toledo and E. Rabani, J. Co...

work page doi:10.1063/5.0106546 2013