Despite the fact that higher-order tensors (HOTs) plague electronic structure methods and severely limits the modeling of interesting chemistry problems, introduction and application of higher-order tensor (HOT) decompositions, specifically the canonical polyadic (CP) decomposition, is fairly limited. The CP decomposition is an incredibly useful sparse tensor factorization that has the ability to disentangle all correlated modes of a tensor. However the complexities associated with CP decomposition have made its application in electronic structure methods difficult. Some of the major issues related to CP decomposition are a product of the mathematics of computing the decomposition: determining the exact CP rank is a non-polynomially hard problem, finding stationary points for rank-R approximations require non-linear optimization techniques, and inexact CP approximations can introduce a large degree of error into tensor networks. While other issues are a result of the construction of computer architectures. For example, computer processing units (CPUs) are organized in a way to maximize the efficiency of dense linear algebra and, thus, the performance of routine tensor algebra kernels, like the Khatri-Rao product, is limited. In this work, we seek to reduce the complexities associated with the CP decomposition and create a route for others to develop reduced-scaling electronic structure theory methods using the CP decomposition. In Chapter 2, we introduce the robust tensor network approximation. This approximation is a way to, in general, eliminate the leading-order error associated with approximated tensors in a network. We utilize the robust network approximation to significantly increase the accuracy of approximating density fitting (DF) integral tensors using rank-deficient CP decompositions in the particle-particle ladder (PPL) diagram of the coupled cluster method with single and double substitutions (CCSD). We show that one can produce results with negligible error in chemically relevant energy differences using a CP rank roughly the same size as the DF fitting basis; which is a significantly smaller rank requirement than found using either a nonrobust approximation or similar grid initialized CP approximations (the pseudospectral (PS) and tensor hypercontraction (THC) approximations). Introduction of the CP approximation, formally, reduces the complexity of the PPL diagram from 𝓞(N⁶) to 𝓞(N⁵) and, using the robust approximation, we are able to observe a cost reduction in CCSD calculations for systems as small as a single water molecule. In Chapter 3, we further demonstrate the utility of the robust network approximation and, in addition, we construct a scheme to optimize a grid-free CP decomposition of the order-four Coulomb integral tensor in 𝓞(N⁴) time. Using these ideas, we reduce the complexity of ten bottleneck contractions from 𝓞(N⁶) to 𝓞(N⁵) in the Laplace transform (LT) formulation of the perturbative triple, (T), correction to CCSD. We show that introducing CP into the LT (T) method with a CP rank roughly the size of the DF fitting basis reduces the cost of computing medium size molecules by a factor of about 2.5 and introduces negligible error into chemically relevant energy differences. Furthermore, we implement these low-cost algorithms using newly developed, optimized tensor algebra kernels in the massively-parallel, block-sparse TiledArray [Calvin, et. al Chemical Reviews 2021 121 (3), 1203-1231] tensor framework. / Doctor of Philosophy / Electronic structure methods and accurate modeling of quantum chemistry have developed alongside the advancements in computer infrastructures. Increasingly large and efficient computers have allowed researchers to model remarkably large chemical systems. Sadly, for as fast as computer infrastructures grow (Moores law predicts that the number of transistors in a computer will double every 18 months) the cost of electronic structure methods grows more quickly. One of the least expensive electronic structure methods, Hartree Fock (HF), grows quartically with molecular size; this means that doubling the size of a molecule increase the number of computer operations by a factor of 16. However, it is known that when chemical systems become sufficiently large, the amount of physical information added to the system grows linearly with system size.[Goedecker, et. al. Comput. Sci. Eng., 2003, 5, (4), 14-21] Unfortunately, standard implementations of electronic structure methods will never achieve linear scaling; the disparity between actual cost and physical scaling of molecules is a result of storing and manipulating data using dense tensors and is known as the curse of dimensionality.[Bellman, Adaptive Control Processes, 1961, 2045, 276] Electronic structure theorists, in their desire to apply accurate methods to increasingly large systems, have known for some time that the cost of conventional algorithms is unreasonably high. These theorists have found that one can reveal sparsity and develop reduced-complexity algorithms using matrix decomposition techniques. However, higher-order tensors (HOTs), tensors with more than two modes, are routinely necessary in algorithm formulations. Matrix decompositions applied to HOTs are not necessarily straight-forward and can have no effect on the limiting behavior of an algorithm. For example, because of the positive definiteness of the Coulomb integral tensor, it is possible to perform a Cholesky decomposition (CD) to reduce the complexity of tensor from an order-4 tensor to a product of order-3 tensors.[Beebe, et. al. Int. J. Quantum Chem., 1977, 12, 683-705] However, using the CD approximated Coulomb integral tensors it is not possible to reduce the complexity of popular methods such as Hartree-Fock or coupled cluster theory. We believe that the next step to reducing the complexity of electronic structure methods is through the accurate application of HOT decompositions. In this work, we only consider a single HOT decomposition: the canonical polyadic (CP) decomposition which represents a tensor as a polyadic sum of products. The CP decomposition disentangles all modes of a tensor by representing an order-N tensor as N order-2 tensors. In this work, we construct the CP decomposition of tensors using algebraic optimization. Our goal, here, is to tackle one of the biggest issues associated with the CP decomposition: accurately approximating tensors and tensor networks. In Chapter 2, we develop a robust formulation to approximate tensor networks, a formulation which removes the leading-order error associated with tensor approximations in a network.[Pierce, et. al. J. Chem. Theory Comput., 2021 17 (4), 2217- 2230] We apply a robust CP approximation to the coupled cluster method with single and double substitutions (CCSD) to reduce the overall cost of the approach. Using this robust CP approximation we can compute CCSD, on average, 2.5-3 times faster and introduce negligibly small error in chemically relevant energy values. Furthermore in Chapter 3, we again use the robust CP network approximation in conjunction with a novel, low cost approach to compute order-four CP decompositions, to reduce the cost of 10 high cost computations in the the perturbative triple, (T), correction to CCSD. By removing these computations, we are able to reduce the cost of (T) by a factor of about 2.5 while introducing significantly small error.

