Addressing first derivative discontinuities in orbital-optimised opposite-spin scaled second-order perturbation theory with regularisation

High-Performance Multi-GPU Analytic RI-MP2 Energy Gradients

This article presents a novel algorithm for the calculation of analytic energy gradients from second-order Møller-Plesset perturbation theory within the Resolution-of-the-Identity approximation (RI-MP2), which is designed to achieve high performance on clusters with multiple graphical processing units (GPUs). The algorithm uses GPUs for all major steps of the calculation, including integral generation, formation of all required intermediate tensors, solution of the Z-vector equation and gradient accumulation. The implementation in the EXtreme Scale Electronic Structure System (EXESS) software package includes a tailored, highly efficient, multistream scheduling system to hide CPU-GPU data transfer latencies and allows nodes with 8 A100 GPUs to operate at over 80% of theoretical peak floating-point performance. Comparative performance analysis shows a significant reduction in computational time relative to traditional multicore CPU-based methods, with our approach achieving up to a 95-fold speedup over the single-node performance of established software such as Q-Chem and ORCA. Additionally, we demonstrate that pairing our implementation with the molecular fragmentation framework in EXESS can drastically lower the computational scaling of RI-MP2 gradient calculations from quintic to subquadratic, enabling further substantial savings in runtime while retaining high numerical accuracy in the resulting gradients.

Reaching High Accuracy for Energetic Properties at Second-Order Perturbation Cost by Merging Self-Consistency and Spin-Opposite Scaling

Quantum chemical methods dealing with challenging systems while retaining low computational costs have attracted attention. In particular, many efforts have been devoted to developing new methods based on second-order perturbation that may be the simplest correlated method beyond Hartree-Fock. We have recently developed a self-consistent perturbation theory named one-body Møller-Plesset second-order perturbation theory (OBMP2) and shown that it can resolve issues caused by the noniterative nature of standard perturbation theory. In this work, we extend the method by introducing spin-opposite scaling to the double-excitation amplitudes, resulting in the O2BMP2 method. We assess the O2BMP2 performance on the triple-bond N2 dissociation, singlet-triplet gaps, and ionization potentials. O2BMP2 performs much better than standard MP2 and reaches the accuracy of coupled-cluster methods in all cases considered in this work.

Optimizing the regularization in size-consistent second-order Brillouin-Wigner perturbation theory

Despite its simplicity and relatively low computational cost, second-order Møller-Plesset perturbation theory (MP2) is well-known to overbind noncovalent interactions between polarizable monomers and some organometallic bonds. In such situations, the pairwise-additive correlation energy expression in MP2 is inadequate. Although energy-gap dependent amplitude regularization can substantially improve the accuracy of conventional MP2 in these regimes, the same regularization parameter worsens the accuracy for small molecule thermochemistry and density-dependent properties. Recently, we proposed a repartitioning of Brillouin-Wigner perturbation theory that is size-consistent to second order (BW-s2), and a free parameter (α) was set to recover the exact dissociation limit of H2 in a minimal basis set. Alternatively α can be viewed as a regularization parameter, where each value of α represents a valid variant of BW-s2, which we denote as BW-s2(α). In this work, we semi-empirically optimize α for noncovalent interactions, thermochemistry, alkane conformational energies, electronic response properties, and transition metal datasets, leading to improvements in accuracy relative to the ab initio parameterization of BW-s2 and MP2. We demonstrate that the optimal α parameter (α = 4) is more transferable across chemical problems than energy-gap-dependent regularization parameters. This is attributable to the fact that the BW-s2(α) regularization strength depends on all of the information encoded in the t amplitudes rather than just orbital energy differences. While the computational scaling of BW-s2(α) is iterative O(N5), this effective and transferable approach to amplitude regularization is a promising route to incorporate higher-order correlation effects at second-order cost.

Repartitioned Brillouin-Wigner perturbation theory with a size-consistent second-order correlation energy

Second-order Møller-Plesset perturbation theory (MP2) often breaks down catastrophically in small-gap systems, leaving much to be desired in its performance for myriad chemical applications such as noncovalent interactions, thermochemistry, and dative bonding in transition metal complexes. This divergence problem has reignited interest in Brillouin-Wigner perturbation theory (BWPT), which is regular at all orders but lacks size consistency and extensivity, severely limiting its application to chemistry. In this work, we propose an alternative partitioning of the Hamiltonian that leads to a regular BWPT perturbation series that, through the second order, is size-extensive, size-consistent (provided its Hartree-Fock reference is also), and orbital invariant. Our second-order size-consistent Brillouin-Wigner (BW-s2) approach can describe the exact dissociation limit of H2 in a minimal basis set, regardless of the spin polarization of the reference orbitals. More broadly, we find that BW-s2 offers improvements relative to MP2 for covalent bond breaking, noncovalent interaction energies, and metal/organic reaction energies, although rivaling coupled-cluster with single and double substitutions for thermochemical properties.

A Similarity Renormalization Group Approach to Green's Function Methods

