{"id":"75ade887-a72d-4fd5-a8f0-da4726790e41","arxiv_id":"1908.05379","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"An intense femtosecond pulse creates depth-dependent electron-hole quasi-temperatures that decay smoothly in silicon and stepwise in 3C-SiC, driven by tunneling versus multiphoton absorption and by band structure.","lead":"This paper simulates how intense femtosecond laser pulses create and heat electron-hole populations near silicon and silicon carbide surfaces. It shows that the resulting quasi-temperatures decay with depth in material-specific ways, which is relevant for modeling laser machining.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The deep-region quasi-temperatures may be artifacts of mapping highly nonthermal carrier distributions onto a single Fermi-Dirac temperature; Eq. (7) supplies an energy moment, not a validated thermal fit.","rationale":"The reader's weakest assumption identifies the same load-bearing point: a single quasi-temperature per carrier type is assumed to represent laser-generated occupation functions that are not demonstrably thermalized on the 48 fs timescale. I sharpen this to a concrete inversion problem: Eq. (7) computes only a first energy moment, and the mapping from that moment to a Fermi-Dirac temperature is unique only if the distribution is already known to be thermal. Figs. 2 and 7 suggest the opposite in the deep regions, where narrow, position-independent peaks dominate. The paper's own wording ('We can assume') flags this step, and no numerical validation is provided. The TDDFT-Maxwell framework itself is established and the carrier-density results may be sound; however, the central scientific claim of the paper is about quasi-temperatures, so the unvalidated projection is decisive. Because the concern does not change the reader's conditional verdict, I recommend keeping the verdict unchanged.","tokens_in":11038,"tokens_out":5073,"duration_ms":57528,"concrete_test":"Recompute Te and Th at each macroscopic layer X by directly fitting the computed occupations O^{k}_{X,i}(Te) to a Fermi-Dirac distribution f(epsilon; mu, T) using the LDA band structure, and compare the fitted temperatures and residual errors with the QTM/STM values from Eq. (7). If the deep-region distributions (X > 0.1 um for Si, X > 0.3 um for 3C-SiC) cannot be fit by any single (mu, T) or the fitted T differs materially from the reported quasi-temperature, then the finite deep-region quasi-temperatures are projection artifacts and the central claim is not supported.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Section III extracts Te and Th from Eq. (7) by computing the band-resolved internal energies U_{X,c} and U_{X,v} and then 'assuming' the QTM/STM quasi-temperature model of ref. [24]. The simulation runs only to 48 fs, and Figs. 2 and 7 show that deep-region occupations are not broad thermal tails: they are narrow, position-independent peaks, e.g., at approximately ±3 eV in 3C-SiC, attributed to four-photon absorption. A Fermi-Dirac temperature is not uniquely determined by the first energy moment; for an arbitrary occupation it is an effective parameter that encodes peak position and carrier density, not a thermalization temperature. At low carrier density, a narrow peak located several eV above the band edge maps to a large 'quasi-temperature' even if the distribution has no thermal width. Hence the central finding that Te and Th remain finite far from the surface may be a mathematical consequence of the projection in Eqs. (5)-(7), not a physical property of the excited carriers. The paper provides no goodness-of-fit check, no convergence test for the 100 (Si) or 84 (3C-SiC) conduction bands, and no comparison with a direct Fermi-Dirac fit, so the band-structure-dependent stepwise Te/Th curves are not yet tied to an actual thermalized population.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper reports multiscale time-dependent density functional theory (TDDFT) plus Maxwell equation simulations of femtosecond laser excitation at silicon and 3C-SiC surfaces. It computes the depth-dependent absorbed energy, electron-hole density, and electron and hole quasi-temperatures extracted from band-resolved occupation probabilities. The main findings are that, near the surface, the quasi-temperatures decrease exponentially with depth, but they remain finite at larger depths because of multi-photon absorption peaks that produce position-independent narrow structures; in silicon the depth profile is smooth and fitted by a double exponential, whereas in 3C-SiC it shows a stepwise decrease attributed to changes in the contributing bands. The paper also connects the depth variation of the Keldysh parameter with the transition between multiphoton and tunneling excitation regimes.","tokens_in":11333,"tokens_out":3036,"duration_ms":31331,"significance":"If the quasi-temperature extraction is valid, the paper provides a first-principles, parameter-free description of the depth-resolved electron-hole distribution in two important semiconductors under femtosecond irradiation, which is useful input for two-temperature-model simulations of laser processing. The simulation setup itself is creditable: the multiscale TDDFT+Maxwell scheme is self-contained, the Keldysh classification is an external criterion, and no target quantity is fitted to force the reported trends. The qualitative consistency between the Keldysh parameter crossing and the dip in the hole quasi-temperature is a genuinely interesting observation that would be valuable if supported by a validated thermalization analysis. The principal weakness is that the central quantitative output, the quasi-temperatures Te and Th, depends on an unvalidated projection onto a quasi-temperature model.","major_comments":[{"comment":"The quasi-temperature extraction is load-bearing and is not validated. Equation (7) defines only the band-resolved internal energies UX,c and UX,v; the text then states 'We can assume the quasi-temperatures' from these moments using the QTM/STM model of reference [24], but does not give the inversion equations, the assumed occupation form, or any consistency check. Figures 2 and 7 show that in the deep region the occupations are narrow, position-independent peaks (for example, around ±3 eV in 3C-SiC, attributed to four-photon absorption), not broad thermal tails. For such a distribution the first energy moment does not uniquely determine a Fermi-Dirac temperature; a narrow peak located several eV above the band edge can map to a large effective 'quasi-temperature' even if the distribution has no thermal width. Therefore the central claim that Te and Th remain finite far from the surface may be an artifact of the projection procedure rather than a property of a thermalized carrier population. The authors should show a direct Fermi-Dirac fit to the occupation functions, report goodness-of-fit or residual measures, and validate the QTM/STM model on nonthermal, low-density distributions before the quantitative temperatures are used to support the conclusions.","section":"Section III, Eq. (7) and text following"},{"comment":"The laser frequency omega_0 is never specified. The pulse envelope and total duration are given (Tp = 10.81 fs, Te = 48.38 fs) and intensities are listed, but the carrier frequency is absent from Eq. (4) and from the parameter discussion. This omission prevents the reader from checking the assignment of the ±3 eV peaks in 3C-SiC to four-photon absorption (which implies a photon energy of about 1.5 eV), and from recomputing the Keldysh parameter plotted in Figs. 4 and 8. The wavelength or angular frequency must be stated explicitly, together with the effective mass and gap values used in the Keldysh formula.","section":"Section II, Eqs. (2) and (4); Figures 4 and 8"},{"comment":"The conclusions about band-structure-dependent stepwise quasi-temperatures are not backed by convergence checks. The silicon calculation uses 16^3 k-points and approximately 100 conduction bands, while 3C-SiC uses 8^3 k-points and 84 conduction bands; the macroscopic mesh spacing is 13 nm and the simulation ends at 48 fs. Since the stepwise features in Fig. 7(d)-(f) are attributed to 'the change of the number of contributing bands,' the result should be shown to be converged with respect to the number of conduction bands, the k-point sampling, and the macroscopic mesh size. A convergence test for Te and Th would materially strengthen the claim that the band structure, rather than a truncation artifact, causes the stepwise behavior.","section":"Section III.A and III.B, numerical parameters"}],"minor_comments":[{"comment":"There are typos in the abstract: 'silicone' should be 'silicon' and 'descrease' should be 'decrease'.","section":"Abstract"},{"comment":"The silicon setup states 'The cubic unit cell containing eight carbon atoms was discretized into grids of 16^3.' This should presumably read 'eight silicon atoms.'","section":"Section III.A, first paragraph"},{"comment":"The sentence '3C-SiC shows a stepwise structure in Fig. 6 (d)-(f)' should refer to Fig. 7, since the stepwise quasi-temperature curves appear in panels (d)-(f) of Fig. 7.","section":"Section III.B, text near Fig. 7"},{"comment":"The sentence 'the deviation between T1 and T1 indicates that a change in the excitation process occurs' should read 'deviation between T1 and T2.'","section":"Section III.A, paragraph on fitting parameters"},{"comment":"Quasi-temperatures are reported in eV but the text also refers to 'temperature'; the authors should state the conversion to Kelvin or explicitly define the units of Te and Th.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within scope for an applied-physics journal and the simulation scheme is interesting. My main concern is that the quantitative central claim rests on the unvalidated quasi-temperature projection from reference [24]; this is fixable in revision by adding direct fits and convergence tests. The missing laser frequency is a simple but important omission that also must be corrected. I would not reject, but the revision needs to be substantive rather than cosmetic."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"You can skip the formalism and go straight to Figs. 3 and 7: the genuinely new output is depth-resolved electron/hole quasi-temperatures for Si and 3C-SiC under femtosecond excitation, computed with the established TDDFT+Maxwell (SALMON) scheme. The observation that the deep-region excitation is dominated by a specific multiphoton channel that leaves quasi-temperatures roughly intensity-independent, and that 3C-SiC shows a stepwise decay where Si is smooth, is new and useful for people who build two-temperature models. The correlation between the hole-temperature dip and the Keldysh parameter crossing is a nice internal consistency check.\n\nThe soft spot is the temperature extraction, and it is load-bearing. Eq. (7) computes band-resolved internal energies, then the QTM/STM model from ref. 24 is assumed to convert those energies (plus carrier density) into a single quasi-temperature per carrier type. But at 48 fs and at low density, the occupations are not Fermi-Dirac-like: Fig. 2 and Fig. 7 show narrow, position-independent peaks (e.g., ±3 eV in SiC, four-photon channel). For such a distribution, the energy moment does not determine a temperature; the QTM output is an effective parameter encoding peak position and density, not a thermalization temperature. So the headline claim that Te and Th stay finite deep in the material may be an artifact of the projection, not a physical property of the carriers.\n\nThe paper also lacks convergence checks on the 100 (Si) or 84 (SiC) conduction bands, gives no direct Fermi-Dirac fit comparison or goodness-of-fit measure, and no code/data deposit. The double-exponential fit (Eq. 8) reproduces the computed Te(X) but is descriptive; interpreting T0 and τ2 as separate physical mechanisms goes beyond the evidence. The LDA band-gap underestimate is acknowledged but not corrected, which matters for quantitative values.\n\nNone of this makes the paper worthless. The calculation is self-contained, the method is established, and the qualitative trends (band-structure dependence, multiphoton vs tunneling regime) are likely robust. But the quantitative quasi-temperatures should be treated as model-dependent effective parameters until the extraction is validated against direct thermal fits or, ideally, experimental data.\n\nI'd send it to peer review—a serious referee could push for the validation and convergence tests. For the reader, this is a conditional accept at best; for the femtosecond-processing community it's a useful pointer, not a reliable input table. I would not cite the temperature values, though I might cite the methodology example.","headline":"New depth-resolved quasi-temperature profiles from TDDFT+Maxwell for Si and 3C-SiC, but the deep-region temperatures are likely artifacts of projecting nonthermal occupations onto a single Fermi-Dirac form.","tokens_in":11863,"tokens_out":2362,"would_cite":false,"duration_ms":24241,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Electron and hole quasi-temperatures stay finite tens of nanometers deep in silicon and 3C-SiC, with profiles set by tunneling versus multiphoton absorption and by band structure.","keywords":["electron-hole quasi-temperature","femtosecond laser processing","silicon","3C-SiC","time-dependent density functional theory","Maxwell's equation","multiphoton absorption","Keldysh parameter"],"falsifier":"A time-resolved photoemission or transient-absorption experiment that measures the carrier energy distribution at depths beyond 100 nm in silicon under 800 nm, 10 fs-scale pulses at ~$10^{13}$ W/$cm^{2}$ would settle the claim: if the deep distribution shows only narrow, intensity-independent multiphoton peaks without a broad ~1 eV thermal tail, the extracted quasi-temperatures are artifacts of the temperature-model projection.","tokens_in":10834,"feed_emoji":"⚡","tokens_out":5515,"duration_ms":49926,"temperature":0.7,"pith_summary":"This paper uses a first-principles simulation that couples time-dependent density functional theory with Maxwell's equations to track how a femtosecond laser pulse excites electrons and holes in silicon and cubic silicon carbide (3C-SiC). It claims that the electron and hole quasi-temperatures, extracted at each depth from the simulated occupation functions, do not simply vanish with distance: they remain finite (around 1 eV in silicon) even hundreds of nanometers inside the material because a multiphoton absorption channel keeps populating well-defined energy states. It further claims that the spatial profile of these temperatures reveals the transition from tunneling excitation at the surface to multiphoton absorption deeper inside, and that the profile depends strongly on the material's band structure. Such depth-resolved carrier distributions are what two-temperature models of laser processing need as input.","feed_headline":"Hot carriers persist deep under femtosecond laser","feed_subtitle":"Simulation links carrier temperature profile to tunneling vs multiphoton absorption and band structure.","key_machinery":"The central machinery is a multiscale implementation that solves the Maxwell wave equation for the laser pulse on a macroscopic depth coordinate X, while at each X a time-dependent Kohn–Sham equation evolves the lattice-periodic electron orbitals in the microscopic coordinate r. Occupations of conduction and valence states at the final time are obtained by projecting the time-dependent orbitals onto the initial ground-state orbitals, giving the change in electron distribution δO at every depth. From those occupations, the reduced internal energies in the valence and conduction bands (Eq. 7) are used with the quasi-temperature model of Ref. 24 to assign electron and hole quasi-temperatures per depth, and the Keldysh parameter γ classifies the local excitation process as tunneling (γ ≪ 1) or multiphoton (γ ≫ 1).","core_discovery":"Under an 800-nm femtosecond pulse at intensities around 5×$10^{12}$ to 2×$10^{13}$ W/$cm^{2}$, the excited electron and hole occupation functions in silicon and 3C-SiC acquire quasi-temperatures whose spatial decay is not monotonic. In silicon, the electron quasi-temperature drops steeply in the first 100 nm, passes through a small dip in the hole temperature where the Keldysh parameter crosses unity, and then saturates at ~1.0 eV independent of laser intensity; the full profile fits a double exponential with a constant offset. In 3C-SiC, the quasi-temperatures instead show stepwise drops as the number of contributing conduction bands changes with depth, and the hole temperature ends up much higher than the electron temperature. The paper concludes that both the excitation mechanism (tunneling versus multiphoton) and the band structure determine the macroscopic electron-hole distribution, and that quantum-mechanical simulation is needed to supply electron-lattice models with accurate initial carrier temperatures.","pith_inferences":["If the deep-region carriers are as hot as the simulation suggests, then the initial electron-lattice coupling in the two-temperature model begins from a hot and spatially extended carrier gas, which would shift predicted damage thresholds and heat-affected zones deeper than in current models.","The stepwise drop in 3C-SiC suggests a testable prediction: changing the laser wavelength to move a different set of conduction bands into resonance should change the depth at which the steps occur.","The simulation stops at 48 fs, before carrier–phonon scattering; one could extend the same method with a longer time window or couple it to molecular dynamics to see how quickly the hot deep-region distribution thermalizes with the lattice."],"forward_implications":["Two-temperature-model simulations of femtosecond laser processing should not assume that carrier temperatures decay monotonically; a finite deep-region temperature must be included.","The value of the saturation temperature in silicon (~1 eV) is set by a multiphoton absorption channel and is nearly independent of laser intensity.","The stepwise temperature profiles in 3C-SiC show that band-structure details, not just absorbed energy, control the carrier distribution; simulations for each material are needed.","The Keldysh parameter alone can locate the crossover depth between tunneling and multiphoton regimes, as the crossover shows up directly as a feature in the hole quasi-temperature."],"supporting_citations":[{"why":"Supplies the multiscale TDDFT+Maxwell implementation that simultaneously propagates the laser field and the electron orbitals.","marker":"[19]"},{"why":"Provides the quasi-temperature model (QTM/STM) used to assign electron and hole temperatures from the simulated occupation functions.","marker":"[24]"},{"why":"Defines the Keldysh parameter used throughout to classify local excitation as tunneling or multiphoton.","marker":"[25]"},{"why":"Gives the projection of time-dependent orbitals onto initial states that yields the electron-hole occupation changes in Eq. (5).","marker":"[23]"},{"why":"Specifies the LDA exchange-correlation functional used in the time-dependent Kohn-Sham equation.","marker":"[22]"},{"why":"Establishes the first-principles method for electron dynamics in periodic solids that the multiscale calculation builds on.","marker":"[20]"},{"why":"Extends the formalism to laser-driven electron dynamics and provides validation for treating the coupled field-electron system.","marker":"[21]"}],"fun_headline_variants":["Carrier heat persists deep inside silicon and 3C-SiC","Stepwise carrier cooling in 3C-SiC under femtosecond laser","Laser-heated carriers cool nonuniformly with depth in Si and SiC","Quasi-temperature profiles reveal band-structure effects in Si and SiC"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result rests on assuming that at each depth the entire nonthermal, laser-generated carrier population can be represented by one quasi-temperature per carrier type; if the carriers are not actually thermalized within the 48 fs simulation time, the reported temperatures are projections of the model rather than physical temperatures.","fun_headline_variants_meta":{"raw":{"variants":["Carrier heat persists deep inside silicon and 3C-SiC","Stepwise carrier cooling in 3C-SiC under femtosecond laser","Laser-heated carriers cool nonuniformly with depth in Si and SiC","Quasi-temperature profiles reveal band-structure effects in Si and SiC"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000297,"raw_usage":{"total_tokens":1680,"prompt_tokens":860,"completion_tokens":820,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":476,"completion_tokens_details":{"reasoning_tokens":740}},"tokens_in":476,"tokens_out":820,"duration_ms":8979,"temperature":1.0,"reasoning_tokens":740,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:16:11.799046+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A time-resolved photoemission or transient-absorption experiment that measures the carrier energy distribution at depths beyond 100 nm in silicon under 800 nm, 10 fs-scale pulses at ~$10^{13}$ W/$cm^{2}$ would settle the claim: if the deep distribution shows only narrow, intensity-independent multiphoton peaks without a broad ~1 eV thermal tail, the extracted quasi-temperatures are artifacts of the temperature-model projection.","supporting_citations":[{"cited_title":"Otobe , author T","cited_arxiv_id":null,"evidence_quote":"Provides the quasi-temperature model (QTM/STM) used to assign electron and hole temperatures from the simulated occupation functions."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the Keldysh parameter used throughout to classify local excitation as tunneling or multiphoton."},{"cited_title":"Yabana , author T","cited_arxiv_id":null,"evidence_quote":"Extends the formalism to laser-driven electron dynamics and provides validation for treating the coupled field-electron system."}],"review_version":1}