Search

We describe our implementation of fermionic tensor network contraction on arbitrary lattices within both a globally ordered and locally ordered formalism. We provide a pedagogical description of these two conventions as implemented for the quimb library. Using hyperoptimized approximate contraction strategies, we present benchmark fermionic projected entangled pair states simulations of finite Hubbard models defined on the three-dimensional diamond lattice and random regular graphs., Comment: 9 pages, 9 figures

The fate of thymine upon excitation by ultraviolet radiation has been the subject of intense debate over the past three decades. Today, it is widely believed that its ultrafast excited state decay stems from a radiationless transition from the bright ${\pi}{\pi}^*$ state to a dark $n{\pi}^*$ state. However, conflicting theoretical predictions have made the experimental data difficult to interpret. Here we simulate the ultrafast dynamics in thymine at the highest level of theory to date, performing wavepacket dynamics with a new coupled cluster method. Our simulation confirms an ultrafast ${\pi}{\pi}^*$ to $n{\pi}^*$ transition (${\tau} = 41 \pm 14$ fs). Furthermore, the predicted oxygen-edge X-ray absorption spectra agree quantitatively with the experimental results. Our simulation also predicts an as-yet uncharacterized photochemical pathway: a ${\pi}{\sigma}^*$ channel that leads to hydrogen dissociation at one of the two N-H bonds in thymine. Similar behavior has been identified in other heteroaromatic compounds, including adenine, and several authors have speculated that a similar pathway may exist in thymine. However, this was never confirmed theoretically or experimentally. This prediction calls for renewed efforts to experimentally identify or exclude the presence of this channel., Comment: 42 pages, 23 figures

Coupled cluster theory is one of the most accurate electronic structure methods for predicting ground and excited state chemistry. However, the presence of numerical artifacts at electronic degeneracies, such as complex energies, has made it difficult to apply it in nonadiabatic dynamics simulations. While it has already been shown that such numerical artifacts can be fully removed by using similarity constrained coupled cluster (SCC) theory [J. Phys. Chem. Lett. 2017, 8, 19, 4801-4807], simulating dynamics requires efficient implementations of gradients and nonadiabatic couplings. Here, we present an implementation of nuclear gradients and nonadiabatic derivative couplings at the similarity constrained coupled cluster singles and doubles (SCCSD) level of theory, thereby making possible nonadiabatic dynamics simulations using a coupled cluster theory that provides a correct description of conical intersections between excited states. We present a few numerical examples that show good agreement with literature values and discuss some limitations of the method., Comment: 40 pages, 3 figures, 2 tables

We present an efficient implementation of analytical non-adiabatic derivative coupling elements for the coupled cluster singles and doubles model. The derivative coupling elements are evaluated in a biorthonormal formulation in which the nuclear derivative acts on the right electronic state, where this state is biorthonormal with respect to the set of left states. This stands in contrast to earlier implementations based on normalized states and a gradient formula for the derivative coupling. As an illustration of the implementation, we determine a minimum energy conical intersection between the n{\pi}* and {\pi}{\pi}* states in the nucleobase thymine., Comment: 8 pages, 4 figures

We present an efficient implementation of ground and excited state CCSD gradients based on Cholesky-decomposed electron repulsion integrals. Cholesky decomposition, like density-fitting, is an inner projection method, and thus similar implementation schemes can be applied for both methods. One well-known advantage of inner projection methods, which we exploit in our implementation, is that one can avoid storing large V3O and V4 arrays by instead considering three-index intermediates. Furthermore, our implementation does not require the formation and storage of Cholesky vector derivatives. The new implementation is shown to perform well, with less than 10% of the time spent calculating the gradients in geometry optimizations. The computational time spent per optimization cycle are furthermore found to be significantly lower compared to other implementations based on an inner projection method. We illustrate the capabilities of the implementation by optimizing the geometry of the retinal molecule (C20H28O) at the CCSD/aug-cc-pVDZ level of theory., Comment: 12 pages, 1 figure

In coupled cluster methods, the electronic states are biorthonormal in the sense that the left states are orthonormal to the right states. Here we present an extension of this formalism to a left and right total molecular wave function. Starting from left and right Born-Huang expansions, we derive projected Schr\"odinger equations for the left and right nuclear wave functions. Observables may be extracted from the resulting wave function pair using standard expressions. The formalism is shown to be invariant under electronic basis transformations, such as normalization of the electronic states. Consequently, the nonadiabatic coupling elements can be expressed with biorthonormal wave functions. Calculating normalization factors that scale as full-CI is therefore not necessary, contrary to claims in the literature. For nuclear dynamics, we therefore need expressions for the vector and scalar couplings in the biorthonormal formalism. We derive these expressions using a Lagrangian formalism., Comment: 41 pages, 1 figure