To follow up on the unexpectedly-good performance of coupled-cluster models with approx- imate inclusion of 3-body clusters [J. Chem. Phys. 151, 064102 (2019)] we performed a more complete assessment of the 3CC method [J. Chem. Phys. 125, 204105 (2006)] for accurate computational thermochemistry in the standard HEAT framework. New spin- integrated implementation of the 3CC method applicable to closed- and open-shell systems utilizes a new automated toolchain for derivation, optimization, and evaluation of operator algebra in many-body electronic structure. We found that with a double-zeta basis set the 3CC correlation energies and their atomization energy contributions are almost always more accurate (with respect to the CCSDTQ reference) than the CCSDT model as well as the standard CCSD(T) model. The mean errors in { 3CC, CCSDT, and CCSD(T) } electronic (per valence electron) and atomization energies were {23, 69, 125} μEh/e and {0.39, 1.92, 2.57} kJ/mol, respectively. The significant and systematic reduction of the error by the 3CC method and its lower cost than CCSDT suggests it as a viable candidate for post-CCSD(T) thermochemistry application. / Doctor of Philosophy / Stepping into the information age, the computing power has rapidly grown over the last half century. Solving chemical problems on computers has improved lives by reducing the cost and time of researching critical technologies. Scientific research is evolving and experimental finding are now supported with a computational model. Doing chemistry on computers requires quantum simulations, which is essentially solving the Schr ̈odinger equation on a computer that simulates a wave function for all the electrons in a system. Different models are built based on how these inter electronic interactions are treated. To predict results with accuracy on par with the experimental findings requires using higher-order wave functions methods.These are computationally expensive and often not practical. The lower-order methods that are easy to implement can be found in all quantum chemistry software packages. On the other hand, the higher-order methods are laborious and error prone to implement manually due to the sheer complexity of theory. Debugging such implementations often requires a lot of effort with the uncertainty in returns. To solve this problem, we implemented a second-quantization toolkit (SeQuant version 2.0) that derives many-body methods, specifically the general-order coupled cluster (CC) model. The CC model is systematically improvable and accurate. One such CC model, the CCSD(T), has been called the gold standard in quantum chemistry. For compactness, these equations are usually derived in their spin-orbital form. The evaluation and storage cost of these methods is reduced up to four-fold by transforming the spin-orbital expressions to a spin-traced form. In this work, the spin-tracing algorithms are described in detail. The general-order coupled cluster approach is used to derive the internally corrected approximate coupled cluster methods. These methods improve the accuracy of a model at a reduced cost. For small molecules, it was observed that the spin-traced evaluation was over three times faster than spin-orbital coupled cluster. To further reduce the cost of calculations, we added explicit correlation to our CC models. These methods improved the quality of our results with a modest increase in the computational cost.

