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.

<p>The vertical electronic excitation energies for the narrow-bandgap polymers LBPF, EP37 and EP62 have been calculated using Density Functional Theory (DFT). Also the vertical excitation energies for the acceptor unit of LBPF have been calculated using the Hartree-Fock (HF), DFT and Coupled Cluster (CC) methods. The calculations cover the visible and infrared wave length region and two strong transitions are obtained, one corresponding to the pi to pi* transition and one corresponding to the pi to Acceptor transition. The excitation energies obtained from DFT are below the corresponding experimental results and attempts have therefore been made to perform bench-marking calculations using a hierarchy of CC methods.</p>

The vertical electronic excitation energies for the narrow-bandgap polymers LBPF, EP37 and EP62 have been calculated using Density Functional Theory (DFT). Also the vertical excitation energies for the acceptor unit of LBPF have been calculated using the Hartree-Fock (HF), DFT and Coupled Cluster (CC) methods. The calculations cover the visible and infrared wave length region and two strong transitions are obtained, one corresponding to the pi to pi* transition and one corresponding to the pi to Acceptor transition. The excitation energies obtained from DFT are below the corresponding experimental results and attempts have therefore been made to perform bench-marking calculations using a hierarchy of CC methods.

Porphyrins are an important class of biomolecules, which are heavily studied, both ex- perimentally and computationally. But, despite the intensive efforts, for many questions we still aren't able to consistently find an agreement between theory and experiment. One of the still unresolved issues is the character of the ground state of the Fe(II)-porphyrin molecule. We used a model of the Fe(II)-porphyrin molecule to study the effects of geometrical changes on the spin states. By carrying out extensive DMRG-CASSCF cal- culations topped with TCCSD correlation treatment we are able to link the effects of these geometrical changes to the experimental results, and predict a quintet ground state for the isolated Fe(II)-porphyrin molecule. Also, using a ligated porphyrin belonging to the iron porphyrin carbene class of molecules, we demonstrate by combining the CASSCF and AC0 methods that geometrical changes outside the porphyrin core cannot be over- looked. 1