{"id":"b54eabc4-cc49-4d05-a749-529740f8d321","arxiv_id":"2607.06414","paper_version":2,"verdict":"ACCEPT","confidence":"HIGH","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":2,"one_line_summary":"Optimizing ring-polymer normal-mode frequencies to exact harmonic radii of gyration produces Eco path integrals that converge like fourth-order Suzuki-Chin methods at pure Trotter cost.","lead":"The authors optimize free normal-mode frequencies of the ring-polymer spring matrix so that harmonic radii of gyration match exact quantum results up to a cutoff frequency. This yields path-integral simulations that converge like fourth-order methods at ordinary second-order cost, with almost no code changes.","discovery_kind":"new_method","skeptic_critique":{"model":"grok-4.5","headline":"No significant objection identified","rationale":"The reader correctly isolates the only modelling assumption that is not rigorously guaranteed (harmonic R^{2} fit improves anharmonic observables). That assumption is transparent, the optimisation is cheap and reproducible (SI code), and the numerical evidence on two condensed-phase systems with machine-learned and empirical potentials is consistent and favourable. Because the paper never claims a general error bound or universality beyond the tested regime, the empirical character of the assumption does not undermine the concrete claim of practical superiority over Trotter and competitive performance versus SC. No stronger technical flaw (loss of positive-definiteness, invalid estimators, non-monotonic convergence, etc.) is present. The verdict therefore remains ACCEPT.","tokens_in":18042,"tokens_out":482,"duration_ms":23367,"concrete_test":"Re-optimize the Eco frequencies for ice with ω_max deliberately set 20 % below and 20 % above the physical O–H stretch (3200 cm⁻¹ and 4800 cm⁻¹) and recompute the P=32 kinetic and potential energies of Fig. 4; if either energy shifts by more than the reported statistical error relative to the original ω_max=4000 cm⁻¹ result, the robustness of the harmonic-fit premise would be weakened.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Eco matches projected-Hessian SC accuracy at pure Trotter cost, reusing estimators and requiring only a few code changes) is internally consistent and empirically supported. The free parameters of any zero-centroid symmetric circulant spring matrix are correctly identified (Sec. II.A, Eqs. 5–13), the least-squares objective for the harmonic R^{2}(ω) (Eqs. 14–17) is well-posed, and the resulting Eco frequencies preserve the structure that validates the standard thermodynamic and virial estimators (Sec. II.C). The transfer of that harmonic fit to anharmonic ice and MOF-5 is only empirical, yet the reported log-log slopes, radial-distribution errors, lattice constants and heat capacities all improve monotonically and match or exceed SC per force evaluation. No hidden inconsistency or counter-example appears in the mathematics or the two non-trivial applications.","agreement_with_reader":"agree"},"referee_report":{"model":"grok-4.5","summary":"The manuscript introduces economised (Eco) path integrals for imaginary-time PIMD. Starting from the observation that any real symmetric circulant spring matrix K with zero centroid eigenvalue is diagonalised by the same bead-to-normal-mode matrix and leaves floor(P/2) free non-zero frequencies, the authors optimise those frequencies by least-squares fitting of the harmonic ring-polymer radius of gyration R^{2}(ω) (and equivalently the quantum energy shift) over 0 ≤ ω ≤ ω_max. The resulting Eco frequencies preserve the structure that validates the standard thermodynamic and centroid-virial estimators, so an Eco calculation re-uses every existing Trotter estimator and requires only a few lines of code change once the frequencies are pre-computed. Numerical tests on hexagonal ice (kinetic/potential energies, RDFs, TePIGS VDOS) and on MOF-5 with a neuroevolution potential (negative thermal expansion, constant-pressure heat capacity) show monotonic convergence that is comparable to projected-Hessian Suzuki-Chin accuracy at pure second-order force cost.","tokens_in":18262,"tokens_out":836,"duration_ms":21065,"significance":"If the reported gains hold more generally, Eco supplies a practical, low-overhead upgrade to the workhorse Trotter discretisation that can become the default for condensed-phase PIMD. The algebraic argument that any zero-centroid circulant K leaves the usual estimators intact is clean and immediately useful; the SI optimisation routine and the GPUMD patch make the method reproducible and immediately deployable. Fair accounting for the extra force evaluations of projected-Hessian SC, together with successful application to a second-derivative property (Cp of MOF-5) that is notoriously hard to converge, strengthens the practical claim. The work therefore sits at the useful intersection of methodological simplicity and demonstrated numerical impact.","major_comments":[],"minor_comments":[{"comment":"Section II.A, first sentence: typographical error “Is is straightforward” should be “It is straightforward”.","section":"II.A"},{"comment":"Figure 5 caption and surrounding text: the comparison with Kapil et al. is valuable, but a brief explicit statement of the machine-learned potential they used (versus q-TIP4P/F here) would help readers judge how close the Trotter baselines really are.","section":"III.A"},{"comment":"Section II.B and Fig. 2: a short remark on how sensitive the optimised frequencies (and subsequent anharmonic observables) are to modest variations in ω_max would be helpful for users who must choose this parameter for a new system.","section":"II.B"},{"comment":"Appendix A: the adiabatic masses m_k = (β_e ω_k,e / β Ω)^{2} m are stated clearly, but a one-line reminder that the same Eco frequencies evaluated at T_e are used would remove any ambiguity for implementers.","section":"Appendix A"},{"comment":"References: the recent Hunt–Althorpe Matsubara-tail paper is properly cited as the conceptual inspiration; a parenthetical note that the present least-squares problem is formally distinct (fitting free ring-polymer frequencies rather than a rational Matsubara tail) would further clarify the relationship.","section":"I"}],"recommendation":"accept","confidential_remarks":"Strong, clean methods paper with open code and a non-trivial materials application. Fits well for a computational-chemistry or chemical-physics journal; I see no novelty or citation concerns."},"author_rebuttal":null,"desk_editor":{"model":"grok-4.5","letter":"This is a practical win for anyone who still runs Trotter PIMD. The free parameters of a zero-centroid symmetric circulant spring matrix have always been there; Zeng and Manolopoulos simply optimise the independent normal-mode frequencies by least-squares matching of the exact harmonic radius of gyration up to a system-dependent ω_max. Once those frequencies are precomputed (seconds of Newton iteration), the rest of the calculation is ordinary Trotter: same estimators, same code, a few lines changed.\n\nWhat they do well is show that the algebraic structure that validates the thermodynamic and centroid-virial estimators survives any such K (Sec. II.C). The ice numbers are clean: Eco energies and RDFs reach the accuracy of a 128-bead Trotter run with roughly half the beads, and the log-log slopes look fourth-order-like. Fair accounting against Kapil et al.’s projected-Hessian SC (extra force evaluations) puts Eco ahead per force call. The MOF-5 GPUMD implementation with a neuroevolution potential is the more impressive test: negative thermal expansion and the low-T heat capacity both converge far faster than Trotter and line up with experiment at P=256 where Trotter is still badly wrong.\n\nThe soft spot is exactly the one the reader flags: the harmonic fit is only empirically transferred to anharmonic condensed-phase observables. It works on ice and MOF-5, but there is no theorem that guarantees it for every potential. ω_max is a free parameter, though a transparent and cheap one. That is a minor caveat, not a load-bearing flaw; the paper is honest about it and the authors themselves flag reaction rates and instantons as the next place to look.\n\nMath is elementary and correct, citations are fair, code and SI are supplied. This is for people who actually run path-integral simulations of materials or water. I would bring it to reading group, I would cite it, and a serious editor should send it to referees without hesitation.","headline":"Clean, immediately usable PIMD acceleration that matches projected-Hessian SC accuracy at pure Trotter cost and reuses every standard estimator.","tokens_in":18832,"tokens_out":517,"would_cite":true,"duration_ms":5868,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"grok-4.5","headline":"Optimising free ring-polymer normal-mode frequencies to exact harmonic radii of gyration produces path integrals that converge like fourth-order methods while keeping second-order Trotter effort and estimators.","keywords":["path integral molecular dynamics","ring polymer","normal mode frequencies","nuclear quantum effects","Trotter discretisation","Suzuki-Chin","MOF-5","hexagonal ice"],"falsifier":"For the same potential, temperature and ω_max, run Eco and Trotter PIMD on hexagonal ice or MOF-5 and check whether Eco kinetic and potential energies, radial distribution peaks or heat capacity still approach the large-P limit faster (roughly fourth-order versus second-order) and with fewer force evaluations than Trotter or projected-Hessian Suzuki-Chin.","tokens_in":18962,"feed_emoji":"⚛️","tokens_out":1063,"duration_ms":15058,"temperature":0.7,"pith_summary":"Standard Trotter path integrals discretise imaginary-time quantum statistics with a fixed spring matrix whose normal-mode frequencies are not free parameters. This paper shows that the remaining free frequencies can be chosen by a simple least-squares fit so that the ring polymer reproduces the exact quantum radius of gyration of every harmonic oscillator up to a chosen maximum frequency. The resulting economised (Eco) path integral uses the same bead-to-mode transform, the same estimators, and only a few line changes in a Trotter code. On hexagonal ice the method reaches the accuracy of projected-Hessian fourth-order Suzuki-Chin calculations at purely second-order cost; on MOF-5 it converges the negative thermal expansion coefficient and the low-temperature heat capacity with far fewer beads than Trotter. A reader who runs path-integral molecular dynamics cares because nuclear quantum effects become affordable for large condensed-phase systems without new estimators or expensive Hessian projections.","feed_headline":"Fewer beads, same accuracy: path integrals get economised","feed_subtitle":"Fit free ring-polymer frequencies to harmonic radii of gyration and fourth-order convergence appears at second-order cost.","key_machinery":"Eco normal-mode frequencies: the floor(P/2) free frequencies obtained by minimising the root-mean-square fractional error between the path-integral and exact harmonic radii of gyration (or equivalently the quantum energy shifts) over 0 ≤ ω ≤ ω_max. Once precomputed they replace the Trotter frequencies inside an otherwise unchanged path-integral molecular-dynamics calculation.","core_discovery":"The Hessian of the free ring-polymer spring potential is a symmetric circulant matrix with a zero centroid eigenvalue; after symmetry constraints, floor(P/2) independent non-zero normal-mode frequencies remain free. Fitting those frequencies by least squares to the exact quantum radii of gyration of harmonic oscillators in the physical frequency window 0 ≤ ω ≤ ω_max produces an economised path integral whose convergence for anharmonic condensed-phase observables is comparable to fourth-order Suzuki-Chin methods while requiring only ordinary Trotter force evaluations and re-using every standard estimator.","pith_inferences":["Because the Eco spring matrix still commutes with the Trotter matrix and preserves cyclic and reversal symmetry, the same optimisation could be inserted into ring-polymer rate theories and instanton calculations that currently rely on Trotter springs.","The appearance of a plateau of nearly degenerate mid-range modes suggests that an even smaller effective number of independent frequencies might be sufficient, opening a route to further compression of the path integral.","If the harmonic fit is performed once at a reference temperature and density, transferability across modest temperature windows would make Eco especially convenient for constant-pressure heat-capacity scans."],"forward_implications":["Existing Trotter PIMD codes can be converted to Eco by precomputing a short table of frequencies and swapping them into the free ring-polymer propagator; no new estimators are required.","Properties that demand hundreds of Trotter beads (heat capacities, negative thermal expansion of frameworks) become accessible at P values comparable to fourth-order methods without projected Hessians.","The same Eco frequencies can be used inside elevated-temperature path-integral ground-state dynamics or other centroid-based spectral methods without altering the dynamical algorithm.","Once ω_max is chosen from a preliminary phonon calculation the method is parameter-free for the remainder of a simulation campaign."],"fun_headline_variants":["Eco path integrals match fourth-order accuracy at Trotter cost","Optimised bead frequencies yield faster path-integral convergence","Fit free ring-polymer modes to harmonic radii of gyration","Economised normal modes reach Suzuki-Chin accuracy cheaply","Least-squares Eco frequencies improve Trotter path integrals"],"cache_read_input_tokens":128,"weakest_assumption_plain":"That a least-squares fit of free normal-mode frequencies to the exact harmonic radius of gyration up to the highest physical frequency of the system will systematically improve convergence for realistic anharmonic condensed-phase observables.","fun_headline_variants_meta":{"raw":{"variants":["Eco path integrals match fourth-order accuracy at Trotter cost","Optimised bead frequencies yield faster path-integral convergence","Fit free ring-polymer modes to harmonic radii of gyration","Economised normal modes reach Suzuki-Chin accuracy cheaply","Least-squares Eco frequencies improve Trotter path integrals"]},"model":"grok-4.5","effort":"low","cost_usd":0.00532,"raw_usage":{"total_tokens":1524,"prompt_tokens":908,"num_sources_used":0,"completion_tokens":66,"cost_in_usd_ticks":53200000,"prompt_tokens_details":{"text_tokens":908,"audio_tokens":0,"image_tokens":0,"cached_tokens":128},"completion_tokens_details":{"audio_tokens":0,"reasoning_tokens":550,"accepted_prediction_tokens":0,"rejected_prediction_tokens":0}},"tokens_in":908,"tokens_out":66,"duration_ms":6097,"temperature":1.0,"reasoning_tokens":550,"cache_read_input_tokens":128,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-07-11T00:33:54.899149+00:00","model_set":{"reader":"grok-4.5"},"falsifier":"For the same potential, temperature and ω_max, run Eco and Trotter PIMD on hexagonal ice or MOF-5 and check whether Eco kinetic and potential energies, radial distribution peaks or heat capacity still approach the large-P limit faster (roughly fourth-order versus second-order) and with fewer force evaluations than Trotter or projected-Hessian Suzuki-Chin.","supporting_citations":[],"review_version":2}