We propose progressive downsampling of wave vectors in the Brillouin zone integrations occurring in the second-order many-body or Møller-Plesset perturbation (MP2) method for extended systems of one-dimensional periodicity. Higher-lying unoccupied and lower-lying occupied Bloch orbitals are subject to downsampling by an exponentially increasing factor (with base n), making the total number of Bloch orbitals included in the MP2 lattice sums grow only logarithmically with respect to the number of basis functions per unit cell. Unlike the mod n downsampling scheme proposed earlier, this log n scheme reduces the scaling of the computational cost and thus achieves a greater speedup as the unit cell size increases. Correct band indexing is essential for accuracy. Two-electron integrals entering the MP2 energy and quasiparticle energy expressions must be multiplied by quadrature weights that are a function of the energy bands involved, and an algorithm to compute the weights is proposed. A combined use of the log n and mod n schemes can speedup the MP2/6-31G** calculation of polyacetylene typically by a factor of 20 with an error in the correlation energy within a few percent relative to the conventional calculation. Similar combinations can reproduce the MP2 quasiparticle energy bands accurately at a fraction of the usual computational cost.
We propose progressive downsampling of wave vectors in the Brillouin zone integrations occurring in the second-order many-body or Møller-Plesset perturbation (MP2) method for extended systems of one-dimensional periodicity. Higher-lying unoccupied and lower-lying occupied Bloch orbitals are subject to downsampling by an exponentially increasing factor (with base n), making the total number of Bloch orbitals included in the MP2 lattice sums grow only logarithmically with respect to the number of basis functions per unit cell. Unlike the mod n downsampling scheme proposed earlier, this log n scheme reduces the scaling of the computational cost and thus achieves a greater speedup as the unit cell size increases. Correct band indexing is essential for accuracy. Two-electron integrals entering the MP2 energy and quasiparticle energy expressions must be multiplied by quadrature weights that are a function of the energy bands involved, and an algorithm to compute the weights is proposed. A combined use of the log n and mod n schemes can speedup the MP2/6-31G** calculation of polyacetylene typically by a factor of 20 with an error in the correlation energy within a few percent relative to the conventional calculation. Similar combinations can reproduce the MP2 quasiparticle energy bands accurately at a fraction of the usual computational cost.
Explicit expressions for electron correlation at the second-order many-body perturbation-theory [MBPT(2)] level are presented and implemented for the total energy per unit cell and for the band structure of extended systems. In the latter case, a formula is presented for a direct evaluation of the band gap rather than obtaining it as a difference of two large numbers. Application is made to alternating trans-polyacetylene. We assess the convergence of MBPT(2) with the number of unit cells (N) included in the lattice summations, the number of k-points (K) taken for the integrals over k in the first Brillouin zone, and the cutoff threshold (10−C) for the two-electron integrals. The MBPT(2) correlation correction to the band structure converges very slowly with N and demands a large K while the MBPT(2) correction to the total energy per unit cell converges much faster with N and needs a much smaller K. Neither MBPT(2) correction is sensitive to the cutoff of the two-electron atomic orbital integrals, 10−C, when C≥5. For polyacetylene, the MBPT(2) band gap is much improved from the SCF result, but does not agree with previous numerical data. Analysis shows that the previous MBPT(2) results were obtained either with too few unit cells such that the convergence with N had not been reached, or that the zeroth-order Hartree–Fock results were inadequately converged. MBPT(2) with a DZP basis improves the Hartree–Fock band gap from 5.57 to 3.22 eV at the experimentally estimated geometry, compared to the measured ∼2 eV peak in the absorption spectrum of the system. We verify our results with three independent programs. We also study the band gap as a function of geometry.
We present a formulation of the explicitly correlated second-order Møller-Plesset (MP2-F12) energy in which all nontrivial post-mean-field steps are formulated with linear computational complexity in system size. The two key ideas are the use of pair-natural orbitals for compact representation of wave function amplitudes and the use of domain approximation to impose the block sparsity. This development utilizes the concepts for sparse representation of tensors described in the context of the domain based local pair-natural orbital-MP2 (DLPNO-MP2) method by us recently [Pinski et al., J. Chem. Phys. 143, 034108 (2015)]. Novel developments reported here include the use of domains not only for the projected atomic orbitals, but also for the complementary auxiliary basis set (CABS) used to approximate the three- and four-electron integrals of the F12 theory, and a simplification of the standard B intermediate of the F12 theory that avoids computation of four-index two-electron integrals that involve two CABS indices. For quasi-1-dimensional systems (n-alkanes), the ON DLPNO-MP2-F12 method becomes less expensive than the conventional ON(5) MP2-F12 for n between 10 and 15, for double- and triple-zeta basis sets; for the largest alkane, C200H402, in def2-TZVP basis, the observed computational complexity is N(∼1.6), largely due to the cubic cost of computing the mean-field operators. The method reproduces the canonical MP2-F12 energy with high precision: 99.9% of the canonical correlation energy is recovered with the default truncation parameters. Although its cost is significantly higher than that of DLPNO-MP2 method, the cost increase is compensated by the great reduction of the basis set error due to explicit correlation.
We present a low-complexity algorithm to calculate the correlation energy of periodic systems in second-order Møller-Plesset (MP2) perturbation theory. In contrast to previous approximation-free MP2 codes, our implementation possesses a quartic scaling, O(N4), with respect to the system size N and offers an almost ideal parallelization efficiency. The general issue that the correlation energy converges slowly with the number of basis functions is eased by an internal basis set extrapolation. The key concept to reduce the scaling is to eliminate all summations over virtual orbitals which can be elegantly achieved in the Laplace transformed MP2 formulation using plane wave basis sets and fast Fourier transforms. Analogously, this approach could allow us to calculate second order screened exchange as well as particle-hole ladder diagrams with a similar low complexity. Hence, the presented method can be considered as a step towards systematically improved correlation energies.
In the Monte Carlo second-order many-body perturbation (MC-MP2) method, the long sum-of-product matrix expression of the MP2 energy, whose literal evaluation may be poorly scalable, is recast into a single high-dimensional integral of functions of electron pair coordinates, which is evaluated by the scalable method of Monte Carlo integration. The sampling efficiency is further accelerated by the redundant-walker algorithm, which allows a maximal reuse of electron pairs. Here, a multitude of graphical processing units (GPUs) offers a uniquely ideal platform to expose multilevel parallelism: fine-grain data-parallelism for the redundant-walker algorithm in which millions of threads compute and share orbital amplitudes on each GPU; coarse-grain instruction-parallelism for near-independent Monte Carlo integrations on many GPUs with few and infrequent interprocessor communications. While the efficiency boost by the redundant-walker algorithm on central processing units (CPUs) grows linearly with the number of electron pairs and tends to saturate when the latter exceeds the number of orbitals, on a GPU it grows quadratically before it increases linearly and then eventually saturates at a much larger number of pairs. This is because the orbital constructions are nearly perfectly parallelized on a GPU and thus completed in a near-constant time regardless of the number of pairs. In consequence, an MC-MP2/cc-pVDZ calculation of a benzene dimer is 2700 times faster on 256 GPUs (using 2048 electron pairs) than on two CPUs, each with 8 cores (which can use only up to 256 pairs effectively). We also numerically determine that the cost to achieve a given relative statistical uncertainty in an MC-MP2 energy increases as O ( n 3 ) or better with system size n, which may be compared with the O ( n 5 ) scaling of the conventional implementation of deterministic MP2. We thus establish the scalability of MC-MP2 with both system and computer sizes.
We present a domain-based local pair natural orbital Møller-Plesset second-order perturbation theory (DLPNO-MP2) for periodic systems, working within a linear combination of atomic orbitals formalism in the Turbomole program package. This approach, Megacell-DLPNO-MP2, embeds a supercell correlation treatment within a megacell and does not involve periodic image summation for the Coulomb integrals. Working in a basis of well-localized Wannier functions, periodicity is instead imposed through rigorous translational symmetry of Hamiltonian integrals and wavefunction parameters. The accuracy of the method is validated through comparison with a complementary periodic DLPNO-MP2 method that employs Born-von Kármán boundary conditions, described in Paper I of this series [Nejad et al., J. Chem. Phys. 163 (2025)]. The PNO approximations are shown to be equivalent in the two approaches and entirely consistent with molecular DLPNO-MP2 calculations. The Megacell-DLPNO-MP2 method displays sub-linear scaling with respect to supercell size at the asymptotic limit, and example calculations are presented with up to 15 000 basis functions in the correlation treatment.