1
|
Elayan IA, Rib L, A Mendes R, Brown A. Beyond Explored Functionals: A Computational Journey of Two-Photon Absorption. J Chem Theory Comput 2024; 20:3879-3893. [PMID: 38648613 DOI: 10.1021/acs.jctc.4c00133] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/25/2024]
Abstract
We present a thorough investigation into the efficacy of 19 density functional theory (DFT) functionals, relative to RI-CC2 results, for computing two-photon absorption (2PA) cross sections (σ2PA) and key dipole moments (|μ00|, |μ11|, |Δμ|, |μ01|) for a series of coumarin dyes in the gas-phase. The functionals include different categories, including local density approximation (LDA), generalized gradient approximation (GGA), hybrid-GGA (H-GGA), range-separated hybrid-GGA (RSH-GGA), meta-GGA (M-GGA), and hybrid M-GGA (HM-GGA), with 14 of them being subjected to analysis for the first time with respect to predicting σ2PA values. Analysis reveals that functionals integrating both short-range (SR) and long-range (LR) corrections, particularly those within the RSH-GGA and HM-GGA classes, outperform the others. Furthermore, the range-separation approach was found more impactful compared to the varying percentages of Hartree-Fock exchange (HF Ex) within different functionals. The functionals traditionally recommended for 2PA do not appear among the top 9 in our study, which is particularly interesting, as these top-performing functionals have not been previously investigated in this context. This list is dominated by M11, QTP variants, ωB97X, ωB97X-V, and M06-2X, surpassing the performance of other functionals, including the commonly used CAM-B3LYP.
Collapse
Affiliation(s)
- Ismael A Elayan
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| | - Laura Rib
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| | - Rodrigo A Mendes
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
- Departamento de Química, Universidade Federal de Mato Grosso, Cuiabá, Mato Grosso 78060-900, Brazil
| | - Alex Brown
- Department of Chemistry, University of Alberta, Edmonton, Alberta T6G 2G2, Canada
| |
Collapse
|
2
|
Mashkovtsev D, Orimoto Y, Aoki Y. Fast and Accurate Calculation of the UV-Vis Spectrum with the Modified Local Excitation Approximation. J Chem Theory Comput 2023; 19:5548-5562. [PMID: 37471461 DOI: 10.1021/acs.jctc.3c00266] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/22/2023]
Abstract
The local excitation approximation (LEA), a method for the calculation of electronic excitations localized in a specific region of a molecule, has been modified with new approaches to enhance the accuracy of the original method. The primary concept behind LEA involves isolating the region of interest as a submolecule from the full molecule using a localization method, followed by calculating electronic excitations solely within this submolecule. In this study, we examined approaches that improve the accuracy in describing the region of interest, particularly its molecular orbital energies. Additionally, the localization method was extended with a new projection technique to accelerate calculations. These approaches were studied in time-dependent density functional theory (TDDFT) calculations applied to four testing systems with a chromophore as the region of interest: two basic linear molecules, acrolein surrounded by 24 water molecules, and a model of a green fluorescent protein. For all studied systems, the results of TDDFT calculations combined with LEA exhibited near-zero error when groups of atoms adjacent to the chromophore were explicitly included in the submolecule. This was achieved with at least a quadratic speedup of the calculation time as a function of the submolecule size.
Collapse
Affiliation(s)
- Denis Mashkovtsev
- Department of Molecular and Material Sciences, Interdisciplinary Graduate School of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| | - Yuuichi Orimoto
- Department of Material Sciences, Faculty of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| | - Yuriko Aoki
- Department of Material Sciences, Faculty of Engineering Sciences, Kyushu University, 6-1 Kasuga-Park, Fukuoka 816-8580, Japan
| |
Collapse
|
3
|
Corzo HH, Hillers-Bendtsen AE, Barnes A, Zamani AY, Pawłowski F, Olsen J, Jørgensen P, Mikkelsen KV, Bykov D. Corrigendum: Coupled cluster theory on modern heterogeneous supercomputers. Front Chem 2023; 11:1256510. [PMID: 37654900 PMCID: PMC10466216 DOI: 10.3389/fchem.2023.1256510] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2023] [Accepted: 07/11/2023] [Indexed: 09/02/2023] Open
Abstract
[This corrects the article DOI: 10.3389/fchem.2023.1154526.].
Collapse
Affiliation(s)
| | | | | | - Abdulrahman Y. Zamani
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, CA, United States
| | - Filip Pawłowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, AL, United States
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University, Aarhus, Denmark
| | - Poul Jørgensen
- Department of Chemistry, Aarhus University, Aarhus, Denmark
| | - Kurt V. Mikkelsen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Dmytro Bykov
- Oak Ridge National Laboratory, Oak Ridge, TN, United States
| |
Collapse
|
4
|
Corzo HH, Hillers-Bendtsen AE, Barnes A, Zamani AY, Pawłowski F, Olsen J, Jørgensen P, Mikkelsen KV, Bykov D. Coupled cluster theory on modern heterogeneous supercomputers. Front Chem 2023; 11:1154526. [PMID: 37388945 PMCID: PMC10303140 DOI: 10.3389/fchem.2023.1154526] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2023] [Accepted: 05/11/2023] [Indexed: 07/01/2023] Open
Abstract
This study examines the computational challenges in elucidating intricate chemical systems, particularly through ab-initio methodologies. This work highlights the Divide-Expand-Consolidate (DEC) approach for coupled cluster (CC) theory-a linear-scaling, massively parallel framework-as a viable solution. Detailed scrutiny of the DEC framework reveals its extensive applicability for large chemical systems, yet it also acknowledges inherent limitations. To mitigate these constraints, the cluster perturbation theory is presented as an effective remedy. Attention is then directed towards the CPS (D-3) model, explicitly derived from a CC singles parent and a doubles auxiliary excitation space, for computing excitation energies. The reviewed new algorithms for the CPS (D-3) method efficiently capitalize on multiple nodes and graphical processing units, expediting heavy tensor contractions. As a result, CPS (D-3) emerges as a scalable, rapid, and precise solution for computing molecular properties in large molecular systems, marking it an efficient contender to conventional CC models.
Collapse
Affiliation(s)
| | | | | | - Abdulrahman Y. Zamani
- Department of Chemistry and Biochemistry and Center for Chemical Computation and Theory, University of California, Merced, CA, United States
| | - Filip Pawłowski
- Department of Chemistry and Biochemistry, Auburn University, Auburn, AL, United States
| | - Jeppe Olsen
- Department of Chemistry, Aarhus University, Aarhus, Denmark
| | - Poul Jørgensen
- Department of Chemistry, Aarhus University, Aarhus, Denmark
| | - Kurt V. Mikkelsen
- Department of Chemistry, University of Copenhagen, Copenhagen, Denmark
| | - Dmytro Bykov
- Oak Ridge National Laboratory, Oak Ridge, TN, United States
| |
Collapse
|
5
|
Sacchetta F, Graf D, Laqua H, Ambroise MA, Kussmann J, Dreuw A, Ochsenfeld C. An effective sub-quadratic scaling atomic-orbital reformulation of the scaled opposite-spin RI-CC2 ground-state model using Cholesky-decomposed densities and an attenuated Coulomb metric. J Chem Phys 2022; 157:104104. [DOI: 10.1063/5.0098719] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An atomic-orbital reformulation of the Laplace-transformed scaled opposite-spin (SOS) coupled cluster singles and doubles (CC2) model within the resolution of the identity (RI) approximation (SOS-RI-CC2) is presented that extends its applicability to molecules with several hundreds of atoms and triple-zeta basis sets. We exploit sparse linear algebra and an attenuated Coulomb metric to decrease the disk space demands and the computational efforts. In this way, an effective sub-quadratic computational scaling is achieved with our ω-SOS-CDD-RI-CC2 model. Moreover, Cholesky decomposition of the ground-state one-electron density matrix reduces the prefactor, allowing for an early crossover with the molecular orbital formulation. The accuracy and performance of the presented method are investigated for various molecular systems.
Collapse
Affiliation(s)
- F. Sacchetta
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - D. Graf
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - H. Laqua
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - M. A. Ambroise
- Chair of Theoretical and Computational Chemistry, Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg, Germany
| | - J. Kussmann
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - A. Dreuw
- Chair of Theoretical and Computational Chemistry, Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg, Germany
| | - C. Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| |
Collapse
|
6
|
D'Cunha R, Crawford TD. Applications of a perturbation-aware local correlation method to coupled cluster linear response properties. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2112627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/15/2022]
Affiliation(s)
- Ruhee D'Cunha
- Department of Chemistry, Virginia Tech, Blacksburg, VA, USA
| | - T. Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, VA, USA
- Molecular Sciences Software Institute, Blacksburg, VA, USA
| |
Collapse
|
7
|
Pant R, Ranga S, Bachhar A, Dutta AK. Pair Natural Orbital Equation-of-Motion Coupled-Cluster Method for Core Binding Energies: Theory, Implementation, and Benchmark. J Chem Theory Comput 2022; 18:4660-4673. [PMID: 35786933 DOI: 10.1021/acs.jctc.2c00165] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the theory and implementation of a lower scaling core-valence separated equation-of-motion coupled-cluster approach based on domain-based local pair natural orbitals for core binding energies. The accuracy of the new method has been compared with that of the standard equation-of-motion coupled-cluster method and experimentally measured results. The use of pair natural orbitals significantly reduces the computation cost and can be applied to large molecules.
Collapse
Affiliation(s)
- Rakesh Pant
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Santosh Ranga
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Arnab Bachhar
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
8
|
Haldar S, Mukhopadhyay T, Dutta AK. A similarity transformed second-order approximate coupled cluster method for the excited states: Theory, implementation, and benchmark. J Chem Phys 2022; 156:014110. [PMID: 34998329 DOI: 10.1063/5.0064889] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present a novel and cost-effective approach of using a second similarity transformation of the Hamiltonian to include the missing higher-order terms in the second-order approximate coupled cluster singles and doubles (CC2) model. The performance of the newly developed ST-EOM-CC2 model has been investigated for the calculation of excitation energies of valence, Rydberg, and charge-transfer excited states. The method shows significant improvement in the excitation energies of Rydberg and charge-transfer excited states as compared to the conventional CC2 method while retaining the good performance of the latter for the valence excited state. This method retains the charge-transfer separability of the charge-transfer excited states, which is a significant advantage over the traditional CC2 method. A second order many-body perturbation theory variant of the new method is also proposed.
Collapse
Affiliation(s)
- Soumi Haldar
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Tamoghna Mukhopadhyay
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
9
|
Ranga S, Dutta AK. A Core-Valence Separated Similarity Transformed EOM-CCSD Method for Core-Excitation Spectra. J Chem Theory Comput 2021; 17:7428-7446. [PMID: 34814683 DOI: 10.1021/acs.jctc.1c00402] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the theory and implementation of a core-valence separated similarity transformed EOM-CCSD (STEOM-CCSD) method for K-edge core excitation spectra. The method can select an appropriate active space using CIS natural orbitals and near "black box" to use. The second similarity transformed Hamiltonian is diagonalized in the space of single excitation. Therefore, the final diagonalization step is free from the convergence problem arising due to the coupling of the core-excited states with the continuum of doubly excited states. Convergence trouble can appear for the preceding core-ionized state calculation in STEOM-CCSD. A core-valence separation (CVS) scheme compatible with the natural orbital based active space selection (CVS-STEOM-CCSD-NO) is implemented to overcome the problem. The CVS-STEOM-CCSD-NO has a similar accuracy to that of the standard CVS-EOM-CCSD method but comes with a lower computational cost. The modification required in the CVS scheme to make use of the CIS natural orbital is highlighted. The suitability of the CVS-STEOM-CCSD-NO method for chemical application is demonstrated by simulating the K-edge spectra of glycine and thymine.
Collapse
Affiliation(s)
- Santosh Ranga
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
10
|
Haldar S, Dutta AK. An efficient Fock space multi-reference coupled cluster method based on natural orbitals: Theory, implementation, and benchmark. J Chem Phys 2021; 155:014105. [PMID: 34241374 DOI: 10.1063/5.0054171] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a natural orbital-based implementation of the intermediate Hamiltonian Fock space coupled-cluster method for the (1, 1) sector of Fock space. The use of natural orbitals significantly reduces the computational cost and can automatically choose an appropriate set of active orbitals. The new method retains the charge transfer separability of the original intermediate Hamiltonian Fock space coupled-cluster method and gives excellent performance for valence, Rydberg, and charge-transfer excited states. It offers significant computational advantages over the popular equation of motion coupled cluster method for excited states dominated by single excitations.
Collapse
Affiliation(s)
- Soumi Haldar
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
11
|
D'Cunha R, Crawford TD. PNO++: Perturbed Pair Natural Orbitals for Coupled Cluster Linear Response Theory. J Chem Theory Comput 2021; 17:290-301. [PMID: 33351627 DOI: 10.1021/acs.jctc.0c01086] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Reduced-scaling methods are needed to make accurate and systematically improvable coupled cluster linear response methods for the calculation of molecular properties tractable for large molecules. In this paper, we examine the perturbed pair natural orbital-based PNO++ approach that creates an orbital space optimized for response properties derived from a lower-cost field-perturbed density matrix. We analyze truncation errors in correlation energies, dynamic polarizabilities, and specific rotations from a coupled cluster singles and doubles (CCSD) reference. We find that incorporating a fixed number of orbitals from the pair natural orbital (PNO) space into the PNO++ method-a new method presented here, the "combined PNO++" approach-recovers accuracy in the CCSD correlation energy while preserving the well-behaved convergence behavior of the PNO++ method for linear response properties.
Collapse
Affiliation(s)
- Ruhee D'Cunha
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - T Daniel Crawford
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States.,Molecular Sciences Software Institute, 1880 Pratt Drive, Suite 1100, Blacksburg, Virginia 24060, United States
| |
Collapse
|
12
|
Folkestad SD, Kjønstad EF, Goletto L, Koch H. Multilevel CC2 and CCSD in Reduced Orbital Spaces: Electronic Excitations in Large Molecular Systems. J Chem Theory Comput 2021; 17:714-726. [PMID: 33417769 PMCID: PMC8016205 DOI: 10.1021/acs.jctc.0c00590] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
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. To treat large systems, we combine MLCC2 and MLCCSD with a reduced-space approach in which 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 and reduces computational cost. The Cholesky factors are constructed directly in the truncated basis, ensuring low storage requirements. Systems for which Hartree-Fock is too expensive 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.
Collapse
Affiliation(s)
- Sarai Dery Folkestad
- Department of Chemistry, Norwegian University of Science and Technology, N-7491, Trondheim, Norway
| | - Eirik F Kjønstad
- Department of Chemistry, Norwegian University of Science and Technology, N-7491, Trondheim, Norway
| | - Linda Goletto
- Department of Chemistry, Norwegian University of Science and Technology, N-7491, Trondheim, Norway
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, N-7491, Trondheim, Norway.,Scuola Normale Superiore, Piazza dei Cavaleri 7, Pisa, 56126, Italy
| |
Collapse
|
13
|
Macetti G, Genoni A. Quantum Mechanics/Extremely Localized Molecular Orbital Embedding Strategy for Excited States: Coupling to Time-Dependent Density Functional Theory and Equation-of-Motion Coupled Cluster. J Chem Theory Comput 2020; 16:7490-7506. [PMID: 33241930 DOI: 10.1021/acs.jctc.0c00956] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The QM/ELMO (quantum mechanics/extremely localized molecular orbital) method is a recently developed embedding technique in which the most important region of the system under examination is treated at fully quantum mechanical level, while the rest is described by means of transferred and frozen extremely localized molecular orbitals. In this paper, we propose the first application of the QM/ELMO approach to the investigation of excited states, and, in particular, we present the coupling of the QM/ELMO philosophy with Time-Dependent Density Functional Theory (TDDFT) and Equation-of-Motion Coupled Cluster with single and double substitutions (EOM-CCSD). The proposed TDDFT/ELMO and EOM-CCSD/ELMO strategies underwent a series of preliminary tests that were already considered for the validation of other embedding methods for excited states. The obtained results showed that the novel techniques allow the accurate description of localized excitations in large systems by only including a relatively small number of atoms in the region treated at fully quantum chemical level. Furthermore, for TDDFT/ELMO, it was also observed that (i) the method enables to avoid the presence of artificial low-lying charge-transfer states that may affect traditional TDDFT calculations, even using functionals that do not take into account long-range corrections, and (ii) the novel approach can be also successfully exploited to investigate local electronic transitions in quite large systems (e.g., reduced model of the Green Fluorescent Protein), and the accuracy of the results can be improved by including a sufficient number of chemically crucial fragments/residues in the quantum mechanical region. Finally, concerning EOM-CCSD/ELMO, it was also seen that, despite the quite crude approximation of an embedding potential given by frozen extremely localized molecular orbitals, the new strategy is able to satisfactorily account for the effects of the environment. This work paves the way to further extensions of the QM/ELMO philosophy for the study of local excitations in extended systems, suggesting the coupling of the QM/ELMO approach with other quantum chemical strategies for excited states, from the simplest ΔSCF techniques to the most advanced and computationally expensive multireferences methods.
Collapse
Affiliation(s)
- Giovanni Macetti
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| | - Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, F-57078 Metz, France
| |
Collapse
|
14
|
Folkestad SD, Koch H. Equation-of-Motion MLCCSD and CCSD-in-HF Oscillator Strengths and Their Application to Core Excitations. J Chem Theory Comput 2020; 16:6869-6879. [PMID: 32955866 PMCID: PMC8011930 DOI: 10.1021/acs.jctc.0c00707] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
We present an implementation of equation-of-motion oscillator strengths for the multilevel CCSD (MLCCSD) model where CCS is used as the lower level method (CCS/CCSD). In this model, the double excitations of the cluster operator are restricted to an active orbital space, whereas the single excitations are unrestricted. Calculated nitrogen K-edge spectra of adenosine, adenosine triphosphate (ATP), and an ATP-water system are used to demonstrate the performance of the model. Projected atomic orbitals (PAOs) are used to partition the virtual space into active and inactive orbital sets. Cholesky decomposition of the Hartree-Fock density is used to partition the occupied orbitals. This Cholesky-PAO partitioning is cheap, scaling as O(N3), and is suitable for the calculation of core excitations, which are localized in character. By restricting the single excitations of the cluster operator to the active space, as well as the double excitations, the CCSD-in-HF model is obtained. A comparison of the two models-MLCCSD and CCSD-in-HF-is presented for the core excitation spectra of the adenosine and ATP systems.
Collapse
Affiliation(s)
- Sarai Dery Folkestad
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim N-7491, Norway
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, Trondheim N-7491, Norway.,Scuola Normale Superiore, Piazza dei Cavaleri 7, Pisa 56126, Italy
| |
Collapse
|
15
|
Ghosh S, Kumar Dutta A, de Souza B, Berraud-Pache R, Izsák R. A new density for transition properties within the similarity transformed equation of motion approach. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1818858] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Soumen Ghosh
- Max-Planck-Institut für Kohlenforschung, Mülheim an der Ruhr, Germany
| | | | | | | | - Róbert Izsák
- Max-Planck-Institut für Kohlenforschung, Mülheim an der Ruhr, Germany
| |
Collapse
|
16
|
Frank MS, Schmitz G, Hättig C. Implementation of the iterative triples model CC3 for excitation energies using pair natural orbitals and Laplace transformation techniques. J Chem Phys 2020; 153:034109. [PMID: 32716174 DOI: 10.1063/5.0012597] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a pair natural orbital (PNO)-based implementation of CC3 excitation energies, which extends our previously published state-specific PNO ansatz for the solution of the excited state eigenvalue problem to methods including connected triple excitations. A thorough analysis of the equations for the excited state triples amplitudes is presented from which we derive a suitable state-specific triple natural orbital basis for the excited state triples amplitudes, which performs equally well for local and non-local excitations. The accuracy of the implementation is evaluated using a large and diverse test set. We find that for states with small contributions from double excitations, a T0 approximation to PNO-CC3 yields accurate results with a mean absolute error (MAE) for TPNO = 10-7 in the range of 0.02 eV. However, for states with larger double excitation contributions, the T0 approximation is found to yield significantly less accurate results, while the Laplace-transformed variant of PNO-CC3 shows a uniform accuracy for singly and doubly excited states (MAE and maximum error of 0.01 eV and 0.07 eV for TPNO = 10-7, respectively). Finally, we apply PNO-CC3 to the calculation of the first excited state of berenil at a S1 minimum geometry, which is shown to be close to a conical intersection. This calculation in the aug-cc-pVTZ basis set (more than 1300 basis functions) is the largest calculation ever performed with CC3 on excitation energies.
Collapse
Affiliation(s)
- Marius S Frank
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44801 Bochum, Germany
| | - Gunnar Schmitz
- Deparment of Chemistry, Aarhus University, DK-8000 Aarhus C, Denmark
| | - Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44801 Bochum, Germany
| |
Collapse
|
17
|
Olsen JMH, Reine S, Vahtras O, Kjellgren E, Reinholdt P, Hjorth Dundas KO, Li X, Cukras J, Ringholm M, Hedegård ED, Di Remigio R, List NH, Faber R, Cabral Tenorio BN, Bast R, Pedersen TB, Rinkevicius Z, Sauer SPA, Mikkelsen KV, Kongsted J, Coriani S, Ruud K, Helgaker T, Jensen HJA, Norman P. Dalton Project: A Python platform for molecular- and electronic-structure simulations of complex systems. J Chem Phys 2020; 152:214115. [PMID: 32505165 DOI: 10.1063/1.5144298] [Citation(s) in RCA: 40] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023] Open
Abstract
The Dalton Project provides a uniform platform access to the underlying full-fledged quantum chemistry codes Dalton and LSDalton as well as the PyFraME package for automatized fragmentation and parameterization of complex molecular environments. The platform is written in Python and defines a means for library communication and interaction. Intermediate data such as integrals are exposed to the platform and made accessible to the user in the form of NumPy arrays, and the resulting data are extracted, analyzed, and visualized. Complex computational protocols that may, for instance, arise due to a need for environment fragmentation and configuration-space sampling of biochemical systems are readily assisted by the platform. The platform is designed to host additional software libraries and will serve as a hub for future modular software development efforts in the distributed Dalton community.
Collapse
Affiliation(s)
- Jógvan Magnus Haugaard Olsen
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Simen Reine
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, N-0315 Oslo, Norway
| | - Olav Vahtras
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Erik Kjellgren
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Karen Oda Hjorth Dundas
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Xin Li
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Janusz Cukras
- Department of Chemistry, University of Warsaw, 02-093 Warsaw, Poland
| | - Magnus Ringholm
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Erik D Hedegård
- Division of Theoretical Chemistry, Lund University, SE-223 62 Lund, Sweden
| | - Roberto Di Remigio
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Nanna H List
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Rasmus Faber
- DTU Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | | | - Radovan Bast
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Thomas Bondo Pedersen
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, N-0315 Oslo, Norway
| | - Zilvinas Rinkevicius
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Stephan P A Sauer
- Department of Chemistry, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark
| | - Kurt V Mikkelsen
- Department of Chemistry, University of Copenhagen, DK-2100 Copenhagen Ø, Denmark
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Kenneth Ruud
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Trygve Helgaker
- Department of Chemistry, Hylleraas Centre for Quantum Molecular Sciences, University of Oslo, N-0315 Oslo, Norway
| | - Hans Jørgen Aa Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Patrick Norman
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| |
Collapse
|
18
|
Haldar S, Dutta AK. A Multilayer Approach to the Equation of Motion Coupled-Cluster Method for the Electron Affinity. J Phys Chem A 2020; 124:3947-3962. [DOI: 10.1021/acs.jpca.0c01793] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Affiliation(s)
- Soumi Haldar
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Achintya Kumar Dutta
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
19
|
Wen X, Graham DS, Chulhai DV, Goodpaster JD. Absolutely Localized Projection-Based Embedding for Excited States. J Chem Theory Comput 2019; 16:385-398. [DOI: 10.1021/acs.jctc.9b00959] [Citation(s) in RCA: 31] [Impact Index Per Article: 6.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Affiliation(s)
- Xuelan Wen
- Department of Chemistry, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455, United States
| | - Daniel S. Graham
- Department of Chemistry, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455, United States
| | - Dhabih V. Chulhai
- Department of Chemistry, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455, United States
| | - Jason D. Goodpaster
- Department of Chemistry, University of Minnesota, 207 Pleasant Street Southeast, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
20
|
Folkestad SD, Koch H. Multilevel CC2 and CCSD Methods with Correlated Natural Transition Orbitals. J Chem Theory Comput 2019; 16:179-189. [DOI: 10.1021/acs.jctc.9b00701] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sarai Dery Folkestad
- Department of Chemistry, Norwegian University of Science and Technology, N-7491 Trondheim, Norway
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, N-7491 Trondheim, Norway
- Scuola Normale Superiore, Piazza dei Cavaleri 7, 56126 Pisa, Italy
| |
Collapse
|
21
|
Mester D, Nagy PR, Kállay M. Reduced-Scaling Correlation Methods for the Excited States of Large Molecules: Implementation and Benchmarks for the Second-Order Algebraic-Diagrammatic Construction Approach. J Chem Theory Comput 2019; 15:6111-6126. [DOI: 10.1021/acs.jctc.9b00735] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, H-1521 Budapest, P.O. Box 91, Hungary
| | - Péter R. Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, H-1521 Budapest, P.O. Box 91, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, H-1521 Budapest, P.O. Box 91, Hungary
| |
Collapse
|
22
|
Izsák R. Single‐reference coupled cluster methods for computing excitation energies in large molecules: The efficiency and accuracy of approximations. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1445] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Róbert Izsák
- Max‐Planck‐Institut für Kohlenforschung Mülheim an der Ruhr Germany
| |
Collapse
|
23
|
Liu J, Sun H, Glover WJ, He X. Prediction of Excited-State Properties of Oligoacene Crystals Using Fragment-Based Quantum Mechanical Method. J Phys Chem A 2019; 123:5407-5417. [PMID: 31187994 DOI: 10.1021/acs.jpca.8b12552] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
A fundamental understanding of the excited-state properties of molecular crystals is of central importance for their optoelectronics applications. In this study, we developed the electrostatically embedded generalized molecular fractionation (EE-GMF) method for the quantitative characterization of the excited-state properties of locally excited molecular clusters. The accuracy of the EE-GMF method is systematically assessed for oligoacene crystals. Our result demonstrates that the EE-GMF method is capable of providing the lowest vertical singlet (S1) and triplet excitation energies (T1), in excellent agreement with the full-system quantum mechanical calculations. Using this method, we also investigated the performance of different density functionals in predicting the excited-state properties of the oligoacene crystals. Among the 13 tested functionals, B3LYP and MN15 give the two lowest overall mean unsigned errors with reference to the experimental S1 and T1 excitation energies. The EE-GMF approach can be readily utilized for studying the excited-state properties of large-scale organic solids at diverse ab initio levels.
Collapse
Affiliation(s)
- Jinfeng Liu
- Department of Basic Medicine and Clinical Pharmacy , China Pharmaceutical University , Nanjing 210009 , China
| | | | - William J Glover
- NYU Shanghai , Shanghai 200122 , China.,NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai 200062 , China.,Department of Chemistry , New York University , New York , New York 10003 , United States
| | - Xiao He
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai , Shanghai 200062 , China
| |
Collapse
|
24
|
Dutta AK, Saitow M, Demoulin B, Neese F, Izsák R. A domain-based local pair natural orbital implementation of the equation of motion coupled cluster method for electron attached states. J Chem Phys 2019; 150:164123. [DOI: 10.1063/1.5089637] [Citation(s) in RCA: 45] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
| | - Masaaki Saitow
- Department of Chemistry, Graduate School of Science, Nagoya University, 1-5 Chikusa-ku, Nagoya, Aichi 464-8602, Japan
| | - Baptiste Demoulin
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Róbert Izsák
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
25
|
Mester D, Kállay M. Reduced-Scaling Approach for Configuration Interaction Singles and Time-Dependent Density Functional Theory Calculations Using Hybrid Functionals. J Chem Theory Comput 2019; 15:1690-1704. [DOI: 10.1021/acs.jctc.8b01199] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2023]
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box
91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box
91, H-1521 Budapest, Hungary
| |
Collapse
|
26
|
Park YC, Perera A, Bartlett RJ. Low scaling EOM-CCSD and EOM-MBPT(2) method with natural transition orbitals. J Chem Phys 2018; 149:184103. [DOI: 10.1063/1.5045340] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Young Choon Park
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Rodney J. Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
27
|
Dutta AK, Saitow M, Riplinger C, Neese F, Izsák R. A near-linear scaling equation of motion coupled cluster method for ionized states. J Chem Phys 2018; 148:244101. [DOI: 10.1063/1.5029470] [Citation(s) in RCA: 40] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Achintya Kumar Dutta
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany
| | - Masaaki Saitow
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, D-45470 Mülheim an der Ruhr, Germany
| | | | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Róbert Izsák
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
28
|
Frank MS, Hättig C. A pair natural orbital based implementation of CCSD excitation energies within the framework of linear response theory. J Chem Phys 2018; 148:134102. [PMID: 29626892 DOI: 10.1063/1.5018514] [Citation(s) in RCA: 39] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a pair natural orbital (PNO)-based implementation of coupled cluster singles and doubles (CCSD) excitation energies that builds upon the previously proposed state-specific PNO approach to the excited state eigenvalue problem. We construct the excited state PNOs for each state separately in a truncated orbital specific virtual basis and use a local density-fitting approximation to achieve an at most quadratic scaling of the computational costs for the PNO construction. The earlier reported excited state PNO construction is generalized such that a smooth convergence of the results for charge transfer states is ensured for general coupled cluster methods. We investigate the accuracy of our implementation by applying it to a large and diverse test set comprising 153 singlet excitations in organic molecules. Already moderate PNO thresholds yield mean absolute errors below 0.01 eV. The performance of the implementation is investigated through the calculations on alkene chains and reveals an at most cubic cost-scaling for the CCSD iterations with the system size.
Collapse
Affiliation(s)
- Marius S Frank
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44801 Bochum, Germany
| | - Christof Hättig
- Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44801 Bochum, Germany
| |
Collapse
|
29
|
Mester D, Nagy PR, Kállay M. Reduced-cost second-order algebraic-diagrammatic construction method for excitation energies and transition moments. J Chem Phys 2018. [DOI: 10.1063/1.5021832] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
- Dávid Mester
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Péter R. Nagy
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
30
|
Dutta AK, Nooijen M, Neese F, Izsák R. Exploring the Accuracy of a Low Scaling Similarity Transformed Equation of Motion Method for Vertical Excitation Energies. J Chem Theory Comput 2017; 14:72-91. [DOI: 10.1021/acs.jctc.7b00802] [Citation(s) in RCA: 72] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Achintya Kumar Dutta
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Marcel Nooijen
- Department of Chemistry, University of Waterloo, Waterloo, Ontario N2L 3G1, Canada
| | - Frank Neese
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Róbert Izsák
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
31
|
Bennie SJ, Curchod BFE, Manby FR, Glowacki DR. Pushing the Limits of EOM-CCSD with Projector-Based Embedding for Excitation Energies. J Phys Chem Lett 2017; 8:5559-5565. [PMID: 29076727 DOI: 10.1021/acs.jpclett.7b02500] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
The calculation of accurate excitation energies using ab initio electronic structure methods such as standard equation of motion coupled cluster singles and doubles (EOM-CCSD) has been cost prohibitive for large systems. In this work, we use a simple projector-based embedding scheme to calculate the EOM-CCSD excitation energies of acrolein solvated in water molecules modeled using density functional theory (DFT). We demonstrate the accuracy of this approach gives excitation energies within 0.01 eV of full EOM-CCSD, but with significantly reduced computational cost. This approach is also shown to be relatively invariant to the choice of functional used in the environment and allows for the description of systems with large numbers of basis functions (>1000) to be treated using state-of-the-art wave function methods. The flexibility of embedding to select orbitals to add to the excited-state method provides insights into the origins of the excitations and can reduce artifacts that could arise in traditional linear response time-dependent DFT (LR-TDDFT).
Collapse
Affiliation(s)
- Simon J Bennie
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Bristol, BS8 1TH, United Kingdom
| | - Basile F E Curchod
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Bristol, BS8 1TH, United Kingdom
| | - Frederick R Manby
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Bristol, BS8 1TH, United Kingdom
| | - David R Glowacki
- Centre for Computational Chemistry, School of Chemistry, University of Bristol , Bristol, BS8 1TH, United Kingdom
- Department of Computer Science, University of Bristol , Bristol, BS8 1UB, United Kingdom
| |
Collapse
|
32
|
Baudin P, Kristensen K. Correlated natural transition orbital framework for low-scaling excitation energy calculations (CorNFLEx). J Chem Phys 2017; 146:214114. [PMID: 28595400 PMCID: PMC5462619 DOI: 10.1063/1.4984820] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2017] [Accepted: 05/18/2017] [Indexed: 11/15/2022] Open
Abstract
We present a new framework for calculating coupled cluster (CC) excitation energies at a reduced computational cost. It relies on correlated natural transition orbitals (NTOs), denoted CIS(D')-NTOs, which are obtained by diagonalizing generalized hole and particle density matrices determined from configuration interaction singles (CIS) information and additional terms that represent correlation effects. A transition-specific reduced orbital space is determined based on the eigenvalues of the CIS(D')-NTOs, and a standard CC excitation energy calculation is then performed in that reduced orbital space. The new method is denoted CorNFLEx (Correlated Natural transition orbital Framework for Low-scaling Excitation energy calculations). We calculate second-order approximate CC singles and doubles (CC2) excitation energies for a test set of organic molecules and demonstrate that CorNFLEx yields excitation energies of CC2 quality at a significantly reduced computational cost, even for relatively small systems and delocalized electronic transitions. In order to illustrate the potential of the method for large molecules, we also apply CorNFLEx to calculate CC2 excitation energies for a series of solvated formamide clusters (up to 4836 basis functions).
Collapse
Affiliation(s)
- Pablo Baudin
- Department of Chemistry, qLEAP Center for Theoretical Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Kasper Kristensen
- Department of Chemistry, qLEAP Center for Theoretical Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|
33
|
Dutta AK, Neese F, Izsák R. A simple scheme for calculating approximate transition moments within the equation of motion expectation value formalism. J Chem Phys 2017; 146:214111. [PMID: 28595413 PMCID: PMC5461178 DOI: 10.1063/1.4984618] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2017] [Accepted: 05/17/2017] [Indexed: 11/15/2022] Open
Abstract
A simple scheme for calculating approximate transition moments within the framework of the equation of motion coupled cluster method is proposed. It relies on a matrix inversion technique to calculate the excited state left eigenvectors and requires no additional cost over that of the excitation energy calculation. The new approximation gives almost identical UV-Vis spectra to that obtained using the standard equation of motion coupled cluster method with single and double excitations for molecules in a standard test set.
Collapse
Affiliation(s)
- Achintya Kumar Dutta
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Frank Neese
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany
| | - Róbert Izsák
- Max-Planck-Institut für Chemische Energiekonversion, Stiftstr. 34-36, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
34
|
Mester D, Nagy PR, Kállay M. Reduced-cost linear-response CC2 method based on natural orbitals and natural auxiliary functions. J Chem Phys 2017; 146:194102. [PMID: 28527453 DOI: 10.1063/1.4983277] [Citation(s) in RCA: 43] [Impact Index Per Article: 6.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A reduced-cost density fitting (DF) linear-response second-order coupled-cluster (CC2) method has been developed for the evaluation of excitation energies. The method is based on the simultaneous truncation of the molecular orbital (MO) basis and the auxiliary basis set used for the DF approximation. For the reduction of the size of the MO basis, state-specific natural orbitals (NOs) are constructed for each excited state using the average of the second-order Møller-Plesset (MP2) and the corresponding configuration interaction singles with perturbative doubles [CIS(D)] density matrices. After removing the NOs of low occupation number, natural auxiliary functions (NAFs) are constructed [M. Kállay, J. Chem. Phys. 141, 244113 (2014)], and the NAF basis is also truncated. Our results show that, for a triple-zeta basis set, about 60% of the virtual MOs can be dropped, while the size of the fitting basis can be reduced by a factor of five. This results in a dramatic reduction of the computational costs of the solution of the CC2 equations, which are in our approach about as expensive as the evaluation of the MP2 and CIS(D) density matrices. All in all, an average speedup of more than an order of magnitude can be achieved at the expense of a mean absolute error of 0.02 eV in the calculated excitation energies compared to the canonical CC2 results. Our benchmark calculations demonstrate that the new approach enables the efficient computation of CC2 excitation energies for excited states of all types of medium-sized molecules composed of up to 100 atoms with triple-zeta quality basis sets.
Collapse
Affiliation(s)
- Dávid Mester
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Péter R Nagy
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- MTA-BME Lendület Quantum Chemistry Research Group, Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
35
|
Baudin P, Kjærgaard T, Kristensen K. CC2 oscillator strengths within the local framework for calculating excitation energies (LoFEx). J Chem Phys 2017; 146:144107. [PMID: 28411600 DOI: 10.1063/1.4979713] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
In a recent work [P. Baudin and K. Kristensen, J. Chem. Phys. 144, 224106 (2016)], we introduced a local framework for calculating excitation energies (LoFEx), based on second-order approximated coupled cluster (CC2) linear-response theory. LoFEx is a black-box method in which a reduced excitation orbital space (XOS) is optimized to provide coupled cluster (CC) excitation energies at a reduced computational cost. In this article, we present an extension of the LoFEx algorithm to the calculation of CC2 oscillator strengths. Two different strategies are suggested, in which the size of the XOS is determined based on the excitation energy or the oscillator strength of the targeted transitions. The two strategies are applied to a set of medium-sized organic molecules in order to assess both the accuracy and the computational cost of the methods. The results show that CC2 excitation energies and oscillator strengths can be calculated at a reduced computational cost, provided that the targeted transitions are local compared to the size of the molecule. To illustrate the potential of LoFEx for large molecules, both strategies have been successfully applied to the lowest transition of the bivalirudin molecule (4255 basis functions) and compared with time-dependent density functional theory.
Collapse
Affiliation(s)
- Pablo Baudin
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Thomas Kjærgaard
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Kasper Kristensen
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| |
Collapse
|
36
|
Høyvik IM, Myhre RH, Koch H. Correlated natural transition orbitals for core excitation energies in multilevel coupled cluster models. J Chem Phys 2017; 146:144109. [DOI: 10.1063/1.4979908] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
|
37
|
Baudin P, Bykov D, Liakh D, Ettenhuber P, Kristensen K. A local framework for calculating coupled cluster singles and doubles excitation energies (LoFEx-CCSD). Mol Phys 2017. [DOI: 10.1080/00268976.2017.1290836] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
Affiliation(s)
- Pablo Baudin
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Aarhus C, Denmark
| | - Dmytro Bykov
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Aarhus C, Denmark
| | - Dmitry Liakh
- Scientific Computing Group, National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN, USA
| | - Patrick Ettenhuber
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Aarhus C, Denmark
| | - Kasper Kristensen
- qLEAP Center for Theoretical Chemistry, Department of Chemistry, Aarhus University, Aarhus C, Denmark
| |
Collapse
|
38
|
Dutta AK, Nooijen M, Neese F, Izsák R. Automatic active space selection for the similarity transformed equations of motion coupled cluster method. J Chem Phys 2017; 146:074103. [DOI: 10.1063/1.4976130] [Citation(s) in RCA: 45] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|