Dynamics of halide ion-water hydrogen bonds in aqueous solutions: dependence on ion size and temperature

Simulating the Competitive Ion Pairing of Hydrated Electrons with Chaotropic Cations

Experiments show that the absorption spectrum of the hydrated electron (ehyd-) blue-shifts in electrolyte solutions compared with what is seen in pure water. This shift has been assigned to the ehyd-'s competitive ion-pairing interactions with the salt cation relative to the salt anion based on the ions' positions on the Hofmeister series. Remarkably, little work has been done investigating the ehyd-'s behavior when the salts have chaotropic cations, which should greatly change the ion-pairing interactions given that the ehyd- is a champion chaotrope. In this work, we remedy this by using mixed quantum/classical simulations to analyze the behavior of two different models of the ehyd- in aqueous RbF and RbI electrolyte solutions as a function of salt concentration. We find that the magnitude of the salt-induced spectral blue-shift is determined by a combination of the number of chaotropic Rb+ cations near the ehyd- and the number of salt anions near those cations so that the spectrum of the ehyd- directly reflects its local environment. We also find that the use of a soft-cavity ehyd- model predicts stronger competitive interactions with Rb+ relative to I- than a more traditional hard cavity model, leading to different predicted spectral shifts that should provide a way to distinguish between the two models experimentally. Our simulations predict that at the same concentration, salts with chaotropic cations should produce larger spectral blue-shifts than salts with kosmotropic cations. We also found that at high salt concentrations with chaotropic cations, the predicted blue-shift is greater when the salt anion is kosmotropic instead of chaotropic. Our goal is for this work to inspire experimentalists to make such measurements, which will help provide a spectroscopic means to distinguish between simulations models that predict different hydration structures for the ehyd-.

Study on the effect of volatile organic compounds on the treatment of high-salt wastewater by low-temperature evaporation

High-salinity wastewater, owing to its intricate composition and challenging treatment requirements, poses a significant hurdle in water environmental governance. In this study, low-temperature evaporation technology is used to tackle wastewater containing the volatile organic compound such as N,N-dimethylacetamide (DMAC). Utilisation of comprehensive approaches involving experimental testing, mathematical modelling, and Aspen Plus software simulations, The influence of DMAC on evaporation efficiency is researched through the following factors which encompassing its effects on boiling point elevation, partial molar activation energy, and the formation of by-products. Additionally, the comparation of the impact of temperature, ionic strength, intermolecular interactions on the evaporation rate and the concentration of the volatile component DMAC in the condensate is also conducted in this study. After conducting a multiple linear regression analysis of evaporation efficiency using the Statistical Product and Service Solutions (SPSS) tool, it was discovered that temperature serves as the primary determinant influencing the evaporation rate. Additionally, ionic strength impacts solution viscosity, intermolecular interactions, and saturated vapour pressure by altering the intermolecular forces, thereby indirectly influencing both the evaporation rate and the quality of condensate water. The comparative analysis of single-effect and double-effect evaporation indicates that the optimal operating condition for double-effect evaporation yields an evaporation rate of 70%, with a remarkable 88% reduction in steam consumption compared to single one. Based on heat and mass balance principles, the mathematical model for double-effect evaporation is established to offer crucial data support for practical industrial applications.

Fluoride-Ion-Responsive Sol-Gel Transition in an L-Cysteine/AgNO3 System: Self-Assembly Peculiarities and Anticancer Activity

Supramolecular hydrogels based on low-molecular-weight compounds are a unique class of so-called "soft" materials, formed by weak non-covalent interactions between precursors at their millimolar concentrations. Due to the variety of structures that can be formed using different low-molecular-weight gelators, they are widely used in various fields of technology and medicine. In this study, we report for the first time an unusual self-assembly process of mixing a hydrosol obtained from L-cysteine and silver nitrate (cysteine-silver sol-CSS) with sodium halides. Modern instrumental techniques such as viscosimetry, UV spectroscopy, dynamic light scattering, zeta potential measurements, SEM and EDS identified that adding fluoride anions to CSS is able to form stable hydrogels of a thixotropic nature, while Cl-, Br- and I- lead to precipitation. The self-assembly process proceeds using a narrow concentration range of F-. An increase in the fluoride anion content in the system leads to a change in the gel network morphology from elongated structures to spherical ones. This fact is reflected in a decrease in the gel viscosity and a number of gel-sol-gel transition cycles. The mechanism of F-'s interaction with hydrosol includes the condensation of anions on the positive surface of the CSS nanoparticles, their binding via electrostatic forces and the formation of a resulting gel carcass. In vitro analysis showed that the hydrogels suppressed human squamous carcinoma cells at a micromolar sample concentration. The obtained soft gels could have potential applications against cutaneous malignancy and as carriers for fluoride anion and other bioactive substance delivery.

