1
|
Lin Z, Liu J, Zhang C, Zheng X, Cheng L. Elucidating Anomalous Intensity Ratios in Chlorine L-Edge X-ray Absorption Spectroscopy: Multiplet Effects and Core Rydberg Transitions. J Phys Chem A 2024; 128:8373-8383. [PMID: 39312206 DOI: 10.1021/acs.jpca.4c04089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/04/2024]
Abstract
A relativistic core-valence-separated equation-of-motion coupled cluster (CVS-EOM-CC) study of chlorine L2,3-edge X-ray absorption near-edge structure (XANES) spectra using CH3Cl and CH2ICl as representative molecules is reported. The nearly identical intensity for the main features in the L2- and L3-edge XANES spectra is attributed to multiplet effects and the overlap between core-valence and core Rydberg transitions. The multiplet effects originating from the interaction between the core hole and the C-Cl σ* orbitals account for around half of the deviation of the L3 and L2 intensity ratio from the 2:1 ratio of the numbers of 2p3/2 and 2p1/2 electrons. The 2p3/2 → 4s core Rydberg transitions are shown to overlap with the 2p1/2 → σ* transitions and contribute to the other half of the intensity anomaly. We demonstrate that triple excitations in CVS-EOM-CC calculations play important roles in accurate simulation of the overlap between the 2p1/2 → σ* and 2p3/2 → 4s transitions.
Collapse
Affiliation(s)
- Zhe Lin
- 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
| | - Chaoqun Zhang
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - 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
|
2
|
Costain TS, Rolston JB, Neville SP, Schuurman MS. A DFT/MRCI Hamiltonian parameterized using only ab initio data. II. Core-excited states. J Chem Phys 2024; 161:114117. [PMID: 39301854 DOI: 10.1063/5.0227385] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2024] [Accepted: 09/05/2024] [Indexed: 09/22/2024] Open
Abstract
A newly parameterized combined density functional theory and multi-reference configuration interaction (DFT/MRCI) Hamiltonian, termed core-valence separation (CVS)-QE12, is defined for the computation of K-shell core-excitation and core-ionization energies. This CVS counterpart to the recently reported QE8 Hamiltonian [Costain et al., J. Chem. Phys, 160, 224106 (2024)] is parameterized by fitting to benchmark quality ab initio data. The definition of the CVS-QE12 and QE8 Hamiltonians differ from previous CVS-DFT/MRCI parameterizations in three primary ways: (i) the replacement of the BHLYP exchange-correlation functional with QTP17 to yield a balanced description of both core and valence excitation energies, (ii) the adoption of a new, three-parameter damping function, and (iii) the introduction of separate scaling of the core-valence and valence-valence Coulombic interactions. Crucially, the parameters of the CVS-QE12 Hamiltonian are obtained via fitting exclusively to highly accurate ab initio vertical core-excitation and ionization energies computed at the CVS-EOM-CCSDT level of theory. The CVS-QE12 Hamiltonian is validated against further benchmark computations and is found to furnish K-edge core vertical excitation and ionization energies exhibiting absolute errors ≤0.5 eV at low computational cost.
Collapse
Affiliation(s)
- Teagan Shane Costain
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
| | - Jibrael B Rolston
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
| | - Simon P Neville
- National Research Council Canada, 100 Sussex Dr., Ottawa, Ontario K1A 0R6, Canada
| | - Michael S Schuurman
- Department of Chemistry and Biomolecular Sciences, University of Ottawa, Ottawa, Ontario K1N 6N5, Canada
- National Research Council Canada, 100 Sussex Dr., Ottawa, Ontario K1A 0R6, Canada
| |
Collapse
|
3
|
Huang M, Evangelista FA. Benchmark Study of Core-Ionization Energies with the Generalized Active Space-Driven Similarity Renormalization Group. J Chem Theory Comput 2024. [PMID: 39271297 PMCID: PMC11428169 DOI: 10.1021/acs.jctc.4c00835] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/15/2024]
Abstract
X-ray photoelectron spectroscopy (XPS) is a powerful experimental technique for probing the electronic structure of molecules and materials; however, interpreting XPS data requires accurate computational methods to model core-ionized states. This work proposes and benchmarks a new approach based on the generalized active space-driven similarity renormalization group (GAS-DSRG) for calculating core-ionization energies and treating correlation effects at the perturbative and nonperturbative levels. We tested the GAS-DSRG across three data sets. First, the vertical core-ionization energies of small molecules containing first-row elements are evaluated. GAS-DSRG achieves mean absolute errors below 0.3 eV, which is comparable to high-level coupled cluster methods. Next, the accuracy of GAS-DSRG is evaluated for larger organic molecules using the CORE65 data set, with the DSRG-MRPT3 level yielding a mean absolute error of only 0.34 eV for 65 core-ionization transitions. Insights are provided into the treatment of static and dynamic correlation, the importance of high-order perturbation theory, and notable differences from density functional theory in the predicted energy ordering of core-ionized states for specific molecules. Finally, vibrationally resolved XPS spectra of diatomic molecules (CO, N2, and O2) are simulated, showing excellent agreement with experimental data.
Collapse
Affiliation(s)
- Meng Huang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
4
|
Park W, Lashkaripour A, Komarov K, Lee S, Huix-Rotllant M, Choi CH. Toward Consistent Predictions of Core/Valence Ionization Potentials and Valence Excitation Energies by MRSF-TDDFT. J Chem Theory Comput 2024; 20:5679-5694. [PMID: 38902891 DOI: 10.1021/acs.jctc.4c00640] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/22/2024]
Abstract
Optimizing exchange-correlation functionals for both core/valence ionization potentials (cIPs/vIPs) and valence excitation energies (VEEs) at the same time in the framework of MRSF-TDDFT is self-contradictory. To overcome the challenge, within the previous "adaptive exact exchange" or double-tuning strategy on Coulomb-attenuating XC functionals (CAM), a new XC functional specifically for cIPs and vIPs was first developed by enhancing exact exchange to both short- and long-range regions. The resulting DTCAM-XI functional achieved remarkably high accuracy in its predictions with errors of less than half eV. An additional concept of "valence attenuation", where the amount of exact exchange for the frontier orbital regions is selectively suppressed, was introduced to consistently predict both VEEs and IPs at the same time. The second functional, DTCAM-XIV, exhibits consistent overall prediction accuracy at ∼0.64 eV. By preferentially optimizing VEEs within the same "valence attenuation" concept, a third functional, DTCAM-VAEE, was obtained, which exhibits improved performance as compared to that of the previous DTCAM-VEE and DTCAM-AEE in the prediction of VEEs, making it an attractive alternative to BH&HLYP. As the combination of "adaptive exchange" and "valence attenuation" is operative, it would be exciting to explore its potential with a more tunable framework in the future.
Collapse
Affiliation(s)
- Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Alireza Lashkaripour
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Konstantin Komarov
- Center for Quantum Dynamics, Pohang University of Science and Technology, Pohang 37673, South Korea
- Department of Chemistry, University of Zürich, Zürich 8057, Switzerland
| | - Seunghoon Lee
- Department of Chemistry, Seoul National University, Seoul 151-747, South Korea
| | | | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
5
|
Fransson T, Pettersson LGM. TDDFT and the x-ray absorption spectrum of liquid water: Finding the "best" functional. J Chem Phys 2024; 160:234105. [PMID: 38884399 DOI: 10.1063/5.0209719] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2024] [Accepted: 06/03/2024] [Indexed: 06/18/2024] Open
Abstract
We investigate the performance of time-dependent density functional theory (TDDFT) for reproducing high-level reference x-ray absorption spectra of liquid water and water clusters. For this, we apply the integrated absolute difference (IAD) metric, previously used for x-ray emission spectra of liquid water [T. Fransson and L. G. M. Pettersson, J. Chem. Theory Comput. 19, 7333-7342 (2023)], in order to investigate which exchange-correlation (xc) functionals yield TDDFT spectra in best agreement to reference, as well as to investigate the suitability of IAD for x-ray absorption spectroscopy spectrum calculations. Considering highly asymmetric and symmetric six-molecule clusters, it is seen that long-range corrected xc-functionals are required to yield good agreement with the reference coupled cluster (CC) and algebraic-diagrammatic construction spectra, with 100% asymptotic Hartree-Fock exchange resulting in the lowest IADs. The xc-functionals with best agreement to reference have been adopted for larger water clusters, yielding results in line with recently published CC theory, but which still show some discrepancies in the relative intensity of the features compared to experiment.
Collapse
Affiliation(s)
- Thomas Fransson
- Department of Physics, AlbaNova University Center, Stockholm University, 10961 Stockholm, Sweden
| | - Lars G M Pettersson
- Department of Physics, AlbaNova University Center, Stockholm University, 10961 Stockholm, Sweden
| |
Collapse
|
6
|
Ehrman J, Shumilov K, Jenkins AJ, Kasper JM, Vitova T, Batista ER, Yang P, Li X. Unveiling Hidden Shake-Up Features in the Uranyl M 4-Edge Spectrum. JACS AU 2024; 4:1134-1141. [PMID: 38559711 PMCID: PMC10976573 DOI: 10.1021/jacsau.3c00838] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/28/2023] [Revised: 02/14/2024] [Accepted: 02/15/2024] [Indexed: 04/04/2024]
Abstract
The M4,5-edge high energy resolution X-ray absorption near-edge structure (HR-XANES) spectra of actinyls offer valuable insights into the electronic structure and bonding properties of heavy-element complexes. To conduct a comprehensive spectral analysis, it is essential to employ computational methods that accurately account for relativistic effects and electron correlation. In this work, we utilize variational relativistic multireference configurational interaction methods to compute and analyze the X-ray M4-edge absorption spectrum of uranyl. By employing these advanced computational techniques, we achieve excellent agreement between the calculated spectral features and experimental observations. Moreover, the calculations unveil significant shake-up features, which arise from the intricate interplay between strongly correlated 3d core-electron and ligand excitations. This research provides important theoretical insights into the spectral characteristics of heavy-element complexes. Furthermore, it establishes the foundation for utilizing M4,5-edge spectroscopy as a means to investigate the chemical activities of such complexes. By leveraging this technique, we can gain a deeper understanding of the bonding behavior and reactivity of heavy-element compounds.
Collapse
Affiliation(s)
- Jordan
N. Ehrman
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Kirill Shumilov
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Andrew J. Jenkins
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Joseph M. Kasper
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Tonya Vitova
- Institute
for Nuclear Waste Disposal (INE), Karlsruhe
Institute of Technology, P.O. Box 3640, Karlsruhe D-76021, Germany
| | - Enrique R. Batista
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Ping Yang
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Xiaosong Li
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
7
|
Song Q, Liu B, Wu J, Zou W, Wang Y, Suo B, Lei Y. GUGA-based MRCI approach with core-valence separation approximation (CVS) for the calculation of the core-excited states of molecules. J Chem Phys 2024; 160:094114. [PMID: 38445728 DOI: 10.1063/5.0189443] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2023] [Accepted: 02/14/2024] [Indexed: 03/07/2024] Open
Abstract
We develop and demonstrate how to use the Graphical Unitary Group Approach (GUGA)-based MRCISD with Core-Valence Separation (CVS) approximation to compute the core-excited states. First, perform a normal Self-Consistent-Field (SCF) or valence MCSCF calculation to optimize the molecular orbitals. Second, rotate the optimized target core orbitals and append to the active space, form an extended CVS active space, and perform a CVS-MCSCF calculation for core-excited states. Finally, construct the CVS-MRCISD expansion space and perform a CVS-MRCISD calculation to optimize the CI coefficients based on the variational method. The CVS approximation with GUGA-based methods can be implemented by flexible truncation of the Distinct Row Table. Eliminating the valence-excited configurations from the CVS-MRCISD expansion space can prevent variational collapse in the Davidson iteration diagonalization. The accuracy of the CVS-MRCISD scheme was investigated for excitation energies and compared with that of the CVS-MCSCF and CVS-CASPT2 methods using the same active space. The results show that CVS-MRCISD is capable of reproducing well-matched vertical core excitation energies that are consistent with experiments by combining large basis sets and a rational reference space. The calculation results also highlight the fact that the dynamic correlation between electrons makes an undeniable contribution in core-excited states.
Collapse
Affiliation(s)
- Qi Song
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, China
| | - Baoyuan Liu
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, China
| | - Junfeng Wu
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, China
| | - Wenli Zou
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, China
| | - Yubin Wang
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, China
| | - Bingbing Suo
- Institute of Modern Physics, Northwest University, Xi'an, Shaanxi 710069, China
| | - Yibo Lei
- College of Chemistry and Materials Science, Northwest University, Xi'an, Shaanxi 710069, People's Republic of China
| |
Collapse
|
8
|
Yuan X, Halbert L, Pototschnig JV, Papadopoulos A, Coriani S, Visscher L, Pereira Gomes AS. Formulation and Implementation of Frequency-Dependent Linear Response Properties with Relativistic Coupled Cluster Theory for GPU-Accelerated Computer Architectures. J Chem Theory Comput 2024; 20:677-694. [PMID: 38193434 DOI: 10.1021/acs.jctc.3c00812] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2024]
Abstract
We present the development and implementation of relativistic coupled cluster linear response theory (CC-LR), which allows the determination of molecular properties arising from time-dependent or time-independent electric, magnetic, or mixed electric-magnetic perturbations (within a common gauge origin for the magnetic properties) as well as taking into account the finite lifetime of excited states in the framework of damped response theory. We showcase our implementation, which is capable to offload the computationally intensive tensor contractions characteristic of coupled cluster theory onto graphical processing units, in the calculation of (a) frequency-(in)dependent dipole-dipole polarizabilities of IIB atoms and selected diatomic molecules, with a particular emphasis on the calculation of valence absorption cross sections for the I2 molecule; (b) indirect spin-spin coupling constants for benchmark systems such as the hydrogen halides (HX, X = F-I) as well the H2Se-H2O dimer as a prototypical system containing hydrogen bonds; and (c) optical rotations at the sodium D line for hydrogen peroxide analogues (H2Y2, Y = O, S, Se, Te). Thanks to this implementation, we are able to show the similarities in performance, but often the significant discrepancies, between CC-LR and approximate methods such as density functional theory. Comparing standard CC response theory with the flavor based upon the equation of motion formalism, we find that for valence properties such as polarizabilities, the two frameworks yield very similar results across the periodic table as found elsewhere in the literature; for properties that probe the core region, such as spin-spin couplings, on the other hand, we show a progressive differentiation between the two as relativistic effects become more important. Our results also suggest that as one goes down the periodic table, it may become increasingly difficult to measure pure optical rotation at the sodium D line due to the appearance of absorbing states.
Collapse
Affiliation(s)
- Xiang Yuan
- Univ. Lille, CNRS, UMR 8523─PhLAM─Physique des Lasers Atomes et Molécules, F-59000 Lille, France
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | - Loïc Halbert
- Univ. Lille, CNRS, UMR 8523─PhLAM─Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | - Johann Valentin Pototschnig
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | - Anastasios Papadopoulos
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | - Sonia Coriani
- DTU Chemistry─Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Lucas Visscher
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, 1081 HV Amsterdam, The Netherlands
| | | |
Collapse
|
9
|
Hait D, Martínez TJ. Predicting the X-ray Absorption Spectrum of Ozone with Single Configuration State Functions. J Chem Theory Comput 2024; 20:873-881. [PMID: 38175153 DOI: 10.1021/acs.jctc.3c01035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2024]
Abstract
X-ray absorption spectra (XAS) of biradicaloid species are often thought to represent a challenge to theoretical methods. This has led to the testing of recently developed multireference techniques on the XAS of ozone, but reproduction of the experimental spectral profile has proven difficult. We utilize a minimal model consisting of a single configuration state function (CSF) per excited state to model core-level excitations of ozone, with the orbitals of each CSF optimized using the restricted open-shell Kohn-Sham (ROKS) method. This protocol leads to semiquantitative agreement with experimental XAS. In fact, we find that low-lying core-hole excited states in biradicaloids can be approximated with individual CSFs, despite the presence of multireference character in the ground state. We also report that the 1s → π* and 1s → σ* transitions have quite distinct widths for O3. This reveals the importance of sampling over a representative range of geometries from the vibrational ground state for properly assessing the accuracy of electronic structure methods against experiments instead of the popular procedure of uniformly broadening stick spectra at the equilibrium geometry.
Collapse
Affiliation(s)
- Diptarka Hait
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
10
|
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
|
11
|
Liao C, Kasper JM, Jenkins AJ, Yang P, Batista ER, Frisch MJ, Li X. State Interaction Linear Response Time-Dependent Density Functional Theory with Perturbative Spin-Orbit Coupling: Benchmark and Perspectives. JACS AU 2023; 3:358-367. [PMID: 36873704 PMCID: PMC9975852 DOI: 10.1021/jacsau.2c00659] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/04/2022] [Revised: 01/17/2023] [Accepted: 01/18/2023] [Indexed: 06/18/2023]
Abstract
Spin-orbit coupling (SOC) is an important driving force in photochemistry. In this work, we develop a perturbative spin-orbit coupling method within the linear response time-dependent density function theory framework (TDDFT-SO). A full state interaction scheme, including singlet-triplet and triplet-triplet coupling, is introduced to describe not only the coupling between the ground and excited states, but also between excited states with all couplings between spin microstates. In addition, expressions to compute spectral oscillator strengths are presented. Scalar relativity is included variationally using the second-order Douglas-Kroll-Hess Hamiltonian, and the TDDFT-SO method is validated against variational SOC relativistic methods for atomic, diatomic, and transition metal complexes to determine the range of applicability and potential limitations. To demonstrate the robustness of TDDFT-SO for large-scale chemical systems, the UV-Vis spectrum of Au25(SR)18 - is computed and compared to experiment. Perspectives on the limitation, accuracy, and capability of perturbative TDDFT-SO are presented via analyses of benchmark calculations. Additionally, an open-source Python software package (PyTDDFT-SO) is developed and released to interface with the Gaussian 16 quantum chemistry software package to perform this calculation.
Collapse
Affiliation(s)
- Can Liao
- Department
of Chemistry, University of Washington, Seattle, Washington98195, United States
| | - Joseph M. Kasper
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
| | - Andrew J. Jenkins
- Department
of Chemistry, University of Washington, Seattle, Washington98195, United States
| | - Ping Yang
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
| | - Enrique R. Batista
- Theoretical
Division, Los Alamos National Laboratory, Los Alamos, New Mexico87545, United States
| | - Michael J. Frisch
- Gaussian
Inc., 340 Quinnipiac Street, Bldg 40, Wallingford, Connecticut06492, United States
| | - Xiaosong Li
- Department
of Chemistry, University of Washington, Seattle, Washington98195, United States
| |
Collapse
|
12
|
Datar A, Wright C, Matthews DA. Theoretical Investigation of the X-ray Stark Effect in Small Molecules. J Phys Chem A 2023; 127:1576-1587. [PMID: 36787229 DOI: 10.1021/acs.jpca.2c08311] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/15/2023]
Abstract
We have studied the Stark effect in the soft x-ray region for various small molecules by calculating the field-dependent x-ray absorption spectra. This effect is explained in terms of the response of molecular orbitals (core and valence), the molecular dipole moment, and the molecular geometry to the applied electric field. A number of consistent trends are observed linking the computed shifts in absorption energies and intensities with specific features of the molecular electronic structure. We find that both the virtual molecular orbitals (valence and/or Rydberg) as well as the core orbitals contribute to observed trends in a complementary fashion. This initial study highlights the potential impact of x-ray Stark spectroscopy as a tool to study electronic structure and environmental perturbations at a submolecular scale.
Collapse
Affiliation(s)
- Avdhoot Datar
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Catherine Wright
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| | - Devin A Matthews
- Department of Chemistry, Southern Methodist University, Dallas, Texas 75275, United States
| |
Collapse
|
13
|
Mester D, Kállay M. Double-Hybrid Density Functional Theory for Core Excitations: Theory and Benchmark Calculations. J Chem Theory Comput 2023; 19:1310-1321. [PMID: 36721871 PMCID: PMC9979613 DOI: 10.1021/acs.jctc.2c01222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
The double-hybrid (DH) time-dependent density functional theory is extended to core excitations. Two different DH formalisms are presented utilizing the core-valence separation (CVS) approximation. First, a CVS-DH variant is introduced relying on the genuine perturbative second-order correction, while an iterative analogue is also presented using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH ansatz. The performance of the new approaches is tested for the most popular DH functionals using the recently proposed XABOOM [J. Chem. Theory Comput.2021, 17, 1618] benchmark set. In order to make a careful comparison, the accuracy and precision of the methods are also inspected. Our results show that the genuine approaches are highly competitive with the more advanced CVS-ADC(2)-based methods if only excitation energies are required. In contrast, as expected, significant differences are observed in oscillator strengths; however, the precision is acceptable for the genuine functionals as well. Concerning the performance of the CVS-DH approaches, the PBE0-2/CVS-ADC(2) functional is superior, while its spin-opposite-scaled variant is also recommended as a cost-effective alternative. For these approaches, significant improvements are realized in the error measures compared with the popular CVS-ADC(2) method.
Collapse
Affiliation(s)
- Dávid Mester
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Müegyetem rkp. 3, H-1111Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111Budapest, Hungary,
| | - Mihály Kállay
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Müegyetem rkp. 3, H-1111Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111Budapest, Hungary,
| |
Collapse
|
14
|
Shaalan Alag A, Jelenfi DP, Tajti A, Szalay PG. Accurate Prediction of Vertical Ionization Potentials and Electron Affinities from Spin-Component Scaled CC2 and ADC(2) Models. J Chem Theory Comput 2022; 18:6794-6801. [PMID: 36269873 PMCID: PMC9890482 DOI: 10.1021/acs.jctc.2c00624] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
The CC2 and ADC(2) wave function models and their spin-component scaled modifications are adopted for predicting vertical ionization potentials (VIPs) and electron affinities (VEAs). The ionic solutions are obtained as electronic excitations in the continuum orbital formalism, making possible the use of existing, widespread quantum chemistry codes with minimal modifications, in full consistency with the treatment of charge transfer excitations. The performance of different variants is evaluated via benchmark calculations on various sets from previous works, containing small- and medium-sized systems, including the nucleobases. It is shown that with the spin-scaled approximate methods, in particular the scaled opposite-spin variant of the ADC(2) method, the accuracy of EOM-CCSD is achievable at a fraction of the computational cost, also outperforming many common electron propagator approaches.
Collapse
Affiliation(s)
- Ahmed Shaalan Alag
- György
Hevesy Doctoral School, Institute of Chemistry,
ELTE Eötvös Loránd University, H-1117Budapest, Hungary
| | - Dávid P. Jelenfi
- György
Hevesy Doctoral School, Institute of Chemistry,
ELTE Eötvös Loránd University, H-1117Budapest, Hungary
| | - Attila Tajti
- Laboratory
of Theoretical Chemistry, Institute of Chemistry,
ELTE Eötvös Loránd University, P.O. Box 32, H-1518Budapest 112, Hungary,E-mail:
| | - Péter G. Szalay
- Laboratory
of Theoretical Chemistry, Institute of Chemistry,
ELTE Eötvös Loránd University, P.O. Box 32, H-1518Budapest 112, Hungary
| |
Collapse
|
15
|
Park W, Alías-Rodríguez M, Cho D, Lee S, Huix-Rotllant M, Choi CH. Mixed-Reference Spin-Flip Time-Dependent Density Functional Theory for Accurate X-ray Absorption Spectroscopy. J Chem Theory Comput 2022; 18:6240-6250. [PMID: 36166346 DOI: 10.1021/acs.jctc.2c00746] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
It is demonstrated that the challenging core-hole particle (CHP) orbital relaxation for core electron spectra can be readily achieved by the mixed-reference spin-flip (MRSF)-time-dependent density functional theory (TDDFT). With the additional scalar relativistic effects on K-edge excitation energies of 24 second- and 17 third-row molecules, the particular ΔCHP-MRSF(R) exhibited near perfect predictions with RMSE ∼0.5 eV, featuring a median value of 0.3 and an interquartile range of 0.4. Overall, the CHP effect is 2-4 times stronger than relativistic ones, contributing more than 20 eV in the cases of sulfur and chlorine third-row atoms. Such high precision allows to explain the splitting and spectral shapes of O, N, and C atom K-edges in the ground state of thymine with atom as well as orbital specific accuracy. The same protocol with a double hole particle relaxation also produced remarkably accurate K-edge spectra of core to valence hole excitation energies from the first (nO8π*) and second (ππ*) excited states of thymine, confirming the assignment of 1s → n excitation for the experimentally observed 526.4 eV peak. Regarding both accuracy and practicality, therefore, MRSF-TDDFT provides a promising protocol for core electron spectra of both ground and excited electronic states alike.
Collapse
Affiliation(s)
- Woojin Park
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Marc Alías-Rodríguez
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, Marseille 13284, France
| | - Daeheum Cho
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| | - Seunghoon Lee
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, United States
| | - Miquel Huix-Rotllant
- Aix-Marseille Univ, CNRS, Institut de Chimie Radicalaire, Marseille 13284, France
| | - Cheol Ho Choi
- Department of Chemistry, Kyungpook National University, Daegu 41566, South Korea
| |
Collapse
|
16
|
Andersen JH, Nanda KD, Krylov AI, Coriani S. Cherry-Picking Resolvents: Recovering the Valence Contribution in X-ray Two-Photon Absorption within the Core-Valence-Separated Equation-of-Motion Coupled-Cluster Response Theory. J Chem Theory Comput 2022; 18:6189-6202. [PMID: 36084326 DOI: 10.1021/acs.jctc.2c00541] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Calculations of first-order response wave functions in the X-ray regime often diverge within correlated frameworks such as equation-of-motion coupled-cluster singles and doubles (EOM-CCSD), a consequence of the coupling with the valence ionization continuum. Here, we extend our strategy of introducing a hierarchy of approximations to the EOM-EE-CCSD resolvent (or, inversely, the model Hamiltonian) involved in the response equations for the calculation of X-ray two-photon absorption (X2PA) cross sections. We exploit the frozen-core core-valence separation (fc-CVS) scheme to first decouple the core and valence Fock spaces, followed by a separate approximate treatment of the valence resolvent. We demonstrate the robust convergence of X-ray response calculations within this framework and compare X2PA spectra of small benchmark molecules with the previously reported density functional theory results.
Collapse
Affiliation(s)
- Josefine H Andersen
- DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| | - Kaushik D Nanda
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Anna I Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, United States
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| |
Collapse
|
17
|
Nascimento DR, Govind N. Computational approaches for XANES, VtC-XES, and RIXS using linear-response time-dependent density functional theory based methods. Phys Chem Chem Phys 2022; 24:14680-14691. [PMID: 35699090 DOI: 10.1039/d2cp01132h] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/23/2023]
Abstract
The emergence of state-of-the-art X-ray light sources has paved the way for novel spectroscopies that take advantage of their atomic specificity to shed light on fundamental physical, chemical, and biological processes both in the static and time domains. The success of these experiments hinges on the ability to interpret and predict core-level spectra, which has opened avenues for theory to play a key role. Over the last two decades, linear-response time-dependent density functional theory (LR-TDDFT), despite various theoretical challenges, has become a computationally attractive and versatile framework to study excited-state spectra including X-ray spectroscopies. In this context, we focus our discussion on LR-TDDFT approaches for the computation of X-ray Near-Edge Structure (XANES), Valence-to-Core X-ray Emission (VtC-XES), and Resonant Inelastic X-ray Scattering (RIXS) spectroscopies in molecular systems with an emphasis on Gaussian basis set implementations. We illustrate these approaches with applications and provide a brief outlook of possible new directions.
Collapse
Affiliation(s)
- Daniel R Nascimento
- Department of Chemistry, The University of Memphis, Memphis, TN, 38152, USA.
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, USA.
| |
Collapse
|
18
|
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
|
19
|
Cunha LA, Hait D, Kang R, Mao Y, Head-Gordon M. Relativistic Orbital-Optimized Density Functional Theory for Accurate Core-Level Spectroscopy. J Phys Chem Lett 2022; 13:3438-3449. [PMID: 35412838 DOI: 10.1021/acs.jpclett.2c00578] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
Core-level spectra of 1s electrons of elements heavier than Ne show significant relativistic effects. We combine advances in orbital-optimized density functional theory (OO-DFT) with the spin-free exact two-component (X2C) model for scalar relativistic effects to study K-edge spectra of third period elements. OO-DFT/X2C is found to be quite accurate at predicting energies, yielding a ∼0.5 eV root-mean-square error versus experiment with the modern SCAN (and related) functionals. This marks a significant improvement over the >50 eV deviations that are typical for the popular time-dependent DFT (TDDFT) approach. Consequently, experimental spectra are quite well reproduced by OO-DFT/X2C, sans empirical shifts for alignment. OO-DFT/X2C combines high accuracy with ground state DFT cost and is thus a promising route for computing core-level spectra of third period elements. We also explored K and L edges of 3d transition metals to identify limitations of the OO-DFT/X2C approach in modeling the spectra of heavier atoms.
Collapse
Affiliation(s)
- Leonardo A Cunha
- 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
| | - 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
| | - Richard Kang
- 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
| | - Yuezhi Mao
- Department of Chemistry, Stanford University, Stanford, California 94305, 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
|
20
|
Yao Y, Golze D, Rinke P, Blum V, Kanai Y. All-Electron BSE@ GW Method for K-Edge Core Electron Excitation Energies. J Chem Theory Comput 2022; 18:1569-1583. [PMID: 35138865 DOI: 10.1021/acs.jctc.1c01180] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We present an accurate computational approach to calculate absolute K-edge core electron excitation energies as measured by X-ray absorption spectroscopy. Our approach employs an all-electron Bethe-Salpeter equation (BSE) formalism based on GW quasiparticle energies (BSE@GW) using numeric atom-centered orbitals (NAOs). The BSE@GW method has become an increasingly popular method for the computation of neutral valence excitation energies of molecules. However, it was so far not applied to molecular K-edge excitation energies. We discuss the influence of different numerical approximations on the BSE@GW calculation and employ in our final setup (i) exact numeric algorithms for the frequency integration of the GW self-energy, (ii) G0W0 and BSE starting points with ∼50% of exact exchange, (iii) the Tamm-Dancoff approximation and (iv) relativistic corrections. We study the basis set dependence and convergence with common Gaussian-type orbital and NAO basis sets. We identify the importance of additional spatially confined basis functions as well as of diffuse augmenting basis functions. The accuracy of our BSE@GW method is assessed for a benchmark set of small organic molecules, previously used for benchmarking the equation-of-motion coupled cluster method [Peng et al., J. Chem. Theory Comput., 2015, 11, 4146], as well as the medium-sized dibenzothiophene (DBT) molecule. Our BSE@GW results for absolute excitation energies are in excellent agreement with the experiment, with a mean average error of only 0.63 eV for the benchmark set and with errors <1 eV for the DBT molecule.
Collapse
Affiliation(s)
- Yi Yao
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Dorothea Golze
- Faculty of Chemistry and Food Chemistry, Technische Universität Dresden, 01062 Dresden, Germany.,Department of Applied Physics, Aalto University, P.O. Box 11100, FI-00076 Aalto, Finland
| | - Patrick Rinke
- Department of Applied Physics, Aalto University, P.O. Box 11100, FI-00076 Aalto, Finland
| | | | - Yosuke Kanai
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
21
|
Yang M, Sissay A, Chen M, Lopata K. Intruder Peak-Free Transient Inner-Shell Spectra Using Real-Time Simulations. J Chem Theory Comput 2022; 18:992-1002. [PMID: 35025498 DOI: 10.1021/acs.jctc.1c00079] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Real-time methods are convenient for simulating core-level absorption spectra but suffer from nonphysical intruder peaks when using atom-centered basis sets. In transient absorption spectra, these peaks exhibit highly nonphysical time-dependent modulations in their energies and oscillator strengths. In this paper, we address the origins of these intruder peaks and propose a straightforward and effective solution based on a filtered dipole operator. In combination with real-time time-dependent density functional theory (RT-TDDFT), we demonstrate how to compute intruder-free attosecond transient X-ray absorption spectra for the aminophenol (C6H7NO) oxygen and nitrogen K-edges and the α-quartz (SiO2) silicon L-edge. Without filtering, the computed spectra are qualitatively wrong. This procedure is suitable for both static and transient inner-shell spectroscopy studies and can easily be implemented in a range of real-time methodologies.
Collapse
Affiliation(s)
- Mengqi Yang
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Adonay Sissay
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Min Chen
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States.,Center for Computation and Technology, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
22
|
de Moura CEV, Sokolov AY. Simulating X-ray photoelectron spectra with strong electron correlation using multireference algebraic diagrammatic construction theory. Phys Chem Chem Phys 2022; 24:4769-4784. [DOI: 10.1039/d1cp05476g] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
A new theoretical approach for the simulations of X-ray photoelectron spectra of strongly correlated molecular systems that combines multireference algebraic diagrammatic construction theory (MR-ADC) with a core–valence separation (CVS) technique.
Collapse
Affiliation(s)
- Carlos E. V. de Moura
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio, 43210, USA
| | - Alexander Yu. Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio, 43210, USA
| |
Collapse
|
23
|
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
|
24
|
Jenkins AJ, Hu H, Lu L, Frisch MJ, Li X. Two-Component Multireference Restricted Active Space Configuration Interaction for the Computation of L-Edge X-ray Absorption Spectra. J Chem Theory Comput 2021; 18:141-150. [PMID: 34908414 DOI: 10.1021/acs.jctc.1c00564] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
X-ray absorption spectroscopy is a powerful probe of local electronic and nuclear structures, providing insights into chemical processes. The theoretical prediction and interpretation of metal L-edge X-ray absorption spectra are complicated by both relativistic effects, including spin-orbit coupling and the multiconfigurational nature of the states involved. This work details an exact two-component multireference restricted active space (RAS) configuration interaction scheme that uses an exact two-component state-averaged complete active space self-consistent-field method, which includes the spin-orbit coupling in a variational manner, for the accurate description of the electronic structure before using a RAS configuration interaction method to describe the core excited states of the X-ray spectrum. Benchmark calculations are presented for a series of iron-containing complexes, with results showing key features of the spectrum being reproduced, including ligand-to-metal charge transfer and shake-up excitations.
Collapse
Affiliation(s)
- Andrew J Jenkins
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hang Hu
- Molecular Engineering and Sciences Institute, University of Washington, Seattle, Washington 98195, United States
| | - Lixin Lu
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Michael J Frisch
- Gaussian Inc., 340 Quinnipiac Street, Building 40, Wallingford, Connecticut 06492, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
25
|
Hohenstein EG, Yu JK, Bannwarth C, List NH, Paul AC, Folkestad SD, Koch H, Martínez TJ. Predictions of Pre-edge Features in Time-Resolved Near-Edge X-ray Absorption Fine Structure Spectroscopy from Hole-Hole Tamm-Dancoff-Approximated Density Functional Theory. J Chem Theory Comput 2021; 17:7120-7133. [PMID: 34623139 DOI: 10.1021/acs.jctc.1c00478] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Time-resolved near-edge X-ray absorption fine structure (TR-NEXAFS) spectroscopy is a powerful technique for studying photochemical reaction dynamics with femtosecond time resolution. In order to avoid ambiguity in TR-NEXAFS spectra from nonadiabatic dynamics simulations, core- and valence-excited states must be evaluated on equal footing and those valence states must also define the potential energy surfaces used in the nonadiabatic dynamics simulation. In this work, we demonstrate that hole-hole Tamm-Dancoff-approximated density functional theory (hh-TDA) is capable of directly simulating TR-NEXAFS spectroscopies. We apply hh-TDA to the excited-state dynamics of acrolein. We identify two pre-edge features in the oxygen K-edge TR-NEXAFS spectrum associated with the S2 (ππ*) and S1 (nπ*) excited states. We show that these features can be used to follow the internal conversion dynamics between the lowest three electronic states of acrolein. Due to the low, O(N2) apparent computational complexity of hh-TDA and our GPU-accelerated implementation, this method is promising for the simulation of pre-edge features in TR-NEXAFS spectra of large molecules and molecules in the condensed phase.
Collapse
Affiliation(s)
- Edward G Hohenstein
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Jimmy K Yu
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States.,Biophysics Program, Stanford University, Stanford, California 94305, United States
| | - Christoph Bannwarth
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Nanna Holmgaard List
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| | - Alexander C Paul
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491 Trondheim, Norway
| | - Sarai D Folkestad
- 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
| | - Todd J Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, United States.,SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, United States
| |
Collapse
|
26
|
Park YC, Perera A, Bartlett RJ. Equation of motion coupled-cluster study of core excitation spectra II: Beyond the dipole approximation. J Chem Phys 2021; 155:094103. [PMID: 34496593 DOI: 10.1063/5.0059276] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present the time-independent (TI) and time-dependent (TD) equation of motion coupled-cluster (EOM-CC) oscillator strengths not limited to those obtained by the dipole approximation. For the conventional TI-EOM-CC, we implement all the terms in the multipole expansion through second order that contributes to the oscillator strength. These include contributions such as magnetic dipole, electric quadrupole, electric octupole, and magnetic quadrupole. In TD-EOM-CC, we only include the quadrupole moment contributions. This augments our previous work [Y. C. Park, A. Perera, and R. J. Bartlett, J. Chem. Phys. 151, 164117 (2019)]. The inclusion of the quadrupole contributions (and all the other contributions through second order in the case of TI-EOM-CCSD) enables us to obtain the intensities for the pre-edge transitions in the metal K-edge spectra, which are dipole inactive. The TI-EOM-CCSD and TD-EOM-CCSD spectra of Ti4+ atoms are used to showcase the implementation of the second-order oscillator strengths. The origin of 1s → e and 1s → t2 in core spectra from iron tetrachloride and titanium tetrachloride is discussed and compared with the experiment.
Collapse
Affiliation(s)
- Young Choon Park
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
27
|
Cooper BC, Koulias LN, Nascimento DR, Li X, DePrince AE. Short Iterative Lanczos Integration in Time-Dependent Equation-of-Motion Coupled-Cluster Theory. J Phys Chem A 2021; 125:5438-5447. [PMID: 34121405 DOI: 10.1021/acs.jpca.1c01102] [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/29/2022]
Abstract
A time-dependent (TD) formulation of equation-of-motion coupled-cluster (EOM-CC) theory can provide excited-state information over an arbitrarily wide energy window with a reduced memory footprint relative to conventional, frequency-domain EOM-CC theory. However, the floating-point costs of the time-integration required by TD-EOM-CC are generally far larger than those of the frequency-domain form of the approach. This work considers the potential of the short iterative Lanczos (SIL) integration scheme [J. Chem. Phys. 1986, 85, 5870-5876] to reduce the floating-point costs of TD-EOM-CC simulations. Low-energy and K-edge absorption features for small molecules are evaluated using TD-EOM-CC with single and double excitations, with the time-integrations carried out via SIL and fourth-order Runge-Kutta (RK4) schemes. Spectra derived from SIL- and RK4-driven simulations are nearly indistinguishable, and with an appropriately chosen subspace dimension, the SIL requires far fewer floating-point operations than are required by RK4. For K-edge spectra, SIL is the more efficient scheme by an average factor of 7.2.
Collapse
Affiliation(s)
- Brandon C Cooper
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Lauren N Koulias
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Daniel R Nascimento
- Department of Chemistry, The University of Memphis, Memphis, Tennessee 38152, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - A Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| |
Collapse
|
28
|
Halbert L, Vidal ML, Shee A, Coriani S, Severo Pereira Gomes A. Relativistic EOM-CCSD for Core-Excited and Core-Ionized State Energies Based on the Four-Component Dirac-Coulomb(-Gaunt) Hamiltonian. J Chem Theory Comput 2021; 17:3583-3598. [PMID: 33944570 DOI: 10.1021/acs.jctc.0c01203] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
We report an implementation of the core-valence separation approach to the four-component relativistic Hamiltonian-based equation-of-motion coupled-cluster with singles and doubles theory (CVS-EOM-CCSD) for the calculation of relativistic core-ionization potentials and core-excitation energies. With this implementation, which is capable of exploiting double group symmetry, we investigate the effects of the different CVS-EOM-CCSD variants and the use of different Hamiltonians based on the exact two-component (X2C) framework on the energies of different core-ionized and -excited states in halogen- (CH3I, HX, and X-, X = Cl-At) and xenon-containing (Xe, XeF2) species. Our results show that the X2C molecular mean-field approach [Sikkema, J.; J. Chem. Phys. 2009, 131, 124116], based on four-component Dirac-Coulomb mean-field calculations (2DCM), is capable of providing core excitations and ionization energies that are nearly indistinguishable from the reference four-component energies for up to and including fifth-row elements. We observe that two-electron integrals over the small-component basis sets lead to non-negligible contributions to core binding energies for the K and L edges for atoms such as iodine or astatine and that the approach based on Dirac-Coulomb-Gaunt mean-field calculations (2DCGM) are significantly more accurate than X2C calculations for which screened two-electron spin-orbit interactions are included via atomic mean-field integrals.
Collapse
Affiliation(s)
- Loïc Halbert
- CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, Université de Lille, F-59000 Lille, France
| | - Marta L Vidal
- DTU Chemistry-Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Avijit Shee
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Sonia Coriani
- DTU Chemistry-Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - André Severo Pereira Gomes
- CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, Université de Lille, F-59000 Lille, France
| |
Collapse
|
29
|
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
|
30
|
Fransson T, Brumboiu IE, Vidal ML, Norman P, Coriani S, Dreuw A. XABOOM: An X-ray Absorption Benchmark of Organic Molecules Based on Carbon, Nitrogen, and Oxygen 1s → π* Transitions. J Chem Theory Comput 2021; 17:1618-1637. [PMID: 33544612 PMCID: PMC8023667 DOI: 10.1021/acs.jctc.0c01082] [Citation(s) in RCA: 34] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Indexed: 01/05/2023]
Abstract
The performance of several standard and popular approaches for calculating X-ray absorption spectra at the carbon, nitrogen, and oxygen K-edges of 40 primarily organic molecules up to the size of guanine has been evaluated, focusing on the low-energy and intense 1s → π* transitions. Using results obtained with CVS-ADC(2)-x and fc-CVS-EOM-CCSD as benchmark references, we investigate the performance of CC2, ADC(2), ADC(3/2), and commonly adopted density functional theory (DFT)-based approaches. Here, focus is on precision rather than on accuracy of transition energies and intensities-in other words, we target relative energies and intensities and the spread thereof, rather than absolute values. The use of exchange-correlation functionals tailored for time-dependent DFT calculations of core excitations leads to error spreads similar to those seen for more standard functionals, despite yielding superior absolute energies. Long-range corrected functionals are shown to perform particularly well compared to our reference data, showing error spreads in energy and intensity of 0.2-0.3 eV and ∼10%, respectively, as compared to 0.3-0.6 eV and ∼20% for a typical pure hybrid. In comparing intensities, state mixing can complicate matters, and techniques to avoid this issue are discussed. Furthermore, the influence of basis sets in high-level ab initio calculations is investigated, showing that reasonably accurate results are obtained with the use of 6-311++G**. We name this benchmark suite as XABOOM (X-ray absorption benchmark of organic molecules) and provide molecular structures and ground-state self-consistent field energies and spectroscopic data. We believe that it provides a good assessment of electronic structure theory methods for calculating X-ray absorption spectra and will become useful for future developments in this field.
Collapse
Affiliation(s)
- Thomas Fransson
- Interdisciplinary
Center for Scientific Computing, Ruprecht-Karls
University, Im Neuenheimer
Feld 205, 69120 Heidelberg, Germany
- Fysikum, Stockholm University, Albanova, 10691 Stockholm, Sweden
| | - Iulia E. Brumboiu
- Department
of Theoretical Chemistry and Biology, KTH
Royal Institute of Technology, 10691 Stockholm, Sweden
- Department
of Chemistry, Korea Advanced Institute of
Science and Technology, 34141 Daejeon, Korea
| | - Marta L. Vidal
- DTU
Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| | - Patrick Norman
- Department
of Theoretical Chemistry and Biology, KTH
Royal Institute of Technology, 10691 Stockholm, Sweden
| | - Sonia Coriani
- DTU
Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
- Department
of Chemistry, NTNU-Norwegian University
of Science and Technology, N-7991 Trondheim, Norway
| | - Andreas Dreuw
- Interdisciplinary
Center for Scientific Computing, Ruprecht-Karls
University, Im Neuenheimer
Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
31
|
Tsuru S, Vidal ML, Pápai M, Krylov AI, Møller KB, Coriani S. An assessment of different electronic structure approaches for modeling time-resolved x-ray absorption spectroscopy. STRUCTURAL DYNAMICS (MELVILLE, N.Y.) 2021; 8:024101. [PMID: 33786337 PMCID: PMC7986275 DOI: 10.1063/4.0000070] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Accepted: 02/11/2021] [Indexed: 05/06/2023]
Abstract
We assess the performance of different protocols for simulating excited-state x-ray absorption spectra. We consider three different protocols based on equation-of-motion coupled-cluster singles and doubles, two of them combined with the maximum overlap method. The three protocols differ in the choice of a reference configuration used to compute target states. Maximum-overlap-method time-dependent density functional theory is also considered. The performance of the different approaches is illustrated using uracil, thymine, and acetylacetone as benchmark systems. The results provide guidance for selecting an electronic structure method for modeling time-resolved x-ray absorption spectroscopy.
Collapse
Affiliation(s)
- Shota Tsuru
- DTU Chemistry, Technical University of Denmark, Kemitorvet Building 207, DK-2800 Kgs. Lyngby, Denmark
| | - Marta L. Vidal
- DTU Chemistry, Technical University of Denmark, Kemitorvet Building 207, DK-2800 Kgs. Lyngby, Denmark
| | - Mátyás Pápai
- DTU Chemistry, Technical University of Denmark, Kemitorvet Building 207, DK-2800 Kgs. Lyngby, Denmark
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, California 90089, USA
| | - Klaus B. Møller
- DTU Chemistry, Technical University of Denmark, Kemitorvet Building 207, DK-2800 Kgs. Lyngby, Denmark
| | - Sonia Coriani
- DTU Chemistry, Technical University of Denmark, Kemitorvet Building 207, DK-2800 Kgs. Lyngby, Denmark
| |
Collapse
|
32
|
Bauman NP, Liu H, Bylaska EJ, Krishnamoorthy S, Low GH, Granade CE, Wiebe N, Baker NA, Peng B, Roetteler M, Troyer M, Kowalski K. Toward Quantum Computing for High-Energy Excited States in Molecular Systems: Quantum Phase Estimations of Core-Level States. J Chem Theory Comput 2021; 17:201-210. [PMID: 33332965 DOI: 10.1021/acs.jctc.0c00909] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
Abstract
This paper explores the utility of the quantum phase estimation (QPE) algorithm in calculating high-energy excited states characterized by the promotion of electrons occupying core-level shells. These states have been intensively studied over the last few decades, especially in supporting the experimental effort at light sources. Results obtained with QPE are compared with various high-accuracy many-body techniques developed to describe core-level states. The feasibility of the quantum phase estimator in identifying classes of challenging shake-up states characterized by the presence of higher-order excitation effects is discussed. We also demonstrate the utility of the QPE algorithm in targeting excitations from specific centers in a molecule. Lastly, we discuss how the lowest-order Trotter formula can be applied to reducing the complexity of the ansatz without affecting the error.
Collapse
Affiliation(s)
- Nicholas P Bauman
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Hongbin Liu
- Microsoft Quantum, Redmond, Washington 98052, United States
| | - Eric J Bylaska
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Sriram Krishnamoorthy
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Guang Hao Low
- Microsoft Quantum, Redmond, Washington 98052, United States
| | | | - Nathan Wiebe
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Nathan A Baker
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Bo Peng
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | | | | | - Karol Kowalski
- Physical Sciences Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| |
Collapse
|
33
|
Moitra T, Madsen D, Christiansen O, Coriani S. Vibrationally resolved coupled-cluster x-ray absorption spectra from vibrational configuration interaction anharmonic calculations. J Chem Phys 2020; 153:234111. [DOI: 10.1063/5.0030202] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Affiliation(s)
- Torsha Moitra
- DTU Chemistry—Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kongens Lyngby, Denmark
| | - Diana Madsen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Ove Christiansen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus C, Denmark
| | - Sonia Coriani
- DTU Chemistry—Department of Chemistry, Technical University of Denmark, Kemitorvet Bldg. 207, DK-2800 Kongens Lyngby, Denmark
| |
Collapse
|
34
|
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
|
35
|
Li X, Govind N, Isborn C, DePrince AE, Lopata K. Real-Time Time-Dependent Electronic Structure Theory. Chem Rev 2020; 120:9951-9993. [DOI: 10.1021/acs.chemrev.0c00223] [Citation(s) in RCA: 64] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195-1700, United States
| | - Niranjan Govind
- Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Christine Isborn
- Department of Chemistry and Chemical Biology, University of California, Merced, California 95343, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306-4390, United States
| | - Kenneth Lopata
- Department of Chemistry, Louisiana State University, Baton Rouge, Louisiana 70803, United States
| |
Collapse
|
36
|
Herbst MF, Fransson T. Quantifying the error of the core-valence separation approximation. J Chem Phys 2020; 153:054114. [PMID: 32770930 DOI: 10.1063/5.0013538] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
For the calculation of core-excited states probed through X-ray absorption spectroscopy, the core-valence separation (CVS) scheme has become a vital tool. This approach allows us to target such states with high specificity, albeit introducing an error. We report the implementation of a post-processing step for CVS excitations obtained within the algebraic-diagrammatic construction scheme for the polarization propagator, which removes this error. Based on this, we provide a detailed analysis of the CVS scheme, identifying its accuracy to be dominated by an error balance between two neglected couplings, one between core and valence single excitations and the other between single and double core excitations. The selection of the basis set is shown to be vital for a proper description of both couplings, with tight polarizing functions being necessary for a good balance of errors. The CVS error is confirmed to be stable across multiple systems, with an element-specific spread for K-edge spectrum calculations of only about ±0.02 eV. A systematic lowering of the CVS error by 0.02 eV-0.03 eV is noted when considering excitations to extremely diffuse states, emulating ionization.
Collapse
Affiliation(s)
- Michael F Herbst
- CERMICS, École des Ponts ParisTech, 6-8 Avenue Blaise Pascal, 77455 Marne-la-Vallée, France; Inria Paris, 75589 Paris Cedex 12, France; and Sorbonne Universitée, Institut des Sciences du Calcul et des Données, ISCD, 75005 Paris, France
| | - Thomas Fransson
- Interdisciplinary Center for Scientific Computing, Heidelberg University, 69120 Heidelberg, Germany and Fysikum, Stockholm University, Albanova, 10691 Stockholm, Sweden
| |
Collapse
|
37
|
Besley NA. Density Functional Theory Based Methods for the Calculation of X-ray Spectroscopy. Acc Chem Res 2020; 53:1306-1315. [PMID: 32613827 DOI: 10.1021/acs.accounts.0c00171] [Citation(s) in RCA: 57] [Impact Index Per Article: 14.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
The availability of new light sources combined with the realization of the unique capabilities of spectroscopy in the X-ray region has driven tremendous advances in the field of X-ray spectroscopy. Currently, these techniques are emerging as powerful analytical tools for the study of a wide range of problems encompassing liquids, materials, and biological systems. Time-resolved measurements add a further dimension to X-ray spectroscopy, opening up the potential to resolve ultrafast chemical processes at an atomic level. X-ray spectroscopy encompasses a range of techniques which provide complementary information, and these include X-ray photoelectron spectroscopy (XPS), X-ray absorption spectroscopy (XAS), X-ray emission spectroscopy (XES), and resonant inelastic X-ray scattering (RIXS). In many studies, the interpretation of the experimental data relies upon calculations to enable the nature of the underlying molecular structure, electronic structure, and bonding to be revealed. Density functional theory (DFT) based methods are some of the most widely used methods for the simulation of X-ray spectra. In this Account, we focus on our recent contributions to the simulation of a range of X-ray spectroscopic techniques using DFT and linear-response time-dependent density functional theory (TDDFT) and show how these methods can provide a computational toolkit for the simulation of X-ray spectroscopy. The importance of the exchange-correlation functional for the calculation of XAS is discussed, and the introduction of short-range corrected functionals is described. The application of these calculations to study large systems through the use of efficient implementations of TDDFT will be highlighted, with the use of these methods illustrated through studies of ionic liquids and transition metal complexes. The extension of TDDFT to calculate XES through the use of a reference determinant for the core-ionized state will be described, and the factors that affect the accuracy of the computed spectra discussed. The application of these approaches will be illustrated through the study of a range of organic molecules and transition metal complexes, which also show how going beyond the dipole approximation in determining the transition intensities can be critical. The application of these approaches to the simulation of the RIXS spectrum of water will also be described, highlighting how ultrafast dynamics on the femtoscale time scale are evident in the measured spectra. In these calculations, the description of the core-ionized and core-excited states becomes increasingly important, and the role of the basis set in accurately describing these states will be explored.
Collapse
Affiliation(s)
- Nicholas A. Besley
- School of Chemistry, University of Nottingham, University Park, Nottingham NG7 2RD, United Kingdom
| |
Collapse
|
38
|
Sarangi R, Vidal ML, Coriani S, Krylov AI. On the basis set selection for calculations of core-level states: different strategies to balance cost and accuracy. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1769872] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
Affiliation(s)
- Ronit Sarangi
- Department of Chemistry, University of Southern California, Los Angeles, CA, USA
| | - Marta L. Vidal
- DTU Chemistry – Department of Chemistry, Technical University of Denmark, Lyngby, Denmark
| | - Sonia Coriani
- DTU Chemistry – Department of Chemistry, Technical University of Denmark, Lyngby, Denmark
| | - Anna I. Krylov
- Department of Chemistry, University of Southern California, Los Angeles, CA, USA
| |
Collapse
|
39
|
Matthews DA. EOM-CC methods with approximate triple excitations applied to core excitation and ionisation energies. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1771448] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2022]
|
40
|
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
|
41
|
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
|
42
|
Kasper JM, Li X. Natural transition orbitals for complex two-component excited state calculations. J Comput Chem 2020; 41:1557-1563. [PMID: 32220083 DOI: 10.1002/jcc.26196] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2019] [Revised: 12/26/2019] [Accepted: 03/09/2020] [Indexed: 01/02/2023]
Abstract
While the natural transition orbital (NTO) method has allowed electronic excitations from time-dependent Hartree-Fock and density functional theory to be viewed in a traditional orbital picture, the extension to multicomponent molecular orbitals such as those used in relativistic two-component methods or generalized Hartree-Fock (GHF) or generalized Kohn-Sham (GKS) is less straightforward due to mixing of spin-components and the inherent inclusion of spin-flip transitions in time-dependent GHF/GKS. An extension of single-component NTOs to the two-component framework is presented, in addition to a brief discussion of the practical aspects of visualizing two-component complex orbitals. Unlike the single-component analog, the method explicitly describes the spin and frequently obtains solutions with several significant orbital pairs. The method is presented using calculations on a mercury atom and a CrO2 Cl2 complex.
Collapse
Affiliation(s)
- Joseph M Kasper
- Department of Chemistry, University of Washington, Seattle, Washington, USA
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington, USA
| |
Collapse
|
43
|
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
|
44
|
Vidal ML, Krylov AI, Coriani S. Dyson orbitals within the fc-CVS-EOM-CCSD framework: theory and application to X-ray photoelectron spectroscopy of ground and excited states. Phys Chem Chem Phys 2020; 22:2693-2703. [DOI: 10.1039/c9cp03695d] [Citation(s) in RCA: 35] [Impact Index Per Article: 8.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Ionization energies and Dyson orbitals within frozen-core core–valence separated equation-of-motion coupled cluster singles and doubles (fc-CVS-EOM-CCSD) enable efficient and reliable calculations of standard XPS and of UV-pump/XPS probe spectra.
Collapse
Affiliation(s)
- Marta L. Vidal
- DTU Chemistry – Department of Chemistry
- Technical University of Denmark
- Kongens Lyngby
- Denmark
| | - Anna I. Krylov
- Department of Chemistry
- University of Southern California
- Los Angeles
- USA
| | - Sonia Coriani
- DTU Chemistry – Department of Chemistry
- Technical University of Denmark
- Kongens Lyngby
- Denmark
| |
Collapse
|
45
|
Nanda KD, Vidal ML, Faber R, Coriani S, Krylov AI. How to stay out of trouble in RIXS calculations within equation-of-motion coupled-cluster damped response theory? Safe hitchhiking in the excitation manifold by means of core–valence separation. Phys Chem Chem Phys 2020; 22:2629-2641. [DOI: 10.1039/c9cp03688a] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
We present a novel approach with robust convergence of the response equations for computing resonant inelastic X-ray scattering (RIXS) cross sections within the equation-of-motion coupled-cluster (EOM-CC) framework.
Collapse
Affiliation(s)
- Kaushik D. Nanda
- Department of Chemistry
- University of Southern California
- Los Angeles
- USA
| | - Marta L. Vidal
- DTU Chemistry – Department of Chemistry
- Technical University of Denmark
- DK-2800
- Denmark
| | - Rasmus Faber
- DTU Chemistry – Department of Chemistry
- Technical University of Denmark
- DK-2800
- Denmark
| | - Sonia Coriani
- DTU Chemistry – Department of Chemistry
- Technical University of Denmark
- DK-2800
- Denmark
| | - Anna I. Krylov
- Department of Chemistry
- University of Southern California
- Los Angeles
- USA
- The Hamburg Centre for Ultrafast Imaging
| |
Collapse
|
46
|
Myhre RH, Coriani S, Koch H. X-ray and UV Spectra of Glycine within Coupled Cluster Linear Response Theory. J Phys Chem A 2019; 123:9701-9711. [PMID: 31549830 DOI: 10.1021/acs.jpca.9b06590] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The coupled cluster models CCSD and CC3 are used to investigate the (core) excited states and ionization energies of glycine in the gas phase. Excited states and ionization energies in the UV spectral range are calculated using a standard coupled cluster linear response, while core-level excited states and ionization potentials are calculated using the core-valence separation approximation. The temperature dependence from different conformers is also assessed.
Collapse
Affiliation(s)
- Rolf H Myhre
- Department of Chemistry , Norwegian University of Science and Technology, NTNU , 7491 Trondheim , Norway.,Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry , University of Oslo , 0315 Oslo , Norway
| | - Sonia Coriani
- DTU Chemistry , Technical University of Denmark , DK-2800 Kongens Lyngby , Denmark.,Aarhus Institute of Advanced Studies , Aarhus University , DK-8000 Aarhus C , Denmark
| | - Henrik Koch
- Department of Chemistry , Norwegian University of Science and Technology, NTNU , 7491 Trondheim , Norway.,Scuola Normale Superiore , Piazza dei Cavalieri 7 , 56126 Pisa , Italy
| |
Collapse
|
47
|
Park YC, Perera A, Bartlett RJ. Equation of motion coupled-cluster for core excitation spectra: Two complementary approaches. J Chem Phys 2019; 151:164117. [PMID: 31675901 DOI: 10.1063/1.5117841] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
This paper presents core excitation spectra from coupled-cluster (CC) theory obtained from both a time-independent and a new time-dependent formalism. The conventional time-independent CC formulation for excited states is the equation-of-motion (EOM-CC) method whose eigenvalues and eigenvectors describe the core excited states. An alternative computational procedure is offered by a time-dependent CC description. In that case, the dipole transition operator is expressed in the CC effective Hamiltonian form and propagated with respect to time. The absorption spectrum is obtained from the CC dipole autocorrelation function via a Fourier transformation. Comparisons are made among the time-dependent results obtained from second-order perturbation theory, to coupled cluster doubles and their linearized forms (CCD and LCCD), to CC singles and doubles (CCSD) and the linearized form (LCCSD). In the time-independent case, considerations of triples (EOM-CCSDT) and quadruples (EOM-CCSDTQ) are used to approach sub-electron volt accuracy. A particular target is the allyl radical, as an example of an open-shell molecule. As the results have to ultimately be the same, the two procedures offer a complementary approach toward analyzing experimental results.
Collapse
Affiliation(s)
- Young Choon Park
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Ajith Perera
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Rodney J Bartlett
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| |
Collapse
|
48
|
Koulias LN, Williams-Young DB, Nascimento DR, DePrince AE, Li X. Relativistic Real-Time Time-Dependent Equation-of-Motion Coupled-Cluster. J Chem Theory Comput 2019; 15:6617-6624. [DOI: 10.1021/acs.jctc.9b00729] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Lauren N. Koulias
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - David B. Williams-Young
- Computational Research Division, Lawrence Berkeley National Laboratory, 1 Cyclotron Road MS 50A-3111, Berkeley, California 94720, United States
| | - Daniel R. Nascimento
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry, Florida State University, Tallahassee, Florida 32306, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
49
|
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
|
50
|
Bokarev SI, Kühn O. Theoretical X‐ray spectroscopy of transition metal compounds. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1433] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Affiliation(s)
| | - Oliver Kühn
- Institut für Physik Universität Rostock Rostock Germany
| |
Collapse
|