Accurate correlated electronic structure methods involve a significant amount of computations and can be only employed to small molecular systems. For example, the coupled-cluster singles, doubles, and perturbative triples model (CCSD(T)), which is known as the ``gold standard" of quantum chemistry for its accuracy, usually can treat molecules with 20-30 atoms. To extend the reach of accurate correlated electronic structure methods to larger molecular systems, we work towards two directions: parallel computing and reduced-cost/scaling approaches. Parallel computing can utilize more computational resources to handle systems that demand more substantial computational efforts. Reduced-cost/scaling approaches, which introduce approximations to the existing electronic structure methods, can significantly reduce the amount of computation and storage requirements. In this work, we introduce a new distributed-memory massively parallel implementation of standard and explicitly correlated (F12) coupled-cluster singles and doubles (CCSD) with canonical bigO{N^6} computational complexity ( C. Peng, J. A. Calvin, F. Pavov{s}evi'c, J. Zhang, and E. F. Valeev, textit{J. Phys. Chem. A} 2016, textbf{120}, 10231.), based on the TiledArray tensor framework. Excellent strong scaling is demonstrated on a multi-core shared-memory computer, a commodity distributed-memory computer, and a national-scale supercomputer. We also present a distributed-memory implementation of the density-fitting (DF) based CCSD(T) method. (C. Peng, J. A. Calvin, and E. F. Valeev, textit{in preparation for submission}) An improved parallel DF-CCSD is presented utilizing lazy evaluation for tensors with more than two unoccupied indices, which makes the DF-CCSD storage requirements always smaller than those of the non-iterative triples correction (T). Excellent strong scaling is observed on both shared-memory and distributed-memory computers equipped with conventional Intel Xeon processors and the Intel Xeon Phi (Knights Landing) processors. With the new implementation, the CCSD(T) energies can be evaluated for systems containing 200 electrons and 1000 basis functions in a few days using a small size commodity cluster, with even more massive computations possible on leadership-class computing resources. The inclusion of F12 correction to the CCSD(T) method makes it converge to basis set limit much more rapidly. The large-scale parallel explicitly correlated coupled-cluster program makes the accurate estimation of the coupled-cluster basis set limit for molecules with 20 or more atoms a routine. Thus, it can be used rigorously to test the emerging reduced-scaling coupled-cluster approaches. Moreover, we extend the pair natural orbital (PNO) approach to excited states through the equation-of-motion coupled cluster singles and doubles (EOM-CCSD) method. (C. Peng, M. C. Clement, and E. F. Valeev, textit{submitted}) We simulate the PNO-EOM-CCSD method using an existing massively parallel canonical EOM-CCSD program. We propose the use of state-averaged PNOs, which are generated from the average of the pair density of excited states, to span the PNO space of all the excited states. The doubles amplitudes in the CIS(D) method are used to compute the state-averaged pair density of excited states. The issue of incorrect states in the state-averaged pair density, caused by an energy reordering of excited states between the CIS(D) and EOM-CCSD, is resolved by simply computing more states than desired. We find that with a truncation threshold of $10^{-7}$, the truncation error for the excitation energy is already below 0.02 eV for the systems tested, while the average number of PNOs is reduced to 50-70 per pair. The accuracy of the PNO-EOM-CCSD method on local, Rydberg and charge transfer states is also investigated. / Ph. D.