The family of Green's function methods based on the GW approximation has gained popularity in the electronic structure theory thanks to its accuracy in weakly correlated systems combined with its cost-effectiveness. Despite this, self-consistent versions still pose challenges in terms of convergence. A recent study [Monino and Loos J. Chem. Phys. 2022, 156, 231101.] has linked these convergence issues to the intruder-state problem. In this work, a perturbative analysis of the similarity renormalization group (SRG) approach is performed on Green's function methods. The SRG formalism enables us to derive, from first-principles, the expression of a naturally static and Hermitian form of the self-energy that can be employed in quasiparticle self-consistent GW (qsGW) calculations. The resulting SRG-based regularized self-energy significantly accelerates the convergence of qsGW calculations, slightly improves the overall accuracy, and is straightforward to implement in existing code.

An in-silico NMR laboratory for nuclear magnetic shieldings computed via finite fields: Exploring nucleus-specific renormalizations of MP2 and MP3

We developed and implemented a method-independent, fully numerical, finite difference approach to calculating nuclear magnetic resonance shieldings, using gauge-including atomic orbitals. The resulting capability can be used to explore non-standard methods, given only the energy as a function of finite-applied magnetic fields and nuclear spins. For example, standard second-order Møller-Plesset theory (MP2) has well-known efficacy for 1H and 13C shieldings and known limitations for other nuclei such as 15N and 17O. It is, therefore, interesting to seek methods that offer good accuracy for 15N and 17O shieldings without greatly increased compute costs, as well as exploring whether such methods can further improve 1H and 13C shieldings. Using a small molecule test set of 28 species, we assessed two alternatives: κ regularized MP2 (κ-MP2), which provides energy-dependent damping of large amplitudes, and MP2.X, which includes a variable fraction, X, of third-order correlation (MP3). The aug-cc-pVTZ basis was used, and coupled cluster with singles and doubles and perturbative triples [CCSD(T)] results were taken as reference values. Our κ-MP2 results reveal significant improvements over MP2 for 13C and 15N, with the optimal κ value being element-specific. κ-MP2 with κ = 2 offers a 30% rms error reduction over MP2. For 15N, κ-MP2 with κ = 1.1 provides a 90% error reduction vs MP2 and a 60% error reduction vs CCSD. On the other hand, MP2.X with a scaling factor of 0.6 outperformed CCSD for all heavy nuclei. These results can be understood as providing renormalization of doubles amplitudes to partially account for neglected triple and higher substitutions and offer promising opportunities for future applications.

Accurate Property Prediction by Second Order Perturbation Theory: The REMP and OO-REMP Hybrids

The prediction of molecular properties such as equilibrium structures or vibrationalwavenumbers is a routine task in computational chemistry. If very high accuracy is required, however, the use of computationally demanding ab initio wavefunction methods is mandatory. We present property calculations utilizing the REMP and OO-REMP hybrid perturbation theories showing that with the latter approach, very accurate results are obtained at second order in perturbation theory. Specifically, equilibrium structures and harmonic vibrational wavenumbers as well as dipole moments of closed and open shell molecules were calculated and compared to the best available experimental results or very accurate calculations.OO-REMP is capable of predicting bond lengths of small closed and open shell molecules with an accuracy of 0.2 pm and 0.5 pm, respectively, often within the range of experimental uncertainty. Equilibrium harmonic vibrational wavenumbers are predicted with an accuracy better than 20 cm−1 . Dipole moments of small closed and open shell molecules are reproduced with a relative error of less than 3 %. Across all investigated properties it turns out that a 20 %:80 % MP:RE mixing ratio consistently provides the best results. This is in line with our previous findings featuring closed and open shell reaction energies.

Unphysical discontinuities, intruder states and regularization in GW methods

By recasting the non-linear frequency-dependent GW quasiparticle equation into a linear eigenvalue problem, we explain the appearance of multiple solutions and unphysical discontinuities in various physical quantities computed within the GW approximation. Considering the GW self-energy as an effective Hamiltonian, it is shown that these issues are key signatures of strong correlation in the (N ± 1)-electron states and can be directly related to the intruder state problem. A simple and efficient regularization procedure inspired by the similarity renormalization group is proposed to avoid such issues and speed up the convergence of partially self-consistent GW calculations.

Do Double-Hybrid Functionals Benefit from Regularization in the PT2 Term? Observations from an Extensive Benchmark

We put to the test a recent suggestion [Shee, J., et al. J. Phys. Chem. Lett. 2021, 12 (50), 12084-12097] that MP2 regularization might improve the performance of double-hybrid density functionals. Using the very large and chemically diverse GMTKN55 benchmark, we find that κ-regularization is indeed beneficial at lower percentages of Hartree-Fock exchange, especially if spin-component scaling is not applied [such as in B2GP-PLYP or ωB97M(2)]. This benefit dwindles for DSD and DOD functionals and vanishes entirely in the ∼70% HF exchange region optimal for them.

UREMP, RO-REMP, and OO-REMP: Hybrid perturbation theories for open-shell electronic structure calculations