Water-in-Salt: Fast Dynamics, Structure, Thermodynamics, and Bulk Properties

We present concentration-dependent dynamics of highly concentrated LiBr solutions and LiCl temperature-dependent dynamics for two high concentrations and compare the results to those of prior LiCl concentration-dependent data. The dynamical data are obtained using ultrafast optical heterodyne-detected optical Kerr effect (OHD-OKE). The OHD-OKE decays are composed of two pairs of biexponentials, i.e., tetra-exponentials. The fastest decay (t1) is the same as pure water's at all concentrations within error, while the second component (t2) slows slightly with concentration. The slower components (t3 and t4), not present in pure water, slow substantially, and their contributions to the decays increase significantly with increasing concentration, similar to LiCl solutions. Simulations of LiCl solutions from the literature show that the slow components arise from large ion/water clusters, while the fast components are from ion/water structures that are not part of large clusters. Temperature-dependent studies (15-95 °C) of two high LiCl concentrations show that decreasing the temperature is equivalent to increasing the room temperature concentration. The LiBr and LiCl concentration dependences and the two LiCl concentrations' temperature dependences all have bulk viscosities that are linearly dependent on τcslow, the correlation time of the slow dynamics (weighted averages of t3 and t4). Remarkably, all four viscosity vs 1/τCslow plots fall on the same line. Application of transition state theory to the temperature-dependent data yields the activation enthalpies and entropies for the dynamics of the large ion/water clusters, which underpin the bulk viscosity.

Effects of Hydrogen Peroxide on the Hydrogen Bonding Structure and Dynamics of Water and Its Influence on the Aqueous Solvation of the Insulin Monomer

The hydrogen bond structure and dynamics of water and hydrogen peroxide (H2O2) in their binary mixtures have been studied at 298 K by classical molecular dynamics simulations. Twelve different concentrations of aqueous-H2O2 solutions are considered for this study. We have analyzed the interactions between water and H2O2 by site-site pair correlation functions and observed that the probability of formation of OW···HP hydrogen bonds are higher compared to OP···HW. The second solvation shell of water is strongly affected by increasing H2O2 concentrations (XP > 0.50), which signifies the destruction of the tetrahedral network structure of water. The translational and rotational dynamics of water and H2O2 do not significantly change up to 25% of H2O2 in aqueous mixtures. The hydrogen bond lifetime of water-water, water-H2O2, and H2O2-H2O2 in the aqueous-H2O2 solutions shows a very minimal change with increasing H2O2 concentrations. In addition to this, we also investigated the effect of H2O2 on the insulin monomer and observed that higher concentrations of H2O2 (XP = 0.10) change the secondary structure. The influence of H2O2 is more on chain-B than that on chain-A in the insulin monomer. The H2O2 occupancy at the protein surface is higher for negatively charged (GLU) and polar (ASN and THR) amino acid residues compared with that for positively charged and neutral residues in the solutions.

Mechanism, Kinetics, and Potential of Mean Force of Evaporation of Water from Aqueous Sodium Chloride Solutions of Varying Concentrations

