Predicting the FCI Energy of Large Systems to Chemical Accuracy from Restricted Active Space Density Matrix Renormalization Group Calculations.

We theoretically derive and validate with large scale simulations a remarkably accurate power law scaling of errors for the restricted active space density matrix renormalization group (DMRG-RAS) method [J. Phys. Chem. A 126, 9709] in electronic structure calculations. This yields a new extrapolation method, DMRG-RAS-X, which reaches chemical accuracy for strongly correlated systems such as the chromium dimer, dicarbon up to a large cc-pVQZ basis and even a large chemical complex such as the FeMoco with significantly lower computational demands than those of previous methods. The method is free of empirical parameters, performed robustly and reliably in all examples we tested, and has the potential to become a vital alternative method for electronic structure calculations in quantum chemistry and more generally for the computation of strong correlations in nuclear and condensed matter physics.

Hilbert space multireference coupled cluster tailored by matrix product states.

In the past decade, the quantum chemical version of the density matrix renormalization group method has established itself as the method of choice for strongly correlated molecular systems. However, despite its favorable scaling, in practice, it is not suitable for computations of dynamic correlation. Several approaches to include that in post-DMRG methods exist; in our group, we focused on the tailored coupled cluster (TCC) approach. This method works well in many situations; however, in exactly degenerate cases (with two or more determinants of equal weight), it exhibits a bias toward the reference determinant representing the Fermi vacuum. Although sometimes it is possible to use a compensation scheme to avoid this bias for energy differences, it is certainly a drawback. In order to overcome this bias of the TCC method, we have developed a Hilbert-space multireference version of tailored CC, which can treat several determinants on an equal footing. We have implemented and compared the performance of three Hilbert-space multireference coupled cluster (MRCC) variants-the state universal one and the Brillouin-Wigner and Mukherjee's state specific ones. We have assessed these approaches on the cyclobutadiene and tetramethyleneethane molecules, which are both diradicals with exactly degenerate determinants at a certain geometry. We have also investigated the sensitivity of the results on the orbital rotation of the highest occupied and lowest unoccupied molecular orbital (HOMO-LUMO) pair, as it is well known that Hilbert-space MRCC methods are not invariant to such transformations.

Toward Large-Scale Restricted Active Space Calculations Inspired by the Schmidt Decomposition.

We present an alternative, memory-efficient, Schmidt decomposition-based description of the inherently bipartite restricted active space (RAS) scheme, which can be implemented effortlessly within the density matrix renormalization group (DMRG) method via the dynamically extended active space procedure. Benchmark calculations are compared against state-of-the-art results of C2 and Cr2, which are notorious for their multireference character. Our results for ground and excited states together with spectroscopic constants demonstrate that the proposed novel approach, dubbed as DMRG-RAS, which is variational and free of uncontrolled method errors, has the potential to outperfom conventional methods for strongly correlated molecules.

Assessing the Accuracy of Tailored Coupled Cluster Methods Corrected by Electronic Wave Functions of Polynomial Cost.

Tailored coupled cluster theory represents a computationally inexpensive way to describe static and dynamical electron correlation effects. In this work, we scrutinize the performance of various coupled cluster methods tailored by electronic wave functions of polynomial cost. Specifically, we focus on frozen-pair coupled cluster (fpCC) methods, which are tailored by pair-coupled cluster doubles (pCCD), and coupled cluster theory tailored by matrix product state wave functions optimized by the density matrix renormalization group (DMRG) algorithm. As test system, we selected a set of various small- and medium-sized molecules containing diatomics (N2, F2, C2, CN+, CO, BN, BO+, and Cr2) and molecules (ammonia, ethylene, cyclobutadiene, benzene, hydrogen chains, rings, and cuboids) for which the conventional single-reference coupled cluster singles and doubles (CCSD) method is not able to produce accurate results for spectroscopic constants, potential energy surfaces, and barrier heights. Most importantly, DMRG-tailored and pCCD-tailored approaches yield similar errors in spectroscopic constants and potential energy surfaces compared to accurate theoretical and/or experimental reference data. Although fpCC methods provide a reliable description for the dissociation pathway of molecules featuring single and quadruple bonds, they fail in the description of triple or hextuple bond-breaking processes or avoided crossing regions.

Orbital entanglement and correlation from pCCD-tailored coupled cluster wave functions.