In dieser Arbeit werden Algorithmen für die Berechnung elektronischer Korrelations- und Anregungsenergien mittels der Coupled-Cluster Methode auf adaptiven Gittern entwickelt und implementiert. Die jeweiligen Funktionen und Operatoren werden adaptiv durch Multiskalenanalyse dargestellt, was eine Basissatz unabängige Beschreibung mit kontrollierter numerischer Genauigkeit ermöglicht. Gleichungen für die Coupled-Cluster Methode werden in einem verallgemeinerten Rahmen, unabhängig von virtuellen Orbitalen und globalen Basissätzen, neu formuliert. Hierzu werden die amplitudengewichteten Anregungen in virtuelle Orbitale ersetzt durch Anregungen in n-Elektronenfunktionen, welche durch Gleichungen im n-Elektronen Ortsraum bestimmt sind. Die erhaltenen Gleichungen können, analog zur Basissatz abh¨angigen Form, mit leicht angepasster Interpretation diagrammatisch dargestellt werden. Aufgrund des singulären Coulomb Potentials werden die Arbeitsgleichungen mit einem explizit korrelierten Ansatz regularisiert. Coupled-Cluster singles mit genäherten doubles (CC2) und ähnliche Modelle werden, für geschlossenschalige Systeme und in regularisierter Form, in die MADNESS Bibliothek (eine allgemeine Bibliothek zur Darstellung von Funktionen und Operatoren mittels Multiskalenanalyse) implementiert. Mit der vorgestellten Methode können elektronische CC2 Paarkorrelationsenergien und Anregungsenergien mit bestimmter numerischer Genauigkeit unabhängig von globalen Basissätzen berechnet werden, was anhand von kleinen Molekülen verifiziert wird / In this work algorithms for the computation of electronic correlation and excitation energies with the Coupled-Cluster method on adaptive grids are developed and implemented. The corresponding functions and operators are adaptively represented with multiresolution analysis allowing a basis-set independent description with controlled numerical accuracy. Equations for the coupled-cluster model are reformulated in a generalized framework independent of virtual orbitals and global basis-sets. For this, the amplitude weighted excitations into virtuals are replaced by excitations into n-electron functions which are determined by projected equations in the n-electron position space. The resulting equations can be represented diagrammatically analogous to basis-set dependent approaches with slightly adjusted rules of interpretation. Due to the singular Coulomb potential, the working equations are regularized with an explicitly correlated ansatz. Coupled-cluster singles with approximate doubles (CC2) and similar models are implemented for closed-shell systems and in regularized form into the MADNESS library (a general library for the representation of functions and operators with multiresolution analysis). With the presented approach electronic CC2 pair-correlation energies and excitation energies can be computed with definite numerical accuracy and without dependence on global basis sets, which is verified on small molecules.

Many-Body Perturbation Theory is put to test as a method for reliable calculations of the electron-electron interaction in two-dimensional quantum dots. We show that second order correlation gives qualitative agreement with experiments on a level which was not found within the Hartree-Fock description. For weaker confinements, the second order correction is shown to be insufficient and higher order contributions must be taken into account. We demonstrate that all order Many-Body Perturbation Theory in the form of the Coupled Cluster Singles and Doubles method yields very reliable results for confinements close to those estimated from experimental data. The possibility to use very large basis sets is shown to be a major advantage compared to Full Configuration Interaction approaches, especially for more than five confined electrons. Also, the possibility to utilize two-electron correlation in combination with tailor made potentials to achieve useful properties is explored. In the case of a two-dimensional quantum dot molecule we vary the interdot distance, and in the case of a two-dimensional quantum ring we vary the ring radius, in order to alter the spectra. In the latter case we demonstrate that correlation in combination with electromagnetic pulses can be used for the realization of quantum logical gates. / At the time of the doctoral defense, the following papers were unpublished and had a status as follows: Paper 5: Manuscript.

