{"id":"5b40b502-a7e7-439e-b48d-f37370a8a1b2","arxiv_id":"2507.14532","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"Uniaxial tensile strain monotonically suppresses the lattice thermal conductivity of cubic boron arsenide by up to about 80%, driven by phonon degeneracy lifting and stronger three- and four-phonon scattering.","lead":"This paper uses first-principles calculations to show that stretching boron arsenide along one direction crushes its ability to conduct heat, cutting thermal conductivity by about 80% at large strain. It explains the drop through broken phonon symmetry and stronger phonon-phonon scattering, and finds the effect is strongest sideways to the stretch.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Strained structures assumed dynamically stable up to 18% strain, but no full-BZ imaginary-mode check is reported; a soft mode would invalidate the central 80% reduction claim.","rationale":"Dynamic stability is the single most load-bearing condition because its failure invalidates the high-strain results outright, whereas other concerns (4ph cutoff radius, q-mesh convergence) would only shift quantitative values. The 80% reduction is already obtained at the 3ph level (b-axis: 2302 to 485 W/mK, Sec. III.C), so the 4ph truncation issue cannot overturn the qualitative central claim. The manuscript's own stress-strain argument addresses homogeneous elastic deformation, not phonon stability, and the 14% dispersion shown is high-symmetry-only. The reader's weakest assumption matches this concern. Since the reader already issued a CONDITIONAL verdict requiring verification, I recommend no change; the proposed phonon-stability test would either resolve the condition or force restricting the claim.","tokens_in":10058,"tokens_out":6834,"duration_ms":81167,"concrete_test":"Recompute the phonon dynamical matrix for the relaxed strained BAs structures at ε = 6%, 10%, 14%, and 18% using the same PBE/520 eV/4×4×4 supercell settings, sampling a uniform q-grid over the full BZ (e.g., 8×8×8 or 16×16×16). If the minimal squared phonon frequency is negative at any strain ≤18%, the structure is dynamically unstable and the κL predictions at and beyond that strain (and the monotonic decrease claim) must be restricted to the stable window; if all frequencies are positive, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim—monotonic suppression of κL under uniaxial [100] tension, reaching ~80% at 18%—requires that the tetragonal strained BAs structures are dynamically stable across the entire 6–18% strain window. The paper justifies the window by the linear stress-strain response (Sec. III.A, Fig. 1a) and notes metallization only at larger strain, but does not report a full-Brillouin-zone imaginary-mode check for the strained cells. Fig. 2c shows phonon dispersions at 14% strain only along high-symmetry paths, which can miss soft modes at generic q-points, and no stability statement is made for 18%. Elastic linearity under homogeneous deformation does not guarantee Born stability or the absence of phonon instabilities. If a soft mode appears within the studied range, the strained reference structure is not the equilibrium one; the reported phonon softening, phase-space enhancement, and 80% suppression would be artifacts of an unstable geometry. This is an internal consistency condition of the harmonic/anharmonic framework, not a disagreement with consensus.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports first-principles DFT+BTE calculations of lattice thermal conductivity (κL) of cubic BAs under uniaxial tensile strain along [100], including both three- and four-phonon scattering. It finds that strain lifts the TA and TO degeneracies and softens the phonon spectrum, increasing scattering phase space and rates, leading to a monotonic reduction of κL with strain and an approximately 80% drop at 18% strain at room temperature. The suppression is stronger perpendicular to the strain direction than along it. The authors attribute this effect to symmetry-breaking that opens new scattering channels, and they support the mechanism with weighted-phase-space decomposition, matrix-element analysis, and small-grain-limit κSG estimates.","tokens_in":10194,"tokens_out":5684,"duration_ms":67652,"significance":"If the predictions hold, this is a useful contribution: it extends strain engineering of κL to anisotropic uniaxial strain, includes four-phonon scattering (important in BAs), identifies a concrete symmetry-based mechanism (degeneracy lifting) with a quantitative decomposition into phase-space and matrix-element contributions, and makes a falsifiable prediction of strong anisotropic suppression. The calculations are parameter-free with respect to the target κL, use an iterative BTE solver with standard DFT settings, and benchmark the stress-strain behavior against diamond. The main caveat is that lattice stability under large strain is asserted rather than directly verified, which affects the reliability of the headline 80% reduction.","major_comments":[{"comment":"The central ~80% reduction claim presumes that all strained structures up to 18% are dynamically stable. The only stability evidence is the linear stress-strain response below 20% (Fig. 1a) and phonon dispersions along high-symmetry paths at 14% (Fig. 2c). Elastic linearity does not rule out imaginary phonon modes at generic q-points, and high-symmetry-path sampling can miss soft-mode instabilities. Please compute full-Brillouin-zone phonon dispersions (or report the lowest phonon frequency on a dense q-grid) for every strained geometry, including 18%, and state explicitly whether any imaginary modes exist. If a soft mode appears, the strained reference is not the equilibrium structure and the predicted κL suppression would need to be reconsidered.","section":"III.A and III.B (Fig. 2c)"},{"comment":"The convergence of the IFC truncation radii and the q-mesh is argued from prior BAs studies at zero strain. Uniaxial strain changes the symmetry to tetragonal and can alter both the range of anharmonic interactions and the required Brillouin-zone sampling. Please report at least one convergence test for a strained cell (e.g., κL versus q-mesh and versus third- and fourth-order interaction radius at ε=14% or 18%) so that the quantitative uncertainty of the ~80% reduction is known.","section":"II (computational parameters)"},{"comment":"The manuscript states that metallization occurs only at larger strain, but it does not report the electronic band gap at the strained geometries. Because the headline number is at ε=18%, please give the band gap at each studied strain and confirm that the 18% structure is still semiconducting; otherwise phonon-electron scattering should be included in the BTE.","section":"III.C"}],"minor_comments":[{"comment":"The sentence 'suggesting that the increased TA scattering rates at this strain level arise from the weakened third-order anharmonicity' is logically inconsistent with the observed decrease in |V(3)|^2; the increased scattering must instead come from the enlarged phase space. Please rephrase.","section":"III.D (inset of Fig. 5e)"},{"comment":"References 19 and 36 are the same publication (B. Wang et al., Phys. Rev. B 106, 184303 (2022)); consolidate them into a single reference.","section":"References"},{"comment":"Please specify whether the strain is engineering strain and whether the stress is the Cauchy stress, since this matters for comparison with experimental loading conditions.","section":"Fig. 1(a)"},{"comment":"The notation for the ± and ±± channel labels is defined only after the equations; consider defining these labels more explicitly at first use to improve readability.","section":"Eqs. (2)-(5)"}],"recommendation":"major_revision","confidential_remarks":"The decisive issue is the missing dynamic-stability check at large strain. If the authors can supply a full-Brillouin-zone imaginary-mode analysis for all strained structures and confirm the semiconducting character at 18%, the paper would be suitable for publication. The duplicated reference and the logical slip in Sec. III.D are trivial to fix. I see no concerns about novelty or citation practices beyond the standard use of methodology references."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Read this one. It is the first systematic uniaxial-strain study of boron arsenide that includes four-phonon scattering, and the central prediction—monotonic suppression of κL under [100] tension, about 80% at 18% strain—is well supported by the internal logic of the calculations. The contrast with the non-monotonic hydrostatic-pressure behavior is striking and makes the paper worth a referee's time.\n\nWhat the paper does well: the DFT+BTE workflow is standard for this material, with a 4x4x4 supercell for IFCs and a 16x16x16 q-mesh for the BTE. The mechanistic decomposition is careful—the authors separate phase space from matrix elements and show that 3ph enhancement comes from both while 4ph enhancement is driven by phase space alone. The branch-resolved analysis tracing the anisotropy to TA degeneracy lifting is persuasive. The comparison to hydrostatic data from the literature is appropriate, and the stress-strain validation against diamond gives some confidence in the elastic regime.\n\nSoft spots, in order of weight. First, the dynamical stability of the strained phases at 6–18% strain is asserted rather than demonstrated. The paper shows phonon dispersions at 14% along high-symmetry paths and notes linear elasticity, but a full-BZ imaginary-mode check is not reported. If a soft mode appears anywhere in that window, the strained reference geometry is not the equilibrium one, and the 80% reduction would be an artifact. That is a verification gap, not a demonstrated error, but it is the load-bearing assumption of the paper and should have been checked. Second, the fourth-order IFC truncation radius (0.35 nm) is not convergence-tested at large strain; it is fixed at a value known for unstrained BAs. Third, no data or code are deposited, which makes the claimed reproducibility hard to verify.\n\nNone of these kill the paper. The physics is internally consistent, and the mechanism is credible. But the headline number should be viewed as provisional until the stability check is done.\n\nWho gets value: phonon transport people and strain-engineering groups. It deserves a serious referee—the claim is important enough to warrant referee time, and the gaps are addressable. If I were handling it, I'd send to review and ask for the stability check, convergence evidence, and a data deposit as revision conditions. I'd cite it once the stability question is settled.","headline":"A solid, mechanistically careful prediction that uniaxial strain suppresses κL in BAs by up to 80%, but the strained structures' dynamical stability is asserted, not checked, and that gap needs closing before the result can be taken as quantitative.","tokens_in":10790,"tokens_out":2937,"would_cite":false,"duration_ms":32790,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["63.20.-e","66.70.-f"],"model":"deepseek-v4-flash","headline":"Applying uniaxial tensile strain along [100] makes the lattice thermal conductivity of cubic boron arsenide fall monotonically, by nearly 80 percent at room temperature under large strain, because the strain lifts TA/TO phonon…","keywords":["boron arsenide","lattice thermal conductivity","uniaxial strain","phonon-phonon scattering","four-phonon scattering","symmetry breaking","phonon degeneracy lifting","first-principles phonon BTE"],"falsifier":"Compute the full phonon dispersion of BAs at 6, 10, 14, and 18 percent uniaxial [100] strain and search the whole Brillouin zone for imaginary frequencies; any soft mode in that range would invalidate the anharmonic force constants and scattering rates that produce the 80 percent suppression.","tokens_in":9781,"feed_emoji":"🔬","tokens_out":7980,"duration_ms":84397,"temperature":0.7,"pith_summary":"This paper argues that stretching cubic boron arsenide along a single crystal direction is a strong and continuously tunable way to lower its exceptionally high lattice thermal conductivity. First-principles phonon Boltzmann transport calculations that include four-phonon scattering show that uniaxial tensile strain along [100] makes $\\kappa_{\\rm L}$ decrease monotonically, with a nearly 80 percent reduction at room temperature under large strain. The mechanism is symmetry breaking: the strain lifts the degeneracy of transverse acoustic and transverse optical phonon branches and softens the phonon spectrum, which opens many new three- and four-phonon scattering channels. The result matters because it points to crystal-symmetry breaking, rather than volume change, as a practical route for engineering heat flow in high-conductivity semiconductors.","feed_headline":"Stretching boron arsenide cuts its heat flow by 80 percent","feed_subtitle":"Tensile strain lifts phonon degeneracies, opening scattering channels that cut lattice heat conduction by 80 percent.","key_machinery":"The load-bearing mechanism is the strain-induced lifting of symmetry-protected phonon degeneracies. In unstrained BAs, two TA branches and two TO branches are degenerate along paths such as $\\Gamma$-$X$ and $\\Gamma$-$L$, protected by the threefold-rotation and mirror symmetries of the cubic lattice. Uniaxial [100] tension breaks those symmetries, splitting the degenerate branches and shifting the acoustic and optical mode frequencies downward. Because phonon-phonon scattering rates are governed by energy and momentum conservation selection rules, the split branches open new scattering channels, especially AAA and AOO processes, and the softened spectrum raises phonon populations, which also enlarges four-phonon phase space. The argument is carried by mode-resolved weighted phase space, a count of available scattering channels satisfying energy and momentum conservation, together with scattering matrix elements that separate the phase-space effect from the change in anharmonic force constants.","core_discovery":"The paper's central claim is that symmetry-breaking uniaxial strain is a distinct and effective control knob for lattice thermal conductivity in cubic boron arsenide. In first-principles phonon Boltzmann transport calculations that include fourth-order anharmonicity, tensile strain applied along the [100] direction makes $\\kappa_{\\rm L}$ decrease monotonically with strain, reaching a nearly 80 percent reduction at 18 percent strain at room temperature along the direction perpendicular to the load. The strain lifts the degeneracy of the two transverse acoustic branches and the two transverse optical branches, which are protected by crystal symmetries, and softens the overall phonon spectrum. These changes enlarge the three- and four-phonon scattering phase space and increase third-order anharmonicity, so both three- and four-phonon scattering rates rise together, in contrast to the competing responses seen under isotropic pressure. The paper also reports that the suppression is markedly anisotropic: heat conduction perpendicular to the strain direction falls far more than along the stretching direction, an effect traced mainly to strain-mediated three-phonon scattering and TA degeneracy lifting.","pith_inferences":["If the degeneracy-lifting picture is right, other cubic semiconductors with degenerate TA/TO branches and a large acoustic-optical gap should show similar strain-driven suppression; screening zinc-blende compounds computationally would reveal which respond most.","Because the strain axis also narrows the electronic band gap, the mechanism offers a way to couple heat-flow suppression with electronic property changes in one device, though the paper stops before metallization.","A membrane or film experiment that measures in-plane and cross-plane thermal conductivity under uniaxial strain could directly test the predicted anisotropy: the perpendicular direction should show the larger drop.","If a phonon instability appears below 18 percent strain, the monotonic suppression would break down; checking the strained phonon dispersion for imaginary modes would settle the practical strain limit."],"forward_implications":["At 18 percent uniaxial tensile strain, the room-temperature lattice thermal conductivity of BAs along the perpendicular direction falls from 1244 to 254 W/mK when four-phonon scattering is included, a reduction of about 80 percent.","The strain response is directional: the perpendicular axis loses far more heat-carrying capacity than the stretching axis, so a single strained sample should display a measurable thermal anisotropy.","Both three- and four-phonon scattering rates grow monotonically with strain, meaning the usual competition between these channels under isotropic pressure is replaced by cooperation under symmetry-breaking strain.","Changes in group velocities and heat capacity contribute at most about a quarter of the drop, so the reduction is primarily an anharmonic-scattering effect rather than a harmonic one."],"supporting_citations":[{"why":"Provides the hydrostatic-pressure baseline showing non-monotonic lattice thermal conductivity in BAs and BSb, which uniaxial strain is contrasted against.","marker":"[14]"},{"why":"Supplies the experimental and computational pressure data for BAs and diamond used to benchmark the stress-strain and pressure calculations.","marker":"[15]"},{"why":"Prior uniaxial-strain study in TiO showing TA splitting and lattice thermal conductivity suppression, the motivation for applying the idea to BAs.","marker":"[20]"},{"why":"Supplies the phonon Boltzmann transport equation framework and third-order anharmonic force constant methodology used to compute thermal conductivity and three-phonon scattering.","marker":"[21]"},{"why":"Provides the four-phonon scattering implementation used to include fourth-order anharmonicity in the thermal conductivity calculations.","marker":"[26]"},{"why":"Establishes that four-phonon scattering significantly affects the lattice thermal conductivity of BAs, justifying its inclusion.","marker":"[33]"},{"why":"Shows that four-phonon processes in BAs are governed by AAOO channels and are not suppressed by the dispersion features that limit three-phonon scattering.","marker":"[34]"},{"why":"Explains the weak three-phonon scattering in BAs from the acoustic-optical gap and acoustic bunching, the baseline that strain modifies.","marker":"[40]"}],"fun_headline_variants":["Stretching boron arsenide quenches heat flow by 80%","Boron arsenide heat flow plummets 80% under uniaxial tension","Uniaxial strain cuts boron arsenide thermal conductivity 80%","Symmetry breaking drives 80% heat-flow drop in boron arsenide","Perpendicular heat flow in boron arsenide falls 80% under tension"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculations assume BAs stays dynamically stable, elastically deformed, and semiconducting up to 18 percent strain, but no imaginary-frequency check of the strained phonon dispersions is reported.","fun_headline_variants_meta":{"raw":{"variants":["Stretching boron arsenide quenches heat flow by 80%","Boron arsenide heat flow plummets 80% under uniaxial tension","Uniaxial strain cuts boron arsenide thermal conductivity 80%","Symmetry breaking drives 80% heat-flow drop in boron arsenide","Perpendicular heat flow in boron arsenide falls 80% under tension"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000634,"raw_usage":{"total_tokens":2943,"prompt_tokens":978,"completion_tokens":1965,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":594,"completion_tokens_details":{"reasoning_tokens":1868}},"tokens_in":594,"tokens_out":1965,"duration_ms":14803,"temperature":1.0,"reasoning_tokens":1868,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:53:30.215102+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the full phonon dispersion of BAs at 6, 10, 14, and 18 percent uniaxial [100] strain and search the whole Brillouin zone for imaginary frequencies; any soft mode in that range would invalidate the anharmonic force constants and scattering rates that produce the 80 percent suppression.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the hydrostatic-pressure baseline showing non-monotonic lattice thermal conductivity in BAs and BSb, which uniaxial strain is contrasted against."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the experimental and computational pressure data for BAs and diamond used to benchmark the stress-strain and pressure calculations."},{"cited_title":"Jin, D.-s","cited_arxiv_id":null,"evidence_quote":"Prior uniaxial-strain study in TiO showing TA splitting and lattice thermal conductivity suppression, the motivation for applying the idea to BAs."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the four-phonon scattering implementation used to include fourth-order anharmonicity in the thermal conductivity calculations."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Shows that four-phonon processes in BAs are governed by AAOO channels and are not suppressed by the dispersion features that limit three-phonon scattering."},{"cited_title":"Lindsay, D","cited_arxiv_id":null,"evidence_quote":"Explains the weak three-phonon scattering in BAs from the acoustic-optical gap and acoustic bunching, the baseline that strain modifies."}],"review_version":1}