Wave functions based on electron-pair states provide inexpensive and reliable models to describe quantum many-body problems containing strongly correlated electrons, given that broken-pair states have been appropriately accounted for by, for instance, a posteriori corrections. In this article, we analyze the performance of electron-pair methods in predicting orbital-based correlation spectra. We focus on the (orbital-optimized) pair-coupled cluster doubles (pCCD) ansatz with a linearized coupled-cluster (LCC) correction. Specifically, we scrutinize how orbital-based entanglement and correlation measures can be determined from a pCCD-tailored CC wave function. Furthermore, we employ the single-orbital entropy, the orbital-pair mutual information, and the eigenvalue spectra of the two-orbital reduced density matrices to benchmark the performance of the LCC correction for the one-dimensional Hubbard model with the periodic boundary condition as well as the N2 and F2 molecules against density matrix renormalization group reference calculations. Our study indicates that pCCD-LCC accurately reproduces the orbital-pair correlation patterns in the weak correlation limit and for molecules close to their equilibrium structure. Hence, we can conclude that pCCD-LCC predicts reliable wave functions in this regime.

DMRG on Top of Plane-Wave Kohn-Sham Orbitals: A Case Study of Defected Boron Nitride.

In this paper, we analyze the numerical aspects of the inherent multireference density matrix renormalization group (DMRG) calculations on top of the periodic Kohn-Sham density functional theory using the complete active space approach. The potential of the framework is illustrated by studying hexagonal boron nitride nanoflakes embedding a charged single boron vacancy point defect by revealing a vertical energy spectrum with a prominent multireference character. We investigate the consistency of the DMRG energy spectrum from the perspective of sample size, basis size, and active space selection protocol. Results obtained from standard quantum chemical atom-centered basis calculations and plane-wave based counterparts show excellent agreement. Furthermore, we also discuss the spectrum of the periodic sheet which is in good agreement with extrapolated data of finite clusters. These results pave the way toward applying the DMRG method in extended correlated solid-state systems, such as point defect qubit in wide band gap semiconductors.

Massively parallel quantum chemical density matrix renormalization group method.

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.

Topologically Protected Correlated End Spin Formation in Carbon Nanotubes.

For most chiralities, semiconducting nanotubes display topologically protected end states of multiple degeneracies. We demonstrate using density matrix renormalization group based quantum chemistry tools that the presence of Coulomb interactions induces the formation of robust end spins. These are the close analogs of ferromagnetic edge states emerging in graphene nanoribbons. The interaction between the two ends is sensitive to the length of the nanotube, its dielectric constant, and the size of the end spins: for S=1/2 end spins, their interaction is antiferromagnetic, while for S>1/2, it changes from antiferromagnetic to ferromagnetic as the nanotube length increases. The interaction between end spins can be controlled by changing the dielectric constant of the environment, thereby providing a possible platform for two-spin quantum manipulations.

Ground state of the Fe(ii)-porphyrin model system corresponds to quintet: a DFT and DMRG-based tailored CC study.

Fe(ii)-porphyrins play an important role in many reactions relevant to material science and biological processes, due to their closely lying spin states. Although the prevalent opinion is that these systems posses the triplet ground state, the recent experiment on Fe(ii)-phthalocyanine under conditions matching those of an isolated molecule points toward the quintet ground state. We present a thorough DFT and DMRG-based tailored CC study of Fe(ii)-porphyrin model, in which we address all previously discussed correlation effects. We examine the importance of geometrical parameters, the Fe-N distances in particular, and conclude that the system possesses the quintet ground state.

Toward DMRG-tailored coupled cluster method in the 4c-relativistic domain.

There are three essential problems in computational relativistic chemistry: Electrons moving at relativistic speeds, close lying states, and dynamical correlation. Currently available quantum-chemical methods are capable of solving systems with one or two of these issues. However, there is a significant class of molecules in which all the three effects are present. These are the heavier transition metal compounds, lanthanides, and actinides with open d or f shells. For such systems, sufficiently accurate numerical methods are not available, which hinders the application of theoretical chemistry in this field. In this paper, we combine two numerical methods in order to address this challenging class of molecules. These are the relativistic versions of coupled cluster methods and the density matrix renormalization group (DMRG) method. To the best of our knowledge, this is the first relativistic implementation of the coupled cluster method externally corrected by DMRG. The method brings a significant reduction of computational costs as we demonstrate on the system of TlH, AsH, and SbH.

Near-Linear Scaling in DMRG-Based Tailored Coupled Clusters: An Implementation of DLPNO-TCCSD and DLPNO-TCCSD(T).

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.

