Computing many-body wave functions with guaranteed precision: The first-order Møller-Plesset wave function for the ground state of helium atom

Scalar Relativistic Effects with Multiwavelets: Implementation and Benchmark

The importance of relativistic effects in quantum chemistry is widely recognized, not only for heavier elements but throughout the periodic table. At the same time, relativistic effects are strongest in the nuclear region, where the description of electrons through a linear combination of atomic orbitals becomes more challenging. Furthermore, the choice of basis sets for heavier elements is limited compared with lighter elements where precise basis sets are available. Thanks to the framework of multiresolution analysis, multiwavelets provide an appealing alternative to overcoming this challenge: they lead to robust error control and adaptive algorithms that automatically refine the basis set description until the desired precision is reached. This allows one to achieve a proper description of the nuclear region. In this work, we extended the multiwavelet-based code MRChem to the scalar zero-order regular approximation framework. We validated our implementation by comparing the total energies for a small set of elements and molecules. To confirm the validity of our implementation, we compared both against a radial numerical code for atoms and the plane-wave-based code EXCITING.

Direct Determination of Optimal Real-Space Orbitals for Correlated Electronic Structure of Molecules

We demonstrate how to determine numerically nearly exact orthonormal orbitals that are optimal for the evaluation of the energy of arbitrary (correlated) states of atoms and molecules by minimization of the energy Lagrangian. Orbitals are expressed in real space using a multiresolution spectral element basis that is refined adaptively to achieve the user-specified target precision while avoiding the ill-conditioning issues that plague AO basis set expansions traditionally used for correlated models of molecular electronic structure. For light atoms, the orbital solver, in conjunction with a variational electronic structure model [selected Configuration Interaction (CI)] provides energies of comparable precision to a state-of-the-art atomic CI solver. The computed electronic energies of atoms and molecules are significantly more accurate than the counterparts obtained with the orbital sets of the same rank expanded in Gaussian AO bases, and can be determined even when linear dependence issues preclude the use of the AO bases. It is feasible to optimize more than 100 fully correlated numerical orbitals on a single computer node, and significant room exists for additional improvement. These findings suggest that real-space orbital representations might be the preferred alternative to AO representations for high-end models of correlated electronic states of molecules and materials.

Reducing Qubit Requirements while Maintaining Numerical Precision for the Variational Quantum Eigensolver: A Basis-Set-Free Approach

We present a basis-set-free approach to the variational quantum eigensolver using an adaptive representation of the spatial part of molecular wave functions. Our approach directly determines system-specific representations of qubit Hamiltonians while fully omitting globally defined basis sets. In this work, we use directly determined pair-natural orbitals on the level of second-order perturbation theory. This results in compact qubit Hamiltonians with high numerical accuracy. We demonstrate initial applications with compact Hamiltonians on up to 22 qubits where conventional representation would for the same systems require 40-100 or more qubits. We further demonstrate reductions in the quantum circuits through the structure of the pair-natural orbitals.

Grid-based diffusion Monte Carlo for fermions without the fixed-node approximation

A diffusion Monte Carlo algorithm is introduced that can determine the correct nodal structure of the wave function of a few-fermion system and its ground-state energy without an uncontrolled bias. This is achieved by confining signed random walkers to the points of a uniform infinite spatial grid, allowing them to meet and annihilate one another to establish the nodal structure without the fixed-node approximation. An imaginary-time propagator is derived rigorously from a discretized Hamiltonian, governing a non-Gaussian, sign-flipping, branching, and mutually annihilating random walk of particles. The accuracy of the resulting stochastic representations of a fermion wave function is limited only by the grid and imaginary-time resolutions and can be improved in a controlled manner. The method is tested for a series of model problems including fermions in a harmonic trap as well as the He atom in its singlet or triplet ground state. For the latter case, the energies approach from above with increasing grid resolution and converge within 0.015E_{h} of the exact basis-set-limit value for the grid spacing of 0.08 a.u. with a statistical uncertainty of 10^{-5}E_{h} without an importance sampling or Jastrow factor.

