pith. sign in

arxiv: 2605.19858 · v1 · pith:IEWT2CJTnew · submitted 2026-05-19 · ⚛️ physics.chem-ph · cond-mat.mtrl-sci· cond-mat.str-el· physics.atm-clus· physics.comp-ph

Accelerated "on-the-fly" coupled-cluster path-integral molecular dynamics: Impact of nuclear quantum effects on an asymmetric proton

Pith reviewed 2026-05-20 02:00 UTC · model grok-4.3

classification ⚛️ physics.chem-ph cond-mat.mtrl-scicond-mat.str-elphysics.atm-clusphysics.comp-ph
keywords path-integral molecular dynamicscoupled-cluster theorynuclear quantum effectsproton transferasymmetric hydrogen bondring-polymer contractionnuclear magnetic shielding
0
0 comments X

The pith

An accelerated coupled-cluster path-integral method shows nuclear quantum effects must be treated together with electron correlation for accurate asymmetric proton positions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops an efficient on-the-fly coupled-cluster path-integral molecular dynamics approach by combining three acceleration techniques. It applies the method to a proton shared between water and formaldehyde and finds that nuclear quantum effects broaden bond-length distributions, move the average proton position closer to the midpoint, and lower the probability it sits nearer to formaldehyde from 81.7 percent to 61.1 percent. The same simulations show that electron correlation and nuclear quantum effects produce shielding changes of similar size that can point in opposite directions. These outcomes indicate that neither effect can be omitted if the goal is to predict the behavior of asymmetric hydrogen bonds reliably.

Core claim

The combination of quantum ring-polymer contraction, second-generation Car-Parrinello-like dynamics of the Hartree-Fock reference, and basis-consistent extrapolation of coupled-cluster amplitudes makes correlated PIMD calculations feasible. Applied to the proton shared by water and formaldehyde, the resulting simulations show that relative to classical nuclei the nuclear quantum effects broaden covalent X-H bond-length distributions, reduce the bias of the shared proton toward formaldehyde, and shift the mean proton-transfer coordinate from 0.206 to 0.135 Å while lowering the probability of finding the proton closer to formaldehyde from 81.7% to 61.1%. The corresponding nuclear magnetic-shi­

What carries the argument

Quantum ring-polymer contraction (qRPC) decomposition, which evaluates the inexpensive Hartree-Fock potential on the full ring polymer while restricting the expensive coupled-cluster correction to the centroid only, augmented by accelerated dynamics and amplitude extrapolation.

If this is right

  • Correlated PIMD becomes computationally practical for finite-temperature studies of proton-sharing systems.
  • Predictive modeling of asymmetric hydrogen bonds requires simultaneous inclusion of correlated electronic structure and nuclear quantum fluctuations.
  • Nuclear quantum effects measurably reduce the classical bias of the shared proton toward one partner molecule.
  • Nuclear magnetic shielding tensors receive contributions from correlation and nuclear quantum motion that are comparable in size and can act in opposing directions.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same acceleration strategy could be transferred to other hydrogen-bonded complexes where classical simulations currently overestimate proton localization.
  • Classical-nucleus models may systematically overstate the stability of certain proton configurations in asymmetric bonds.
  • Direct comparison of the simulated shielding tensors with high-resolution gas-phase NMR data would test whether the opposing-effect observation holds in experiment.

Load-bearing premise

Evaluating the coupled-cluster correction only at the average position of the ring-polymer beads still captures the proton-transfer coordinate and shielding tensors accurately enough that the added acceleration layers do not introduce large uncontrolled errors.

What would settle it

An experimental measurement of the proton-transfer coordinate histogram or the average nuclear magnetic shielding tensor for the water-formaldehyde complex that lies well outside the range produced by the quantum-nuclear simulation.

Figures

Figures reproduced from arXiv: 2605.19858 by Hossam Elgabarty, Thomas D. K\"uhne, Thomas Spura.

