1
|
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
|
2
|
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
|
3
|
Arias-Martinez JE, Wu H, Head-Gordon M. Generalization of One-Center Nonorthogonal Configuration Interaction Singles to Open-Shell Singlet Reference States: Theory and Application to Valence-Core Pump-Probe States in Acetylacetone. J Chem Theory Comput 2024; 20:752-766. [PMID: 38164934 DOI: 10.1021/acs.jctc.3c01139] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2024]
Abstract
We formulate a one-center nonorthogonal configuration interaction singles (1C-NOCIS) theory for the computation of core excited states of an initial singlet state with two unpaired electrons. This model, which we refer to as 1C-NOCIS two-electron open-shell (2eOS), is appropriate for computing the K-edge near-edge X-ray absorption spectra (NEXAS) of the valence excited states of closed-shell molecules relevant to pump-probe time-resolved (TR) NEXAS experiments. With the inclusion of core-hole relaxation effects and explicit spin adaptation, 1C-NOCIS 2eOS requires mild shifts to match experiment, is free of artifacts due to spin contamination, and can capture the high-energy region of the spectrum beyond the transitions into the singly occupied molecular orbitals (SOMOs). Calculations on water and thymine illustrate the different key features of excited-state NEXAS, namely, the core-to-SOMO transitions as well as shifts and spin-splittings in the transitions analogous to those of the ground state. Simulations of the TR-NEXAS of acetylacetone after excitation to its π → π* singlet excited state at the carbon K-edge, an experiment carried out recently, showcase the ability of 1C-NOCIS 2eOS to efficiently simulate NEXAS based on nonadiabatic molecular dynamics simulations.
Collapse
Affiliation(s)
- Juan E Arias-Martinez
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Hamlin Wu
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
4
|
Haasler M, Maier TM, Kaupp M. Toward a correct treatment of core properties with local hybrid functionals. J Comput Chem 2023; 44:2461-2477. [PMID: 37635647 DOI: 10.1002/jcc.27211] [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: 07/04/2023] [Revised: 08/01/2023] [Accepted: 08/03/2023] [Indexed: 08/29/2023]
Abstract
In local hybrid functionals (LHs), a local mixing function (LMF) determines the position-dependent exact-exchange admixture. We report new LHs that focus on an improvement of the LMF in the core region while retaining or partly improving upon the high accuracy in the valence region exhibited by the LH20t functional. The suggested new pt-LMFs are based on a Padé form and modify the previously used ratio between von Weizsäcker and Kohn-Sham local kinetic energies by different powers of the density to enable flexibly improved approximations to the correct high-density and iso-orbital limits relevant for the innermost core region. Using TDDFT calculations for a set of K-shell core excitations of second- and third-period systems including accurate state-of-the-art relativistic orbital corrections, the core part of the LMF is optimized, while the valence part is optimized as previously reported for test sets of atomization energies and reaction barriers (Haasler et al., J Chem Theory Comput 2020, 16, 5645). The LHs are completed by a calibration function that minimizes spurious nondynamical correlation effects caused by the gauge ambiguities of exchange-energy densities, as well as by B95c meta-GGA correlation. The resulting LH23pt functional relates to the previous LH20t functional but specifically improves upon the core region.
Collapse
Affiliation(s)
- Matthias Haasler
- Technische Universität Berlin, Institute of Chemistry Theoretical Chemistry/Quantum Chemistry, Berlin, Germany
| | - Toni M Maier
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Braunschweig, Germany
| | - Martin Kaupp
- Technische Universität Berlin, Institute of Chemistry Theoretical Chemistry/Quantum Chemistry, Berlin, Germany
| |
Collapse
|
5
|
Fransson T, Pettersson LGM. Calibrating TDDFT Calculations of the X-ray Emission Spectrum of Liquid Water: The Effects of Hartree-Fock Exchange. J Chem Theory Comput 2023; 19:7333-7342. [PMID: 37787584 PMCID: PMC10601479 DOI: 10.1021/acs.jctc.3c00728] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Indexed: 10/04/2023]
Abstract
The structure and dynamics of liquid water continue to be debated, with insight provided by, among others, X-ray emission spectroscopy (XES), which shows a split in the high-energy 1b1 feature. This split is yet to be reproduced by theory, and it remains unclear if these difficulties are related to inaccuracies in dynamics simulations, spectrum calculations, or both. We investigate the performance of different methods for calculating XES of liquid water, focusing on the ability of time-dependent density functional theory (TDDFT) to reproduce reference spectra obtained by high-level coupled cluster and algebraic-diagrammatic construction scheme calculations. A metric for evaluating the agreement between theoretical spectra termed the integrated absolute difference (IAD), which considers the integral of shifted difference spectra, is introduced and used to investigate the performance of different exchange-correlation functionals. We find that computed spectra of symmetric and asymmetric model water structures are strongly and differently influenced by the amount of Hartree-Fock exchange, with best agreement to reference spectra for ∼40-50%. Lower percentages tend to yield high density of contributing states, resulting in too broad features. The method introduced here is useful also for other spectrum calculations, in particular where the performance for ensembles of structures are evaluated.
Collapse
Affiliation(s)
- Thomas Fransson
- Department of Physics, AlbaNova
University Center, Stockholm University, 109 61 Stockholm, Sweden
| | - Lars G. M. Pettersson
- Department of Physics, AlbaNova
University Center, Stockholm University, 109 61 Stockholm, Sweden
| |
Collapse
|
6
|
Jana S, Herbert JM. Fractional-Electron and Transition-Potential Methods for Core-to-Valence Excitation Energies Using Density Functional Theory. J Chem Theory Comput 2023; 19:4100-4113. [PMID: 37312236 DOI: 10.1021/acs.jctc.3c00202] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
Methods for computing X-ray absorption spectra based on a constrained core hole (possibly containing a fractional electron) are examined. These methods are based on Slater's transition concept and its generalizations, wherein core-to-valence excitation energies are determined using Kohn-Sham orbital energies. Methods examined here avoid promoting electrons beyond the lowest unoccupied molecular orbital, facilitating robust convergence. Variants of these ideas are systematically tested, revealing a best-case accuracy of 0.3-0.4 eV (with respect to experiment) for K-edge transition energies. Absolute errors are much larger for higher-lying near-edge transitions but can be reduced below 1 eV by introducing an empirical shift based on a charge-neutral transition-potential method, in conjunction with functionals such as SCAN, SCAN0, or B3LYP. This procedure affords an entire excitation spectrum from a single fractional-electron calculation, at the cost of ground-state density functional theory and without the need for state-by-state calculations. This shifted transition-potential approach may be especially useful for simulating transient spectroscopies or in complex systems where excited-state Kohn-Sham calculations are challenging.
Collapse
Affiliation(s)
- Subrata Jana
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
7
|
Yang KR, Kyro GW, Batista VS. The landscape of computational approaches for artificial photosynthesis. NATURE COMPUTATIONAL SCIENCE 2023; 3:504-513. [PMID: 38177419 DOI: 10.1038/s43588-023-00450-1] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Accepted: 04/11/2023] [Indexed: 01/06/2024]
Abstract
Artificial photosynthesis is an attractive strategy for converting solar energy into fuels, largely because the Earth receives enough solar energy in one hour to meet humanity's energy needs for an entire year. However, developing devices for artificial photosynthesis remains difficult and requires computational approaches to guide and assist the interpretation of experiments. In this Perspective, we discuss current and future computational approaches, as well as the challenges of designing and characterizing molecular assemblies that absorb solar light, transfer electrons between interfaces, and catalyze water-splitting and fuel-forming reactions.
Collapse
Affiliation(s)
- Ke R Yang
- Department of Chemistry, Yale University, New Haven, CT, USA
- Energy Sciences Institute, Yale University, West Haven, CT, USA
| | - Gregory W Kyro
- Department of Chemistry, Yale University, New Haven, CT, USA
| | - Victor S Batista
- Department of Chemistry, Yale University, New Haven, CT, USA.
- Energy Sciences Institute, Yale University, West Haven, CT, USA.
| |
Collapse
|
8
|
Troß J, Carter-Fenk K, Cole-Filipiak NC, Schrader P, Word M, McCaslin LM, Head-Gordon M, Ramasesha K. Excited-State Dynamics during Primary C-I Homolysis in Acetyl Iodide Revealed by Ultrafast Core-Level Spectroscopy. J Phys Chem A 2023; 127:4103-4114. [PMID: 37103479 DOI: 10.1021/acs.jpca.3c01414] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/28/2023]
Abstract
In typical carbonyl-containing molecules, bond dissociation events follow initial excitation to nπC═O* states. However, in acetyl iodide, the iodine atom gives rise to electronic states with mixed nπC═O* and nσC-I* character, leading to complex excited-state dynamics, ultimately resulting in dissociation. Using ultrafast extreme ultraviolet (XUV) transient absorption spectroscopy and quantum chemical calculations, we present an investigation of the primary photodissociation dynamics of acetyl iodide via time-resolved spectroscopy of core-to-valence transitions of the I atom after 266 nm excitation. The probed I 4d-to-valence transitions show features that evolve on sub-100-fs time scales, reporting on excited-state wavepacket evolution during dissociation. These features subsequently evolve to yield spectral signatures corresponding to free iodine atoms in their spin-orbit ground and excited states with a branching ratio of 1.1:1 following dissociation of the C-I bond. Calculations of the valence excitation spectrum via equation-of-motion coupled cluster with single and double substitutions (EOM-CCSD) show that initial excited states are of spin-mixed character. From the initially pumped spin-mixed state, we use a combination of time-dependent density functional theory (TDDFT)-driven nonadiabatic ab initio molecular dynamics and EOM-CCSD calculations of the N4,5 edge to reveal a sharp inflection point in the transient XUV signal that corresponds to rapid C-I homolysis. By examining the molecular orbitals involved in the core-level excitations at and around this inflection point, we are able to piece together a detailed picture of C-I bond photolysis in which d → σ* transitions give way to d → p excitations as the bond dissociates. We also report theoretical predictions of short-lived, weak 4d → 5d transitions in acetyl iodide, validated by weak bleaching in the experimental transient XUV spectra. This joint experimental-theoretical effort has thus unraveled the detailed electronic structure and dynamics of a strongly spin-orbit coupled system.
Collapse
Affiliation(s)
- Jan Troß
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, United States
| | - Kevin Carter-Fenk
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Neil C Cole-Filipiak
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, United States
| | - Paul Schrader
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, United States
| | - Mi'Kayla Word
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, United States
| | - Laura M McCaslin
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, 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
| | - Krupa Ramasesha
- Combustion Research Facility, Sandia National Laboratories, Livermore, California 94550, United States
| |
Collapse
|
9
|
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
|