pith. sign in

arxiv: 2602.10267 · v2 · submitted 2026-02-10 · ⚛️ physics.chem-ph

Constrained nuclear-electronic orbital second-order Moller-Plesset perturbation theory

Pith reviewed 2026-05-16 02:13 UTC · model grok-4.3

classification ⚛️ physics.chem-ph
keywords constrained nuclear-electronic orbitalMP2 perturbation theorynuclear quantum effectsvibrational averagingmolecular propertiesmulticomponent methodselectronic-nuclear correlation
0
0 comments X

The pith

CNEO-MP2 captures nuclear vibrational effects on molecular properties in one calculation by adding electronic-nuclear correlation.

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

The paper derives a multicomponent second-order Møller-Plesset perturbation theory method inside the constrained nuclear-electronic orbital framework, starting from a multicomponent generalization of the Hylleraas functional. This CNEO-MP2 approach includes both electronic-nuclear and nuclear correlation to compute vibrationally averaged properties such as internuclear distances, bond angles, potential energy surfaces, and vibrational frequencies. A sympathetic reader would care because the method obtains nuclear quantum effects including zero-point energy and isotopic shifts directly during a single geometry optimization or energy evaluation, removing the need for separate higher-order force-constant calculations required by many existing techniques.

Core claim

Derived from a multicomponent generalization of the Hylleraas functional, the CNEO-MP2 method includes electronic-nuclear and nuclear correlation to calculate vibrationally averaged molecular properties within the constrained nuclear-electronic orbital framework, as demonstrated by benchmarks on diatomic and small polyatomic molecules and ions.

What carries the argument

Multicomponent generalization of the Hylleraas functional, used to generate the second-order perturbation expansion inside the CNEO framework.

If this is right

  • CNEO-MP2 produces internuclear distances and bond angles that already incorporate vibrational averaging.
  • It generates potential energy surfaces and vibrational frequencies that reflect nuclear quantum effects directly.
  • Nuclear quantum effects such as zero-point energy and isotopic shifts are obtained without additional post-processing force-constant calculations.
  • The method applies to both diatomic and small polyatomic molecules and ions in a single computation.

Where Pith is reading between the lines

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

  • The single-step workflow could lower the cost of including nuclear quantum effects when screening larger sets of molecules.
  • Routine use in quantum chemistry codes might make vibrational averaging a default option rather than an optional add-on.
  • The approach might be extended to compute averaged properties along reaction paths where zero-point corrections matter.

Load-bearing premise

The multicomponent generalization of the Hylleraas functional supplies a valid starting point for a second-order perturbation expansion inside the constrained nuclear-electronic orbital framework.

What would settle it

Benchmark calculations on the test set of molecules and ions showing that CNEO-MP2 internuclear distances or vibrational frequencies deviate significantly from experimental values or from established methods that include nuclear quantum effects would falsify the central claim.

Figures

Figures reproduced from arXiv: 2602.10267 by Gabrielle B. Tucker, Kurt R. Brorsen.