The mechanism, kinetics, and potential of mean force of evaporation of water from aqueous NaCl solutions are investigated through both unbiased molecular dynamics simulations and also biased simulations using the umbrella sampling method. The results are obtained for aqueous solutions of three different NaCl concentrations ranging from 0.6 to 6.0 m and also for pure water. The rate of evaporation is found to decrease in the presence of ions. It is found that the process of evaporation of a surface water molecule from ionic solutions can be triggered through its collision with another water or chloride ion. Such collisions provide the additional kinetic energy that is required for evaporation. However, when the collision takes place with a Cl^- ion, the evaporation of the escaping water also involves a collision with water in the vicinity of the ion at the same time along with the ion-water collision. These two collisions together provide the required kinetic energy for escape of the evaporating water molecule. Thus, the mechanism of evaporation process of ionic solutions can be more complex than that of pure water. The potential of mean force (PMF) of evaporation is found to be positive and it increases with increasing ion concentration. Also, no barrier in the PMF is found to be present for the condensation of water from vapor phase to the surfaces of the solutions. A detailed analysis of the unsuccessful evaporation attempts by surface water molecules is also made in the current study.

Residence Times from Molecular Dynamics Simulations

In this work, efficient methods are proposed for the calculation, from molecular dynamics trajectories, of residence times (RTs) and related quantities. One of these was designed to obtain RT distributions, from which mean residence times (MRTs), residual time distributions, and mean residual times can be computed. This method does not require the assumptions and approximations made by the most commonly used methods. Its link to the most popular method in the literature is identified. It is shown how the much faster new method can be used as a replacement for the latter and the advantages of doing so. Also, a prescription for estimating the uncertainty in the MRTs obtained though the proposed method is provided. Another even faster method for the calculation of the MRTs, their uncertainties, and the mean residual times is also proposed. It yields exactly the same results as the first one but does not allow to obtain the mentioned distributions. Another very popular method, based on autocorrelation functions, for computing MRTs is analyzed in terms of its assumptions and approximations. An alternative, also based on autocorrelation functions, which is faster and requires fewer assumptions is presented. A prescription for the calculation of the uncertainty of the MRTs obtained with the latter method is also provided. In the literature, there are a few methods to discard short transient escapes. Here, an algorithm is suggested that is much faster than the most used one and allows a more integral treatment of the process. Also, it is more widely applicable because it is a preprocessing step that can be used in conjunction with any of the proposed methods mentioned above. The main disadvantage of these two approaches to discard brief escapes is that the maximum duration allowed for an escape to be considered transient appears as a parameter to be determined for the particular system under study. As an alternative, a parameter-free method of a similar character is also proposed to estimate the mean residence time of particles that reached a constant probability of leaving the region of interest.

Decontamination of water co-polluted by copper, toluene and tetrahydrofuran using lauric acid

Co-contamination by organic solvents (e.g., toluene and tetrahydrofuran) and metal ions (e.g., Cu2+) is common in industrial wastewater and in industrial sites. This manuscript describes the separation of THF from water in the absence of copper ions, as well as the treatment of water co-polluted with either THF and copper, or toluene and copper. Tetrahydrofuran (THF) and water are freely miscible in the absence of lauric acid. Lauric acid separates the two solvents, as demonstrated by proton nuclear magnetic resonance (1H NMR) and Attenuated Total Reflection-Fourier Transform Infrared Spectroscopy (ATR-FTIR). The purity of the water phase separated from 3:7 (v/v) THF:water mixtures using 1 M lauric acid is ≈87%v/v. Synchrotron small angle X-Ray scattering (SAXS) indicates that lauric acid forms reverse micelles in THF, which swell in the presence of water (to host water in their interior) and ultimately lead to two free phases: 1) THF-rich and 2) water-rich. Deprotonated lauric acid (laurate ions) also induces the migration of Cu2+ ions in either THF (following separation from water) or in toluene (immiscible in water), enabling their removal from water. Laurate ions and copper ions likely interact through physical interactions (e.g., electrostatic interactions) rather than chemical bonds, as shown by ATR-FTIR. Inductively coupled plasma-optical emission spectrometry (ICP-OES) demonstrates up to 60% removal of Cu2+ ions from water co-polluted by CuSO4 or CuCl2 and toluene. While lauric acid emulsifies water and toluene in the absence of copper ions, copper salts destabilize emulsions. This is beneficial, to avoid that copper ions are re-entrained in the water phase alongside with toluene, following their migration in the toluene phase. The effect of copper ions on emulsion stability is explained based on the decreased interfacial activity and compressional rigidity of interfacial films, probed using a Langmuir trough. In wastewater treatment, lauric acid (a powder) can be mixed directly in the polluted water. In the context of groundwater remediation, lauric acid can be solubilized in canola oil to enable its injection to treat aquifers co-polluted by organic solvents and Cu2+. In this application, injectable filters obtained by injecting cationic hydroxyethylcellulose (HEC +) would impede the flow of toluene and copper ions partitioned in it, protecting downstream receptors. Co-contaminants can be subsequently extracted upstream of the filters (using pumping wells), to enable their simultaneous removal from aquifers.