An accurate description of the electron correlation energy in closed- and open-shell molecules is shown to be obtained by a second-order perturbation theory (PT) termed REMP. REMP is a hybrid of the Retaining the Excitation degree (RE) and the Møller–Plesset (MP) PTs. It performs particularly encouragingly in an orbital-optimized variant (OO-REMP) where the reference wavefunction is given by an unrestricted Slater determinant whose spin orbitals are varied such that the total energy becomes a minimum. While the approach generally behaves less satisfactorily with unrestricted Hartree–Fock references, reasonable performance is observed for restricted Hartree–Fock and restricted open-shell Hartree–Fock references. Inclusion of single excitations to OO-REMP is investigated and found—as in similar investigations—to be dissatisfying as it deteriorates performance. For the non-multireference subset of the accurate W4-11 benchmark set of Karton et al. [Chem. Phys. Lett. 510, 165–178 (2011)], OO-REMP predicts most atomization and reaction energies with chemical accuracy (1 kcal mol−1) if complete-basis-set extrapolation with augmented and core-polarized basis sets is used. For the W4-11 related test-sets, the error estimates obtained with the OO-REMP method approach those of coupled-cluster with singles, doubles and perturbative triples [CCSD(T)] within 20%–35%. The best performance of OO-REMP is found for a mixing ratio of 20%:80% MP:RE, which is essentially independent of whether radical stabilization energies, barrier heights, or reaction energies are investigated. Orbital optimization is shown to improve the REMP approach for both closed and open shell cases and outperforms coupled-cluster theory with singles and doubles (CCSD), spin-component scaled Møller-Plesset theory at second order (SCS-MP2), and density functionals, including double hybrids in all the cases considered.

Can second-order perturbation theory accurately predict electron density of open-shell molecules? The importance of self-consistency

Electron density plays an essential role in predicting molecular properties.

Regularized Second-Order Møller-Plesset Theory: A More Accurate Alternative to Conventional MP2 for Noncovalent Interactions and Transition Metal Thermochemistry for the Same Computational Cost

Second-order Møller-Plesset theory (MP2) notoriously breaks down for π-driven dispersion interactions and dative bonds in transition metal complexes. Herein, we investigate three physically justified forms of single-parameter, energy-gap dependent regularization which can yield high and transferable accuracy for a variety of noncovalent interactions (including S22, S66, and L7 test sets) and (mostly closed shell) transition metal thermochemistry. Regularization serves to damp overestimated pairwise additive contributions, renormalizing first-order amplitudes such that the effects of higher-order correlations are incorporated. The optimal parameter values for the noncovalent and transition metal sets are 1.1, 0.7, and 0.4 for κ, σ, and σ² regularizers, respectively. However, such regularization slightly degrades the accuracy of conventional MP2 for some small-molecule test sets, most of which have relatively large average frontier energy gaps. Our results suggest that appropriately regularized MP2 models may improve double hybrid density functionals, at no additional cost over conventional MP2.

Improving Perturbation Theory for Open-Shell Molecules via Self-Consistency

We present an extension of our one-body Møller-Plesset second-order perturbation (OBMP2) method for open-shell systems. We derived the OBMP2 Hamiltonian through the canonical transformation followed by the cumulant approximation to reduce many-body operators into one-body ones. The resulting Hamiltonian consists of an uncorrelated Fock (unperturbed Hamiltonian) and a one-body correlation potential (perturbed Hamiltonian) composed of only double excitations. Molecular orbitals and associated energy levels are then relaxed via self-consistency, similar to Hartree-Fock, in the presence of the correlation at the MP2 level. We demonstrate the OBMP2 performance by considering two examples well-known for requiring orbital optimization: bond breaking and isotropic hyperfine coupling constants. In contrast to noniterative MP2, we show that OBMP2 can yield a smooth transition through the unrestriction point and accurately predict isotropic hyperfine coupling constants.

Third-Order Møller-Plesset Perturbation Theory Made Useful? Choice of Orbitals and Scaling Greatly Improves Accuracy for Thermochemistry, Kinetics, and Intermolecular Interactions

We develop and test methods that include second- and third-order perturbation theory (MP3) using orbitals obtained from regularized orbital-optimized second-order perturbation theory, κ-OOMP2, denoted as MP3:κ-OOMP2. Testing MP3:κ-OOMP2 shows RMS errors that are 1.7-5 times smaller than those of MP3 across 7 data sets. To do still better, empirical training of the scaling factors for the second- and third-order correlation energies and the regularization parameter on one of those data sets led to an unregularized scaled (c₂ = 1.0; c₃ = 0.8) denoted as MP2.8:κ-OOMP2. MP2.8:κ-OOMP2 yields significant additional improvement over MP3:κ-OOMP2 in 4 of 6 test data sets on thermochemistry, kinetics, and noncovalent interactions. Remarkably, these two methods outperform coupled cluster with singles and doubles in 5 of the 7 data sets considered, at greatly reduced cost (no O(N6) iterations).

Distinguishing artificial and essential symmetry breaking in a single determinant: approach and application to the C60, C36, and C20 fullerenes

The existence of a generalized Hartree–Fock solution in C₆₀ has led to controversy on whether C₆₀ is polyradicaloid (or strongly correlated). We attempt to end the controversy with κ-OOMP2 which removes the illusion of this artificial symmetry breaking. We conclude that C₆₀ is not strongly correlated.