Using Higher-Order Singular Value Decomposition To Define Weakly

Aug 22, 2017 - Using Higher-Order Singular Value Decomposition To Define Weakly Coupled and Strongly Correlated Clusters: The n-Body Tucker Approximat...
3 downloads 8 Views 2MB Size
Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES

Article

Using higher-order singular value decomposition to define weakly coupled and strongly correlated cluster states: the n-body Tucker approximation Nicholas J. Mayhall J. Chem. Theory Comput., Just Accepted Manuscript • DOI: 10.1021/acs.jctc.7b00696 • Publication Date (Web): 22 Aug 2017 Downloaded from http://pubs.acs.org on August 22, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

Journal of Chemical Theory and Computation is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 14

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Journal of Chemical Theory and Computation

Using higher-order singular value decomposition to define weakly coupled and strongly correlated clusters: the n-body Tucker approximation Nicholas J. Mayhall∗ Department of Chemistry, Virginia Tech, Blacksburg, VA 24060, USA An approximate wavefunction ansatz is presented which describes low-energy states of a highly clustered molecular system as a linear combination of multiple reduced-rank tensors. Using the Tucker decomposition as a way to obtain local clusters states, the exact solution is solved for in the space spanned by a small number of states on each cluster, with complete correlation occuring between limited numbers of clusters at a time. In this initial study, we report the implementation for a Heisenberg spin Hamiltonian with numerical examples of regular grid spin lattices, and ab initio-derived spin-Hamiltonians used to analyze the approximation. From these results, we find that the proposed method works well when the Hamiltonian interactions within a cluster are larger than between a cluster, and when this is not true, the method is not effective.

I.

mode product:

INTRODUCTION

Tensors, or multilinear arrays, are useful mathematical objects for storing and manipulating higher-dimensional datasets where data entries carry several attributes. However, due to the fact that each additional dimension significantly increases the number of possible data entries, mathematical techniques for efficient storage and analysis of tensors is becoming increasingly important. One particularly important group of techniques are referred to as tensor decompositions,1–4 which has been applied to many data intensive problems, i.e., data mining,5,6 facial recognition,7 signal processing.8 These are attempts to generalize the techniques of matrix factorizations, such as the eigen-decomposition or singular value decomposition (SVD), so as to obtain a rank-exposing orthogonal subspace representation of the tensor. Although, there is no complete SVD generalization to n-dimensional tensors which retains all of the desirable properties of the SVD, there are a multiple useful decompositions which have been studied. Most of these fall into one of two categories, 1) PARAFAC/CANDECOMP3 or 2) Tucker decompositions (Higher-Order SVD, HOSVD).4 In PARAFAC, one attempts to find a compact representation of a d-mode input tensor T as a sum of rank-1 tensors (i.e., a sum of kronecker products of d vectors), T ≈

R X

(2) (d) λr u(1) r ⊗ ur ⊗ · · · ⊗ ur ,

(1)

r

where R is the tensor rank, and λr is a scalar coefficient. (d) Whereas with the Tucker decomposition, the vectors ur are replaced with matrices and the coefficient becomes another tensor, C: T ≈ C ×1 U (1) ×2 U (2) · · · ×d U (d) ,

(2)

where C is a smaller “core tensor”, U (d) is a unitary matrix in the vector space associated with mode d of T called a “Tucker factor”, and C ×d U (d) denotes the d-

C ×n U (n) ⇔Ci1 ,i2 ,...,in−1 ,in ,in+1 ,... Uin ,j = Ci1 ,i2 ,...,in−1 ,j,in+1 ,...

(3)

The Tucker decomposition can also be expressed with index notation: Ti,j,...,d ≈ Cα,β,...,γ Ui,α Uj,β · · · Ud,γ

(4)

As is easily seen in Eq. 4, the Tucker decomposition is essentially a change of basis, from i to α, from j to β, etc. Without detailing the various computational algorithms for performing PARAFAC and Tucker decompositions (which is outside the scope of this paper), we simply emphasize that it is often desirable to have accurate approximations to the original tensor, but which contain much less data, resulting in more efficient storage, algebraic manipulation, and analysis. Although much of the literature pertaining to tensor decompositions exist outside of the physical sciences, there are many areas of science which benefit from efficient tensor decomposition strategies. Due to the high dimensionality of the associated quantities, quantum chemistry is an ideal example of such an area of science. During the past several years, tensor decomposition techniques have become increasingly valuable for increasing the efficiency of several quantum chemistry algorithms. Most notably perhaps, has been the significant acceleration of the two-electron repulsion integrals based on Resolution of the Identity (RI) approximations or Cholesky Decompositions (CD).9–16,16–22 More recently, the Tensor HyperContraction (THC) methodology23–26 (and related methods),27 have been developed to decompose the 4-mode two-electron integral tensor, which can reduce both storage requirements and overall scaling of methods like MP2. The MP2 or coupled cluster amplitude tensors have also been factorized via SVD,28,29 , PARAFAC,30,31 , THC,24 and Tucker decompositions,32 to understand the effect correlation has on orbitals, and to improve computational efficiency. Compression techniques have also been used to obtain optimized orbital active spaces for aggressively truncating configuration expansions in corre-

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 2 of 14 2

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