Time-resolved terahertz-Raman spectroscopy reveals that cations and anions distinctly modify intermolecular interactions of water

The solvation of ions changes the physical, chemical and thermodynamic properties of water, and the microscopic origin of this behaviour is believed to be ion-induced perturbation of water's hydrogen-bonding network. Here we provide microscopic insights into this process by monitoring the dissipation of energy in salt solutions using time-resolved terahertz-Raman spectroscopy. We resonantly drive the low-frequency rotational dynamics of water molecules using intense terahertz pulses and probe the Raman response of their intermolecular translational motions. We find that the intermolecular rotational-to-translational energy transfer is enhanced by highly charged cations and is drastically reduced by highly charged anions, scaling with the ion surface charge density and ion concentration. Our molecular dynamics simulations reveal that the water-water hydrogen-bond strength between the first and second solvation shells of cations increases, while it decreases around anions. The opposite effects of cations and anions on the intermolecular interactions of water resemble the effects of ions on the stabilization and denaturation of proteins.

Ion permeation, selectivity, and electronic polarization in fluoride channels

Fluoride channels (Flucs) export toxic F^- from the cytoplasm. Crystallography and mutagenesis have identified several conserved residues crucial for fluoride transport, but the permeation mechanism at the molecular level has remained elusive. Herein, we have applied constant-pH molecular dynamics and free-energy-sampling methods to investigate fluoride permeation through a Fluc protein from Escherichia coli. We find that fluoride is facile to permeate in its charged form, i.e., F^-, by traversing through a non-bonded network. The extraordinary F^- selectivity is gained by the hydrogen-bonding capability of the central binding site and the Coulombic filter at the channel entrance. The F^- permeation rate calculated using an electronically polarizable force field is significantly more accurate compared with the experimental value than that calculated using a more standard additive force field, suggesting an essential role for electronic polarization in the F^--Fluc interactions.

Bridging Structure, Dynamics, and Thermodynamics: An Example Study on Aqueous Potassium Halides

Aqueous salt systems are ubiquitous in all areas of life. The ions in these solutions impose important structural and dynamic perturbations to water. In this study, we employ a combined neutron scattering, nuclear magnetic resonance, and computational modeling approach to deconstruct ion-specific perturbations to water structure and dynamics and shed light on the molecular origins of bulk thermodynamic properties of the solutions. Our approach uses the atomistic scale resolution offered to us by neutron scattering and computational modeling to investigate how the properties of particular short-ranged microenvironments within aqueous systems can be related to bulk properties of the system. We find that by considering only the water molecules in the first hydration shell of the ions that the enthalpy of hydration can be determined. We also quantify the range over which ions perturb water structure by calculating the average enthalpic interaction between a central halide anion and the surrounding water molecules as a function of distance and find that the favorable anion-water enthalpic interactions only extend to ∼4 Å. We further validate this by showing that ions induce structure in their solvating water molecules by examining the distribution of dipole angles in the first hydration shell of the ions but that this perturbation does not extend into the bulk water. We then use these structural findings to justify mathematical models that allow us to examine perturbations to rotational and diffusive dynamics in the first hydration shell around the potassium halide ions from NMR measurements. This shows that as one moves down the halide series from fluorine to iodine, and ionic charge density is therefore reduced, that the enthalpy of hydration becomes less negative. The first hydration shell also becomes less well structured, and rotational and diffusive motions of the hydrating water molecules are increased. This reduction in structure and increase in dynamics are likely the origin of the previously observed increased entropy of hydration as one moves down the halide series. These results also suggest that simple monovalent potassium halide ions induce mostly local perturbations to water structure and dynamics.

A Molecular Interpretation of the Dynamics of Diffusive Mass Transport of Water within a Glassy Polyetherimide