Figure 1
Figure 1. Figure 1: Single-component MP2 potential energy surface of the bifluoride anion (FHF− ) as a function of the position of the hydrogen atom in the xz-plane. The fluorine atoms are fixed at 1.146 and -1.146 ˚A [PITH_FULL_IMAGE:figures/full_fig_p024_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Single-component MP2 potential energy surface of the deuterium bifluoride anion (FDF− ) as a function of the position of the deuterium atom in the xz-plane. The fluorine atoms are fixed at 1.146 and -1.146 ˚A [PITH_FULL_IMAGE:figures/full_fig_p024_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Multicomponent CNEO-MP2 potential energy surface of the bifluoride anion (FHF− ) as a function of the position of the hydrogen atom in the xz-plane. Here the fluorine atoms are fixed at 1.146 and -1.146 ˚A. 24 [PITH_FULL_IMAGE:figures/full_fig_p024_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Multicomponent CNEO-MP2 potential energy surface of the deuterium bifluoride anion (FDF− ) as a function of the position of the deuterium atom in the xz-plane. Here the fluorine atoms are fixed at 1.146 and -1.146 ˚A. 5.4 The Zundel Cation: Geometry and Frequencies The Zundel cation (H5O+ 2 ), is an important system for studying proton transfer and other effects of nuclear vibration due to the relatively f… view at source ↗
Figure 5
Figure 5. Figure 5: Geometry of the Zundel (H5O+ 2 ) cation, calculated with CNEO-MP2, using the aug-cc-pVTZ electronic basis set, and the PB4-D nuclear basis set for the shared proton. The optimized structure has C2 symmetry. This visualization was generated with Avogadro 2.0 [132]. The blue arrow corresponds to the positive z-axis, the red arrow to the positive x-axis, and the green arrow to the positive y-axis. 26 [PITH_F… view at source ↗
read the original abstract

A multicomponent second-order M{\o}ller-Plesset perturbation theory (MP2) method is derived and implemented within the constrained nuclear-electronic orbital (CNEO) framework from a multicomponent generalization of the Hylleraas functional. The CNEO-MP2 method includes electronic-nuclear and nuclear correlation in the calculation of vibrationally averaged molecular properties. Nuclear quantum effects like vibrational averaging, isotopic effects, and zero-point energy can be captured in a single calculation or geometry optimization with CNEO-MP2, eliminating the need to perform costly subsequent calculations to determine higher order force constants as required with many existing methods used to determine vibrational effects upon molecular properties. The CNEO-MP2 method is benchmarked on a test set of diatomic and small polyatomic molecules and ions. Herein, we present internuclear distances, bond angles, potential energy surfaces, and vibrational frequencies calculated with the CNEO-MP2 method to demonstrate that it correctly captures the effects of nuclear vibrational motion upon molecular properties.

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 derives and implements a multicomponent second-order Møller-Plesset perturbation theory (CNEO-MP2) within the constrained nuclear-electronic orbital (CNEO) framework, starting from a multicomponent generalization of the Hylleraas functional. The method incorporates electronic-nuclear and nuclear correlation to compute vibrationally averaged molecular properties such as internuclear distances, bond angles, potential energy surfaces, and vibrational frequencies. Benchmarks on diatomic and small polyatomic molecules and ions are presented to demonstrate capture of nuclear quantum effects including vibrational averaging, isotopic shifts, and zero-point energy in a single calculation or geometry optimization, without requiring separate higher-order force constant evaluations.

Significance. If the central derivation is sound, the CNEO-MP2 approach offers a meaningful advance by embedding nuclear quantum effects directly into a perturbative electronic-structure framework, potentially reducing the computational overhead associated with post hoc vibrational averaging or anharmonic corrections in conventional methods. This could prove useful for systems where nuclear motion significantly influences properties, such as hydrogen-bonded complexes or isotopically sensitive reactions.

major comments (2)
  1. [§2.3, Eq. (12)] §2.3, Eq. (12): The multicomponent Hylleraas functional is generalized to CNEO, but the text does not explicitly introduce Lagrange multipliers for the nuclear-electronic orbital constraints or demonstrate that the first-order stationarity conditions remain satisfied after the MP2 amplitudes are solved. Without this step, it is unclear whether the second-order energy correction remains strictly within the constrained manifold or admits unphysical contributions that do not correspond to vibrational averaging.
  2. [§4.2, Table 2] §4.2, Table 2: The reported mean absolute deviations for vibrational frequencies (approximately 15 cm⁻¹ across the test set) are presented without direct comparison to CCSD(T) or experimental values with propagated uncertainties; this weakens the claim that CNEO-MP2 “correctly captures” the effects, as the magnitude of improvement over CNEO-HF or other baselines is not quantified.
minor comments (2)
  1. [§2.1] The definition of the nuclear-electronic orbital basis functions in §2.1 is introduced without an accompanying table of symbols, making it difficult to track the distinction between electronic and nuclear orbital indices throughout the amplitude equations.
  2. [Figure 3] Figure 3 caption states that the potential energy surface is “vibrationally averaged,” but the figure itself shows only the raw CNEO-MP2 curve; a second panel or inset comparing the averaged versus unaveraged surface would improve clarity.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the positive assessment of the CNEO-MP2 method and its potential utility for incorporating nuclear quantum effects. We address the two major comments below with clarifications and planned revisions.

read point-by-point responses
  1. Referee: [§2.3, Eq. (12)] §2.3, Eq. (12): The multicomponent Hylleraas functional is generalized to CNEO, but the text does not explicitly introduce Lagrange multipliers for the nuclear-electronic orbital constraints or demonstrate that the first-order stationarity conditions remain satisfied after the MP2 amplitudes are solved. Without this step, it is unclear whether the second-order energy correction remains strictly within the constrained manifold or admits unphysical contributions that do not correspond to vibrational averaging.

    Authors: The CNEO framework enforces the nuclear-electronic orbital constraints at the mean-field (CNEO-HF) level via Lagrange multipliers that define the constrained orbitals. The multicomponent Hylleraas functional is then constructed and minimized within this fixed constrained orbital basis, so the MP2 amplitudes are solved in the space orthogonal to the constraints. This construction ensures the second-order correction remains strictly on the constrained manifold by design. We acknowledge that the manuscript does not explicitly restate the stationarity conditions after the amplitudes are obtained; we will add a short paragraph and an auxiliary equation in §2.3 of the revised manuscript to demonstrate that the first-order stationarity conditions continue to hold. revision: partial

  2. Referee: [§4.2, Table 2] §4.2, Table 2: The reported mean absolute deviations for vibrational frequencies (approximately 15 cm⁻¹ across the test set) are presented without direct comparison to CCSD(T) or experimental values with propagated uncertainties; this weakens the claim that CNEO-MP2 “correctly captures” the effects, as the magnitude of improvement over CNEO-HF or other baselines is not quantified.

    Authors: We agree that direct comparisons and quantification of improvement would strengthen the presentation. In the revised manuscript we will augment Table 2 (and the accompanying text) with CCSD(T) reference values for the subset of molecules where they are available, experimental vibrational frequencies, and the corresponding mean absolute deviations for CNEO-HF. This will allow explicit quantification of the improvement provided by the MP2 correlation terms. revision: yes

Circularity Check

0 steps flagged

CNEO-MP2 derivation from multicomponent Hylleraas generalization is self-contained with no reduction to inputs

full rationale

The paper derives the CNEO-MP2 method explicitly from a multicomponent generalization of the Hylleraas functional inside the CNEO framework, then benchmarks it on diatomic and polyatomic test sets for properties such as internuclear distances and vibrational frequencies. No equations or steps in the abstract or described derivation reduce a prediction to a fitted parameter, self-definition, or self-citation chain; the Hylleraas starting point is treated as an external functional whose stationarity conditions are assumed to hold under CNEO constraints, and the MP2 amplitudes are obtained by standard perturbation theory. The central claim therefore retains independent content from the functional generalization and the subsequent implementation, yielding a normal non-circular finding.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

Review performed on abstract only; no explicit free parameters, axioms, or invented entities are stated in the provided text.

pith-pipeline@v0.9.0 · 5476 in / 1060 out tokens · 123526 ms · 2026-05-16T02:13:11.538134+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

138 extracted references · 138 canonical work pages

  1. [1]

    D. M. Bishop. Molecular motion in electric fields.Rev. Mod. Phys.1990,62(2), 343–374

  2. [2]

    R. P. A. Bettens. Bound State Potential Energy Surface Construction: Ab Initio Zero-point Energies and Vibrationally Averaged Rotational Constants.J. Am. Chem. Soc.2003,125(2), 584–587

  3. [3]

    Puzzarini and J

    C. Puzzarini and J. F. Stanton. Connections between the accuracy of rotational constants and equi- librium molecular structures.Phys. Chem. Chem. Phys.2023,25(3), 1421–1429. 30

  4. [4]

    Solomon and J

    P. Solomon and J. M. Schulman. The Effects of Vibrational Averaging on the Calculated Nuclear Spin–Spin Coupling Constants of Ammonia.J. Am. Chem. Soc.1977,99(24), 7776–7779

  5. [5]

    Egidi, J

    F. Egidi, J. Bloino, C. Cappelli, V. Barone, and J. Tomasi. Tuning of NMR and EPR parameters by vibrational averaging and environmental effects: an integrated computational approach.Mol. Phys. 2013,111(9-11), 1345–1354

  6. [6]

    Grigoleit and M

    S. Grigoleit and M. B¨ uhl. Thermal Effects and Vibrational Corrections to Transition Metal NMR Chemical Shifts.Chem. Eur. J.2004,10(21), 5541–5552

  7. [7]

    T. Tuttle. Averaging Semiempirical NMR Chemical Shifts: Dynamic Effects on the Subpicosecond Time Scale.J. Phys. Chem. A2009,113(43), 11723–11733

  8. [8]

    Draˇ c´ ınsk´ y and P

    M. Draˇ c´ ınsk´ y and P. Hodgkinson. A molecular dynamics study of the effects of fast molecular motions on solid-state NMR parameters.Cryst. Eng. Comm.2013,15(43), 8705–8712

  9. [9]

    Benzi, O

    C. Benzi, O. Crescenzi, M. Pavone, and V. Barone. Reliable NMR chemical shifts for molecules in solution by methods rooted in density functional theory.Magn. Reson. Chem.2004,42(S1), S57–S67

  10. [10]

    Ishii, T

    Y. Ishii, T. Terao, and S. Hayashi. Theory and simulation of vibrational effects on structural mea- surements by solid-state nuclear magnetic resonance.J. Chem. Phys.1997,107(8), 2760–2774

  11. [11]

    Toyama, T

    M. Toyama, T. Oka, and Y. Morino. Effect of vibration and rotation on the internuclear distance.J. Mol. Spec.1964,13(1-4), 193–213

  12. [12]

    X. Chen, Z. Rinkevicius, Z. Cao, K. Ruud, and H. ˚Agren. Zero-point vibrational corrections to isotropic hyperfine coupling constants in polyatomic molecules.Phys. Chem. Chem. Phys.2011,13 (2), 696–707

  13. [13]

    B. C. Mort and J. Autschbach. Magnitude of Zero-Point Vibrational Corrections to Optical Rotation in Rigid Organic Molecules: A Time-Dependent Density Functional Study.J. Phys. Chem. A2005, 109(38), 8617–8623

  14. [14]

    J. M. Bowman. The Self-Consistent-Field Approach to Polyatomic Vibrations.Acc. Chem. Res.1986, 19(7), 202–208

  15. [15]

    Christiansen

    O. Christiansen. Vibrational structure theory: new vibrational wave function methods for calculation 31 of anharmonic vibrational energies and vibrational contributions to molecular properties.Phys. Chem. Chem. Phys.2007,9(23), 2942–2953

  16. [16]

    Shiga, M

    M. Shiga, M. Tachikawa, and S. Miura. A unified scheme forab initiomolecular orbital theory and path integral molecular dynamics.J. Chem. Phys.2001,115(20), 9149–9159

  17. [17]

    K. R. Glaesemann and L. E. Fried. A path integral approach to molecular thermochemistry.J. Chem. Phys.2003,118(4), 1596–1603

  18. [18]

    L´ opez-Ciudad, R

    T. L´ opez-Ciudad, R. Ram´ ırez, J. Schulte, M. C. B¨ ohm. Anharmonic effects on the structural and vibrational properties of the ethyl radical: A path integral Monte Carlo study.J. Chem. Phys.2003, 119(8), 4328–4338

  19. [19]

    H. H. Nielsen, The Vibration-Rotation Energies of Molecules.Rev. Mod. Phys.1951,23(2), 90–136

  20. [20]

    I. M. Mills. Vibrational perturbation theory.J. Mol. Spec.1961,5(1-6), 334–340

  21. [21]

    Neugebauer and B

    J. Neugebauer and B. A. Hess. Fundamental vibrational frequencies of small polyatomic molecules from density-functional calculations and vibrational perturbation theory.J. Chem. Phys2003,118 (16), 7215–7225

  22. [22]

    Bloino and V

    J. Bloino and V. Barone. A second-order perturbation theory route to vibrational averages and transition properties of molecules: General formulation and application to infrared and vibrational circular dichroism spectroscopies.J. Chem. Phys.2012,136(12), 124108-1–124108-15

  23. [23]

    P. R. Franke, J. F. Stanton, and G. E. Douberly. How to VPT2: Accurate and Intuitive Simula- tions of CH Stretching Infrared Specta Using VPT2+K with Large Effective Hamiltonian Resonance Treatments.J. Phys. Chem. A2021,125(6), 1301–1324

  24. [24]

    R. C. Fortenberry. Picking up Good Vibrations through Quartic Force Fields and Vibrational Per- turbation Theory.J. Phys. Chem. Lett2024,15(25), 6528–6537

  25. [25]

    M´ atyus

    E. M´ atyus. Pre-Born-Oppenheimer molecular structure theory.Mol. Phys.2019,117(5), 590–609

  26. [26]

    Tachikawa, K

    M. Tachikawa, K. Mori, H. Nakai, K. Iguchi. An extension of ab initio molecular orbital theory to nuclear motion.Chem. Phys. Lett.1998,290(4-6), 437–442. 32

  27. [27]

    S. Webb, T. Iordanov, and S. Hammes-Schiffer. Multiconfigurational nuclear-electronic orbital ap- proach: Incorporation of nuclear quantum effects in electronic structure calculations.J. Chem. Phys. 2002,117(9) 4106–4118

  28. [28]

    Reyes, M

    A. Reyes, M. V. Pak, and S. Hammes-Schiffer. Investigation of isotope effects with the nuclear- electronic orbital approach.J. Chem. Phys.2005,123(6), 064104-1–064104-8

  29. [29]

    R. U. Khan and R. Tonner-Zech. Geometrical Isotope Effects on Chemical Bonding in Hydrogen Bonded Systems: Combining Nuclear-Electronic Orbital DFT and Energy Decomposition Analysis. J. Comput. Chem.2025,46(24), 1–12

  30. [30]

    L. E. Smith, V. Briega-Martos, Y. Yang, and S. Hammes-Schiffer. Isotope Effects for Water at Pt(111) Computed with Nuclear-Electronic Orbital Theory.J. Chem. Phys. Lett.2025,16(51), 13054–13061

  31. [31]

    I. L. Thomas. Protonic Structure of Molecules. I. Ammonia Molecules.Phys. Rev.1969,185(1), 90–94

  32. [32]

    I. L. Thomas. Protonic Structure of Molecules. II. Methodology, Center-of-Mass Transformation, and the Structure of Methane, Ammonia, and Water.Phys. Rev. A1970,2(4), 1200–1208

  33. [33]

    H. Nakai. Simultaneous Determination of Nuclear and Electronic Wave Functions Without Born- Oppenheimer Approximation: Ab Initio NO+MO/HF Theory.Int. J. Quant. Chem.2002,86(6), 511–517

  34. [34]

    H. Nakai. Nuclear Orbital Plus Molecular Orbital Theory: Simultaneous Determination of Nuclear and Electronic Wave Functions Without Born–Oppenheimer Approximation.Int. J. Quant. Chem. 2007,107(14), 2849–2869

  35. [35]

    Pavoˇ sevi´ c, T

    F. Pavoˇ sevi´ c, T. Culpitt, and S. Hammes-Schiffer. Multicomponent Quantum Chemistry: Integrating Electronic and Nuclear Quantum Effects via the Nuclear–Electronic Orbital Method.Chem. Rev.2020, 120(9), 4222–4253

  36. [36]

    M. V. Pak, A. Chakraborty, and S. Hammes-Schiffer. Density Functional Theory Treatment of Elec- tron Correlation in the Nuclear-Electronic Orbital Approach.J. Phys. Chem. A2007,111(20), 4522–4526

  37. [37]

    Chakraborty, M

    A. Chakraborty, M. V. Pak, and S. Hammes-Schiffer. Properties of the exact universal functional in multicomponent density functional theory.J. Chem. Phys.2009,131(12), 124115-1–124115-8. 33

  38. [38]

    Chakraborty, M

    A. Chakraborty, M. V. Pak, and S. Hammes-Schiffer. Development of Electron-Proton Density Func- tionals for Multicomponent Density Functional Theory.Phys. Rev. Lett.2008,101(15), 153001- 1–153001-4; Erratum:Phys. Rev. Lett.2011,106(16), 169902-1

  39. [39]

    Y. Yang, K. R. Brorsen, T. Culpitt, M. V. Pak, and S. Hammes-Schiffer. Development of a prac- tical multicomponent density functional for electron-proton correlation to produce accurate proton densities.J. Chem. Phys.2017,147(11), 114113-1–114113-5

  40. [40]

    K. R. Brorsen, Y. Yang, and S. Hammes-Schiffer. Multicomponent Density Functional Theory: Im- pact of Nuclear Quantum Effects on Proton Affinities and Geometries.J. Chem. Phys. Lett.2017,8 (15), 3488–3493

  41. [41]

    K. R. Brorsen, P. E. Schneider, and S. Hammes-Schiffer. Alternative forms and transferability of electron-proton correlation functionals in nuclear-electronic orbital density functional theory.J. Chem. Phys.2018,149(4), 044110-1 –044110-7

  42. [42]

    Y. Yang, T. Culpitt, and S. Hammes-Schiffer. Multicomponent Time-Dependent Density Functional Theory: Proton and Electron Excitation Energies.J. Phys. Chem. Lett.2018,9(7), 1765–1770

  43. [43]

    Culpitt, Y

    T. Culpitt, Y. Yang, F. Pavoˇ sevi´ c, Z. Tao, and S. Hammes-Schiffer. Enhancing the applicability of multicomponent time-dependent density functional theory.J. Chem. Phys.2019,150(20), 201101- 1–201101-6

  44. [44]

    Z. Tao, Y. Yang, and S. Hammes-Schiffer. Multicomponent density functional theory: Including the density gradient in the electron-proton correlation functional for hydrogen and deuterium.J. Chem. Phys.2019,151(12), 124102-1–124102-7

  45. [45]

    Hammes-Schiffer

    Q.Yu and S. Hammes-Schiffer. Nuclear-Electronic Orbital Multistate Density Functional Theory.J. Phys. Chem. Lett.2020,11(23), 10106–10113

  46. [46]

    Hasecke and R

    L. Hasecke and R. A. Mata. Multicomponent Double-Hybrid Density Functional Theory.J. Chem. Theory Comput.2025,21(22), 11509–11520

  47. [47]

    Gimferrer, L

    M. Gimferrer, L. Hasecke, M. B¨ odecker, and R. A. Mata. Accurate vibrational hydrogen-bond shift predictions with multicomponent DFT.Chem. Sci.2025,16(24), 11002–11011. 34

  48. [48]

    Nakai, K

    H. Nakai, K. Sodeyama, and M. Hoshino. Non-Born–Oppenheimer theory for simultaneous determi- nation of vibrational and electronic excited states: ab initio NO+MO/CIS theory.Chem. Phys. Lett. 2001,345(1-2), 118–124

  49. [49]

    Nakai and K

    H. Nakai and K. Sodeyama. Many-body effects in nonadiabatic molecular theory for simultaneous determination of nuclear and electronic wave functions:Ab initioNOMO/MBPT and CC methods. J. Chem. Phys.2003,118(3), 1119–1127

  50. [50]

    Iordanov and S

    T. Iordanov and S. Hammes-Schiffer. Vibrational analysis for the nuclear–electronic orbital method. J. Chem. Phys.2003,118(21), 9489–9496

  51. [51]

    Swalina, M

    C. Swalina, M. V. Pak, and S. Hammes-Schiffer. Alternative formulation of many-body perturbation theory for electron–proton correlation.Chem. Phys. Lett.2005,404(4–6), 394–399

  52. [52]

    J. H. Skone, M. V. Pak, and S. Hammes-Schiffer. Nuclear-electronic orbital nonorthogonal configu- ration interaction approach.J. Chem. Phys.2005,123(13), 134108-1–134108-8

  53. [53]

    Auer and S

    B. Auer and S. Hammes-Schiffer. Localized Hartree product treatment of multiple protons in the nuclear-electronic orbital framework.J. Chem. Phys.2010,132(8), 084110-1–084110-8

  54. [54]

    Pavoˇ sevi´ c, T

    F. Pavoˇ sevi´ c, T. Culpitt, and S. Hammes-Schiffer. Multicomponent Coupled Cluster Singles and Doubles Theory within the Nuclear-Electronic Orbital Framework.J. Chem. Theory Comp.2019,15 (1), 338–347

  55. [55]

    Pavoˇ sevi´ c and S

    F. Pavoˇ sevi´ c and S. Hammes-Schiffer. Multicomponent coupled cluster singles and doubles and Brueckner doubles methods: Proton densities and energies.J. Chem. Phys.2019,151(7), 074104- 1–074104-7

  56. [56]

    Pavoˇ sevi´ c, B

    F. Pavoˇ sevi´ c, B. J. G. Rousseau, and S. Hammes-Schiffer. Multicomponent Orbital-Optimized Per- turbation Theory Methods: Approaching Coupled Cluster Accuracy at Lower Cost.J. Phys. Chem. Lett.2020,11(4), 1578–1583

  57. [57]

    O. J. Fajen and K. R. Brorsen. Separation of electron-electron and electron-proton correlation in mul- ticomponent orbital-optimized perturbation theory.J. Chem. Phys2020,152(19), 194107-1–194107-8

  58. [58]

    Pavoˇ sevi´ c, Z

    F. Pavoˇ sevi´ c, Z. Tao, and S. Hammes-Schiffer. Multicomponent Coupled Cluster Singles and Doubles with Density Fitting: Protonated Water Tetramers with Quantized Protons.J. Phys. Chem. Lett. 2021,12(6), 1631–1637. 35

  59. [59]

    Alaal and K

    N. Alaal and K. R. Brorsen. Multicomponent heat-bath configuration interaction with the perturba- tive correction for the calculation of protonic excited states.J. Chem. Phys.2021,155(23), 234107- 1–234107-11

  60. [60]

    O. J. Fajen, and K. R. Brorsen. Multicomponent MP4 and the inclusion of triple excitations in multicomponent many-body methods.J. Chem. Phys.2021,155(23), 234108-1–234108-9

  61. [61]

    Fowler and K

    D. Fowler and K. R. Brorsen. (T) Correction for Multicomponent Coupled-Cluster Theory for a Single Quantum Proton.J. Chem. Theory Comput.2022,18(12), 7298–7305. Erratum:J. Chem. Theory Comput.2025,21(9), 5005–5007

  62. [62]

    Hasecke and R

    L. Hasecke and R. A. Mata. Nuclear Quantum Effects Made Accessible: Local Density Fitting in Multicomponent Methods.J. Chem. Theory Comput.2023,19(22), 8223–8233

  63. [63]

    Hasecke and R

    L. Hasecke and R. A. Mata. Local Electronic Correlation in Multicomponent Møller-Plesset Pertur- bation Theory.J. Chem. Theory Comput.2023,20(22), 9928–9938

  64. [64]

    R. J. Goudy, F. Pavoˇ sevi´ c, and S. Hammes-Schiffer. Triple excitations in nuclear–electronic orbital coupled cluster theory for multiple quantum protons.J. Chem. Phys.2025,163(22), 224119-1–224119- 8

  65. [65]

    J. H. Fetherholf, F. Pavoˇ sevi´ c, and S. Hammes-Schiffer. Nuclear-electronic orbital second-order cou- pled cluster for excited states.arXiv Preprint2025, 1–10

  66. [66]

    R. M. Nieminen, E. Boro´ nski, and L. J. Lantto. Two-component density functional theory: Applica- tion to positron states.Phys. Rev. B1985,32(2), 1377–1379

  67. [67]

    Wiktor, G

    J. Wiktor, G. Jomard, and M. Torrent. Two-component density functional theory within the projector augmented-wave approach: Accurate self-consistent computations of positron lifetimes and momentum distributions.Phys. Rev. B2015,92(23), 125113-1–125113-15

  68. [68]

    Boro´ nski and R

    E. Boro´ nski and R. M. Nieminen. Electron-positron density-functional theory.Phys. Rev. B1986, 34(6), 3820–3831

  69. [69]

    D. V. Makhov and L. J. Lewis. Two-component density functional theory calculations of positron lifetimes for small vacancy clusters in silicon.Phys. Rev. B2005,71(20), 205215-1–205215-6. 36

  70. [70]

    Goli and S

    M. Goli and S. Shahbazian. Two-component density functional theory for muonic molecules: Inclusion of the electron–positive muon correlation functional.J. Chem. Phys.2022,156(4), 044104-1–044104- 16

  71. [71]

    N. S. Riyahi, M. Goli, and S. Shahbazian. Quantifying errors of electron-proton/muon correlation functionals through the Kohn-Sham inversion of a two-component model system.Phys. Rev. B2023, 108(24), 245155-1–245155-20

  72. [72]

    Udagawa, T

    T. Udagawa, T. Tsuneda, and M. Tachikawa. Electron-nucleus correlation functional for multicom- ponent density functional thoery.Phys. Rev. A2014,89(5), 052519-1–052519-6

  73. [73]

    Z. Pan, L. Peng, Z. Pan, Y. Yuan, H. Zhang, and B. Ye. F-µbond length andµSR depolarization spectrum calculation for fluoride using two-component density functional theory.Chinese Phys. B 2023,32(8), 087602-1–087602-4

  74. [74]

    Rayka, M

    M. Rayka, M. Goli, and S. Shahbazian. Effective electronic-only Kohn-Sham equations for the muonic molecules.Phys. Chem. Chem. Phys.2018,20(13), 8802–8811

  75. [75]

    Goli and S

    M. Goli and S. Shahbazian. Developing effective electronic-only coupled-cluster and Møller-Plesset perturbation theories for the muonic molecules.Phys. Chem. Chem. Phys.2018,20(24), 16749–16760

  76. [76]

    Rayka, M

    M. Rayka, M. Goli, and S. Shahbazian. Toward a muon-specific electronic structure theory: effective electronic Hartree-Fock equations for muonic molecules.Phys. Chem. Chem. Phys.2018,20(6), 4466–4477

  77. [77]

    B. H. Ellis, S. Aggarway, and A. Chakraborty. Development of the Multicomponent Coupled-Cluster Theory for Investigation of Multiexcitonic Interactions.J. Chem. Theory Comput.2016,12(1), 188–200

  78. [78]

    Nakai, M

    H. Nakai, M. Hoshino, K. Miyamoto, and S. Hyodo. Elimination of translational and rotational motions in nuclear orbital plus molecular orbital theory.J. Chem. Phys.2005,122(16), 164101- 1–164101-10

  79. [79]

    Miyamoto, M

    K. Miyamoto, M. Hoshino, and H. Nakai. Elimination of Translational and Rotational Motions in Nu- clear Orbital Plus Molecular Orbital Theory: Contribution of the First-Order Rovibration Coupling. J. Chem. Theory Comput.2006,2(6), 1544–1550. 37

  80. [80]

    A. D. Bochevearov, E. F. Valeev, and C. D. Sherrill. The electron and nuclear orbitals model: current challenges and future prospectsMol. Phys.2006,102(1), 111–123

Showing first 80 references.