lated calculations.33–41 A nice overview of the role of tensor decompositions in quantum chemistry can be found in Ref. 42. Tensor decompositions also have intimate connections to several emerging techniques for low-entanglement wavefunctions for strongly correlated systems.43 These are in contrast to the conventional methods for strong correlation, such as complete active space self consistent field (CASSCF),44 which this is effectively a bruteforce procedure and only possible when the number of strongly correlated orbitals is small. Although originally developed as a way to exploit scale separation in spin lattices,45–47 density matrix renormalization group (DMRG) soon became recognized as an algorithm which variationally optimizes a wavefunction having a specific decomposed tensor format known as matrix product state (MPS).43,48–52 Efficient codes have since demonstrated that DMRG can solve quantum chemistry problems which would otherwise be intractable.53–70 Because the MPS is a one-dimensional factorization of the FCI coefficient tensor, ideal performance only occurs for molecular systems with quasi one-dimensional geometry, although this has been found to not be overly restricting in many cases. Several new techniques such as complete graph tensor network states (CGTN)71 , tree tensor networks (TTN)72 , and multiscale renormalization ansatz (MERA)73 have been developed to exploit low-entanglement in higher dimensional systems, and many other strong correlation techniques are also being explored.50,74–84 While the approaches mentioned above attempt to find low-rank approximations by disentangling the indices of the FCI tensor, other approaches seek to find local subspaces which result in compact representations. Some of the approaches most relevant to this current work are the Active-Space Decomposition (ASD) of Shiozaki and coworkers,85–87 Block-Correlated Coupled Cluster (BCCC) of Li and coworkers,88–90 the cluster Mean-Field (cMF) method of Scuceria and coworkers,91 the Renormalized Exciton Method (REM) of Ma and coworkers,92–95 and the ab initio Frenkel Davydov Exciton Model (AIFDEM) of Herbert and coworkers.96–98 In ASD, the FCI space of a clustered system is decomposed as a direct product of low-energy states on each cluster. This is essentially an approximate Tucker decomposition of the FCI tensor. It is approximate in the sense that the “Tucker factors” are chosen to be eigenfunctions of the local cluster Hamiltonian, as opposed to the eigenvectors of a reduced density matrix (i.e., HOSVD). Parker et al. have shown86 that a benzene dimer’s lowenergy spectrum can be nearly converged by treating only tens of states per cluster. However, for larger numbers of clusters, even tens of states per cluster will quickly become intractable. In order to avoid this exponential increase in dimension, the DMRG algorithm was then used to approximate the ASD method for arbitrary numbers of clusters,85 with the increased correlation within a cluster leading to a decreased number of important renormalized

states compared to conventional DMRG. The BCCC method is a generalization of the coupled cluster ansatz to the case of a direct product reference wavefunction. The direct product wavefunction can capture important strong correlation effects within a block. Then, in order to capture the missing dynamical correlation, Li and coworkers88–90 use genaralized creation/annihilation operators which act to destroy a cluster’s reference state and replace it with a higher energy state (states which were discarded in the original definition of the direct product reference state). Similar to conventional coupled-cluster, these generalized excitation operators are then exponentiated and applied to the reference, with the resulting matrix elements derived in a similar (but vastly more complicated) manner as that used in single-reference coupled cluster. This has been very successful in quantum chemistry applications using both a CAS-CI reference,89,90,99,100 and as a perturbative correction to a GVB reference.101 In the cMF method, a direct product wavefunction is again assumed to describe the basic qualitative features of a given strongly correlated state. However, rather than choosing an arbitrary direct product state, Scuceria and Jim´enez-Hoyos variationally optimized the molecular orbitals by minimizing the energy of the direct product wavefunction.91 This orbital optimization was found to be critical for obtaining accurate results. After optimizing the orbitals, a second-order correction (referred to as cPT2) was included which recovered a significant amount of the missing correlation for a wide range of correlation strengths of a Hubbard model. The REM approach of Malrieu and coworkers,92 can also be considered as a tensor decomposed wavefunction. Here, the large system’s wavefunction is approximated as a linear combination of block excitations, in which only a single block is excited at a time, while all others remain in their ground states. Optimizing the wavefunction in this space requires the Hamiltonian matrix elements between the various block excited states. However, instead of using the bare Hamiltonian, Malrieu and coworkers use Bloch’s effective Hamiltonian theory102,103 to obtain an effective interaction between the block excited states. This implicitly folds in the effects of both higher energy block states and cooperative multi-block excitations into the model space Hamiltonian matrix. This has been shown to be effective on gapped spin-lattices,92,104 in addition to the excited states of chromophore aggregates using the ab initio Hamiltonian.93–95 Targeting similar problems as REM, the ab initio Frenkel Davydov Exciton Model (AIFDEM) of Herbert and coworkers also seeks to approximate the low-lying singly-excited states of aggregates via a direct product wavefunction ansatz. However, unlike REM, AIFDEM uses the bare Hamiltonian in the evaluation of the exciton couplings. Although this prevents the implicit inclusion of multiblock states, (i.e., partial hybridization between Frenkel and charge-transfer excitations), AIFDEM does include some explicit 2-block charge-transfer states by

ACS Paragon Plus Environment

Page 3 of 14

Journal of Chemical Theory and Computation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 4 of 14 4

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The exact coefficient tensor can therefore be written as:     ˆ A PˆB + λQ ˆ B PˆC + λQ ˆC C =C PˆA + λQ (9)

increase the number of degrees of freedom, it is necessary to avoid the construction of the full Hamiltonian matrix by only evaluating the action of the Hamiltonian onto each of the Tucker block core tensors (σ vector). This is =C PˆA PˆB PˆC 0th Order essentially the same approach taken in conventional CI ˆ A PˆB PˆC + C PˆA Q ˆ B Pˆf C + C PˆA PˆB Q ˆ C 1st Order calculations using Davidson’s algorithm.106 The contri+ CQ ˆAQ ˆ B PˆC + C Q ˆ A PˆB Q ˆ C + C PˆA Q ˆB Q ˆ C 2nd Order bution to the σ vector will include a sum over all one + CQ and two cluster Hamiltonian terms. For example, the ˆAQ ˆB Q ˆC + CQ 3rd Order operator coupling clusters A and B will generate the following contribution to the σ vector, thus leading to a definition of the exact coefficient tensor k ˆ AB |Aα′ Bβ ′ Cγ ′ i Cαk′ β ′ γ ′ (13) σαβγ ⇐ hAα Bβ Cγ | H in terms of a sum of decreasingly important compressed  th X X tensors. Our aim is to truncate this sum at n order and Jab hAα Bβ | Sˆa+ Sˆb− + Sˆa− Sˆb+ = solve for the associated core tensors. a∈A b∈B As a simple example, suppose we were to truncate 