We present efficient implementations of the multilevel CC2 (MLCC2) and multilevel CCSD (MLCCSD) models. As the system size increases, MLCC2 and MLCCSD exhibit the scaling of the lower level coupled cluster model. In order to treat large systems, we combine MLCC2 and MLCCSD with a reduced orbital space approach where the multilevel coupled cluster calculation is performed in a significantly truncated molecular orbital basis. The truncation scheme is based on the selection of an active region of the molecular system and the subsequent construction of localized Hartree-Fock orbitals. These orbitals are used in the multilevel coupled cluster calculation. The electron repulsion integrals are Cholesky decomposed using a screening protocol that guarantees accuracy in the truncated molecular orbital basis. The Cholesky factors are constructed directly in the truncated basis, ensuring low storage requirements. Even larger systems can be treated by using a multilevel Hartree-Fock reference. With the reduced space approach, we can handle systems with more than a thousand atoms. This is demonstrated for paranitroaniline in aqueous solution., Comment: 42 pages, 9 figures

We present an ab initio correlated approach to study molecules that interact strongly with quantum fields in an optical cavity. Quantum electrodynamics coupled cluster theory provides a non-perturbative description of cavity-induced effects in ground and excited states. Using this theory, we show how quantum fields can be used to manipulate charge transfer and photochemical properties of molecules. We propose a strategy to lift electronic degeneracies and induce modifications in the ground state potential energy surface close to a conical intersection., Comment: 18 pages, 12 figures

We introduce a multimodel approach to solve coupled cluster equations, employing a quasi Newton algorithm for the ground state and an Olsen algorithm for the excited states. In these algorithms, both of which can be viewed as Newton algorithms, the Jacobian matrix of a lower level coupled cluster model is used in Newton equations associated with the target model. Improvements in convergence then implies savings for sufficiently large molecular systems, since the computational cost of macroiterations scales more steeply with system size than the cost of microiterations. The multimodel approach is suitable when there is a lower level Jacobian matrix that is much more accurate than the zeroth order approximation. Applying the approach to the CC3 equations, using the CCSD approximation of the Jacobian, we show that the time spent to determine the ground and valence excited states can be significantly reduced. We also find improved convergence for core excited states, indicating that similar savings will be obtained with an explicit implementation of the core-valence separated CCSD Jacobian transformation., Comment: 22 pages and 3 figures

The eT program is an open source electronic structure package with emphasis on coupled cluster and multilevel methods. It includes efficient spin adapted implementations of ground and excited singlet states, as well as equation of motion oscillator strengths, for CCS, CC2, CCSD, and CC3. Furthermore, eT provides unique capabilities such as multilevel Hartree-Fock and multilevel CC2, real-time propagation for CCS and CCSD, and efficient CC3 oscillator strengths. With a coupled cluster code based on an efficient Cholesky decomposition algorithm for the electronic repulsion integrals, eT has similar advantages as codes using density fitting, but with strict error control. Here we present the main features of the program and demonstrate its performance through example calculations. Because of its availability, performance, and unique capabilities, we expect eT to become a valuable resource to the electronic structure community., Comment: 31 pages and 10 figures - supplementary information included in the uploaded files

We present an algorithm where only the Cholesky basis is determined in the decomposition procedure. This allows for improved screening and a partitioned matrix decomposition scheme, both of which significantly reduce memory usage and computational cost. After the basis has been determined, an inner projection technique is used to construct the Cholesky vectors. The algorithm extends the application range of the methodology and is well suited for multilevel methods. We apply the algorithm to systems with up to 80000 atomic orbitals., Comment: 17 pages, 6 figures, submitted to The Journal of Chemical Physics

We derive the crossing conditions at conical intersections between electronic states in coupled cluster theory, and show that if the coupled cluster Jacobian matrix is nondefective, two (three) independent conditions are correctly placed on the nuclear degrees of freedom for an inherently real (complex) Hamiltonian. Calculations using coupled cluster theory on an $2 {^{1}}A' / 3 {^{1}}A'$ conical intersection in hypofluorous acid illustrate the nonphysical artifacts associated with defects at accidental same-symmetry intersections. In particular, the observed intersection seam is folded about a space of the correct dimensionality, indicating that minor modifications to the theory are required for it to provide a correct description of conical intersections in general. We find that an accidental symmetry allowed $1 {^{1}}A" / 2 {^{1}}A"$ intersection in hydrogen sulfide is properly described, showing no artifacts as well as linearity of the energy gap to first order in the branching plane., Comment: 9 pages and 4 figures

The motion of electrons and nuclei in photochemical events often involve conical intersections, degeneracies between electronic states. They serve as funnels for nuclear relaxation - on the femtosecond scale - in processes where the electrons and nuclei couple nonadiabatically. Accurate ab initio quantum chemical models are essential for interpreting experimental measurements of such phenomena. In this paper we resolve a long-standing problem in coupled cluster theory, presenting the first formulation of the theory that correctly describes conical intersections between excited electronic states of the same symmetry. This new development demonstrates that the highly accurate coupled cluster theory can be applied to describe dynamics on excited electronic states involving conical intersections., Comment: 8 pages and 3 figures and including supporting information (with corrections and improved notation)