1
|
Wang C, Tian W, Zhou K. Ab Initio Simulation of Liquid Water without Artificial High Temperature. J Chem Theory Comput 2024. [PMID: 39219067 DOI: 10.1021/acs.jctc.4c00650] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/04/2024]
Abstract
Comprehending the structure and dynamics of water is crucial in various fields, such as water desalination, ion separation, electrocatalysis, and biochemical processes. While reported works show that the ab initio molecular dynamics (AIMD) can accurately portray water's structure, the artificial high temperature (AHT) from 120 to 30 K is needed to mimic the quantum nature of hydrogen-bond network from GGA, metaGGA to hybrid functionals. The AHT proves to be an inadequate approach for systems involving aqueous multiphase mixtures, such as water-solid interfaces and aqueous solutions. This is due to the activation of additional phonons in other phases, which can lead to an overestimation of the dynamics of nearby water molecules. In this work, we find that the regularized SCAN (rSCAN) functional effectively captures both the structure and dynamics of liquid water at ambient conditions without AHT. Moreover, rSCAN closely matches experimental results for the hydration structures of alkali, alkali earth, and halide ions. We anticipate that the versatile and accurate rSCAN functional will emerge as a key tool based on ab initio simulation for investigating chemical processes in aqueous environments.
Collapse
Affiliation(s)
- Chenyu Wang
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou 215006, China
| | - Wei Tian
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou 215006, China
| | - Ke Zhou
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou 215006, China
| |
Collapse
|
2
|
Knight B, Mondal R, Han N, Pietra NF, Hall BA, Edgar KJ, Vaissier Welborn V, Madsen LA, De Yoreo JJ, Dove PM. Kinetics of Calcite Nucleation onto Sulfated Chitosan Derivatives and Implications for Water-Polysaccharide Interactions during Crystallization of Sparingly Soluble Salts. CRYSTAL GROWTH & DESIGN 2024; 24:6338-6353. [PMID: 39131446 PMCID: PMC11311137 DOI: 10.1021/acs.cgd.4c00602] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/01/2024] [Revised: 07/01/2024] [Accepted: 07/02/2024] [Indexed: 08/13/2024]
Abstract
Anionic macromolecules are found at sites of CaCO3 biomineralization in diverse organisms, but their roles in crystallization are not well-understood. We prepared a series of sulfated chitosan derivatives with varied positions and degrees of sulfation, DS(SO3 -), and measured calcite nucleation rate onto these materials. Fitting the classical nucleation theory model to the kinetic data reveals the interfacial free energy of the calcite-polysaccharide-solution system, γnet, is lowest for nonsulfated controls and increases with DS(SO3 -). The kinetic prefactor also increases with DS(SO3 -). Simulations of Ca2+-H2O-chitosan systems show greater water structuring around sulfate groups compared to uncharged substituents, independent of sulfate location. Ca2+-SO3 - interactions are solvent-separated by distances that are inversely correlated with DS(SO3 -) of the polysaccharide. The simulations also predict SO3 - and NH3 + groups affect the solvation waters and HCO3 - ions associated with Ca2+. Integrating the experimental and computational evidence suggests sulfate groups influence nucleation by increasing the difficulty of displacing near-surface water, thereby increasing γnet. By correlating γnet and net charge per monosaccharide for diverse polysaccharides, we suggest the solvent-separated interactions of functional groups with Ca2+ influence thermodynamic and kinetic components to crystallization by similar solvent-dominated processes. The findings reiterate the importance of establishing water structure and properties at macromolecule-solution interfaces.
Collapse
Affiliation(s)
- Brenna
M. Knight
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department
of Geosciences, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Ronnie Mondal
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
- Macromolecules
Innovation Institute, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Nizhou Han
- Department
of Geosciences, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Nicholas F. Pietra
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
- Macromolecules
Innovation Institute, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Brady A. Hall
- GlycoMIP, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Kevin J. Edgar
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
- Macromolecules
Innovation Institute, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department
of Sustainable Biomaterials, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Valerie Vaissier Welborn
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
- Macromolecules
Innovation Institute, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Louis A. Madsen
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
- Macromolecules
Innovation Institute, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - James J. De Yoreo
- Physical
Sciences Division, Physical and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99352, United States
- Department
of Materials Science and Engineering, University
of Washington, Seattle, Washington 98195, United States
| | - Patricia M. Dove
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department
of Geosciences, Virginia Tech, Blacksburg, Virginia 24061, United States
- Macromolecules
Innovation Institute, Virginia Tech, Blacksburg, Virginia 24061, United States
- Department
of Materials Science and Engineering, Virginia
Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
3
|
Tang B, Bai Q, Fang YG, Francisco JS, Zhu C, Fang WH. Mechanistic Insights into N 2O 5-Halide Ions Chemistry at the Air-Water Interface. J Am Chem Soc 2024; 146:21742-21751. [PMID: 39074151 DOI: 10.1021/jacs.4c05850] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/31/2024]
Abstract
The activation of halogens (X = Cl, Br, I) by N2O5 is linked to NOx sources, ozone concentrations, NO3 reactivity, and the chemistry of halide-containing aerosol particles. However, a detailed chemical mechanism is still lacking. Herein, we explored the chemistry of the N2O5···X- systems at the air-water interface. Two different reaction pathways were identified for the reaction of N2O5 with X- at the air-water interface: the formation of XNO2 or XONO, along with NO3-. In the case of the Cl- system, the ClNO2 generation pathway is more favorable, while for the Br- and I- systems, the formation of BrONO and IONO is barrierless, making them the predominant products. Furthermore, the mechanisms of formation of X2 from XNO2 and XONO were also investigated. The high energy barriers of reactions and the high free energies of the products compared to those of the reactants indicate that ClNO2 is stable at the air-water interface. Contrary to the widely held belief regarding X2 producing from the reaction of XNO2 with X-, our calculations demonstrate that BrONO and IONO initially form stable BrONO···Br- and IONO···I- complexes, which then subsequently react with Br- and I- to form Br3- and I3-, respectively. Finally, Br3- and I3- decompose to form Br2 and I2. These findings have significant implications for experimental interpretation and offer new insights into halogen cycling in the atmosphere.
Collapse
Affiliation(s)
- Bo Tang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Qi Bai
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Ye-Guang Fang
- Laboratory of Theoretical and Computational Nanoscience, CAS Key Laboratory of Nanosystem and Hierarchical Fabrication, CAS Center for Excellence in Nanoscience, National Center for Nanoscience and Technology, Beijing 100190, China
| | - Joseph S Francisco
- Department of Chemistry, University of Pennsylvania, Philadelphia, Pennsylvania 19104, United States
| | - Chongqin Zhu
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Wei-Hai Fang
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| |
Collapse
|
4
|
Kastinen T, Batys P, Tolmachev D, Laasonen K, Sammalkorpi M. Ion-Specific Effects on Ion and Polyelectrolyte Solvation. Chemphyschem 2024; 25:e202400244. [PMID: 38712639 DOI: 10.1002/cphc.202400244] [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: 03/05/2024] [Revised: 05/02/2024] [Accepted: 05/02/2024] [Indexed: 05/08/2024]
Abstract
Ion-specific effects on aqueous solvation of monovalent counter ions, Na+ ${^+ }$ , K+ ${^+ }$ , Cl- ${^- }$ , and Br- ${^- }$ , and two model polyelectrolytes (PEs), poly(styrene sulfonate) (PSS) and poly(diallyldimethylammonium) (PDADMA) were here studied with ab initio molecular dynamics (AIMD) and classical molecular dynamics (MD) simulations based on the OPLS-aa force-field which is an empirical fixed point-charge force-field. Ion-specific binding to the PE charge groups was also characterized. Both computational methods predict similar response for the solvation of the PEs but differ notably in description of ion solvation. Notably, AIMD captures the experimentally observed differences in Cl- ${^- }$ and Br- ${^- }$ anion solvation and binding with the PEs, while the classical MD simulations fail to differentiate the ion species response. Furthermore, the findings show that combining AIMD with the computationally less costly classical MD simulations allows benefiting from both the increased accuracy and statistics reach.
Collapse
Affiliation(s)
- Tuuva Kastinen
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076, Aalto, Finland
- Academy of Finland Center of Excellence in Life-Inspired Hybrid Materials (LIBER), Aalto University, P.O. Box 16100, 00076, Aalto, Finland
- Faculty of Engineering and Natural Sciences, Tampere University, P.O. Box 541, 33014, Tampere University, Finland
| | - Piotr Batys
- Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, PL-30239, Krakow, Poland
| | - Dmitry Tolmachev
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076, Aalto, Finland
- Academy of Finland Center of Excellence in Life-Inspired Hybrid Materials (LIBER), Aalto University, P.O. Box 16100, 00076, Aalto, Finland
| | - Kari Laasonen
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076, Aalto, Finland
| | - Maria Sammalkorpi
- Department of Chemistry and Materials Science, Aalto University, P.O. Box 16100, 00076, Aalto, Finland
- Academy of Finland Center of Excellence in Life-Inspired Hybrid Materials (LIBER), Aalto University, P.O. Box 16100, 00076, Aalto, Finland
- Department of Bioproducts and Biosystems, Aalto University, P.O. Box 16100, 00076, Aalto, Finland
| |
Collapse
|
5
|
Awgchew H, Beyene S, Kifilu A. Potassium adsorption capacity and desorption kinetics in soils of Qenberenaweti sub-watershed, central highlands of Ethiopia. Heliyon 2024; 10:e31336. [PMID: 38803861 PMCID: PMC11129098 DOI: 10.1016/j.heliyon.2024.e31336] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2023] [Revised: 05/14/2024] [Accepted: 05/14/2024] [Indexed: 05/29/2024] Open
Abstract
Determining the supply and uptake of K nutrient and the dynamic equilibrium (adsorption-desorption) reactions among the K forms in the soils are not commonly addressed in the highlands of Ethiopia. A study was therefore initiated to determine the adsorption capacity of the exchangeable K and the release kinetics of the non-exchangeable K in the soils of the Qenberenaweti Sub-watershed. Twelve disturbed surface (0-20 cm) soil sub-samples were collected from every farmland which was representative of each pre-identified soil type (Vertic Cambisols, Pellic Vertisols, Pisoplinthic Luvisols, Relictistagnic Cambisols, Pisoplinthic Cambisols, and Plinthofractic Cambisols). A composite sample was made in duplicate for the determination of K adsorption capacity and desorption kinetics per soil type. The mean maximum (69.47 ± 4.31 %) and minimum (56.16 ± 6.04 %) K adsorption rates were obtained from the Plinthofractic Cambisols and Vertic Cambisols, respectively. Among the tested isotherm models, the goodness of the Freundlich was better fit the data of all experimental soils; hence, a modified equation of this model (qe = aCeb/a) could be used to describe the theoretical doses of K fertilizers required to develop K levels in soil solutions. The highest constant K releases from the Plinthofractic Cambisols (47 mg kg-1), Pisiopllintic Cambisols (46 mg kg-1), and Pisoplinthic Luvisols (44 mg kg-1) were attained at the 9th extraction. In comparison, it was noticed at the 7th and 11th extractions of the Relictistagnic Cambisols (45 mg kg-1) and both Pellic Vertisols (48 mg kg-1) and Vertic Cambisols (42 mg kg-1), respectively. The equation of power function was the best to successfully describe the released K+ from all the experimental soils. Eventually, determining the adsorption capacity and release kinetics of K at a site-specific level helps to know the relative potential of the soils to supply K and also plan for an effective K fertilization strategy.
Collapse
Affiliation(s)
- Haymanot Awgchew
- School of Plant and Horticultural Science, Hawassa University, Hawassa, Ethiopia
| | - Sheleme Beyene
- School of Plant and Horticultural Science, Hawassa University, Hawassa, Ethiopia
| | - Alemayehu Kifilu
- School of Plant and Horticultural Science, Hawassa University, Hawassa, Ethiopia
| |
Collapse
|
6
|
Kacenauskaite L, Moncada Cohen M, Van Wyck SJ, Fayer MD. Fast Structural Dynamics in Concentrated HCl Solutions: From Proton Hopping to the Bulk Viscosity. J Am Chem Soc 2024; 146:12355-12364. [PMID: 38682723 DOI: 10.1021/jacs.3c11620] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/01/2024]
Abstract
Concentrated acid solutions, particularly HCl, have been studied extensively to examine the proton hopping and infrared spectral signatures of hydronium ions. Much less attention has been given to the structural dynamics of concentrated HCl solutions. Here, we apply optical heterodyne detected-optical Kerr effect (OHD-OKE) measurements to examine HCl concentration-dependent dynamics from moderate (0.8 m) to very high (15.5 m) concentrations and compare the results to the dynamics of NaCl solutions, as Na+ is similar in size to the hydronium cation. Both HCl and NaCl OHD-OKE signals decay as triexponentials at all concentrations, in contrast to pure water, which decays as a biexponential. Two remarkable features of the HCl dynamics are the following: (1) the bulk viscosity is linearly related to the slowest decay constant, t3, and (2) the concentration-dependent proton hopping times, determined by ab initio MD simulations and 2D IR chemical exchange experiments, both obtained from the literature, fall on the same line as the slowest structural dynamics relaxation time, t3, within experimental error. The structural dynamics of hydronium/chloride/water clusters, with relaxation times t3, are responsible for the concentration dependence of microscopic property of proton hopping and the macroscopic bulk viscosity. The slowest time constant (t3), which does not have a counterpart in pure water, is 3 ps at 0.8 m and increases by a factor of ∼2 by 15.5 m. The two fastest HCl decay constants, t1 and t2, are similar to those of pure water and increase mildly with the concentration.
Collapse
Affiliation(s)
- Laura Kacenauskaite
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
- Nano-Science Center & Department of Chemistry, University of Copenhagen, Copenhagen 2100, Denmark
| | - Max Moncada Cohen
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Stephen J Van Wyck
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| | - Michael D Fayer
- Department of Chemistry, Stanford University, Stanford, California 94305, United States
| |
Collapse
|
7
|
Moison H, Aufort J, Benoit M, Méheut M. On Local Structure Equilibration of Ca 2+ in Solution by Ab Initio Molecular Dynamics. J Phys Chem B 2024; 128:3167-3181. [PMID: 38525554 DOI: 10.1021/acs.jpcb.3c07308] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/26/2024]
Abstract
Analyzing the stable isotopic ratio of Ca offers valuable insights into a wide range of applications from climate reconstruction to bone cancer diagnosis and agricultural nutrient improvement. While the first hydration shell of Ca in solution is expected to play a major role in its fractionation properties, the coordination of Ca in water remains a subject of debate. In this work, Ca2+ in water has been modeled by means of ab initio molecular dynamics simulations using various exchange and correlation functionals and at different temperatures. Results show a significant effect of the functional on the average Ca2+ coordination, depending on its tendency to over- or understructure liquid water. The BLYP functional with Grimme-D2 correction was judged as the most accurate among those tested based on its accuracy to reproduce water structural and diffusion properties. Using this functional, the effect of temperature has been systematically investigated, focusing on means to limit the uncertainty in our assessments of the average coordination of Ca2+ ions by (1) estimating the number of water exchanges in the simulations and (2) implementing a statistical approach based on Markov chains. The findings indicate, especially, that our simulations at 300, 350, and 400 K do not yield converged results due to potential equilibration problems. These observations impose substantial constraints on the trustworthiness of numerous estimates in the existing literature that depend on trajectories with insufficient exchanges. We estimate Ca2+ coordination values of 6.8 ± 0.1, 6.8 ± 0.1, 6.7 ± 0.2, and 6.7 ± 0.2 at 600, 550, 500, and 450 K respectively. At lower temperatures (300, 350, and 400 K), while obtaining definitive values for Ca2+ coordination remains challenging, our research does indicate a potential temperature-related influence on coordination with an average Ca2+ coordination at 300 K as low as 6.2.
Collapse
Affiliation(s)
- Hugo Moison
- GET, OMP, Université Paul Sabatier, 14, avenue Édouard Belin, 31400 Toulouse, France
| | - Julie Aufort
- GET, OMP, Université Paul Sabatier, 14, avenue Édouard Belin, 31400 Toulouse, France
- IMPMC, Sorbonne Université, 4, place Jussieu, 75252 Paris Cedex 5, France
| | - Magali Benoit
- CEMES CNRS and Université Toulouse, 39 rue Jeanne Marvig, 31055 Toulouse, France
| | - Merlin Méheut
- GET, OMP, Université Paul Sabatier, 14, avenue Édouard Belin, 31400 Toulouse, France
| |
Collapse
|
8
|
Giri P, Barath V S, Dhurua S, Maity S, Gazi R, Jana M. Boron-containing fullerene-based salts with cyclic carbonate solvents as electrolytes for Li-ion batteries and beyond. Phys Chem Chem Phys 2024; 26:9317-9328. [PMID: 38444289 DOI: 10.1039/d3cp04706g] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/07/2024]
Abstract
Replacement of carbon atoms by a heteroatom in fullerene is a promising route that enhances the electronic properties of fullerenes and results in hetero fullerene-based effective agents ensuring applications in vivid fields of the solar cell, cathode materials for batteries, etc. Towards the development of new electrolyte salts, attention has been paid to facilitating ion mobility in particular and moderate stability of the anions in addition. From the atomistic molecular dynamics simulation studies, for the first time, we uncover that the boron-containing hetero fullerene, C59B- anion-based LiC59B, and NaC59B salts in cyclic carbonate solvents can act as efficient electrolytes by improving the transport phenomenon of the metal ions in solution, importantly for Li+ and satisfactorily for Na+ as compared to their commonly used BF4- anion based salts. Additionally, our study revealed that apart from LiC59B, and NaC59B salts, C58B22- based MgC58B2 salt can facilitate the ionic conductivity of the electrolyte. The properties of the proposed electrolyte under an electric field and different temperatures were investigated. Some of the bulk properties of the used electrolytes to some extent were found to be improved in the presence of these salts. The first principle-based electrochemical calculations further justify the stability of the proposed anions. The initial investigation from the Reactive force-field (ReaxFF) based atomistic simulations study elucidates that LiC59B reduces the decomposition of the EC solvent compared to LiBF4 and facilitates solvent stability.
Collapse
Affiliation(s)
- Piyusaranjan Giri
- Molecular Simulation Laboratory, Department of Chemistry, National Institute of Technology, Rourkela-769008, India.
| | - Sourab Barath V
- Molecular Simulation Laboratory, Department of Chemistry, National Institute of Technology, Rourkela-769008, India.
| | - Shakuntala Dhurua
- Molecular Simulation Laboratory, Department of Chemistry, National Institute of Technology, Rourkela-769008, India.
| | - Sankar Maity
- Molecular Simulation Laboratory, Department of Chemistry, National Institute of Technology, Rourkela-769008, India.
| | - Rabiul Gazi
- Molecular Simulation Laboratory, Department of Chemistry, National Institute of Technology, Rourkela-769008, India.
| | - Madhurima Jana
- Molecular Simulation Laboratory, Department of Chemistry, National Institute of Technology, Rourkela-769008, India.
| |
Collapse
|
9
|
Palos E, Caruso A, Paesani F. Consistent density functional theory-based description of ion hydration through density-corrected many-body representations. J Chem Phys 2023; 159:181101. [PMID: 37947509 DOI: 10.1063/5.0174577] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2023] [Accepted: 10/23/2023] [Indexed: 11/12/2023] Open
Abstract
Delocalization error constrains the accuracy of density functional theory in describing molecular interactions in ion-water systems. Using Na+ and Cl- in water as model systems, we calculate the effects of delocalization error in the SCAN functional for describing ion-water and water-water interactions in hydrated ions, and demonstrate that density-corrected SCAN (DC-SCAN) predicts n-body and interaction energies with an accuracy approaching coupled cluster theory. The performance of DC-SCAN is size-consistent, maintaining an accurate description of molecular interactions well beyond the first solvation shell. Molecular dynamics simulations at ambient conditions with many-body MB-SCAN(DC) potentials, derived from the many-body expansion, predict the solvation structure of Na+ and Cl- in quantitative agreement with reference data, while simultaneously reproducing the structure of liquid water. Beyond rationalizing the accuracy of density-corrected models of ion hydration, our findings suggest that our unified density-corrected MB formalism holds great promise for efficient DFT-based simulations of condensed-phase systems with chemical accuracy.
Collapse
Affiliation(s)
- Etienne Palos
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Alessandro Caruso
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
| | - Francesco Paesani
- Department of Chemistry and Biochemistry, University of California San Diego, La Jolla, California 92093, USA
- Materials Science and Engineering, University of California San Diego, La Jolla, California 92093, USA
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093, USA
| |
Collapse
|
10
|
Huang M, Jiang Y, Luo Z, Wang J, Ding Z, Guo X, Liu X, Wang Y. Transition metal doped WSi 2N 4monolayer for water splitting electrocatalysts: a first-principles study. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2023; 35:485001. [PMID: 37665141 DOI: 10.1088/1361-648x/acf263] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/14/2023] [Accepted: 08/21/2023] [Indexed: 09/05/2023]
Abstract
High-performance water splitting electrocatalysts are urgently needed in the face of the environmental degradation and energy crisis. The first principles method was used in this study to systematically examine the electronic characteristics of transition metal (Sc, Ti, V, Cr, Mn, Fe, and Ru) doped WSi2N4(TM@WSi2N4) and its potential as oxygen evolution reaction (OER) catalysts. Our study shows that the doping of TM atoms significantly improves the catalytic performance of TM@WSi2N4, especially Fe@WSi2N4shows a low overpotential (ηOER= 470 mV). Interestingly, we found that integrated-crystal orbital Hamilton population and d-band center can be used as descriptors to explain the high catalytic activity of Fe@WSi2N4. Subsequently, Fe@WSi2N4exhibits the best hydrogen evolution reaction (HER) activity with a universal overpotential of 47 mV on N1sites. According to our research, Fe@WSi2N4offers a promising substitute for precious metals as a catalyst for overall water splitting with low OER and HER overpotentials.
Collapse
Affiliation(s)
- Mengya Huang
- College of Big Data and Information Engineering, Guizhou University, Guiyang 550025, People's Republic of China
- Key Laboratory of Micro-Nano-Electronics of Guizhou Province, Guiyang 550025, People's Republic of China
- College of Big Health, Guizhou Medical University, Guiyang 550025, People's Republic of China
| | - Yan Jiang
- College of Big Data and Information Engineering, Guizhou University, Guiyang 550025, People's Republic of China
- Key Laboratory of Micro-Nano-Electronics of Guizhou Province, Guiyang 550025, People's Republic of China
| | - Zijiang Luo
- School of Information, Guizhou University of Finance and Economics, Guiyang 550025, People's Republic of China
| | - Jihong Wang
- College of Big Data and Information Engineering, Guizhou University, Guiyang 550025, People's Republic of China
| | - Zhao Ding
- College of Big Data and Information Engineering, Guizhou University, Guiyang 550025, People's Republic of China
- Power Semiconductor Device Reliability Research Center of the Ministry of Education, Guizhou University, Guiyang 550025, People's Republic of China
- Key Laboratory of Micro-Nano-Electronics of Guizhou Province, Guiyang 550025, People's Republic of China
| | - Xiang Guo
- College of Big Data and Information Engineering, Guizhou University, Guiyang 550025, People's Republic of China
- Power Semiconductor Device Reliability Research Center of the Ministry of Education, Guizhou University, Guiyang 550025, People's Republic of China
- Key Laboratory of Micro-Nano-Electronics of Guizhou Province, Guiyang 550025, People's Republic of China
| | - Xuefei Liu
- School of Physics and Electronic Science, Guizhou Normal University, Guiyang 550025, People's Republic of China
| | - Yi Wang
- College of Big Data and Information Engineering, Guizhou University, Guiyang 550025, People's Republic of China
- Power Semiconductor Device Reliability Research Center of the Ministry of Education, Guizhou University, Guiyang 550025, People's Republic of China
- Key Laboratory of Micro-Nano-Electronics of Guizhou Province, Guiyang 550025, People's Republic of China
| |
Collapse
|
11
|
Zhou K, Qian C, Liu Y. Quantifying the Structure of Water and Hydrated Monovalent Ions by Density Functional Theory-Based Molecular Dynamics. J Phys Chem B 2022; 126:10471-10480. [PMID: 36451081 DOI: 10.1021/acs.jpcb.2c05330] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/05/2022]
Abstract
The accurate description of the structures of water and hydrated ions is important in electrochemical desalination, ion separation, and supercapacitors. In this work, we present an ab initio atomistic simulation-based study to explore the structure of water and hydrated monovalent ions (Li+, Na+, K+, Rb+, F-, and Cl-) at ambient conditions using generalized gradient approximation (GGA)-based methods with and without van der Waals correction (PBE, PBE + D3, and revPBE + D3) and recently developed strongly constrained and appropriately normed (SCAN) meta-GGA. We find that both revPBE + D3 and SCAN can well capture the structure of bulk water with +30 K artificial high temperature in contrast to overstructuring water using PBE and PBE + D3. However, being the same as PBE + D3, revPBE + D3 overestimates the structure of the hydration shell, especially for monovalent cations. Surprisingly, SCAN can well match the experimental results of hydrated monovalent ions. Detailed structure analyzes of entropy reveal that the hydration shell under the level of PBE + D3 and revPBE + D3 is more disordered and looser than SCAN. The successful prediction of the flexible SCAN functional could facilitate the exploration of complex ionic processes in the aqueous phase, the interactions of hydrated ions with surfaces, and solvation states in nanopores at an accurate, efficient, predictive, and ab initio level.
Collapse
Affiliation(s)
- Ke Zhou
- College of Energy, Soochow Institute for Energy and Materials InnovationS (SIEMIS), Jiangsu Provincial Key Laboratory for Advanced Carbon Materials and Wearable Energy Technologies, Soochow University, Suzhou215006, China.,Laboratory for Multiscale Mechanics and Medical Science, SV LAB, School of Aerospace, Xi'an Jiaotong University, Xi'an710049, China
| | - Chen Qian
- Department of Mechanical Engineering, Zhejiang University, Hangzhou310058, China
| | - Yilun Liu
- Laboratory for Multiscale Mechanics and Medical Science, SV LAB, School of Aerospace, Xi'an Jiaotong University, Xi'an710049, China
| |
Collapse
|
12
|
Kundu A, Mamatkulov SI, Brünig FN, Bonthuis DJ, Netz RR, Elsaesser T, Fingerhut BP. Short-Range Cooperative Slow-down of Water Solvation Dynamics Around SO 4 2--Mg 2+ Ion Pairs. ACS PHYSICAL CHEMISTRY AU 2022; 2:506-514. [PMID: 36465835 PMCID: PMC9706802 DOI: 10.1021/acsphyschemau.2c00034] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 08/01/2022] [Revised: 09/15/2022] [Accepted: 09/15/2022] [Indexed: 06/17/2023]
Abstract
The presence of ions affects the structure and dynamics of water on a multitude of length and time scales. In this context, pairs of Mg2+ and SO4 2- ions in water constitute a prototypical system for which conflicting pictures of hydration geometries and dynamics have been reported. Key issues are the molecular pair and solvation shell geometries, the spatial range of electric interactions, and their impact on solvation dynamics. Here, we introduce asymmetric SO4 2- stretching vibrations as new and most specific local probes of solvation dynamics that allow to access ion hydration dynamics at the dilute concentration (0.2 M) of a native electrolyte environment. Highly sensitive heterodyne 2D-IR spectroscopy in the fingerprint region of the SO4 2- ions around 1100 cm-1 reveals a specific slow-down of solvation dynamics for hydrated MgSO4 and for Na2SO4 in the presence of Mg2+ ions, which manifests as a retardation of spectral diffusion compared to aqueous Na2SO4 solutions in the absence of Mg2+ ions. Extensive molecular dynamics and density functional theory QM/MM simulations provide a microscopic view of the observed ultrafast dephasing and hydration dynamics. They suggest a molecular picture where the slow-down of hydration dynamics arises from the structural peculiarities of solvent-shared SO4 2--Mg2+ ion pairs.
Collapse
Affiliation(s)
- Achintya Kundu
- Max-Born-Institut
Für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin12489, Germany
| | - Shavkat I. Mamatkulov
- Institute
of Material Sciences of Uzbekistan Academy of Sciences, Tashkent100084, Uzbekistan
| | | | - Douwe Jan Bonthuis
- Institute
of Theoretical and Computational Physics, Graz University of Technology, Graz8010, Austria
| | - Roland R. Netz
- Fachbereich
Physik, Freie Universität Berlin, Berlin14195, Germany
| | - Thomas Elsaesser
- Max-Born-Institut
Für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin12489, Germany
| | - Benjamin P. Fingerhut
- Max-Born-Institut
Für Nichtlineare Optik und Kurzzeitspektroskopie, Berlin12489, Germany
| |
Collapse
|
13
|
Xi C, Zheng F, Gao G, Song Z, Zhang B, Dong C, Du XW, Wang LW. Ion Solvation Free Energy Calculation Based on Ab Initio Molecular Dynamics Using a Hybrid Solvent Model. J Chem Theory Comput 2022; 18:6878-6891. [PMID: 36253911 DOI: 10.1021/acs.jctc.1c01298] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Free energy calculation of small molecules or ion species in aqueous solvent is one of the most important problems in electrochemistry study. Although there are many previous approaches to calculate such free energies, they are based on ab initio methods and suffer from various limitations and approximations. In the current work, we developed a hybrid approach based on ab initio molecular dynamics (AIMD) simulations to calculate the ion solvation energy. In this approach, a small water cluster surrounding the central ion is used, and implicit solvent model is applied outside the water cluster. A dynamic potential well is used during AIMD to keep the water cluster together. Quasi-harmonic approximation is used to calculate the entropy contribution, while the total energy average is used to calculate the enthalpy term. The obtained solvation voltages of the bulk metal agree with experiments within 0.3 eV, and the simulation results for the solvation energies of gaseous ions are close to the experimental observations. Besides the free energies, radial pair distribution functions and coordination numbers of hydrated cations are also obtained. The remaining challenges of this method are also discussed.
Collapse
Affiliation(s)
- Cong Xi
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States.,Institute of New Energy Materials, School of Materials Science and Engineering, Tianjin University, Tianjin30072, People's Republic of China
| | - Fan Zheng
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Guoping Gao
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Zhigang Song
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Buyu Zhang
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| | - Cunku Dong
- Institute of New Energy Materials, School of Materials Science and Engineering, Tianjin University, Tianjin30072, People's Republic of China
| | - Xi-Wen Du
- Institute of New Energy Materials, School of Materials Science and Engineering, Tianjin University, Tianjin30072, People's Republic of China
| | - Lin-Wang Wang
- Materials Science Division, Lawrence Berkeley National Laboratory, Berkeley, California94720, United States
| |
Collapse
|
14
|
Khanal R, Irle S. Quantum chemical investigation of the effect of alkali metal ions on the dynamic structure of water in aqueous solutions. RSC Adv 2022; 12:25500-25510. [PMID: 36275866 PMCID: PMC9480497 DOI: 10.1039/d2ra04563j] [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: 07/22/2022] [Accepted: 08/26/2022] [Indexed: 11/21/2022] Open
Abstract
We report quantum chemical molecular dynamics (MD) simulations based on the density-functional tight-binding (DFTB) method to investigate the effect of K+, Na+, and Mg2+ ions in aqueous solutions on the static and dynamic structure of bulk water at room temperature and with various concentrations. The DFTB/MD simulations were validated for the description of ion solvation in aqueous ionic solutions by comparing static pair distribution functions (PDFs) as well as the cation solvation shell between experimental and available ab initio DFT data. The effect of the cations on the water structure, as well as relative differences between K+, Na+, and Mg2+ cations, were analyzed in terms of atomically resolved PDFs as well as time-dependent Van Hove correlation functions (VHFs). The investigation of the VHFs reveals that salt ions generally slow down the dynamic decay of the pair correlations in the water solvation sphere, irrespective of the cation size or charge. The analysis of partial metal-oxygen VHFs indicates that there are long-lived correlations between water and Na+ over long distances, in contrast to K+ and Mg2+.
Collapse
Affiliation(s)
- Rabi Khanal
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory Oak Ridge Tennessee 37831 USA
| | - Stephan Irle
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory Oak Ridge Tennessee 37831 USA
| |
Collapse
|
15
|
Saha S, Dutta B, Ghosh M, Chowdhury J. Adsorption of 4-Mercapto Pyridine with Gold Nanoparticles Embedded in the Langmuir-Blodgett Film Matrix of Stearic Acid: SERS, XPS Studies Aided by Born-Oppenheimer on the Fly Dynamics, Time-Resolved Wavelet Transform Theory, and DFT. ACS OMEGA 2022; 7:27818-27830. [PMID: 35990435 PMCID: PMC9386704 DOI: 10.1021/acsomega.1c07321] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/15/2023]
Abstract
This paper reports the adsorptive behavior of the 4-mercaptopyridine (4MPy) molecule soaked in gold nanoparticles (AuNPs) that remain embedded in the bilayer Langmuir-Blodgett (LB) film matrix of stearic acid (SA) for various soaking times (STs). The as-fabricated substrate proved to be an efficient SERS sensing platform that can sense the analyte 4MPy molecules at trace concentrations of ∼1.0 × 10-9 M. The XPS study not only reveals the adsorption of 4Mpy molecules with AuNPs via a sulfur atom but also suggests partial degradation of the analyte molecule upon adsorption. This observation is further substantiated from the SERS spectral profile, which shows unusual broadening of the enhanced Raman bands of the molecule at higher STs. The experimental observations are supported by Born-Oppenheimer on-the-fly molecular dynamics (BO-OF-MD), time-resolved wavelet transform theory (WT), and the DFT calculations based on adcluster models. Selective enhancements of Raman bands in the SERS spectra further suggest the involvement of charge transfer (CT) interaction to the overall enhancements of Raman bands of the analyte molecule. The molecule → CT contribution has been estimated from electron density difference calculations and the corresponding CT distance; the amount of CT is also envisaged.
Collapse
Affiliation(s)
- Somsubhra Saha
- Department of Physics, Jadavpur University, 188, Raja S.C. Mallick Rd, Kolkata 700032, India
| | - Bipan Dutta
- Department of Physics, Sammilani Mahavidyalaya, Baghajatin Station, E. M. Bypass, Kolkata 700094, India
| | - Manash Ghosh
- Department of Spectroscopy, Indian Association for the Cultivation of Science, Jadavpur, Kolkata 700032, India
| | - Joydeep Chowdhury
- Department of Physics, Jadavpur University, 188, Raja S.C. Mallick Rd, Kolkata 700032, India
| |
Collapse
|
16
|
Smirnov PR, Grechin OV. Structure of the Immediate Environment of Ions in Aqueous Solutions of Calcium Chloride, Based on Data from X-ray Diffraction Analysis. RUSSIAN JOURNAL OF PHYSICAL CHEMISTRY A 2022. [DOI: 10.1134/s0036024422040288] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/23/2022]
|
17
|
Kröger J, Podjaski F, Savasci G, Moudrakovski I, Jiménez-Solano A, Terban MW, Bette S, Duppel V, Joos M, Senocrate A, Dinnebier R, Ochsenfeld C, Lotsch BV. Conductivity Mechanism in Ionic 2D Carbon Nitrides: From Hydrated Ion Motion to Enhanced Photocatalysis. ADVANCED MATERIALS (DEERFIELD BEACH, FLA.) 2022; 34:e2107061. [PMID: 34870342 DOI: 10.1002/adma.202107061] [Citation(s) in RCA: 19] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/06/2021] [Revised: 11/20/2021] [Indexed: 05/12/2023]
Abstract
Carbon nitrides are among the most studied materials for photocatalysis; however, limitations arise from inefficient charge separation and transport within the material. Here, this aspect is addressed in the 2D carbon nitride poly(heptazine imide) (PHI) by investigating the influence of various counterions, such as M = Li+ , Na+ , K+ , Cs+ , Ba2+ , NH4 + , and tetramethyl ammonium, on the material's conductivity and photocatalytic activity. These ions in the PHI pores affect the stacking of the 2D layers, which further influences the predominantly ionic conductivity in M-PHI. Na-containing PHI outperforms the other M-PHIs in various relative humidity (RH) environments (0-42%RH) in terms of conductivity, likely due to pore-channel geometry and size of the (hydrated) ion. With increasing RH, the ionic conductivity increases by 4-5 orders of magnitude (for Na-PHI up to 10-5 S cm-1 at 42%RH). At the same time, the highest photocatalytic hydrogen evolution rate is observed for Na-PHI, which is mirrored by increased photogenerated charge-carrier lifetimes, pointing to efficient charge-carrier stabilization by, e.g., mobile ions. These results indicate that also ionic conductivity is an important parameter that can influence the photocatalytic activity. Besides, RH-dependent ionic conductivity is of high interest for separators, membranes, or sensors.
Collapse
Affiliation(s)
- Julia Kröger
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569, Stuttgart, Germany
- Department of Chemistry, University of Munich, LMU, Butenandtstr. 5-13, 81377, Munich, Germany
- Cluster of Excellence E-Conversion, Lichtenbergstr. 4a, 85748, Garching, Germany
| | - Filip Podjaski
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569, Stuttgart, Germany
- Cluster of Excellence E-Conversion, Lichtenbergstr. 4a, 85748, Garching, Germany
| | - Gökcen Savasci
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569, Stuttgart, Germany
- Department of Chemistry, University of Munich, LMU, Butenandtstr. 5-13, 81377, Munich, Germany
- Cluster of Excellence E-Conversion, Lichtenbergstr. 4a, 85748, Garching, Germany
| | - Igor Moudrakovski
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569, Stuttgart, Germany
| | - Alberto Jiménez-Solano
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569, Stuttgart, Germany
| | - Maxwell W Terban
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569, Stuttgart, Germany
| | - Sebastian Bette
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569, Stuttgart, Germany
| | - Viola Duppel
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569, Stuttgart, Germany
| | - Markus Joos
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569, Stuttgart, Germany
| | - Alessandro Senocrate
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569, Stuttgart, Germany
| | - Robert Dinnebier
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569, Stuttgart, Germany
| | - Christian Ochsenfeld
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569, Stuttgart, Germany
- Department of Chemistry, University of Munich, LMU, Butenandtstr. 5-13, 81377, Munich, Germany
- Cluster of Excellence E-Conversion, Lichtenbergstr. 4a, 85748, Garching, Germany
| | - Bettina V Lotsch
- Max Planck Institute for Solid State Research, Heisenbergstr. 1, 70569, Stuttgart, Germany
- Department of Chemistry, University of Munich, LMU, Butenandtstr. 5-13, 81377, Munich, Germany
- Cluster of Excellence E-Conversion, Lichtenbergstr. 4a, 85748, Garching, Germany
| |
Collapse
|
18
|
Hydration and aggregation in aqueous xylitol solutions in the wide temperature range. J Mol Liq 2021. [DOI: 10.1016/j.molliq.2021.116838] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
19
|
Midgley S, Di Tommaso D, Fleitmann D, Grau-Crespo R. Sulfate and Molybdate Incorporation at the Calcite-Water Interface: Insights from Ab Initio Molecular Dynamics. ACS EARTH & SPACE CHEMISTRY 2021; 5:2066-2073. [PMID: 34476320 PMCID: PMC8397469 DOI: 10.1021/acsearthspacechem.1c00131] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Revised: 06/29/2021] [Accepted: 07/13/2021] [Indexed: 06/13/2023]
Abstract
Sulfur and molybdenum trace impurities in speleothems (stalagmites and stalactites) can provide long and continuous records of volcanic activity, which are important for past climatic and environmental reconstructions. However, the chemistry governing the incorporation of the trace element-bearing species into the calcium carbonate phases forming speleothems is not well understood. Our previous work has shown that substitution of tetrahedral oxyanions [XO4]2- (X = S and Mo) replacing [CO3]2- in CaCO3 bulk phases (except perhaps for vaterite) is thermodynamically unfavorable with respect to the formation of competing phases, due to the larger size and different shape of the [XO4]2- tetrahedral anions in comparison with the flat [CO3]2- anions, which implied that most of the incorporation would happen at the surface rather than at the bulk of the mineral. Here, we present an ab initio molecular dynamics study, exploring the incorporation of these impurities at the mineral-water interface. We show that the oxyanion substitution at the aqueous calcite (10.4) surface is clearly favored over bulk incorporation, due to the lower structural strain on the calcium carbonate solid. Incorporation at surface step sites is even more favorable for both oxyanions, thanks to the additional interface space afforded by the surface line defect to accommodate the tetrahedral anion. Differences between sulfate and molybdate substitutions can be mostly explained by the size of the anions. The molybdate oxyanion is more difficult to incorporate in the calcite bulk than the smaller sulfate oxyanion. However, when molybdate is substituted at the surface, the elastic cost is avoided because the oxyanion protrudes out of the surface and gains stability via the interaction with water at the interface, which in balance results in more favorable surface substitution for molybdate than for sulfate. The detailed molecular-level insights provided by our calculations will be useful to understand the chemical basis of S- and Mo-based speleothem records.
Collapse
Affiliation(s)
- Scott
D. Midgley
- Department
of Chemistry, University of Reading, Whiteknights, Reading RG6 6DX, U.K.
| | - Devis Di Tommaso
- Department
of Chemistry, School of Biological and Chemical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, U.K.
| | - Dominik Fleitmann
- Department
of Archaeology, University of Reading, Whiteknights, Reading RG6 6AB, U.K.
- Department
of Environmental Sciences, University of
Basel, Bernoullistrasse 32, Basel 4056, Switzerland
| | - Ricardo Grau-Crespo
- Department
of Chemistry, University of Reading, Whiteknights, Reading RG6 6DX, U.K.
| |
Collapse
|
20
|
Hayamizu K, Chiba Y, Haishi T. Dynamic ionic radius of alkali metal ions in aqueous solution: a pulsed-field gradient NMR study. RSC Adv 2021; 11:20252-20257. [PMID: 35479919 PMCID: PMC9033755 DOI: 10.1039/d1ra02301b] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/23/2021] [Accepted: 05/19/2021] [Indexed: 01/26/2023] Open
Abstract
The dynamic behavior of alkali metal ions, Li+, Na+, K+, Rb+ and Cs+ in aqueous solutions is one of the most important topics in solution chemistry. Since these alkali metals contain nuclear magnetic resonance (NMR) active nuclei, it is possible to directly measure the diffusion constants of the alkali metal ions using the pulsed field gradient (PFG) NMR method. In this paper, the 7Li, 23Na, 87Rb, 133Cs and 1H resonances are observed for diffusion constants in aqueous solution and the solvent H2O. Until now, the values of the diffusion constant have been lacking when discussing hydration effects around alkali metal ions. It is known that the static ionic radius (Rion) increases with increasing the atomic number, and the experimental diffusion constants also increase with increasing the atomic number, which is opposite to the Stokes–Einstein (SE) relation. It suggests that alkali metal ions diffuse through a space of 10−6 m accompanying the hydrated spheres with a time interval of 10−3 s. For each alkali metal ion, the dynamic ionic radius is evaluated. Stokes radius (dynamic ionic radius) of the alkali metal ions versus the ionic radius (Rion) at 303 K. The dotted line is a guide for the 1 : 1 relation.![]()
Collapse
Affiliation(s)
- Kikuko Hayamizu
- Institute of Applied Physics, Tsukuba University Tennodai Tsukuba 305-8573 Japan
| | - Yusuke Chiba
- Graduate School of Pure and Applied Sciences and Tsukuba Research Center for Energy, Materials Science (TREMS), University of Tsukuba Tennodai Tsukuba 305-8573 Japan
| | | |
Collapse
|
21
|
Zhan C, Sun Y, Aydin F, Wang YM, Pham TA. Confinement effects on the solvation structure of solvated alkaline metal cations in a single-digit 1T-MoS 2 nanochannel: A first-principles study. J Chem Phys 2021; 154:164706. [PMID: 33940836 DOI: 10.1063/5.0047554] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022] Open
Abstract
Confinement plays an important role in determining ion transport in porous materials, which, in turn, may influence the performance of many energy storage and desalination devices. In this work, we combined density functional theory (DFT) with an implicit solvation model and ab initio molecular dynamics (AIMD) to investigate the effects of nanoconfinement on several solvated alkaline metal cations in a single-digit 1T-MoS2 nanochannel. Our DFT calculations with a solvation model indicated that cations with stronger hydration energy introduce a higher number of co-intercalated water molecules into the channel, consistent with early experimental observation obtained for MXene (2D transition metal carbide) channels. The predicted optimal water numbers for the cations were then used for AIMD simulations that explicitly include the effects of the solvent. When compared with the cations in bulk solution, our simulations showed that the hydration structure and coordination number (CN) of the solvated cations confined in the MoS2 channel can be significantly altered. We found that larger cations with weaker hydration energy (K+, Rb+, and Cs+) exhibited a distinctive CN decrease under confinement, while smaller cations (Li+ and Na+) retained a similar hydration shell as in the bulk solution. More specifically, the hydration shell of large cations (K+, Rb+, and Cs+) in MoS2 showed similar features of the coordination angle to the bulk, which suggests the partially broken hydration shell with no geometry change under confinement. Our simulations provided insights into the change of the hydration structure of alkaline metal cations under confinement, which may have important implications on their transport in the 1T-MoS2 channel.
Collapse
Affiliation(s)
- Cheng Zhan
- Material Science Division, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Yangyunli Sun
- Department of Chemistry, University of California, Riverside, California 92521, USA
| | - Fikret Aydin
- Material Science Division, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| | - Y Morris Wang
- Department of Materials Science and Engineering, University of California, Log Angeles, California 90095, USA
| | - Tuan Anh Pham
- Material Science Division, Lawrence Livermore National Laboratory, Livermore, California 94550, USA
| |
Collapse
|
22
|
Assessing the Intrinsic Strengths of Ion–Solvent and Solvent–Solvent Interactions for Hydrated Mg2+ Clusters. INORGANICS 2021. [DOI: 10.3390/inorganics9050031] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/17/2022] Open
Abstract
Information resulting from a comprehensive investigation into the intrinsic strengths of hydrated divalent magnesium clusters is useful for elucidating the role of aqueous solvents on the Mg2+ ion, which can be related to those in bulk aqueous solution. However, the intrinsic Mg–O and intermolecular hydrogen bond interactions of hydrated magnesium ion clusters have yet to be quantitatively measured. In this work, we investigated a set of 17 hydrated divalent magnesium clusters by means of local vibrational mode force constants calculated at the ωB97X-D/6-311++G(d,p) level of theory, where the nature of the ion–solvent and solvent–solvent interactions were interpreted from topological electron density analysis and natural population analysis. We found the intrinsic strength of inner shell Mg–O interactions for [Mg(H2O)n]2+ (n = 1–6) clusters to relate to the electron density at the bond critical point in Mg–O bonds. From the application of a secondary hydration shell to [Mg(H2O)n]2+ (n = 5–6) clusters, stronger Mg–O interactions were observed to correspond to larger instances of charge transfer between the lp(O) orbitals of the inner hydration shell and the unfilled valence shell of Mg. As the charge transfer between water molecules of the first and second solvent shell increased, so did the strength of their intermolecular hydrogen bonds (HBs). Cumulative local vibrational mode force constants of explicitly solvated Mg2+, having an outer hydration shell, reveal a CN of 5, rather than a CN of 6, to yield slightly more stable configurations in some instances. However, the cumulative local mode stretching force constants of implicitly solvated Mg2+ show the six-coordinated cluster to be the most stable. These results show that such intrinsic bond strength measures for Mg–O and HBs offer an effective way for determining the coordination number of hydrated magnesium ion clusters.
Collapse
|
23
|
Smirnov PR. Structural Parameters of the Nearest Surrounding of Group II
Metal Ions in Oxygen-Containing Solvents. RUSS J GEN CHEM+ 2021. [DOI: 10.1134/s1070363221030129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/22/2022]
|