this expansion after first order terms. The wavefunction (14) +2Sˆaz Sˆbz |Aα′ Bβ ′ i Cγ Cγ′ Cαk′ β ′ γ ′ would then be written as,  XX + − − + Jab Sa,αα = ′ Sb,ββ ′ + Sa,αα′ Sb,ββ ′ |ψi =Cijk |Ai Bj Ck i (10) a∈A b∈B  k z z ≈ (Cαβγ Uiα Ujβ Ukγ + Cαβγ V U U (11) ¯ iα ¯ jβ kγ (15) + 2Sa,αα ′ Sb,ββ ′ Cα′ β ′ γ  k |A B C i + Cαβγ U V U + C U U V ¯ i j k iα j β¯ kγ αβ¯ γ iα jβ k¯ γ (16) =Hαβα′ β ′ Cα′ β ′ γ E E ˜ ˜ ˜ ˜ ˜ ˜ ≈Cαβγ Aα Bβ Cγ + Cαβγ Aα¯ Bβ Cγ ¯ where a and b correspond to lattice sites in clusters A E E and B, respectively. ˜ ˜ ˜ ˜ ˜ ˜ + Cαβγ (12) ¯ Aα Bβ¯ Cγ + Cαβ¯ γ Aα Bβ Cγ ¯ Because the Sˆz operator is diagonal we can easily use this quantity to project each Tucker block onto single Ms where Uiα and Viα¯ are the vectors in the P and Q spaces, block. In such a case, the dimer Hamiltonian in Eq. 16 respectively. Because these projection operators are decan be collapsed onto the given subspace before contractfined in terms of a HOSVD, we will refer to each of ing with the trial vectors. these sets Eof basis vectors as Tucker blocks. For instance, In the cases presented here, the non-interacting Hamil ˜ ˜ ˜ ˆ 1B = L H ˆ A , is used to precondition the Aα Bβ Cγ defines a Tucker block defined by the states tonian, H A Hamiltonian matrix in the Davidson algorithm. This spanned by α, β, γ. seems to be effective under the same conditions required In the simplest approach (and the one explored in this for the n-body Tucker method to be accurate (“highly initial paper) one could then use the variational principle clustered” systems). However, as inter-cluster couplings to optimize these coefficients. Although a variational opincrease relative to the gaps between the P and Q spaces timization does not provide a purely size-extensive theof the clusters, the preconditioner becomes less wellory, it does guarantee size-consistency in certain cases. suited, and other directions will need to be explored. For instance, a system comprised of an arbitrary number of clusters computed at infinity, is exactly equal to the sum of the individual cluster energies. This is a result B. Self-consistent cluster state optimization of the direct product reference. However, if one had four clusters arranged in two groups, with infinite distance beBy defining the cluster states in terms of the HOSVD tween the groups and non-zero inter-cluster interactions of the exact tensor, it is obvious that one does not diwithin a group, a 2-body Tucker approximation would rectly arrive at a computationally useful algorithm. The not yield size consistent results. In future work we will inexact solution is a prerequisite, so we need to define a vestigate the use of perturbation theory to determine the coefficient tensors, which would be exactly size-extensive. bottom-up strategy for obtaining accurate cluster states. In order to optimize the wavefunction, we will need If we were to use only the zeroth-order Tucker block, to evaluate matrix elements in this representation. In an HOSVD of the resulting tensor will always return the same states, since no mixing with the Q space is included. this first paper, we explore the implementation of a However, because the higher order Tucker blocks mix the model spin Hamiltonian to study the convergence properP and Q spaces, the sum of multiple tensors will allow the ties and behavior of this wavefunction ansatz, although occupation numbers to spill into the Q space, thus giving fully ab initio implementation will be the next step in us a new set of singular vectors, which changes the variathis work. Our implementation is inspired by Shiozaki’s tional space accessible within our truncated wavefunction Active-space Decomposition (ASD)85–87 in that we comapproximation. This naturally defines a self-consistency pute matrix elements using operators represented diloop between the cluster state partitioning and an eigenrectly in the associated local basis. However, because state of the full system Hamiltonian. In order to ensure the addition of Q space contributions will significantly

ACS Paragon Plus Environment

Page 5 of 14

Journal of Chemical Theory and Computation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 6 of 14

Page 7 of 14

Journal of Chemical Theory and Computation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 8 of 14 8

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

puted individually and added together. This also has an interesting interpretation from an quantum embedding perspective. For example, in an incremental formulation, a 2-body interaction between clusters A and C would be obtained by diagonalizing the Hamiltonian in the basis spanned by all states on A and C, while the rest of the clusters are compressed into their respective P spaces. As such, each resulting n-body group of clusters would be computed with a quantum bath, enabling direct entanglement between the n-body collection and the rest of the system. Preliminary data is promising for this approach, and will be the focus of future work.

B.

H4 Clusters

In real systems, the interactions are not only between nearest neighbor sites as in the above 2-d spin lattices. More generally, all sites can interact with one another with a magnitude which decays with increasing distance. In order to explore the performance of the nbody Tucker expansion with systems which contain interactions between all sites, we have constructed an ab initio model consisting of 24 hydrogen atoms in space by starting with a methane cluster and deleting the carbon atoms. These hydrogen atoms are thus arranged into near-tetrahedral clusters, which themselves are arranged in a near-octahedron. In this study, each neartetrahedral group of four hydrogen atoms defines a cluster for the n-body Tucker decomposition. A potential energy surface scan is explored in which one of the clusters is pushed into the center of the octahedron. This system has no exact symmetry, so no interactions are necessarily zero. This ab initio system is then mapped onto a spin lattice, using the spin-flip CI plus Bloch effective Hamiltonian strategy, which has been presented previously,109,110 and then solved with the currently proposed n-body Tucker method. At each position of the PES scan, an ab initio spin-flip calculation is performed, the exchange coupling constants are extracted, and the resulting spin Hamiltonian is solved. Each ab initio spinflip calculation has been performed with a development version of QChem.111 These results are presented in Fig. 6. This ab initio parameterization of a spin Hamiltonian has been demonstrated to provide qualitative agreement with ab initio results for many cases. In this particular case, we do not expect quantitative accuracy, due to the strong coupling of the hydrogen atoms when they become close, for which the Bloch effective Hamiltonian projection has a norm of only 0.747 for the most compact geometry. Nonetheless, this does provide a physically motivated (and qualitatively meaningful) parameterization of a spin lattice to better understand this methodology. The extracted model Hamiltonian parameters are available in the supplementary information for each point on the PES. Moving from right to left on this PES in Fig. 6(b),

one can see that all approaches work well for the system when they are separated. However, as the rightmost cluster moves closer to the others, the defined clustering pattern becomes non-ideal. The errors increase quickly for the 0-body Tucker curve and 2-body Tucker curve, while the 4-body Tucker curve remains nearly flat, with the largest error of 0.004 eV for the point at −0.65 ˚ A. The reason the 2-body Tucker method is unable to fully recover this interaction energy is easily seen by the exchange interaction graphs in Fig. 6(a). For this compressed geometry, the largest interactions in the system are between the translated cluster, and two other clusters. Thus, in order to contain the strongest interactions simultaneously, groups of at least three clusters would be needed, which are absent in the 2-body Tucker approximation. One way to improve the 2-body Tucker results would be to use a geometry dependent clustering. For example, if at r= −0.65 ˚ A a new clustering is chosen such that the strongest interactions are within a cluster (still with only four sites per cluster), the NB2 error is only 0.13 eV, compared to 1.72 eV with the original clustering. We plan on exploring better ways of clustering in the future, drawing inspiration from pair entanglement measures like the “mutual information” index from Reiher and Stein.112 Alternatively, one could increase the number of states in the P space for the clusters. If the three strongly coupled clusters are given four eigenstates each (one singlet and one triply degenerate triplet), a 1-body Tucker approximation can be used instead, providing an error of only 0.05 eV. In Fig. 6(c), the convergence of the cluster states is examined using the 2-body Tucker approximation for scan distance r=−0.65 ˚ A, which is the most compact geometry point. By comparing the energy convergence behavior of the DIIS extrapolated results, and the nonextrapolated results (simple iterations similar to taking Roothaan steps in Hartree-Fock), it is clearly seen that the uncoupled DIIS extrapolation is very effective in this case for improving the convergence of the cluster states (or Tucker factors), with the number of iterations required decreasing from over 1000 to 23, to obtain convergence in the energy change to within 10−9 eV. Without the DIIS extrapolations, the iterations only slowly change the energy for the first 100 iterations before rapidly decreasing the energy, before ultimately plateauing at the same point obtained with the DIIS extrapolation. The norms of the DIIS error vectors for each cluster is also shown in this figure, each represented with a grey line. Although DIIS is known to be very effective at accelerating self-consistency optimizations, the impressive performance here is perhaps surprising, due to the fact that the we are extrapolating each cluster’s reduced density matrix, which are local quantities, instead of a global quantity. As described above, this means that each cluster carries it’s own DIIS subspace, and that each extrapolation occurs independently of one another. While this seems to work well for the cases investigated so far, we do expect to find situations where this is less than ideal,

ACS Paragon Plus Environment

Page 9 of 14

Journal of Chemical Theory and Computation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 10 of 14 10

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

IV.

CONCLUSIONS

gether using an incremental approach for expectation values, and finally an ab initio implementation for fermionic systems to allow direct quantum chemistry calculations to be performed with this methodology.

In this paper, we propose a new wavefunction approximation for strongly correlated but highly clustered molecular systems, which combines Tucker tensor decompositions with more commonly used n-body organizations. The exact coefficient tensor is approximated as a sum of low-rank tensors (in the Tucker sense), and the resulting compressed coefficient tensors are determined using the variational principle. To explore the suitability of this approximation under different conditions, we have implemented this method to solve the Heisenberg spin Hamiltonian, which is a good model of molecules in the strongly correlated limit. Numerical calculations are performed for regular grid spin lattices, in addition to ab initio-derived spin lattices of different dimensionalities. From the data in this paper, we see that the n-body Tucker method works well for well-clustered systems, and incurs increasingly large errors when the clusters become less and less defined. Future work will focus on using perturbation theory and coupled-electron pairbased methods instead of the variational principle to parameterize the wavefunction, (which will recover full size-extensivity), avoiding the full-system expansion alto-

The author would like to thank an anonymous reviewer for pointing out the similarity of the presented n-body Tucker decomposition with the REM92–95 and AIFDEM methods.96,97 Support for this work was provided by the Department of Chemistry at Virginia Tech. The authors acknowledge Advanced Research Computing at Virginia Tech for providing computational resources and technical support that have contributed to the results reported within this paper. URL: http://www.arc.vt.edu

Electronic address: [email protected] Cattell, R. B. Parallel proportional profiles and other principles for determining the choice of factors by rotation. Psychometrika 1944, 9, 267–283. Appellof, C. J.; Davidson, E. R. Strategies for analyzing data from video fluorometric monitoring of liquid chromatographic effluents. Anal. Chem. 1981, 53, 2053–2056. Hitchcock, F. L. The Expression of a Tensor or a Polyadic as a Sum of Products. J. Math. Phys. 1927, 6, 164–189. Tucker, L. R. Some mathematical notes on three-mode factor analysis. Psychometrika 1966, 31, 279–311. Kobayashi, M.; Aono, M. Vector Space Models for Search and Cluster Mining. In Surv. Text Min.; Springer New York: New York, NY, 2004; pp 103–122. Bader, B. W.; Berry, M. W.; Browne, M. Discussion Tracking in Enron Email Using PARAFAC. In Surv. Text Min. II ; Springer London: London, 2008; pp 147–163. Vasilescu, M. A. O.; Terzopoulos, D. Multilinear Analysis of Image Ensembles: TensorFaces. In Multilinear Analysis of Image Ensembles: TensorFaces; Springer, Berlin, Heidelberg, 2002; pp 447–460. De Lathauwer, L.; Vandewalle, J. Dimensionality reduction in higher-order signal processing and rank(R1,R2,. . . ,RN) reduction in multilinear algebra. Linear Algebra Appl. 2004, 391, 31–55. Wilson, S. Universal basis sets and Cholesky decomposition of the two-electron integral matrix. Comput. Phys. Commun. 1990, 58, 71–81. Sch¨ utz, M.; Manby, F. R.; Sch¨ utz, M.; Ahlrichs, R. Linear scaling local coupled cluster theory with density fitting. Part I: 4-external integrals. Phys. Chem. Chem. Phys. 2003, 5, 3349–3358. Werner, H.-J.; Manby, F. R.; Knowles, P. J. Fast linear

scaling second-order Møller-Plesset perturbation theory (MP2) using local and density fitting approximations. J. Chem. Phys. 2003, 118, 8149–8160. Manby, F. R. Density fitting in second-order linear-r12 MøllerPlesset perturbation theory. J. Chem. Phys. 2003, 119, 4607–4613. Aquilante, F.; Pedersen, T. B. Quartic scaling evaluation of canonical scaled opposite spin second-order MøllerPlesset correlation energy using Cholesky decompositions. Chem. Phys. Lett. 2007, 449, 354–357. Aquilante, F.; Pedersen, T. B.; Lindh, R. Low-cost evaluation of the exchange Fock matrix from Cholesky and density fitting representations of the electron repulsion integrals. J. Chem. Phys. 2007, 126, 194106. Boman, L.; Koch, H.; S´ anchez de Mer´ as, A. Method specific Cholesky decomposition: Coulomb and exchange energies. J. Chem. Phys. 2008, 129, 134107. Weigend, F.; Kattannek, M.; Ahlrichs, R. Approximated electron repulsion integrals: Cholesky decomposition versus resolution of the identity methods. J. Chem. Phys. 2009, 130, 164106. Aquilante, F.; Gagliardi, L.; Pedersen, T. B.; Lindh, R. Atomic Cholesky decompositions: a route to unbiased auxiliary basis sets for density fitting approximation with tunable accuracy and efficiency. J. Chem. Phys. 2009, 130, 154107. Zienau, J.; Clin, L.; Doser, B.; Ochsenfeld, C. Choleskydecomposed densities in Laplace-based second-order MøllerPlesset perturbation theory. J. Chem. Phys. 2009, 130, 204112. Manzer, S.; Horn, P. R.; Mardirossian, N.; HeadGordon, M. Fast, accurate evaluation of exact exchange: The occ-RI-K algorithm. J. Chem. Phys. 2015, 143,

∗ 1

2

3

4

5

6

7

8

9

10

11

Supporting Information

All Cartesian coordinates for all structures used in this paper can be found in the Supporting Information.

Acknowledgement

12

13

14

15

16

17

18

19

ACS Paragon Plus Environment

Page 11 of 14

Journal of Chemical Theory and Computation 11

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

20

21

22

23

24

25

26

27

28

29

30

31

32

33

34

35

36

024113. Manzer, S. F.; Epifanovsky, E.; Head-Gordon, M. Efficient Implementation of the Pair Atomic Resolution of the Identity Approximation for Exact Exchange for Hybrid and Range-Separated Density Functionals. J. Chem. Theory Comput. 2015, 11, 518–527. Hollman, D. S.; Schaefer, H. F.; Valeev, E. F. A tight distance-dependent estimator for screening three-center Coulomb integrals over Gaussian basis functions. J. Chem. Phys. 2015, 142, 154106. Hollman, D. S.; Schaefer, H. F.; Valeev, E. F. Semi-exact concentric atomic density fitting: Reduced cost and increased accuracy compared to standard density fitting. J. Chem. Phys. 2014, 140, 064109. Hohenstein, E. G.; Parrish, R. M.; Mart´ınez, T. J. Tensor hypercontraction density fitting. I. Quartic scaling second- and third-order Møller-Plesset perturbation theory. J. Chem. Phys. 2012, 137, 044103. Hohenstein, E. G.; Parrish, R. M.; Sherrill, C. D.; Mart´ınez, T. J. Communication: Tensor hypercontraction. III. Least-squares tensor hypercontraction for the determination of correlated wavefunctions. J. Chem. Phys. 2012, 137, 221101. Parrish, R. M.; Hohenstein, E. G.; Mart´ınez, T. J.; Sherrill, C. D. Tensor hypercontraction. II. Least-squares renormalization. J. Chem. Phys. 2012, 137, 224106. Parrish, R. M.; Hohenstein, E. G.; Schunck, N. F.; Sherrill, C. D.; Martinez, T. J. Tensor hypercontraction: A universal technique for the resolution of matrix elements of local, finite-range $N$-body potentials in many-body quantum problems 2013, 4. Schmitz, G.; Madsen, N. K.; Christiansen, O. Atomicbatched tensor decomposed two-electron repulsion integrals. The Journal of Chemical Physics 2017, 146, 134112. Kinoshita, T.; Hino, O.; Bartlett, R. J. Singular value decomposition approach for the approximate coupledcluster method. J. Chem. Phys. 2003, 119, 7756. Beran, G. J. O.; Head-Gordon, M. Extracting dominant pair correlations from many-body wave functions. J. Chem. Phys. 2004, 121, 78–88. Benedikt, U.; Auer, A. A.; Espig, M.; Hackbusch, W. Tensor decomposition in post-Hartree-Fock methods. I. Twoelectron integrals and MP2. J. Chem. Phys. 2011, 134, 054118. Benedikt, U.; Bohm, K.-H.; Auer, A. A. Tensor decomposition in post-HartreeFock methods. II. CCD implementation. J. Chem. Phys. 2013, 139, 224101. Bell, F.; Lambrecht, D.; Head-Gordon, M. Higher order singular value decomposition in quantum chemistry. Mol. Phys. 2010, 108, 2759–2773. Pito´ ak, M.; Holka, F.; Neogr´ ady, P.; Urban, M. Optimized virtual orbitals for correlated calculations: Towards large scale CCSD(T) calculations of molecular dipole moments and polarizabilities. J. Mol. Struct. THEOCHEM 2006, 768, 79–89. Adamowicz, L.; Bartlett, R. J.; Sadlej, A. J. Optimized virtual orbital space for highlevel correlated calculations. II. Electric properties. J. Chem. Phys. 1988, 88, 5749– 5758. Adamowicz, L.; Bartlett, R. J. Optimized virtual orbital space for highlevel correlated calculations. J. Chem. Phys. 1987, 86, 6314–6324. Ahlrichs, R.; Driessler, F.; Lischka, H.; Staemmler, V.;

37

38

39

40

41

42

43

44

45

46

47

48

49

50

51

52

53

54

55

56

Kutzelnigg, W. PNO-CI (pair natural orbital configuration interaction) and CEPA-PNO (coupled electron pair approximation with pair natural orbitals) calculations of molecular systems. II. The molecules BeH2 , BH, BH3 , CH4 , CH− . J. Chem. Phys. 1975, 62, 1235–1247. Ahlrichs, R.; Driessler, F. Direct determination of pair natural orbitals. Theor. Chim. Acta 1975, 36, 275–287. Meyer, W. PNOCI Studies of electron correlation effects. I. Configuration expansion by means of nonorthogonal orbitals, and application to the ground state and ionized states of methane. J. Chem. Phys. 1973, 58, 1017–1035. Meyer, W. Ionization energies of water from PNO-CI calculations. Int. J. Quantum Chem. 2009, 5, 341–348. Neese, F.; Wennmohs, F.; Hansen, A. Efficient and accurate local approximations to coupled-electron pair approaches: An attempt to revive the pair natural orbital method. J. Chem. Phys. 2009, 130, 114108. Taylor, P. R. A rapidly convergent CI expansion based on several reference configurations, using optimized correlating orbitals. J. Chem. Phys. 1981, 74, 1256–1270. Benedikt, U. PhD, Chemnitz University of Technology, 2013. Chan, G. K. L. Low entanglement wavefunctions. Wiley Interdiscip. Rev. Comput. Mol. Sci. 2012, 2, 907–920. Roos, B. O. The complete active space SCF method in a fock-matrix-based super-CI formulation. Int. J. Quantum Chem. 1980, 18, 175–189. White, S.; Noack, R. Real-space quantum renormalization groups. Phys. Rev. Lett. 1992, 68, 3487–3490. White, S. R. Density matrix formulation for quantum renormalization groups. Phys. Rev. Lett. 1992, 69, 2863– 2866. White, S. R. Density-matrix algorithms for quantum renormalization groups. Phys. Rev. B 1993, 48, 10345– 10356. Chan, G. K.-L.; Dorando, J. J.; Ghosh, D.; Hachmann, J.; Neuscamman, E.; Wang, H.; Yanai, T. An Introduction to the Density Matrix Renormalization Group Ansatz in Quantum Chemistry. In Front. Quantum Syst. Chem. Phys.; Springer Netherlands: Dordrecht, 2008; pp 49–65. Hachmann, J.; Cardoen, W.; Chan, G. K.-L. Multireference correlation in long molecules with the quadratic scaling density matrix renormalization group. J. Chem. Phys. 2006, 125, 144101. Chan, G. K.-L.; Zgid, D. The Density Matrix Renormalization Group in Quantum Chemistry. Annu. Rep. Comput. Chem. 2009, 5, 149–162. Wouters, S.; Van Neck, D. The density matrix renormalization group for ab initio quantum chemistry. Eur. Phys. J. D 2014, 68, 272. Schollw¨ ock, U. The density-matrix renormalization group in the age of matrix product states. Ann. Phys. (N. Y). 2011, 326, 96–192. White, S. R. Density-matrix algorithms for quantum renormalization groups. Phys. Rev. B 1993, 48, 10345– 10356. White, S. R.; Martin, R. L. Ab initio quantum chemistry using the density matrix renormalization group. J. Chem. Phys. 1999, 110, 4127. Chan, G. K.-L.; Head-Gordon, M. Highly correlated calculations with a polynomial cost algorithm: A study of the density matrix renormalization group. J. Chem. Phys. 2002, 116, 4462. Chan, G. K.-L.; Keselman, A.; Nakatani, N.; Li, Z.;

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

Page 12 of 14 12

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

57

58

59

60

61

62

63

64

65

66

67

68

69

70

71

72

73

74

White, S. R. Matrix Product Operators, Matrix Product States, and ab initio Density Matrix Renormalization Group algorithms 2016. Hachmann, J.; Dorando, J. J.; Avil´es, M.; Chan, G. K.L. The radical character of the acenes: a density matrix renormalization group study. J. Chem. Phys. 2007, 127, 134309. Chan, G. K.-L.; Sharma, S. The density matrix renormalization group in quantum chemistry. Annu. Rev. Phys. Chem. 2011, 62, 465–81. Sharma, S.; Chan, G. K.-L. Communication: A flexible multi-reference perturbation theory by minimizing the Hylleraas functional with matrix product states. J. Chem. Phys. 2014, 141, 111101. Sharma, S.; Alavi, A. Multireference linearized coupled cluster theory for strongly correlated systems using matrix product states. J. Chem. Phys. 2015, 143, 102815. Sharma, S.; Sivalingam, K.; Neese, F.; Chan, G. K.-L. Low-energy spectrum of iron-sulfur clusters directly from many-particle quantum mechanics. Nat. Chem. 2014, 6, 927–33. Hedeg˚ ard, E. D.; Reiher, M. Polarizable Embedding Density Matrix Renormalization Group. J. Chem. Theory Comput. 2016, 12, 4242–4253. Marti, K. H.; Reiher, M. The Density Matrix Renormalization Group Algorithm in Quantum Chemistry. Zeitschrift fur Phys. Chemie 2010, 224, 583–599. Stein, C. J.; von Burg, V.; Reiher, M. The Delicate Balance of Static and Dynamic Electron Correlation. J. Chem. Theory Comput. 2016, 12, 3764–3773. Marti, K. H.; Reiher, M. New electron correlation theories for transition metal chemistry. Phys. Chem. Chem. Phys. 2011, 13, 6750–9. Hedeg˚ ard, E. D.; Knecht, S.; Kielberg, J. S.; Jensen, H. J. A.; Reiher, M. Density matrix renormalization group with efficient dynamical electron correlation through range separation. J. Chem. Phys. 2015, 142, 224108. Mizukami, W.; Kurashige, Y.; Yanai, T. More π Electrons Make a Difference: Emergence of Many Radicals on Graphene Nanoribbons Studied by Ab Initio DMRG Theory. J. Chem. Theory Comput. 2013, 9, 401–407. Kurashige, Y.; Yanai, T. Second-order perturbation theory with a density matrix renormalization group selfconsistent field reference function: Theory and application to the study of chromium dimer. J. Chem. Phys. 2011, 135, 094104. Kurashige, Y.; Chan, G. K.-L.; Yanai, T. Entangled quantum electronic wavefunctions of the Mn4CaO5 cluster in photosystem II. Nat. Chem. 2013, 5, 660–666. Kurashige, Y.; Yanai, T. High-performance ab initio density matrix renormalization group method: Applicability to large-scale multireference problems for metal compounds. J. Chem. Phys. 2009, 130, 234114. Marti, K. H.; Bauer, B.; Reiher, M.; Troyer, M.; Verstraete, F. Complete-graph tensor network states: a new fermionic wave function ansatz for molecules. New J. Phys. 2010, 12, 103008. Shi, Y.-Y.; Duan, L.-M.; Vidal, G. Classical simulation of quantum many-body systems with a tree tensor network. Phys. Rev. A 2006, 74, 022320. Vidal, G. Entanglement Renormalization. Phys. Rev. Lett. 2007, 99, 220405. Neuscamman, E. Size Consistency Error in the Antisymmetric Geminal Power Wave Function can be Completely

75

76

77

78

79

80

81

82

83

84

85

86

87

88

89

90

91

Removed. Phys. Rev. Lett. 2012, 109, 203001. Mazziotti, D. Anti-Hermitian Contracted Schr¨ odinger Equation: Direct Determination of the Two-Electron Reduced Density Matrices of Many-Electron Molecules. Phys. Rev. Lett. 2006, 97, 143002. Parkhill, J. A.; Head-Gordon, M. A truncation hierarchy of coupled cluster models of strongly correlated systems based on perfect-pairing references: the singles+doubles models. J. Chem. Phys. 2010, 133, 124102. Parkhill, J. A.; Lawler, K.; Head-Gordon, M. The perfect quadruples model for electron correlation in a valence active space. J. Chem. Phys. 2009, 130, 084101. Small, D. W.; Head-Gordon, M. A fusion of the closedshell coupled cluster singles and doubles method and valence-bond theory for bond breaking. J. Chem. Phys. 2012, 137, 114103. Shao, Y.; Head-Gordon, M.; Krylov, A. I. The spinflip approach within time-dependent density functional theory: Theory and applications to diradicals. J. Chem. Phys. 2003, 118, 4807. Krylov, A. Spin-flip configuration interaction: an electronic structure model that is both variational and sizeconsistent. Chem. Phys. Lett. 2001, 350, 522–530. Evangelista, F. A. A driven similarity renormalization group approach to quantum many-body problems. J. Chem. Phys. 2014, 141, 054109. Zhang, T.; Evangelista, F. A. A Deterministic Projector Configuration Interaction Approach for the Ground State of Quantum Many-Body Systems. J. Chem. Theory Comput. 2016, 12, 4326–4337. Fosso-Tande, J.; Nguyen, T.-S.; Gidofalvi, G.; DePrince, A. E. Large-Scale Variational Two-Electron Reduced-Density-Matrix-Driven Complete Active Space Self-Consistent Field Methods. J. Chem. Theory Comput. 2016, 12, 2260–2271. DePrince, A. E. Variational optimization of the twoelectron reduced-density matrix under pure-state Nrepresentability conditions. J. Chem. Phys. 2016, 145, 164109. Parker, S. M.; Shiozaki, T. Communication: Active space decomposition with multiple sites: density matrix renormalization group algorithm. J. Chem. Phys. 2014, 141, 211102. Parker, S. M.; Seideman, T.; Ratner, M. A.; Shiozaki, T. Communication: Active-space decomposition for molecular dimers. J. Chem. Phys. 2013, 139, 021108. Parker, S. M.; Shiozaki, T. Quasidiabatic States from Active Space Decomposition. J. Chem. Theory Comput. 2014, 140715170536002. Li, S. Block-correlated coupled cluster theory: The general formulation and its application to the antiferromagnetic Heisenberg model. J. Chem. Phys. 2004, 120, 5017. Fang, T.; Shen, J.; Li, S. Block correlated coupled cluster method with a complete-active-space self-consistent-field reference function: The formula for general active spaces and its applications for multibond breaking systems. J. Chem. Phys. 2008, 128, 224107. Shen, J.; Li, S. Block correlated coupled cluster method with the complete active-space self-consistent-field reference function: Applications for low-lying electronic excited states. J. Chem. Phys. 2009, 131, year. Jim´enez-Hoyos, C. A.; Scuseria, G. E. Cluster-based mean-field and perturbative description of strongly correlated fermion systems: Application to the one- and

ACS Paragon Plus Environment

Page 13 of 14

Journal of Chemical Theory and Computation 13

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

92

93

94

95

96

97

98

99

100

101

102

103

104

105

106

107

108

two-dimensional Hubbard model. Phys. Rev. B 2015, 92, 085101. Al Hajj, M.; Malrieu, J.-P.; Guih´ery, N. Renormalized excitonic method in terms of block excitations: Application to spin lattices. Phys. Rev. B 2005, 72, 224412. Ma, Y.; Ma, H. Calculating excited states of molecular aggregates by the renormalized excitonic method. J. Phys. Chem. A 2013, 117, 3655–3665. Ma, Y.; Liu, Y.; Ma, H. A new fragment-based approach for calculating electronic excitation energies of large systems. J. Chem. Phys. 2012, 136, 024113. Zhang, H.; Malrieu, J.-P.; Ma, H.; Ma, J. Implementation of renormalized excitonic method at ab initio level. J. Comput. Chem. 2011, 33, 34–43. Morrison, A. F.; Herbert, J. M. Low-Scaling Quantum Chemistry Approach to Excited-State Properties via an ab Initio Exciton Model: Application to Excitation Energy Transfer in a Self-Assembled Nanotube. J. Phys. Chem. Lett. 2015, 6, 4390–4396. Morrison, A. F.; You, Z.-Q.; Herbert, J. M. Ab Initio Implementation of the Frenkel-Davydov Exciton Model: A Naturally Parallelizable Approach to Computing Collective Excitations in Crystals and Aggregates. J. Chem. Theory Comput. 2014, 10, 5366–76. Herbert, J. M.; Zhang, X.; Morrison, A. F.; Liu, J. Beyond Time-Dependent Density Functional Theory Using Only Single Excitations: Methods for Computational Studies of Excited States in Complex Systems. Accounts of Chemical Research 2016, 49, 931–941. Fang, T.; Shen, J.; Li, S. Block correlated coupled cluster method with a complete-active-space self-consistent-field reference function: The implementation for low-lying excited states. J. Chem. Phys. 2008, 129, 234106. Fang, T.; Li, S. Block correlated coupled cluster theory with a complete active-space self-consistent-field reference function: The formulation and test applications for single bond breaking. J. Chem. Phys. 2007, 127, 204108. Xu, E.; Li, S. Block correlated second order perturbation theory with a generalized valence bond reference function. J. Chem. Phys. 2013, 139, 174111. Bloch, C. Sur la th´eorie des perturbations des ´etats li´es. Nucl. Phys. 1958, 6, 329–347. des Cloizeaux, J. Extension d’une formule de Lagrange ` a des probl`emes de valeurs propres. Nucl. Phys. 1960, 20, 321–346. Al Hajj, M.; Malrieu, J.-P. On the spin gaps of conjugated hydrocarbon polymers. The Journal of Chemical Physics 2007, 127, 144902. Morrison, A. F.; Herbert, J. M. Analytic derivative couplings and first-principles exciton/phonon coupling constants for an ab initio Frenkel-Davydov exciton model: Theory, implementation, and application to compute triplet exciton mobility parameters for crystalline tetracene. The Journal of Chemical Physics 2017, 146, 224110. Davidson, E. R. The iterative calculation of a few of the lowest eigenvalues and corresponding eigenvectors of large real-symmetric matrices. J. Comput. Phys. 1975, 17, 87– 94. Pulay, P. ImprovedSCF convergence acceleration. J. Comput. Chem. 1982, 3, 556–560. Pulay, P. Convergence acceleration of iterative sequences. the case of scf iteration. Chem. Phys. Lett. 1980, 73, 393– 398.

109

110

111

112

Mayhall, N. J. From Model Hamiltonians to ab Initio Hamiltonians and Back Again: Using Single Excitation Quantum Chemistry Methods To Find Multiexciton States in Singlet Fission Materials. J. Chem. Theory Comput. 2016, 12, 4263–4273. Mayhall, N. J.; Head-Gordon, M. Computational Quantum Chemistry for Multiple-Site Heisenberg Spin Couplings Made Simple: Still Only One SpinFlip Required. J. Phys. Chem. Lett. 2015, 6, 1982–1988. Shao, Y. et al. Advances in molecular quantum chemistry contained in the Q-Chem 4 program package. Mol. Phys. 2014, 113, 184–215. Stein, C. J.; Reiher, M. Automated Selection of Active Orbital Spaces. J. Chem. Theory Comput. 2016.

ACS Paragon Plus Environment

Journal of Chemical Theory and Computation

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 ACS Paragon Plus Environment

Page 14 of 14