1
|
Fedorov DG. Use of caps in the auxiliary basis set formulation of the fragment molecular orbital method. J Comput Chem 2024; 45:1540-1551. [PMID: 38490813 DOI: 10.1002/jcc.27345] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2024] [Revised: 02/19/2024] [Accepted: 02/29/2024] [Indexed: 03/17/2024]
Abstract
An auxiliary polarization formulation of the fragment molecular orbital (FMO) method is developed, combining a basis set correction computed for capped isolated fragments with a polarization obtained from uncapped fragments. For a set of organic and inorganic test systems, it is shown that the total energy and atomic charges are accurately reproduced with respect to full unfragmented calculations. It is demonstrated that the method is accurate for computing electronic excited states. The developed approach is applied to rank the isomers of chignolin from experimental NMR data (PDB: 1UAO) according to their relative energy. Contributions of polarization and basis set effects to pair interactions between fragments are elucidated.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Japan
| |
Collapse
|
2
|
Tuckman H, Neuscamman E. Aufbau Suppressed Coupled Cluster Theory for Electronically Excited States. J Chem Theory Comput 2024; 20:2761-2773. [PMID: 38502102 DOI: 10.1021/acs.jctc.3c01285] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/20/2024]
Abstract
We introduce an approach to improve single-reference coupled cluster theory in settings where the Aufbau determinant is absent from or plays only a small role in the true wave function. Using a de-excitation operator that can be efficiently hidden within a similarity transform, we create a coupled cluster wave function in which de-excitations work to suppress the Aufbau determinant and produce wave functions dominated by other determinants. Thanks to an invertible and fully exponential form, the approach is systematically improvable, size consistent, size extensive, and, interestingly, size intensive in a granular way that should make the adoption of some ground state techniques, such as local correlation, relatively straightforward. In this initial study, we apply the general formalism to create a state-specific method for orbital-relaxed, singly excited states. We find that this approach matches the accuracy of similar-cost equation-of-motion methods in valence excitations while offering improved accuracy for charge transfer states. We also find the approach to be more accurate than excited-state-specific perturbation theory in both types of states.
Collapse
Affiliation(s)
- Harrison Tuckman
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
3
|
Watrous AG, Davis MC, Fortenberry RC. Performance of EOM-CCSD(T)(a)*-Based Quartic Force Fields in Computing Fundamental, Anharmonic Vibrational Frequencies of Molecular Electronically Excited States with Application to the Ã1A″ State of :CCH 2 (Vinylidene). J Phys Chem A 2024; 128:2150-2161. [PMID: 38466814 DOI: 10.1021/acs.jpca.3c08168] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/13/2024]
Abstract
Highly accurate anharmonic vibrational frequencies of electronically excited states are not as easily computed as their ground electronic state counterparts, but recently developed approximate triple excited state methods may be changing that. One emerging excited state method is equation of motion coupled cluster theory at the singles and doubles level with perturbative triples computed via the (a)* formalism, EOMEE-CCSD(T)(a)*. One of the most employed means for the ready computation of vibrational anharmonic frequencies for ground electronic states is second-order vibrational perturbation theory (VPT2), a theory based on quartic force fields (QFFs),fourth-order Taylor series expansions of the potential portion of the internuclear Watson Hamiltonian. The combination of these two is herein benchmarked for its performance for use as a means of computing rovibrational spectra of electronically excited states. Specifically, the EOMEE-CCSD(T)(a)* approach employing a complete basis set extrapolation along with core electron inclusion and relativity (the so-called "CcCR" approach) defining the QFF produces anharmonic fundamental vibrational frequencies within 2.83%, on the average, of reported gas-phase experimentally assigned values for the test set including the A ~ 1 A ″ states of HCF, HCCl, HSiF, HNO, and HPO. However, some states have exceptional accuracy in the fundamentals, most notably for ν2 of A ~ 1 A ″ HCCl in which the CcCR QFF value is within 1.8 cm-1 at 927.9 cm-1 (or 0.2%) of the experiment. Additionally, this approach produces rotational constants to, on the absolute average, within 0.41% of available experimental data, showcasing notable accuracy in the computation of rovibronic spectral data. Furthermore, utilizing a hybrid approach composed of harmonic CcCR force constants along with a set of simple EOMEE-CCSD(T)(a)*/aug-cc-pVQZ QFF cubic and quartic force constants is faster than using pure CcCR and better represents those modes that suffer from numerical instability in the anharmonic portion of the QFF, implying that this so-called "CcCR + QZ" QFF approach may be the best for future applications. Finally, complete, rovibrational spectral data are provided for A ~ 1 A 2 :CCH2, a molecule of potential astrochemical interest, in order to aid in its potential future experimental rovibronic characterization.
Collapse
Affiliation(s)
- Alexandria G Watrous
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| | - Megan C Davis
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| | - Ryan C Fortenberry
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| |
Collapse
|
4
|
Fedorov DG. Analysis of Site Energies and Excitonic Couplings: The Role of Symmetry and Polarization. J Phys Chem A 2024; 128:1154-1162. [PMID: 38302431 DOI: 10.1021/acs.jpca.3c06293] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2024]
Abstract
An excitonic coupling model is developed based on an equation-of-motion coupled cluster combined with the fragment molecular orbital method. The effects of polarization and excitonic coupling on the splitting of quasi-degenerate levels in systems containing multiple chromophores are elucidated on dimers of formaldehyde, water, formic acid, hydrogen fluoride, and carbon monoxide. It is shown that the level structure is mainly determined by the mutual polarization of chromophores and to a lesser extent by the excitonic coupling. The role of symmetry in excitonic coupling in dimers is discussed. The excitonic coupling between all residues in the photoactive yellow protein (PDB: 2PHY) is analyzed.
Collapse
Affiliation(s)
- Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Central 2, Umezono 1-1-1, Tsukuba 305-8568, Japan
| |
Collapse
|
5
|
Tuckman H, Neuscamman E. Excited-State-Specific Pseudoprojected Coupled-Cluster Theory. J Chem Theory Comput 2023; 19:6160-6171. [PMID: 37676752 DOI: 10.1021/acs.jctc.3c00194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/09/2023]
Abstract
We present an excited-state-specific coupled-cluster approach in which both the molecular orbitals and cluster amplitudes are optimized for an individual excited state. The theory is formulated via a pseudoprojection of the traditional coupled-cluster wavefunction that allows correlation effects to be introduced atop an excited-state mean field starting point. The approach shares much in common with ground-state CCSD, including size extensivity and an N6 cost scaling. Preliminary numerical tests show that, when augmented with N5 cost perturbative corrections for key terms, the method can improve over excited-state-specific second-order perturbation theory in valence, charge transfer, and Rydberg states.
Collapse
Affiliation(s)
- Harrison Tuckman
- Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
6
|
Wang X, Hasan M, Fan L, Wang Y, Li H, Slaughter DS, Centurion M. Mass-selected ion-molecule cluster beam apparatus for ultrafast photofragmentation studies. THE REVIEW OF SCIENTIFIC INSTRUMENTS 2023; 94:095111. [PMID: 37724931 DOI: 10.1063/5.0148194] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Accepted: 08/31/2023] [Indexed: 09/21/2023]
Abstract
We describe an apparatus for investigating the excited-state dissociation dynamics of mass-selected ion-molecule clusters by mass-resolving and detecting photofragment-ions and neutrals, in coincidence, using an ultrafast laser operating at high repetition rates. The apparatus comprises a source that generates ion-molecule clusters, a time-of-flight spectrometer, and a mass filter that selects the desired anions, and a linear-plus-quadratic reflectron mass spectrometer that discriminates the fragment anions after the femtosecond laser excites the clusters. The fragment neutrals and anions are then captured by two channeltron detectors. The apparatus performance is tested by measuring the photofragments: I-, CF3I-, and neutrals from photoexcitation of the ion-molecule cluster CF3I·I- using femtosecond UV laser pulses with a wavelength of 266 nm. The experimental results are compared with our ground state and excited state electronic structure calculations as well as the existing results and calculations, with particular attention to the generation mechanism of the anion fragments and dissociation channels of the ion-molecule cluster CF3I·I- in the charge-transfer excited state.
Collapse
Affiliation(s)
- Xiaojun Wang
- Department of Physics and Astronomy, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA
| | - Mahmudul Hasan
- Department of Physics and Astronomy, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA
| | - Lin Fan
- Department of Physics and Astronomy, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA
| | - Yibo Wang
- Department of Physics and Astronomy, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA
| | - Hui Li
- Department of Chemistry, Nebraska Center for Materials and Nanoscience, and Center for Integrated Biomolecular Communication, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA
| | - Daniel S Slaughter
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Rd., Berkeley, California 94720, USA
| | - Martin Centurion
- Department of Physics and Astronomy, University of Nebraska-Lincoln, Lincoln, Nebraska 68588, USA
| |
Collapse
|
7
|
Gururangan K, Piecuch P. Converging high-level coupled-cluster energetics via adaptive selection of excitation manifolds driven by moment expansions. J Chem Phys 2023; 159:084108. [PMID: 37610021 DOI: 10.1063/5.0162873] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2023] [Accepted: 08/03/2023] [Indexed: 08/24/2023] Open
Abstract
A novel approach to rapidly converging high-level coupled-cluster (CC) energetics in an automated fashion is proposed. The key idea is an adaptive selection of excitation manifolds defining higher--than--two-body components of the cluster operator inspired by CC(P;Q) moment expansions. The usefulness of the resulting methodology is illustrated by molecular examples where the goal is to recover the electronic energies obtained using the CC method with a full treatment of singly, doubly, and triply excited clusters (CCSDT) when the noniterative triples corrections to CCSD fail.
Collapse
Affiliation(s)
- Karthik Gururangan
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
- Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|
8
|
Magoulas I, Evangelista FA. Unitary Coupled Cluster: Seizing the Quantum Moment. J Phys Chem A 2023; 127:6567-6576. [PMID: 37523485 PMCID: PMC10424243 DOI: 10.1021/acs.jpca.3c02781] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2023] [Revised: 07/08/2023] [Indexed: 08/02/2023]
Abstract
Shallow, CNOT-efficient quantum circuits are crucial for performing accurate computational chemistry simulations on current noisy quantum hardware. Here, we explore the usefulness of noniterative energy corrections, based on the method of moments of coupled-cluster theory, for accelerating convergence toward full configuration interaction. Our preliminary numerical results relying on iteratively constructed ansätze suggest that chemically accurate energies can be obtained with substantially more compact circuits, implying enhanced resilience to gate and decoherence noise.
Collapse
Affiliation(s)
- Ilias Magoulas
- Department of Chemistry and
Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Francesco A. Evangelista
- Department of Chemistry and
Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
9
|
Shee A, Yeh CN, Peng B, Kowalski K, Zgid D. Triple Excitations in Green's Function Coupled Cluster Solver for Studies of Strongly Correlated Systems in the Framework of Self-Energy Embedding Theory. J Phys Chem Lett 2023; 14:2416-2424. [PMID: 36856741 DOI: 10.1021/acs.jpclett.2c03616] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/18/2023]
Abstract
Embedding theories became important approaches used for accurate calculations of both molecules and solids. In these theories, a small chosen subset of orbitals is treated with an accurate method, called an impurity solver, capable of describing higher correlation effects. Ideally, such a chosen fragment should contain multiple orbitals responsible for the chemical and physical behavior of the compound. Handling a large number of chosen orbitals presents a very significant challenge for the current generation of solvers used in the physics and chemistry community. Here, we develop a Green's function coupled cluster singles doubles and triples (GFCCSDT) solver that can be used for a quantitative description in both molecules and solids. This solver allows us to treat orbital spaces that are inaccessible to other accurate solvers. At the same time, GFCCSDT maintains high accuracy of the resulting self-energy. Moreover, in conjunction with the GFCCSD solver, it allows us to test the systematic convergence of computational studies. Developing the CC family of solvers paves the road to fully systematic Green's function embedding calculations in solids. In this paper, we focus on the investigation of GFCCSDT self-energies for a strongly correlated problem of SrMnO3 solid. Subsequently, we apply this solver to solid MnO showing that an approximate variant of GFCCSDT is capable of yielding a high accuracy orbital resolved spectral function.
Collapse
Affiliation(s)
- Avijit Shee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Chia-Nan Yeh
- Department of Physics, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Bo Peng
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Karol Kowalski
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Dominika Zgid
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
10
|
Nowak A, Boguslawski K. A configuration interaction correction on top of pair coupled cluster doubles. Phys Chem Chem Phys 2023; 25:7289-7301. [PMID: 36810525 DOI: 10.1039/d2cp05171k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Abstract
Numerous numerical studies have shown that geminal-based methods are a promising direction to model strongly correlated systems with low computational costs. Several strategies have been introduced to capture the missing dynamical correlation effects, which typically exploit a posteriori corrections to account for correlation effects associated with broken-pair states or inter-geminal correlations. In this article, we scrutinize the accuracy of the pair coupled cluster doubles (pCCD) method extended by configuration interaction (CI) theory. Specifically, we benchmark various CI models, including, at most double excitations against selected CC corrections as well as conventional single-reference CC methods. A simple Davidson correction is also tested. The accuracy of the proposed pCCD-CI approaches is assessed for challenging small model systems such as the N2 and F2 dimers and various di- and triatomic actinide-containing compounds. In general, the proposed CI methods considerably improve spectroscopic constants compared to the conventional CCSD approach, provided a Davidson correction is included in the theoretical model. At the same time, their accuracy lies between those of the linearized frozen pCCD and frozen pCCD variants.
Collapse
Affiliation(s)
- Artur Nowak
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland.
| | - Katharina Boguslawski
- Institute of Physics, Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Toruń, Grudziadzka 5, 87-100 Toruń, Poland.
| |
Collapse
|
11
|
Ning J, Truhlar DG. Electronic Excitation of ortho-Fluorothiophenol. J Phys Chem A 2023; 127:1469-1474. [PMID: 36749710 DOI: 10.1021/acs.jpca.2c08674] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
Abstract
ortho-Fluorothiophenol (o-FTP) photodissociates through the well-known πσ* process. The fluorine atom of o-FTP introduces a feature in the photodissociation of o-FTP that does not occur in most other πσ* processes because the fluorine atom can form a hydrogen bond with the hydrogen atom of the SH group. Theoretical computations can serve as a good way to study these reactions because they usually proceed very quickly, and the current spectroscopies cannot probe the details of the processes as thoroughly as theory can. Here we use completely renormalized equation-of-motion coupled cluster theory with single and double excitations and a quasiperturbative treatment of connected triple excitations (CR-EOM-CCSD(T)) and quasidegenerate perturbation theory, in particular extended multistate complete-active-space second-order perturbation theory (XMS- CASPT2), to calculate the four lowest singlet states of o-FTP and hybrid density functional theory to optimize the geometries of the two lowest singlet states. We find that ten active electrons in nine active orbitals are sufficient to provide a good reference function for all four states. We find that the ground electronic state and the first excited singlet state both exhibit strongly bent hydrogen bonds. We also use density functional theory with the Tamm-Dancoff approximation and the SMD solvation model to successfully simulate the electronic spectrum of o-FTP in n-hexane solvent.
Collapse
Affiliation(s)
- Jiaxin Ning
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| | - Donald G Truhlar
- Department of Chemistry, Chemical Theory Center, and Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, United States
| |
Collapse
|
12
|
Köse ME. Estimation of Excited-State Geometries of Benzene and Fluorobenzene through Vibronic Analyses of Absorption Spectra. ACS OMEGA 2022; 7:32764-32774. [PMID: 36120020 PMCID: PMC9476181 DOI: 10.1021/acsomega.2c04615] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/21/2022] [Accepted: 08/26/2022] [Indexed: 06/15/2023]
Abstract
The parameters used in theoretical modeling of vibrational patterns within Franck-Condon (FC) approximation can be adjusted to match the vibrationally well-resolved experimental absorption spectrum of molecules. These simulation parameters can then be used to reveal the structural changes occurring between the initial and final states assuming the harmonic oscillator approximation holds for both states. Such a theoretical approach has been applied to benzene and fluorobenzene to disclose the first excited-state geometries of both compounds. The carbon-carbon bond length of benzene in the 1B2u state has been calculated as 1.430 Å, which is in very good agreement with the experimental bond length of 1.432 Å. The FC spectral fit method has been exploited to reveal the 1B2 state of fluorobenzene as well. Commonly employed density functional theory (DFT) and time-dependent DFT methods have been used to calculate the ground- and excited-state geometries of both compounds, respectively. The comparison of geometrical parameters and vibrational frequencies at the relevant states shows that frequently used hybrid functionals perform quite well in the ground state, whereas their performances drop considerably while predicting the excited-state properties. Among the hybrid functionals studied, TD-B3LYP with 6-31+G(d) basis set can be chosen to calculate the excited-state properties of molecules, albeit with much less anticipation of accuracy from the performance that B3LYP usually shows at the ground state.
Collapse
|
13
|
Fu Y, Niu Y, Liang X, Liu H, Zou W. High-Level Theoretical Study of the C 1Π States of GaCl, GaBr, InCl, and InBr. J Phys Chem A 2022; 126:5565-5573. [PMID: 35948043 DOI: 10.1021/acs.jpca.2c04240] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
For the first singlet excited states C1Π of IIIA group monohalides GaCl, GaBr, InCl, and InBr, a very shallow potential well may exist in the flat potential energy curve, which poses a challenge to the theoretical accuracy. In this study, high-level theoretical calculations have been performed through the Feller-Peterson-Dixon composite approach to study the C1Π states, and the obtained spectroscopic constants agree well with the known experimental ones. It is found that the C1Π states are sensitive to the size of basis functions, spin-orbit coupling, and strong correlations mainly due to triple excitations. The final results show that the C1Π states of GaCl and InCl are quasi-bound with one (v' = 0) and four (v' = 0-3) vibrational levels, respectively, being consistent with the experimental findings, whereas the C1Π states of GaBr and InBr are repulsive. Our conclusions deny the existence of higher vibrational levels v' = 1 for GaCl, v' > 3 for InCl, and v' ≥ 0 for InBr in previous experimental and theoretical studies of C1Π.
Collapse
Affiliation(s)
- Yingxiang Fu
- Institute of Modern Physics, Northwest University, and Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi'an 710127, Shaanxi, P. R. China
| | - Yue Niu
- Institute of Modern Physics, Northwest University, and Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi'an 710127, Shaanxi, P. R. China
| | - Ximin Liang
- Institute of Modern Physics, Northwest University, and Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi'an 710127, Shaanxi, P. R. China
| | - Haitao Liu
- Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing, 100088, P. R. China
| | - Wenli Zou
- Institute of Modern Physics, Northwest University, and Shaanxi Key Laboratory for Theoretical Physics Frontiers, Xi'an 710127, Shaanxi, P. R. China
| |
Collapse
|
14
|
Paul AC, Folkestad SD, Myhre RH, Koch H. Oscillator Strengths in the Framework of Equation of Motion Multilevel CC3. J Chem Theory Comput 2022; 18:5246-5258. [PMID: 35921447 PMCID: PMC9476665 DOI: 10.1021/acs.jctc.2c00164] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
We present an efficient implementation of the equation
of motion
oscillator strengths for the closed-shell multilevel coupled cluster
singles and doubles with perturbative triples method (MLCC3) in the
electronic structure program eT. The orbital space is split into an active part treated with
CC3 and an inactive part computed at the coupled cluster singles and
doubles (CCSD) level of theory. Asymptotically, the CC3 contribution
scales as floating-point operations, where nV is the total number of virtual orbitals while nv and no are the
number of active virtual and occupied orbitals, respectively. The
CC3 contribution, thus, only scales linearly with the full system
size and can become negligible compared to the cost of CCSD. We demonstrate
the capabilities of our implementation by calculating the ultraviolet–visible
spectrum of azobenzene and a core excited state of betaine 30 with
more than 1000 molecular orbitals.
Collapse
Affiliation(s)
- Alexander C Paul
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway
| | - Sarai Dery Folkestad
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway
| | - Rolf H Myhre
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway.,Scuola Normale Superiore, Piazza dei Cavaleri 7, 56126 Pisa, Italy
| |
Collapse
|
15
|
Liu J, Matthews DA, Cheng L. Quadratic Unitary Coupled-Cluster Singles and Doubles Scheme: Efficient Implementation, Benchmark Study, and Formulation of an Extended Version. J Chem Theory Comput 2022; 18:2281-2291. [PMID: 35312299 DOI: 10.1021/acs.jctc.1c01210] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
An efficient implementation of the quadratic unitary coupled-cluster singles and doubles (qUCCSD) scheme for calculations of electronic ground and excited states using an unrestricted molecular spin-orbital formulation and an efficient tensor contraction library is reported. The accuracy of the qUCCSD scheme and the efficiency of the present implementation are demonstrated using extensive benchmark calculations of excitation energies and an application to S0 → S1 vertical excitation energies for cis- and trans-4a,4b-dihydrotriphenylene. The qUCCSD scheme has been shown to provide improved excitation energies compared with the UCC3 scheme formulated based on perturbation theory. A UCC truncation scheme that can provide excitation energies correct through the fourth order is also presented to further improve the accuracy of the qUCCSD scheme.
Collapse
Affiliation(s)
- Junzi Liu
- School of Chemistry and Biological Engineering, University of Science and Technology Beijing, Beijing, 100083, China
| | - Devin A Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
16
|
Sold S, Mummaneni BC, Michenfelder NC, Peng Y, Powell AK, Unterreiner AN, Lefkidis G, Hübner W. Experimental and Theoretical Study of the Ultrafast Dynamics of a Ni 2 Dy 2 -Compound in DMF After UV/Vis Photoexcitation. ChemistryOpen 2021; 11:e202100153. [PMID: 34931474 PMCID: PMC9059312 DOI: 10.1002/open.202100153] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2021] [Revised: 11/04/2021] [Indexed: 01/30/2023] Open
Abstract
We present a combined experimental and theoretical study of the ultrafast transient absorption spectroscopy results of a {Ni2Dy2}‐compound in DMF, which can be considered as a prototypic molecule for single molecule magnets. We apply state‐of‐the‐art ab initio quantum chemistry to quantitatively describe the optical properties of an inorganic complex system comprising ten atoms to form the chromophoric unit, which is further stabilized by surrounding ligands. Two different basis sets are used for the calculations to specifically identify two dominant peaks in the ground state. Furthermore, we theoretically propagate the compound's correlated many‐body wavefunction under the influence of a laser pulse as well as relaxation processes and compare against the time‐resolved absorption spectra. The experimental data can be described with a time constant of several hundreds of femtoseconds attributed to vibrational relaxation and trapping into states localized within the band gap. A second time constant is ascribed to the excited state while trap states show lifetimes on a longer timescale. The theoretical propagation is performed with the density‐matrix formalism and the Lindblad superoperator, which couples the system to a thermal bath, allowing us to extract relaxation times from first principles.
Collapse
Affiliation(s)
- S Sold
- Department of Physics and Research Center OPTIMAS, Technische Universität Kaiserslautern, P.O. Box 3049, 67653, Kaiserslautern, Germany
| | - B C Mummaneni
- Department of Physics and Research Center OPTIMAS, Technische Universität Kaiserslautern, P.O. Box 3049, 67653, Kaiserslautern, Germany
| | - N C Michenfelder
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Fritz-Haber-Weg 2, 76131, Karlsruhe, Germany
| | - Y Peng
- Institute of Inorganic Chemistry, Karlsruhe Institute of Technology, Engesserstrasse 15, 76131, Karlsruhe, Germany.,Institute of Nanotechnology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344, Eggenstein-Leopoldshafen, Germany
| | - A K Powell
- Institute of Inorganic Chemistry, Karlsruhe Institute of Technology, Engesserstrasse 15, 76131, Karlsruhe, Germany.,Institute of Nanotechnology, Karlsruhe Institute of Technology, Hermann-von-Helmholtz-Platz 1, 76344, Eggenstein-Leopoldshafen, Germany
| | - A-N Unterreiner
- Institute of Physical Chemistry, Karlsruhe Institute of Technology, Fritz-Haber-Weg 2, 76131, Karlsruhe, Germany
| | - G Lefkidis
- Department of Physics and Research Center OPTIMAS, Technische Universität Kaiserslautern, P.O. Box 3049, 67653, Kaiserslautern, Germany.,School of Mechanics, Civil Engineering and Architecture, Northwestern Polytechnical University, Xi'an, 710072, China
| | - W Hübner
- Department of Physics and Research Center OPTIMAS, Technische Universität Kaiserslautern, P.O. Box 3049, 67653, Kaiserslautern, Germany
| |
Collapse
|
17
|
Capistran BA, Yuwono SH, Moemeni M, Maity S, Vahdani A, Borhan B, Jackson JE, Piecuch P, Dantus M, Blanchard GJ. Excited-State Dynamics of a Substituted Fluorene Derivative. The Central Role of Hydrogen Bonding Interactions with the Solvent. J Phys Chem B 2021; 125:12242-12253. [PMID: 34726920 PMCID: PMC9254887 DOI: 10.1021/acs.jpcb.1c06474] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Substituted fluorene structures have demonstrated unusual photochemical properties. Previous reports on the substituted fluorene Schiff base FR0-SB demonstrated super photobase behavior with a ΔpKb of ∼14 upon photoexcitation. In an effort to understand the basis for this unusual behavior, we have examined the electronic structure and relaxation dynamics of the structural precursor of FR0-SB, the aldehyde FR0, in protic and aprotic solvents using time-resolved fluorescence spectroscopy and quantum chemical calculations. The calculations show three excited singlet states in relatively close energetic proximity. The spectroscopic data are consistent with relaxation dynamics from these electronic states that depend on the presence and concentration of solvent hydroxyl functionality. These results underscore the central role of solvent hydrogen bonding to the FR0 aldehyde oxygen in mediating the relaxation dynamics within this molecule.
Collapse
Affiliation(s)
- Briana A Capistran
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Stephen H Yuwono
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Mehdi Moemeni
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Soham Maity
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Aria Vahdani
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Babak Borhan
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - James E Jackson
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, United States
| | - Marcos Dantus
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, United States
| | - G J Blanchard
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
18
|
Liu J, Cheng L. Unitary coupled-cluster based self-consistent polarization propagator theory: A quadratic unitary coupled-cluster singles and doubles scheme. J Chem Phys 2021; 155:174102. [PMID: 34742195 DOI: 10.1063/5.0062090] [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/14/2022] Open
Abstract
The development of a quadratic unitary coupled-cluster singles and doubles (qUCCSD) based self-consistent polarization propagator method is reported. We present a simple strategy for truncating the commutator expansion of the unitary version of coupled-cluster transformed Hamiltonian H̄. The qUCCSD method for the electronic ground state includes up to double commutators for the amplitude equations and up to cubic commutators for the energy expression. The qUCCSD excited-state eigenvalue equations include up to double commutators for the singles-singles block of H̄, single commutators for the singles-doubles and doubles-singles blocks, and the bare Hamiltonian for the doubles-doubles block. Benchmark qUCCSD calculations of the ground-state properties and excitation energies for representative molecules demonstrate significant improvement of the accuracy and robustness over the previous UCC3 scheme derived using Møller-Plesset perturbation theory.
Collapse
Affiliation(s)
- Junzi Liu
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
19
|
Park W, Shen J, Lee S, Piecuch P, Filatov M, Choi CH. Internal Conversion between Bright (1 1Bu+) and Dark (2 1Ag-) States in s- trans-Butadiene and s- trans-Hexatriene. J Phys Chem Lett 2021; 12:9720-9729. [PMID: 34590847 DOI: 10.1021/acs.jpclett.1c02707] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Internal conversion (IC) between the two lowest singlet excited states, 11Bu+ and 21Ag-, of s-trans-butadiene and s-trans-hexatriene is investigated using a series of single- and multi- reference wave function and density functional theory (DFT) methodologies. Three independent types of the equation-of-motion coupled-cluster (EOMCC) theory capable of providing an accurate and balanced description of one- as well as two-electron transitions, abbreviated as δ-CR-EOMCC(2,3), DIP-EOMCC(4h2p){No}, and DEA-EOMCC(4p2h){Nu} or DEA-EOMCC(3p1h,4p2h){Nu}, consistently predict that the 11Bu+/21Ag- crossing in both molecules occurs along the bond length alternation coordinate. However, the analogous 11Bu+ and 21Ag- potentials obtained with some multireference approaches, such as CASSCF and MRCIS(D), as well as with the linear-response formulation of time-dependent DFT (TDDFT), do not cross. Hence, caution needs to be exercised when studying the low-lying singlet excited states of polyenes with conventional multiconfigurational methods and TDDFT. The multistate many-body perturbation theory methods, such as XMCQDPT2, do correctly reproduce the curve crossing. Among the simplest and least expensive computational methodologies, the DFT approaches that incorporate the contributions of doubly excited configurations, abbreviated as MRSF (mixed reference spin-flip) TDDFT and SSR(4,4), accurately reproduce our best EOMCC results. This is highly promising for nonadiabatic molecular dynamics simulations in larger systems.
Collapse
Affiliation(s)
- Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, United States
| | - Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
20
|
Casanova-Páez M, Goerigk L. Time-Dependent Long-Range-Corrected Double-Hybrid Density Functionals with Spin-Component and Spin-Opposite Scaling: A Comprehensive Analysis of Singlet-Singlet and Singlet-Triplet Excitation Energies. J Chem Theory Comput 2021; 17:5165-5186. [PMID: 34291643 DOI: 10.1021/acs.jctc.1c00535] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
Following the work on spin-component and spin-opposite scaled (SCS/SOS) global double hybrids for singlet-singlet excitations by Schwabe and Goerigk [ J. Chem. Theory Comput. 2017, 13, 4307-4323] and our own works on new long-range corrected (LC) double hybrids for singlet-singlet and singlet-triplet excitations [ J. Chem. Theory Comput. 2019, 15, 4735-4744 and J. Chem. Phys. 2020, 153, 064106], we present new LC double hybrids with SCS/SOS that demonstrate further improvement over previously published results and methods. We introduce new unscaled and scaled versions of different global and LC double hybrids based on Becke88 or PBE exchange combined with LYP, PBE, or P86 correlation. For singlet-singlet excitations, we cross-validate them on six benchmark sets that cover small to medium-sized chromophores with different excitation types (local-valence, Rydberg, and charge transfer). For singlet-triplet excitations, we perform the cross-validation on three different benchmark sets following the same analysis as in our previous work in 2020. In total, 203 excitations are analyzed. Our results confirm and extend those of Schwabe and Goerigk regarding the superior performance of SCS and SOS variants compared to their unscaled parents by decreasing mean absolute deviations, root-mean-square deviations, or error spans by more than half and bringing absolute mean deviations closer to zero. Our SCS/SOS variants are shown to be highly efficient and robust for the computation of vertical excitation energies, which even outperform specialized double hybrids that also contain an LC in their perturbative part. In particular, our new SCS/SOS-ωPBEPP86 and SCS/SOS-ωB88PP86 functionals are four of the most accurate and robust methods tested in this work, and we fully recommend them for future applications. However, if the relevant SCS and SOS algorithms are not available to the user, we suggest ωPBEPP86 as the best unscaled method in this work.
Collapse
Affiliation(s)
- Marcos Casanova-Páez
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| | - Lars Goerigk
- School of Chemistry, The University of Melbourne, Parkville, Victoria 3010, Australia
| |
Collapse
|
21
|
Mester D, Kállay M. Spin-Scaled Range-Separated Double-Hybrid Density Functional Theory for Excited States. J Chem Theory Comput 2021; 17:4211-4224. [PMID: 34152771 PMCID: PMC8280718 DOI: 10.1021/acs.jctc.1c00422] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2021] [Indexed: 11/28/2022]
Abstract
Our recently presented range-separated (RS) double-hybrid (DH) time-dependent density functional approach [J. Chem. Theory Comput. 17, 927 (2021)] is combined with spin-scaling techniques. The proposed spin-component-scaled (SCS) and scaled-opposite-spin (SOS) variants are thoroughly tested for almost 500 excitations including the most challenging types. This comprehensive study provides useful information not only about the new approaches but also about the most prominent methods in the DH class. The benchmark calculations confirm the robustness of the RS-DH ansatz, while several tendencies and deficiencies are pointed out for the existing functionals. Our results show that the SCS variant consistently improves the results, while the SOS variant preserves the benefits of the original RS-DH method reducing its computational expenses. It is also demonstrated that, besides our approaches, only the nonempirical functionals provide balanced performance for general applications, while particular methods are only suggested for certain types of excitations.
Collapse
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
|
22
|
Davis MC, Fortenberry RC. (T)+EOM Quartic Force Fields for Theoretical Vibrational Spectroscopy of Electronically Excited States. J Chem Theory Comput 2021; 17:4374-4382. [PMID: 34165980 DOI: 10.1021/acs.jctc.1c00307] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
(T)+EOM quartic force fields (QFFs) are proposed for ab initio rovibrational properties of electronically excited states of small molecules. The (T)+EOM method is a simple treatment of the potential surface of the excited state using a composite energy from the CCSD(T) energy for the ground-state configuration and the EOM-CCSD excitation energy for the target state. The method is benchmarked with two open-shell species, HOO and HNF, and two closed-shell species, HNO and HCF. A (T)+EOM QFF with a complete basis set extrapolation (C) and corrections for core correlation (cC) and scalar relativity (R), dubbed (T)+EOM/CcCR, achieves a mean absolute error (MAE) as low as 1.6 cm-1 for the à 2A' state of HOO versus an established benchmark QFF with CCSD(T)-F12b/cc-pVTZ-F12 (F12-TZ) for this variationally accessible electronically excited state. The MAE for anharmonic frequencies for (T)+EOM/CcCR versus F12-TZ for HNF is 7.5 cm-1. The closed-shell species are compared directly with the experiment, where a simpler (T)+EOM QFF using the aug-cc-pVTZ basis set compares more favorably than the more costly (T)+EOM/CcCR, suggesting a possible influence of decreasing accuracy with basis set size. Scans along internal coordinates are also provided which show reasonable modeling of the potential surface by (T)+EOM compared to benchmark QFFs computed for variationally accessible electronic states. The agreement between (T)+EOM/CcCR with F12-TZ and CcCR benchmarks is also shown to be quite accurate for rotational constants and geometries, with an MAE of 0.008 MHz for the rotational constants of (T)+EOM/CcCR versus CcCR for à 2A' HOO and agreement within 0.003 Šfor bond lengths.
Collapse
Affiliation(s)
- Megan C Davis
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| | - Ryan C Fortenberry
- Department of Chemistry & Biochemistry, University of Mississippi, University, Mississippi 38677-1848, United States
| |
Collapse
|
23
|
Mester D, Kállay M. A Simple Range-Separated Double-Hybrid Density Functional Theory for Excited States. J Chem Theory Comput 2021; 17:927-942. [PMID: 33400872 PMCID: PMC7884002 DOI: 10.1021/acs.jctc.0c01135] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2020] [Indexed: 01/12/2023]
Abstract
A simple and robust range-separated (RS) double-hybrid (DH) time-dependent density functional approach is presented for the accurate calculation of excitation energies of molecules within the Tamm-Dancoff approximation. The scheme can be considered as an excited-state extension of the ansatz proposed by Toulouse and co-workers [J. Chem. Phys. 2018, 148, 164105], which is based on the two-parameter decomposition of the Coulomb potential, for which both the exchange and correlation contributions are range-separated. A flexible and efficient implementation of the new scheme is also presented, which facilitates its extension to any combination of exchange and correlation functionals. The performance of the new approximation is tested for singlet excitations on several benchmark compilations and thoroughly compared to that of representative DH, RS hybrid, and RS DH functionals. The one-electron basis set dependence and computation times are also assessed. Our results show that the new approach improves on standard DHs in most cases, and it can provide a more robust and accurate alternative. In addition, on average, it noticeably surpasses the existing RS hybrid and RS DH functionals.
Collapse
Affiliation(s)
- Dávid Mester
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, Budapest, H-1521, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and
Materials Science, Budapest University of
Technology and Economics, P.O. Box 91, Budapest, H-1521, Hungary
| |
Collapse
|
24
|
Abstract
We present a new and efficient implementation of the closed shell coupled cluster singles and doubles with perturbative triples method (CC3) in the electronic structure program eT. Asymptotically, a ground state calculation has an iterative cost of 4nV4nO3 floating point operations (FLOP), where nV and nO are the number of virtual and occupied orbitals, respectively. The Jacobian and transpose Jacobian transformations, required to iteratively solve for excitation energies and transition moments, both require 8nV4nO3 FLOP. We have also implemented equation of motion (EOM) transition moments for CC3. The EOM transition densities require recalculation of triples amplitudes, as nV3nO3 tensors are not stored in memory. This results in a noniterative computational cost of 10nV4nO3 FLOP for the ground state density and 26nV4nO3 FLOP per state for the transition densities. The code is compared to the CC3 implementations in CFOUR, DALTON, and PSI4. We demonstrate the capabilities of our implementation by calculating valence and core excited states of l-proline.
Collapse
Affiliation(s)
- Alexander
C. Paul
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
| | - Rolf H. Myhre
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
| | - Henrik Koch
- Department
of Chemistry, Norwegian University of Science
and Technology, NTNU, 7491 Trondheim, Norway
- Scuola
Normale Superiore, Piazza dei Cavaleri 7, 56126 Pisa, Italy
| |
Collapse
|
25
|
Goerigk L, Casanova-Paéz M. The Trip to the Density Functional Theory Zoo Continues: Making a Case for Time-Dependent Double Hybrids for Excited-State Problems. Aust J Chem 2021. [DOI: 10.1071/ch20093] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
This account is written for general users of time-dependent density functional theory (TD-DFT) methods as well as chemists who are unfamiliar with the field. It includes a brief overview of conventional TD-DFT approaches and recommendations for applications to organic molecules based on our own experience. The main emphasis of this work, however, lies in providing the first in-depth review of time-dependent double-hybrid density functionals. They were first established in 2007 with very promising follow-up studies in the subsequent four years before developments or applications became scarce. The topic has regained more interest since 2017, and this account reviews those latest developments led by our group. These developments have shown unprecedented robustness for a variety of different types of electronic excitations when compared to more conventional TD-DFT methods. In particular, time-dependent double hybrids do not suffer from artificial ghost states and are able to reproduce exciton-coupled absorption spectra. Our latest methods include range separation and belong to the currently best TD-DFT methods for singlet-singlet excitations in organic molecules. While there is still room for improvement and further development in this space, we hope that this account encourages users to adjust their computational protocols to such new methods to provide more real-life testing and scenarios.
Collapse
|
26
|
Vilarrubias P. Electronic spectroscopy of some small anions containing S, N and O using CR-EOM-CCSD(T) method. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1797915] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Pere Vilarrubias
- Generalitat de Catalunya, Departament d’Ensenyament, IES Castelló d’Empúries, Girona, Spain
| |
Collapse
|
27
|
Braïda B, Chen Z, Wu W, Hiberty PC. Valence Bond Alternative Yielding Compact and Accurate Wave Functions for Challenging Excited States. Application to Ozone and Sulfur Dioxide. J Chem Theory Comput 2020; 17:330-343. [PMID: 33319998 DOI: 10.1021/acs.jctc.0c00598] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A novel state-averaged version of ab initio nonorthogonal valence bond method is described, for the sake of accurate theoretical studies of excited states in the valence bond framework. With respect to standard calculations in the molecular orbital framework, the state-averaged breathing-orbital valence bond (BOVB) method has the advantage to be free from the penalizing constraint for the ground and excited state(s) to share the same unique set of orbitals. The ability of the BOVB method to faithfully describe excited states and to compute accurate transition energies from the ground state is tested on the five lowest-lying singlet electronic states of ozone and sulfur dioxide, among which 11B2 and 21A1 are the challenging ones. As the 11A2, 11B1, and 11B2 states are of different symmetries than the ground state, they can be calculated at the state-specific BOVB level. On the other hand, the 21A1 states and the 11A1 ground states, which are of like symmetry, are calculated with the state-averaged BOVB technique. In all cases, the calculated vertical energies are close to the experimental values when available, and at par with the most sophisticated calculations in the molecular framework, despite the extreme compactness of the BOVB wave functions, made of no more than 5-9 valence bond structures in all cases. The features that allow the combination of compactness and accuracy in challenging cases are analyzed. For the "ionic" 11B2 states, which are the site of important charge fluctuations, it is because of the built-in dynamic correlation inherent to the BOVB method. For the 21A1 ones, this is the fact that these states have the degree of freedom of having different orbitals than the ground states, even though they are of like symmetry and calculated simultaneously using the newly implemented state-average BOVB algorithm. Finally, the description of the excited states in terms of Lewis structures is insightful, rationalizing the fast ring closure for the 21A1 state of ozone and predicting some diradical character in the so-called "ionic" 11B2 states.
Collapse
Affiliation(s)
- Benoît Braïda
- Laboratoire de Chimie Théorique, Sorbonne Université, UMR7616 CNRS, Paris 75252 France
| | - Zhenhua Chen
- College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and the State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, Fujian 361005, China
| | - Wei Wu
- College of Chemistry and Chemical Engineering, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, and the State Key Laboratory of Physical Chemistry of Solid Surfaces, Xiamen University, Xiamen, Fujian 361005, China
| | - Philippe C Hiberty
- Université Paris-Saclay, CNRS, Institut de Chimie Physique UMR8000, Orsay 91405, France
| |
Collapse
|
28
|
Lahiri J, Moemeni M, Kline J, Magoulas I, Yuwono SH, Laboe M, Shen J, Borhan B, Piecuch P, Jackson JE, Blanchard GJ, Dantus M. Isoenergetic two-photon excitation enhances solvent-to-solute excited-state proton transfer. J Chem Phys 2020; 153:224301. [PMID: 33317305 PMCID: PMC7725536 DOI: 10.1063/5.0020282] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2020] [Accepted: 11/15/2020] [Indexed: 01/05/2023] Open
Abstract
Two-photon excitation (TPE) is an attractive means for controlling chemistry in both space and time. Since isoenergetic one- and two-photon excitations (OPE and TPE) in non-centrosymmetric molecules are allowed to reach the same excited state, it is usually assumed that they produce similar excited-state reactivity. We compare the solvent-to-solute excited-state proton transfer of the super photobase FR0-SB following isoenergetic OPE and TPE. We find up to 62% increased reactivity following TPE compared to OPE. From steady-state spectroscopy, we rule out the involvement of different excited states and find that OPE and TPE spectra are identical in non-polar solvents but not in polar ones. We propose that differences in the matrix elements that contribute to the two-photon absorption cross sections lead to the observed enhanced isoenergetic reactivity, consistent with the predictions of our high-level coupled-cluster-based computational protocol. We find that polar solvent configurations favor greater dipole moment change between ground and excited states, which enters the probability for TPE as the absolute value squared. This, in turn, causes a difference in the Franck-Condon region reached via TPE compared to OPE. We conclude that a new method has been found for controlling chemical reactivity via the matrix elements that affect two-photon cross sections, which may be of great utility for spatial and temporal precision chemistry.
Collapse
Affiliation(s)
- Jurick Lahiri
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Mehdi Moemeni
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Jessica Kline
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Ilias Magoulas
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Stephen H. Yuwono
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Maryann Laboe
- Department of Chemical Engineering and Material Science, Michigan State University, East Lansing, Michigan 48824, USA
| | - Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Babak Borhan
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Piotr Piecuch
- Authors to whom correspondence should be addressed: , Tel.: +1-517-353-0501; , Tel.: +1-517-353-1151; , Tel.: +1-517-353-0504; , Tel.: +1-517-353-1105; and , Tel.: +1-517-353-1191
| | - James E. Jackson
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - G. J. Blanchard
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Marcos Dantus
- Authors to whom correspondence should be addressed: , Tel.: +1-517-353-0501; , Tel.: +1-517-353-1151; , Tel.: +1-517-353-0504; , Tel.: +1-517-353-1105; and , Tel.: +1-517-353-1191
| |
Collapse
|
29
|
Yuwono SH, Chakraborty A, Emiliano Deustua J, Shen J, Piecuch P. Accelerating convergence of equation-of-motion coupled-cluster computations using the semi-stochastic CC(P;Q) formalism. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1817592] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Stephen H. Yuwono
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
| | - Arnab Chakraborty
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
| | | | - Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, MI, USA
- Department of Physics & Astronomy, Michigan State University, East Lansing, MI, USA
| |
Collapse
|
30
|
Clune R, Shea JAR, Neuscamman E. N5-Scaling Excited-State-Specific Perturbation Theory. J Chem Theory Comput 2020; 16:6132-6141. [DOI: 10.1021/acs.jctc.0c00308] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Rachel Clune
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Jacqueline A. R. Shea
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
31
|
Matthews DA, Cheng L, Harding ME, Lipparini F, Stopkowicz S, Jagau TC, Szalay PG, Gauss J, Stanton JF. Coupled-cluster techniques for computational chemistry: The CFOUR program package. J Chem Phys 2020; 152:214108. [DOI: 10.1063/5.0004837] [Citation(s) in RCA: 214] [Impact Index Per Article: 53.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Affiliation(s)
- Devin A. Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, USA
| | - Lan Cheng
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Michael E. Harding
- Institut für Physikalische Chemie, Karlsruher Institut für Technologie (KIT), Kaiserstr. 12, D-76131 Karlsruhe, Germany
| | - Filippo Lipparini
- Dipartimento di Chimica e Chimica Industriale, Università di Pisa, Via G. Moruzzi 13, I-56124 Pisa, Italy
| | - Stella Stopkowicz
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - Thomas-C. Jagau
- Department of Chemistry, University of Munich (LMU), Butenandtstr. 5-13, D-81377 Munich, Germany
| | - Péter G. Szalay
- ELTE Eötvös Loránd University, Institute of Chemistry, Pázmány Péter sétány 1/A, H-1117 Budapest, Hungary
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, D-55128 Mainz, Germany
| | - John F. Stanton
- Quantum Theory Project, Departments of Chemistry and Physics, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
32
|
Aprà E, Bylaska EJ, de Jong WA, Govind N, Kowalski K, Straatsma TP, Valiev M, van Dam HJJ, Alexeev Y, Anchell J, Anisimov V, Aquino FW, Atta-Fynn R, Autschbach J, Bauman NP, Becca JC, Bernholdt DE, Bhaskaran-Nair K, Bogatko S, Borowski P, Boschen J, Brabec J, Bruner A, Cauët E, Chen Y, Chuev GN, Cramer CJ, Daily J, Deegan MJO, Dunning TH, Dupuis M, Dyall KG, Fann GI, Fischer SA, Fonari A, Früchtl H, Gagliardi L, Garza J, Gawande N, Ghosh S, Glaesemann K, Götz AW, Hammond J, Helms V, Hermes ED, Hirao K, Hirata S, Jacquelin M, Jensen L, Johnson BG, Jónsson H, Kendall RA, Klemm M, Kobayashi R, Konkov V, Krishnamoorthy S, Krishnan M, Lin Z, Lins RD, Littlefield RJ, Logsdail AJ, Lopata K, Ma W, Marenich AV, Martin Del Campo J, Mejia-Rodriguez D, Moore JE, Mullin JM, Nakajima T, Nascimento DR, Nichols JA, Nichols PJ, Nieplocha J, Otero-de-la-Roza A, Palmer B, Panyala A, Pirojsirikul T, Peng B, Peverati R, Pittner J, Pollack L, Richard RM, Sadayappan P, Schatz GC, Shelton WA, Silverstein DW, Smith DMA, Soares TA, Song D, Swart M, Taylor HL, Thomas GS, Tipparaju V, Truhlar DG, Tsemekhman K, Van Voorhis T, Vázquez-Mayagoitia Á, Verma P, Villa O, Vishnu A, Vogiatzis KD, Wang D, Weare JH, Williamson MJ, Windus TL, Woliński K, Wong AT, Wu Q, Yang C, Yu Q, Zacharias M, Zhang Z, Zhao Y, Harrison RJ. NWChem: Past, present, and future. J Chem Phys 2020; 152:184102. [PMID: 32414274 DOI: 10.1063/5.0004997] [Citation(s) in RCA: 302] [Impact Index Per Article: 75.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Abstract
Specialized computational chemistry packages have permanently reshaped the landscape of chemical and materials science by providing tools to support and guide experimental efforts and for the prediction of atomistic and electronic properties. In this regard, electronic structure packages have played a special role by using first-principle-driven methodologies to model complex chemical and materials processes. Over the past few decades, the rapid development of computing technologies and the tremendous increase in computational power have offered a unique chance to study complex transformations using sophisticated and predictive many-body techniques that describe correlated behavior of electrons in molecular and condensed phase systems at different levels of theory. In enabling these simulations, novel parallel algorithms have been able to take advantage of computational resources to address the polynomial scaling of electronic structure methods. In this paper, we briefly review the NWChem computational chemistry suite, including its history, design principles, parallel tools, current capabilities, outreach, and outlook.
Collapse
Affiliation(s)
- E Aprà
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - E J Bylaska
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - W A de Jong
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - N Govind
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - K Kowalski
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - T P Straatsma
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - M Valiev
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - H J J van Dam
- Brookhaven National Laboratory, Upton, New York 11973, USA
| | - Y Alexeev
- Argonne Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - J Anchell
- Intel Corporation, Santa Clara, California 95054, USA
| | - V Anisimov
- Argonne Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - F W Aquino
- QSimulate, Cambridge, Massachusetts 02139, USA
| | - R Atta-Fynn
- Department of Physics, The University of Texas at Arlington, Arlington, Texas 76019, USA
| | - J Autschbach
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - N P Bauman
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - J C Becca
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - D E Bernholdt
- Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | | | - S Bogatko
- 4G Clinical, Wellesley, Massachusetts 02481, USA
| | - P Borowski
- Faculty of Chemistry, Maria Curie-Skłodowska University in Lublin, 20-031 Lublin, Poland
| | - J Boschen
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| | - J Brabec
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, 18223 Prague 8, Czech Republic
| | - A Bruner
- Department of Chemistry and Physics, University of Tennessee at Martin, Martin, Tennessee 38238, USA
| | - E Cauët
- Service de Chimie Quantique et Photophysique (CP 160/09), Université libre de Bruxelles, B-1050 Brussels, Belgium
| | - Y Chen
- Facebook, Menlo Park, California 94025, USA
| | - G N Chuev
- Institute of Theoretical and Experimental Biophysics, Russian Academy of Science, Pushchino, Moscow Region 142290, Russia
| | - C J Cramer
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - J Daily
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - M J O Deegan
- SKAO, Jodrell Bank Observatory, Macclesfield SK11 9DL, United Kingdom
| | - T H Dunning
- Department of Chemistry, University of Washington, Seattle, Washington 98195, USA
| | - M Dupuis
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - K G Dyall
- Dirac Solutions, Portland, Oregon 97229, USA
| | - G I Fann
- Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - S A Fischer
- Chemistry Division, U. S. Naval Research Laboratory, Washington, DC 20375, USA
| | - A Fonari
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - H Früchtl
- EaStCHEM and School of Chemistry, University of St. Andrews, St. Andrews KY16 9ST, United Kingdom
| | - L Gagliardi
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - J Garza
- Departamento de Química, División de Ciencias Básicas e Ingeniería, Universidad Autónoma Metropolitana-Iztapalapa, Col. Vicentina, Iztapalapa, C.P. 09340 Ciudad de México, Mexico
| | - N Gawande
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - S Ghosh
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 5545, USA
| | - K Glaesemann
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - A W Götz
- San Diego Supercomputer Center, University of California, San Diego, La Jolla, California 92093, USA
| | - J Hammond
- Intel Corporation, Santa Clara, California 95054, USA
| | - V Helms
- Center for Bioinformatics, Saarland University, D-66041 Saarbrücken, Germany
| | - E D Hermes
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94551, USA
| | - K Hirao
- Next-generation Molecular Theory Unit, Advanced Science Institute, RIKEN, Saitama 351-0198, Japan
| | - S Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, USA
| | - M Jacquelin
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - L Jensen
- Department of Chemistry, The Pennsylvania State University, University Park, Pennsylvania 16802, USA
| | - B G Johnson
- Acrobatiq, Pittsburgh, Pennsylvania 15206, USA
| | - H Jónsson
- Faculty of Physical Sciences, University of Iceland, Reykjavík, Iceland and Department of Applied Physics, Aalto University, FI-00076 Aalto, Espoo, Finland
| | - R A Kendall
- Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - M Klemm
- Intel Corporation, Santa Clara, California 95054, USA
| | - R Kobayashi
- ANU Supercomputer Facility, Australian National University, Canberra, Australia
| | - V Konkov
- Chemistry Program, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - S Krishnamoorthy
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - M Krishnan
- Facebook, Menlo Park, California 94025, USA
| | - Z Lin
- Department of Physics, University of Science and Technology of China, Hefei, China
| | - R D Lins
- Aggeu Magalhaes Institute, Oswaldo Cruz Foundation, Recife, Brazil
| | | | - A J Logsdail
- Cardiff Catalysis Institute, School of Chemistry, Cardiff University, Cardiff, Wales CF10 3AT, United Kingdom
| | - K Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | - W Ma
- Institute of Software, Chinese Academy of Sciences, Beijing, China
| | - A V Marenich
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - J Martin Del Campo
- Departamento de Física y Química Teórica, Facultad de Química, Universidad Nacional Autónoma de México, México City, Mexico
| | - D Mejia-Rodriguez
- Quantum Theory Project, Department of Physics, University of Florida, Gainesville, Florida 32611, USA
| | - J E Moore
- Intel Corporation, Santa Clara, California 95054, USA
| | - J M Mullin
- DCI-Solutions, Aberdeen Proving Ground, Maryland 21005, USA
| | - T Nakajima
- Computational Molecular Science Research Team, RIKEN Center for Computational Science, Kobe, Hyogo 650-0047, Japan
| | - D R Nascimento
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - J A Nichols
- Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, USA
| | - P J Nichols
- Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA
| | - J Nieplocha
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - A Otero-de-la-Roza
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, 33006 Oviedo, Spain
| | - B Palmer
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - A Panyala
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - T Pirojsirikul
- Department of Chemistry, Prince of Songkla University, Hat Yai, Songkhla 90112, Thailand
| | - B Peng
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - R Peverati
- Chemistry Program, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - J Pittner
- J. Heyrovský Institute of Physical Chemistry, Academy of Sciences of the Czech Republic, v.v.i., 18223 Prague 8, Czech Republic
| | - L Pollack
- StudyPoint, Boston, Massachusetts 02114, USA
| | | | - P Sadayappan
- School of Computing, University of Utah, Salt Lake City, Utah 84112, USA
| | - G C Schatz
- Department of Chemistry, Northwestern University, Evanston, Illinois 60208, USA
| | - W A Shelton
- Cain Department of Chemical Engineering, Louisiana State University, Baton Rouge, Louisiana 70803, USA
| | | | - D M A Smith
- Intel Corporation, Santa Clara, California 95054, USA
| | - T A Soares
- Dept. of Fundamental Chemistry, Universidade Federal de Pernambuco, Recife, Brazil
| | - D Song
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - M Swart
- ICREA, 08010 Barcelona, Spain and Universitat Girona, Institut de Química Computacional i Catàlisi, Campus Montilivi, 17003 Girona, Spain
| | - H L Taylor
- CD-adapco/Siemens, Melville, New York 11747, USA
| | - G S Thomas
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - V Tipparaju
- Cray Inc., Bloomington, Minnesota 55425, USA
| | - D G Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | | | - T Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Á Vázquez-Mayagoitia
- Argonne Leadership Computing Facility, Argonne National Laboratory, Argonne, Illinois 60439, USA
| | - P Verma
- 1QBit, Vancouver, British Columbia V6E 4B1, Canada
| | - O Villa
- NVIDIA, Santa Clara, California 95051, USA
| | - A Vishnu
- Pacific Northwest National Laboratory, Richland, Washington 99352, USA
| | - K D Vogiatzis
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996, USA
| | - D Wang
- College of Physics and Electronics, Shandong Normal University, Jinan, Shandong 250014, China
| | - J H Weare
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - M J Williamson
- Department of Chemistry, Cambridge University, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - T L Windus
- Department of Chemistry, Iowa State University and Ames Laboratory, Ames, Iowa 50011, USA
| | - K Woliński
- Faculty of Chemistry, Maria Curie-Skłodowska University in Lublin, 20-031 Lublin, Poland
| | - A T Wong
- Qwil, San Francisco, California 94107, USA
| | - Q Wu
- Brookhaven National Laboratory, Upton, New York 11973, USA
| | - C Yang
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Q Yu
- AMD, Santa Clara, California 95054, USA
| | - M Zacharias
- Department of Physics, Technical University of Munich, 85748 Garching, Germany
| | - Z Zhang
- Stanford Research Computing Center, Stanford University, Stanford, California 94305, USA
| | - Y Zhao
- State Key Laboratory of Silicate Materials for Architectures, International School of Materials Science and Engineering, Wuhan University of Technology, Wuhan 430070, China
| | - R J Harrison
- Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, USA
| |
Collapse
|
33
|
Barca GMJ, Bertoni C, Carrington L, Datta D, De Silva N, Deustua JE, Fedorov DG, Gour JR, Gunina AO, Guidez E, Harville T, Irle S, Ivanic J, Kowalski K, Leang SS, Li H, Li W, Lutz JJ, Magoulas I, Mato J, Mironov V, Nakata H, Pham BQ, Piecuch P, Poole D, Pruitt SR, Rendell AP, Roskop LB, Ruedenberg K, Sattasathuchana T, Schmidt MW, Shen J, Slipchenko L, Sosonkina M, Sundriyal V, Tiwari A, Galvez Vallejo JL, Westheimer B, Włoch M, Xu P, Zahariev F, Gordon MS. Recent developments in the general atomic and molecular electronic structure system. J Chem Phys 2020; 152:154102. [PMID: 32321259 DOI: 10.1063/5.0005188] [Citation(s) in RCA: 518] [Impact Index Per Article: 129.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A discussion of many of the recently implemented features of GAMESS (General Atomic and Molecular Electronic Structure System) and LibCChem (the C++ CPU/GPU library associated with GAMESS) is presented. These features include fragmentation methods such as the fragment molecular orbital, effective fragment potential and effective fragment molecular orbital methods, hybrid MPI/OpenMP approaches to Hartree-Fock, and resolution of the identity second order perturbation theory. Many new coupled cluster theory methods have been implemented in GAMESS, as have multiple levels of density functional/tight binding theory. The role of accelerators, especially graphical processing units, is discussed in the context of the new features of LibCChem, as it is the associated problem of power consumption as the power of computers increases dramatically. The process by which a complex program suite such as GAMESS is maintained and developed is considered. Future developments are briefly summarized.
Collapse
Affiliation(s)
- Giuseppe M J Barca
- Research School of Computer Science, Australian National University, Canberra, ACT 2601, Australia
| | - Colleen Bertoni
- Argonne Leadership Computing Facility, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Laura Carrington
- EP Analytics, 12121 Scripps Summit Dr. Ste. 130, San Diego, California 92131, USA
| | - Dipayan Datta
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Nuwan De Silva
- Department of Physical and Biological Sciences, Western New England University, Springfield, Massachusetts 01119, USA
| | - J Emiliano Deustua
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Dmitri G Fedorov
- Research Center for Computational Design of Advanced Functional Materials (CD-FMat), National Institute of Advanced Industrial Science and Technology (AIST), Umezono 1-1-1, Tsukuba 305-8568, Japan
| | - Jeffrey R Gour
- Microsoft, 15590 NE 31st St., Redmond, Washington 98052, USA
| | - Anastasia O Gunina
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Emilie Guidez
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217, USA
| | - Taylor Harville
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Stephan Irle
- Computational Science and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37830, USA
| | - Joe Ivanic
- Advanced Biomedical Computational Science, Frederick National Laboratory for Cancer Research, Frederick, Maryland 21702, USA
| | - Karol Kowalski
- Physical Sciences Division, Battelle, Pacific Northwest National Laboratory, K8-91, P.O. Box 999, Richland, Washington 99352, USA
| | - Sarom S Leang
- EP Analytics, 12121 Scripps Summit Dr. Ste. 130, San Diego, California 92131, USA
| | - Hui Li
- Department of Chemistry, University of Nebraska, Lincoln, Nebraska 68588, USA
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, People's Republic of China
| | - Jesse J Lutz
- Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87185, USA
| | - Ilias Magoulas
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Joani Mato
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Vladimir Mironov
- Department of Chemistry, Lomonosov Moscow State University, Leninskie Gory 1/3, Moscow 119991, Russian Federation
| | - Hiroya Nakata
- Kyocera Corporation, Research Institute for Advanced Materials and Devices, 3-5-3 Hikaridai Seika-cho, Souraku-gun, Kyoto 619-0237, Japan
| | - Buu Q Pham
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - David Poole
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Spencer R Pruitt
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Alistair P Rendell
- Research School of Computer Science, Australian National University, Canberra, ACT 2601, Australia
| | - Luke B Roskop
- Cray Inc., a Hewlett Packard Enterprise Company, 2131 Lindau Ln #1000, Bloomington, Minnesota 55425, USA
| | - Klaus Ruedenberg
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | | | - Michael W Schmidt
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Lyudmila Slipchenko
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Masha Sosonkina
- Department of Computational Modeling and Simulation Engineering, Old Dominion University, Norfolk, Virginia 23529, USA
| | - Vaibhav Sundriyal
- Department of Computational Modeling and Simulation Engineering, Old Dominion University, Norfolk, Virginia 23529, USA
| | - Ananta Tiwari
- EP Analytics, 12121 Scripps Summit Dr. Ste. 130, San Diego, California 92131, USA
| | - Jorge L Galvez Vallejo
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Bryce Westheimer
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Marta Włoch
- 530 Charlesina Dr., Rochester, Michigan 48306, USA
| | - Peng Xu
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Federico Zahariev
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Mark S Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| |
Collapse
|
34
|
Yuwono SH, Magoulas I, Piecuch P. Quantum computation solves a half-century-old enigma: Elusive vibrational states of magnesium dimer found. SCIENCE ADVANCES 2020; 6:eaay4058. [PMID: 32284970 PMCID: PMC7124942 DOI: 10.1126/sciadv.aay4058] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/15/2019] [Accepted: 01/09/2020] [Indexed: 05/29/2023]
Abstract
The high-lying vibrational states of the magnesium dimer (Mg2), which has been recognized as an important system in studies of ultracold and collisional phenomena, have eluded experimental characterization for half a century. Until now, only the first 14 vibrational states of Mg2 have been experimentally resolved, although it has been suggested that the ground-state potential may support five additional levels. Here, we present highly accurate ab initio potential energy curves based on state-of-the-art coupled-cluster and full configuration interaction computations for the ground and excited electronic states involved in the experimental investigations of Mg2. Our ground-state potential unambiguously confirms the existence of 19 vibrational levels, with ~1 cm-1 root mean square deviation between the calculated rovibrational term values and the available experimental and experimentally derived data. Our computations reproduce the latest laser-induced fluorescence spectrum and provide guidance for the experimental detection of the previously unresolved vibrational levels.
Collapse
Affiliation(s)
- Stephen H. Yuwono
- Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA
| | - Ilias Magoulas
- Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, MI 48824, USA
- Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA
| |
Collapse
|
35
|
|
36
|
Ma F, Wang Z, Guo M, Wang F. Approximate equation-of-motion coupled-cluster methods for electron affinities of closed-shell molecules. J Chem Phys 2020; 152:124111. [PMID: 32241115 DOI: 10.1063/1.5142736] [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 investigate performance of the equation-of-motion coupled-cluster method at the single and doubles level (EOM-CCSD) and a series of approximate methods based on EOM-CCSD on electron affinities (EA) of closed-shell cations and neutral molecules with positive and negative EAs in this work. Our results confirm that P-EOM-MBPT2 can provide reasonable EAs when molecules with significant multireference character are not considered and its mean absolute error on EAs of these molecules is around or less than 0.2 eV. Its accuracy is comparable to that of the more expensive EOM-CCSD(2) method. Results of EOM-CCSD(2), P-EOM-MBPT2, and CIS(D∞) indicate that the [[H, ac +], T2] term in the 1h2p-1h block is more important on EAs than the term neglected in the 1h2p-1h2p block in P-EOM-MBPT2. We proposed an economical method where EAs from CIS(D∞) are corrected by treating this [[H, ac +], T2] term in the 1h2p-1h block perturbatively [corr-CIS(D∞)]. EAs with corr-CIS(D∞) agree very well with those of P-EOM-MBPT2 with a difference of less than 0.02 eV. Computational scaling of this method is N4 for the iterative part and N5 for some non-iterative steps. Its storage requirement is only of OV3. Corr-CIS(D∞) is an economical and reliable method on EAs, and it can be applied to EAs of large molecules.
Collapse
Affiliation(s)
- Fengjiao Ma
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610065, People's Republic of China
| | - Zhifan Wang
- School of Electronic Engineering, Chengdu Technological University, Chengdu 611730, People's Republic of China
| | - Minggang Guo
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610065, People's Republic of China
| | - Fan Wang
- Institute of Atomic and Molecular Physics, Key Laboratory of High Energy Density Physics and Technology, Ministry of Education, Sichuan University, Chengdu 610065, People's Republic of China
| |
Collapse
|
37
|
Lewis AM, Berkelbach TC. Ab Initio Linear and Pump-Probe Spectroscopy of Excitons in Molecular Crystals. J Phys Chem Lett 2020; 11:2241-2246. [PMID: 32109074 DOI: 10.1021/acs.jpclett.0c00031] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Linear and nonlinear spectroscopies are powerful tools used to investigate the energetics and dynamics of electronic excited states of both molecules and crystals. While highly accurate ab initio calculations of molecular spectra can be performed relatively routinely, extending these calculations to periodic systems is challenging. Here, we present calculations of the linear absorption spectrum and pump-probe two-photon photoemission spectra of the naphthalene crystal using equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD). Molecular acene crystals are of interest due to the low-energy multiexciton singlet states they exhibit, which have been studied extensively as intermediates involved in singlet fission. Our linear absorption spectrum is in good agreement with experiment, predicting a first exciton absorption peak at 4.4 eV, and our two-photon photoemission spectra capture the qualitative behavior of multiexciton states, whose double-excitation character cannot be captured by current methods. The simulated pump-probe spectra provide support for existing interpretations of two-photon photoemission experiments in closely related acene crystals such as tetracene and pentacene.
Collapse
Affiliation(s)
- Alan M Lewis
- Department of Chemistry and James Franck Institute, University of Chicago, Chicago, Illinois 60637, United States
| | - Timothy C Berkelbach
- Department of Chemistry, Columbia University, New York, New York 10027, United States
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| |
Collapse
|
38
|
Scaling exchange and correlation in the on-top density functional of multiconfiguration pair-density functional theory: effect on electronic excitation energies and bond energies. Theor Chem Acc 2020. [DOI: 10.1007/s00214-019-2539-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/25/2022]
|
39
|
Pokhilko P, Izmodenov D, Krylov AI. Extension of frozen natural orbital approximation to open-shell references: Theory, implementation, and application to single-molecule magnets. J Chem Phys 2020; 152:034105. [DOI: 10.1063/1.5138643] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Affiliation(s)
- Pavel Pokhilko
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| | - Daniil Izmodenov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
- Department of Chemistry, Lomonosov Moscow State University, Moscow 119991, Russia
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089-0482, USA
| |
Collapse
|
40
|
Jahangiri S, Arrazola JM, Quesada N, Delgado A. Quantum algorithm for simulating molecular vibrational excitations. Phys Chem Chem Phys 2020; 22:25528-25537. [DOI: 10.1039/d0cp03593a] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
We introduce a quantum algorithm for simulating molecular vibrational excitations during vibronic transitions. The algorithm is used to simulate vibrational excitations of pyrrole and butane during photochemical and mechanochemical excitations.
Collapse
|
41
|
Hong Y, Hou B, Jiang H, Zhang J. Machine learning and artificial neural network accelerated computational discoveries in materials science. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1450] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Yang Hong
- Department of Chemistry University of Nebraska‐Lincoln Lincoln Nebraska
| | - Bo Hou
- Department of Engineering University of Cambridge Cambridge UK
| | - Hengle Jiang
- Holland Computing Center University of Nebraska‐Lincoln Lincoln Nebraska
| | - Jingchao Zhang
- Holland Computing Center University of Nebraska‐Lincoln Lincoln Nebraska
| |
Collapse
|
42
|
Musiał M, Meissner L, Cembrzynska J. The intermediate Hamiltonian Fock-space coupled-cluster method with approximate evaluation of the three-body effects. J Chem Phys 2019; 151:184102. [DOI: 10.1063/1.5124806] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Monika Musiał
- Institute of Chemistry, University of Silesia in Katowice, Szkolna 9, 40-006 Katowice, Poland
| | - Leszek Meissner
- Institute of Physics, Nicholaus Copernicus University, Grudziadzka 5/7, 87-100 Toruń, Poland
| | - Justyna Cembrzynska
- Institute of Physics, Nicholaus Copernicus University, Grudziadzka 5/7, 87-100 Toruń, Poland
| |
Collapse
|
43
|
Zhang L, Truhlar DG, Sun S. Full-dimensional three-state potential energy surfaces and state couplings for photodissociation of thiophenol. J Chem Phys 2019; 151:154306. [DOI: 10.1063/1.5124870] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- Linyao Zhang
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People’s Republic of China
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Donald G. Truhlar
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Shaozeng Sun
- School of Energy Science and Engineering, Harbin Institute of Technology, Harbin 150001, People’s Republic of China
| |
Collapse
|
44
|
Lahiri J, Moemeni M, Kline J, Borhan B, Magoulas I, Yuwono SH, Piecuch P, Jackson JE, Dantus M, Blanchard GJ. Proton Abstraction Mediates Interactions between the Super Photobase FR0-SB and Surrounding Alcohol Solvent. J Phys Chem B 2019; 123:8448-8456. [PMID: 31532676 DOI: 10.1021/acs.jpcb.9b06580] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We report on the motional and proton transfer dynamics of the super photobase FR0-SB in the series of normal alcohols C1 (methanol) through C8 (n-octanol) and ethylene glycol. Steady-state and time-resolved fluorescence data reveal that the proton abstraction dynamics of excited FR0-SB depend on the identity of the solvent and that the transfer of the proton from solvent to FR0-SB*, forming FR0-HSB+*, fundamentally alters the nature of interactions between the excited molecule and its surroundings. In its unprotonated state, solvent interactions with FR0-SB* are consistent with slip limit behavior, and in its protonated form, intermolecular interactions are consistent with a much stronger interaction of FR0-HSB+* with the deprotonated solvent RO-. We understand the excited-state population dynamics in the context of a kinetic model involving a transition state wherein FR0-HSB+* is still bound to the negatively charged alkoxide, prior to solvation of the two charged species. Data acquired in ethylene glycol confirm the hypothesis that the rotational diffusion dynamics of FR0-SB* are largely mediated by solvent viscosity while proton transfer dynamics are mediated by the lifetime of the transition state. Taken collectively, our results demonstrate that FR0-SB* extracts solvent protons efficiently and in a predictable manner, consistent with a ca. 3-fold increase in dipole moment upon photoexcitation as determined by ab initio calculations based on the equation-of-motion coupled-cluster theory.
Collapse
|
45
|
Wang T, Brédas J. Quantum‐Chemical Evaluation of Impact of Chlorination versus Fluorination on the Electronic Properties of Donors and Acceptors for Organic Solar Cells. ADVANCED THEORY AND SIMULATIONS 2019. [DOI: 10.1002/adts.201900136] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Tonghui Wang
- School of Chemistry and Biochemistry & Center for Organic Photonics and Electronics Georgia Institute of Technology Atlanta GA 30332–0400 USA
| | - Jean‐Luc Brédas
- School of Chemistry and Biochemistry & Center for Organic Photonics and Electronics Georgia Institute of Technology Atlanta GA 30332–0400 USA
| |
Collapse
|
46
|
Cheng L. A study of non-iterative triples contributions in relativistic equation-of-motion coupled-cluster calculations using an exact two-component Hamiltonian with atomic mean-field spin-orbit integrals: Application to uranyl and other heavy-element compounds. J Chem Phys 2019; 151:104103. [DOI: 10.1063/1.5113796] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
47
|
Casanova-Páez M, Dardis MB, Goerigk L. ωB2PLYP and ωB2GPPLYP: The First Two Double-Hybrid Density Functionals with Long-Range Correction Optimized for Excitation Energies. J Chem Theory Comput 2019; 15:4735-4744. [PMID: 31298850 DOI: 10.1021/acs.jctc.9b00013] [Citation(s) in RCA: 86] [Impact Index Per Article: 17.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Double-hybrid density functionals are currently the most accurate density functionals for ground-state properties and electronic excitations. Nevertheless, the lack of a long-range correction scheme makes them unreliable when it comes to long-range excitations. For this reason, we propose the first two time-dependent double-hybrid functionals with correct asymptotic long-range behavior named ωB2PLYP and ωB2GPPLYP. Herein, we demonstrate their excellent performance and show that they are the most accurate and robust time-dependent density functional theory methods for electronic excitation energies. They provide a balanced description of local-valence, Rydberg, and charge-transfer states. They are also able to tackle the difficult first two transitions in polycyclic aromatic hydrocarbons and show very promising results in a preliminary study on transition-metal compounds, exemplified for titanium dioxide clusters. This work shows that double hybrids can be systematically improved also for excitation energies, and further work in this field is warranted.
Collapse
Affiliation(s)
| | - Michael B Dardis
- School of Chemistry , The University of Melbourne , Victoria 3010 , Australia
| | - Lars Goerigk
- School of Chemistry , The University of Melbourne , Victoria 3010 , Australia
| |
Collapse
|
48
|
Zhang QH, Liu YZ, Yin WY, Zhu RS, Yan YH, Ding PF, Bumaliya A. Physical properties and stepwise dissociation of Halon 2402 under external electric field. Chem Phys Lett 2019. [DOI: 10.1016/j.cplett.2019.03.037] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
49
|
Pitzer M, Ozga C, Küstner-Wetekam C, Reiß P, Knie A, Ehresmann A, Jahnke T, Giuliani A, Nahon L. State-Dependent Fragmentation of Protonated Uracil and Uridine. J Phys Chem A 2019; 123:3551-3557. [PMID: 30943036 DOI: 10.1021/acs.jpca.9b01822] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Action spectroscopy using photon excitation in the VUV range (photon energy 4.5-9 eV) was performed on protonated uracil (UraH+) and uridine (UrdH+). The precursor ions with m/ z 113 and m/ z 245, respectively, were produced by an electrospray ionization source and accumulated inside a quadrupole ion trap mass spectrometer. After irradiation with tunable synchrotron radiation, product ion mass spectra were obtained. Fragment yields as a function of excitation energy show several maxima that can be attributed to the photoexcitation into different electronic states. For uracil, vertically excited states were calculated using the equation-of-motion coupled cluster approach and compared to the observed maxima. This allows to establish correlations between electronic states and the resulting fragment masses and can thus help to disentangle the complex de-excitation and fragmentation pathways of nucleic acid building blocks. Photofragmentation of the nucleoside uridine shows a significantly lower variety of fragments, indicating stabilization of the nucleobase by the attached sugar.
Collapse
Affiliation(s)
- Martin Pitzer
- Department of Chemical and Biological Physics , Weizmann Institute of Science , Rehovot 7610001 , Israel.,Institute of Physics and Center for Interdisciplinary Nanostructure Science and Technology (CINSaT) , University of Kassel , 34132 Kassel , Germany
| | - Christian Ozga
- Institute of Physics and Center for Interdisciplinary Nanostructure Science and Technology (CINSaT) , University of Kassel , 34132 Kassel , Germany
| | - Catmarna Küstner-Wetekam
- Institute of Physics and Center for Interdisciplinary Nanostructure Science and Technology (CINSaT) , University of Kassel , 34132 Kassel , Germany
| | - Philipp Reiß
- Institute of Physics and Center for Interdisciplinary Nanostructure Science and Technology (CINSaT) , University of Kassel , 34132 Kassel , Germany
| | - André Knie
- Institute of Physics and Center for Interdisciplinary Nanostructure Science and Technology (CINSaT) , University of Kassel , 34132 Kassel , Germany
| | - Arno Ehresmann
- Institute of Physics and Center for Interdisciplinary Nanostructure Science and Technology (CINSaT) , University of Kassel , 34132 Kassel , Germany
| | - Till Jahnke
- Institute for Nuclear Physics , Goethe-University , 60438 Frankfurt , Germany
| | | | | |
Collapse
|
50
|
Deustua JE, Yuwono SH, Shen J, Piecuch P. Accurate excited-state energetics by a combination of Monte Carlo sampling and equation-of-motion coupled-cluster computations. J Chem Phys 2019; 150:111101. [DOI: 10.1063/1.5090346] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- J. Emiliano Deustua
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Stephen H. Yuwono
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Jun Shen
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
| | - Piotr Piecuch
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, USA
- Department of Physics and Astronomy, Michigan State University, East Lansing, Michigan 48824, USA
| |
Collapse
|