{"id":"6beaa7e9-2088-4ecb-87e2-d0207ddc1ce8","arxiv_id":"2604.01471","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":5.5,"correctness_risk":"low","formal_verification":"none","parameter_count":1,"one_line_summary":"TUNA is a specialized open-source quantum chemistry program for atoms and diatomics that derives all properties from numerical energy derivatives via a simple CLI.","lead":"TUNA is an open-source quantum chemistry code restricted to atoms and diatomics that exposes many electronic-structure methods through one uniform command line. It is built so that every property is obtained by numerical differentiation of the energy, making it a compact teaching and method-development platform.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified","rationale":"The reader correctly identifies the strongest claim (Abstract + §2) and the only non-trivial practical assumption (finite-difference stability). That assumption is real but secondary: the paper’s value proposition is the uniform numerical architecture and the open, documented implementation, not a guarantee of machine-precision derivatives for every conceivable diatomic. No internal inconsistency, circular argument or unsupported theoretical leap appears. Significance is modest (pedagogical/software), novelty moderate, correctness risk low, and reproducibility high—exactly the profile that justifies ACCEPT for a methods/software preprint. No adjustment to the verdict is warranted.","tokens_in":9088,"tokens_out":488,"duration_ms":5166,"concrete_test":"Clone the public repository, install via pip, and re-run the exact command given in §2 (TUNA BDE : H H 1.0 : CCSD cc-pVQZ : ZPE VPT2 EXTRAPOLATE). Confirm that the reported BDE is 103.32 kcal mol^{-1} (or within 0.1 kcal mol^{-1}) and that the finite-difference path completes without SCF/CC divergence. This single end-to-end check verifies both reproducibility and the practical adequacy of the chosen stencils for a high-order correlated method.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The paper’s central claim is architectural and pedagogical rather than a new physical derivation: once an energy evaluator exists, diatomic symmetry makes numerical differentiation a uniform route to structures, frequencies, response properties and dynamics. That claim is directly supported by the released MIT-licensed code, the documented CLI, the explicit finite-difference stencils (§2), and a concrete numerical example (CCSD/cc-pVQZ BDE of H2 matching experiment to 0.05 kcal mol^{-1}). The reader’s weakest assumption—that optimised step sizes h plus tightened convergence remain stable for every method and property—is a legitimate practical caveat, but it is not load-bearing for the claim as stated. The paper never asserts that a single universal h works for all regimes; it only asserts that the numerical-differentiation principle yields a transparent, uniform program. Because the contribution is software infrastructure with open code and documentation, residual numerical-edge-case risk does not undermine the architectural argument.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"TUNA is an open-source Python quantum chemistry program restricted to atoms and diatomics. Its design principle is that once an electronic-structure method can evaluate the energy, structural, vibrational, response and dynamical properties are obtained uniformly by numerical differentiation, exploiting diatomic symmetry. The manuscript describes a compact CLI covering SPE, OPT, FREQ, ANHARM, SCAN, MD, BDE, IP and EA; a broad method set (HF, DFT through double hybrids, MP2–MP4 and variants, CCD through CCSDTQ, QCI, CEPA, CIS); many Gaussian basis sets; finite-difference stencils for gradients and higher derivatives; and a concrete H2 BDE example (CCSD/cc-pVQZ + ZPE/VPT2 + extrapolation) that matches experiment to 0.05 kcal mol−1. The code is MIT-licensed, installable from PyPI, and accompanied by theoretical documentation.","tokens_in":9289,"tokens_out":712,"duration_ms":6414,"significance":"If the released code and documentation match the description, TUNA fills a genuine niche: a transparent, uniform teaching and prototyping environment for electronic structure on the simplest chemically interesting systems. The architectural claim—energy evaluation plus numerical differentiation under diatomic symmetry—is sound and is supported by open code, explicit stencils (§2), and a falsifiable numerical demonstration. Strengths that should be credited include the MIT-licensed Python/NumPy/Cython implementation, the consistent CLI, reference implementations of high-order CC and double-hybrid DFT, and the integrated plotting of orbitals, densities and vibrational wavefunctions. Residual numerical caveats about finite-difference step sizes do not undermine the pedagogical and developmental value of the contribution.","major_comments":[],"minor_comments":[{"comment":"§2 states that h has been optimised and that tightened convergence is used for higher derivatives, but no numerical values of h, stencil widths, or convergence thresholds are given. A short table or appendix listing the defaults would improve reproducibility without changing the architectural claim.","section":null},{"comment":"The H2 BDE example is persuasive but solitary. One additional sentence or table entry for a second diatomic (e.g. N2 or HF frequency or BDE) would better illustrate breadth across methods.","section":null},{"comment":"Figure 2 captions and axis labels appear partially corrupted by Unicode/encoding artefacts in the manuscript source; these should be cleaned for the published version.","section":null},{"comment":"§5 asserts competitiveness with established codes via symmetry and Cython integrals; a single wall-time comparison (even for H2 or N2 at a fixed level) would make that claim concrete.","section":null},{"comment":"Table 2 lists exchange and correlation functionals separately; a brief note on which combinations are actually exposed as named XC functionals (and any double-hybrid defaults) would help users.","section":null}],"recommendation":"accept","confidential_remarks":"This is a software/methods paper rather than a new physical result. For a computational-chemistry or software-focused venue it is a clear accept; if the journal’s bar is primarily novel theory, the editor may still find the pedagogical and reference-implementation value sufficient. The reader’s and skeptic’s assessments align with mine: no load-bearing technical flaw."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"Harry Brough’s TUNA is a focused software paper: a Python/NumPy/SciPy (plus Cython integrals) package restricted to atoms and diatomics that implements a wide method stack—HF, DFT through double hybrids, MP2–4 and SCS/DLPNO variants, CCD through CCSDTQ, QCI, CEPA, CIS(+D)—behind a single CLI of the form TUNA [calc] : atoms dist : method basis : [keywords]. The design claim is explicit and honest: once you can evaluate an energy, every property (opt, force, harmonic/anharmonic freq, BDE/IP/EA, finite-field response, AIMD, scans) is obtained by numerical differentiation, exploiting 1-D symmetry so that stencils stay practical. That is what is actually new—the deliberate scope cut plus the uniform numerical layer—not the electronic-structure methods themselves.\n\nWhat it does well is execution and packaging. The H2 BDE example (CCSD/cc-pVQZ + ZPE/VPT2 + extrapolation → 103.32 vs exp 103.27 kcal mol⁻¹) is a clean external check. Orbitals, densities and PES plots are built-in; the manual is described as theory-plus-pasteable inputs; code is MIT, on GitHub and pip-installable as quantumtuna. Citations to Psi4NumPy, HarPy, Crawford projects, PyDFT and ORCA’s simple-input style are appropriate and non-circular. No free parameters are being fitted to data; the only free knobs are the usual finite-difference h and stencil widths, which the paper says were optimised and paired with tighter convergence.\n\nSoft spots are real but proportionate. The paper is a software description, not a methods paper with broad validation tables; numerical stability of the higher stencils across every method and every diatomic is asserted rather than exhaustively demonstrated. Planned features (relativistic, magnetic response, CASSCF) are listed as future work. None of that undercuts the architectural claim or the released artefact.\n\nThis is for chemical-education groups, method developers who want a readable sandbox, and anyone who needs quick, consistent benchmarks on diatomics. It is not a theory advance and will not change industrial practice. I would send it to peer review as a methods/software contribution; the code and documentation already make the central claim checkable. Worth engaging if you teach or prototype on small systems; otherwise a useful reference to keep on the shelf.","headline":"Clean, open diatomic-only QC package whose real novelty is the uniform numerical-diff architecture and one-line CLI, not new theory.","tokens_in":9874,"tokens_out":611,"would_cite":true,"duration_ms":5978,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"One principle turns any energy method into structures, vibrations and response properties for atoms and diatomics.","keywords":["quantum chemistry software","diatomic molecules","numerical differentiation","electronic structure methods","coupled cluster","density functional theory","teaching platform","response properties"],"falsifier":"Run the same high-order coupled-cluster or double-hybrid calculation for a polarisability or anharmonic frequency both with TUNA’s default finite-difference scheme and with a well-validated analytic or much denser numerical reference; a systematic discrepancy larger than chemical or spectroscopic accuracy would falsify the numerical-stability premise.","tokens_in":9979,"feed_emoji":"⚛️","tokens_out":934,"duration_ms":7458,"temperature":0.7,"pith_summary":"TUNA is an open-source quantum chemistry program restricted to atoms and diatomic molecules. The author argues that this narrow scope, combined with a single design rule, produces a program that is simultaneously a teaching tool, a benchmarking environment and a development sandbox. The rule is simple: once a method can evaluate the energy, every structural, vibrational and field-response property is obtained by numerical differentiation on the one-dimensional bond coordinate. Because diatomic symmetry reduces many observables to one or two independent components, carefully chosen finite-difference stencils become practical rather than a last resort. The result is a uniform command-line interface that runs Hartree–Fock, density-functional, Møller–Plesset and high-order coupled-cluster methods, and immediately yields optimisations, anharmonic frequencies, bond-dissociation energies, polarisabilities and molecular-dynamics trajectories for every one of them. Reference implementations and detailed theory documentation make the code an accessible starting point for students and for researchers who want to prototype new algorithms on the simplest instructive molecules.","feed_headline":"One rule turns any energy method into full diatomic properties","feed_subtitle":"TUNA shows that numerical differentiation plus diatomic symmetry yields a transparent teaching and benchmarking code","key_machinery":"Numerical differentiation of the energy under diatomic symmetry: central differences for gradients, five- to nine-point stencils for higher derivatives, with step sizes and convergence thresholds chosen once for the whole program so that every implemented electronic-structure method inherits the full suite of properties automatically.","core_discovery":"The paper establishes that a deliberately restricted molecular domain (atoms and diatomics) plus the uniform principle “energy evaluation implies all properties by numerical differentiation” yields a transparent, consistent electronic-structure program. Once any method supplies an energy, optimised geometries, harmonic and non-perturbative anharmonic frequencies, finite-field response properties, scans and ab initio trajectories become available without method-specific analytic derivatives.","pith_inferences":["The same numerical-differentiation philosophy could be stress-tested on a small set of linear triatomics to measure how quickly the symmetry advantage erodes.","Because every property is a numerical derivative, TUNA provides a ready-made platform for quantifying finite-difference error bars as a function of method and basis set—an analysis rarely published for production codes.","The planned addition of magnetic response and CASSCF will test whether the single-principle design continues to scale once analytic derivatives become far more complicated than the energy itself."],"forward_implications":["Any newly implemented energy method automatically gains geometry optimisations, frequencies, response properties and dynamics without additional analytic-derivative code.","Students can compare restricted Hartree–Fock, DFT, MP2–MP4 and CCSDTQ on the same diatomic with identical syntax and immediately plot orbitals, densities and potential-energy curves.","Benchmarking campaigns on bond breaking, spin contamination and basis-set extrapolation become uniform across a wide ladder of methods.","Developers can prototype improved algorithms or new functionals on the simplest instructive systems and obtain full property suites at once.","Basis-set extrapolation applied to the energy automatically propagates to approximate complete-basis response properties."],"fun_headline_variants":["One energy rule unlocks all diatomic properties by numerical differentiation","Energy evaluation alone yields geometries frequencies and trajectories","Diatomic symmetry plus numerical differentiation builds complete properties","Any energy method now supplies full properties for atoms and diatomics","Restricted domain turns simple energy codes into teaching benchmarking tools"],"cache_read_input_tokens":128,"weakest_assumption_plain":"The claim rests on the premise that a single set of carefully chosen finite-difference step sizes and tightened energy thresholds remains numerically stable and accurate for every method and every property across the diatomics the program targets.","fun_headline_variants_meta":{"raw":{"variants":["One energy rule unlocks all diatomic properties by numerical differentiation","Energy evaluation alone yields geometries frequencies and trajectories","Diatomic symmetry plus numerical differentiation builds complete properties","Any energy method now supplies full properties for atoms and diatomics","Restricted domain turns simple energy codes into teaching benchmarking tools"]},"model":"grok-4.5","effort":"low","cost_usd":0.00679,"raw_usage":{"total_tokens":1663,"prompt_tokens":701,"num_sources_used":0,"completion_tokens":81,"cost_in_usd_ticks":67900000,"prompt_tokens_details":{"text_tokens":701,"audio_tokens":0,"image_tokens":0,"cached_tokens":256},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":881,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":701,"tokens_out":81,"duration_ms":7548,"temperature":1.0,"reasoning_tokens":881,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-14T19:56:25.179292+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"Run the same high-order coupled-cluster or double-hybrid calculation for a polarisability or anharmonic frequency both with TUNA’s default finite-difference scheme and with a well-validated analytic or much denser numerical reference; a systematic discrepancy larger than chemical or spectroscopic accuracy would falsify the numerical-stability premise.","supporting_citations":[],"review_version":1}