1
|
Mandal A, Berquist EJ, Herbert JM. A new parameterization of the DFT/CIS method with applications to core-level spectroscopy. J Chem Phys 2024; 161:044114. [PMID: 39051834 DOI: 10.1063/5.0220535] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2024] [Accepted: 07/08/2024] [Indexed: 07/27/2024] Open
Abstract
Time-dependent density functional theory (TD-DFT) within a restricted excitation space is an efficient means to compute core-level excitation energies using only a small subset of the occupied orbitals. However, core-to-valence excitation energies are significantly underestimated when standard exchange-correlation functionals are used, which is partly traceable to systemic issues with TD-DFT's description of Rydberg and charge-transfer excited states. To mitigate this, we have implemented an empirically modified combination of configuration interaction with single substitutions (CIS) based on Kohn-Sham orbitals, which is known as "DFT/CIS." This semi-empirical approach is well-suited for simulating x-ray near-edge spectra, as it contains sufficient exact exchange to model charge-transfer excitations yet retains DFT's low-cost description of dynamical electron correlation. Empirical corrections to the matrix elements enable semi-quantitative simulation of near-edge x-ray spectra without the need for significant a posteriori shifts; this should be useful in complex molecules and materials with multiple overlapping x-ray edges. Parameter optimization for use with a specific range-separated hybrid functional makes this a black-box method intended for both core and valence spectroscopy. Results herein demonstrate that realistic K-edge absorption and emission spectra can be obtained for second- and third-row elements and 3d transition metals, with promising results for L-edge spectra as well. DFT/CIS calculations require absolute shifts that are considerably smaller than what is typical in TD-DFT.
Collapse
Affiliation(s)
- Aniket Mandal
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | | | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
2
|
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
|
3
|
Folkestad SD, Paul AC, Paul Née Matveeva R, Coriani S, Odelius M, Iannuzzi M, Koch H. Understanding X-ray absorption in liquid water using triple excitations in multilevel coupled cluster theory. Nat Commun 2024; 15:3551. [PMID: 38670938 PMCID: PMC11053016 DOI: 10.1038/s41467-024-47690-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2023] [Accepted: 04/03/2024] [Indexed: 04/28/2024] Open
Abstract
X-ray absorption (XA) spectroscopy is an essential experimental tool to investigate the local structure of liquid water. Interpretation of the experiment poses a significant challenge and requires a quantitative theoretical description. High-quality theoretical XA spectra require reliable molecular dynamics simulations and accurate electronic structure calculations. Here, we present the first successful application of coupled cluster theory to model the XA spectrum of liquid water. We overcome the computational limitations on system size by employing a multilevel coupled cluster framework for large molecular systems. Excellent agreement with the experimental spectrum is achieved by including triple excitations in the wave function and using molecular structures from state-of-the-art path-integral molecular dynamics. We demonstrate that an accurate description of the electronic structure within the first solvation shell is sufficient to successfully model the XA spectrum of liquid water within the multilevel framework. Furthermore, we present a rigorous charge transfer analysis of the XA spectrum, which is reliable due to the accuracy and robustness of the electronic structure methodology. This analysis aligns with previous studies regarding the character of the prominent features of the XA spectrum of liquid water.
Collapse
Affiliation(s)
- Sarai Dery Folkestad
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491, Trondheim, Norway
| | - Alexander C Paul
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491, Trondheim, Norway
| | - Regina Paul Née Matveeva
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491, Trondheim, Norway
| | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, DTU, 2800, Kongens Lyngby, Denmark
| | - Michael Odelius
- Department of Physics, Stockholm University, 10691, Stockholm, Sweden
| | - Marcella Iannuzzi
- Department of Chemistry, University of Zurich, 8057, Zürich, Switzerland
| | - Henrik Koch
- Department of Chemistry, Norwegian University of Science and Technology, NTNU, 7491, Trondheim, Norway.
| |
Collapse
|
4
|
Caricato M. A Perspective on the Simulation of Electronic Circular Dichroism and Circularly Polarized Luminescence Spectra in Chiral Solid Materials. J Phys Chem A 2024; 128:1197-1206. [PMID: 38295762 DOI: 10.1021/acs.jpca.3c08095] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/23/2024]
Abstract
Chiral materials have shown tremendous potential for many technological applications, such as optoelectronics, sensing, magnetism, information technology, and imaging. Characterization of these materials is mostly based on chiroptical spectroscopies, such as electronic circular dichroism (ECD) and circularly polarized luminescence (CPL). These experimental measurements would greatly benefit from theoretical simulations for interpretation of the spectra as well as predictions on new materials. While ECD and CPL simulations are well established for molecular systems, they are not for materials. In this Perspective, we describe the theoretical quantities necessary to simulate ECD and CPL spectra in oriented systems. Then, we discuss the approximate strategies currently used to perform these calculations, what computational machinery is already available to develop more general approaches, and some of the open challenges for the simulation of ECD and CPL spectra in solid materials. When methods that are as reliable and computationally efficient as those for molecules are developed, these simulations will provide invaluable insight and guidance for the rational design of optically active materials.
Collapse
Affiliation(s)
- Marco Caricato
- Department of Chemistry, University of Kansas, 1567 Irving Hill Road, Lawrence, Kansas 66045, United States
| |
Collapse
|
5
|
Herbert JM, Zhu Y, Alam B, Ojha AK. Time-Dependent Density Functional Theory for X-ray Absorption Spectra: Comparing the Real-Time Approach to Linear Response. J Chem Theory Comput 2023; 19:6745-6760. [PMID: 37708349 DOI: 10.1021/acs.jctc.3c00673] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/16/2023]
Abstract
We simulate X-ray absorption spectra at elemental K-edges using time-dependent density functional theory (TDDFT) in both its conventional linear-response implementation and its explicitly time-dependent or "real-time" formulation. Real-time TDDFT simulations enable broadband spectra calculations without the need to invoke frozen occupied orbitals ("core/valence separation"), but we find that these spectra are often contaminated by transitions to the continuum that originate from lower-energy core and semicore orbitals. This problem becomes acute in triple-ζ basis sets, although it is sometimes sidestepped in double-ζ basis sets. Transitions to the continuum acquire surprisingly large dipole oscillator strengths, leading to spectra that are difficult to interpret. Meaningful spectra can be recovered by means of a filtering technique that decomposes the spectrum into contributions from individual occupied orbitals, and the same procedure can be used to separate L- and K-edge spectra arising from different elements within a given molecule. In contrast, conventional linear-response TDDFT requires core/valence separation but is free of these artifacts. It is also significantly more efficient than the real-time approach, even when hundreds of individual states are needed to reproduce near-edge absorption features and even when Padé approximants are used to reduce the real-time simulations to just 2-4 fs of time propagation. Despite the cost, the real-time approach may be useful to examine the validity of the core/valence separation approximation.
Collapse
Affiliation(s)
- John M Herbert
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
- Chemical Physics Graduate Program, The Ohio State University, Columbus, Ohio 43210, United States
| | - Ying Zhu
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
- Chemical Physics Graduate Program, The Ohio State University, Columbus, Ohio 43210, United States
| | - Bushra Alam
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Avik Kumar Ojha
- Department of Chemistry & Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
6
|
He ZH, Huang YY, Ji GF, Chen J, Wu Q. The Structure Properties of Carbon Materials Formed in 2,4,6-Triamino-1,3,5-Trinitrobenzene Detonation: A Theoretical Insight for Nucleation of Diamond-like Carbon. Int J Mol Sci 2023; 24:12568. [PMID: 37628750 PMCID: PMC10454052 DOI: 10.3390/ijms241612568] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2023] [Revised: 08/03/2023] [Accepted: 08/06/2023] [Indexed: 08/27/2023] Open
Abstract
The structure and properties of nano-carbon materials formed in explosives detonation are always a challenge, not only for the designing and manufacturing of these materials but also for clearly understanding the detonation performance of explosives. Herein, we study the dynamic evolution process of condensed-phase carbon involved in 2,4,6-Triamino-1,3,5-trinitrobenzene (TATB) detonation using the quantum-based molecular dynamics method. Various carbon structures such as, graphene-like, diamond-like, and "diaphite", are obtained under different pressures. The transition from a C sp2- to a sp3-hybrid, driven by the conversion of a hexatomic to a non-hexatomic ring, is detected under high pressure. A tightly bound nucleation mechanism for diamond-like carbon dominated by a graphene-like carbon layer is uncovered. The graphene-like layer is readily constructed at the early stage, which would connect with surrounding carbon atoms or fragments to form the tetrahedral structure, with a high fraction of sp3-hybridized carbon. After that, the deformed carbon layers further coalesce with each other by bonding between carbon atoms within the five-member ring, to form the diamond-like nucleus. The complex "diaphite" configuration is detected during the diamond-like carbon nucleation, which illustrates that the nucleation and growth of detonation nano-diamond would accompany the intergrowth of graphene-like layers.
Collapse
Affiliation(s)
- Zheng-Hua He
- National Key Laboratory of Shock Wave and Detonation Physics, Institute of Fluid Physics, China Academy of Engineering Physics, Mianyang 621900, China; (Z.-H.H.)
| | - Yao-Yao Huang
- National Key Laboratory of Shock Wave and Detonation Physics, Institute of Fluid Physics, China Academy of Engineering Physics, Mianyang 621900, China; (Z.-H.H.)
| | - Guang-Fu Ji
- National Key Laboratory of Shock Wave and Detonation Physics, Institute of Fluid Physics, China Academy of Engineering Physics, Mianyang 621900, China; (Z.-H.H.)
| | - Jun Chen
- National Key Laboratory of Computational Physics, Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
| | - Qiang Wu
- National Key Laboratory of Shock Wave and Detonation Physics, Institute of Fluid Physics, China Academy of Engineering Physics, Mianyang 621900, China; (Z.-H.H.)
| |
Collapse
|
7
|
Carter-Fenk K, Cunha LA, Arias-Martinez JE, Head-Gordon M. Electron-Affinity Time-Dependent Density Functional Theory: Formalism and Applications to Core-Excited States. J Phys Chem Lett 2022; 13:9664-9672. [PMID: 36215404 DOI: 10.1021/acs.jpclett.2c02564] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
The lack of particle-hole attraction and orbital relaxation within time-dependent density functional theory (TDDFT) lead to extreme errors in the prediction of K-edge X-ray absorption spectra (XAS). We derive a linear-response formalism that uses optimized orbitals of the n - 1-electron system as the reference, building orbital relaxation and a proper hole into the initial density. Our approach is an exact generalization of the static-exchange approximation that ameliorates the particle-hole interaction error associated with the adiabatic approximation and reduces errors in TDDFT XAS by orders of magnitude. With a statistical performance of just 0.5 eV root-mean-square error and the same computational scaling as TDDFT under the core-valence separation approximation, we anticipate that this approach will be of great utility in XAS calculations of large systems.
Collapse
Affiliation(s)
- Kevin Carter-Fenk
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California94720, United States
| | - Leonardo A Cunha
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California94720, United States
| | - Juan E Arias-Martinez
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| |
Collapse
|
8
|
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: 1] [Impact Index Per Article: 0.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
|
9
|
Mališ M, Vandaele E, Luber S. Spin-Orbit Couplings for Nonadiabatic Molecular Dynamics at the ΔSCF Level. J Chem Theory Comput 2022; 18:4082-4094. [PMID: 35666703 DOI: 10.1021/acs.jctc.1c01046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
A procedure for the calculation of spin-orbit coupling (SOC) at the delta self-consistent field (ΔSCF) level of theory is presented. Singlet and triplet excited electronic states obtained with the ΔSCF method are expanded into a linear combination of singly excited Slater determinants composed of ground electronic state Kohn-Sham orbitals. This alleviates the nonorthogonality between excited and ground electronic states and introduces a framework, similar to the auxiliary wave function at the time-dependent density functional theory (TD-DFT) level, for the calculation of observables. The ΔSCF observables of the formaldehyde system were compared to reference TD-DFT values. Our procedure gives all components (energies, gradients, nonadiabatic couplings, and SOC terms) at the ΔSCF level of theory for conducting efficient, full-atomistic nonadiabatic molecular dynamics with intersystem crossing, particularly in condensed phase systems.
Collapse
Affiliation(s)
- Momir Mališ
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Eva Vandaele
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| | - Sandra Luber
- Department of Chemistry, University of Zürich, Winterthurerstrasse 190, 8057 Zürich, Switzerland
| |
Collapse
|
10
|
Liang J, Feng X, Hait D, Head-Gordon M. Revisiting the Performance of Time-Dependent Density Functional Theory for Electronic Excitations: Assessment of 43 Popular and Recently Developed Functionals from Rungs One to Four. J Chem Theory Comput 2022; 18:3460-3473. [PMID: 35533317 DOI: 10.1021/acs.jctc.2c00160] [Citation(s) in RCA: 54] [Impact Index Per Article: 27.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023]
Abstract
In this paper, the performance of more than 40 popular or recently developed density functionals is assessed for the calculation of 463 vertical excitation energies against the large and accurate QuestDB benchmark set. For this purpose, the Tamm-Dancoff approximation offers a good balance between computational efficiency and accuracy. The functionals ωB97X-D and BMK are found to offer the best performance overall with a root-mean square error (RMSE) of around 0.27 eV, better than the computationally more demanding CIS(D) wave function method with a RMSE of 0.36 eV. The results also suggest that Jacob's ladder still holds for time-dependent density functional theory excitation energies, though hybrid meta generalized-gradient approximations (meta-GGAs) are not generally better than hybrid GGAs. Effects of basis set convergence, gauge invariance correction to meta-GGAs, and nonlocal correlation (VV10) are also studied, and practical basis set recommendations are provided.
Collapse
Affiliation(s)
- Jiashu Liang
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California at Berkeley, Berkeley, California 94720, United States
| | - Xintian Feng
- Q-Chem Inc., Pleasanton, California 94588, United States
| | - Diptarka Hait
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California at Berkeley, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California at Berkeley, Berkeley, California 94720, United States.,Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
11
|
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
|
12
|
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
|
13
|
Carter-Fenk K, Head-Gordon M. On the choice of reference orbitals for linear-response calculations of solution-phase K-edge X-ray absorption spectra. Phys Chem Chem Phys 2022; 24:26170-26179. [DOI: 10.1039/d2cp04077h] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
X-ray absorption spectra of liquids calculated with linear-response theories like TDDFT and CIS are dramatically improved with core-ion reference orbitals.
Collapse
Affiliation(s)
- Kevin Carter-Fenk
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, CA 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA
| |
Collapse
|
14
|
Taherivardanjani S, Elfgen R, Reckien W, Suarez E, Perlt E, Kirchner B. Benchmarking the Computational Costs and Quality of Vibrational Spectra from Ab Initio Simulations. ADVANCED THEORY AND SIMULATIONS 2021. [DOI: 10.1002/adts.202100293] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Shima Taherivardanjani
- Mulliken Center for Theoretical Chemistry Institute for Physical and Theoretical Chemistry Beringstr. 4 Bonn D‐53115 Germany
| | - Roman Elfgen
- Mulliken Center for Theoretical Chemistry Institute for Physical and Theoretical Chemistry Beringstr. 4 Bonn D‐53115 Germany
| | - Werner Reckien
- Mulliken Center for Theoretical Chemistry Institute for Physical and Theoretical Chemistry Beringstr. 4 Bonn D‐53115 Germany
| | - Estela Suarez
- Institute for Advanced Simulation Jülich Supercomputing Centre, Forschungszentrum Jülich GmbH Wilhelm‐Johnen‐Straße Jülich D‐52425 Germany
| | - Eva Perlt
- Otto Schott Institute of Materials Research Faculty of Physics and Astronomy Friedrich‐Schiller‐Universität Jena Löbdergraben 32 Jena D‐07743 Germany
| | - Barbara Kirchner
- Mulliken Center for Theoretical Chemistry Institute for Physical and Theoretical Chemistry Beringstr. 4 Bonn D‐53115 Germany
| |
Collapse
|
15
|
Bussy A, Hutter J. First-principles correction scheme for linear-response time-dependent density functional theory calculations of core electronic states. J Chem Phys 2021; 155:034108. [PMID: 34293885 DOI: 10.1063/5.0058124] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Linear-response time-dependent density functional theory (LR-TDDFT) for core level spectroscopy using standard local functionals suffers from self-interaction error and a lack of orbital relaxation upon creation of the core hole. As a result, LR-TDDFT calculated x-ray absorption near edge structure spectra needed to be shifted along the energy axis to match experimental data. We propose a correction scheme based on many-body perturbation theory to calculate the shift from first-principles. The ionization potential of the core donor state is first computed and then substituted for the corresponding Kohn-Sham orbital energy, thus emulating Koopmans's condition. Both self-interaction error and orbital relaxation are taken into account. The method exploits the localized nature of core states for efficiency and integrates seamlessly in our previous implementation of core level LR-TDDFT, yielding corrected spectra in a single calculation. We benchmark the correction scheme on molecules at the K- and L-edges as well as for core binding energies and report accuracies comparable to higher order methods. We also demonstrate applicability in large and extended systems and discuss efficient approximations.
Collapse
Affiliation(s)
- Augustin Bussy
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| | - Jürg Hutter
- Department of Chemistry, University of Zurich, Winterthurerstrasse 190, CH-8057 Zürich, Switzerland
| |
Collapse
|