1
|
Arias-Martinez JE, Wu H, Head-Gordon M. Generalization of One-Center Nonorthogonal Configuration Interaction Singles to Open-Shell Singlet Reference States: Theory and Application to Valence-Core Pump-Probe States in Acetylacetone. J Chem Theory Comput 2024; 20:752-766. [PMID: 38164934 DOI: 10.1021/acs.jctc.3c01139] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2024]
Abstract
We formulate a one-center nonorthogonal configuration interaction singles (1C-NOCIS) theory for the computation of core excited states of an initial singlet state with two unpaired electrons. This model, which we refer to as 1C-NOCIS two-electron open-shell (2eOS), is appropriate for computing the K-edge near-edge X-ray absorption spectra (NEXAS) of the valence excited states of closed-shell molecules relevant to pump-probe time-resolved (TR) NEXAS experiments. With the inclusion of core-hole relaxation effects and explicit spin adaptation, 1C-NOCIS 2eOS requires mild shifts to match experiment, is free of artifacts due to spin contamination, and can capture the high-energy region of the spectrum beyond the transitions into the singly occupied molecular orbitals (SOMOs). Calculations on water and thymine illustrate the different key features of excited-state NEXAS, namely, the core-to-SOMO transitions as well as shifts and spin-splittings in the transitions analogous to those of the ground state. Simulations of the TR-NEXAS of acetylacetone after excitation to its π → π* singlet excited state at the carbon K-edge, an experiment carried out recently, showcase the ability of 1C-NOCIS 2eOS to efficiently simulate NEXAS based on nonadiabatic molecular dynamics simulations.
Collapse
Affiliation(s)
- Juan E Arias-Martinez
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Hamlin Wu
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
2
|
Hirao K, Nakajima T, Chan B. Core-Level 2s and 2p Binding Energies of Third-Period Elements (P, S, and Cl) Calculated by Hartree-Fock and Kohn-Sham Δ SCF Theory. J Phys Chem A 2023; 127:7954-7963. [PMID: 37703090 DOI: 10.1021/acs.jpca.3c04783] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/14/2023]
Abstract
In the present study, we investigate the use of the ΔSCF method and Slater's transition state (STS) theory to calculate the binding energies of the 2s and 2p electrons of third-period elements (P, S, and Cl). Both the Hartree-Fock (HF) and Kohn-Sham (KS) approximations are examined. The STS approximation performs well in reproducing the ΔSCF values. However, for the ΔSCF method itself, while the binding energy of the 2p electrons is accurately predicted, the results for 2s are fairly sensitive to the functional, exhibiting significant variations due to self-interaction errors (SIE). Nonetheless, the variations in chemical shifts between different species remain relatively small, and the values agree with experiments due to the cancellation of SIE. A notable observation is that the chemical shifts of the 2s and 2p electrons are similar, indicating a perturbation caused by the valence electrons. The error in the absolute binding energy of KS ΔSCF against the experiment is nearly constant for the same element in different molecules, and it depends largely on the functional owing to SIE. A shifting scheme previously developed can be employed to reproduce the experimental 2s and 2p binding energies, with dependence on the functional and atom but not on the molecule even for 2s KS ΔSCF binding energies. Upon obtaining the corrected binding energies, we find that the gap between 2s and 2p binding energy is nearly independent of chemical environment for a given element: 57.5, 63.9, and 70.9 eV for the elements P, S, and Cl, respectively.
Collapse
Affiliation(s)
- Kimihiko Hirao
- Fukui Institute for Fundamental Chemistry, Kyoto University, Takano, Nishihiraki-cho 34-4, Sakyo-ku, Kyoto 606-8103, Japan
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan
| | - Takahito Nakajima
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan
| | - Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| |
Collapse
|
3
|
Huang M, Evangelista FA. A study of core-excited states of organic molecules computed with the generalized active space driven similarity renormalization group. J Chem Phys 2023; 158:124112. [PMID: 37003756 DOI: 10.1063/5.0137096] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/03/2023] Open
Abstract
This work examines the accuracy and precision of x-ray absorption spectra computed with a multireference approach that combines generalized active space (GAS) references with the driven similarity renormalization group (DSRG). We employ the x-ray absorption benchmark of organic molecule (XABOOM) set, consisting of 116 transitions from mostly organic molecules [Fransson et al., J. Chem. Theory Comput. 17, 1618 (2021)]. Several approximations to a full-valence active space are examined and benchmarked. Absolute excitation energies and intensities computed with the GAS-DSRG truncated to second-order in perturbation theory are found to systematically underestimate experimental and reference theoretical values. Third-order perturbative corrections significantly improve the accuracy of GAS-DSRG absolute excitation energies, bringing the mean absolute deviation from experimental values down to 0.32 eV. The ozone molecule and glyoxylic acid are particularly challenging for second-order perturbation theory and are examined in detail to assess the importance of active space truncation and intruder states.
Collapse
Affiliation(s)
- Meng Huang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, USA
| |
Collapse
|
4
|
Hirao K, Nakajima T, Chan B, Lee HJ. The core ionization energies calculated by delta SCF and Slater's transition state theory. J Chem Phys 2023; 158:064112. [PMID: 36792520 DOI: 10.1063/5.0140032] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
Abstract
The core ionization energies of the second-period and third-period elements are studied by ΔSCF and Slater's transition state (STS) theory by using Hartree-Fock (HF) and Kohn-Sham (KS) approximations. Electron correlation increases the estimated core ionization energies, while the self-interaction error (SIE) decreases them, especially for the third-period elements and is a more significant factor. As a result, while HF lacks electron correlation, it is free of SIE and reasonably predicts the core ionization energies. The core ionization energies calculated by HF STS are very close to those calculated by HF ΔSCF, showing that STS reasonably describes the relaxation of the core hole. The core ionization energies calculated by KS are particularly sensitive to the SIE of the functional used, with functionals having less SIE yielding more accurate ΔSCF core ionization energies. Consequently, BH&HLYP gives better results than B3LYP and LC-BOP since BH&HLYP is the hybrid functional with high proportion of the exact HF exchange. Although the core ionization energies are underestimated by ΔSCF due to SIE, STS gives larger core ionization energies than ΔSCF due to a concave behavior of the error curves of STS, which is also related to SIE. The mean absolute deviations of STS relative to ΔSCF, and relative to the experiment, are almost constant regardless of the nuclei among the element in the second period, and likewise among those in the third period. The systematic nature suggests that shifting the STS core ionization energies may be useful. We propose the shifted STS (1) for reproducing ΔSCF values, and the shifted STS (2) to reproduce the observed ones for KS calculations. Both schemes work quite well. The calculated results of KS ΔSCF and STS vary depending on the functional. However, the variation of each species' shifted STS (2) is very small, and all shifted STS (2) values are close to the observed ones. As the shifted STS require only one SCF calculation, they are simple and practical for predicting the core ionization energies.
Collapse
Affiliation(s)
- Kimihiko Hirao
- Fukui Institute for Fundamental Chemistry, Kyoto University, Takano, Nishihiraki-cho 34-4, Sakyo-ku, Kyoto 606-8103, Japan
| | - Takahito Nakajima
- RIKEN Center for Computational Science, 7-1-26, Minatojima-minami-machi, Chuo-ku, Kobe 650-0047, Japan
| | - Bun Chan
- Graduate School of Engineering, Nagasaki University, Bunkyo 1-14, Nagasaki 852-8521, Japan
| | - Ho-Jin Lee
- Department of Natural Sciences, Southwest Tennessee Community College, Memphis, Tennessee 38015, USA
| |
Collapse
|
5
|
Wibowo M, Huynh BC, Cheng CY, Irons TJP, Teale AM. Understanding ground and excited-state molecular structure in strong magnetic fields using the maximum overlap method. Mol Phys 2022. [DOI: 10.1080/00268976.2022.2152748] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Affiliation(s)
- Meilani Wibowo
- School of Chemistry, University of Nottingham, University Park, Nottingham, UK
| | - Bang C. Huynh
- School of Chemistry, University of Nottingham, University Park, Nottingham, UK
| | - Chi Y. Cheng
- School of Chemistry, University of Nottingham, University Park, Nottingham, UK
| | - Tom J. P. Irons
- School of Chemistry, University of Nottingham, University Park, Nottingham, UK
| | - Andrew M. Teale
- School of Chemistry, University of Nottingham, University Park, Nottingham, UK
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Blindern, Oslo, Norway
| |
Collapse
|
6
|
Burton HGA. Generalized nonorthogonal matrix elements. II: Extension to arbitrary excitations. J Chem Phys 2022; 157:204109. [DOI: 10.1063/5.0122094] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022] Open
Abstract
Electronic structure methods that exploit nonorthogonal Slater determinants face the challenge of efficiently computing nonorthogonal matrix elements. In a recent publication [H. G. A. Burton, J. Chem. Phys. 154, 144109 (2021)], I introduced a generalized extension to the nonorthogonal Wick’s theorem that allows matrix elements to be derived between excited configurations from a pair of reference determinants with a singular nonorthogonal orbital overlap matrix. However, that work only provided explicit expressions for one- and two-body matrix elements between singly- or doubly-excited configurations. Here, this framework is extended to compute generalized nonorthogonal matrix elements between higher-order excitations. Pre-computing and storing intermediate values allows one- and two-body matrix elements to be evaluated with an [Formula: see text] scaling relative to the system size, and the LIBGNME computational library is introduced to achieve this in practice. These advances make the evaluation of all nonorthogonal matrix elements almost as easy as their orthogonal counterparts, facilitating a new phase of development in nonorthogonal electronic structure theory.
Collapse
Affiliation(s)
- Hugh G. A. Burton
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
7
|
Carter-Fenk K, Cunha LA, Arias-Martinez JE, Head-Gordon M. Electron-Affinity Time-Dependent Density Functional Theory: Formalism and Applications to Core-Excited States. J Phys Chem Lett 2022; 13:9664-9672. [PMID: 36215404 DOI: 10.1021/acs.jpclett.2c02564] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
The lack of particle-hole attraction and orbital relaxation within time-dependent density functional theory (TDDFT) lead to extreme errors in the prediction of K-edge X-ray absorption spectra (XAS). We derive a linear-response formalism that uses optimized orbitals of the n - 1-electron system as the reference, building orbital relaxation and a proper hole into the initial density. Our approach is an exact generalization of the static-exchange approximation that ameliorates the particle-hole interaction error associated with the adiabatic approximation and reduces errors in TDDFT XAS by orders of magnitude. With a statistical performance of just 0.5 eV root-mean-square error and the same computational scaling as TDDFT under the core-valence separation approximation, we anticipate that this approach will be of great utility in XAS calculations of large systems.
Collapse
Affiliation(s)
- Kevin Carter-Fenk
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California94720, United States
| | - Leonardo A Cunha
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California94720, United States
| | - Juan E Arias-Martinez
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| |
Collapse
|
8
|
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
|
9
|
Straatsma TP, Broer R, Sánchez-Mansilla A, Sousa C, de Graaf C. GronOR: Scalable and Accelerated Nonorthogonal Configuration Interaction for Molecular Fragment Wave Functions. J Chem Theory Comput 2022; 18:3549-3565. [PMID: 35640094 DOI: 10.1021/acs.jctc.2c00266] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
GronOR is a program package for nonorthogonal configuration interaction calculations. Electronic wave functions are constructed in terms of antisymmetrized products of multiconfiguration molecular fragment wave functions. The computational complexity of the nonorthogonal methodologies implemented in GronOR applied to large molecular assemblies requires a design that takes full advantage of massively parallel supercomputer architectures and accelerator technologies. This work describes the implementation strategy and resulting performance characteristics. In addition to parallelization and acceleration, the software development strategy includes aspects of fault resiliency and heterogeneous computing. The program was designed for large-scale supercomputers but also runs effectively on small clusters and workstations for small molecular systems. GronOR is available as open source to the scientific community.
Collapse
Affiliation(s)
- T P Straatsma
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6373, United States.,Department of Chemistry and Biochemistry, University of Alabama, Tuscaloosa, Alabama 35487-0336, United States
| | - R Broer
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, 9747 AG Groningen, The Netherlands
| | - A Sánchez-Mansilla
- Department of Physical and Inorganic Chemistry, Universitat Rovira i Virgili, C. Marcel·lí Domingo 1, 43007 Tarragona, Spain
| | - C Sousa
- Department of Physical Chemistry and Institut de Química Teòrica i Computacional, Universitat de Barcelona, 08028 Barcelona, Spain
| | - C de Graaf
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, 9747 AG Groningen, The Netherlands.,Department of Physical and Inorganic Chemistry, Universitat Rovira i Virgili, C. Marcel·lí Domingo 1, 43007 Tarragona, Spain.,ICREA, Pg. Lluís Companys 23, 08010 Barcelona, Spain
| |
Collapse
|
10
|
Hait D, Oosterbaan KJ, Carter-Fenk K, Head-Gordon M. Computing x-ray absorption spectra from linear-response particles atop optimized holes. J Chem Phys 2022; 156:201104. [PMID: 35649868 DOI: 10.1063/5.0092987] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Abstract
State specific orbital optimized density functional theory (OO-DFT) methods, such as restricted open-shell Kohn-Sham (ROKS), can attain semiquantitative accuracy for predicting x-ray absorption spectra of closed-shell molecules. OO-DFT methods, however, require that each state be individually optimized. In this Communication, we present an approach to generate an approximate core-excited state density for use with the ROKS energy ansatz, which is capable of giving reasonable accuracy without requiring state-specific optimization. This is achieved by fully optimizing the core-hole through the core-ionized state, followed by the use of electron-addition configuration interaction singles to obtain the particle level. This hybrid approach can be viewed as a DFT generalization of the static-exchange (STEX) method and can attain ∼0.6 eV rms error for the K-edges of C-F through the use of local functionals, such as PBE and OLYP. This ROKS(STEX) approach can also be used to identify important transitions for full OO ROKS treatment and can thus help reduce the computational cost of obtaining OO-DFT quality spectra. ROKS(STEX), therefore, appears to be a useful technique for the efficient prediction of x-ray absorption spectra.
Collapse
Affiliation(s)
- Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Katherine J Oosterbaan
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Kevin Carter-Fenk
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
11
|
Zheng X, Zhang C, Jin Z, Southworth SH, Cheng L. Benchmark relativistic delta-coupled-cluster calculations of K-edge core-ionization energies of third-row elements. Phys Chem Chem Phys 2022; 24:13587-13596. [PMID: 35616685 DOI: 10.1039/d2cp00993e] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
A benchmark computational study of K-edge core-ionization energies of third-row elements using relativistic delta-coupled-cluster (ΔCC) methods and a revised core-valence separation (CVS) scheme is reported. High-level relativistic (HLR) corrections beyond the spin-free exact two-component theory in its one-electron variant (SFX2C-1e), including the contributions from two-electron picture-change effects, spin-orbit coupling, the Breit term, and quantum electrodynamics effects, have been taken into account and demonstrated to play an important role. Relativistic ΔCC calculations are shown to provide accurate results for core-ionization energies of third-row elements. The SFX2C-1e-CVS-ΔCC results augmented with HLR corrections show a maximum deviation of less than 0.5 eV with respect to experimental values.
Collapse
Affiliation(s)
- Xuechen Zheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, MD 21218, USA.
| | - Chaoqun Zhang
- Department of Chemistry, The Johns Hopkins University, Baltimore, MD 21218, USA.
| | - Zheqi Jin
- Department of Chemistry, University College London, London, WC1E 6BT, UK
| | - Stephen H Southworth
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, MD 21218, USA.
| |
Collapse
|
12
|
Sánchez-Mansilla A, Sousa C, Kathir RK, Broer R, Straatsma TP, de Graaf C. On the role of dynamic electron correlation in non-orthogonal configuration interaction with fragments. Phys Chem Chem Phys 2022; 24:11931-11944. [PMID: 35521680 DOI: 10.1039/d2cp00772j] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Two different approaches have been implemented to include the effect of dynamic electron correlation in the Non-Orthogonal Configuration Interaction for Fragments (NOCI-F) method. The first is based on shifting the diagonal matrix elements of the NOCI matrix, while the second incorporates the dynamic correlation explicitly in the fragment wave functions used to construct the many-electron basis functions of the NOCI. The two approaches are illustrated for the calculation of the electronic coupling relevant in singlet fission and the coupling of spin moments in organic radicals. Comparison of the calculated diabatic couplings, the NOCI energies and wave functions shows that dynamic electron correlation is not only efficiently but also effectively incorporated by the shifting approach and can largely affect the coupling between electronic states. Also, it brings the NOCI coupling of the spin moments in close agreement with benchmark calculations.
Collapse
Affiliation(s)
- A Sánchez-Mansilla
- Departament de Química Física i Inorgànica, Universitat Rovira i Virgili, Tarragona, Spain
| | - C Sousa
- Departament de Ciència de Materials i Química Física and Institut de Química Teòrica i Computacional, Universitat de Barcelona, Spain.
| | - R K Kathir
- Zernike Institute of Advanced Materials, University of Groningen, The Netherlands
| | - R Broer
- Zernike Institute of Advanced Materials, University of Groningen, The Netherlands
| | - T P Straatsma
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, TN 37831-6373, USA.,Department of Chemistry and Biochemistry, University of Alabama, Tuscaloosa, AL 35487-0336, USA
| | - C de Graaf
- Departament de Química Física i Inorgànica, Universitat Rovira i Virgili, Tarragona, Spain.,Zernike Institute of Advanced Materials, University of Groningen, The Netherlands.,ICREA, Pg. Lluís Companys 23, Barcelona, Spain.
| |
Collapse
|
13
|
Carter-Fenk K, Head-Gordon M. On the choice of reference orbitals for linear-response calculations of solution-phase K-edge X-ray absorption spectra. Phys Chem Chem Phys 2022; 24:26170-26179. [DOI: 10.1039/d2cp04077h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
X-ray absorption spectra of liquids calculated with linear-response theories like TDDFT and CIS are dramatically improved with core-ion reference orbitals.
Collapse
Affiliation(s)
- Kevin Carter-Fenk
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| |
Collapse
|
14
|
Huang M, Li C, Evangelista FA. Theoretical Calculation of Core-Excited States along Dissociative Pathways beyond Second-Order Perturbation Theory. J Chem Theory Comput 2021; 18:219-233. [PMID: 34964628 DOI: 10.1021/acs.jctc.1c00884] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We extend the multireference driven similarity renormalization (MR-DSRG) method to compute core-excited states by combining it with a GASSCF treatment of orbital relaxation and static electron correlation effects. We consider MR-DSRG treatments of dynamical correlation truncated at the level of perturbation theory (DSRG-MRPT2/3) and iterative linearized approximations with one- and two-body operators [MR-LDSRG(2)] in combination with a spin-free exact-two-component (X2C) one-electron treatment of scalar relativistic effects. This approach is calibrated and tested on a series of 16 core-excited states of five closed- and open-shell diatomic molecules containing first-row elements (C, N, and O). All GASSCF-MR-DSRG theories show excellent agreement with experimental adiabatic transitions energies, with mean absolute errors ranging between 0.17 and 0.35 eV, even for the challenging partially doubly excited states of the N2+ molecule. The vibrational structure of all these transitions, obtained from using a full potential energy scan, shows a mean absolute error as low as 25 meV for DSRG-MRPT2 and 12/13 meV for DSRG-MRPT3 and MR-LDSRG(2). We generally find that a treatment of dynamical correlation that goes beyond the second-order level in perturbation theory improves the accuracy of the potential energy surface, especially in the bond-dissociation region.
Collapse
Affiliation(s)
- Meng Huang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chenyang Li
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States.,Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
15
|
David G, Irons TJP, Fouda AEA, Furness JW, Teale AM. Self-Consistent Field Methods for Excited States in Strong Magnetic Fields: a Comparison between Energy- and Variance-Based Approaches. J Chem Theory Comput 2021; 17:5492-5508. [PMID: 34517708 DOI: 10.1021/acs.jctc.1c00236] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Self-consistent field methods for excited states offer an attractive low-cost route to study not only excitation energies but also properties of excited states. Here, we present the generalization of two self-consistent field methods, the maximum overlap method (MOM) and the σ-SCF method, to calculate excited states in strong magnetic fields and investigate their stability and accuracy in this context. These methods use different strategies to overcome the well-known variational collapse of energy-based optimizations to the lowest solution of a given symmetry. The MOM tackles this problem in the definition of the orbital occupations to constrain the self-consistent field procedure to converge on excited states, while the σ-SCF method is based on the minimization of the variance instead of the energy. To overcome the high computational cost of the variance minimization, we present a new implementation of the σ-SCF method with the resolution of identity approximation, allowing the use of large basis sets, which is an important requirement for calculations in strong magnetic fields. The accuracy of these methods is assessed by comparison with the benchmark literature data for He, H2, and CH+. The results reveal severe limitations of the variance-based scheme, which become more acute in large basis sets. In particular, many states are not accessible using variance optimization. Detailed analysis shows that this is a general feature of variance optimization approaches due to the masking of local minima in the optimization. In contrast, the MOM shows promising performance for computing excited states under these conditions, yielding results consistent with available benchmark data for a diverse range of electronic states.
Collapse
Affiliation(s)
- Grégoire David
- School of Chemistry, University of Nottingham, Nottingham NG7 2RD, U.K
| | - Tom J P Irons
- School of Chemistry, University of Nottingham, Nottingham NG7 2RD, U.K
| | - Adam E A Fouda
- School of Chemistry, University of Nottingham, Nottingham NG7 2RD, U.K.,Chemical Sciences and Engineering Division, Argonne National Laboratory, 9700 S Cass Avenue, Lemont, Illinois 60439, United States
| | - James W Furness
- Department of Physics and Engineering Physics, Tulane University, New Orleans, Louisiana 70118, United States
| | - Andrew M Teale
- School of Chemistry, University of Nottingham, Nottingham NG7 2RD, U.K.,Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, Oslo N-0315, Norway
| |
Collapse
|
16
|
Epifanovsky E, Gilbert ATB, Feng X, Lee J, Mao Y, Mardirossian N, Pokhilko P, White AF, Coons MP, Dempwolff AL, Gan Z, Hait D, Horn PR, Jacobson LD, Kaliman I, Kussmann J, Lange AW, Lao KU, Levine DS, Liu J, McKenzie SC, Morrison AF, Nanda KD, Plasser F, Rehn DR, Vidal ML, You ZQ, Zhu Y, Alam B, Albrecht BJ, Aldossary A, Alguire E, Andersen JH, Athavale V, Barton D, Begam K, Behn A, Bellonzi N, Bernard YA, Berquist EJ, Burton HGA, Carreras A, Carter-Fenk K, Chakraborty R, Chien AD, Closser KD, Cofer-Shabica V, Dasgupta S, de Wergifosse M, Deng J, Diedenhofen M, Do H, Ehlert S, Fang PT, Fatehi S, Feng Q, Friedhoff T, Gayvert J, Ge Q, Gidofalvi G, Goldey M, Gomes J, González-Espinoza CE, Gulania S, Gunina AO, Hanson-Heine MWD, Harbach PHP, Hauser A, Herbst MF, Hernández Vera M, Hodecker M, Holden ZC, Houck S, Huang X, Hui K, Huynh BC, Ivanov M, Jász Á, Ji H, Jiang H, Kaduk B, Kähler S, Khistyaev K, Kim J, Kis G, Klunzinger P, Koczor-Benda Z, Koh JH, Kosenkov D, Koulias L, Kowalczyk T, Krauter CM, Kue K, Kunitsa A, Kus T, Ladjánszki I, Landau A, Lawler KV, Lefrancois D, Lehtola S, Li RR, Li YP, Liang J, Liebenthal M, Lin HH, Lin YS, Liu F, Liu KY, Loipersberger M, Luenser A, Manjanath A, Manohar P, Mansoor E, Manzer SF, Mao SP, Marenich AV, Markovich T, Mason S, Maurer SA, McLaughlin PF, Menger MFSJ, Mewes JM, Mewes SA, Morgante P, Mullinax JW, Oosterbaan KJ, Paran G, Paul AC, Paul SK, Pavošević F, Pei Z, Prager S, Proynov EI, Rák Á, Ramos-Cordoba E, Rana B, Rask AE, Rettig A, Richard RM, Rob F, Rossomme E, Scheele T, Scheurer M, Schneider M, Sergueev N, Sharada SM, Skomorowski W, Small DW, Stein CJ, Su YC, Sundstrom EJ, Tao Z, Thirman J, Tornai GJ, Tsuchimochi T, Tubman NM, Veccham SP, Vydrov O, Wenzel J, Witte J, Yamada A, Yao K, Yeganeh S, Yost SR, Zech A, Zhang IY, Zhang X, Zhang Y, Zuev D, Aspuru-Guzik A, Bell AT, Besley NA, Bravaya KB, Brooks BR, Casanova D, Chai JD, Coriani S, Cramer CJ, Cserey G, DePrince AE, DiStasio RA, Dreuw A, Dunietz BD, Furlani TR, Goddard WA, Hammes-Schiffer S, Head-Gordon T, Hehre WJ, Hsu CP, Jagau TC, Jung Y, Klamt A, Kong J, Lambrecht DS, Liang W, Mayhall NJ, McCurdy CW, Neaton JB, Ochsenfeld C, Parkhill JA, Peverati R, Rassolov VA, Shao Y, Slipchenko LV, Stauch T, Steele RP, Subotnik JE, Thom AJW, Tkatchenko A, Truhlar DG, Van Voorhis T, Wesolowski TA, Whaley KB, Woodcock HL, Zimmerman PM, Faraji S, Gill PMW, Head-Gordon M, Herbert JM, Krylov AI. Software for the frontiers of quantum chemistry: An overview of developments in the Q-Chem 5 package. J Chem Phys 2021; 155:084801. [PMID: 34470363 PMCID: PMC9984241 DOI: 10.1063/5.0055522] [Citation(s) in RCA: 412] [Impact Index Per Article: 137.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023] Open
Abstract
This article summarizes technical advances contained in the fifth major release of the Q-Chem quantum chemistry program package, covering developments since 2015. A comprehensive library of exchange-correlation functionals, along with a suite of correlated many-body methods, continues to be a hallmark of the Q-Chem software. The many-body methods include novel variants of both coupled-cluster and configuration-interaction approaches along with methods based on the algebraic diagrammatic construction and variational reduced density-matrix methods. Methods highlighted in Q-Chem 5 include a suite of tools for modeling core-level spectroscopy, methods for describing metastable resonances, methods for computing vibronic spectra, the nuclear-electronic orbital method, and several different energy decomposition analysis techniques. High-performance capabilities including multithreaded parallelism and support for calculations on graphics processing units are described. Q-Chem boasts a community of well over 100 active academic developers, and the continuing evolution of the software is supported by an "open teamware" model and an increasingly modular design.
Collapse
Affiliation(s)
- Evgeny Epifanovsky
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | | | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yuezhi Mao
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Pavel Pokhilko
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Alec F. White
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Marc P. Coons
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Adrian L. Dempwolff
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Zhengting Gan
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Diptarka Hait
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Paul R. Horn
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Leif D. Jacobson
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | | - Jörg Kussmann
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Adrian W. Lange
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Daniel S. Levine
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Simon C. McKenzie
- Research School of Chemistry, Australian National University, Canberra, Australia
| | | | - Kaushik D. Nanda
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Dirk R. Rehn
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Marta L. Vidal
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | | | - Ying Zhu
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Bushra Alam
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Benjamin J. Albrecht
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | - Ethan Alguire
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Josefine H. Andersen
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | - Vishikh Athavale
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Dennis Barton
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Khadiza Begam
- Department of Physics, Kent State University, Kent, Ohio 44242, USA
| | - Andrew Behn
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Nicole Bellonzi
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Yves A. Bernard
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Hugh G. A. Burton
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Abel Carreras
- Donostia International Physics Center, 20080 Donostia, Euskadi, Spain
| | - Kevin Carter-Fenk
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | | - Alan D. Chien
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | | | - Vale Cofer-Shabica
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Saswata Dasgupta
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Marc de Wergifosse
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Jia Deng
- Research School of Chemistry, Australian National University, Canberra, Australia
| | | | - Hainam Do
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Beringstr. 4, 53115 Bonn, Germany
| | - Po-Tung Fang
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | | | - Qingguo Feng
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Triet Friedhoff
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - James Gayvert
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Qinghui Ge
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Gergely Gidofalvi
- Department of Chemistry and Biochemistry, Gonzaga University, Spokane, Washington 99258, USA
| | - Matthew Goldey
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Joe Gomes
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Sahil Gulania
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Anastasia O. Gunina
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Phillip H. P. Harbach
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Andreas Hauser
- Institute of Experimental Physics, Graz University of Technology, Graz, Austria
| | | | - Mario Hernández Vera
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Manuel Hodecker
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Zachary C. Holden
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Shannon Houck
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Xunkun Huang
- Department of Chemistry, Xiamen University, Xiamen 361005, China
| | - Kerwin Hui
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Bang C. Huynh
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Maxim Ivanov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Ádám Jász
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Hyunjun Ji
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Hanjie Jiang
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Benjamin Kaduk
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Sven Kähler
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Kirill Khistyaev
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Jaehoon Kim
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Gergely Kis
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | | | - Zsuzsanna Koczor-Benda
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Joong Hoon Koh
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Dimitri Kosenkov
- Department of Chemistry, Purdue University, West Lafayette, Indiana 47907, USA
| | - Laura Koulias
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | | | - Caroline M. Krauter
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Karl Kue
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - Alexander Kunitsa
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Thomas Kus
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | | | - Arie Landau
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Keith V. Lawler
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Daniel Lefrancois
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | | | - Run R. Li
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Yi-Pei Li
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Jiashu Liang
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Marcus Liebenthal
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Hung-Hsuan Lin
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - You-Sheng Lin
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Fenglai Liu
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | | | - Arne Luenser
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Aaditya Manjanath
- Institute of Chemistry, Academia Sinica, 128, Academia Road Section 2, Nangang District, Taipei 11529, Taiwan
| | - Prashant Manohar
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Erum Mansoor
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Sam F. Manzer
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Shan-Ping Mao
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | | | - Thomas Markovich
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Stephen Mason
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Simon A. Maurer
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - Peter F. McLaughlin
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | | | - Jan-Michael Mewes
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Stefanie A. Mewes
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Pierpaolo Morgante
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - J. Wayne Mullinax
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | | | | | - Alexander C. Paul
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Suranjan K. Paul
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Fabijan Pavošević
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | - Zheng Pei
- School of Electrical and Computer Engineering, University of Oklahoma, Norman, Oklahoma 73019, USA
| | - Stefan Prager
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Emil I. Proynov
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Ádám Rák
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Eloy Ramos-Cordoba
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Bhaskar Rana
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Alan E. Rask
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Adam Rettig
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Ryan M. Richard
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Fazle Rob
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Elliot Rossomme
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Tarek Scheele
- Institute for Physical and Theoretical Chemistry, University of Bremen, Bremen, Germany
| | - Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Matthias Schneider
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Nickolai Sergueev
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Shaama M. Sharada
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Wojciech Skomorowski
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - David W. Small
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Christopher J. Stein
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Yu-Chuan Su
- Department of Physics, National Taiwan University, Taipei 10617, Taiwan
| | - Eric J. Sundstrom
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Zhen Tao
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, USA
| | - Jonathan Thirman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Gábor J. Tornai
- Stream Novation Ltd., Práter utca 50/a, H-1083 Budapest, Hungary
| | - Takashi Tsuchimochi
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Norm M. Tubman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | - Oleg Vydrov
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Jan Wenzel
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Jon Witte
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Atsushi Yamada
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Kun Yao
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Sina Yeganeh
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Shane R. Yost
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Alexander Zech
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Geneva 4, Switzerland
| | - Igor Ying Zhang
- Department of Chemistry, Fudan University, Shanghai 200433, China
| | - Xing Zhang
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Yu Zhang
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Dmitry Zuev
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Alán Aspuru-Guzik
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, USA
| | - Alexis T. Bell
- Department of Chemical Engineering, University of California, Berkeley, California 94720, USA
| | - Nicholas A. Besley
- School of Chemistry, University of Nottingham, Nottingham, United Kingdom
| | - Ksenia B. Bravaya
- Department of Chemistry, Boston University, Boston, Massachusetts 02215, USA
| | - Bernard R. Brooks
- Laboratory of Computational Biophysics, National Institute of Health, Bethesda, Maryland 20892, USA
| | - David Casanova
- Donostia International Physics Center, 20080 Donostia, Euskadi, Spain
| | | | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kgs Lyngby, Denmark
| | | | | | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, USA
| | - Robert A. DiStasio
- Department of Chemistry and Chemical Biology, Cornell University, Ithaca, New York 14853, USA
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Barry D. Dunietz
- Department of Chemistry and Biochemistry, Kent State University, Kent, Ohio 44240, USA
| | - Thomas R. Furlani
- Department of Chemistry, University at Buffalo, State University of New York, Buffalo, New York 14260, USA
| | - William A. Goddard
- Materials and Process Simulation Center, California Institute of Technology, Pasadena, California 91125, USA
| | | | - Teresa Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | | | | | | | - Yousung Jung
- Graduate School of Energy, Environment, Water and Sustainability (EEWS), Korea Advanced Institute of Science and Technology (KAIST), Daejeon 34141, Republic of Korea
| | - Andreas Klamt
- COSMOlogic GmbH & Co. KG, Imbacher Weg 46, D-51379 Leverkusen, Germany
| | - Jing Kong
- Q-Chem, Inc., 6601 Owens Drive, Suite 105, Pleasanton, California 94588, USA
| | - Daniel S. Lambrecht
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | | | | | - C. William McCurdy
- Department of Chemistry, University of California, Davis, California 95616, USA
| | - Jeffrey B. Neaton
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Christian Ochsenfeld
- Department of Chemistry, Ludwig Maximilian University, Butenandtstr. 7, D-81377 München, Germany
| | - John A. Parkhill
- Department of Chemistry and Biochemistry, University of Notre Dame, Notre Dame, Indiana 46556, USA
| | - Roberto Peverati
- Department of Chemistry, Florida Institute of Technology, Melbourne, Florida 32901, USA
| | - Vitaly A. Rassolov
- Department of Chemistry and Biochemistry, University of South Carolina, Columbia, South Carolina 29208, USA
| | | | | | | | - Ryan P. Steele
- Department of Chemistry and Henry Eyring Center for Theoretical Chemistry, University of Utah, Salt Lake City, Utah 84112, USA
| | - Joseph E. Subotnik
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, USA
| | - Alex J. W. Thom
- Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511 Luxembourg, Luxembourg
| | - Donald G. Truhlar
- Department of Chemistry, University of Minnesota, Minneapolis, Minnesota 55455, USA
| | - Troy Van Voorhis
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Tomasz A. Wesolowski
- Department of Physical Chemistry, University of Geneva, 30, Quai Ernest-Ansermet, CH-1211 Geneva 4, Switzerland
| | - K. Birgitta Whaley
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - H. Lee Woodcock
- Department of Chemistry, University of South Florida, Tampa, Florida 33620, USA
| | - Paul M. Zimmerman
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Shirin Faraji
- Zernike Institute for Advanced Materials, University of Groningen, 9774AG Groningen, The Netherlands
| | | | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA,Author to whom correspondence should be addressed:
| |
Collapse
|
17
|
Zhao R, Grofe A, Wang Z, Bao P, Chen X, Liu W, Gao J. Dynamic-then-Static Approach for Core Excitations of Open-Shell Molecules. J Phys Chem Lett 2021; 12:7409-7417. [PMID: 34328742 DOI: 10.1021/acs.jpclett.1c02039] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
Delta self-consistent-field methods are widely used in studies of electronically excited states. However, the nonaufbau determinants are generally spin-contaminated. Here, we describe a general approach for spin-coupling interactions of open-shell molecules, making use of multistate density functional theory (MSDFT). In particular, the effective exchange integrals that determine spin coupling are obtained by enforcing the multiplet degeneracy of the S+1 state in the MS = S manifold. Consequently, they are consistent with the energy of the high-spin state that is adequately treated by Kohn-Sham density functional theory (DFT) and, thereby, free of double counting of correlation. The method was applied to core excitations of open-shell molecules and compared with those by spin-adapted time-dependent DFT. An excellent agreement with experiment was found employing the BLYP functional and aug-cc-pCVQZ basis set. Overall, MSDFT provides an effective combination of the strengths of DFT and wave function theory to achieve efficiency and accuracy.
Collapse
Affiliation(s)
- Ruoqi Zhao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, Guangdong, China
- Institute of Theoretical Chemistry, Jilin University Changchun 130023, Jilin, China
| | - Adam Grofe
- Institute of Theoretical Chemistry, Jilin University Changchun 130023, Jilin, China
| | - Zikuan Wang
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, Shandong, China
| | - Peng Bao
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Xin Chen
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, Guangdong, China
| | - Wenjian Liu
- Qingdao Institute for Theoretical and Computational Sciences, Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao 266237, Shandong, China
| | - Jiali Gao
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, Guangdong, China
- Department of Chemistry and Supercomputing Institute, University of Minnesota, Minneapolis 55455, Minnesota, United States
| |
Collapse
|
18
|
Fouda AEA, Ho PJ. Site-specific generation of excited state wavepackets with high-intensity attosecond x rays. J Chem Phys 2021; 154:224111. [PMID: 34241215 DOI: 10.1063/5.0050891] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
High-intensity attosecond x rays can produce coherent superpositions of valence-excited states through two-photon Raman transitions. The broad-bandwidth, high-field nature of the pulses results in a multitude of accessible excited states. Multiconfigurational quantum chemistry with the time-dependent Schrödinger equation is used to examine population transfer dynamics in stimulated x-ray Raman scattering of the nitric oxide oxygen and nitrogen K-edges. Two pulse schemes initiate wavepackets of different characters and demonstrate how chemical differences between core-excitation pathways affect the dynamics. The population transfer to valence-excited states is found to be sensitive to the electronic structure and pulse conditions, highlighting complexities attributed to the Rabi frequency. The orthogonally polarized two-color-pulse setup has increased selectivity while facilitating longer, less intense pulses than the one-pulse setup. Population transfer in the 1s → Rydberg region is more effective but less selective at the nitrogen K-edge; the selectivity is reduced by double core-excited states. Result interpretation is aided by resonant inelastic x-ray scattering maps.
Collapse
Affiliation(s)
- Adam E A Fouda
- Chemical Sciences and Engineering Division, Argonne National Laboratory, 9700 S Cass Avenue, Lemont, Illinois 60439, USA
| | - Phay J Ho
- Chemical Sciences and Engineering Division, Argonne National Laboratory, 9700 S Cass Avenue, Lemont, Illinois 60439, USA
| |
Collapse
|
19
|
Besley NA. Density Functional Theory Calculations of Core-Electron Binding Energies at the K-Edge of Heavier Elements. J Chem Theory Comput 2021; 17:3644-3651. [PMID: 34042434 DOI: 10.1021/acs.jctc.1c00171] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
The capability to determine core-electron binding energies (CEBEs) is vital in the analysis of X-ray photoelectron spectroscopy, and the continued development of light sources has made inner shell spectroscopy of heavier elements increasingly accessible. Density functional theory is widely used to determine CEBEs of lighter elements (boron-fluorine). It is shown that good performance of exchange-correlation functionals for these elements does not necessarily translate to the calculation of CEBEs for the heavier elements from the next row of the periodic table, and in general, larger errors are observed. Two strategies are explored that improve the accuracy of the calculated CEBEs. The first is to apply element and functional dependent energy corrections, and the second is a reparametrization of a short-range corrected functional. This functional is able to reproduce experimental phosphorus and sulfur K-edge CEBEs with an average error of 0.15 eV demonstrating the importance of reducing the self-interaction error associated with the core electrons and represents progress toward a density functional theory calculation that performs equally well for ionization at the K-edge of all elements.
Collapse
Affiliation(s)
- Nicholas A Besley
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, U.K
| |
Collapse
|
20
|
Burton HGA. Generalized nonorthogonal matrix elements: Unifying Wick's theorem and the Slater-Condon rules. J Chem Phys 2021; 154:144109. [PMID: 33858143 DOI: 10.1063/5.0045442] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Matrix elements between nonorthogonal Slater determinants represent an essential component of many emerging electronic structure methods. However, evaluating nonorthogonal matrix elements is conceptually and computationally harder than their orthogonal counterparts. While several different approaches have been developed, these are predominantly derived from the first-quantized generalized Slater-Condon rules and usually require biorthogonal occupied orbitals to be computed for each matrix element. For coupling terms between nonorthogonal excited configurations, a second-quantized approach such as the nonorthogonal Wick's theorem is more desirable, but this fails when the two reference determinants have a zero many-body overlap. In this contribution, we derive an entirely generalized extension to the nonorthogonal Wick's theorem that is applicable to all pairs of determinants with nonorthogonal orbitals. Our approach creates a universal methodology for evaluating any nonorthogonal matrix element and allows Wick's theorem and the generalized Slater-Condon rules to be unified for the first time. Furthermore, we present a simple well-defined protocol for deriving arbitrary coupling terms between nonorthogonal excited configurations. In the case of overlap and one-body operators, this protocol recovers efficient formulas with reduced scaling, promising significant computational acceleration for methods that rely on such terms.
Collapse
Affiliation(s)
- Hugh G A Burton
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford OX1 3QZ, United Kingdom
| |
Collapse
|
21
|
Bilalbegović G, Maksimović A, Valencic LA, Lehtola S. Sulfur Molecules in Space by X-rays: A Computational Study. ACS EARTH & SPACE CHEMISTRY 2021; 5:436-448. [PMID: 33842801 PMCID: PMC8028330 DOI: 10.1021/acsearthspacechem.0c00238] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/29/2020] [Revised: 02/12/2021] [Accepted: 02/15/2021] [Indexed: 06/12/2023]
Abstract
X-ray astronomy lacks high resolution spectra of interstellar dust analogues and molecules, severely hampering interstellar medium studies based on upcoming X-ray missions. Various theoretical approaches may be used to address this problem, but they must first be shown to reproduce reliable spectra compared to the experiment. In this work, we calculate the sulfur K edge X-ray absorption spectra of H2S, SO2, and OCS, whose spectra are already known from X-ray experiments and predict the X-ray spectrum of CS, which as far as we are aware has not been measured, thereby hampering its detection by X-ray telescopes. We chose these four molecules as the astrochemistry of sulfur is an unsolved problem and as the four molecules are already known to exist in space. We consider three types of methods for modeling the X-ray spectra: more accurate calculations with the algebraic-diagrammatic construction (ADC) and the CC2, CCSD, and CC3 coupled cluster (CC) approaches as well as more affordable ones with transition potential density functional theory (TP-DFT). A comparison of our computational results to previously reported experimental spectra shows that the core-valence separation (CVS) approaches CVS-ADC(2)-x and CVS-CC3 generally yield a good qualitative level of agreement with the experiment, suggesting that they can be used for interpreting measured spectra, while the TP-DFT method is not reliable for these molecules. However, quantitative agreement with the experiment is still outside the reach of the computational methods studied in this work.
Collapse
Affiliation(s)
- Goranka Bilalbegović
- Department
of Physics, Faculty of Science, University
of Zagreb, Bijenička
cesta 32, 10000 Zagreb, Croatia
| | - Aleksandar Maksimović
- Center
of Excellence for Advanced Materials and Sensing Devices, Rudjer Bošković Institute, Bijenička cesta 54, 10000 Zagreb, Croatia
| | - Lynne A. Valencic
- NASA
Goddard Space Flight Center, Greenbelt, 20771 Maryland, United States
- Department
of Physics & Astronomy, The Johns Hopkins
University, 366 Bloomberg
Center, 3400 N. Charles Street, Baltimore, 21218 Maryland, United States
| | - Susi Lehtola
- Department
of Chemistry, University of Helsinki, P.O. Box 55, A. I. Virtasen aukio
1, FI-00014 Helsinki, Finland
| |
Collapse
|
22
|
Besley NA. Modeling of the spectroscopy of core electrons with density functional theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2021. [DOI: 10.1002/wcms.1527] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Affiliation(s)
- Nicholas A. Besley
- School of Chemistry, University of Nottingham University Park Nottingham UK
| |
Collapse
|
23
|
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
|
24
|
Garner SM, Neuscamman E. A variational Monte Carlo approach for core excitations. J Chem Phys 2020; 153:144108. [DOI: 10.1063/5.0020310] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Affiliation(s)
- Scott M. Garner
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Eric Neuscamman
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
25
|
Hait D, Haugen EA, Yang Z, Oosterbaan KJ, Leone SR, Head-Gordon M. Accurate prediction of core-level spectra of radicals at density functional theory cost via square gradient minimization and recoupling of mixed configurations. J Chem Phys 2020; 153:134108. [DOI: 10.1063/5.0018833] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Diptarka Hait
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Eric A. Haugen
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Zheyue Yang
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Katherine J. Oosterbaan
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Stephen R. Leone
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Department of Physics, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
26
|
Fouda AEA, Seitz LC, Hauschild D, Blum M, Yang W, Heske C, Weinhardt L, Besley NA. Observation of Double Excitations in the Resonant Inelastic X-ray Scattering of Nitric Oxide. J Phys Chem Lett 2020; 11:7476-7482. [PMID: 32787301 DOI: 10.1021/acs.jpclett.0c01981] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
The nitrogen K-edge resonant inelastic X-ray scattering (RIXS) map of nitric oxide (NO) has been measured and simulated to provide a detailed analysis of the observed features. High-resolution experimental RIXS maps were collected using an in situ gas flow cell and a high-transmission soft X-ray spectrometer. Accurate descriptions of the ground, excited, and core-excited states are based upon restricted active space self-consistent-field calculations using second order multiconfigurational perturbation theory. The nitrogen K-edge RIXS map of NO shows a range of features that can be assigned to intermediate states arising from 1s → π* and 1s → Rydberg excitations; additional bands are attributed to doubly excited intermediate states comprising 1s → π* and π → π* excitations. These results provide a detailed picture of RIXS for an open-shell molecule and an extensive description of the core-excited electronic structure of NO, an important molecule in many chemical and biological processes.
Collapse
Affiliation(s)
- Adam E A Fouda
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, U.K
| | - Linsey C Seitz
- Institute for Photon Science and Synchrotron Radiation (IPS), Karlsruhe Institute of Technology (KIT), Hermann-v.-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- Department of Chemical and Biological Engineering, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, United States
| | - Dirk Hauschild
- Institute for Photon Science and Synchrotron Radiation (IPS), Karlsruhe Institute of Technology (KIT), Hermann-v.-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- Institute for Chemical Technology and Polymer Chemistry (ITCP), Karlsruhe Institute of Technology (KIT), Engesserstrasse 18/20, 76128 Karlsruhe, Germany
- Department of Chemistry and Biochemistry, University of Nevada, Las Vegas (UNLV), 4505 Maryland Parkway, Las Vegas, Nevada 89154-4003, United States
| | - Monika Blum
- Department of Chemistry and Biochemistry, University of Nevada, Las Vegas (UNLV), 4505 Maryland Parkway, Las Vegas, Nevada 89154-4003, United States
- Advanced Light Source, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, California 94720, United States
| | - Wanli Yang
- Advanced Light Source, Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, California 94720, United States
| | - Clemens Heske
- Institute for Photon Science and Synchrotron Radiation (IPS), Karlsruhe Institute of Technology (KIT), Hermann-v.-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- Institute for Chemical Technology and Polymer Chemistry (ITCP), Karlsruhe Institute of Technology (KIT), Engesserstrasse 18/20, 76128 Karlsruhe, Germany
- Department of Chemistry and Biochemistry, University of Nevada, Las Vegas (UNLV), 4505 Maryland Parkway, Las Vegas, Nevada 89154-4003, United States
| | - Lothar Weinhardt
- Institute for Photon Science and Synchrotron Radiation (IPS), Karlsruhe Institute of Technology (KIT), Hermann-v.-Helmholtz-Platz 1, 76344 Eggenstein-Leopoldshafen, Germany
- Institute for Chemical Technology and Polymer Chemistry (ITCP), Karlsruhe Institute of Technology (KIT), Engesserstrasse 18/20, 76128 Karlsruhe, Germany
- Department of Chemistry and Biochemistry, University of Nevada, Las Vegas (UNLV), 4505 Maryland Parkway, Las Vegas, Nevada 89154-4003, United States
| | - Nicholas A Besley
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, U.K
| |
Collapse
|
27
|
Dong X, Mahler AD, Kempfer-Robertson EM, Thompson LM. Global Elucidation of Self-Consistent Field Solution Space Using Basin Hopping. J Chem Theory Comput 2020; 16:5635-5644. [DOI: 10.1021/acs.jctc.0c00488] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xinju Dong
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40205, United States
| | - Andrew D. Mahler
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40205, United States
| | | | - Lee M. Thompson
- Department of Chemistry, University of Louisville, Louisville, Kentucky 40205, United States
| |
Collapse
|
28
|
Burton HGA, Thom AJW. Reaching Full Correlation through Nonorthogonal Configuration Interaction: A Second-Order Perturbative Approach. J Chem Theory Comput 2020; 16:5586-5600. [DOI: 10.1021/acs.jctc.0c00468] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Affiliation(s)
- Hugh G. A. Burton
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| | - Alex J. W. Thom
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, U.K
| |
Collapse
|
29
|
Ehlert C, Klamroth T. PSIXAS: A Psi4 plugin for efficient simulations of X-ray absorption spectra based on the transition-potential and Δ-Kohn-Sham method. J Comput Chem 2020; 41:1781-1789. [PMID: 32394459 DOI: 10.1002/jcc.26219] [Citation(s) in RCA: 18] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2020] [Revised: 04/21/2020] [Accepted: 04/21/2020] [Indexed: 01/25/2023]
Abstract
Near edge X-ray absorption fine structure (NEXAFS) spectra and their pump-probe extension (PP-NEXAFS) offer insights into valence- and core-excited states. We present PSIXAS, a recent implementation for simulating NEXAFS and PP-NEXAFS spectra by means of the transition-potential and the Δ-Kohn-Sham method. The approach is implemented in form of a software plugin for the Psi4 code, which provides access to a wide selection of basis sets as well as density functionals. We briefly outline the theoretical foundation and the key aspects of the plugin. Then, we use the plugin to simulate PP-NEXAFS spectra of thymine, a system already investigated by others and us. It is found that larger, extended basis sets are needed to obtain more accurate absolute resonance positions. We further demonstrate that, in contrast to ordinary NEXAFS simulations, where the choice of the density functional plays a minor role for the shape of the spectrum, for PP-NEXAFS simulations the choice of the density functional is important. Especially hybrid functionals (which could not be used straightforwardly before to simulate PP-NEXAFS spectra) and their amount of "Hartree-Fock like" exact exchange affects relative resonance positions in the spectrum.
Collapse
Affiliation(s)
- Christopher Ehlert
- Heidelberg Institute for Theoretical Studies (HITS gGmbH), Heidelberg, Germany.,Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg, Germany
| | | |
Collapse
|
30
|
Kathir RK, de Graaf C, Broer R, Havenith RWA. Reduced Common Molecular Orbital Basis for Nonorthogonal Configuration Interaction. J Chem Theory Comput 2020; 16:2941-2951. [PMID: 32279493 PMCID: PMC7222100 DOI: 10.1021/acs.jctc.9b01144] [Citation(s) in RCA: 19] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023]
Abstract
Electron and charge transfers are part of many vital processes in nature and technology. Ab initio descriptions of these processes provide useful insights that can be utilized for applications. A combination of the embedded cluster material model and nonorthogonal configuration interaction (NOCI), in which the cluster wave functions are expanded in many-electron basis functions (MEBFs) consisting of spin-adapted, antisymmetrized products of multiconfigurational wave functions of fragments (which are usually molecules) in the cluster, appears to provide a compromise between accuracy and calculation time. Additional advantages of this NOCI-Fragments approach are the chemically convenient interpretation of the wave function in terms of molecular states, and the direct accessibility of electronic coupling between diabatic states to describe energy and electron transfer processes. Bottlenecks in this method are the large number of two-electron integrals that have to be handled for the calculation of an electronic coupling matrix element and the enormous number of matrix elements over determinant pairs that have to be evaluated for the calculation of one matrix element between the MEBFs. We show here how we created a reduced common molecular orbital basis that is utilized to significantly reduce the number of two-electron integrals that need to be handled. The results obtained with this basis do not show any loss of accuracy in relevant quantities like electronic couplings and vertical excitation energies. We also show a significant reduction in computation time without loss in accuracy when matrix elements over determinant pairs with small weights are neglected in the NOCI. These improvements in the methodology render NOCI-Fragments to be also applicable to treat clusters of larger molecular systems with larger atomic basis sets and larger active spaces, as the computation time becomes dependent on the number of occupied orbitals and less dependent on the size of the active space.
Collapse
Affiliation(s)
- R K Kathir
- Theoretical Chemistry, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Coen de Graaf
- Department of Physical and Inorganic Chemistry, Universitat Rovira i Virgili, C. Marcel·lı́ Domingo 1, 43007 Tarragona, Spain.,Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands.,ICREA, Passeig Lluı́s Companys 23, 08010 Barcelona, Spain
| | - Ria Broer
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands
| | - Remco W A Havenith
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands.,Stratingh Institute for Chemistry, University of Groningen, Nijenborgh 4, 9747 AG Groningen, The Netherlands.,Ghent Quantum Chemistry Group, Department of Inorganic and Physical Chemistry, Ghent University, 9000 Gent, Belgium
| |
Collapse
|
31
|
Zheng X, Liu J, Doumy G, Young L, Cheng L. Hetero-site Double Core Ionization Energies with Sub-electronvolt Accuracy from Delta-Coupled-Cluster Calculations. J Phys Chem A 2020; 124:4413-4426. [DOI: 10.1021/acs.jpca.0c00901] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xuechen Zheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Junzi Liu
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Gilles Doumy
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Linda Young
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|
32
|
Straatsma TP, Broer R, Faraji S, Havenith RWA, Suarez LEA, Kathir RK, Wibowo M, de Graaf C. GronOR: Massively parallel and GPU-accelerated non-orthogonal configuration interaction for large molecular systems. J Chem Phys 2020; 152:064111. [PMID: 32061226 DOI: 10.1063/1.5141358] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
GronOR is a program package for non-orthogonal configuration interaction calculations for an electronic wave function built in terms of anti-symmetrized products of multi-configuration molecular fragment wave functions. The two-electron integrals that have to be processed may be expressed in terms of atomic orbitals or in terms of an orbital basis determined from the molecular orbitals of the fragments. The code has been specifically designed for execution on distributed memory massively parallel and Graphics Processing Unit (GPU)-accelerated computer architectures, using an MPI+OpenACC/OpenMP programming approach. The task-based execution model used in the implementation allows for linear scaling with the number of nodes on the largest pre-exascale architectures available, provides hardware fault resiliency, and enables effective execution on systems with distinct central processing unit-only and GPU-accelerated partitions. The code interfaces with existing multi-configuration electronic structure codes that provide optimized molecular fragment orbitals, configuration interaction coefficients, and the required integrals. Algorithm and implementation details, parallel and accelerated performance benchmarks, and an analysis of the sensitivity of the accuracy of results and computational performance to thresholds used in the calculations are presented.
Collapse
Affiliation(s)
- T P Straatsma
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6373, USA
| | - R Broer
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, Groningen, The Netherlands
| | - S Faraji
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, Groningen, The Netherlands
| | - R W A Havenith
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, Groningen, The Netherlands
| | - L E Aguilar Suarez
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, Groningen, The Netherlands
| | - R K Kathir
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, Groningen, The Netherlands
| | - M Wibowo
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, Groningen, The Netherlands
| | - C de Graaf
- Theoretical Chemistry Group, Zernike Institute for Advanced Materials, University of Groningen, Groningen, The Netherlands
| |
Collapse
|
33
|
Hait D, Head-Gordon M. Highly Accurate Prediction of Core Spectra of Molecules at Density Functional Theory Cost: Attaining Sub-electronvolt Error from a Restricted Open-Shell Kohn-Sham Approach. J Phys Chem Lett 2020; 11:775-786. [PMID: 31917579 DOI: 10.1021/acs.jpclett.9b03661] [Citation(s) in RCA: 74] [Impact Index Per Article: 18.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/03/2023]
Abstract
We present the use of the recently developed square gradient minimization (SGM) algorithm for excited-state orbital optimization to obtain spin-pure restricted open-shell Kohn-Sham (ROKS) energies for core excited states of molecules. The SGM algorithm is robust against variational collapse and offers a reliable route to converging orbitals for target excited states at only 2-3 times the cost of ground-state orbital optimization (per iteration). ROKS/SGM with the modern SCAN/ωB97X-V functionals is found to predict the K-edge of C, N, O, and F to a root mean squared error of ∼0.3 eV. ROKS/SGM is equally effective at predicting L-edge spectra of third period elements, provided a perturbative spin-orbit correction is employed. This high accuracy can be contrasted with traditional time-dependent density functional theory (TDDFT), which typically has greater than 10 eV error and requires translation of computed spectra to align with experiment. ROKS is computationally affordable (having the same scaling as ground-state DFT and a slightly larger prefactor) and can be applied to geometry optimizations/ab initio molecular dynamics of core excited states, as well as condensed phase simulations. ROKS can also model doubly excited/ionized states with one broken electron pair, which are beyond the ability of linear response based methods.
Collapse
Affiliation(s)
- Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry , University of California , Berkeley , California 94720 , United States
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry , University of California , Berkeley , California 94720 , United States
- Chemical Sciences Division , Lawrence Berkeley National Laboratory , Berkeley , California 94720 , United States
| |
Collapse
|
34
|
Oosterbaan KJ, White AF, Hait D, Head-Gordon M. Generalized single excitation configuration interaction: an investigation into the impact of the inclusion of non-orthogonality on the calculation of core-excited states. Phys Chem Chem Phys 2020; 22:8182-8192. [DOI: 10.1039/c9cp06592j] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
In this paper, we investigate different non-orthogonal generalizations of the configuration interaction with single substitutions (CIS) method and their impact on the calculation of core-excited states.
Collapse
Affiliation(s)
| | - Alec F. White
- Division of Chemistry and Chemical Engineering
- California Institute of Technology
- Pasadena
- USA
| | - Diptarka Hait
- Department of Chemistry
- University of California
- Berkeley
- USA
- Chemical Sciences Division
| | - Martin Head-Gordon
- Department of Chemistry
- University of California
- Berkeley
- USA
- Chemical Sciences Division
| |
Collapse
|
35
|
Nite J, Jiménez-Hoyos CA. Low-Cost Molecular Excited States from a State-Averaged Resonating Hartree–Fock Approach. J Chem Theory Comput 2019; 15:5343-5351. [DOI: 10.1021/acs.jctc.9b00579] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Jacob Nite
- Department of Chemistry, Wesleyan University, Middletown, Connecticut 06459, United States
| | | |
Collapse
|
36
|
Burton HGA, Thom AJW. General Approach for Multireference Ground and Excited States Using Nonorthogonal Configuration Interaction. J Chem Theory Comput 2019; 15:4851-4861. [DOI: 10.1021/acs.jctc.9b00441] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Affiliation(s)
- Hugh G. A. Burton
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, U.K
| | - Alex J. W. Thom
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge, CB2 1EW, U.K
| |
Collapse
|
37
|
Zheng X, Cheng L. Performance of Delta-Coupled-Cluster Methods for Calculations of Core-Ionization Energies of First-Row Elements. J Chem Theory Comput 2019; 15:4945-4955. [DOI: 10.1021/acs.jctc.9b00568] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Xuechen Zheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Lan Cheng
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| |
Collapse
|