Despite an immense progress in recent decades, a precise treatment of strongly corre- lated molecular systems still remains a challenge as of today. To help solve this problem, we have developed a massively parallel implementation of DMRG, called MOLMPS. For nonrelativistic systems requiring accurate treatment of both static and dynamical corre- lation, we have extended MOLMPS by the means of the almost linear scaling DLPNO- TCCSD method. In relativistic domain, we are the first who employed the 4c-CCSD to add dynamical correlation on top of DMRG, yielding the 4c-TCCSD method. When tested on benchmarks like the π-conjugated anthracene tetramer with CAS(63,63) and the FeMoco cofactor with CAS(113,76). We showed a good parallel performance on up to about 2000 CPUs. On the example of Iron(II)-Porphyrin model, we showed that the DLPNO-TCCSD captures 99.9% of TCCSD correlation energy. Our spectroscopic study on heavy diatomics showed that the 4c-TCCSD approach increases the precision of underlying CCSD to the order of CCSD(T) and that it is a promising approach. The thesis discusses three different implementations of quantum chemical methods based on QC-DMRG. 1
The central problem in the modern electronic structure theory is the calculation of cor- relation energy, possibly by an approach that would account for both static and dynamic correlation in an efficient, balanced and accurate way. In this thesis, I present a collection of methods that combine the effective treatment of dynamic correlation by the coupled cluster theory with density matrix renormalization group, a well-established technique for calculations of strongly correlated systems. The connection between them is achieved via the tailored coupled clusters (TCC) ansatz, which conveniently does not impose any ad- ditional computational costs. After the successful initial assessment, we developed more efficient implementations of these methods by employing the local approaches based on pair natural orbitals. This way, we extended the range of possible applications to larger systems with thousands of basis functions. To assess the accuracy of TCC as well as its local counterparts, we performed a variety of benchmark calculations ranging from small, yet challenging systems such as the nitrogen molecule or tetramethyleneethane diradical, to larger molecules like oxo-Mn(Salen) or Fe(II)-porphyrin model. 1
Recent developments in quantum embedding theories have provided attractive approaches to correlated calculations for large systems. In this work, we extend our previous work [ J. Chem. Theory Comput. 2019, 15, 4497–4506; J. Phys. Chem. Lett. 2019, 10, 6368–6374] on bootstrap embedding (BE) to enable correlated ab initio calculations at the coupled cluster with singles and doubles (CCSD) level for large molecules. We introduce several new algorithmic developments that significantly reduce the computational cost of BE, while maintaining its accuracy. The resulting implementation scales as O ( N 3 ) for the integral transform and O ( N ) for the CCSD calculation. Numerical results on a series of conjugated molecules suggest that BE with reasonably sized fragments can recover more than 99.5% of the total correlation energy of a full CCSD calculation, while the required computational resources (time and storage) compare favorably to one popular local correlation scheme: domain localized pair natural orbital (DLPNO). The largest BE calculation in this work involves ∼2900 basis functions and can be performed on a single node with 16 CPU cores and 64 GB of memory in a few days. We anticipate that these developments represent an important step toward the application of BE to solve practical problems.
Computational chemistry has made remarkable progress in the last couple of decades due to the availability of better and more powerful computers. However, there is definitely more room for improvement. There are problems in material science and biology that require 102 to 104 atoms to be considered and the importance of such systems cannot be over-emphasized. Addressing such systems with ab-initio methods is still a quantum chemist’s dream. ¨ Brute force solutions of the Schrodinger equation are not possible for reasonably interesting and important systems due to the large number of Slater determinants that have to be stored in a computer’s memory. Thus electronic structure theorists have made chemically and mathematically intuitive approximations to simplify the problems. A particularly difficult class of systems that is very interesting to the theoretical and experimental chemists is the strongly correlated systems. High Tc superconductors, nanotubes, graphene sheets, transition metal complexes and photosynthetic materials all fall under this category. These systems are too large to be treated by traditional quantum chemical methods. Another reason why such systems have not been studied sufficiently by ab-initio method is that they comprise of complicated strong electron-electron correlation (and therefore the name strongly correlated), thus, making them impossible to be qualitatively defined by a molecular orbital picture, which has been a quantum chemist’s fa- vorite tool for a long time now. Therefore, in this thesis, we have tried to develop, improve and apply methods for treating such strongly correlated systems. The thesis is broadly divided into two parts. In the first part, we discuss the orbital optimized DMRG method developed to treat static correlation in large strongly correlated systems. In the next part we discuss the two methods that can treat dynamic correlation in such complicated systems.
We present a new implementation of density matrix renormalization group based tailored coupled clusters method (TCCSD), which employs the domain-based local pair natural orbital approach (DLPNO). Compared to the previous local pair natural orbital (LPNO) version of the method, the new implementation is more accurate, offers more favorable scaling, and provides more consistent behavior across the variety of systems. On top of the singles and doubles, we include the perturbative triples correction (T), which is able to retrieve even more dynamic correlation. The methods were tested on three systems: tetramethyleneethane, oxo-Mn(Salen), and iron(II)-porphyrin model. The first two were revisited to assess the performance with respect to LPNO-TCCSD. For oxo-Mn(Salen), we retrieved between 99.8 and 99.9% of the total canonical correlation energy which is an improvement of 0.2% over the LPNO version in less than 63% of the total LPNO runtime. Similar results were obtained for iron(II)-porphyrin. When the perturbative triples correction was employed, irrespective of the active space size or system, the obtained energy differences between two spin states were within the chemical accuracy of 1 kcal/mol using the default DLPNO settings.
We present, to the best of our knowledge, the first attempt to exploit the super-computer platform for quantum chemical density matrix renormalization group (QC-DMRG) calculations. We have developed the parallel scheme based on the in-house MPI global memory library, which combines operator and symmetry sector parallelisms, and tested its performance on three different molecules, all typical candidates for QC-DMRG calculations. In case of the largest calculation, which is the nitrogenase FeMo cofactor cluster with the active space comprising 113 electrons in 76 orbitals and bond dimension equal to 6000, our parallel approach scales up to approximately 2000 CPU cores.
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.