We present the development of a real-space and projected congruent transformation method for treating electron correlation in chemical systems. This method uses an explicitly correlated function for performing congruent transformation on the electronic Hamiltonian. As a result of this transformation, the electronic Hamiltonian is transformed into a sum of two-, three-, four-, five-, and six-particle operators. Efficient computational implementation of these many-particle operators continues to be challenging for application of the congruent transformation approach for many-electron systems. In this work, we present a projected congruent transformed Hamiltonian (PCTH) approach to avoid computation of integrals involving operators that couple more than two particles. The projected congruent transformation becomes identical to the real-space congruent transformation in the limit of infinite basis size. However, for practical calculations, the projection is always performed on a finite-dimensional space. We show that after representing the contributing expressions of the PCTH in terms of diagrams, it is possible to identify a subset of diagrams that can be summed up to infinite order. This technique, denoted as partial infinite-order summation (PIOS), partly alleviates the limitation from the finite-basis representation of the PCTH method. The PCTH and PCTH-PIOS methods were applied to an isoelectronic series of 10-electron systems $(\mathrm{Ne},\mathrm{HF},{\mathrm{H}}_{2}\mathrm{O},{\mathrm{NH}}_{3},{\mathrm{CH}}_{4})$ and results were compared with configuration interaction (CISD) calculations. The results indicate that the PCTH-PIOS method can treat electron-electron correlations while avoiding explicit construction and diagonalization of the Hamiltonian matrix.
We present the development of a real-space and projected congruent transformation method for treating electron correlation in chemical systems. This method uses an explicitly correlated function for performing congruent transformation on the electronic Hamiltonian. As a result of this transformation, the electronic Hamiltonian is transformed into a sum of two-, three-, four-, five-, and six-particle operators. Efficient computational implementation of these many-particle operators continues to be challenging for application of the congruent transformation approach for many-electron systems. In this work, we present a projected congruent transformed Hamiltonian (PCTH) approach to avoid computation of integrals involving operators that couple more than two particles. The projected congruent transformation becomes identical to the real-space congruent transformation in the limit of infinite basis size. However, for practical calculations, the projection is always performed on a finite-dimensional space. We show that after representing the contributing expressions of the PCTH in terms of diagrams, it is possible to identify a subset of diagrams that can be summed up to infinite order. This technique, denoted as partial infinite-order summation (PIOS), partly alleviates the limitation from the finite-basis representation of the PCTH method. The PCTH and PCTH-PIOS methods were applied to an isoelectronic series of 10-electron systems (Ne,HF,H2O,NH3,CH4) and results were compared with configuration interaction (CISD) calculations. The results indicate that the PCTH-PIOS method can treat electron-electron correlations while avoiding explicit construction and diagonalization of the Hamiltonian matrix.
An implementation of the single-site density matrix renormalization group (DMRG) approach for the quantum-chemical application to calculate approximate as well as exact electronic ground-state energies and wave functions of atoms and molecules is presented.The algorithm is formulated in terms of the variational and sizeconsistent class of wave functions, matrix product states, and their operator analogues, matrix product operators.A procedure to construct quantum-chemical Hamiltonian operators with long-range Coulomb interactions as matrix product operators is developed.The many-electron Hamiltonian operators are represented as sums and products of one-electron operators.For the efficient matrix product operator construction of quantum-chemical Hamiltonians, different basis spaces including discrete position and momentum space as well as molecular and atomic orbitals are investigated.Discrete position space yields the best matrix product operator size scaling, but introduces artificial long-range couplings in the kinetic energy term of the Hamiltonian which are a consequence of the DMRG algorithm in terms of matrix product states.The huge amount of discrete space points needed for an accurate description poses a problem because of the electron-electron interaction.Therefore, a momentum space or Fourier series expansion of the electron-electron interaction is shown which makes it possible to compress its matrix product operator representation.DMRG applied to quantum-chemical Hamiltonians in matrix product operator form can serve as means to directly (without mean-field treatment in advance) investigate quantum-chemical problems and find the analytically exact ground-state solution of small atoms and molecules.The next basis space analyzed in detail is the atomic orbital basis.The matrix product operator construction and the direct application of the DMRG method in this basis is examined.A physical compression scheme based on a singular value decomposition of the electron-electron interaction in the atomic orbital basis is demonstrated.Furthermore, DMRG calculations in the atomic orbital basis for certain systems where mean-field theory yields qualitatively wrong results and where higher-order approximative correction methods like truncated configuration interaction and coupled cluster cannot fully capture strong electron correlation effects are compared to results obtained in the molecular orbital basis.The comparison shows that atomic orbitals describe the physical situation well and thus lead to a tremendous reduction of the matrix product state dimension and the computational effort.i
One of the major obstacles in modern quantum chemistry is the need to employ exceedingly large one-electron basis sets in electron correlation methods to obtain reliably accurate results. The family of explicitly correlated wave function methods circumvents this issue by introducing terms that explicitly depend on the interelectronic distance, thereby enforcing the correct electron–electron cusp behavior, albeit at the cost of additional complicated many-electron integrals when evaluated without further approximation. The main focus of the present work is the efficient, low-scaling evaluation of explicitly correlated approaches, examined here through the explicitly correlated correction to second-order Møller–Plesset perturbation theory (MP2-F12) as a representative example. While these methods are highly effective in reducing the basis set incompleteness error, the cost associated with the corresponding correction often far exceeds that of the underlying correlation method. To overcome these limitations, several techniques and strategies are introduced to significantly reduce computational demand and resolve the main bottlenecks that arise. A particularly straightforward approach is provided by the scaled MP2-F12 ansatz, which introduces a single empirically determined factor that scales the direct-type geminal–geminal terms to closely reproduce the MP2-F12 correction, thereby omitting the most expensive exchange-type contributions. Furthermore, the computation of Fock matrix elements spanning multiple orbital spaces is drastically accelerated by extending techniques originally developed for Hartree–Fock and density functional theory. The direct contributions are evaluated using a modified resolution-of-the-identity Coulomb (RI-J) approach employing the J-engine ansatz for integral construction, while the exchange contributions are treated via three-dimensional numerical quadrature in the form of the seminumerical sn-LinK method. These approaches reduce the computational cost by more than three orders of magnitude compared to the approximation-free evaluation while maintaining high accuracy, with additional speedups achieved by exploiting the massive parallelism of graphics processing units. The central focus of this work is the efficient decomposition of complicated exchange-type multi-electron integrals, where numerical quadrature plays a key role in improving performance. The complementary auxiliary basis set resolution of the identity (CABS-RI) reduces the most expensive exchange-type integral to a three-electron integral, which is then efficiently decomposed via numerical quadrature. In this context, a highly efficient batch-wise, distance-dependent integral screening exploiting the short-range behavior of F12 operators is introduced, significantly improving performance and enabling near-linear-scaling evaluations while retaining high accuracy. These strategies are extended to all exchange-type contributions in F12 theory, leading to an alternative formulation of integrals better suited for efficient factorization. A novel combination of numerical quadrature with density fitting is presented for the efficient computation of products of four-center two-electron integrals, alongside algorithms covering all exchange-type contributions, including advantageous batching strategies. In general, the formal computational cost with system size M is reduced from O(M⁵) to O(M⁴), with further improvements appearing feasible through the use of localized molecular orbitals combined with advanced integral screening techniques.
This thesis is devoted to the development of numerical methods for correlated electrons in realistic systems based on perturbation theories using Feynman diagrams. Electron correlations in realistic systems give rise to a wide range of interesting physical properties. Diagrammatic expansions and its low-order approximations are key building blocks of theoretical and numerical methods for studying correlated electrons. Recent developments of the diagrammatic quantum Monte Carlo (DiagMC) methods enable efficient stochastic evaluations of diagrams and does not suffer from numerical sign problem which worsens with increasing system size. It is therefore a promising direction to apply DiagMC techniques to realistic electron systems in the pursuit of ab initio calculations of correlated materials. This thesis starts by reviewing basic concepts and definitions of realistic electron systems, the perturbation theory based on Feynman diagrams, Monte Carlo evaluations of the series, as well as self-consistent diagrammatic methods. Low-order diagrammatic methods are applied to realistic calculations and benchmarked against large collections of state-of-the-art many-body methods. The thesis then proposes a diagrammatic Monte Carlo solver for realistic molecules and impurities based on a recursive evaluation of diagrams, which is tested on realistic molecules and impurities with varying temperature, system size, and interaction strength. We observe that this method is ideal for problems with many orbitals with moderate correlations, but may suffer from series divergence at strong correlations. Next, we propose the inchworm Monte Carlo algorithm for interaction expansion, which overcomes limitations of divergent diagrammatic series by replacing the standard perturbation series with a sequence of incremental series expansions. A prototype implementation is presented for system with on-site Hubbard interactions. Finally, we discuss technical aspects of numerical representations of the Green's function, which has significant impact on the efficiency of diagrammatic methods. We review several common numerical representations, and introduce our development of the sparse sampling method, which provides compact representations of the Green's function in both imaginary time and frequency as well as fast transformation formula between the two.
Closed-form expressions for all matrix elements required for variational calculation of the electronic structure of periodic solids have been derived using a basis of explicitly correlated Gaussians (ECGs). Periodic basis functions are constructed by summing shifted correlated Gaussians over all composite lattice translations, where a generalized unfolding theorem reduces the resulting double lattice sum to a single sum through a unified computational framework for overlap, kinetic energy, and Coulomb potential operators. The formalism has been validated through application to an infinite one-dimensional hydrogen chain, where the ground-state energy per atom computed in the thermodynamic limit is shown to agree with finite-chain results extrapolated by other many-body methods.
Quantum simulations of electronic structure with a transformed Hamiltonian that includes some electron correlation effects are demonstrated. The transcorrelated Hamiltonian used in this work is efficiently constructed classically, at polynomial cost, by an approximate similarity transformation with an explicitly correlated two-body unitary operator. This Hamiltonian is Hermitian, includes no more than two-particle interactions, and is free of electron-electron singularities. We investigate the effect of such a transformed Hamiltonian on the accuracy and computational cost of quantum simulations by focusing on a widely used solver for the Schrödinger equation, namely the variational quantum eigensolver method, based on the unitary coupled cluster with singles and doubles (q-UCCSD) Ansatz. Nevertheless, the formalism presented here translates straightforwardly to other quantum algorithms for chemistry. Our results demonstrate that a transcorrelated Hamiltonian, paired with extremely compact bases, produces explicitly correlated energies comparable to those from much larger bases. For the chemical species studied here, explicitly correlated energies based on an underlying 6-31G basis had cc-pVTZ quality. The use of the very compact transcorrelated Hamiltonian reduces the number of CNOT gates required to achieve cc-pVTZ quality by up to two orders of magnitude, and the number of qubits by a factor of three.