Modern Electronic Structure Theory using Tensor Product States

Strongly correlated systems have been a major challenge for a long time in the field of theoretical chemistry. For such systems, the relevant portion of the Hilbert space scales exponentially, preventing efficient simulation on large systems. However, in many cases, the Hilbert space can be partitioned into clusters on the basis of strong and weak interactions. In this work, we mainly focus on an approach where we partition the system into smaller orbital clusters in which we can define many-particle cluster states and use traditional many-body methods to capture the rest of the inter-cluster correlations.

This dissertation can be mainly divided into two parts. In the first part of this dissertation, the clustered ansatz, termed as tensor product states (TPS), is used to study large strongly correlated systems. In the second part, we study a particular type of strongly correlated system, correlated triplet pair states that arise in singlet fission.

The many-body expansion (MBE) is an efficient tool that has a long history of use for calculating interaction energies, binding energies, lattice energies, and so on. We extend the incremental full configuration interaction originally proposed for a Slater determinant to a tensor product state (TPS) based wavefunction. By partitioning the active space into smaller orbital clusters, our approach starts from a cluster mean-field reference TPS configuration and includes the correlation contribution of the excited TPSs using a many-body expansion. This method, named cluster many-body expansion (cMBE), improves the convergence of MBE at lower orders compared to directly doing a block-based MBE from an RHF reference. The performance of the cMBE method is also tested on a graphene nano-sheet with a very large active space of 114 electrons in 114 orbitals, which would require 1066 determinants for the exact FCI solution.

Selected CI (SCI) using determinants becomes intractable for large systems with strong correlation. We introduce a method for SCI algorithms using tensor product states which exploits local molecular structure to significantly reduce the number of SCI variables. We demonstrate the potential of this method, called tensor product selected configuration interaction (TPSCI), using a few model Hamiltonians and molecular examples. These numerical results show that TPSCI can be used to significantly reduce the number of SCI variables in the variational space, and thus paving a path for extending these deterministic and variational SCI approaches to a wider range of physical systems.

The extension of the TPSCI algorithm for excited states is also investigated. TPSCI with perturbative corrections provides accurate excitation energies for low-lying triplet states with respect to extrapolated results. In the case of traditional SCI methods, accurate excitation energies are obtained only after extrapolating calculations with large variational dimensions compared to TPSCI. We provide an intuitive connection between lower triplet energy mani- folds with Hückel molecular orbital theory, providing a many-body version of Hückel theory for excited triplet states.

The n-body Tucker ansatz (which is a truncated TPS wavefunction) developed in our group provides a good approximation to the low-lying states of a clusterable spin system. In this approach, a Tucker decomposition is used to obtain local cluster states which can be truncated to prune the full Hilbert space of the system. As a truncated variational approach, it has been observed that the self-consistently optimized n-body Tucker method is not size- extensive, a property important for many-body methods. We explore the use of perturbation theory and linearized coupled-cluster methods to obtain a robust yet efficient approximation. Perturbative corrections to the n-body Tucker method have been implemented for the Heisenberg Hamiltonian and numerical data for various lattices and molecular systems has been presented to show the applicability of the method.

In the second part of this dissertation, we focus on studying a particular type of strongly correlated states that occurs in singlet fission material. The correlated triplet pair state 1(TT) is a key intermediate in the singlet fission process, and understanding the mechanism by which it separates into two independent triplet states is critical for leveraging singlet fission for improving solar cell efficiency. This separation mechanism is dominated by two key interactions: (i) the exchange interaction (K) between the triplets which leads to the spin splitting of the biexciton state into 1(TT),3(TT) and 5(TT) states, and (ii) the triplet-triplet energy transfer integral (t) which enables the formation of the spatially separated (but still spin entangled) state 1(T...T). We develop a simple ab initio technique to compute both the triplet-triplet exchange (K) and triplet-triplet energy transfer coupling (t). Our key findings reveal new conditions for successful correlated triplet pair state dissociation. The biexciton exchange interaction needs to be ferromagnetic or negligible compared to the triplet energy transfer for favorable dissociation. We also explore the effect of chromophore packing to reveal geometries where these conditions are achieved for tetracene.

We also provide a simple connectivity rule to predict whether the through-bond coupling will be stabilizing or destabilizing for the (TT) state in covalently linked singlet fission chromophores. By drawing an analogy between the chemical system and a simple spin-lattice, one is able to determine the ordering of the multi-exciton spin state via a generalized usage of Ovchinnikov's rule. In the case of meta connectivity, we predict 5(TT) to be formed and this is later confirmed by experimental techniques like time-resolved electron spin resonance (TR-ESR).