Toward the efficient local tailored coupled cluster approximation and the peculiar case of oxo-Mn(Salen).

We introduce a new implementation of the coupled cluster method with single and double excitations tailored by the matrix product state wave functions (DMRG-TCCSD), which employs the local pair natural orbital (LPNO) approach. By exploiting locality in the coupled cluster stage of the calculation, we were able to remove some of the limitations that hindered the application of the canonical version of the method to larger systems and/or with larger basis sets. We assessed the accuracy of the approximation using two systems: tetramethyleneethane (TME) and oxo-Mn(Salen). Using the default cut-off parameters, we were able to recover over 99.7% and 99.8% of the canonical correlation energy for the triplet and singlet state of TME, respectively. In the case of oxo-Mn(Salen), we found that the amount of retrieved canonical correlation energy depends on the size of the complete active space (CAS)-we retrieved over 99.6% for the larger 27 orbital CAS and over 99.8% for the smaller 22 orbital CAS. The use of LPNO-TCCSD allowed us to perform these calculations up to quadruple-ζ basis set, amounting to 1178 basis functions. Moreover, we examined dependence of the ground state of oxo-Mn(Salen) on the CAS composition. We found that the inclusion of 4dxy orbital plays an important role in stabilizing the singlet state at the DMRG-CASSCF level via double-shell effect. However, by including dynamic correlation, the ground state was found to be triplet regardless of the size of the basis set or the composition of CAS, which is in agreement with previous findings by canonical DMRG-TCCSD in smaller basis.

Quantum information-based analysis of electron-deficient bonds.

Recently, the correlation theory of the chemical bond was developed, which applies concepts of quantum information theory for the characterization of chemical bonds, based on the multiorbital correlations within the molecule. Here, for the first time, we extend the use of this mathematical toolbox for the description of electron-deficient bonds. We start by verifying the theory on the textbook example of a molecule with three-center two-electron bonds, namely, diborane(6). We then show that the correlation theory of the chemical bond is able to properly describe the bonding situation in more exotic molecules which have been synthesized and characterized only recently, in particular, the diborane molecule with four hydrogen atoms [diborane(4)] and a neutral zerovalent s-block beryllium complex, whose surprising stability was attributed to a strong three-center two-electron π bond stretching across the C-Be-C core. Our approach is of high importance especially in the light of a constant chase after novel compounds with extraordinary properties where the bonding is expected to be unusual.

Numerical and Theoretical Aspects of the DMRG-TCC Method Exemplified by the Nitrogen Dimer.

In this article, we investigate the numerical and theoretical aspects of the coupled-cluster method tailored by matrix-product states. We investigate formal properties of the used method, such as energy size consistency and the equivalence of linked and unlinked formulation. The existing mathematical analysis is here elaborated in a quantum chemical framework. In particular, we highlight the use of what we have defined as a complete active space-external space gap describing the basis splitting between the complete active space and the external part generalizing the concept of a HOMO-LUMO gap. Furthermore, the behavior of the energy error for an optimal basis splitting, i.e., an active space choice minimizing the density matrix renormalization group-tailored coupled-cluster singles doubles error, is discussed. We show numerical investigations on the robustness with respect to the bond dimensions of the single orbital entropy and the mutual information, which are quantities that are used to choose a complete active space. Moreover, the dependence of the ground-state energy error on the complete active space has been analyzed numerically in order to find an optimal split between the complete active space and external space by minimizing the density matrix renormalization group-tailored coupled-cluster error.

Elucidating cation-cation interactions in neptunyl dications using multi-reference ab initio theory.

Understanding the binding mechanism in neptunyl clusters formed due to cation-cation interactions is of crucial importance in nuclear waste reprocessing and related areas of research. Since experimental manipulations with such species are often rather limited, we have to rely on quantum-chemical predictions of their electronic structures and spectroscopic parameters. In this work, we present a state-of-the-art quantum chemical study of the T-shaped and diamond-shaped neptunyl(v) and neptunyl(vi) dimers. Specifically, we scrutinize their molecular structures, (implicit and explicit) solvation effects, the interplay of static and dynamical correlation, and the influence of spin-orbit coupling on the ground state and lowest-lying excited states for different total spin states and total charges of the neptunyl dications. Furthermore, we use the picture of interacting orbitals (quantum entanglement and correlation analysis) to identify strongly correlated orbitals in the cation-cation complexes that should be included in complete active space calculations. Most importantly, our study highlights the complex interplay of correlation effects and relativistic corrections in the description of the ground and lowest-lying excited states of neptunyl dications.

