Implementation and applications of density-fitted symmetry-adapted perturbation theory

Noncovalent interactions play a vital role throughout much of chemistry. The understanding and characterization of these interactions is an area where theoretical chemistry can provide unique insight. While many methods have been developed to study noncovalent interactions, symmetry-adapted perturbation theory (SAPT) stands out as one of the most robust. In addition to providing energetic information about an interaction, it provides insight into the underlying physics of the interaction by decomposing the energy into electrostatics, exchange, induction and dispersion. Therefore, SAPT is capable of not only answering questions about how strongly a complex is bound, but also why it is bound. This proves to be an invaluable tool for the understanding of noncovalent interactions in complex systems. The wavefunction-based formulation of SAPT can provide qualitative results for large systems as well as quantitative results for smaller systems. In order to extend the applicability of this method, approximations to the two-electron integrals must be introduced. At low-order, the introduction of density fitting approximations allows SAPT computations to be performed on systems with up to 220 atoms and 2850 basis functions. Higher-orders of SAPT, which boasts accuracy rivaling the best theoretical methods, can be applied to systems with over 40 atoms. Higher-order SAPT also benefits from approximations that attempt to truncate unneccesary unoccupied orbitals.PhDCommittee Chair: Sherrill, C. David; Committee Member: Bongiorno, Angelo; Committee Member: Brown, Kenneth R.; Committee Member: Harvey, Stephen C.; Committee Member: Liotta, Charles L.; Committee Member: Whetten, Robert L

Tensor hypercontraction: A universal technique for the resolution of matrix elements of local, finite-range NN-body potentials in many-body quantum problems

Configuration-space matrix elements of N-body potentials arise naturally and ubiquitously in the Ritz-Galerkin solution of many-body quantum problems. For the common specialization of local, finite-range potentials, we develop the eXact Tensor HyperContraction (X-THC) method, which provides a quantized renormalization of the coordinate-space form of the N-body potential, allowing for a highly separable tensor factorization of the configuration-space matrix elements. This representation allows for substantial computational savings in chemical, atomic, and nuclear physics simulations, particularly with respect to difficult "exchange-like" contractions.Comment: Third version of the manuscript after referee's comments. In press in PRL. Main text: 4 pages, 2 figures, 1 table; Supplemental material (also included): 14 pages, 2 figures, 2 table

Improvement of the coupled-cluster singles and doubles method via scaling same- and opposite-spin components of the double excitation correlation energy

© 2008 American Institute of Physics. The electronic version of this article is the complete one and can be found at: http://dx.doi.org/10.1063/1.2883974DOI: 10.1063/1.2883974There has been much interest in cost-free improvements to second-order Moller-Plesset perturbation theory (MP2) via scaling the same- and opposite-spin components of the correlation energy (spin-component scaled MP2). By scaling the same- and opposite-spin components of the double excitation correlation energy from the coupled-cluster of single and double excitations (CCSD) method, similar improvements can be achieved. Optimized for a set of 48 reaction energies, scaling factors were determined to be 1.13 and 1.27 for the same- and opposite-spin components, respectively. Preliminary results suggest that the spin-component scaled CCSD (SCS-CCSD) method will outperform all MP2 type methods considered for describing intermolecular interactions. Potential energy curves computed with the SCS-CCSD method for the sandwich benzene dimer and methane dimer reproduce the benchmark CCSD(T) potential curves with errors of only a few hundredths of 1 kcal mol⁻¹ for the minima. The performance of the SCS-CCSD method suggests that it is a reliable, lower cost alternative to the CCSD(T) method

Ultrafast internal conversion and photochromism in gas-phase salicylideneaniline