Combining Pseudopotential and All Electron Density Functional Theory for the Efficient Calculation of Core Spectra Using a Multiresolution Approach

Broadly speaking, the calculation of core spectra such as electron energy loss spectra (EELS) at the level of density functional theory (DFT) usually relies on one of two approaches: conceptually more complex but computationally efficient projector augmented wave based approaches or more straightforward but computationally more intensive all electron (AE) based approaches. In this work we present an alternative method, which aims to find a middle ground between the two. Specifically, we have implemented an approach in the multiwavelet madness molecular DFT code that permits a combination of atoms treated at the AE and pseudopotential (PSP) level. Atoms for which one wishes to calculate the core edges are thus treated at an AE level, while the remainder can be treated at the PSP level. This is made possible thanks to the multiresolution approach of madness, which permits accurate and efficient calculations at both the AE and PSP level. Through examples of a small molecule and a carbon nanotube, we demonstrate the potential applications of our approach.

Optimized Pair Natural Orbitals for the Coupled Cluster Methods

We present the coupled-cluster singles and doubles method formulated in terms of truncated pair natural orbitals (PNO) that are optimized to minimize the effect of truncation. Compared to the standard ground-state PNO coupled-cluster approaches, in which truncated PNOs derived from first-order Møller-Plesset (MP1) amplitudes are used to compress the CC wave operator, the iteratively optimized PNOs ("iPNOs") offer moderate improvement for small PNO ranks but rapidly increase their effectiveness for large PNO ranks. The error introduced by PNO truncation in the CCSD energy is reduced by orders of magnitude in the asymptotic regime, with an insignificant increase in PNO ranks. The effect of PNO optimization is particularly effective when combined with Neese's perturbative correction for the PNO incompleteness of the CCSD energy. The use of the perturbative correction in combination with the PNO optimization procedure seems to produce the most precise approximation to the canonical CCSD energies for small and large PNO ranks. For the standard benchmark set of noncovalent binding energies, remarkable improvements with respect to the standard PNO approach range from a factor of 3 with PNO truncation threshold τ_PNO = 10^-6 (with the maximum PNO truncation error in the binding energy of only 0.1 kcal/mol) to more than 2 orders of magnitude with τ_PNO = 10^-9.

Lowering of the complexity of quantum chemistry methods by choice of representation

The complexity of the standard hierarchy of quantum chemistry methods is not invariant to the choice of representation. This work explores how the scaling of common quantum chemistry methods can be reduced using real-space, momentum-space, and time-dependent intermediate representations without introducing approximations. We find the scalings of exact Gaussian basis Hartree-Fock theory, second-order Møller-Plesset perturbation theory, and coupled cluster theory (specifically, linearized coupled cluster doubles and the distinguishable cluster approximation with doubles) to be O(N³), O(N³), and O(N⁵), respectively, where N denotes the system size. These scalings are not asymptotic and hold over all ranges of N.

Magnetic properties with multiwavelets and DFT: the complete basis set limit achieved

Multiwavelets are emerging as an attractive alternative to traditional basis sets such as Gaussian-type orbitals and plane waves. One of their distinctive properties is the ability to reach the basis set limit (often a chimera for traditional approaches) reliably and consistently by fixing the desired precision ε. We present our multiwavelet implementation of the linear response formalism, applied to static magnetic properties, at the self-consistent field level of theory (both for Hartree-Fock and density functional theories). We demonstrate that the multiwavelets consistently improve the accuracy of the results when increasing the desired precision, yielding results that have four to five digits precision, thus providing a very useful benchmark which could otherwise only be estimated by extrapolation methods. Our results show that magnetizabilities obtained with the augmented quadruple-ζ basis (aug-cc-pCVQZ) are practically at the basis set limit, whereas absolute nuclear magnetic resonance shielding tensors are more challenging: even by making use of a standard extrapolation method, the accuracy is not substantially improved. In contrast, our results provide a benchmark that: (1) confirms the validity of the extrapolation ansatz; (2) can be used as a reference to achieve a property-specific extrapolation scheme, thus providing a means to obtain much better extrapolated results; (3) allows us to separate functional-specific errors from basis-set ones and thus to assess the level of cancellation between basis set and functional errors often exploited in density functional theory.