Extensive theoretical studies on the ground and excited state properties of systems containing heavy atoms have shown that accurate prediction of transition energies and related properties requires the incorporation of both relativistic and higher order correlation and relaxation effects as these effects are strongly inter- wined. The relativistic and dynamical electron correlation effects can be incor- porated in many-electron systems through a variety of many-body methods like configuration interaction (CI), coupled cluster method (CCM) etc. which are very powerful and effective tool for high precision description of electron correlation in many-electron systems. In this thesis, we investigate the relativistic and correlation effects in heavy atomic and molecular systems using these two highly correlated many-body methods. It is well recognized that, heavy polar diatomic molecules such as BaF, YbF, TlF, PbO, etc. are the leading experimental candidates for the search of violation of Parity (P ) and Time-reversal (T ) symmetry. The experimental detection of such P,T-odd effects in atoms and molecules has important consequences for the theory of fundamental interactions or for physics beyond the standard model (SM). For instance, a series of experiments on TlF have already been reported which provide the tightest limit available on the tensor coupling constant C , proton electric dipole moment (EDM) dp , etc. Experiments on YbF and BaF molecules are also of fundamental significance to the study of symmetry violation in nature, as these experiments have the potential to detect effects due to the electron EDMde. It is therefore imperative that high precession calculations are necessary to interpret these ongoing (and perhaps forthcoming) experimental outcome. For example, the knowledge of the effective electric field E(characterized by Wd) at the unpaired electron is required to link the experimentally determined P,T-odd frequency shift with the electron EDM de. We begin with a brief review of P,T-odd effects in heavy atoms and polar diatomics and the possible mechanisms which can give rise to such effects, in particular, the one arises due to the intrinsic electron EDM de. The P,T-odd interaction constant Wd is computed for the ground (2∑ ) state of YbF and BaF molecules using all-electron DF orbitals at the restricted active space (RAS) CI level. The RASCI space used for both systems in this calculation is sufficiently large to incorporate important core-core, core-valence, and valence-valence electron correlation effects. In addition to Wd, we also report the dipole moment (µe ) for these systems to assess the reliability of the method. The basis set dependency of Wd is also analyzed. The single reference coupled cluster (SRCC) method, developed by the cluster expansion of a single determinant reference function, is one of the most sophisticated method for treating dynamical correlation effects in a size-extensive manner. The non-uniqueness of the exponential nature of the wave operator diversifies the methods in multi-reference context. The multi-reference coupled cluster (MRCC) strategies fall within two broad classes: (a) State-Universal (SU), a Hilbert-space approach and (b) Valence-Universal (VU), a Fock-space approach. In this thesis, we shall be mainly concerned with the VU-MRCC which unlike SU-MRCC uses a single wave operator that not only correlates the reference functions, but also all the lower valence (or the so called subdued) sectors, obtained by deleting the occupancies systematically. The linear response theory (LRT) or equation of motion (EOM) method is another possible alternative which is nowadays extensively used to compute the atomic and molecular properties. Although, the CCLRT or EOM-CC method is not fully extensive in nature, this method has some distinct advantages over the traditional VU-MRCC theory. Further, for one-valence problem like ionization processes, the CCLRT/EOM-CC is formally equivalent to VU-MRCC, and hence, size-extensive. In this thesis, the core-extensive CCLRT and core-valence extensive (all electron) VU-MRCC methods are applied to compute the ground and excited state properties of various atomic and molecular systems (HCl, CuH, Ag, Sr, Yb and Hg) using nonrelativistic and relativistic (for heavy atoms) spinors. The similarities and differences in the structure of these two formalisms are also addressed. We also investigate the ground and excited state properties of HCN which is a system of astrophysical importance. This system has raised interest among the astrophysicists due to its detection in the atmosphere of Titan and Carbon stars. HCN has also been identified via radio-techniques in both comets and interstellar atmosphere. In the ash-photolysis of oxazole, iso-oxazole, and thiozole a transient band system was observed in the region 2500-3050 Å. This band system was attributed to a meta-stable form of HCN, i.e, either HNC or triplet HCN. We carry out detailed theoretical investigations using CCLRT and complete active space self-consistent field (CASSCF) method to characterize this unidentified band and other experimentally observed transitions.

