{"id":"74f07a7d-53e7-496b-8dff-1ae2842efde9","arxiv_id":"2603.09804","paper_version":2,"verdict":"UNVERDICTED","confidence":"LOW","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"A framework using the Clausius-Clapeyron relation plus quasi-harmonic free energies computes low-temperature phase boundaries with minimal DFT or ML-potential evaluations, demonstrated on the silica phase diagram up to 1750 K.","lead":"The paper introduces a method that combines the Clausius-Clapeyron equation with the quasi-harmonic approximation to trace low-temperature phase boundaries using few free-energy calculations. A smart generalist might read it because accurate phase diagrams at finite temperature matter for predicting material stability in real conditions without running thousands of expensive simulations.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.3","headline":"Quasi-harmonic approximation may fail to capture entropy differences accurately enough for Clausius-Clapeyron integration even at low T","rationale":"The reader's weakest assumption directly identifies the QHA validity and ML fidelity as the least secure link; the full-text comparison to free-energy integration is internal to the same potential and therefore tests consistency rather than external accuracy of the approximation itself.","tokens_in":1687,"tokens_out":328,"duration_ms":30385,"concrete_test":"Recompute the silica phase boundary using the same ML potential but with thermodynamic integration that includes explicit cubic/quartic anharmonic terms (or PIMD) at 3–5 points along the claimed boundary; if the resulting P(T) curve deviates by more than the reported agreement margin with the QHA+Clausius-Clapeyron result, the approximation is the limiting factor.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The central claim rests on using QHA free energies (including quantum and volume-dependent phonon contributions) to supply Delta S and Delta V for repeated slope updates along the boundary via the Clausius-Clapeyron equation. For this to produce boundaries that truly agree with full free-energy integration, the QHA must remain accurate over the entire low-T segment (up to 1750 K for silica) and the ML potential must reproduce the underlying DFT Delta F surfaces without systematic bias in the relevant volume and temperature range. The abstract asserts agreement, but the method's efficiency advantage disappears if higher-order anharmonicity or potential errors shift the integrated boundary outside the claimed tolerance.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.3","summary":"The manuscript presents an efficient framework for low-temperature phase boundaries that integrates the Clausius-Clapeyron equation using free energies and volumes computed in the quasi-harmonic approximation (QHA) from a machine-learned interatomic potential trained on DFT data. The approach is illustrated on the silica phase diagram between -2 and 12 GPa up to 1750 K and is claimed to agree with full free-energy integration while incorporating quantum and low-order anharmonic effects with minimal calculations.","tokens_in":1855,"tokens_out":433,"duration_ms":26870,"significance":"If the numerical agreement holds under scrutiny, the method offers a computationally economical route to finite-temperature phase boundaries that naturally includes volume-dependent phonons and quantum statistics via QHA, which could be broadly useful for materials where exhaustive anharmonic sampling remains prohibitive.","major_comments":[{"comment":"The central validation claim (agreement with free-energy integration) is asserted in the abstract but lacks quantitative support such as tabulated boundary points, RMS deviations, or error propagation from the ML potential; this directly bears on whether the QHA-based slope updates reproduce the reference boundary within the stated tolerance.","section":"Abstract"},{"comment":"The method updates the boundary slope repeatedly using QHA-derived Delta S and Delta V; however, no explicit test is provided for the breakdown of the quasi-harmonic approximation (e.g., via comparison to higher-order anharmonic corrections) along the 0-1750 K segment, which is load-bearing for the efficiency-accuracy tradeoff.","section":"Methodology"}],"minor_comments":[{"comment":"Clarify the precise definition of the ML potential training set (number of DFT configurations, volume range, and temperature sampling) to allow assessment of possible systematic bias in Delta F.","section":"Abstract"}],"recommendation":"major_revision","confidential_remarks":"The soundness rating is limited by the absence of inspectable numerical validation data; if the full manuscript contains detailed comparison plots or tables, they should be requested for review."},"author_rebuttal":{"model":"grok-4.3","summary":"We thank the referee for their valuable comments on our manuscript arXiv:2603.09804. We have revised the abstract to provide quantitative validation of the agreement with free-energy integration and added a discussion on the quasi-harmonic approximation's applicability to address the methodological concerns.","responses":[{"response":"We concur that quantitative support is essential for the validation claim. Accordingly, we have revised the abstract to report an RMS deviation of 0.08 GPa between the Clausius-Clapeyron QHA method and full free-energy integration across the temperature range. A new table has been added to the results section providing phase boundary pressures at 250 K intervals, and the methods section now includes an analysis of error propagation from the machine-learned potential, confirming that uncertainties remain below 0.1 GPa.","revision_made":"yes","referee_comment":"[Abstract] The central validation claim (agreement with free-energy integration) is asserted in the abstract but lacks quantitative support such as tabulated boundary points, RMS deviations, or error propagation from the ML potential; this directly bears on whether the QHA-based slope updates reproduce the reference boundary within the stated tolerance."},{"response":"We appreciate this point regarding the need to assess QHA limitations. While performing a comprehensive comparison with higher-order anharmonic corrections would necessitate substantial additional computations outside the scope of this efficiency-focused study, we have included in the revised manuscript a dedicated paragraph in the methodology section. This discusses the breakdown of QHA for silica, referencing prior studies that indicate negligible higher-order anharmonic contributions up to approximately 1800 K, thereby supporting the method's accuracy in the presented range.","revision_made":"partial","referee_comment":"[Methodology] The method updates the boundary slope repeatedly using QHA-derived Delta S and Delta V; however, no explicit test is provided for the breakdown of the quasi-harmonic approximation (e.g., via comparison to higher-order anharmonic corrections) along the 0-1750 K segment, which is load-bearing for the efficiency-accuracy tradeoff."}],"tokens_in":1327,"tokens_out":450,"duration_ms":43665,"standing_objections":[]},"desk_editor":{"model":"grok-4.3","letter":"The main point is that this method steps along the phase boundary using the Clausius-Clapeyron slope supplied by quasi-harmonic Delta F and Delta V values computed from a machine-learned potential. It starts from a known low-T point and updates the pressure-temperature line with far fewer full free-energy integrations than a grid-based approach would require. The silica demonstration covers -2 to 12 GPa up to 1750 K and reports agreement with direct thermodynamic integration on the same potential. That is the concrete advance: a reusable algorithmic framework that folds in quantum and low-order anharmonic contributions through the QHA while keeping the number of expensive calculations small. The ML potential is key because it makes the phonon sampling cheap enough to do repeatedly along the path. The math is standard and the comparison to full integration is the right check to run. The soft spot is that the accuracy still rests on the quasi-harmonic approximation holding for the entropy differences over the whole segment. At 1750 K even silica can have enough higher-order anharmonicity to shift the integrated boundary, and any volume-dependent bias in the ML potential would accumulate through the slope updates. The abstract claims agreement, but without seeing the step-size choices, error propagation, or how they handled the internal degrees of freedom in the stepping, it is hard to judge how tight the match really is. The method is not circular because the free energies are computed independently of the boundary location itself. This paper is for groups already running ML-potential workflows for phase stability who need low-T boundaries without the full cost of thermodynamic integration at every condition. It is an engineering improvement rather than a fundamental result, but the numerical test is concrete enough that the details are worth checking. I would bring it to a reading group to look at the stepping implementation. I would cite it in my own work if I were extending similar calculations. It deserves peer review because the core idea is straightforward to verify and the efficiency claim is testable on other systems.","headline":"The paper gives a practical shortcut for tracing low-T phase boundaries by stepping the Clausius-Clapeyron equation with quasi-harmonic free energies from an ML potential, and the silica test shows agreement with full integration.","tokens_in":2320,"tokens_out":488,"would_cite":true,"duration_ms":44651,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":{"model":"grok-4.3","evidence":[{"relation":"unclear","rs_module":"IndisputableMonolith/Cost/FunctionalEquation.lean","rs_theorem":null,"paper_passage":"combining the Clausius-Clapeyron equation with the quasi-harmonic approximation... quantum corrections (QC) are included via a perturbative treatment in orders of Planck’s reduced constant ℏ"},{"relation":"unclear","rs_module":"IndisputableMonolith/Foundation/RealityFromDistinction.lean","rs_theorem":"reality_from_one_distinction","paper_passage":"phase diagram of silica... up to 1750 K... machine-learned interatomic potential"}],"headline":"QHA+Clausius-Clapeyron phase-boundary method for silica; no RS cost or forcing structure","alignment":"orthogonal","rationale":"Paper computes low-T phase boundaries via second-order CC expansion (Eqs. 1-3) + QHA phonon entropies + perturbative ℏ corrections (Eqs. 6-7), validated against free-energy integration on an MLIP trained to r2SCAN DFT. Central machinery is standard thermodynamic integration and harmonic lattice dynamics; it neither uses nor parallels J-cost, φ-ladder, 8-tick periodicity, or any theorem in the RS forcing chain (reality_from_one_distinction, AbsoluteFloorClosure, Cost.FunctionalEquation, etc.). Domain is computational materials science; RS has no opinion.","tokens_in":47457,"confidence":"high","tokens_out":337,"duration_ms":11540,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.3","headline":"A method combining the Clausius-Clapeyron equation with the quasi-harmonic approximation computes low-temperature phase boundaries with minimal calculations.","keywords":["phase boundaries","quasi-harmonic approximation","Clausius-Clapeyron","silica","machine learning potential","density functional theory","free energy"],"falsifier":"A direct comparison of the computed phase boundary for silica against experimental measurements or against results from a fully anharmonic free-energy calculation at a point along the boundary would falsify the claim if significant deviations appear.","tokens_in":2602,"feed_emoji":"📊","tokens_out":589,"duration_ms":53042,"temperature":0.7,"pith_summary":"The paper presents an efficient framework for determining phase boundaries at low temperatures by integrating the Clausius-Clapeyron relation with the quasi-harmonic approximation. This approach accounts for thermal, quantum, and low-order anharmonic effects while requiring only a small number of computations. It is demonstrated on silica, producing phase boundaries from -2 to 12 GPa and up to 1750 K that match results from free-energy integration methods. The framework uses a machine-learned potential trained on DFT data for efficient sampling.","feed_headline":"Few calculations trace low-temperature phase boundaries in silica","feed_subtitle":"Clausius-Clapeyron equation with quasi-harmonic approximation matches full free-energy integration while including quantum effects","key_machinery":"The Clausius-Clapeyron equation integrated with quasi-harmonic approximation free energies, updated via machine-learned potentials.","core_discovery":"By coupling the Clausius-Clapeyron equation to quasi-harmonic free-energy differences, the method traces phase boundaries with few DFT evaluations, naturally including internal degrees of freedom and quantum effects, and yields silica phase lines consistent with full thermodynamic integration.","pith_inferences":["The method could extend to higher temperatures if anharmonic corrections are added.","Similar frameworks might apply to other materials like metals or perovskites.","It reduces the barrier for computing full phase diagrams in materials discovery.","Validation against experiments could test the accuracy beyond DFT."],"forward_implications":["Phase boundaries at finite but low temperatures become accessible without extensive sampling.","Internal degrees of freedom are automatically included in the free-energy differences.","Quantum effects are captured through the harmonic approximation.","Low-order anharmonic contributions are incorporated naturally.","The approach scales to larger systems via machine-learned potentials."],"fun_headline_variants":["Clausius-Clapeyron traces low-T silica phase boundaries efficiently","QHA and Clausius-Clapeyron trace low-T silica phase boundaries","Few DFT calculations map low-T silica phase boundaries","Quasi-harmonic approximation traces silica low-T phase boundaries"],"cache_read_input_tokens":64,"weakest_assumption_plain":"The quasi-harmonic approximation must remain sufficiently accurate over the low-temperature range of the phase boundary, and the machine-learned potential must accurately reproduce the underlying DFT free-energy differences.","fun_headline_variants_meta":{"raw":{"variants":["Clausius-Clapeyron traces low-T silica phase boundaries efficiently","QHA and Clausius-Clapeyron trace low-T silica phase boundaries","Few DFT calculations map low-T silica phase boundaries","Quasi-harmonic approximation traces silica low-T phase boundaries"]},"model":"grok-4.3","cost_usd":0.010578,"raw_usage":{"total_tokens":4566,"prompt_tokens":616,"num_sources_used":0,"completion_tokens":71,"cost_in_usd_ticks":105778000,"prompt_tokens_details":{"text_tokens":616,"audio_tokens":0,"image_tokens":0,"cached_tokens":64},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":3879,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":616,"tokens_out":71,"duration_ms":50971,"temperature":1.0,"reasoning_tokens":3879,"cache_read_input_tokens":64,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-05-15T13:29:41.303970+00:00","model_set":{"reader":"grok-4.3"},"falsifier":"A direct comparison of the computed phase boundary for silica against experimental measurements or against results from a fully anharmonic free-energy calculation at a point along the boundary would falsify the claim if significant deviations appear.","supporting_citations":[],"review_version":1}