The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.
Our overarching goal is the development of wavefunction-based quantum chemistry methods that give good results for both weakly- and strongly-correlated systems, that preserve as many fundamental properties of the exact wavefunction as possible, and that are computationally affordable and black-box for the end user. Many of our contributions in this area are based on the archetypal quantum chemistry method, coupled-cluster theory, and they have resulted from our efforts to understand and remedy the failure of symmetry-adapted coupled-cluster theory in the presence of strong correlation. In this thesis, we recount these investigations and discuss the key insights and results that they have produced. We describe in detail the novel wavefunction approaches we have devised, as well as the approximations and considerations one must make in order to arrive at useful equations. We believe that this work could be important for the continued development of affordable methods capable of accurately describing a wide variety of quantum chemical systems.
The calculation of non-covalent interactions of highly polarizable molecules is an outstanding problem in chemistry and physics. Difficulties arise both in accurate treatments of dispersion interactions between molecules and in the evaluation of surfaces of sufficient size for use in further computations. While ``gold-standard" calculations using coupled cluster theory with singles, doubles and perturbative triples provides a reliable method for evaluating interactions the computational cost involved in calculating more than a few points of an interaction surface becomes prohibitive for even medium sized molecules. Expanding the long range interaction into a van der Waals series reduces the cost of evaluating a surface to the computation of a few parameters. We describe here the implementation of a new computer program for calculating van der Waals coefficients for arbitrary molecules using the sum over states method. The laboratory-frame transformation of the computed van der Waals surface and the inclusion of rotational state dressing on the surface is derived. Analytic approximations for the interactions of two linear molecules in the presence of a small DC electric field are also derived. Recent achievements in the formation and manipulation of ultracold polar molecules have opened the door to exciting new studies in cold chemical reactions. To characterize the energetics and reaction pathways of trimer and tetramer formation and reactions, we have computed the structure and thermochemistry of the model trimer Li
3 and every X
2Y
2 alkali tetramer through Cs. Related interest in the control of molecular ion reactions is also investigated with recent theoretical results for the rubidium hydroxide reaction presented here.
Perhaps the most important approximations to the electronic structure problem in quantum chemistry are those based on coupled cluster and density functional theories. Coupled cluster theory has been called the ``gold standard'' of quantum chemistry due to the high accuracy that it achieves for weakly correlated systems. Kohn--Sham density functionals based on semilocal approximations are, without a doubt, the most widely used methods in chemistry and material science because of their high accuracy/cost ratio. The root of the success of coupled cluster and density functionals is their ability to efficiently describe the dynamic part of the electron correlation. However, both traditional coupled cluster and density functional approximations may fail catastrophically when substantial static correlation is present. This severely limits the applicability of these methods to a plethora of important chemical and physical problems such as, e.g., the description of bond breaking, transition states, transition metal-, lanthanide- and actinide-containing compounds, and superconductivity. In an attempt to tackle this problem, nonstandard (single-reference) coupled cluster-based techniques that aim to describe static correlation have been recently developed: pair coupled cluster doubles (pCCD) and singlet-paired coupled cluster doubles (CCD0). The ability to describe static correlation in pCCD and CCD0 comes, however, at the expense of important amounts of dynamic correlation so that the high accuracy of standard coupled cluster becomes unattainable. Thus, the reliable and efficient description of static and dynamic correlation in a simultaneous manner remains an open problem for quantum chemistry and many-body theory in general. In this thesis, different ways to combine pCCD and CCD0 with density functionals in order to describe static and dynamic correlation simultaneously (and efficiently) are explored. The combination of wavefunction and density functional methods has a long history in quantum chemistry (practical implementations have appeared in the literature since the 1970s). However, this kind of techniques have not achieved widespread use due to problems such as double counting of correlation and the symmetry dilemma---the fact that wavefunction methods respect the symmetries of Hamiltonian, while modern functionals are designed to work with broken symmetry densities. Here, particular mathematical features of pCCD and CCD0 are exploited to avoid these problems in an efficient manner. The two resulting families of approximations, denoted as pCCD+DFT and CCD0+DFT, are shown to be able to describe static and dynamic correlation in standard benchmark calculations. Furthermore, it is also shown that CCD0+DFT lends itself to combination with correlation from the direct random phase approximation (dRPA). Inclusion of dRPA in the long-range via the technique of range-separation allows for the description of dispersion correlation, the remaining part of the correlation. Thus, when combined with the dRPA, CCD0+DFT can account for all three-types of electron correlation that are necessary to accurately describe molecular systems. Lastly, applications of CCD0+DFT to actinide chemistry are considered in this work. The accuracy of CCD0+DFT for predicting equilibrium geometries and vibrational frequencies of actinide molecules and ions is assessed and compared to that of well-established quantum chemical methods. For this purpose, the
f0 actinyl series (UO
22+, NpO
23+, PuO
24+), the isoelectronic NUN, and Thorium (ThO, ThO
2+) and Nobelium (NoO, NoO
2) oxides are studied. It is shown that the CCD0+DFT description of these species agrees with available experimental data and is comparable with the results given by the highest-level calculations that are possible for such heavy compounds while being, at least, an order of magnitude lower in computational cost.