Figure 1
Figure 1. Figure 1: Mean number of iterations required to converge the cluster and [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Energy and force deviations along a CC-MD trajectory propagated with two iterations [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Partial C–H PCF obtained from the CC-MD and CC-PIMD simulations. The first [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Free-energy distribution, in kcal mol−1 , of the shared proton in our CC-MD and CC-PIMD simulations as a function of the intermolecular O–O distance and the proton reaction coordinate ν. 8 [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Partial O–H+ PCF obtained from the CC-MD and CC-PIMD simulations, decomposed into the contributions from the formaldehyde oxygen OCH2O and the water oxygen OH2O to the shared proton H+. 4.4 Nuclear magnetic resonance (NMR) response of the asymmetric hydrogen bond Nuclear magnetic shielding tensors are highly sensitive to small changes in nuclear geometry and electronic structure. A bond-length change of on… view at source ↗
Figure 6
Figure 6. Figure 6: Distribution of the isotropic nuclear magnetic shielding [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Distribution of the isotropic shielding of the hydrogens bound to water and carbon. [PITH_FULL_IMAGE:figures/full_fig_p011_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Decomposition of the NMR response of the shared proton and the formaldehyde C=O [PITH_FULL_IMAGE:figures/full_fig_p011_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Impact of nuclear quantum effects on the isotropic nuclear magnetic shielding [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Impact of electron correlation on the isotropic nuclear magnetic shielding [PITH_FULL_IMAGE:figures/full_fig_p014_10.png] view at source ↗
read the original abstract

We present an accelerated ``on-the-fly'' coupled-cluster path-integral molecular dynamics (PIMD) method for finite-temperature simulations in which electron correlation and nuclear quantum effects are treated simultaneously. The approach is based on our quantum ring-polymer contraction (qRPC) technique, in which the inexpensive Hartree-Fock potential is evaluated on the full ring-polymer, while the expensive coupled-cluster correction is evaluated on the centroid only. This qRPC decomposition is combined with a second-generation Car-Parrinello-like dynamics of the Hartree-Fock reference and a basis-consistent extrapolation of the coupled-cluster and de-excitation amplitudes. The combination of all three acceleration layers is essential for making correlated PIMD calculations feasible. We apply the method to a proton shared by water and formaldehyde. Relative to classical nuclei, nuclear quantum effects broaden covalent X--H bond-length distributions, reduce the bias of the shared proton toward formaldehyde, and shift the mean proton-transfer coordinate from 0.206 to 0.135A. The probability of finding the proton closer to formaldehyde decreases from 81.7$\%$ to 61.1$\%$. The corresponding nuclear magnetic shielding tensors show that electron correlation and nuclear quantum effects are of comparable magnitude and can act in opposite directions. These results demonstrate that predictive simulations of asymmetric hydrogen bonds require a simultaneous treatment of correlated electronic structure and nuclear quantum fluctuations.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript presents an accelerated on-the-fly coupled-cluster path-integral molecular dynamics (PIMD) method based on quantum ring-polymer contraction (qRPC), in which the Hartree-Fock potential is evaluated on the full ring-polymer while the coupled-cluster correction is evaluated only on the centroid. This is combined with second-generation Car-Parrinello-like dynamics of the Hartree-Fock reference and basis-consistent extrapolation of coupled-cluster amplitudes. Applied to an asymmetric proton shared between water and formaldehyde, the simulations report that nuclear quantum effects broaden X-H bond distributions, shift the mean proton-transfer coordinate from 0.206 Å to 0.135 Å, and reduce the probability of the proton being closer to formaldehyde from 81.7% to 61.1%. Electron correlation and nuclear quantum effects are found to be of comparable magnitude and can act in opposite directions on the nuclear magnetic shielding tensors.

Significance. If the qRPC approximation and acceleration layers preserve the required accuracy, the work is significant for enabling feasible correlated PIMD on systems with asymmetric hydrogen bonds. It provides concrete numerical evidence that both electron correlation and nuclear quantum fluctuations must be treated simultaneously for predictive modeling, with opposing effects on observables such as shielding tensors. The combination of qRPC, Car-Parrinello-like HF dynamics, and extrapolation is a practical strength that addresses computational bottlenecks in the field.

major comments (2)
  1. [Results and Methods sections] The central claim that simultaneous treatment of correlation and NQE is required rests on the accuracy of the qRPC decomposition for the proton-transfer coordinate and shielding tensors. No direct benchmark comparing qRPC-CC-PIMD statistics to a reference full-CC ring-polymer calculation (on the target system or a smaller proxy) is provided, leaving open the possibility of systematic under- or over-correction in the centroid-only CC term.
  2. [Results] The reported shifts (mean coordinate 0.206 Å to 0.135 Å; probability 81.7% to 61.1%) and the statement of comparable-magnitude opposing effects on shielding lack accompanying error bars, convergence tests with respect to ring-polymer size, or statistical uncertainties. These are load-bearing for the quantitative conclusions about NQE impact.
minor comments (2)
  1. [Abstract] In the abstract, the unit '0.135A' should include a space and proper angstrom symbol for consistency with standard notation.
  2. [Abstract] The LaTeX artifact '%$' in the abstract should be corrected to a properly rendered percentage symbol.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their careful reading of the manuscript and for the constructive comments. We address each major comment below, indicating the revisions we will make where appropriate.

read point-by-point responses
  1. Referee: [Results and Methods sections] The central claim that simultaneous treatment of correlation and NQE is required rests on the accuracy of the qRPC decomposition for the proton-transfer coordinate and shielding tensors. No direct benchmark comparing qRPC-CC-PIMD statistics to a reference full-CC ring-polymer calculation (on the target system or a smaller proxy) is provided, leaving open the possibility of systematic under- or over-correction in the centroid-only CC term.

    Authors: We acknowledge that a direct benchmark against full coupled-cluster ring-polymer calculations would provide stronger validation of the qRPC approximation for the quantities of interest. Such calculations remain computationally prohibitive for the target system and even for smaller proxies at the CCSD(T) level, which motivated the development of the qRPC approach in the first place. In earlier work introducing qRPC, we performed exact comparisons on smaller model systems demonstrating that the centroid-only correlation correction introduces only minor errors for structural observables. We will expand the Methods section with a more detailed discussion of these prior validations, the theoretical basis for applying the correction at the centroid, and an estimate of the expected residual error for the proton-transfer coordinate and shielding tensors in the present application. revision: partial

  2. Referee: [Results] The reported shifts (mean coordinate 0.206 Å to 0.135 Å; probability 81.7% to 61.1%) and the statement of comparable-magnitude opposing effects on shielding lack accompanying error bars, convergence tests with respect to ring-polymer size, or statistical uncertainties. These are load-bearing for the quantitative conclusions about NQE impact.

    Authors: We agree that the quantitative results require accompanying statistical uncertainties and convergence information. In the revised manuscript we will report error bars on the mean proton-transfer coordinate, the reported probabilities, and the shielding tensor components, obtained via block averaging of the production trajectories. We will also add explicit convergence tests with respect to ring-polymer size (e.g., results for 8, 16, and 32 beads) to confirm that the observed shifts are stable. These additions will be placed in the Results section. revision: yes

standing simulated objections not resolved
  • Direct benchmark of qRPC-CC-PIMD statistics against a full-CC ring-polymer reference calculation on the target system or a suitable smaller proxy, which remains computationally infeasible.

Circularity Check

0 steps flagged

No significant circularity; method accelerations and results are independent of target observables

full rationale

The paper's claimed chain introduces qRPC decomposition (HF on full ring-polymer, CC on centroid), second-generation Car-Parrinello-like HF dynamics, and basis-consistent extrapolation to enable correlated PIMD. These accelerations are technical and not defined circularly in terms of the proton-transfer coordinate, probability shifts, or shielding tensors. Reported results (mean coordinate 0.206 Å → 0.135 Å; probability 81.7% → 61.1%) are direct simulation outputs. Self-citation to prior qRPC work supports the base technique but is not load-bearing for the new combination or the conclusion that simultaneous treatment is required; the derivation remains self-contained against external molecular-dynamics benchmarks with no reductions by construction, fitted inputs renamed as predictions, or imported uniqueness theorems.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central feasibility claim rests on the accuracy of the qRPC decomposition and the two additional acceleration layers; these are domain assumptions rather than derived quantities. No explicit free parameters or new invented entities are named in the abstract.

axioms (1)
  • domain assumption The Hartree-Fock potential evaluated on the full ring-polymer plus coupled-cluster correction on the centroid accurately captures the relevant potential energy surface for proton transfer.
    This is the core of the qRPC technique invoked to make the calculation tractable.

pith-pipeline@v0.9.0 · 5810 in / 1405 out tokens · 32005 ms · 2026-05-20T02:00:11.778660+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

76 extracted references · 76 canonical work pages

  1. [1]

    Feynman and Albert R

    Richard P. Feynman and Albert R. Hibbs.Quantum Mechanics and Path Integrals. McGraw- Hill, New York, 1965

  2. [2]

    Parrinello and A

    M. Parrinello and A. Rahman. Study of an F center in molten KCl.J. Chem. Phys., 80(2):860, 1984

  3. [3]

    David Chandler and Peter G. Wolynes. Exploiting the isomorphism between quantum theory and classical statistical mechanics of polyatomic fluids.J. Chem. Phys., 74(7):4078, 1981

  4. [4]

    Manolopoulos, Thomas E

    Scott Habershon, David E. Manolopoulos, Thomas E. Markland, and Thomas F. Miller. Ring-polymer molecular dynamics: Quantum effects in chemical dynamics from classical trajectories in an extended phase space.Annu. Rev. Phys. Chem., 64(1):387–413, 2013

  5. [5]

    Ab initio path-integral molecular dynamics.Z

    Dominik Marx and Michele Parrinello. Ab initio path-integral molecular dynamics.Z. Phys. B Cond. Mat., 95(2):143–144, 1994. 15

  6. [6]

    Tuckerman, Dominik Marx, Michael L

    Mark E. Tuckerman, Dominik Marx, Michael L. Klein, and Michele Parrinello. Efficient and general algorithms for path integral Car–Parrinello molecular dynamic.J. Chem. Phys., 104(14):5579–5588, 1996

  7. [7]

    Berard, and Gregory A

    Marc Pavese, Daniel R. Berard, and Gregory A. Voth. Ab initio centroid molecular dynamics: a fully quantum method for condensed-phase dynamics simulations.Chem. Phys. Lett., 300(1-2):93–98, 1999

  8. [8]

    A unified scheme for ab initio molecular orbital theory and path integral molecular dynamics.J

    Motoyuki Shiga, Masanori Tachikawa, and Shinichi Miura. A unified scheme for ab initio molecular orbital theory and path integral molecular dynamics.J. Chem. Phys., 115(20):9149, 2001

  9. [9]

    Tuckerman, Dominik Marx, Michael L

    Mark E. Tuckerman, Dominik Marx, Michael L. Klein, and Michele Parrinello. On the quantum nature of the shared proton in hydrogen bonds.Science, 275(5301):817–820, 1997

  10. [10]

    An ab initio path integral monte carlo simulation method for molecules and clusters: Application to Li4 and Li+ 5 .J

    Ruben O Weht, Jorge Kohanoff, Darıo A Estrin, and Charusita Chakravarty. An ab initio path integral monte carlo simulation method for molecules and clusters: Application to Li4 and Li+ 5 .J. Chem. Phys., 108(21):8848–8858, 1998

  11. [11]

    Tunnelling and zero-point motion in high-pressure ice.Nature, 392(6673):258–261, 1998

    Magali Benoit, Dominik Marx, and Michele Parrinello. Tunnelling and zero-point motion in high-pressure ice.Nature, 392(6673):258–261, 1998

  12. [12]

    Tuckerman, Jürg Hutter, and Michele Parrinello

    Dominik Marx, Mark E. Tuckerman, Jürg Hutter, and Michele Parrinello. The nature of the hydrated excess proton in water.Nature, 397(6720):601–604, 1999

  13. [13]

    Tuckerman, Dominik Marx, and Michele Parrinello

    Mark E. Tuckerman, Dominik Marx, and Michele Parrinello. The nature and transport mechanism of hydrated hydroxide ions in aqueous solution.Nature, 417(6892):925–929, 2002

  14. [14]

    Nuclear quantum effects in water.Phys

    Joseph Morrone and Roberto Car. Nuclear quantum effects in water.Phys. Rev. Lett., 101(1), 2008

  15. [15]

    On-the-Fly

    Anna Kaczmarek, Motoyuki Shiga, and Dominik Marx. Quantum Effects on Vibrational and Electronic Spectra of Hydrazine Studied by “On-the-Fly” ab Initio Ring Polymer Molecular Dynamics.J. Phys. Chem. A, 113(10):1985–1994, 2009

  16. [16]

    Xin-Zheng Li, Matthew I. J. Probert, Ali Alavi, and Angelos Michaelides. Quantum nature of the proton in water-hydroxyl overlayers on metal surfaces.Phys. Rev. Lett., 104(6), 2010

  17. [17]

    X.-Z. Li, B. Walker, and A. Michaelides. Quantum nature of the hydrogen bond.Proc. Nat. Acad. Sci. USA, 108(16):6369–6373, 2011

  18. [18]

    Ceriotti, J

    M. Ceriotti, J. Cuny, M. Parrinello, and D. E. Manolopoulos. Nuclear quantum effects and hydrogen bond fluctuations in water.Proc. Nat. Acad. Sci. USA, 110(39):15591–15596, 2013

  19. [19]

    Giberti, A

    F. Giberti, A. A. Hassanali, M. Ceriotti, and M. Parrinello. The role of quantum effects on structural and electronic fluctuations in neat and charged water.J. Phys. Chem. B, 118:13226–13235, 2014

  20. [20]

    Hassanali, and Thomas D

    Jan Kessler, Hossam Elgabarty, Thomas Spura, Kristof Karhan, Pouya Partovi-Azar, Ali A. Hassanali, and Thomas D. Kühne. Structure and dynamics of the instantaneous water/vapor interface revisited by path-integral and ab initio molecular dynamics simulations.J. Phys. Chem. B, 119(31):10079–10086, 2015

  21. [21]

    Thomas Spura, Christopher John, Scott Habershon, and Thomas D. Kühne. Nuclear quantum effects in liquid water from path-integral simulations using an ab initio force- matching approach.Mol. Phys., 113(8):808–822, 2015. 16

  22. [22]

    Azadi, R

    S. Azadi, R. Singh, and T. D. Kühne. Nuclear quantum effects induce metallization of dense solid molecular hydrogen.J. Comp. Chem., 39:262–268, 2018

  23. [23]

    Timothy Clark, Julian Joachim Heske, and Thomas D. Kühne. Opposing electronic and nuclear quantum effects on hydrogen bonds in H2O and D2O.ChemPhysChem, 20(19):2461– 2465, 2019

  24. [24]

    Kühne and Rustam Z

    Thomas D. Kühne and Rustam Z. Khaliullin. Electronic signature of the instantaneous asymmetry in the first coordination shell of liquid water.Nature Communications, 4:1450, 2013

  25. [25]

    Kühne and Rustam Z

    Thomas D. Kühne and Rustam Z. Khaliullin. Nature of the asymmetry in the hydrogen- bond networks of hexagonal ice and liquid water.J. Am. Chem. Soc., 136(9):3395–3399, 2014

  26. [26]

    Kühne, Matthias Krack, and Michele Parrinello

    Thomas D. Kühne, Matthias Krack, and Michele Parrinello. Static and Dynamical Properties of Liquid Water from First Principles by a Novel Car- Parrinello-like Approach.J. Chem. Theory Comput., 5(2):235–241, 2009

  27. [27]

    John A. Pople. Quantum chemical models (nobel lecture).Angewandte Chemie International Edition, 38(13-14):1894–1902, 1999

  28. [28]

    Wiley, Chichester, 2000

    Trygve Helgaker, Poul Jørgensen, and Jeppe Olsen.Molecular Electronic-Structure Theory. Wiley, Chichester, 2000

  29. [29]

    On-the-fly

    Thomas Spura, Hossam Elgabarty, and Thomas D. Kühne. “On-the-fly” coupled cluster path-integral molecular dynamics: Impact of nuclear quantum effects on the protonated water dimer.Phys. Chem. Chem. Phys., 17:14355 – 14359, 2015

  30. [30]

    Christopher John, Thomas Spura, Scott Habershon, and Thomas D. Kühne. Quantum ring-polymer contraction method: Including nuclear quantum effects at no additional computational cost in comparison to ab initio molecular dynamics.Phys. Rev. E, 93:043305, 2016

  31. [31]

    Markland and David E

    Thomas E. Markland and David E. Manolopoulos. A refined ring polymer contraction scheme for systems with electrostatic interactions.Chem. Phys. Lett., 464(4-6):256–261, 2008

  32. [32]

    Markland, and David E

    Scott Habershon, Thomas E. Markland, and David E. Manolopoulos. Competing quantum effects in the dynamics of a flexible water model.J. Chem. Phys., 131(2):024501, 2009

  33. [33]

    Kühne, Matthias Krack, Fawzi R

    Thomas D. Kühne, Matthias Krack, Fawzi R. Mohamed, and Michele Parrinello. Efficient and accurate car-parrinello-like approach to born-oppenheimer molecular dynamics.Phys. Rev. Lett., 98(6):066401, 2007

  34. [34]

    Thomas D. Kühne. Second generation car-parrinello molecular dynamics.WIREs: Comput. Mol. Sci., 4(4):391–406, 2014

  35. [35]

    On the correlation problem in atomic and molecular systems

    Jiří Čížek. On the correlation problem in atomic and molecular systems. calculation of wavefunction components in ursell-type expansion using quantum-field theoretical methods. J. Chem. Phys., 45:4256, 1966

  36. [36]

    Coupled-cluster theory

    Jürgen Gauss. Coupled-cluster theory. In P. von R. Schleyer, editor,Encyclopedia of Computational Chemistry, volume 1, pages 615–636. Wiley, Chichester, 1998

  37. [37]

    Bartlett and Monika Musiał

    Rodney J. Bartlett and Monika Musiał. Coupled-cluster theory in quantum chemistry.Rev. Mod. Phys., 79:291–352, 2007. 17

  38. [38]

    On the concepts of connectivity, separability, and consistency: An illustration by partitioned diagrams and numerical probing.Chem

    Michael Hanrath. On the concepts of connectivity, separability, and consistency: An illustration by partitioned diagrams and numerical probing.Chem. Phys., 356(1-3):31–38, 2009

  39. [39]

    Ab initio methods for the calculation of NMR shielding and indirect spin–spin coupling constants.Chem

    Trygve Helgaker, Michał Jaszuński, and Kenneth Ruud. Ab initio methods for the calculation of NMR shielding and indirect spin–spin coupling constants.Chem. Rev., 99(1):293–352, 1999

  40. [40]

    Bartlett and George D

    Rodney J. Bartlett and George D. Purvis. Many-body perturbation theory, coupled-pair many-electron theory, and the importance of quadruple excitations for the correlation problem.Int. J. Quantum Chem., 14(5):561–581, 1978

  41. [41]

    Adamowicz, W

    L. Adamowicz, W. D. Laidig, and R. J. Bartlett. Analytical gradients for the coupled-cluster method.Int. J. Quantum Chem., 26(S18):245–254, 1984

  42. [42]

    Coupled cluster response functions.J

    Henrik Koch and Poul Jørgensen. Coupled cluster response functions.J. Chem. Phys., 93(5):3333, 1990

  43. [43]

    Analytic energy derivatives for coupled-cluster methods describing excited states: General formulas and comparison of computational costs.Int

    PG Szalay. Analytic energy derivatives for coupled-cluster methods describing excited states: General formulas and comparison of computational costs.Int. J. Quantum Chem., 55(2):151–163, 1995

  44. [44]

    Variational principles and linked-cluster exp s expansions for static and dynamic many-body problems.Annals of Physics, 151(2):311–382, 1983

    Jouko Arponen. Variational principles and linked-cluster exp s expansions for static and dynamic many-body problems.Annals of Physics, 151(2):311–382, 1983

  45. [45]

    Mihály Kállay, Jürgen Gauss, and Péter G. Szalay. Analytic first derivatives for general coupled-cluster and configuration interaction models.J. Chem. Phys., 119(6):2991–3004, 2003

  46. [46]

    On the non-orthogonality problem connected with the use of atomic wave functions in the theory of molecules and crystals.J

    Per-Olov Löwdin. On the non-orthogonality problem connected with the use of atomic wave functions in the theory of molecules and crystals.J. Chem. Phys., 18(3):365, 1950

  47. [47]

    J. C. Slater and G. F. Koster. Simplified LCao method for the periodic potential problem. Phys. Rev., 94(6):1498–1524, 1954

  48. [48]

    Time-reversible always stable predictor-corrector method for molecular dynamics of polarizable molecules.J

    Jiří Kolafa. Time-reversible always stable predictor-corrector method for molecular dynamics of polarizable molecules.J. Comput. Chem., 25(3):335–342, 2004

  49. [49]

    Markland and David E

    Thomas E. Markland and David E. Manolopoulos. An efficient ring polymer contraction scheme for imaginary time path integral simulations.J. Chem. Phys., 129(2):024105, 2008

  50. [50]

    Manolopoulos

    Michele Ceriotti, Joshua More, and David E. Manolopoulos. i-pi: A python interface for ab initio path integral molecular dynamics simulations.Comp. Phys. Comm., 185(3):1019–1026, 2014

  51. [51]

    Matthews, Lan Cheng, Michael E

    Devin A. Matthews, Lan Cheng, Michael E. Harding, Filippo Lipparini, Stella Stopkowicz, Thomas-C. Jagau, Péter G. Szalay, Jürgen Gauss, and John F. Stanton. Coupled-cluster techniques for computational chemistry: The CFOUR program package.J. Chem. Phys., 152(21):214108, 2020

  52. [52]

    Scheiner, Gustavo E

    Andrew C. Scheiner, Gustavo E. Scuseria, Julia E. Rice, Timothy J. Lee, and Henry F. Schaefer. Analytic evaluation of energy gradients for the single and double excitation coupled cluster (CCSD) wave function: Theory and application.J. Chem. Phys., 87(9):5361, 1987

  53. [53]

    Stanton, Jürgen Gauss, John D

    John F. Stanton, Jürgen Gauss, John D. Watts, and Rodney J. Bartlett. A direct prod- uct decomposition approach for symmetry exploitation in many-body methods. i. energy calculations.J. Chem. Phys., 94(6):4334, 1991. 18

  54. [54]

    Stanton, and Rodney J

    Jürgen Gauss, John F. Stanton, and Rodney J. Bartlett. Coupled-cluster open-shell analytic gradients: Implementation of the direct product decomposition approach in energy gradient calculations.J. Chem. Phys., 95(4):2623, 1991

  55. [55]

    Stanton, John D Watts, and Rodney J Bartlett

    Jürgen Gauss, Walter J Lauderdale, John F. Stanton, John D Watts, and Rodney J Bartlett. Analytic energy gradients for open-shell coupled-cluster singles and doubles (CCSD) calculations using restricted open-shell hartree-fock (ROHF) reference functions. Chem. Phys. Lett., 182(3):207–215, 1991

  56. [56]

    Thom H. Dunning. Gaussian basis sets for use in correlated molecular calculations. i. the atoms boron through neon and hydrogen.J. Chem. Phys., 90(2):1007, 1989

  57. [57]

    Turney, Andrew C

    Justin M. Turney, Andrew C. Simmonett, Robert M. Parrish, Edward G. Hohenstein, Francesco A. Evangelista, Justin T. Fermann, Benjamin J. Mintz, Lori A. Burns, Jeremiah J. Wilke, Micah L. Abrams, and et al. Psi4: an open-source ab initio electronic structure program.WIREs: Comput. Mol. Sci., 2(4):556–565, 2011

  58. [58]

    Daniel Crawford and Henry F

    T. Daniel Crawford and Henry F. Schaefer. An introduction to coupled cluster theory for computational chemists.Rev. Comp. Chem., 14:33–136, 2000

  59. [59]

    Jürgen Gauss and John F. Stanton. Perturbative treatment of triple excitations in coupled- cluster calculations of nuclear magnetic shielding constants.J. Chem. Phys., 104(7):2574, 1996

  60. [60]

    Jürgen Gauss and John F. Stanton. Electron-correlated methods for the calculation of NMR chemical shifts. In M. Kaupp, M. Bühl, and V. G. Malkin, editors,Calculation of NMR and EPR Parameters: Theory and Applications, pages 123–139. Wiley-VCH, Weinheim, 2004

  61. [61]

    Harding, Thorsten Metzroth, Jürgen Gauss, and Alexander A

    Michael E. Harding, Thorsten Metzroth, Jürgen Gauss, and Alexander A. Auer. Parallel calculation of CCSD and CCSD(T) analytic first and second derivatives.J. Chem. Theory Comput., 4(1):64–74, 2008

  62. [62]

    Convergence acceleration of iterative sequences

    Péter Pulay. Convergence acceleration of iterative sequences. The case of SCF iteration. Chem. Phys. Lett., 73(2):393–398, 1980

  63. [63]

    P. Pulay. Improved SCF convergence acceleration.J. Comput. Chem., 3(4):556–560, 1982

  64. [64]

    Shielding tensor data and structure: the bond-related chemical shift concept

    AR Grimmer. Shielding tensor data and structure: the bond-related chemical shift concept. In JA Tossell, editor,Nuclear magnetic shieldings and molecular structure, volume 386, pages 191–201, 1993

  65. [65]

    Rovibrational effects, temperature dependence, and isotope effects on the nuclear shielding tensors of water: A new 17O absolute shielding scale.J

    Juha Vaara, Juhani Lounila, Kenneth Ruud, and Trygve Helgaker. Rovibrational effects, temperature dependence, and isotope effects on the nuclear shielding tensors of water: A new 17O absolute shielding scale.J. Chem. Phys., 109(19):8388–8397, 1998

  66. [66]

    Relative importance of first and second derivatives of nuclear magnetic resonance chemical shifts and spin-spin coupling constants for vibrational averaging.J

    Martin Dračínský, Jakub Kaminský, and Petr Bouř. Relative importance of first and second derivatives of nuclear magnetic resonance chemical shifts and spin-spin coupling constants for vibrational averaging.J. Chem. Phys., 130(9):094106, 2009

  67. [67]

    Böhm, Joachim Schulte, and Rafael Ramírez

    Michael C. Böhm, Joachim Schulte, and Rafael Ramírez. On the influence of nuclear fluctuations on calculated NMR shieldings of benzene and ethylene: a Feynman path integral-ab initio investigation.Int. J. Quantum Chem., 86(3):280–296, 2001

  68. [68]

    Joachim Schulte, Rafael Ramírez, and Michael C. Böhm. Influence of nuclear fluctuations on the NMR parameters of bullvalene: A Feynman path integral – Ab initio study.Chem. Phys. Lett., 432(4-6):579–584, 2006. 19

  69. [69]

    The chemical shift of deproto- nated water dimer: Ab initio path integral simulation.J

    Motoyuki Shiga, Kimichi Suzuki, and Masanori Tachikawa. The chemical shift of deproto- nated water dimer: Ab initio path integral simulation.J. Chem. Phys., 132(11):114104, 2010

  70. [70]

    Effects of Quantum Nuclear Delocalisation on NMR Parameters from Path Integral Molecular Dynamics.Chem

    Martin Dračínský and Paul Hodgkinson. Effects of Quantum Nuclear Delocalisation on NMR Parameters from Path Integral Molecular Dynamics.Chem. Eur. J., 20(8):2201–2207, 2014

  71. [71]

    Hassanali, Jérôme Cuny, Michele Ceriotti, Chris J

    Ali A. Hassanali, Jérôme Cuny, Michele Ceriotti, Chris J. Pickard, and Michele Par- rinello. The fuzzy quantum proton in the hydrogen chloride hydrates.J. Am. Chem. Soc., 134(20):8557–8569, 2012

  72. [72]

    Jürgen Gauss and John F. Stanton. Coupled-cluster calculations of nuclear magnetic resonance chemical shifts.J. Chem. Phys., 103(9):3561, 1995

  73. [73]

    J. S. Waugh, L. M. Huber, and U. Haeberlen. Approach to high-resolution nmr in solids. Phys. Rev. Lett., 20(5):180–182, 1968

  74. [74]

    S. G. Kukolich. Proton magnetic shielding tensors from spin-rotation measurements on formaldehyde and ammonia.J. Am. Chem. Soc., 97(20):5704–5707, 1975

  75. [75]

    The proton nuclear magnetic shielding tensors in biphenyl: Experiment and theory.J

    Frank Schönborn, Heike Schmitt, Herbert Zimmermann, Ulrich Haeberlen, Clémence Corminboeuf, Gisbert Großmann, and Thomas Heine. The proton nuclear magnetic shielding tensors in biphenyl: Experiment and theory.J. Magn. Reson., 175(1):52–64, 2005

  76. [76]

    Proton magnetic shielding tensor in liquid water.J

    Kristofer Modig and Bertil Halle. Proton magnetic shielding tensor in liquid water.J. Am. Chem. Soc., 124(40):12031–12041, 2002. 20