In dieser Arbeit wird eine unrestricted Coupled-Cluster CC2 Response-Methode für die Berechnung von Eigenschaften erster und zweiter Ordnung, mit dem elektronischen g-Tensor als Schwerpunkt, präsentiert. Lokale Korrelations- und Dichtefittingnäherungen wurden verwendet. Die fundamentalen Konzepte notwendig für das Verständnis von Coupled-Cluster-Theorie, Dichtefitting, lokaler Korrelation, allgemeinen Coupled-Cluster Eigenschaften und dem elektronischen g-Tensor werden diskutiert. Die berechneten g-Tensoren werden mit denen durch Coupled-Cluster Singles and Doubles, Dichtefunktionaltheorie und Experiment erhaltenen verglichen. Effizienz und Genauigkeit der Näherung wird untersucht. Ein detailierter Anhang beschreibt die diagrammatische Coupled-Cluster-Theorie sowie ihre Anwendung zur Herleitung der verwendeten Arbeitsgleichungen. Die in dieser Arbeit entwickelte Methode ermöglicht es, den elektronischen g-Tensor von ausgedehnten Systemen mit einer Methode, die nicht auf Dichtefunktionaltheorie basiert, quantitativ vorherzusagen. Damit ist sie ein wichtiger Schritt hin zur Entwicklung von niedrig skalierenden Coupled-Cluster-Methoden höherer Ordnung für diese Art von Problem. / This work presents an unrestricted coupled-cluster CC2 response method using local correlation and density fitting approximations for the calculation of first and second order properties with particular focus on the electronic g-tensor. The fundamental concepts related to coupled-cluster theory, density fitting, local correlation, general coupled-cluster properties and the electronic g-tensor are discussed. The calculated g-tensors are benchmarked against those obtained from coupled-cluster singles and doubles, density functional theory and experiment. Efficiency and accuracy of the approximations is investigated. A detailed appendix covers the fundamentals of diagrammatic coupled-cluster and its application to the derivation of the working equations. The method presented in this thesis enables the quantitative prediction of the electronic g-tensor of extended systems with a method other than density functional theory. It represents an important step towards the development of low-scaling higher order coupled-cluster methods for this type of problem.

No description available.

Many-body perturbation theory and coupled-cluster theory, combined with carefully constructed basis sets, can be used to accurately compute the properties of small molecules. We applied a series of methods and basis sets aimed at reaching the ab initio limit to determine the barrier to planarity for ethylene cation. For potential energy surfaces corresponding to bond dissociation, a single Slater determinant is no longer an appropriate reference, and the single-reference hierarchy breaks down. We computed full configuration interaction benchmark data for calibrating new and existing quantum chemical methods for the accurate description of potential energy surfaces. We used the data to calibrate single-reference configuration interaction, perturbation theory, and coupled-cluster theory and multi-reference configuration interaction and perturbation theory, using various types of molecular orbitals, for breaking single and multiple bonds on ground-state and excited-state surfaces. We developed a determinant-based method which generalizes the formulation of many-body wave functions and energy expectation values. We used the method to calibrate single-reference and multi-reference configuration interaction and coupled-cluster theories, using different types of molecular orbitals, for the symmetric dissociation of water. We extended the determinant-based method to work with general configuration lists, enabling us to study, for the first time, arbitrarily truncated coupled-cluster wave functions. We used this new capability to study the importance of configurations in configuration interaction and coupled-cluster wave functions at different regions of a potential energy surface.

The work presented in this thesis is dedicated to developing inexpensive quantum-chemical models that are able to produce smooth and physically correct potential energy curves for the dissociation of single covalent bonds. It is well known that the energies produced by many ab initio theories scaling as the fifth order with the system size (for instance, second-order Moller-Plesset (MP2) and Epstein-Nesbet perturbation theories) diverge at large interatomic separations. We show that the divergent behavior of such perturbation schemes is due to a small number of terms in the energy expressions. Then, we demonstrate that the self-consistent replacement of these terms by their analogs from the coupled cluster theory (such as CCSD) allows one to redress the erroneous behavior of the perturbation theories without the damage to the overall scaling. We also investigate the accuracy of these hybrid perturbation theory-coupled cluster theories near equilibrium geometry. Judging from the computed spectroscopic constants and shapes of the potential energy curves, one such model, denoted MP2-CCSD(II) in this work, performs consistently better than the MP2 theory at essentially the same computational cost.