The diffusion process of water molecules within a polyetherimide (PEI) glassy matrix has been analyzed by combining the experimental analysis of water sorption kinetics performed by FTIR spectroscopy with theoretical information gathered from Molecular Dynamics simulations and with the expression of water chemical potential provided by a non-equilibrium lattice fluid model able to describe the thermodynamics of glassy polymers. This approach allowed us to construct a convincing description of the diffusion mechanism of water in PEI providing molecular details of the process related to the effects of the cross- and self-hydrogen bonding established in the system on the dynamics of water mass transport.

Hydrogen-Bond Structure and Low-Frequency Dynamics of Electrolyte Solutions: Hydration Numbers from ab Initio Water Reorientation Dynamics and Dielectric Relaxation Spectroscopy

We present an atomistic simulation scheme for the determination of the hydration number (h) of aqueous electrolyte solutions based on the calculation of the water dipole reorientation dynamics. In this methodology, the time evolution of an aqueous electrolyte solution generated from ab initio molecular dynamics simulations is used to compute the reorientation time of different water subpopulations. The value of h is determined by considering whether the reorientation time of the water subpopulations is retarded with respect to bulk-like behavior. The application of this computational protocol to magnesium chloride (MgCl2 ) solutions at different concentrations (0.6-2.8 mol kg-1 ) gives h values in excellent agreement with experimental hydration numbers obtained using GHz-to-THz dielectric relaxation spectroscopy. This methodology is attractive because it is based on a well-defined criterion for the definition of hydration number and provides a link with the molecular-level processes responsible for affecting bulk solution behavior. Analysis of the ab initio molecular dynamics trajectories using radial distribution functions, hydrogen bonding statistics, vibrational density of states, water-water hydrogen bonding lifetimes, and water dipole reorientation reveals that MgCl2 has a considerable influence on the hydrogen bond network compared with bulk water. These effects have been assigned to the specific strong Mg-water interaction rather than the Cl-water interaction.

Solvation Shell of the Nitrite Ion in Water: An Ab Initio Molecular Dynamics Study

We performed ab initio molecular dynamics simulation of a nitrite ion in water to investigate the structural and dynamical properties of its hydration shell. The nitrite ion is found to exhibit strong asymmetry toward hydrogen bonding due to its two different types of hydrogen bond acceptor sites. This difference is better captured through further partitioning of the hydration shell into its proximal and distal regions. The frequency shifts of the stretch modes of hydration shell water reveal that the nitrogen site forms a stronger hydrogen bond than its oxygen sites with the latter forming hydrogen bonds, which are similar in strength to that between a pair of water molecules. The escape dynamics of water from the hydration shell is found to be rather slow, which seems to classify the nitrite ion as a structure-maker. However, the dynamics of orientational and hydrogen bond relaxation reveal a faster mobility of water molecules in the hydration shell than bulk water in spite of strong ion-water interactions. It is found that the nitrite ion can hold water molecules in its solvation shell and still make them rotate fast in its vicinity through switching of their hydrogen bonds between its nitrogen and oxygen acceptor sites. The dipole moment of the solute in water is also calculated in the present study.

Density functional theory based molecular dynamics study of solution composition effects on the solvation shell of metal ions