An Isolated Molecule of Iron(II) Phthalocyanin Exhibits Quintet Ground-State: A Nexus between Theory and Experiment.

Iron(II) phthalocyanine (FePc) is an important member of the phthalocyanines family with potential applications in the fields of electrocatalysis, magnetic switching, electrochemical sensing, and phototheranostics. Despite the importance of electronic properties of FePc in these applications, a reliable determination of its ground-state is still challenging. Here we present combined state of the art computational methods and experimental approaches, that is, Mössbauer spectroscopy and Superconducting Quantum Interference Device (SQUID) magnetic measurements to identify the ground state of FePc. While the nature of the ground state obtained with density functional theory (DFT) depends on the functional, giving mostly the triplet state, multi-reference complete active space second-order perturbation theory (CASPT2) and density matrix renormalization group (DMRG) methods assign quintet as the FePc ground-state in gas-phase. This has been confirmed by the hyperfine parameters obtained from 57 Fe Mössbauer spectroscopy performed in frozen monochlorobenzene. The use of monochlorobenzene guarantees an isolated nature of the FePc as indicated by a zero Weiss temperature. The results open doors for exploring the ground state of other metal porphyrin molecules and their controlled spin transitions via external stimuli.

The Intricate Case of Tetramethyleneethane: A Full Configuration Interaction Quantum Monte Carlo Benchmark and Multireference Coupled Cluster Studies.

We have performed a full configuration interaction (FCI) quality benchmark calculation for the tetramethyleneethane molecule in the cc-pVTZ basis set employing a subset of complete active space second order perturbation theory, CASPT2(6,6), natural orbitals for the FCI quantum Monte Carlo calculation. The results are in an excellent agreement with the previous large scale diffusion Monte Carlo calculations by Pozun et al. and available experimental results. Our computations verified that there is a maximum on the potential energy surface (PES) of the ground singlet state (1A) 45° torsional angle, and the corresponding vertical singlet-triplet energy gap is 0.01 eV. We have employed this benchmark for the assessment of the accuracy of Mukherjee's coupled clusters with up to triple excitations (MkCCSDT) and CCSD tailored by the density matrix renormalization group method (DMRG). Multireference MkCCSDT with CAS(2,2) model space, though giving good values for the singlet-triplet energy gap, is not able to properly describe the shape of the multireference singlet PES. Similarly, DMRG(24,25) is not able to correctly capture the shape of the singlet surface, due to the missing dynamic correlation. On the other hand, the DMRG-tailored CCSD method describes the shape of the ground singlet state with excellent accuracy but for the correct ordering requires computation of the zero-spin-projection component of the triplet state (3B1).

T3NS: Three-Legged Tree Tensor Network States.

We present a new variational tree tensor network state (TTNS) ansatz, the three-legged tree tensor network state (T3NS). Physical tensors are interspersed with branching tensors. Physical tensors have one physical index and at most two virtual indices, as in the matrix product state (MPS) ansatz of the density matrix renormalization group (DMRG). Branching tensors have no physical index, but up to three virtual indices. In this way, advantages of DMRG, in particular a low computational cost and a simple implementation of symmetries, are combined with advantages of TTNS, namely incorporating more entanglement. Our code is capable of simulating quantum chemical Hamiltonians, and we present several proof-of-principle calculations on LiF, N2, and the bis(µ-oxo) and µ-η2:η2 peroxo isomers of [Cu2O2]2+.

The correlation theory of the chemical bond.

The quantum mechanical description of the chemical bond is generally given in terms of delocalized bonding orbitals, or, alternatively, in terms of correlations of occupations of localised orbitals. However, in the latter case, multiorbital correlations were treated only in terms of two-orbital correlations, although the structure of multiorbital correlations is far richer; and, in the case of bonds established by more than two electrons, multiorbital correlations represent a more natural point of view. Here, for the first time, we introduce the true multiorbital correlation theory, consisting of a framework for handling the structure of multiorbital correlations, a toolbox of true multiorbital correlation measures, and the formulation of the multiorbital correlation clustering, together with an algorithm for obtaining that. These make it possible to characterise quantitatively, how well a bonding picture describes the chemical system. As proof of concept, we apply the theory for the investigation of the bond structures of several molecules. We show that the non-existence of well-defined multiorbital correlation clustering provides a reason for debated bonding picture.