{"id":"effc69a6-d6cd-40d3-90e8-deb560ced0fc","arxiv_id":"1908.07648","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"The nonlinear expansion of the photoluminescence spot in h-BN encapsulated WSe2 monolayers is dominated by kinetic-energy relaxation of hot excitons, saturating at high density due to Auger heating balancing phonon cooling.","lead":"This paper shows that the fast, time-dependent spreading of light emission in atomically thin WSe2 comes from hot excitons cooling down, not from an intrinsic anomalous diffusion. The result explains a long-standing puzzle in 2D semiconductor physics and could guide design of exciton-based devices.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The fit underpinning the central claim assumes tau* equals the 0.23 ns PL lifetime with no independent cooling-time measurement; if tau* is shorter, the extracted excess temperatures and saturation in Fig. 4 may be fitting artifacts.","rationale":"The reader correctly identified the unvalidated tau* = lifetime assumption as a load-bearing weakness. I agree that this is the primary concern: the model in Eq. (3) cannot yield reliable T* values unless tau* is known independently, and the saturation that is the paper's central observable is extracted from that same fit. I also note a second, more basic issue that the reader did not emphasize: the manuscript states MSD = 2 D(t) t for a time-dependent diffusivity, whereas the diffusion equation requires MSD = 2∫ D(s) ds. If the former expression was used, the fitted T* values are mathematically biased, independent of the tau* question. This does not automatically falsify the qualitative hot-transport picture, because the photon-energy-dependent saturation is an independent cross-check; however, that cross-check still passes through the same model. The appropriate verdict remains CONDITIONAL: the central claim is plausible and well-motivated, but it should not be accepted until the model is tested with a free tau* (or a direct measurement of the cooling time) and the corrected integral relation. The proposed re-fit is a single, concrete check that would settle whether the saturation is physical or a fitting artifact. I therefore keep the reader's verdict unchanged, while noting that the integration issue strengthens the case for conditional acceptance rather than full acceptance.","tokens_in":7727,"tokens_out":6532,"duration_ms":158712,"concrete_test":"Re-fit the MSD transients behind Fig. 1 using the correct integrated relation MSD(t) = 2∫_0^t D(s) ds and treat tau* as a free parameter (or determine it independently, e.g., from time-resolved PL spectra probing the exciton temperature); recompute the Fig. 4 excess-temperature saturation for both 405 nm and 520 nm excitation. If the saturation does not persist for both photon energies with a physically reasonable tau*, the central claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim rests on extracting the excess exciton temperature T* by fitting the MSD transient with the model in Eq. (3), where the kinetic-energy relaxation time tau* is set equal to the measured PL lifetime (0.23 ns) on the assumption that the gas completely relaxes before recombination. No independent measurement of the hot-exciton cooling time is provided, yet tau* controls the shape of the entire transient: if tau* is much shorter than 0.23 ns, as is typical for exciton-phonon cooling in TMDs at room temperature, then T* must be much larger to reproduce the early fast rise, and the density dependence of T*—including the claimed saturation in Fig. 4—is not identifiable from the MSD alone. A separate, compounding issue is that the text writes MSD = 2 D(t) t for a time-dependent diffusivity; the correct relation is MSD = 2∫_0^t D(s) ds. If the former was used in the fitting, the extracted T* values are biased in a way that could create or mask the saturation. Because both issues enter the same fit that produces Fig. 4, the central claim is not independently supported without a direct test of tau* or a corrected, free-tau fit.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript reports spatiotemporally resolved photoluminescence measurements of exciton transport in h-BN encapsulated WSe2 monolayers at room temperature. The authors observe that the mean squared displacement (MSD) of the exciton cloud evolves nonlinearly at early times, with the initial expansion rate increasing with excitation density and saturating at high densities. They attribute the fast early expansion to hot exciton transport, model the exciton temperature as T(t)=T_L+T*exp(-t/tau*), and extract an excess temperature T* by fitting the MSD transients. They also compare two excitation photon energies and find that the saturation density is higher for lower photon energy, which they interpret as evidence for a balance between Auger-assisted hot-exciton generation and phonon-assisted cooling. They argue that Auger broadening of the PL profile is negligible at the relevant densities based on independent time-resolved PL measurements.","tokens_in":8012,"tokens_out":4129,"duration_ms":480568,"significance":"If the central claim is correct, the work provides a useful experimental demonstration of hot-exciton transport at room temperature in a technologically relevant TMD monolayer and offers a mechanism for the previously reported time-dependent exciton diffusivity. The strengths of the paper are the direct spatiotemporal measurement, the careful control for Auger broadening via independent TRPL and integrated-PL analysis, and the photon-energy comparison that provides a physically motivated cross-check. However, the main model relies on an unmeasured assumption that the kinetic-energy relaxation time equals the exciton lifetime, and the possible use of an incorrect MSD-to-diffusivity relation could bias the extracted temperatures. These issues affect the central quantitative claim and require additional analysis or experiments.","major_comments":[{"comment":"The kinetic-energy relaxation time tau* is set equal to the measured PL lifetime (0.23 ns) solely on the assumption that the exciton gas completely relaxes before recombining; no independent measurement of the exciton cooling time is provided. Because tau* controls the time-decay of the instantaneous diffusivity in the model, the fitted T* values and their density dependence, including the saturation shown in Fig. 4, are identifiable only if tau* is known. If the actual cooling time is substantially shorter than 0.23 ns, as is common for exciton-phonon relaxation in TMDs, the fitted T* values would be correspondingly different and the reported saturation could be an artifact of the assumed tau*. The manuscript should either provide an independent determination of tau* (for example, from time-resolved PL linewidth or transient absorption) or, failing that, present fits with tau* treated as a free parameter, with confidence intervals, and show that the saturation and the photon-energy trend survive.","section":"Eq. (3), Fig. 4"},{"comment":"The text writes the MSD as <Δr^2(t)> = 2D(t)t for a time-dependent diffusivity, but the correct relation is <Δr^2(t)> = 2∫_0^t D(s) ds. If the non-integrated expression was used in the fits, the extracted excess temperatures are biased in a way that could either create or mask the saturation. The authors should state explicitly which expression was actually used in the fitting, and if the instantaneous form was used, the analysis should be repeated with the integrated form.","section":"MSD definition before Eq. (3)"},{"comment":"The saturation of T* with increasing excitation density is a fitted trend rather than an independent prediction, since T* is the amplitude of the exponential term in Eq. (3) that is adjusted to match the early-time MSD rise. The photon-energy comparison provides a useful consistency check, but it involves the same model and the same fitted T*, so it does not by itself resolve the degeneracy between tau* and T*. A quantitative test of the proposed Auger-heating/phonon-cooling balance, or at least a sensitivity analysis showing that the density-dependent saturation is robust to the model assumptions, is needed before the abstract's mechanistic claim can be considered established.","section":"Fig. 4 and abstract"}],"minor_comments":[{"comment":"Many inline equations use placeholder-like symbols (e.g., \"〈Δσ%(')〉\", \"8(')\", \"DEF\") that appear to be OCR artifacts; the manuscript should be typeset with clean mathematical notation.","section":"Throughout"},{"comment":"References 14 and 37 are the same Bouchaud and Georges paper, and reference 55 (Najafi et al., Nat. Commun. 8, 15177) is duplicated by reference 59; the reference list should be deduplicated.","section":"References"},{"comment":"The Auger constant is estimated using an assumed 11.5% absorption at 405 nm; the authors should provide the resulting uncertainty in C_A and in the density threshold below which Auger broadening is claimed to be negligible.","section":"Fig. 2 and supplementary material"},{"comment":"Figure 1 and Figure 4 were obtained on different samples; the manuscript should state this explicitly and discuss how sample-to-sample variation is controlled when comparing the two data sets.","section":"Figures 1 and 4"},{"comment":"The text refers to supplementary Figures S3 and S4 and to supplementary fabrication details, but no supplementary material is included with the manuscript; these should be provided or the references to them should be amended.","section":"Supplementary material"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Quick read for you. The paper claims that the nonlinear MSD in WSe2 monolayers comes from the relaxation of hot exciton kinetic energy, not from Auger broadening or disorder. That claim is new for TMDs, and I think it is probably right in its qualitative form. The best part is the control experiment: they vary the excitation photon energy (405 vs 520 nm) and find the saturation of the initial expansion shifts to higher density for lower photon energy. That is a directional prediction of the hot exciton picture, and it goes the right way.\n\nThe Auger-broadening accounting is also solid. They fit TRPL with the rate equation, get a consistent Auger constant from two methods, and show it cannot explain the saturation. For the low-density regime they quantify the broadening correction. Credit where due.\n\nThe soft spots are in the model, and they are load-bearing for the quantitative conclusions. They assume the kinetic energy relaxation time tau* equals the 0.23 ns PL lifetime, with no independent measurement of the cooling time. If tau* is actually shorter (which is typical for exciton-phonon cooling at room temperature), the extracted excess temperatures are inflated and the density-dependent saturation in Fig. 4 could be a fitting artifact. Second, the model writes MSD = 2 D(t) t with time-dependent D(t). The correct relation is MSD = 2∫_0^t D(s) ds. Using the product form changes the shape of the fit and biases T*. Both issues enter the same fit that produces Fig. 4, so I would not trust the absolute T* values or the exact saturation density yet.\n\nThat said, the central conclusion has independent support from the photon-energy comparison, which is not a fitted artifact in the same way. The authors need to re-fit with free tau* (or measure it), use the integral formula, and show the raw MSD data with error bars and deposited code/data. Then this would be a solid paper.\n\nThis is for the community working on exciton transport in 2D materials. I would keep it in the pipeline. Send it to peer review, but with referees who will insist on the tau* measurement and the corrected MSD relation. I wouldn't cite the numbers as they stand, but the qualitative demonstration is worth knowing.","headline":"Hot exciton transport as the explanation for nonlinear MSD in WSe2 is plausible and has a nice photon-energy cross-check, but the paper's quantitative fit leans on an unverified tau*-equals-lifetime assumption and a wrong MSD=D(t)t integral relation.","tokens_in":8553,"tokens_out":4754,"would_cite":false,"duration_ms":144853,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"In h-BN-encapsulated WSe2 monolayers, the nonlinear spread of photoexcited excitons is caused by the cooling of a hot exciton gas, and the saturation of that spread at high density comes from Auger heating balancing phonon cooling.","keywords":["hot exciton transport","WSe2 monolayer","mean squared displacement","anomalous diffusion","Auger heating","exciton-phonon relaxation","time-resolved photoluminescence","h-BN encapsulation"],"falsifier":"Measure the exciton cooling time independently—for example by time-resolved photoluminescence of the exciton emission energy or by transient absorption after non-resonant pumping at the same densities—and compare it with 0.23 ns. If the observed spectral relaxation or hot-carrier cooling time is much shorter than the lifetime, the MSD fits assign too much weight to a slow exponential temperature decay and the reported excess temperatures and saturation would not survive.","tokens_in":7521,"feed_emoji":"🔬","tokens_out":9409,"duration_ms":79131,"temperature":0.7,"pith_summary":"This paper reports room-temperature, time-resolved photoluminescence imaging of exciton spreading in h-BN encapsulated WSe2 monolayers and argues that the apparent anomalous diffusion—a fast early expansion that slows to a constant rate—is the spatial signature of a hot exciton gas cooling toward the lattice temperature. The initial expansion rate grows with excitation density and saturates at high density, and the density at which saturation appears depends on the excitation photon energy. The authors conclude that the early fast motion is hot exciton transport, with Auger-assisted exciton generation heating the gas and phonon-assisted scattering cooling it until a dynamical balance caps the expansion rate. This matters because time-dependent exciton diffusivities have been seen in TMD monolayers before without a physical explanation, and distinguishing hot from cold exciton transport is directly relevant for room-temperature excitonic devices.","feed_headline":"Hot excitons, not disorder, set WSe2 expansion rate","feed_subtitle":"The gas cools from fast hot-exciton motion to slow diffusion; Auger heating caps the initial speed.","key_machinery":"The central object is the time-dependent exciton diffusivity $D(t)=\\mu_h k_B T(t)/q$, tied to the instantaneous effective temperature $T(t)$ of the exciton gas through the Einstein relation, with $T(t)=T_L+T^*\\exp(-t/\\tau^*)$. This turns a spatial measurement—the mean squared displacement of the Gaussian photoluminescence profile—into a thermometer for the kinetic energy of the hot exciton gas. The supporting machinery is the Auger rate equation $n'(t)=-n/\\tau-C_A n^2$, used to measure the Auger constant and to show that Auger broadening contributes negligibly below roughly $4\\times10^{11}\\,\\mathrm{cm^{-2}}$, plus the photon-energy-dependent saturation measurement that separates hot-exciton transport from Auger artifacts.","core_discovery":"The central claim is that the nonlinear evolution of the mean squared displacement of a non-resonantly excited exciton gas in h-BN encapsulated WSe2 monolayers is dominated by the relaxation of the gas's excess kinetic energy, not by Auger broadening or by disorder-induced anomalous diffusion. Using a time-varying diffusivity $D(t)=\\mu_h k_B T(t)/q$ with an exponentially relaxing temperature $T(t)=T_L+T^*\\exp(-t/\\tau^*)$, the authors fit the measured MSD and extract initial excess temperatures that increase with excitation density and saturate at high densities. The saturation is interpreted as a balance between Auger-assisted hot exciton generation, which heats the gas, and phonon-assisted relaxation, which cools it. A control experiment lowers the excitation photon energy from 3.1 eV to 2.4 eV and finds that saturation shifts to higher excitation densities, as expected if the initial kinetic energy of the hot gas controls the effect.","pith_inferences":["The same cooling-gas picture should apply to other monolayer TMDs and van der Waals heterostructures, predicting stronger apparent time-dependent diffusivity for high photon energy, high density, and low lattice temperature; the paper's control experiment already shows the photon-energy trend.","If the model is correct, the initial expansion speed should scale roughly as $\\sqrt{T^*}$ and therefore depend on the exciton effective mass, which could be tested by comparing different monolayer materials.","The density at which saturation occurs could serve as a quantitative probe of the Auger heating rate: fitting the density dependence of $T^*$ should constrain the ratio of Auger-assisted hot-exciton generation to phonon cooling.","A direct time-resolved measurement of the exciton gas temperature, for example through emission linewidth or phonon-sideband spectroscopy, would test the exponential relaxation assumed in the model without relying on the lifetime approximation."],"forward_implications":["Time-dependent diffusivities measured in TMD monolayers can reflect the cooling of a hot exciton gas rather than anomalous diffusion from disorder.","Initial diffusivities can reach roughly $4\\,\\mathrm{cm^2\\,s^{-1}}$ at high excitation density, far above the cold-exciton value, so density-dependent transport measurements must account for excess kinetic energy.","At high densities, Auger-assisted hot exciton generation sets an upper bound on the gas temperature and hence on the early expansion speed; raising the pump photon energy lowers the density needed to reach that bound.","In h-BN encapsulated monolayers, Auger broadening is small below about $4\\times10^{11}\\,\\mathrm{cm^{-2}}$, so the fast early expansion can be read as genuine transport rather than profile flattening.","Room-temperature excitonic devices can in principle use a fast hot-exciton transport regime before cooling and a slower cold-exciton regime afterward."],"supporting_citations":[{"why":"Supplies the spatiotemporally resolved photoluminescence mapping method used to measure the MSD evolution in WSe2 monolayers.","marker":"7"},{"why":"Establishes the time-correlated single-photon counting scanning technique for building temporally and spatially resolved photoluminescence maps.","marker":"27"},{"why":"Provides the exciton transport visualization approach that the measurement technique extends to TMD monolayers.","marker":"28"},{"why":"Defines Auger broadening and the integrated-PL relation used to cross-check the Auger constant and rule out profile flattening.","marker":"3"},{"why":"Supplies the rate equation $n'=-n/\\tau-C_A n^2$ and the Auger recombination framework used to fit the time-resolved photoluminescence.","marker":"32"},{"why":"Provides the hot-exciton framework in which non-resonant excitation creates excitons with excess kinetic energy that relax via phonon scattering.","marker":"29"},{"why":"Shows exciton-exciton annihilation is strongly suppressed in h-BN encapsulated TMD monolayers, supporting the small Auger-broadening contribution.","marker":"45"},{"why":"Introduces Auger heating as the process by which non-radiative exciton recombination transfers energy to nearby excitons and raises the gas temperature.","marker":"46"},{"why":"Supplies the relation between instantaneous diffusivity and effective carrier temperature used as the basis for the time-varying diffusivity model.","marker":"59"},{"why":"Documents that phonon scattering rates increase with carrier or exciton energy, the premise behind the saturation balance between Auger heating and phonon cooling.","marker":"60"}],"fun_headline_variants":["Hot exciton cooling, not disorder, sets WSe2 expansion","Exciton gas expansion traced to kinetic energy relaxation","WSe2 hot excitons: expansion speed set by cooling, not disorder","Auger heating caps hot exciton expansion in WSe2","Exciton gas expansion slows as excess kinetic energy relaxes"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the kinetic-energy relaxation time $\\tau^*$ of the hot exciton gas can be set equal to the measured exciton lifetime (0.23 ns) in the MSD fit; no independent measurement of the cooling time is provided, and if true cooling is much faster, the extracted excess temperatures and the reported saturation could be fitting artifacts.","fun_headline_variants_meta":{"raw":{"variants":["Hot exciton cooling, not disorder, sets WSe2 expansion","Exciton gas expansion traced to kinetic energy relaxation","WSe2 hot excitons: expansion speed set by cooling, not disorder","Auger heating caps hot exciton expansion in WSe2","Exciton gas expansion slows as excess kinetic energy relaxes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.00073,"raw_usage":{"total_tokens":3212,"prompt_tokens":830,"completion_tokens":2382,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":446,"completion_tokens_details":{"reasoning_tokens":2296}},"tokens_in":446,"tokens_out":2382,"duration_ms":15466,"temperature":1.0,"reasoning_tokens":2296,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T12:00:32.374701+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the exciton cooling time independently—for example by time-resolved photoluminescence of the exciton emission energy or by transient absorption after non-resonant pumping at the same densities—and compare it with 0.23 ns. If the observed spectral relaxation or hot-carrier cooling time is much shorter than the lifetime, the MSD fits assign too much weight to a slow exponential temperature decay and the reported excess temperatures and saturation would not survive.","supporting_citations":[],"review_version":1}