We present an ab initio molecular dynamics study of the alkali metal ions Li+, Na+, K+ and Cs+, and of the alkaline earth metal ions Mg2+ and Ca2+ in both pure water and electrolyte solutions containing the counterions Cl- and SO42-. Simulations were conducted using different density functional theory methods (PBE, BLYP and revPBE), with and without the inclusion of dispersion interactions (-D3). Analysis of the ion-water structure and interaction strength, water exchange between the first and second hydration shell, and hydrogen bond network and low-frequency reorientation dynamics around the metal ions have been used to characterise the influence of solution composition on the ionic solvation shell. Counterions affect the properties of the hydration shell not only when they are directly coordinated to the metal ion, but also when they are at the second coordination shell. Chloride ions reduce the sodium hydration shell and expand the calcium hydration shell by stabilizing under-coordinated hydrated Na(H2O)5+ complexes and over-coordinated Ca(H2O)72+. The same behaviour is observed in CaSO4(aq), where Ca2+ and SO42- form almost exclusively solvent-shared ion pairs. Water exchange between the first and second hydration shell around Ca2+ in CaSO4(aq) is drastically decelerated compared with the simulations of the hydrated metal ion (single Ca2+, no counterions). Velocity autocorrelation function analysis, used to probe the strength of the local ion-water interaction, shows a smoother decay of Mg2+ in MgCl2(aq), which is a clear indication of a looser inter-hexahedral vibration in the presence of chloride ions located in the second coordination shell of Mg2+. The hydrogen bond statistics and orientational dynamics in the ionic solvation shell show that the influence on the water-water network cannot only be ascribed to the specific cation-water interaction, but also to the subtle interplay between the level of hydration of the ions, and the interactions between ions, especially those of opposite charge. As many reactive processes involving solvated metal ions occur in environments that are far from pure water but rich in ions, this computational study shows how the solution composition can result in significant differences in behaviour and function of the ionic solvation shell.

Halogenide anions as halogen and hydrogen bond acceptors in iodopyridinium halogenides

Structures of iodopyridinium halogenides have demonstrated why iodide, the weakest halogen bond acceptor among the halogenides, preferentially forms halogen bonds.

Chloride ions as integral parts of hydrogen bonded networks in aqueous salt solutions: the appearance of solvent separated anion pairs

Characteristic hydrogen bonded motifs, including solvent separated anion pairs, in concentrated aqueous LiCl solutions.

Heterogeneous Occupancy and Vibrational Dynamics of Spatially Patterned Water Molecules

We performed first-principles molecular dynamics simulations of relatively dilute aqueous solutions of sulfate and thiosulfate dianions to analyze the structure, dynamics, and vibrational spectral properties of water molecules around the solute, especially the spatially patterned solvent molecules in the first solvation layer and the extended layers. This study also involves the investigation of dynamics of dangling OH groups in these layers and their role in patterning the water molecules around the dianions. Structural evaluation of the systems is carried out by radial distribution functions, number integrals, and spatial distribution functions. The lifetime of dangling OH groups inside the solvation shell is compared more to that of the bulk. By constructing the O-H groups in three ensembles (S1, S2, and S3) around the anion, we show that the frequency distribution of OH modes in the S1 ensemble show red-shifting for both sulfate and thiosulfate. The O-H groups in the S2 ensemble of the sulfate-water system show red-shifting by 10 cm^-1, while in the case of thiosulfate-water, these O-H groups show blue-shifting by 8 cm^-1. The water molecules in S1 and S2 subensembles have slower dynamics compared to those in the bulk (S3). The dynamics of various kinds of hydrogen bonds were characterized by hydrogen bond population correlation functions. The spectral diffusion of solvation shell O-H modes was performed through a frequency-time correlation function. We find a significant amount of orientational retardation of water molecules in the S1 layer and moderate retardation in the S2 layer as compared to that in the bulk, S3 layer. All these findings, the red shift of the OH stretching frequency in S1 and S2 layers, slowing down of the orientational dynamics of OH vectors in S1 and S2 layers, and less diffusivity of water in S1 and S2 layers, show the long-range kosmotropic effect of multivalent sulfate and thiosulfate oxyanions. Due to the long-range effect, heterogeneous occupancy of water molecules is observed, and the water molecules are found to arrange in a patterned manner in the vicinity of anions with varied local density.

Vibration Spectral Dynamics of Weakly Coordinating Water Molecules near an Anion: FPMD Simulations of an Aqueous Solution of Tetrafluoroborate

