{"id":"5ade2419-8f14-4ab3-ada5-a36c949cca9d","arxiv_id":"2607.18085","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"First-principles KKR-CPA calculations give the first full Te,Tl-dependent electronic transport parameter set for Ti-6Al-4V; lower alloy thermal conductivity is nearly offset by weaker electron-phonon coupling in two-temperature-model lattice heating.","lead":"Researchers computed the electronic heat transport and electron–lattice coupling of titanium and the alloy Ti-6Al-4V from first principles, for use in laser-ablation models. The alloy's heat conduction is much lower, but it barely changes the simulated lattice heating because two opposing effects cancel.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 19% functional-form claim is anchored at lattice temperatures ≈8–9.6 kK, while the κe map stops at 2.5 kK; unspecified high-Tl extrapolation could change the conclusion.","rationale":"The reader's weakest_assumption already targets the high-Tl extrapolation, and my assessment converges on the same point. The central claim rests on the 19% shift in peak lattice temperature, which occurs at Tl values roughly four times the upper bound of the calculated κ_e map. Because the extrapolation rule is not stated, the magnitude of the shift is not independently reproducible, and it is exactly the kind of missing support the reviewing rule asks me to flag. I considered whether the neglect of electron–electron scattering in the Kubo–Greenwood transport calculation is a more fundamental issue; it is plausible and worth checking, but the text is less explicit about it, and the extrapolation gap alone is sufficient to keep the verdict conditional. I do not find a mathematical or logical error in the cancellation argument for Ti vs Ti-6Al-4V. The proposed concrete test—running the same TTM with several explicit high-Tl continuations—would settle whether the headline number is robust. Since the reader already issued CONDITIONAL for this reason, my recommendation is UNCHANGED, not a new verdict.","tokens_in":20603,"tokens_out":5905,"duration_ms":309668,"concrete_test":"Obtain the κ_e(Te,Tl) maps (authors state data are available on request) and rerun the τ_p = 300 fs phase-explosion TTM with three explicit Tl treatments above 2500 K: (i) clamp κ_e to its 2500 K values; (ii) linear decrease to 50% at 2Tm; (iii) use a liquid-Ti Wiedemann–Franz/Drude estimate above 2500 K. Compare T_l,max for Ti, Ti-6Al-4V, and the low-T model. If the low-T overshoot remains ≈19±4% across all three extrapolations, the headline survives; if it varies by more than that, the 19% claim is unverified.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Fig. 4(b) reports T_l,max ≈8.1 kK (Ti) and ≈9.6 kK (low-T model) in the phase-explosion regime, and the central 19% conclusion is drawn there. But the first-principles κ_e(Te,Tl) map is generated from σ(E) at discrete Tl up to 2100 K, linearly extrapolated only to 2500 K (Sec. II B), and Sec. II C explicitly disclaims validity above ≈1.3Tm. The manuscript never states how κ_e is continued to Tl ≈8–9.6 kK. If the map is clamped at its 2500 K values, decays, or is replaced by a liquid-metal estimate above Tm, the computed peak lattice temperature—and hence the 19% overshoot—can shift materially. Since the comparison model κ_e = κ_e0 Te/Tl has its own strong Tl dependence, the discrepancy at high Tl is not fixed without specifying both models consistently. The conclusion that the functional form matters more than the elemental/alloy distinction depends on a number that may be an artifact of an unstated extrapolation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes the electronic transport properties of hcp Ti and the Ti-6Al-4V alloy using the KKR-CPA method with the alloy-analogy model and the Kubo-Greenwood formalism, obtaining temperature-dependent electronic thermal conductivity, electron-phonon coupling, and electronic heat capacity. These parameters are then used in one-dimensional two-temperature-model (TTM) simulations of ultrashort-pulse laser excitation. The principal claims are: (i) the Ti resistivity agrees with ab initio electron-phonon calculations and experiment, including the high-temperature saturation; (ii) Ti-6Al-4V and elemental Ti produce almost identical peak lattice temperatures in TTM despite large differences in transport parameters; and (iii) replacing the first-principles thermal conductivity with the commonly used low-temperature Drude form κ_e = κ_e0 T_e/T_l shifts the peak lattice temperature by up to 19%, so the functional form of κ_e matters more than the elemental/alloy distinction.","tokens_in":20910,"tokens_out":6800,"duration_ms":78125,"significance":"If the central 19% claim is correct, the paper makes a practically important point for modeling ultrafast laser ablation of Ti-6Al-4V: the choice of the κ_e(T_e,T_l) functional form is more consequential than whether one uses elemental Ti or alloy-specific parameters. The paper has clear strengths: the Ti resistivity benchmark against independent ABINIT electron-phonon calculations and experiment (Fig. 1a) is convincing; the comparison between the full generalized transport-coefficient result and the Drude-type approximations (Fig. 3b) is transparent; and the pulse-duration-dependent sensitivity analysis (Fig. 4c,d) provides useful physical insight. However, the quantitative 19% headline is currently not fully supported because it is obtained in a regime of lattice temperature far above the range for which the first-principles κ_e map is defined or documented. The alloy κ_e is also not independently validated, which weakens the 1.4% cancellation claim.","major_comments":[{"comment":"The central 19% claim is anchored in the phase-explosion regime, where the TTM reaches T_l,max ≈ 8.1 kK for Ti and ≈ 9.6 kK for the low-T model. The first-principles κ_e map, however, is computed only up to T_l = 2500 K (Sec. II B: “computed for a discrete set of lattice temperatures up to T_l = 2100 K and then linearly extrapolated to T_l = 2500 K”), and Sec. II C explicitly disclaims strict validity above ≈ 1.3 T_m. The manuscript never states how κ_e(T_e,T_l) is continued from 2500 K to 8–9.6 kK during the TTM integration. Since the comparison model κ_e = κ_e0 T_e/T_l has a strong 1/T_l dependence, the discrepancy between the two models at high T_l is not fixed without specifying the continuation of both consistently. The 19% overshoot could change materially if the first-principles map is clamped at 2500 K, extrapolated with a different functional form, or replaced by a liquid-metal","section":"Sec. II B, II C and Fig. 4(b)"},{"comment":"The conclusion that Ti and Ti-6Al-4V give nearly identical lattice temperatures (1.4% difference) relies on the alloy-specific κ_e, whose peak value is a factor of 6.4 lower than in Ti. Unlike the elemental Ti resistivity, the alloy κ_e is not benchmarked against any independent calculation or measurement. The equilibrium resistivity comparison in Fig. 2(d) is a useful check of the CPA disorder model, but κ_e is an energy-integrated Onsager quantity and its alloy suppression depends on the detailed σ(E) shape, not just the residual resistivity. I recommend at least one independent verification of the alloy transport map (e.g., a Wiedemann–Franz check at low T_e, or a second method) or an explicit estimate of the uncertainty in the alloy κ_e before the 1.4% cancellation is used quantitatively.","section":"Sec. III C and Fig. 3(d)"},{"comment":"The phase-explosion threshold is defined as T_l,max = 0.9 T_crit with T_crit = 8980 K for titanium, and the simulated peak lattice temperatures in Fig. 4(b) are 8.1–9.6 kK. For the low-T model, T_l,max = 9638 K exceeds T_crit, while no phase change or liquid-state transport is included in the TTM. The authors acknowledge that mechanical response and phase change are not considered, and use T_l,max only as a proxy, but for the specific claim about the phase-explosion regime it would be important to clarify what physical meaning is assigned to T_l,max above T_crit and whether the liquid-state κ_e and C_l should not be used there. The first-principles κ_e from the crystalline hcp calculation is formally not applicable in this regime, and the assertion that the electronic DOS changes weakly upon melting [57] does not directly ensure that σ(E) or κ_e behaves the same way.","section":"Sec. II C and TTM setup"}],"minor_comments":[{"comment":"The abstract says “low-temperature Drude limit”, while the TTM comparison in Sec. III D replaces the first-principles conductivity with the linear form κ_e = κ_e0 T_e/T_l, which is the low-temperature limit of Eq. (13), not the full Drude scattering model including electron-electron scattering. Please use consistent terminology.","section":"Abstract and Sec. III D"},{"comment":"State explicitly that “f-character” means angular momentum l = 3, and consider giving the angular-momentum cutoff convergence test in the main text rather than only in the Supplemental Material.","section":"Sec. II A"},{"comment":"The statement “data available from the corresponding authors upon request” is suboptimal for a first-principles dataset that is meant to be reusable. A repository link (e.g., Zenodo or a materials database) would improve reproducibility.","section":"Data availability"},{"comment":"The caption says the κ_e maps are “interpolated ... from GTC calculations at four lattice temperatures”, but the interpolation method is not specified. State whether linear or spline interpolation is used and how the maps behave outside the 300–2500 K T_l range.","section":"Fig. 3(c,d)"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal and has a valuable, usable dataset for the laser-ablation community. The Ti benchmark is solid and the sensitivity analysis is informative. The main reason for major revision, rather than rejection, is that the headline quantitative claim (19%) depends on an unspecified and formally invalid high-T_l extrapolation of κ_e; this is fixable by documenting the continuation and testing its sensitivity. The alloy κ_e validation should also be strengthened, but this is secondary. I have no concerns about citation practice or novelty disclosure."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the genuinely new thing here is the alloy-specific electronic transport parameter set for Ti-6Al-4V under electron-phonon nonequilibrium — κe, Ce, G — computed with KKR-CPA-AAM and Kubo–Greenwood, and the observation that G and κe changes cancel in the TTM, leaving Ti and Ti-6Al-4V with nearly identical peak lattice temperatures. That cancellation is a credible, useful insight, and the sensitivity analysis explains the mechanism.\n\nThe paper does several things well. The Ti resistivity cross-check against independent abinit electron–phonon calculations and experiment is solid, including the high-temperature saturation, which they trace to mean-free-path saturation in the CPA–AAM. The GTC vs Drude comparison at low Te is transparent and appropriate. They also flag their own domain-of-validity limit above ~1.3Tm, which is honest.\n\nThe biggest soft spot is the continuation of κe(Te,Tl) beyond the computed map. The map stops at Tl = 2500 K, yet the phase-explosion TTM run reaches lattice temperatures of 8–9.6 kK. The paper never states how κe is continued there — clamped, extrapolated, liquid-metal estimate? Their disclaimer in Sec. II C (“extend the present treatment beyond its strict domain of validity”) shows awareness, but not which curve actually entered Fig. 4(b). Because the headline 19% overshoot of the low-T Drude model is drawn from exactly that regime, the number is not pinned down until the continuation is specified and sensitivity to it is checked. This is fixable and does not invalidate the method.\n\nSecond, the alloy κe map is not validated against any independent calculation or experiment. The Ti cross-check gives indirect confidence, but CPA is most stressed in a disordered alloy. Error bars would help. Third, no data or code are shipped; “available upon request” is weak for a paper whose main product is a parameter set. The ablation-threshold comparison across different experiments is treated cautiously, which is fair.\n\nOverall, the computational pipeline is sound, the dataset is useful to anyone modeling ultrashort-pulse laser ablation of Ti-6Al-4V, and the compensation story is likely right. The 19% functional-form claim is the one quantitative conclusion that could shift with better handling of the high-Tl regime. Fix that, supply the data, and this is a solid paper. I’d send it to a serious referee.","headline":"Useful alloy-specific transport parameter set for Ti-6Al-4V with a credible G-κe compensation story, but the headline 19% claim sits on an unspecified high-temperature extrapolation.","tokens_in":21478,"tokens_out":4532,"would_cite":true,"duration_ms":46087,"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":"In simulations of ultrashort-pulse ablation of Ti-6Al-4V, the functional form of the electronic thermal conductivity shifts peak lattice temperature by up to 19% — more than the elemental-Ti versus alloy-specific distinction.","keywords":["Ti-6Al-4V","ultrashort-pulse laser ablation","two-temperature model","electronic thermal conductivity","electron-phonon coupling","Kubo-Greenwood","coherent potential approximation","Mott-Ioffe-Regel limit"],"falsifier":"Time-resolved measurement of the transient lattice temperature of Ti and Ti-6Al-4V after a 300-fs excitation pulse (for example, ultrafast electron or X-ray diffraction), compared side by side with two-temperature-model predictions using the first-principles κe map versus the Drude form; if the Drude form matches the measured temperature better, the central 19% claim would be contradicted.","tokens_in":20465,"feed_emoji":"⚡","tokens_out":6145,"duration_ms":87073,"temperature":0.7,"pith_summary":"This paper tries to establish that the dominant source of error in two-temperature-model predictions of ultrashort-pulse laser ablation of Ti-6Al-4V is not the common practice of substituting elemental-titanium parameters, but the functional form chosen for the electronic thermal conductivity. It computes temperature-dependent electronic transport for both hcp Ti and the alloy from first principles, treating chemical and thermal disorder on equal footing. The central numerical discovery is that alloy-specific parameters and elemental-titanium parameters give nearly the same peak lattice temperatures (within about 1.4%), because reduced electron-phonon coupling offsets a 6.4-fold reduction in thermal conductivity. Replacing the first-principles conductivity with the widespread low-temperature Drude form shifts the peak lattice temperature by up to 19%, making the transport-model choice more consequential than the material-property choice. A sympathetic reader would care because most existing simulations use that Drude form, meaning their predictions carry a systematic error larger than the material-property uncertainty they are trying to correct.","feed_headline":"Conductivity model choice shifts ablation predictions by 19%","feed_subtitle":"For the workhorse titanium alloy, the conductivity model—not the material parameters—dominates the error.","key_machinery":"The central object is the electronic thermal conductivity map κe(Te, Tl), computed via Kubo-Greenwood linear response within a Korringa-Kohn-Rostoker Green's function scheme that combines the coherent potential approximation for chemical disorder with an alloy-analogy model for thermal disorder. From the energy-resolved conductivity σ(E), generalized Onsager transport coefficients yield κe(Te, Tl) including the Seebeck correction. These maps feed a two-temperature model whose logarithmic sensitivity analysis isolates the influence of each parameter; the key mechanism is that the electron-phonon coupling G and the conductivity prefactor κe0 have opposite-sign sensitivities that cancel in the","core_discovery":"The paper shows that for the workhorse titanium alloy Ti-6Al-4V, the accuracy of two-temperature-model simulations of ultrashort-pulse ablation is governed more by the functional form of the electronic thermal conductivity than by whether the parameters are taken from elemental Ti or from the alloy itself. First-principles Kubo-Greenwood calculations yield thermal conductivity maps that saturate and then decrease with electronic temperature, peaking near 2.97 kW m−1 K−1 for Ti but only 0.47 kW m−1 K−1 for Ti-6Al-4V. Despite that factor-of-6.4 difference, the alloy's lower electron-phonon coupling compensates, so peak lattice temperatures from the two parameter sets differ by only about 1.4%.","pith_inferences":["If the 19% shift is representative, past and existing two-temperature-model studies of Ti-6Al-4V that adopt the linear Te/Tl conductivity may carry a systematic overestimate of peak lattice temperature, although relative comparisons between materials could remain valid since both are affected similarly.","The same first-principles machinery could be applied to other commercial alloys; a useful test would be whether the G-κe compensation observed here is a general feature of strongly chemically disordered metals or specific to Ti-6Al-4V.","Because the paper uses peak lattice temperature as a proxy for ablation threshold rather than an explicit spallation or phase-explosion criterion, combining the new parameter set with a hydrodynamic or molecular-dynamics ablation model and comparing crater-depth data would directly test the proxy's validity.","A clean falsifying experiment: measure the ablation-threshold ratio of Ti to Ti-6Al-4V at nanosecond pulses; if it does not approach the predicted ~1.09, the breakdown of G-κe compensation at long pulse durations is wrong."],"forward_implications":["A complete first-principles parameter set for hcp Ti-6Al-4V—electronic heat capacity, electron-phonon coupling, and κe(Te, Tl)—is now available for two-temperature-model simulations without substituting elemental-titanium values.","Using alloy-specific parameters instead of elemental Ti changes predicted peak lattice temperatures by only about 1.4% at 300 fs, consistent with experimental ablation thresholds that differ by about 3%.","Replacing first-principles κe with the low-temperature Drude form raises the peak lattice temperature by 7% at the melting threshold and 19% at phase explosion, so the functional form dominates predictive error.","The near-identical lattice response of Ti and Ti-6Al-4V stems from opposite-sign sensitivities of G and κe that cancel; this compensation is strongest for sub-picosecond pulses and weakens at longer durations, predicting about a 9% higher ablation threshold for Ti in the nanosecond limit.","The experimentally observed high-temperature resistivity saturation of Ti is reproduced by the calculation and traced to loss of quasiparticle coherence, with mean free paths approaching the Mott-Ioffe-Regel limit."],"fun_headline_variants":["Model form beats titanium alloy choice in ablation code","Thermal conductivity model drives laser ablation predictions","Alloy vs element matters less than transport model shape","For Ti-6Al-4V, conductivity model dominates ablation error"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The transport map is trusted at lattice temperatures above about 1.3 times the melting point, where the hcp crystal formalism formally breaks down; the paper linearly extrapolates the conductivity input to 2500 K and then uses it up to the phase-explosion regime near 8000–9600 K.","fun_headline_variants_meta":{"raw":{"variants":["Model form beats titanium alloy choice in ablation code","Thermal conductivity model drives laser ablation predictions","Alloy vs element matters less than transport model shape","For Ti-6Al-4V, conductivity model dominates ablation error"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000185,"raw_usage":{"total_tokens":1245,"prompt_tokens":916,"completion_tokens":329,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":660,"completion_tokens_details":{"reasoning_tokens":265}},"tokens_in":660,"tokens_out":329,"duration_ms":3780,"temperature":1.0,"reasoning_tokens":265,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T16:02:34.358970+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Time-resolved measurement of the transient lattice temperature of Ti and Ti-6Al-4V after a 300-fs excitation pulse (for example, ultrafast electron or X-ray diffraction), compared side by side with two-temperature-model predictions using the first-principles κe map versus the Drude form; if the Drude form matches the measured temperature better, the central 19% claim would be contradicted.","supporting_citations":[],"review_version":1}