1
|
Borrelli W, Mei KJ, Park SJ, Schwartz BJ. Partial Molar Solvation Volume of the Hydrated Electron Simulated Via DFT. J Phys Chem B 2024; 128:2425-2431. [PMID: 38422045 PMCID: PMC10945486 DOI: 10.1021/acs.jpcb.3c05091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2023] [Revised: 02/03/2024] [Accepted: 02/14/2024] [Indexed: 03/02/2024]
Abstract
Different simulation models of the hydrated electron produce different solvation structures, but it has been challenging to determine which simulated solvation structure, if any, is the most comparable to experiment. In a recent work, Neupane et al. [J. Phys. Chem. B 2023, 127, 5941-5947] showed using Kirkwood-Buff theory that the partial molar volume of the hydrated electron, which is known experimentally, can be readily computed from an integral over the simulated electron-water radial distribution function. This provides a sensitive way to directly compare the hydration structure of different simulation models of the hydrated electron with experiment. Here, we compute the partial molar volume of an ab-initio-simulated hydrated electron model based on density-functional theory (DFT) with a hybrid functional at different simulated system sizes. We find that the partial molar volume of the DFT-simulated hydrated electron is not converged with respect to the system size for simulations with up to 128 waters. We show that even at the largest simulation sizes, the partial molar volume of DFT-simulated hydrated electrons is underestimated by a factor of 2 with respect to experiment, and at the standard 64-water size commonly used in the literature, DFT-based simulations underestimate the experimental solvation volume by a factor of ∼3.5. An extrapolation to larger box sizes does predict the experimental partial molar volume correctly; however, larger system sizes than those explored here are currently intractable without the use of machine-learned potentials. These results bring into question what aspects of the predicted hydrated electron radial distribution function, as calculated by DFT-based simulations with the PBEh-D3 functional, deviate from the true solvation structure.
Collapse
Affiliation(s)
| | | | - Sanghyun J. Park
- Department of Chemistry and
Biochemistry, University of California,
Los Angeles, Los Angeles, California 90095-1569, United States
| | - Benjamin J. Schwartz
- Department of Chemistry and
Biochemistry, University of California,
Los Angeles, Los Angeles, California 90095-1569, United States
| |
Collapse
|
2
|
Sopena Moros A, Li S, Li K, Doumy G, Southworth SH, Otolski C, Schaller RD, Kumagai Y, Rubensson JE, Simon M, Dakovski G, Kunnus K, Robinson JS, Hampton CY, Hoffman DJ, Koralek J, Loh ZH, Santra R, Inhester L, Young L. Tracking Cavity Formation in Electron Solvation: Insights from X-ray Spectroscopy and Theory. J Am Chem Soc 2024; 146:3262-3269. [PMID: 38270463 PMCID: PMC10859959 DOI: 10.1021/jacs.3c11857] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2024]
Abstract
We present time-resolved X-ray absorption spectra of ionized liquid water and demonstrate that OH radicals, H3O+ ions, and solvated electrons all leave distinct X-ray-spectroscopic signatures. Particularly, this allows us to characterize the electron solvation process through a tool that focuses on the electronic response of oxygen atoms in the immediate vicinity of a solvated electron. Our experimental results, supported by ab initio calculations, confirm the formation of a cavity in which the solvated electron is trapped. We show that the solvation dynamics are governed by the magnitude of the random structural fluctuations present in water. As a consequence, the solvation time is highly sensitive to temperature and to the specific way the electron is injected into water.
Collapse
Affiliation(s)
- Arturo Sopena Moros
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestraße 85, Hamburg 22607, Germany
| | - Shuai Li
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Kai Li
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
- Department of Physics and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
| | - Gilles Doumy
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Stephen H Southworth
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Christopher Otolski
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Richard D Schaller
- Center for Nanoscale Materials, Argonne National Laboratory, Lemont, Illinous 60439, United States
- Department of Chemistry, Northwestern University, 2145 N. Sheridan Rd., Evanston, Illinois 60208, United States
| | - Yoshiaki Kumagai
- Department of Applied Physics, Tokyo University of Agriculture and Technology, Tokyo 184-8588, Japan
| | - Jan-Erik Rubensson
- Department of Physics and Astronomy, Uppsala University, Box 516, Uppsala SE-75120, Sweden
| | - Marc Simon
- Laboratoire de Chimie Physique-Matière et Rayonnement, LCPMR, Sorbonne Université, CNRS, Paris F-75005, France
| | | | | | | | | | | | - Jake Koralek
- LCLS, SLAC, Menlo Park, California 94025, United States
| | - Zhi-Heng Loh
- School of Chemistry, Chemical Engineering and Biotechnology, and School of Physical and Mathematical Sciences, Nanyang Technological University, Singapore 637371, Singapore
| | - Robin Santra
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestraße 85, Hamburg 22607, Germany
- Department of Physics, Universität Hamburg, Notkestraße 9, Hamburg 22607, Germany
| | - Ludger Inhester
- Center for Free-Electron Laser Science CFEL, Deutsches Elektronen-Synchrotron DESY, Notkestraße 85, Hamburg 22607, Germany
| | - Linda Young
- Chemical Sciences and Engineering Division, Argonne National Laboratory, Lemont, Illinois 60439, United States
- Department of Physics and James Franck Institute, The University of Chicago, Chicago, Illinois 60637, United States
| |
Collapse
|
3
|
Muchova E, Hollas D, Holland DMP, Bacellar C, Leroy L, Barillot TR, Longetti L, Coreno M, de Simone M, Grazioli C, Chergui M, Ingle RA. Jahn-Teller effects in initial and final states: high-resolution X-ray absorption, photoelectron and Auger spectroscopy of allene. Phys Chem Chem Phys 2023; 25:6733-6745. [PMID: 36799466 DOI: 10.1039/d2cp05299g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/18/2023]
Abstract
Carbon K-edge resonant Auger spectra of gas-phase allene following excitation of the pre-edge 1s → π* transitions are presented and analysed with the support of EOM-CCSD/cc-pVTZ calculations. X-Ray absorption (XAS), X-ray photoelectron (XPS), valence band and non-resonant Auger spectra are also reanalysed with a series of computational approaches. The results presented demonstrate the importance of including nuclear ensemble effects for simulating X-ray observables and as an effective strategy for capturing Jahn-Teller effects in spectra.
Collapse
Affiliation(s)
- Eva Muchova
- Department of Physical Chemistry, University of Chemistry and Technology, Technická 5, 166 28 Prague, Czech Republic
| | - Daniel Hollas
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, BS8 1TS, UK
| | | | - Camila Bacellar
- Laboratoire de Spectroscopie Ultrarapide, Ecole Polytechnique Fédérale de Lausanne, ISIC, FSB, CH-1015 Lausanne, Switzerland
| | - Ludmila Leroy
- Laboratoire de Spectroscopie Ultrarapide, Ecole Polytechnique Fédérale de Lausanne, ISIC, FSB, CH-1015 Lausanne, Switzerland
| | - Thomas R Barillot
- Laboratoire de Spectroscopie Ultrarapide, Ecole Polytechnique Fédérale de Lausanne, ISIC, FSB, CH-1015 Lausanne, Switzerland
| | - Luca Longetti
- Laboratoire de Spectroscopie Ultrarapide, Ecole Polytechnique Fédérale de Lausanne, ISIC, FSB, CH-1015 Lausanne, Switzerland
| | - Marcello Coreno
- ISM-CNR, Istituto di Struttura dei Materiali, LD2 Unit, 34149 Trieste, Italy
| | | | - Cesare Grazioli
- IOM-CNR, Istituto Officina dei Materiali, 34149 Trieste, Italy
| | - Majed Chergui
- Laboratoire de Spectroscopie Ultrarapide, Ecole Polytechnique Fédérale de Lausanne, ISIC, FSB, CH-1015 Lausanne, Switzerland
| | - Rebecca A Ingle
- Department of Chemistry, 20 Gordon Street, London, WC1H 0AJ, UK.
| |
Collapse
|
4
|
Carter-Fenk K, Johnson BA, Herbert JM, Schenter GK, Mundy CJ. Birth of the Hydrated Electron via Charge-Transfer-to-Solvent Excitation of Aqueous Iodide. J Phys Chem Lett 2023; 14:870-878. [PMID: 36657160 DOI: 10.1021/acs.jpclett.2c03460] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
A primary means to generate hydrated electrons in laboratory experiments is excitation to the charge-transfer-to-solvent (CTTS) state of a solute such as I-(aq), but this initial step in the genesis of e-(aq) has never been simulated directly using ab initio molecular dynamics. We report the first such simulations, combining ground- and excited-state simulations of I-(aq) with a detailed analysis of fluctuations in the Coulomb potential experienced by the nascent solvated electron. What emerges is a two-step picture of the evolution of e-(aq) starting from the CTTS state: I-(aq) + hν → I-*(aq) → I•(aq) + e-(aq). Notably, the equilibrated ground state of e-(aq) evolves from I-*(aq) without any nonadiabatic transitions, simply as a result of solvent reorganization. The methodology used here should be applicable to other photochemical electron transfer processes in solution, an important class of problems directly relevant to photocatalysis and energy transfer.
Collapse
Affiliation(s)
- Kevin Carter-Fenk
- Physical Science Division, Pacific Northwest National Laboratory, Richland, Washington99352, United States
- Department of Chemistry, University of California, Berkeley, California94720, United States
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio43210, United States
| | - Britta A Johnson
- Physical Science Division, Pacific Northwest National Laboratory, Richland, Washington99352, United States
| | - John M Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio43210, United States
| | - Gregory K Schenter
- Physical Science Division, Pacific Northwest National Laboratory, Richland, Washington99352, United States
| | - Christopher J Mundy
- Physical Science Division, Pacific Northwest National Laboratory, Richland, Washington99352, United States
- Department of Chemical Engineering, University of Washington, Seattle, Washington98195, United States
| |
Collapse
|
5
|
Park SJ, Narvaez WA, Schwartz BJ. Ab Initio Studies of Hydrated Electron/Cation Contact Pairs: Hydrated Electrons Simulated with Density Functional Theory Are Too Kosmotropic. J Phys Chem Lett 2023; 14:559-566. [PMID: 36630724 DOI: 10.1021/acs.jpclett.2c03705] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
We have performed the first DFT-based ab initio MD simulations of a hydrated electron (eaq-) in the presence of Na+, a system chosen because ion-pairing behavior in water depends sensitively on the local hydration structure. Experiments show that eaq-'s interact weakly with Na+; the eaq-'s spectrum blue shifts by only a few tens of meV upon ion pairing without changing shape. We find that the spectrum of the DFT-simulated eaq- red shifts and changes shape upon interaction with Na+, in contrast with experiment. We show that this is because the hydration structure of the DFT-simulated eaq- is too ordered or kosmotropic. Conversely, simulations that produce eaq-'s with a less ordered or chaotropic hydration structure form weaker ion pairs with Na+, yielding predicted spectral blue shifts in better agreement with experiment. Thus, ab initio simulations based on hybrid GGA DFT functionals fail to produce the correct solvation structure for the hydrated electron.
Collapse
Affiliation(s)
- Sanghyun J Park
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| | - Wilberth A Narvaez
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| | - Benjamin J Schwartz
- Department of Chemistry and Biochemistry, University of California, Los Angeles, Los Angeles, California 90095-1569, United States
| |
Collapse
|