{"status":"ok","message-type":"work","message-version":"1.0.0","message":{"indexed":{"date-parts":[[2026,10,1]],"date-time":"2026-10-01T10:59:46Z","timestamp":1790852386533,"version":"4.1.0"},"reference-count":38,"publisher":"AIP Publishing","issue":"13","license":[{"start":{"date-parts":[[2026,10,1]],"date-time":"2026-10-01T00:00:00Z","timestamp":1790812800000},"content-version":"vor","delay-in-days":0,"URL":"https:\/\/creativecommons.org\/licenses\/by\/4.0\/"}],"funder":[{"DOI":"10.13039\/501100000288","name":"Royal Society","doi-asserted-by":"publisher","award":["NIF\/R1\/254124"],"award-info":[{"award-number":["NIF\/R1\/254124"]}],"id":[{"id":"10.13039\/501100000288","id-type":"DOI","asserted-by":"publisher"}]}],"content-domain":{"domain":["pubs.aip.org"],"crossmark-restriction":true},"short-container-title":[],"published-print":{"date-parts":[[2026,10,7]]},"abstract":"<jats:p>The Hessian of the ring polymer spring potential in the standard Trotter path integral is a P \u00d7 P symmetric circulant matrix with a centroid eigenvalue of zero. All such matrices commute and are diagonalized by the same bead-to-normal mode transformation matrix, and their eigenvalues contain \u2308P\/2\u2309 \u2212 1 degenerate pairs by symmetry. However, this still leaves some freedom to improve on the Trotter approximation: one can optimize the remaining \u230aP\/2\u230b independent non-zero normal mode frequencies to fit the exact quantum mechanical radii of gyration of harmonic ring polymers with frequencies in the range 0 \u2264 \u03c9 \u2264 \u03c9max, where \u03c9max is the maximum physical frequency in the problem of interest. The optimization involves solving a simple least squares problem for the optimum (economized or \u201cEco\u201d) internal mode frequencies. The remainder of the calculation then proceeds in the same way as a Trotter path integral calculation. An example application to hexagonal ice shows that the convergence of the Eco path integral is comparable to that of the fourth order Suzuki\u2013Chin path integral, but with purely second order Trotter effort. There is no need to calculate the projected Hessians that arise in the Suzuki\u2013Chin method by finite differences; there is no need to develop any new estimators for observables, and once the Eco frequencies have been calculated, the implementation of the Eco path integral involves changing just a few lines of a Trotter path integral code. To provide a more impressive example, we have implemented the Eco method in GPUMD and used it to converge the (negative) thermal expansion coefficient and the constant pressure heat capacity of MOF-5 with a machine-learned neuroevolution potential.<\/jats:p>","DOI":"10.1063\/5.0351653","type":"journal-article","created":{"date-parts":[[2026,10,1]],"date-time":"2026-10-01T10:48:13Z","timestamp":1790851693000},"update-policy":"https:\/\/doi.org\/10.1063\/aip-crossmark-policy-page","source":"Crossref","is-referenced-by-count":0,"title":["Economized path integrals"],"prefix":"10.1063","volume":"165","author":[{"given":"Zezhu","family":"Zeng","sequence":"first","affiliation":[{"name":"Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford , South Parks Road, Oxford OX1 3QZ,","place":["United Kingdom"]}],"role":[{"vocabulary":"crossref","role":"author"}]},{"ORCID":"https:\/\/orcid.org\/0000-0002-7111-0763","authenticated-orcid":false,"given":"David E.","family":"Manolopoulos","sequence":"additional","affiliation":[{"name":"Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford , South Parks Road, Oxford OX1 3QZ,","place":["United Kingdom"]}],"role":[{"vocabulary":"crossref","role":"author"}]}],"member":"317","published-online":{"date-parts":[[2026,10,1]]},"reference":[{"key":"2026100106480561900_c1","doi-asserted-by":"publisher","first-page":"146401","DOI":"10.1103\/physrevlett.98.146401","article-title":"Generalized neural-network representation of high-dimensional potential-energy surfaces","volume":"98","year":"2007","journal-title":"Phys. Rev. Lett."},{"key":"2026100106480561900_c2","doi-asserted-by":"publisher","first-page":"136403","DOI":"10.1103\/physrevlett.104.136403","article-title":"Gaussian approximation potentials: The accuracy of quantum mechanics, without the electrons","volume":"104","year":"2010","journal-title":"Phys. Rev. Lett."},{"key":"2026100106480561900_c3","doi-asserted-by":"publisher","first-page":"860","DOI":"10.1063\/1.446740","article-title":"Study of an F center in molten KCl","volume":"80","year":"1984","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c4","volume-title":"Quantum Mechanics and Path Integrals","year":"1965"},{"key":"2026100106480561900_c5","doi-asserted-by":"publisher","first-page":"4078","DOI":"10.1063\/1.441588","article-title":"Exploiting the isomorphism between quantum theory and classical statistical mechanics of polyatomic fluids","volume":"74","year":"1981","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c6","doi-asserted-by":"publisher","first-page":"0109","DOI":"10.1038\/s41570-017-0109","article-title":"Nuclear quantum effects enter the mainstream","volume":"2","year":"2018","journal-title":"Nat. Rev. Chem."},{"key":"2026100106480561900_c7","doi-asserted-by":"publisher","first-page":"3765","DOI":"10.1143\/jpsj.53.3765","article-title":"Monte Carlo calculation of quantum systems. II. Higher order correction","volume":"53","year":"1984","journal-title":"J.\u00a0Phys. Soc. Jpn."},{"key":"2026100106480561900_c8","doi-asserted-by":"publisher","first-page":"425","DOI":"10.1016\/0375-9601(95)00266-6","article-title":"Hybrid exponential product formulas for unbounded operators with possible applications to Monte Carlo simulations","volume":"201","year":"1995","journal-title":"Phys. Lett. A"},{"key":"2026100106480561900_c9","doi-asserted-by":"publisher","first-page":"344","DOI":"10.1016\/s0375-9601(97)00003-0","article-title":"Symplectic integrators from composite operator factorizations","volume":"226","year":"1997","journal-title":"Phys. Lett. A"},{"key":"2026100106480561900_c10","doi-asserted-by":"publisher","first-page":"234103","DOI":"10.1063\/1.4971438","article-title":"High-order path integrals made easy","volume":"145","year":"2016","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c11","doi-asserted-by":"publisher","first-page":"084104","DOI":"10.1063\/1.3556661","article-title":"Accelerating the convergence of path integral dynamics with a generalized Langevin equation","volume":"134","year":"2011","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c12","doi-asserted-by":"publisher","first-page":"100604","DOI":"10.1103\/physrevlett.109.100604","article-title":"Efficient first-principles calculation of the quantum kinetic energy and momentum distribution of nuclei","volume":"109","year":"2012","journal-title":"Phys. Rev. Lett."},{"key":"2026100106480561900_c13","doi-asserted-by":"publisher","first-page":"1368","DOI":"10.1039\/c5sc03443d","article-title":"Modeling quantum nuclei with perturbed path integral molecular dynamics","volume":"7","year":"2016","journal-title":"Chem. Sci."},{"key":"2026100106480561900_c14","doi-asserted-by":"publisher","first-page":"084113","DOI":"10.1063\/5.0314385","article-title":"Exploiting the path-integral radius of gyration in open quantum dynamics","volume":"164","year":"2026","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c15","volume-title":"Circulant Matrices","year":"1994","edition":"2"},{"key":"2026100106480561900_c16","doi-asserted-by":"publisher","first-page":"124104","DOI":"10.1063\/1.3489925","article-title":"Efficient stochastic thermostatting of path integral molecular dynamics","volume":"133","year":"2010","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c17","doi-asserted-by":"publisher","first-page":"5150","DOI":"10.1063\/1.442815","article-title":"On path integral Monte Carlo simulations","volume":"76","year":"1982","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c18","doi-asserted-by":"publisher","first-page":"024501","DOI":"10.1063\/1.3167790","article-title":"Competing quantum effects in the dynamics of a flexible water model","volume":"131","year":"2009","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c19","doi-asserted-by":"publisher","first-page":"1518","DOI":"10.1063\/1.473300","article-title":"Unit cells for the simulation of hexagonal ice","volume":"106","year":"1997","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c20","doi-asserted-by":"publisher","first-page":"014101","DOI":"10.1063\/1.2408420","article-title":"Canonical sampling through velocity rescaling","volume":"126","year":"2007","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c22","doi-asserted-by":"publisher","first-page":"234116","DOI":"10.1063\/1.4883861","article-title":"How to remove the spurious resonances from ring polymer molecular dynamics","volume":"140","year":"2014","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c23","doi-asserted-by":"publisher","first-page":"054109","DOI":"10.1063\/1.5100587","article-title":"Path-integral dynamics of water using curvilinear centroids","volume":"151","year":"2019","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c24","doi-asserted-by":"publisher","first-page":"231101","DOI":"10.1063\/5.0076704","article-title":"Fast quasi-centroid molecular dynamics","volume":"155","year":"2021","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c25","doi-asserted-by":"publisher","first-page":"181102","DOI":"10.1063\/5.0120386","article-title":"Quantum dynamics using path integral coarse-graining","volume":"157","year":"2022","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c26","doi-asserted-by":"publisher","first-page":"9172","DOI":"10.1021\/acs.jpcb.3c05028","article-title":"Fast quasi-centroid molecular dynamics for water and ice","volume":"127","year":"2023","journal-title":"J.\u00a0Phys. Chem. B"},{"key":"2026100106480561900_c27","doi-asserted-by":"publisher","first-page":"2028","DOI":"10.1063\/1.470959","article-title":"Adiabatic path integral molecular dynamics methods. II. Algorithms","volume":"104","year":"1996","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c28","doi-asserted-by":"publisher","first-page":"031101","DOI":"10.1063\/1.3290958","article-title":"Communications: On artificial frequency shifts in infrared spectra obtained from centroid molecular dynamics: Quantum liquid water","volume":"132","year":"2010","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c29","doi-asserted-by":"publisher","first-page":"114801","DOI":"10.1063\/5.0106617","article-title":"GPUMD: A package for constructing accurate machine-learned potentials and performing highly efficient atomistic simulations","volume":"157","year":"2022","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c30","doi-asserted-by":"publisher","first-page":"064109","DOI":"10.1063\/5.0241006","article-title":"Highly efficient path-integral molecular dynamics simulations with GPUMD using neuroevolution potentials: Case studies on thermal properties of materials","volume":"162","year":"2025","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c31","doi-asserted-by":"publisher","first-page":"e2415664122","DOI":"10.1073\/pnas.2415664122","article-title":"Lattice distortion leads to glassy thermal transport in crystalline Cs3Bi2I6Cl3","volume":"122","year":"2025","journal-title":"Proc. Natl. Acad. Sci. U. S. A."},{"key":"2026100106480561900_c32","doi-asserted-by":"publisher","first-page":"104309","DOI":"10.1103\/physrevb.104.104309","article-title":"Neuroevolution machine learning potentials: Combining high accuracy and low cost in atomistic simulations and application to heat transport","volume":"104","year":"2021","journal-title":"Phys. Rev. B"},{"key":"2026100106480561900_c33","unstructured":"Z.\n            Zeng\n           and D. E.Manolopoulos, A patch to GPUMD for economised path-integrals (2026); https:\/\/github.com\/ZengZezhu\/Eco-Path-Integrals-in-GPUMD."},{"key":"2026100106480561900_c34","doi-asserted-by":"publisher","first-page":"18","DOI":"10.1038\/s41524-024-01205-w","article-title":"Machine learned force-fields for an ab-initio quality description of metal-organic frameworks","volume":"10","year":"2024","journal-title":"npj Comput. Mater."},{"key":"2026100106480561900_c35","doi-asserted-by":"publisher","first-page":"16181","DOI":"10.1021\/jp103212z","article-title":"Elucidating negative thermal expansion in MOF-5","volume":"114","year":"2010","journal-title":"J.\u00a0Phys. Chem. C"},{"key":"2026100106480561900_c36","doi-asserted-by":"publisher","first-page":"1","DOI":"10.1016\/j.micromeso.2015.05.047","article-title":"Specific heat capacities of MOF-5, Cu-BTC, Fe-BTC, MOF-177 and MIL-53 (Al) over wide temperature ranges: Measurements and application of empirical group contribution method","volume":"217","year":"2015","journal-title":"Microporous Mesoporous Mater."},{"key":"2026100106480561900_c37","doi-asserted-by":"publisher","first-page":"2348","DOI":"10.1021\/acs.jctc.5c02073","article-title":"Graph-based internal coordinate analysis for transition state characterization","volume":"22","year":"2026","journal-title":"J.\u00a0Chem. Theory Comput."},{"key":"2026100106480561900_c38","doi-asserted-by":"publisher","first-page":"104101","DOI":"10.1063\/1.2013257","article-title":"Path-integral virial estimator based on the scaling of fluctuation coordinates: Application to quantum clusters","volume":"123","year":"2005","journal-title":"J.\u00a0Chem. Phys."},{"key":"2026100106480561900_c39","doi-asserted-by":"publisher","first-page":"204102","DOI":"10.1063\/5.0300048","article-title":"Vibrational spectra of materials and molecules from partially adiabatic elevated-temperature centroid molecular dynamics","volume":"163","year":"2025","journal-title":"J.\u00a0Chem. Phys."}],"container-title":["The Journal of Chemical Physics"],"original-title":[],"language":"en","link":[{"URL":"https:\/\/pubs.aip.org\/aip\/jcp\/article-pdf\/doi\/10.1063\/5.0351653\/21157157\/134104_1_5.0351653.pdf","content-type":"application\/pdf","content-version":"vor","intended-application":"syndication"},{"URL":"https:\/\/pubs.aip.org\/aip\/jcp\/article-pdf\/doi\/10.1063\/5.0351653\/21157157\/134104_1_5.0351653.pdf","content-type":"unspecified","content-version":"vor","intended-application":"similarity-checking"}],"deposited":{"date-parts":[[2026,10,1]],"date-time":"2026-10-01T10:48:15Z","timestamp":1790851695000},"score":1,"resource":{"primary":{"URL":"https:\/\/pubs.aip.org\/jcp\/article\/165\/13\/134104\/3405671\/Economized-path-integrals"}},"subtitle":[],"short-title":[],"issued":{"date-parts":[[2026,10,1]]},"references-count":38,"journal-issue":{"issue":"13","published-print":{"date-parts":[[2026,10,7]]}},"URL":"https:\/\/doi.org\/10.1063\/5.0351653","relation":{},"ISSN":["0021-9606","1089-7690"],"issn-type":[{"value":"0021-9606","type":"print"},{"value":"1089-7690","type":"electronic"}],"subject":[],"published-other":{"date-parts":[[2026,10,7]]},"published":{"date-parts":[[2026,10,1]]},"article-number":"134104"}}