Coupled-Cluster in Real Space. 1. CC2 Ground State Energies Using Multiresolution Analysis

A framework to calculate CC2 approximated coupled-cluster ground state correlation energies in a multiresolution basis is derived and implemented into the MADNESS library. The CC2 working equations are formulated in first quantization which makes them suitable for real-space methods. The first quantized equations can be interpreted diagrammatically using the usual diagrams from second quantization with adjusted interpretation rules. Singularities arising from the nuclear and electronic potentials are regularized by explicitly taking the nuclear and electronic cusps into account. The regularized three- and six-dimensional cluster functions are represented directly on an adaptive grid. The resulting equations are free of singularities and virtual orbitals, which results in a low intrinsic scaling. Correlation energies close to the basis set limit are computed for small molecules. This work is the first step toward CC2 excitation energies in a multiresolution basis.

Analytic second nuclear derivatives of Hartree-Fock and DFT using multi-resolution analysis

We present the formalism, implementation, and numerical results for the computation of second derivatives with respect to nuclear displacements of molecules in the formalism of multi-resolution analysis. The highly singular nuclear potentials are partially regularized to improve the numerical stability. Vibrational frequencies are well reproduced to within an RMS of a few cm^-1 compared to large basis set LCAO (linear combination of atomic orbitals) calculations. Intermolecular modes, hindered rotations, and heavy atoms may lead to loss of precision. Tight precision thresholds are therefore necessary to converge to numerically stable results.

Tensor numerical methods in quantum chemistry: from Hartree–Fock to excitation energies

We resume the recent successes of the grid-based tensor numerical methods and discuss their prospects in real-space electronic structure calculations.

Multiresolution quantum chemistry in multiwavelet bases: excited states from time-dependent Hartree–Fock and density functional theory via linear response

A fully numerical method for the time-dependent Hartree–Fock and density functional theory (TD-HF/DFT) with the Tamm–Dancoff (TD) approximation is presented in a multiresolution analysis (MRA) approach.

The grid-based fast multipole method – a massively parallel numerical scheme for calculating two-electron interaction energies

A grid-based fast multipole method has been developed for calculating two-electron interaction energies for non-overlapping charge densities.

Real-space numerical grid methods in quantum chemistry

This themed issue reports on recent progress in the fast developing field of real-space numerical grid methods in quantum chemistry.

Adaptive multiconfigurational wave functions

A method is suggested to build simple multiconfigurational wave functions specified uniquely by an energy cutoff Λ. These are constructed from a model space containing determinants with energy relative to that of the most stable determinant no greater than Λ. The resulting Λ-CI wave function is adaptive, being able to represent both single-reference and multireference electronic states. We also consider a more compact wave function parameterization (Λ+SD-CI), which is based on a small Λ-CI reference and adds a selection of all the singly and doubly excited determinants generated from it. We report two heuristic algorithms to build Λ-CI wave functions. The first is based on an approximate prescreening of the full configuration interaction space, while the second performs a breadth-first search coupled with pruning. The Λ-CI and Λ+SD-CI approaches are used to compute the dissociation curve of N2 and the potential energy curves for the first three singlet states of C2. Special attention is paid to the issue of energy discontinuities caused by changes in the size of the Λ-CI wave function along the potential energy curve. This problem is shown to be solvable by smoothing the matrix elements of the Hamiltonian. Our last example, involving the Cu2O2(2+) core, illustrates an alternative use of the Λ-CI method: as a tool to both estimate the multireference character of a wave function and to create a compact model space to be used in subsequent high-level multireference coupled cluster computations.