Salicylideneaniline (SA) is an archetypal system for excited-state intramolecular proton transfer (ESIPT) in non-planar systems. Multiple channels for relaxation involving both the keto and enol forms have been proposed after excitation to S1 with near-UV light. Here, we present transient absorption measurements of hot gas-phase SA, jet-cooled SA, and SA in Ar clusters using cavity-enhanced transient absorption spectroscopy (CE-TAS). Assignment of the spectra is aided by simulated TAS spectra, computed by applying time-dependent complete active space configuration interaction (TD-CASCI) to structures drawn from nonadiabatic molecular dynamics simulations. We find prompt ESIPT in all conditions followed by the rapid generation of the trans keto metastable photochrome state and fluorescent keto state in parallel. Increasing the internal energy increases the photochrome yield and decreases the fluorescent yield and fluorescent state lifetime observed in TAS. In Ar clusters, internal conversion of SA is severely hindered, but the photochrome yield is unchanged. Taken together, these results are consistent with the photochrome being produced via the vibrationally excited keto population after ESIPT

Ultrafast internal conversion and photochromism in gas-phase salicylideneaniline

Salicylidenaniline (SA) is an archetypal system for excited-state intramolecular proton transfer (ESIPT) in non-planar systems. Multiple channels for relaxation involving both the keto and enol forms have been proposed after excitation to S$_1$ with near-UV light. Here we present transient absorption measurements of hot gas-phase SA, jet-cooled SA, and SA in Ar clusters using cavity-enhanced transient absorption spectroscopy (CE-TAS). Assignment of the spectra is aided by simulated TAS spectra, computed by applying time-dependent complete active space configuration interaction (TD-CASCI) to structures drawn from nonadiabatic molecular dynamics simulations. We find prompt ESIPT in all conditions followed by the rapid parallel generation of the trans-keto metastable photochrome state and fluorescent keto state in parallel. Increasing the internal energy increases the photochrome yield and decreases the fluorescent yield and fluorescent state lifetime observed in TAS. In Ar clusters, internal conversion of SA is severely hindered but the photochrome yield is unchanged. Taken together, these results are consistent with the photochrome being produced via the vibrationally excited keto population after ESIPT

Efficient Quantum Analytic Nuclear Gradients with Double Factorization

Efficient representations of the Hamiltonian such as double factorization drastically reduce circuit depth or number of repetitions in error corrected and noisy intermediate scale quantum (NISQ) algorithms for chemistry. We report a Lagrangian-based approach for evaluating relaxed one- and two-particle reduced density matrices from double factorized Hamiltonians, unlocking efficiency improvements in computing the nuclear gradient and related derivative properties. We demonstrate the accuracy and feasibility of our Lagrangian-based approach to recover all off-diagonal density matrix elements in classically-simulated examples with up to 327 quantum and 18470 total atoms in QM/MM simulations, with modest-sized quantum active spaces. We show this in the context of the variational quantum eigensolver (VQE) in case studies such as transition state optimization, ab initio molecular dynamics simulation and energy minimization of large molecular systems.Comment: 22 pages, 5 figure

Atomistic non-adiabatic dynamics of the LH2 complex with a GPU-accelerated ab initio exciton model

We present GPU-accelerated ab initio molecular dynamics simulations of nonadiabatic dynamics in the LH2 complex in full atomistic detail.</p

Psi4

Psi4 is an ab initio electronic structure program providing methods such as Hartree-Fock, density functional theory, configuration interaction, and coupled-cluster theory. The 1.1 release represents a major update meant to automate complex tasks, such as geometry optimization using complete-basis-set extrapolation or focal-point methods. Conversion of the top-level code to a Python module means that Psi4 can now be used in complex workflows alongside other Python tools. Several new features have been added with the aid of libraries providing easy access to techniques such as density fitting, Cholesky decomposition, and Laplace denominators. The build system has been completely rewritten to simplify interoperability with independent, reusable software components for quantum chemistry. Finally, a wide range of new theoretical methods and analyses have been added to the code base, including functional-group and open-shell symmetry adapted perturbation theory, density-fitted coupled cluster with frozen natural orbitals, orbital-optimized perturbation and coupled-cluster methods (e.g., OO-MP2 and OO-LCCD), density-fitted multiconfigurational self-consistent field, density cumulant functional theory, algebraic-diagrammatic construction excited states, improvements to the geometry optimizer, and the "X2C" approach to relativistic corrections, among many other improvements