1
|
Markmann V, Pan J, Hansen BL, Haubro ML, Nimmrich A, Lenzen P, Levantino M, Katayama T, Adachi SI, Gorski-Bilke S, Temps F, Dohn AO, Møller KB, Nielsen MM, Haldrup K. Real-time structural dynamics of the ultrafast solvation process around photo-excited aqueous halides. Chem Sci 2024; 15:11391-11401. [PMID: 39055005 PMCID: PMC11268492 DOI: 10.1039/d4sc01912a] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2024] [Accepted: 06/11/2024] [Indexed: 07/27/2024] Open
Abstract
This work investigates and describes the structural dynamics taking place following charge-transfer-to-solvent photo-abstraction of electrons from I- and Br- ions in aqueous solution following single- and 2-photon excitation at 202 nm and 400 nm, respectively. A Time-Resolved X-ray Solution Scattering (TR-XSS) approach with direct sensitivity to the structure of the surrounding solvent as the water molecules adopt a new equilibrium configuration following the electron-abstraction process is utilized to investigate the structural dynamics of solvent shell expansion and restructuring in real-time. The structural sensitivity of the scattering data enables a quantitative evaluation of competing models for the interaction between the nascent neutral species and surrounding water molecules. Taking the I0-O distance as the reaction coordinate, we find that the structural reorganization is delayed by 0.1 ps with respect to the photoexcitation and completes on a time scale of 0.5-1 ps. On longer time scales we determine from the evolution of the TR-XSS difference signal that I0: e- recombination takes place on two distinct time scales of ∼20 ps and 100 s of picoseconds. These dynamics are well captured by a simple model of diffusive evolution of the initial photo-abstracted electron population where the charge-transfer-to-solvent process gives rise to a broad distribution of electron ejection distances, a significant fraction of which are in the close vicinity of the nascent halogen atoms and recombine on short time scales.
Collapse
Affiliation(s)
- Verena Markmann
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Jaysree Pan
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Bianca L Hansen
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Morten L Haubro
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Amke Nimmrich
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
- Department of Chemistry and Molecular Biology, University of Gothenburg Gothenburg Sweden
| | - Philipp Lenzen
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Matteo Levantino
- European Synchrotron Radiation Facility CS40220 Grenoble 38043 Cedex 9 France
| | - Tetsuo Katayama
- Japan Synchrotron Radiation Research Institute Kouto 1-1-1, Sayo Hyogo 679-5198 Japan
- RIKEN SPring-8 Center 1-1-1 Kouto, Sayo Hyogo 679-5148 Japan
| | - Shin-Ichi Adachi
- Institute of Materials Structure Science, High Energy Accelerator Research Organization (KEK) 1-1 Oho, Tsukuba Ibaraki 305-0801 Japan
- Department of Materials Structure Science, School of High Energy Accelerator Science 1-1 Oho, Tsukuba Ibaraki 305-0801 Japan
| | | | - Friedrich Temps
- Christian-Albrechts-University Kiel Olshausenstr. 40 24098 Kiel Germany
| | - Asmus O Dohn
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
- Science Institute, University of Iceland 107 Reykjavík Iceland
| | - Klaus B Møller
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Martin M Nielsen
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| | - Kristoffer Haldrup
- Technical University of Denmark Anker Engelunds Vej 1 2800 Lyngby Denmark
| |
Collapse
|
2
|
Reidelbach M, Bai M, Schneeberger M, Zöllner MS, Kubicek K, Kirchberg H, Bressler C, Thorwart M, Herrmann C. Solvent Dynamics of Aqueous Halides before and after Photoionization. J Phys Chem B 2023; 127:1399-1413. [PMID: 36728132 DOI: 10.1021/acs.jpcb.2c07992] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
Electron transfer reactions can be strongly influenced by solvent dynamics. We study the photoionization of halides in water as a model system for such reactions. There are no internal nuclear degrees of freedom in the solute, allowing the dynamics of the solvent to be uniquely identified. We simulate the equilibrium solvent dynamics for Cl-, Br-, I-, and their respective neutral atoms in water, comparing quantum mechanical/molecular mechanical (QM/MM) and classical molecular dynamics (MD) methods. On the basis of the obtained configurations, we calculate the extended X-ray absorption fine structure (EXAFS) spectra rigorously based on the MD snapshots and compare them in detail with other theoretical and experimental results available in the literature. We find our EXAFS spectra based on QM/MM MD simulations in good agreement with their experimental counterparts for the ions. Classical MD simulations for the ions lead to EXAFS spectra that agree equally well with the experiment when it comes to the oscillatory period of the signal, even though they differ from the QM/MM radial distribution functions extracted from the MD. The amplitude is, however, considerably overestimated. This suggests that to judge the reliability of theoretical simulation methods or to elucidate fine details of the atomistic dynamics of the solvent based on EXAFS spectra, the amplitude as well as the oscillatory period need to be considered. If simulations fail qualitatively, as does the classical MD for the aqueous neutral halogen atoms, the resulting EXAFS will also be strongly affected in both oscillatory period and amplitude. The good reliability of QM/MM-based EXAFS simulations, together with clear qualitative differences in the EXAFS spectra found between halides and their atomic counterparts, suggests that a combined theory and experimental EXAFS approach is suitable for elucidating the nonequilibrium solvent dynamics in the photoionization of halides and possibly also for electron transfer reactions in more complex systems.
Collapse
Affiliation(s)
- Marco Reidelbach
- Department of Chemistry, Universität Hamburg, Harbor Bldg. 610, Luruper Chaussee 149, 22761Hamburg, Germany.,The Hamburg Centre of Ultrafast Imaging, Luruper Chaussee 149, 22761Hamburg, Germany
| | - Mei Bai
- The Hamburg Centre of Ultrafast Imaging, Luruper Chaussee 149, 22761Hamburg, Germany.,I. Institut für Theoretische Physik, Universität Hamburg, Notkestr. 9, 22607Hamburg, Germany
| | - Michaela Schneeberger
- Department of Chemistry, Universität Hamburg, Harbor Bldg. 610, Luruper Chaussee 149, 22761Hamburg, Germany.,The Hamburg Centre of Ultrafast Imaging, Luruper Chaussee 149, 22761Hamburg, Germany
| | - Martin Sebastian Zöllner
- Department of Chemistry, Universität Hamburg, Harbor Bldg. 610, Luruper Chaussee 149, 22761Hamburg, Germany.,The Hamburg Centre of Ultrafast Imaging, Luruper Chaussee 149, 22761Hamburg, Germany
| | - Katharina Kubicek
- The Hamburg Centre of Ultrafast Imaging, Luruper Chaussee 149, 22761Hamburg, Germany.,Department of Physics, Universität Hamburg, Notkestr. 85, 22607Hamburg, Germany.,European XFEL, Holzkoppel 4, 22869Schenefeld, Germany
| | - Henning Kirchberg
- The Hamburg Centre of Ultrafast Imaging, Luruper Chaussee 149, 22761Hamburg, Germany.,I. Institut für Theoretische Physik, Universität Hamburg, Notkestr. 9, 22607Hamburg, Germany
| | - Christian Bressler
- The Hamburg Centre of Ultrafast Imaging, Luruper Chaussee 149, 22761Hamburg, Germany.,Department of Physics, Universität Hamburg, Notkestr. 85, 22607Hamburg, Germany.,European XFEL, Holzkoppel 4, 22869Schenefeld, Germany
| | - Michael Thorwart
- The Hamburg Centre of Ultrafast Imaging, Luruper Chaussee 149, 22761Hamburg, Germany.,I. Institut für Theoretische Physik, Universität Hamburg, Notkestr. 9, 22607Hamburg, Germany
| | - Carmen Herrmann
- Department of Chemistry, Universität Hamburg, Harbor Bldg. 610, Luruper Chaussee 149, 22761Hamburg, Germany.,The Hamburg Centre of Ultrafast Imaging, Luruper Chaussee 149, 22761Hamburg, Germany
| |
Collapse
|
3
|
Gomez DT, Pratt LR, Asthagiri DN, Rempe SB. Hydrated Anions: From Clusters to Bulk Solution with Quasi-Chemical Theory. Acc Chem Res 2022; 55:2201-2212. [PMID: 35829622 PMCID: PMC9386901 DOI: 10.1021/acs.accounts.2c00078] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023]
Abstract
The interactions of hydrated ions with molecular and macromolecular solution and interface partners are strong on a chemical energy scale. Here, we recount the foremost ab initio theory for the evaluation of the hydration free energies of ions, namely, quasi-chemical theory (QCT). We focus on anions, particularly halides but also the hydroxide anion, because they have been outstanding challenges for all theories. For example, this work supports understanding the high selectivity for F- over Cl- in fluoride-selective ion channels despite the identical charge and the size similarity of these ions. QCT is built by the identification of inner-shell clusters, separate treatment of those clusters, and then the integration of those results into the broader-scale solution environment. Recent work has focused on a close comparison with mass-spectrometric measurements of ion-hydration equilibria. We delineate how ab initio molecular dynamics (AIMD) calculations on ion-hydration clusters, elementary statistical thermodynamics, and electronic structure calculations on cluster structures sampled from the AIMD calculations obtain just the free energies extracted from the cluster experiments. That theory-experiment comparison has not been attempted before the work discussed here, but the agreement is excellent with moderate computational effort. This agreement reinforces both theory and experiment and provides a numerically accurate inner-shell contribution to QCT. The inner-shell complexes involving heavier halides display strikingly asymmetric hydration clusters. Asymmetric hydration structures can be problematic for the evaluation of the QCT outer-shell contribution with the polarizable continuum model (PCM). Nevertheless, QCT provides a favorable setting for the exploitation of PCM when the inner-shell material shields the ion from the outer solution environment. For the more asymmetrically hydrated, and thus less effectively shielded, heavier halide ions clustered with waters, the PCM is less satisfactory. We therefore investigate an inverse procedure in which the inner-shell structures are sampled from readily available AIMD calculations on the bulk solutions. This inverse procedure is a remarkable improvement; our final results are in close agreement with a standard tabulation of hydration free energies, and the final composite results are independent of the coordination number on the chemical energy scale of relevance, as they should be. Finally, a comparison of anion hydration structure in clusters and bulk solutions from AIMD simulations emphasize some differences: the asymmetries of bulk solution inner-shell structures are moderated compared with clusters but are still present, and inner hydration shells fill to slightly higher average coordination numbers in bulk solution than in clusters.
Collapse
Affiliation(s)
- Diego T. Gomez
- Department
of Chemical & Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, United States,
| | - Lawrence R. Pratt
- Department
of Chemical & Biomolecular Engineering, Tulane University, New Orleans, Louisiana 70118, United States,
| | - Dilipkumar N. Asthagiri
- Department
of Chemical and Biomolecular Engineering, Rice University, Houston, Texas 77005, United States,
| | - Susan B. Rempe
- Center
for Integrated Nanotechnologies, Sandia
National Laboratories, Albuquerque, New Mexico 87185, United States,
| |
Collapse
|
4
|
Orabi EA, Öztürk TN, Bernhardt N, Faraldo-Gómez JD. Corrections in the CHARMM36 Parametrization of Chloride Interactions with Proteins, Lipids, and Alkali Cations, and Extension to Other Halide Anions. J Chem Theory Comput 2021; 17:6240-6261. [PMID: 34516741 DOI: 10.1021/acs.jctc.1c00550] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The nonpolarizable CHARMM force field is one of the most widely used energy functions for all-atom biomolecular simulations. Chloride is the only halide ion included in the latest version, CHARMM36m, and is used widely in simulation studies, often as an electrolyte ion but also as the biological substrate of transport proteins and enzymes. Here, we find that existing parameters systematically underestimate the interaction of Cl- with proteins and lipids. Accordingly, when examined in solution, little to no Cl-association can be observed with most components of the protein, including backbone, polar side chains and aromatic rings. The strength of the interaction with cationic side chains and with alkali ions is also incongruent with experimental measurements, specifically osmotic coefficients of concentrated solutions. Consistent with these findings, a 4-μs trajectory of the Cl--specific transport protein CLC-ec1 shows irreversible Cl- dissociation from the so-called Scen binding site, even in a 150 mM NaCl buffer. To correct for these deficiencies, we formulate a series of pair-specific Lennard-Jones parameters that override those resulting from the conventional Lorentz-Berthelot combination rules. These parameters, referred to as NBFIX, are systematically calibrated against available experimental data as well as ab initio geometry optimizations and energy evaluations, for a wide set of binary and ternary Cl- complexes with protein and lipid analogs and alkali cations. Analogously, we also formulate parameter sets for the other three biological halide ions, namely, fluoride, bromide, and iodide. The resulting parameters are used to calculate the potential of mean force defining the interaction of each anion and each of the protein and lipid analogues in bulk water, revealing association free energies in the range of -0.3 to -3.3 kcal/mol, with the F- complexes being the least stable. The NBFIX corrections also preserve the Cl- occupancy of CLC-ec1 in a second 4-μs trajectory. We posit that these optimized molecular-mechanics models provide a more realistic foundation for all-atom simulation studies of processes entailing changes in hydration, recognition, or transport of halide anions.
Collapse
Affiliation(s)
- Esam A Orabi
- Theoretical Molecular Biophysics Laboratory, National Heart, Lung and Blood Institute, National Institutes of Health, 10 Center Drive, Bethesda, Maryland 20814, United States
| | - Tuǧba N Öztürk
- Theoretical Molecular Biophysics Laboratory, National Heart, Lung and Blood Institute, National Institutes of Health, 10 Center Drive, Bethesda, Maryland 20814, United States.,Department of Biochemistry and Molecular Biophysics, Washington University School of Medicine, St. Louis, Missouri 63110, United States
| | - Nathan Bernhardt
- Theoretical Molecular Biophysics Laboratory, National Heart, Lung and Blood Institute, National Institutes of Health, 10 Center Drive, Bethesda, Maryland 20814, United States
| | - José D Faraldo-Gómez
- Theoretical Molecular Biophysics Laboratory, National Heart, Lung and Blood Institute, National Institutes of Health, 10 Center Drive, Bethesda, Maryland 20814, United States
| |
Collapse
|
5
|
Ravikumar B, Mynam M, Repaka S, Rai B. Solvation shell dynamics explains charge transport characteristics of LIB electrolytes. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116613] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
6
|
Paul SK, Herbert JM. Probing Interfacial Effects on Ionization Energies: The Surprising Banality of Anion-Water Hydrogen Bonding at the Air/Water Interface. J Am Chem Soc 2021; 143:10189-10202. [PMID: 34184532 DOI: 10.1021/jacs.1c03131] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Liquid microjet photoelectron spectroscopy is an increasingly common technique to measure vertical ionization energies (VIEs) of aqueous solutes, but the interpretation of these experiments is subject to questions regarding sensitivity to bulk versus interfacial solvation environments. We have computed aqueous-phase VIEs for a set of inorganic anions, using a combination of molecular dynamics simulations and electronic structure calculations, with results that are in excellent agreement with experiment regardless of whether the simulation data are restricted to ions at the air/water interface or to those in bulk aqueous solution. Although the computed VIEs are sensitive to ion-water hydrogen bonding, we find that the short-range solvation structure is sufficiently similar in both environments that it proves impossible to discriminate between the two on the basis of the VIE, a conclusion that has important implications for the interpretation of liquid-phase photoelectron spectroscopy. More generally, analysis of the simulation data suggests that the surface activity of soft anions is largely a second or third solvation shell effect, arising from disruption of water-water hydrogen bonds and not from significant changes in first-shell anion-water hydrogen bonding.
Collapse
Affiliation(s)
- Suranjan K Paul
- 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
|
Borkowski AK, Piskulich ZA, Thompson WH. Examining the Hofmeister Series through Activation Energies: Water Diffusion in Aqueous Alkali-Halide Solutions. J Phys Chem B 2021; 125:350-359. [PMID: 33382267 DOI: 10.1021/acs.jpcb.0c09965] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The effect of ions on the properties of aqueous solutions is often categorized in terms of the Hofmeister series that ranks them from chaotropes ("structure-breakers"), which weaken the surrounding hydrogen-bond network to kosmotropes ("structure-makers"), which enhance it. Here, we investigate the Hofmeister series in ∼1 M sodium-halide solutions using molecular dynamics simulations to calculate the effect of the identity and proximity of the halide anion on both the water diffusion coefficient and its activation energy. A recently developed method for calculating the activation energy from a single-temperature simulation is used, which also permits a rigorous decomposition into contributions from different interactions and motions. The mechanisms of the salt effects on the water dynamics are explored by separately considering water molecules based on their location relative to the ions. The results show that water diffusion is accelerated moving down the halide group from F- to I-. The behavior of the diffusion activation energy, Ea, is more complex, indicating a significant role for entropic effects. However, water molecules in the first or second solvation shell of an ion exhibit a decrease in Ea moving down the halide series and Na+ exhibits a larger effect than any of the anions. The Ea for water molecules within the second solvation shell of an ion are modest, indicating a short-ranged nature of the ion influence.
Collapse
Affiliation(s)
- Ashley K Borkowski
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Zeke A Piskulich
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| | - Ward H Thompson
- Department of Chemistry, University of Kansas, Lawrence, Kansas 66045, United States
| |
Collapse
|
8
|
Muralidharan A, Pratt L, Chaudhari M, Rempe S. Quasi-chemical theory for anion hydration and specific ion effects: Cl-(aq) vs. F-(aq). ACTA ACUST UNITED AC 2019. [DOI: 10.1016/j.cpletx.2019.100037] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022]
|
9
|
Zhang Y, Zeng W, Huang L, Liu W, Jia E, Zhao Y, Wang F, Zhu Z. In Situ Liquid Secondary Ion Mass Spectrometry: A Surprisingly Soft Ionization Process for Investigation of Halide Ion Hydration. Anal Chem 2019; 91:7039-7046. [PMID: 30950268 DOI: 10.1021/acs.analchem.8b05804] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
The understanding of ion solvation phenomena is of significance due to their influences on many important chemical, biological, and environmental processes. Mass spectrometry (MS)-based methods have been used to investigate this topic with molecular insights. As ion-solvent interactions are weak, ionization processes should be as soft as possible in order to retain solvation structures. An in situ liquid secondary ion MS (SIMS) approach developed in our group has been recently utilized in investigations of Li ion solvation in nonaqueous solution, and it detected a series of solvated Li ions. As traditionally SIMS has long been recognized as a hard ionization process with strong damage occurring at the sputtering interface, it is very interesting to study further how soft in situ liquid SIMS can be. In this work, we used halide ion hydration as an example to compare the ionization performance of the in situ liquid SIMS approach with regular electrospray ionization MS (ESI-MS). Results show that, although ESI has been recognized as a soft ionization method, nearly no solvated halide ions were detected by ESI-MS analysis, which acquired only strong signals of salt ion clusters. In contrast, in liquid SIMS spectra, a series of obvious hydrated halide ion compositions could be observed. We further evaluated the hydration numbers of halide ions and revealed the effects of the ion size, charge density, and polarizability on the hydration phenomenon. Our findings demonstrated that the in situ liquid SIMS approach is surprisingly soft, and it is expected to have very broad applications on investigation of various ion-solvent interactions and many other interesting chemical processes (e.g., the initial nucleation of nanoparticle formation) in liquid environment.
Collapse
Affiliation(s)
- Yanyan Zhang
- Beijing National Laboratory for Molecular Sciences, National Centre for Mass Spectrometry in Beijing, CAS Key Laboratory of Analytical Chemistry for Living Biosystems , Institute of Chemistry, Chinese Academy of Sciences , Beijing 100190 , China.,Environmental Molecular Sciences Laboratory , Pacific Northwest National Laboratory , Richland , Washington 99354 , United States.,University of Chinese Academy of Sciences , Beijing 100049 , China
| | - Wenjuan Zeng
- Beijing National Laboratory for Molecular Sciences, National Centre for Mass Spectrometry in Beijing, CAS Key Laboratory of Analytical Chemistry for Living Biosystems , Institute of Chemistry, Chinese Academy of Sciences , Beijing 100190 , China.,University of Chinese Academy of Sciences , Beijing 100049 , China
| | - Liuqin Huang
- Environmental Molecular Sciences Laboratory , Pacific Northwest National Laboratory , Richland , Washington 99354 , United States
| | - Wen Liu
- Environmental Molecular Sciences Laboratory , Pacific Northwest National Laboratory , Richland , Washington 99354 , United States
| | - Endong Jia
- University of Chinese Academy of Sciences , Beijing 100049 , China.,The Key Laboratory of Solar Thermal Energy and Photovoltaic System , Institute of Electrical Engineering, Chinese Academy of Sciences , Beijing 100190 , China
| | - Yao Zhao
- Beijing National Laboratory for Molecular Sciences, National Centre for Mass Spectrometry in Beijing, CAS Key Laboratory of Analytical Chemistry for Living Biosystems , Institute of Chemistry, Chinese Academy of Sciences , Beijing 100190 , China
| | - Fuyi Wang
- Beijing National Laboratory for Molecular Sciences, National Centre for Mass Spectrometry in Beijing, CAS Key Laboratory of Analytical Chemistry for Living Biosystems , Institute of Chemistry, Chinese Academy of Sciences , Beijing 100190 , China.,University of Chinese Academy of Sciences , Beijing 100049 , China
| | - Zihua Zhu
- Environmental Molecular Sciences Laboratory , Pacific Northwest National Laboratory , Richland , Washington 99354 , United States
| |
Collapse
|
10
|
Holden ZC, Rana B, Herbert JM. Analytic gradient for the QM/MM-Ewald method using charges derived from the electrostatic potential: Theory, implementation, and application to ab initio molecular dynamics simulation of the aqueous electron. J Chem Phys 2019; 150:144115. [DOI: 10.1063/1.5089673] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/30/2023] Open
Affiliation(s)
- Zachary C. Holden
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Bhaskar Rana
- 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
|
11
|
Fedkin MV, Shin YK, Dasgupta N, Yeon J, Zhang W, van Duin D, van Duin ACT, Mori K, Fujiwara A, Machida M, Nakamura H, Okumura M. Development of the ReaxFF Methodology for Electrolyte-Water Systems. J Phys Chem A 2019; 123:2125-2141. [PMID: 30775922 DOI: 10.1021/acs.jpca.8b10453] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A new ReaxFF reactive force field has been developed for water-electrolyte systems including cations Li+, Na+, K+, and Cs+ and anions F-, Cl-, and I-. The reactive force field parameters have been trained against quantum mechanical (QM) calculations related to water binding energies, hydration energies and energies of proton transfer. The new force field has been validated by applying it to molecular dynamics (MD) simulations of the ionization of different electrolytes in water and comparison of the results with experimental observations and thermodynamics. Radial distribution functions (RDF) determined for most of the atom pairs (cation or anion with oxygen and hydrogen of water) show a good agreement with the RDF values obtained from DFT calculations. On the basis of the applied force field, the ReaxFF simulations have described the diffusion constants for water and electrolyte ions in alkali metal hydroxide and chloride salt solutions as a function of composition and electrolyte concentration. The obtained results open opportunities to advance ReaxFF methodology to a wide range of applications involving electrolyte ions and solutions.
Collapse
Affiliation(s)
- Mark V Fedkin
- Department of Mechanical and Nuclear Engineering , The Pennsylvania State University , University Park , Pennsylvania 16802 , United States
| | - Yun Kyung Shin
- Department of Mechanical and Nuclear Engineering , The Pennsylvania State University , University Park , Pennsylvania 16802 , United States
| | - Nabankur Dasgupta
- Department of Engineering Science and Mechanics , The Pennsylvania State University , University Park , Pennsylvania 16802 , United States
| | - Jejoon Yeon
- Department of Mechanical and Nuclear Engineering , The Pennsylvania State University , University Park , Pennsylvania 16802 , United States.,Center for Composite Materials , University of Delaware , Newark , Delaware 19716 , United States
| | - Weiwei Zhang
- Department of Mechanical and Nuclear Engineering , The Pennsylvania State University , University Park , Pennsylvania 16802 , United States
| | - Diana van Duin
- Department of Mechanical and Nuclear Engineering , The Pennsylvania State University , University Park , Pennsylvania 16802 , United States
| | - Adri C T van Duin
- Department of Mechanical and Nuclear Engineering , The Pennsylvania State University , University Park , Pennsylvania 16802 , United States
| | | | - Atsushi Fujiwara
- Materials Science Department , MOLSIS Inc. , 3-19-9, Hatchobori , Chuo-ku, Tokyo 104-0032 , Japan
| | | | | | - Masahiko Okumura
- Center for Computational Science & e-Systems , Japan Atomic Energy Agency , 178-4-4 Wakashiba , Kashiwa , Chiba 277-0871 , Japan
| |
Collapse
|
12
|
Réal F, Gomes ASP, Guerrero Martínez YO, Ayed T, Galland N, Masella M, Vallet V. Structural, dynamical, and transport properties of the hydrated halides: How do At(-) bulk properties compare with those of the other halides, from F(-) to I(-)? J Chem Phys 2016; 144:124513. [PMID: 27036467 DOI: 10.1063/1.4944613] [Citation(s) in RCA: 24] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022] Open
Abstract
The properties of halides from the lightest, fluoride (F(-)), to the heaviest, astatide (At(-)), have been studied in water using a polarizable force-field approach based on molecular dynamics (MD) simulations at the 10 ns scale. The selected force-field explicitly treats the cooperativity within the halide-water hydrogen bond networks. The force-field parameters have been adjusted to ab initio data on anion/water clusters computed at the relativistic Möller-Plesset second-order perturbation theory level of theory. The anion static polarizabilities of the two heaviest halides, I(-) and At(-), were computed in the gas phase using large and diffuse atomic basis sets, and taking into account both electron correlation and spin-orbit coupling within a four-component framework. Our MD simulation results show the solvation properties of I(-) and At(-) in aqueous phase to be very close. For instance, their first hydration shells are structured and encompass 9.2 and 9.1 water molecules at about 3.70 ± 0.05 Å, respectively. These values have to be compared to the F(-), Cl(-), and Br(-) ones, i.e., 6.3, 8.4, and 9.0 water molecules at 2.74, 3.38, and 3.55 Å, respectively. Moreover our computations predict the solvation free energy of At(-) in liquid water at ambient conditions to be 68 kcal mol(-1), a value also close the I(-) one, about 70 kcal mol(-1). In all, our simulation results for I(-) are in excellent agreement with the latest neutron- and X-ray diffraction studies. Those for the At(-) ion are predictive, as no theoretical or experimental data are available to date.
Collapse
Affiliation(s)
- Florent Réal
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | - André Severo Pereira Gomes
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | | | - Tahra Ayed
- CEISAM UMR CNRS 6230, Université de Nantes, 2 Rue de la Houssinière, BP 92208 F-44322 Nantes Cedex 3, France
| | - Nicolas Galland
- CEISAM UMR CNRS 6230, Université de Nantes, 2 Rue de la Houssinière, BP 92208 F-44322 Nantes Cedex 3, France
| | - Michel Masella
- Laboratoire de Biologie Structurale et Radiobiologie, Service de Bioénergétique, Biologie Structurale et Mécanismes, Institut de Biologie et de Technologies de Saclay, CEA Saclay, F-91191 Gif sur Yvette Cedex, France
| | - Valérie Vallet
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| |
Collapse
|
13
|
Antalek M, Pace E, Hedman B, Hodgson KO, Chillemi G, Benfatto M, Sarangi R, Frank P. Solvation structure of the halides from x-ray absorption spectroscopy. J Chem Phys 2016; 145:044318. [PMID: 27475372 PMCID: PMC4967075 DOI: 10.1063/1.4959589] [Citation(s) in RCA: 29] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2016] [Accepted: 07/11/2016] [Indexed: 11/14/2022] Open
Abstract
Three-dimensional models for the aqueous solvation structures of chloride, bromide, and iodide are reported. K-edge extended X-ray absorption fine structure (EXAFS) and Minuit X-ray absorption near edge (MXAN) analyses found well-defined single shell solvation spheres for bromide and iodide. However, dissolved chloride proved structurally distinct, with two solvation shells needed to explain its strikingly different X-ray absorption near edge structure (XANES) spectrum. Final solvation models were as follows: iodide, 8 water molecules at 3.60 ± 0.13 Å and bromide, 8 water molecules at 3.40 ± 0.14 Å, while chloride solvation included 7 water molecules at 3.15 ± 0.10 Å, and a second shell of 7 water molecules at 4.14 ± 0.30 Å. Each of the three derived solvation shells is approximately uniformly disposed about the halides, with no global asymmetry. Time-dependent density functional theory calculations simulating the chloride XANES spectra following from alternative solvation spheres revealed surprising sensitivity of the electronic state to 6-, 7-, or 8-coordination, implying a strongly bounded phase space for the correct structure during an MXAN fit. MXAN analysis further showed that the asymmetric solvation predicted from molecular dynamics simulations using halide polarization can play no significant part in bulk solvation. Classical molecular dynamics used to explore chloride solvation found a 7-water solvation shell at 3.12 (-0.04/+0.3) Å, supporting the experimental result. These experiments provide the first fully three-dimensional structures presenting to atomic resolution the aqueous solvation spheres of the larger halide ions.
Collapse
Affiliation(s)
- Matthew Antalek
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Stanford University, Menlo Park, California 94025, USA
| | - Elisabetta Pace
- Laboratori Nazionali di Frascati-INFN, P.O. Box 13, 00044 Frascati, Italy
| | - Britt Hedman
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Stanford University, Menlo Park, California 94025, USA
| | - Keith O Hodgson
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Giovanni Chillemi
- CINECA, SCAI-SuperComputing Applications and Innovation Department, Via dei Tizii 6, 00185 Roma, Italy
| | - Maurizio Benfatto
- Laboratori Nazionali di Frascati-INFN, P.O. Box 13, 00044 Frascati, Italy
| | - Ritimukta Sarangi
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Stanford University, Menlo Park, California 94025, USA
| | - Patrick Frank
- Stanford Synchrotron Radiation Lightsource, SLAC National Accelerator Laboratory, Stanford University, Menlo Park, California 94025, USA
| |
Collapse
|
14
|
Sun L, Li X, Tu Y, Ågren H. Origin of ion selectivity at the air/water interface. Phys Chem Chem Phys 2016; 17:4311-8. [PMID: 25574599 DOI: 10.1039/c4cp03338h] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
Among many characteristics of ions, their capability to accumulate at air/water interfaces is a particular issue that has been the subject of much research attention. For example, the accumulation of halide anions (Cl(-), Br(-), I(-)) at the water surface is of great importance to heterogeneous reactions that are of environmental concern. However, the actual mechanism that drives anions towards the air/water interface remains unclear. In this work, we have performed atomistic simulations using polarizable models to mimic ionic behavior under atmospheric conditions. We find that larger anions are abundant at the water surface and that the cations are pulled closer to the surface by the counterions. We propose that polarization effects stabilize the anions with large radii when approaching the surface. This energetically more favorable situation is caused by the fact that the more polarized anions at the surface attract water molecules more strongly. Of relevance is also the ordering of the surface water molecules with their hydrogen atoms pointing outwards which induce an external electronic field that leads to a different surface behavior of anions and cations. The water-water interaction is weakened by the distinct water-ion attraction, a point contradicting the proposition that F(-) is a kosmotrope. The simulation results thus allow us to obtain a more holistic understanding of the interfacial properties of ionic solutions and atmospheric aerosols.
Collapse
Affiliation(s)
- Lu Sun
- Department of Theoretical Chemistry and Biology, School of Biotechnology, Royal Institute of Technology, S-10691 Stockholm, Sweden.
| | | | | | | |
Collapse
|
15
|
Karmakar A, Chandra A. Water in Hydration Shell of an Iodide Ion: Structure and Dynamics of Solute-Water Hydrogen Bonds and Vibrational Spectral Diffusion from First-Principles Simulations. J Phys Chem B 2015; 119:8561-72. [DOI: 10.1021/jp510714e] [Citation(s) in RCA: 30] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Anwesa Karmakar
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, India 208016
| | - Amalendu Chandra
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur, India 208016
| |
Collapse
|
16
|
Karmakar A, Chandra A. Ab initio molecular dynamics studies of hydrogen bonded structure, molecular motion, and frequency fluctuations of water in the vicinity of azide ions. J Chem Phys 2015; 142:164505. [DOI: 10.1063/1.4918579] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Anwesa Karmakar
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| | - Amalendu Chandra
- Department of Chemistry, Indian Institute of Technology Kanpur, Kanpur 208016, India
| |
Collapse
|
17
|
Karmakar A, Chandra A. Effects of dispersion interaction on vibrational spectral diffusion in aqueous NaBr solutions: An ab initio molecular dynamics study. Chem Phys 2015. [DOI: 10.1016/j.chemphys.2014.11.021] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/24/2022]
|
18
|
An ab initio molecular dynamics study of the hydrogen bonded structure, dynamics and vibrational spectral diffusion of water in the ion hydration shell of a superoxide ion. Chem Phys 2014. [DOI: 10.1016/j.chemphys.2014.10.001] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
19
|
Soniat M, Rick SW. Charge transfer effects of ions at the liquid water/vapor interface. J Chem Phys 2014; 140:184703. [DOI: 10.1063/1.4874256] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
|
20
|
Smirnov PR. Structural parameters of the nearest surrounding of halide ions in the aqueous electrolyte solutions. RUSS J GEN CHEM+ 2013. [DOI: 10.1134/s107036321308001x] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
21
|
Badu S, Truflandier L, Autschbach J. Quadrupolar NMR Spin Relaxation Calculated Using Ab Initio Molecular Dynamics: Group 1 and Group 17 Ions in Aqueous Solution. J Chem Theory Comput 2013; 9:4074-86. [PMID: 26592401 DOI: 10.1021/ct400419s] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Electric field gradient (EFG) fluctuations for the monoatomic ions (7)Li(+), (23)Na(+), (35)Cl(-), (81)Br(-), and (127)I(-) in aqueous solution are studied using Car-Parrinello ab initio molecular dynamics (aiMD) simulations based on density functional theory. EFG calculations are typically performed with 1024 ion-solvent configurations from the aiMD simulation, using the Zeroth Order Regular Approximation (ZORA) relativistic Hamiltonian. Autocorrelation functions for the spherical EFG tensor elements are computed, transformed into the corresponding spectral densities (under the extreme narrowing condition), and subsequently converted into NMR quadrupolar relaxation rates for the ions. The relaxation rates are compared with experimental data. The order of magnitude is correctly predicted by the simulations. The computational protocol is tested in detail for (81)Br(-).
Collapse
Affiliation(s)
- Shyam Badu
- Department of Chemistry, University at Buffalo-State University of New York , Buffalo, New York 14260-3000, United States
| | - Lionel Truflandier
- Institut des Sciences Moléculaires, Université Bordeaux I , 351 Cours de la Libration, 33405 Talence, France
| | - Jochen Autschbach
- Department of Chemistry, University at Buffalo-State University of New York , Buffalo, New York 14260-3000, United States
| |
Collapse
|
22
|
Allolio C, Salas-Illanes N, Desmukh YS, Hansen MR, Sebastiani D. H-Bonding Competition and Clustering in Aqueous LiI. J Phys Chem B 2013; 117:9939-46. [DOI: 10.1021/jp4033468] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/16/2023]
Affiliation(s)
- Christoph Allolio
- Department of Chemistry, Martin-Luther Universität Halle-Wittenberg, von-Danckelmann-Platz
4, 06120 Halle/Saale, Germany
| | - Nora Salas-Illanes
- Department of Chemistry, Martin-Luther Universität Halle-Wittenberg, von-Danckelmann-Platz
4, 06120 Halle/Saale, Germany
| | - Yogesh S. Desmukh
- Department of Chemical
Engineering and Chemistry, Eindhoven University of Technology, P.O. Box 513, Eindhoven, Netherlands
- Dutch Polymer Institute (DPI), P.O. Box 902, 5600AX Eindhoven, The Netherlands
| | - Michael Ryan Hansen
- Dutch Polymer Institute (DPI), P.O. Box 902, 5600AX Eindhoven, The Netherlands
- Max Planck Institute for Polymer Research, Ackermannweg 10, 55128 Mainz,
Germany
| | - Daniel Sebastiani
- Department of Chemistry, Martin-Luther Universität Halle-Wittenberg, von-Danckelmann-Platz
4, 06120 Halle/Saale, Germany
| |
Collapse
|
23
|
Ahmed M, Singh AK, Mondal JA, Sarkar SK. Water in the Hydration Shell of Halide Ions Has Significantly Reduced Fermi Resonance and Moderately Enhanced Raman Cross Section in the OH Stretch Regions. J Phys Chem B 2013; 117:9728-33. [DOI: 10.1021/jp403618x] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Affiliation(s)
- Mohammed Ahmed
- Radiation & Photochemistry Division, Bhabha Atomic Research Centre, Trombay, Mumbai-400085, India
| | - Ajay K. Singh
- Radiation & Photochemistry Division, Bhabha Atomic Research Centre, Trombay, Mumbai-400085, India
| | - Jahur A. Mondal
- Radiation & Photochemistry Division, Bhabha Atomic Research Centre, Trombay, Mumbai-400085, India
| | - Sisir K. Sarkar
- Radiation & Photochemistry Division, Bhabha Atomic Research Centre, Trombay, Mumbai-400085, India
| |
Collapse
|
24
|
A first principles simulation study of vibrational spectral diffusion in aqueous NaBr solutions: Dynamics of water in ion hydration shells. Chem Phys 2013. [DOI: 10.1016/j.chemphys.2012.11.007] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
|
25
|
Baer MD, Mundy CJ. An ab initio approach to understanding the specific ion effect. Faraday Discuss 2013; 160:89-101; discussion 103-20. [DOI: 10.1039/c2fd20113e] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
26
|
Matsuo T, Arata T, Oda T, Fujiwara S. Difference in hydration structures between F-actin and myosin subfragment-1 detected by small-angle X-ray and neutron scattering. Biophysics (Nagoya-shi) 2013; 9:99-106. [PMID: 27493547 PMCID: PMC4629667 DOI: 10.2142/biophysics.9.99] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/05/2013] [Accepted: 06/13/2013] [Indexed: 12/01/2022] Open
Abstract
Hydration structures around F-actin and myosin subfragment-1 (S1), which play central roles as counterparts in muscle contraction, were investigated by small-angle X-ray scattering (SAXS) and small-angle neutron scattering (SANS). The radius of gyration of chymotryptic S1 was evaluated to be 41.3±1.1 Å for SAXS, 40.1±3.0 Å for SANS in H2O, and 37.8±0.8 Å for SANS in D2O, respectively. The values of the cross-sectional radius of gyration of F-actin were 25.4±0.03 Å for SAXS, 23.4±2.4 Å for SANS in H2O, and 22.6 ± 0.6 Å for SANS in D2O, respectively. These differences arise from different contributions of the hydration shell to the scattering curves. Analysis by model calculations showed that the hydration shell of S1 has the average density 10–15% higher than bulk water, being the typical hydration shell. On the other hand, the hydration shell of F-actin has the average density more than 19% higher than bulk water, indicating that F-actin has a denser, unusual hydration structure. The results indicate a difference in the hydration structures around F-actin and S1. The unusual hydration structure around F-actin may have the structural property of so-called “hyper-mobile water” around F-actin.
Collapse
Affiliation(s)
- Tatsuhito Matsuo
- Quantum Beam Science Directorate, Japan Atomic Energy Agency, Tokai-mura, Naka-gun, Ibaraki 319-1195, Japan
| | - Toshiaki Arata
- Department of Biological Sciences, Graduate School of Science, Osaka University, Toyonaka, Osaka 560-0043, Japan
| | - Toshiro Oda
- RIKEN SPring-8 Center, RIKEN Harima Institute, Sayo-gun, Hyogo 679-5148, Japan
| | - Satoru Fujiwara
- Quantum Beam Science Directorate, Japan Atomic Energy Agency, Tokai-mura, Naka-gun, Ibaraki 319-1195, Japan
| |
Collapse
|
27
|
Zidi ZS. On the stability of ion water clusters at atmospheric conditions: Open system Monte Carlo simulation. J Chem Phys 2012; 137:124107. [DOI: 10.1063/1.4754528] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
|
28
|
Piatkowski L, Bakker HJ. Vibrational dynamics of the bending mode of water interacting with ions. J Chem Phys 2012; 135:214509. [PMID: 22149804 DOI: 10.1063/1.3664866] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
We studied the vibrational relaxation dynamics of the bending mode (ν(2)) of the H(2)O water molecules in the presence of different salts (LiCl, LiBr, LiI, NaI, CsI, NaClO(4), and NaBF(4)). The linear and nonlinear spectra of the bending mode show distinct responses of water molecules hydrating the anions. We observe that the bending mode of water molecules that are hydrogen-bonded to an anion exhibits much slower relaxation rates (T(1)~1ps) than water molecules that are hydrogen-bonded to other water molecules (T(1)=400 fs). We find that the effect of the anion on the absorption spectrum and relaxation time constant of the water bending mode is not only determined by the strength of the hydrogen-bond interaction but also by the shape of the anion.
Collapse
Affiliation(s)
- L Piatkowski
- FOM Institute AMOLF, Science Park 104, 1098 XG Amsterdam, The Netherlands.
| | | |
Collapse
|
29
|
Trumm M, Martínez YOG, Réal F, Masella M, Vallet V, Schimmelpfennig B. Modeling the hydration of mono-atomic anions from the gas phase to the bulk phase: The case of the halide ions F−, Cl−, and Br−. J Chem Phys 2012; 136:044509. [DOI: 10.1063/1.3678294] [Citation(s) in RCA: 34] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
|
30
|
Boisson J, Stirnemann G, Laage D, Hynes JT. Water reorientation dynamics in the first hydration shells of F- and I-. Phys Chem Chem Phys 2011; 13:19895-901. [PMID: 21915404 DOI: 10.1039/c1cp21834d] [Citation(s) in RCA: 50] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Molecular dynamics and analytic theory results are presented for the reorientation dynamics of first hydration shell water molecules around fluoride and iodide anions. These ions represent the extremes of the (normal) halide series in terms of their size and conventional structure-making and -breaking categorizations. The simulated reorientation times are consistent with NMR and ultrafast IR experimental results. They are also in good agreement with the theoretical predictions of the analytic Extended Jump Model. Analysis through this model shows that while sudden, large amplitude jumps (in which the reorienting water exchanges hydrogen-bond partners) are the dominant reorientation pathway for the I(-) case, they are comparatively less important for the F(-) case. In particular, the diffusive reorientation of an intact F(-)···H(2)O hydrogen-bonded pair is found to be most important for the reorientation time, a feature related to the greater hydrogen-bond strength for the F(-)···H(2)O pair. The dominance of this effect for e.g. multiply charged ions is suggested.
Collapse
Affiliation(s)
- Jean Boisson
- Ecole Normale Supérieure, Chemistry Department, UMR CNRS-ENS-UPMC 8640, Paris, France
| | | | | | | |
Collapse
|
31
|
Mallik BS, Chandra A. An ab initio molecular dynamics study of supercritical aqueous ionic solutions: Hydrogen bonding, rotational dynamics and vibrational spectral diffusion. Chem Phys 2011. [DOI: 10.1016/j.chemphys.2011.06.031] [Citation(s) in RCA: 12] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
32
|
Heisler IA, Mazur K, Meech SR. Low-Frequency Modes of Aqueous Alkali Halide Solutions: An Ultrafast Optical Kerr Effect Study. J Phys Chem B 2011; 115:1863-73. [DOI: 10.1021/jp111239v] [Citation(s) in RCA: 54] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ismael A. Heisler
- School of Chemistry, University of East Anglia, Norwich NR4 7TJ, U.K
| | - Kamila Mazur
- School of Chemistry, University of East Anglia, Norwich NR4 7TJ, U.K
| | - Stephen R. Meech
- School of Chemistry, University of East Anglia, Norwich NR4 7TJ, U.K
| |
Collapse
|
33
|
Kirchner B, di Dio PJ, Hutter J. Real-world predictions from ab initio molecular dynamics simulations. Top Curr Chem (Cham) 2011; 307:109-53. [PMID: 21842358 DOI: 10.1007/128_2011_195] [Citation(s) in RCA: 66] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
In this review we present the techniques of ab initio molecular dynamics simulation improved to its current stage where the analysis of existing processes and the prediction of further chemical features and real-world processes are feasible. For this reason we describe the relevant developments in ab initio molecular dynamics leading to this stage. Among them, parallel implementations, different basis set functions, density functionals, and van der Waals corrections are reported. The chemical features accessible through AIMD are discussed. These are IR, NMR, as well as EXAFS spectra, sampling methods like metadynamics and others, Wannier functions, dipole moments of molecules in condensed phase, and many other properties. Electrochemical reactions investigated by ab initio molecular dynamics methods in solution, on surfaces as well as complex interfaces, are also presented.
Collapse
Affiliation(s)
- Barbara Kirchner
- Wilhelm-Ostwald Institute of Physical and Theoretical Chemistry, University of Leipzig, Linnéstr. 2, 04103 Leipzig, Germany.
| | | | | |
Collapse
|
34
|
Fulton JL, Schenter GK, Baer MD, Mundy CJ, Dang LX, Balasubramanian M. Probing the Hydration Structure of Polarizable Halides: A Multiedge XAFS and Molecular Dynamics Study of the Iodide Anion. J Phys Chem B 2010; 114:12926-37. [DOI: 10.1021/jp106378p] [Citation(s) in RCA: 66] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- John L. Fulton
- Chemical and Materials Sciences Division, Pacific Northwest National Laboratory, Richland, Washington, 99354, Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany, and Advanced Photon Source, Argonne National Laboratory, Argonne, Illinois 60439
| | - Gregory K. Schenter
- Chemical and Materials Sciences Division, Pacific Northwest National Laboratory, Richland, Washington, 99354, Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany, and Advanced Photon Source, Argonne National Laboratory, Argonne, Illinois 60439
| | - Marcel D. Baer
- Chemical and Materials Sciences Division, Pacific Northwest National Laboratory, Richland, Washington, 99354, Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany, and Advanced Photon Source, Argonne National Laboratory, Argonne, Illinois 60439
| | - Christopher J. Mundy
- Chemical and Materials Sciences Division, Pacific Northwest National Laboratory, Richland, Washington, 99354, Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany, and Advanced Photon Source, Argonne National Laboratory, Argonne, Illinois 60439
| | - Liem X. Dang
- Chemical and Materials Sciences Division, Pacific Northwest National Laboratory, Richland, Washington, 99354, Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany, and Advanced Photon Source, Argonne National Laboratory, Argonne, Illinois 60439
| | - Mahalingam Balasubramanian
- Chemical and Materials Sciences Division, Pacific Northwest National Laboratory, Richland, Washington, 99354, Lehrstuhl für Theoretische Chemie, Ruhr-Universität Bochum, D-44780 Bochum, Germany, and Advanced Photon Source, Argonne National Laboratory, Argonne, Illinois 60439
| |
Collapse
|
35
|
Terrier C, Vitorge P, Gaigeot MP, Spezia R, Vuilleumier R. Density functional theory based molecular dynamics study of hydration and electronic properties of aqueous La3+. J Chem Phys 2010; 133:044509. [DOI: 10.1063/1.3460813] [Citation(s) in RCA: 33] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
36
|
The solvent shell structure of aqueous iodide: X-ray absorption spectroscopy and classical, hybrid QM/MM and full quantum molecular dynamics simulations. Chem Phys 2010. [DOI: 10.1016/j.chemphys.2010.03.023] [Citation(s) in RCA: 46] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022]
|
37
|
D’Angelo P, Migliorati V, Guidoni L. Hydration Properties of the Bromide Aqua Ion: the Interplay of First Principle and Classical Molecular Dynamics, and X-ray Absorption Spectroscopy. Inorg Chem 2010; 49:4224-31. [DOI: 10.1021/ic9025574] [Citation(s) in RCA: 72] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Affiliation(s)
- Paola D’Angelo
- Dipartimento di Chimica, Università di Roma “La Sapienza”, P.le A. Moro 5, 00185 Roma, Italy
| | - Valentina Migliorati
- Dipartimento di Chimica, Università di Roma “La Sapienza”, P.le A. Moro 5, 00185 Roma, Italy
| | - Leonardo Guidoni
- Dipartimento di Chimica, Ingegneria Chimica e Materiali, Università degli Studi dell’Aquila, via Campo di Pile, zona industriale di Pile, 67100, L’Aquila, Italy
| |
Collapse
|
38
|
|
39
|
Tongraar A, Hannongbua S, Rode BM. QM/MM MD Simulations of Iodide Ion (I−) in Aqueous Solution: A Delicate Balance between Ion−Water and Water−Water H-Bond Interactions. J Phys Chem A 2010; 114:4334-9. [DOI: 10.1021/jp910435d] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Anan Tongraar
- School of Chemistry, Institute of Science, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand, and Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 52a, A-6020 Innsbruck, Austria
| | - Supot Hannongbua
- School of Chemistry, Institute of Science, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand, and Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 52a, A-6020 Innsbruck, Austria
| | - Bernd Michael Rode
- School of Chemistry, Institute of Science, Suranaree University of Technology, Nakhon Ratchasima 30000, Thailand, Department of Chemistry, Faculty of Science, Chulalongkorn University, Bangkok 10330, Thailand, and Theoretical Chemistry Division, Institute of General, Inorganic and Theoretical Chemistry, University of Innsbruck, Innrain 52a, A-6020 Innsbruck, Austria
| |
Collapse
|
40
|
Carnevale V, Raugei S. Structural aspects of the solvation shell of lysine and acetylated lysine: A Car–Parrinello and classical molecular dynamics investigation. J Chem Phys 2009; 131:225103. [DOI: 10.1063/1.3268703] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
|
41
|
Osuna SÃ, Swart M, Baerends EJ, Bickelhaupt FM, Solà M. Homolytic versus Heterolytic Dissociation of Alkalimetal Halides: The Effect of Microsolvation. Chemphyschem 2009; 10:2955-65. [DOI: 10.1002/cphc.200900480] [Citation(s) in RCA: 13] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
42
|
Lin YS, Auer BM, Skinner JL. Water structure, dynamics, and vibrational spectroscopy in sodium bromide solutions. J Chem Phys 2009; 131:144511. [DOI: 10.1063/1.3242083] [Citation(s) in RCA: 126] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/26/2023] Open
|
43
|
Guàrdia E, Skarmoutsos I, Masia M. On Ion and Molecular Polarization of Halides in Water. J Chem Theory Comput 2009; 5:1449-53. [DOI: 10.1021/ct900096n] [Citation(s) in RCA: 70] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Affiliation(s)
- Elvira Guàrdia
- Departament de Física i Enginyeria Nuclear, Universitat Politècnica de Catalunya, Campus Nord B4-B5, Barcelona 08034, Spain
| | - Ioannis Skarmoutsos
- Departament de Física i Enginyeria Nuclear, Universitat Politècnica de Catalunya, Campus Nord B4-B5, Barcelona 08034, Spain
| | - Marco Masia
- Dipartimento di Chimica, Università degli Studi di Sassari, Sardinian Laboratory for Computational Materials Science SLACS (INFM-CNR) and INSTM, Via Vienna 2, 07100 Sassari, Italy
| |
Collapse
|
44
|
Perera PN, Browder B, Ben-Amotz D. Perturbations of Water by Alkali Halide Ions Measured using Multivariate Raman Curve Resolution. J Phys Chem B 2009; 113:1805-9. [PMID: 19199691 DOI: 10.1021/jp808732s] [Citation(s) in RCA: 65] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
|
45
|
Mallik BS, Semparithi A, Chandra A. A first principles theoretical study of vibrational spectral diffusion and hydrogen bond dynamics in aqueous ionic solutions: D2O in hydration shells of Cl(-) ions. J Chem Phys 2009; 129:194512. [PMID: 19026071 DOI: 10.1063/1.3006032] [Citation(s) in RCA: 98] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A theoretical study of vibrational spectral diffusion and hydrogen bond dynamics in aqueous ionic solutions is presented from first principles without employing any empirical potential models. The present calculations are based on ab initio molecular dynamics for trajectory generation and wavelet analysis of the simulated trajectories for time dependent frequency calculations. Results are obtained for two different deuterated aqueous solutions: the first one is a relatively dilute solution of a single Cl(-) ion and the second one is a concentrated solution of NaCl ( approximately 3M) dissolved in liquid D(2)O. It is found that the frequencies of OD bonds in the anion hydration shell, i.e., those which are hydrogen bonded to the chloride ion, have a higher stretch frequency than those in the bulk water. Also, on average, the frequencies of hydration shell OD modes are found to increase with increase in the anion-water hydrogen bond distance. On the dynamical side, when the vibrational spectral diffusion is calculated exclusively for the hydration shell water molecules in the first solution, the dynamics reveals three time scales: a short-time relaxation ( approximately 200 fs) corresponding to the dynamics of intact ion-water hydrogen bonds, a slower relaxation ( approximately 3 ps) corresponding to the lifetimes of chloride ion-water hydrogen bonds, and another longer-time constant ( approximately 20 ps) corresponding to the escape dynamics of water from the anion hydration shell. Existence of such three time scales for hydration shell water molecules was also reported earlier for water containing a single iodide ion using classical molecular dynamics [B. Nigro et al., J. Phys. Chem. A 110, 11237 (2006)]. Hence, the present study confirms the basic results of this earlier work using a different methodology. However, when the vibrational spectral diffusion is calculated over all the OD modes, only two time scales of approximately 150 fs and approximately 2.7 ps are found without the slowest component of approximately 20 ps. This is likely because of the very small weight that the hydration shell water molecules carry to the overall spectral diffusion in the solution containing a single ion. For the concentrated solution also, the slowest component of approximately 20 ps is not found in the spectral diffusion of all water molecules because a distinct separation between the hydration shell and bulk water in terms of their stretch frequencies does not hold at this high concentration regime. The present first principles results are compared with those of the available experiments and classical simulations.
Collapse
Affiliation(s)
- Bhabani S Mallik
- Department of Chemistry, Indian Institute of Technology, Kanpur 208016, India
| | | | | |
Collapse
|
46
|
Shevkunov SV. Structural transition in the OH−(H2O) n cluster in water vapors. COLLOID JOURNAL 2008. [DOI: 10.1134/s1061933x08060161] [Citation(s) in RCA: 11] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
47
|
Glover WJ, Larsen RE, Schwartz BJ. The roles of electronic exchange and correlation in charge-transfer-to-solvent dynamics: Many-electron nonadiabatic mixed quantum/classical simulations of photoexcited sodium anions in the condensed phase. J Chem Phys 2008; 129:164505. [DOI: 10.1063/1.2996350] [Citation(s) in RCA: 23] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
48
|
Affiliation(s)
- H. J. Bakker
- FOM Institute for Atomic and Molecular Physics, Kruislaan 407, 1098 SJ Amsterdam, The Netherlands
| |
Collapse
|
49
|
|
50
|
Nigro B, Re S, Laage D, Rey R, Hynes JT. On the ultrafast infrared spectroscopy of anion hydration shell hydrogen bond dynamics. J Phys Chem A 2007; 110:11237-43. [PMID: 17004732 DOI: 10.1021/jp064846m] [Citation(s) in RCA: 62] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Molecular Dynamics simulations are used to examine the title issue for the I-/HOD/D2O solution system in connection with recent ultrafast infrared spectroscopic experiments. It is argued that the long "modulation time" associated with the spectral diffusion of the OH frequency, extracted in these experiments, should be interpreted as reflecting the escape time of an HOD from the first hydration shell of the I- ion, i.e., the residence time of an HOD in this solvation shell. Shorter time features related to the oscillation of the OH ...I- hydrogen bond and the breaking and making of this bond are also discussed.
Collapse
Affiliation(s)
- Bruno Nigro
- Département de Chimie, UMR 8640 PASTEUR, Ecole Normale Supérieure, 24 rue Lhomond, 75231 Paris, France.
| | | | | | | | | |
Collapse
|