{"id":"7605f3cd-2752-4d6c-bec1-cde51dd45e6f","arxiv_id":"2607.01673","paper_version":1,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Stoichiometric LiNbO3 exhibits intrinsically low room-temperature thermal conductivity because stronger anharmonicity and larger scattering phase space suppress phonon lifetimes by 1-2 orders of magnitude relative to high-conductivity benchmarks like cBAs.","lead":"The paper measures thermal conductivity of stoichiometric lithium niobate with laser pump-probe experiments and machine-learned simulations, reporting values orders of magnitude below silicon due to strong phonon anharmonicity. Smart generalists might read it to understand heat management limits in lithium niobate electro-optic devices.","discovery_kind":"new_application","skeptic_critique":{"model":"grok-4.3","headline":"Accuracy of ML potential for anharmonic interactions lacks direct DFT benchmarks independent of κ match","rationale":"The reader's weakest_assumption correctly isolates the single point where the argument is least secure: the MLIP must faithfully reproduce anharmonic physics for the origin claim to hold. No other internal inconsistency appears in the abstract or described results; the experimental-simulation agreement on total κ is consistent but insufficient to validate the decomposition.","tokens_in":1866,"tokens_out":343,"duration_ms":21334,"concrete_test":"Extract third-order interatomic force constants from the published MLIP for a 2×2×2 supercell of sLN; recompute them with direct DFT (same functional and cutoff) on identical displacements; if the root-mean-square difference in |Φ_3| exceeds 20% or alters the dominant scattering channels, the anharmonicity-based explanation is unreliable.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim attributes low κ to stronger anharmonicity and larger scattering phase space (suppressing lifetimes 1-2 orders of magnitude) by comparing harmonic properties (C_v, v_g) favorably to cBAs while showing enhanced anharmonicity in sLN. This decomposition is obtained entirely from the MLIP-based simulations. The measured κ and T^{-α} (α≈1) scaling confirm intrinsic scattering but do not independently verify the microscopic breakdown. If the MLIP systematically overestimates third-order force constants or phase space (even while reproducing total κ), the origin attribution fails. The paper states the potential is not tuned to κ, yet no section details quantitative validation of anharmonic quantities against direct DFT.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript reports FDTR measurements of thermal conductivity in stoichiometric LiNbO3 (sLN) that agree with MLIP-based atomistic simulations at room temperature. The temperature scaling κ ~ T^{-α} with α ≈ 1 indicates intrinsic phonon-phonon scattering. By comparing harmonic properties (C_v, v_g) to cBAs and showing enhanced anharmonicity plus larger scattering phase space in sLN, the authors attribute the low κ to 1-2 order-of-magnitude suppression of phonon lifetimes, yielding a maximum MFP of ~140 nm and size effects in films thinner than 1 μm.","tokens_in":1992,"tokens_out":538,"duration_ms":20074,"significance":"The independent experimental FDTR data confirming intrinsically low κ in sLN, together with the T^{-α} scaling, establish a useful baseline for thermal management in LN-based electro-optic devices. The decomposition separating harmonic from anharmonic contributions is a strength if the MLIP is shown to be accurate for third-order interactions; the work would then supply falsifiable microscopic predictions for size-dependent transport.","major_comments":[{"comment":"The manuscript provides no quantitative, independent DFT benchmarks of anharmonic quantities (third-order force constants, mode-resolved scattering rates, or phase-space volumes) from the MLIP. Because the central claim that low κ originates from stronger anharmonicity and larger phase space (rather than harmonic properties) rests entirely on the MLIP-derived lifetimes and MFPs, the absence of such benchmarks makes the microscopic origin attribution load-bearing and unverified.","section":"Methods (MLIP training/validation) or Results (phonon lifetime comparison)"},{"comment":"The statement that harmonic properties are 'not responsible' for the low κ relies on the MLIP reproducing C_v and v_g that are comparable to or higher than cBAs; however, no sensitivity analysis or error propagation from the MLIP fitting is shown to confirm that small errors in these quantities would not alter the conclusion.","section":"Results (comparison to cBAs)"}],"minor_comments":[{"comment":"The abstract claims agreement between measured and simulated κ but does not report the numerical values or relative difference; this should be stated explicitly.","section":"Abstract"},{"comment":"The maximum MFP of approximately 140 nm is stated without reference to the specific phonon mode or cumulative MFP plot from which it is extracted.","section":"Results (MFP analysis)"}],"recommendation":"major_revision","confidential_remarks":null},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for the constructive comments and for recognizing the significance of the experimental confirmation of intrinsically low thermal conductivity in stoichiometric LiNbO3. We address each major comment below.","responses":[{"response":"We agree that direct DFT benchmarks of third-order quantities would provide additional reassurance. The MLIP was trained on DFT-derived energies, forces, and stresses from configurations that sample anharmonic behavior, and its accuracy for thermal transport—including anharmonic effects—is validated by quantitative agreement with independent FDTR measurements of both the room-temperature κ value and the T^{-α} scaling. Discrepancies in third-order interactions would have produced clear mismatches with experiment. We will expand the Methods section with further details on the training set composition and any available cross-checks for second-order properties, but full mode-resolved third-order DFT benchmarks remain computationally prohibitive for the supercell sizes required.","revision_made":"partial","referee_comment":"The manuscript provides no quantitative, independent DFT benchmarks of anharmonic quantities (third-order force constants, mode-resolved scattering rates, or phase-space volumes) from the MLIP. Because the central claim that low κ originates from stronger anharmonicity and larger phase space (rather than harmonic properties) rests entirely on the MLIP-derived lifetimes and MFPs, the absence of such benchmarks makes the microscopic origin attribution load-bearing and unverified."},{"response":"We acknowledge that an explicit sensitivity analysis would strengthen the presentation. The harmonic quantities (heat capacity and group velocities) are obtained from the dynamical matrix and are less sensitive to fitting errors than lifetimes. The factor of 10–100 suppression in lifetimes dominates the conductivity difference; even a 20% uncertainty in C_v or v_g would not change the conclusion that anharmonicity and phase space are the primary origin. In the revised manuscript we will add a short paragraph discussing the magnitude of MLIP fitting errors on second-order properties and confirming that the lifetime contrast remains decisive.","revision_made":"yes","referee_comment":"The statement that harmonic properties are 'not responsible' for the low κ relies on the MLIP reproducing C_v and v_g that are comparable to or higher than cBAs; however, no sensitivity analysis or error propagation from the MLIP fitting is shown to confirm that small errors in these quantities would not alter the conclusion."}],"tokens_in":1577,"tokens_out":500,"duration_ms":21039,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The punchline is that this work delivers the first combined experimental and simulation study of thermal transport in stoichiometric lithium niobate, showing low intrinsic conductivity due to anharmonicity.\n\nThe new part is the FDTR measurement that fills the data gap, plus the comparison to cBAs that rules out harmonic properties as the culprit. The simulations match the measured room-temperature value and the temperature scaling with alpha near 1 supports phonon-phonon scattering dominance. The size effect in thin films is a useful practical takeaway.\n\nIt does well on establishing the low kappa as intrinsic and on using the benchmark material to isolate the mechanism.\n\nThe soft spot is in the ML potential. The claim about stronger anharmonicity and larger phase space suppressing lifetimes by 1-2 orders rests on that potential. The stress-test concern holds: there is no direct DFT validation of the anharmonic interactions shown, so if the potential overestimates those, the origin story could be off even while total kappa matches. The paper notes it's not tuned to kappa, but independent checks would strengthen it.\n\nThis paper is for people in integrated photonics and nonlinear optics who care about thermal management in LN devices. The experimental data is the core value, and the analysis adds context.\n\nIt deserves serious refereeing because it addresses a real scarcity in the literature with both measurement and theory.","headline":"New FDTR data shows low intrinsic kappa in sLN with anharmonicity as the driver via ML sims, but the potential lacks independent anharmonic DFT checks.","tokens_in":2490,"tokens_out":354,"would_cite":true,"duration_ms":27120,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"Stoichiometric lithium niobate conducts heat poorly because its phonons scatter far more strongly than in high-conductivity materials.","keywords":["lithium niobate","thermal conductivity","phonon anharmonicity","stoichiometric LiNbO3","frequency-domain thermoreflectance","phonon lifetimes","mean free path","size effects"],"falsifier":"A direct measurement of phonon lifetimes in stoichiometric lithium niobate that are not one to two orders of magnitude shorter than those in cubic boron arsenide would falsify the proposed microscopic origin.","tokens_in":2767,"feed_emoji":"","tokens_out":504,"duration_ms":20020,"temperature":0.7,"pith_summary":"The paper measures the thermal conductivity of stoichiometric lithium niobate with frequency-domain thermoreflectance and obtains values orders of magnitude below those of silicon and other semiconductors. Atomistic simulations reproduce the measurements and show that phonon heat capacity and group velocities are comparable to or higher than in ultrahigh-conductivity cubic boron arsenide. The decisive difference is stronger anharmonicity together with a larger scattering phase space, which shorten phonon lifetimes by one to two orders of magnitude. This produces a maximum mean free path near 140 nm and makes conductivity in thin films fall to half the bulk value already at 10 nm thickness.","feed_headline":"Short phonon lifetimes cause low thermal conductivity in lithium niobate","feed_subtitle":"Values are orders of magnitude below silicon; conductivity halves at 10 nm film thickness due to strong anharmonicity.","key_machinery":"stronger anharmonicity and larger scattering phase space that suppress phonon lifetimes relative to cubic boron arsenide","core_discovery":"The low thermal conductivity of stoichiometric LiNbO3 originates from substantially stronger anharmonicity and larger scattering phase space, which suppress phonon lifetimes by 1-2 orders of magnitude, leading to a maximum phonon mean free path of approximately 140 nm. Harmonic properties such as heat capacity and group velocities are not responsible, as they are either higher than or comparable to those in cubic boron arsenide. The temperature dependence follows a T to the minus one scaling, confirming intrinsic phonon-phonon scattering dominance, and size effects appear in films below 1 micrometer.","pith_inferences":[],"forward_implications":[],"fun_headline_variants":["Anharmonicity suppresses phonon lifetimes in stoichiometric LiNbO3","Larger scattering space dominates low kappa in lithium niobate","LiNbO3 heat capacity rivals cBAs but kappa stays low","Phonon paths limited to 140 nm in stoichiometric LiNbO3"],"cache_read_input_tokens":2112,"weakest_assumption_plain":"The machine-learned interatomic potential accurately reproduces the anharmonic phonon interactions of stoichiometric lithium niobate.","fun_headline_variants_meta":{"raw":{"variants":["Anharmonicity suppresses phonon lifetimes in stoichiometric LiNbO3","Larger scattering space dominates low kappa in lithium niobate","LiNbO3 heat capacity rivals cBAs but kappa stays low","Phonon paths limited to 140 nm in stoichiometric LiNbO3"]},"model":"grok-4.3","cost_usd":0.006995,"raw_usage":{"total_tokens":3321,"prompt_tokens":831,"num_sources_used":0,"completion_tokens":73,"cost_in_usd_ticks":69949500,"prompt_tokens_details":{"text_tokens":831,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":2417,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":831,"tokens_out":73,"duration_ms":19590,"temperature":1.0,"reasoning_tokens":2417,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-03T10:23:04.136116+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A direct measurement of phonon lifetimes in stoichiometric lithium niobate that are not one to two orders of magnitude shorter than those in cubic boron arsenide would falsify the proposed microscopic origin.","supporting_citations":[],"review_version":1}