The extent to which the ions affect the nearby water molecules will decide the structure-making or breaking nature of those ions in aqueous solutions. The effects of a weakly coordinating anion on the structure, dynamics, and vibrational properties of water molecules are not so significant as compared to an anion capable of making strong ion-water hydrogen bonds. The present work deals with the first-principles molecular dynamics study of an aqueous solution of such a weakly coordinating anion, tetrafluoroborate (BF₄^-), using dispersion-corrected DFT-based first-principles molecular dynamics (FPMD) simulations. Various structural, dynamical, and spectral properties, such as radial distribution functions (RDFs), rotational dynamics, vibrational density of states (VDOS), hydrogen bond as well as dangling OH autocorrelation functions, and residence dynamics, were calculated to investigate the effects of the anion on nearby water molecules. The process of spectral diffusion was assessed through a time series wavelet transformation of trajectories obtained from FPMD simulations. The first ion-water solvation shell extends up to 5.5 Å, containing around 20 water molecules. The lifetime of the ion-water hydrogen bond is found to be 1.19 ps, whereas the water-water hydrogen bond lifetime is found to be 1.13 ps. Inside the solvation shell, the persistence time of dangling OH chromophores and the average frequency of OH modes inside the solvation shell are found to be more compared to bulk. Three time scales are found for solvation shell OH modes from the frequency-frequency correlation function. A very short time scale is found for the intact ion-water interaction; the short time scale is for the ion-water hydrogen bond, and the long time scale is for escape dynamics of water molecules from the ion solvation shell. From the mean squared displacement, it is found that solvation water molecules diffuse slower than the bulk. However, solvation shell water molecules show faster relaxation from the analysis of rotational anisotropy. Within the longer time scale of spectral diffusion, this process (which is related to various dynamics of the molecules) is not yet complete, as compared to fast anisotropic decay. This fact is similar to the experimental finding of spectral diffusion and anisotropy time scales in the aqueous solution of borohydride anion. The calculated results are also compared with available experimental data wherever possible.

Born-Oppenheimer Molecular Dynamics Simulations of a Bromate Ion in Water Reveal Its Dual Kosmotropic and Chaotropic Behavior

The solvation structure and dynamics of a bromate (BrO₃^-) ion in water are studied by means of Born-Oppenheimer molecular dynamics simulations at two different temperatures using the Becke-Lee-Yang-Parr functional with Grimme D3 dispersion corrections. The bromate ion possesses a pyramidal structure, and it has two types of solvation sites, namely, the bromine and oxygen atoms. We have looked at different radial and orientational distributions of water molecules around the bromate ion and also investigated their hydrogen bonding properties. The solvation structure of the bromate ion is also compared with that of the iodate (IO₃^-) ion, which is structurally rather similar to the bromate ion and was found to have some unusual solvation properties in water. It is found that the bromate ion follows a similar trend as that followed by the iodate ion as far as the solvation structure is concerned. However, the effect of the former on surrounding water is found to be much weaker than that of the latter. On the dynamical side, we have looked at diffusion, residence dynamics, and also the orientational and hydrogen bond relaxation of water molecules around the BrO₃^- ion and compared them with those of the bulk. Dynamical results are presented for both H₂O and D₂O around the BrO₃^- ion. Interpretation of the dynamical results in terms of structure-making (kosmotropic)/-breaking (chaotropic) properties of the BrO₃^- ion reveals that the bromine atom of this ion acts as a water structure breaker, whereas the three oxygens act as water structure makers. Thus, in spite of being a single ion, the bromate ion has dual characteristics and the experimentally observed kosmotropic ability of this ion is actually a trade-off between a chaotropic site (the bromine atom) and three kosmotropic sites (three oxygen atoms) that are present in the ion.

Interaction between Water and Alkali Metal Ions and Its Temperature Dependence Revealed by Oxygen K-Edge X-ray Absorption Spectroscopy

Interaction between water molecules and alkali metal ions in aqueous salt solutions has been studied by the oxygen K-edge soft X-ray absorption spectroscopy (XAS) in transmission mode. In the measurement of several alkali halide aqueous solutions with different alkali chlorides (Li, Na, and K) and different sodium halides (Cl, Br, and I), the pre-edge component arising from the hydration water molecules shows a blue shift in peak energy as strongly depending on cations but not on anions. In the temperature dependent measurement, the pre-edge component arising from water molecules beyond the first hydration shell shows the same behavior as that of pure liquid water. On the other hand, the pre-edge component arising from water molecules in the first hydration shell of Li⁺ ions is not evidently dependent on the temperature, indicating that the hydration water molecules are more strongly bound with Li⁺ ions than the other water molecules. These experimental results are supported by the results of radial distribution functions of the first hydration shell and their temperature dependence, evaluated by molecular dynamics simulations.