Water in Hydration Shell of an Iodide Ion: Structure and Dynamics of Solute-Water Hydrogen Bonds and Vibrational Spectral Diffusion from First-Principles Simulations

All-Round Ionic Liquids for Shuttle-Free Zinc-Iodine Battery

The practical implementation of aqueous zinc-iodine batteries (ZIBs) is hindered by the rampant Zn dendrites growth, parasite corrosion, and polyiodide shuttling. In this work, ionic liquid EMIM[OAc] is employed as an all-round solution to mitigate challenges on both the Zn anode and the iodine cathode side. First, the EMIM+ embedded lean-water inner Helmholtz plane (IHP) and inert solvation sheath modulated by OAc- effectively repels H2 O molecules away from the Zn anode surface. The preferential adsorption of EMIM+ on Zn metal facilitates uniform Zn nucleation via a steric hindrance effect. Second, EMIM+ can reduce the polyiodide shuttling by hindering the iodine dissolution and forming an EMIM+ -I3 - dominated phase. These effects holistically enhance the cycle life, which is manifested by both Zn || Zn symmetric cells and Zn-I2 full cells. ZIBs with EAc deliver a capacity decay rate of merely 0.01 ‰ per cycle after over 18,000 cycles at 4 A g-1 , and lower self-discharge and better calendar life than the ZIBs without ionic liquid EAc additive.

Role of Surface Features on the Initial Dissolution of CH3NH3PbI3 Perovskite in Liquid Water: An Ab Initio Molecular Dynamics Study

The degradation of CH3NH3PbI3 (MAPbI3) hybrid organic inorganic perovskite (HOIP) by water has been the major issue hampering its use in commercial perovskites solar cells (PSCs), as MAPbI3 HOIP has been known to easily degrade in the presence of water. Even though there have been numerous studies investigating this phenomenon, there is still no consensus on the mechanisms of the initial stages of dissolution. Here, we attempt to consolidate differing mechanistic interpretations previously reported in the literature through the use of the first-principles constrained ab initio molecular dynamics (AIMD) to study both the energetics and mechanisms that accompany the degradation of MAPbI3 HOIP in liquid water. By comparing the dissolution free energy barrier between surface species of different surficial types, we find that the dominant dissolution mechanisms of surface species vary widely based on the specific surface features. The high sensitivity of the dissolution mechanism to surface features has contributed to the many dissolution mechanisms proposed in the literature. In contrast, the dissolution free energy barriers are mainly determined by the dissolving species rather than the type of surfaces, and the type of surfaces the ions are dissolving from is inconsequential toward the dissolution free energy barrier. However, the presence of surface defects such as vacancy sites is found to significantly lower the dissolution free energy barriers. Based on the estimated dissolution free energy barriers, we propose that the dissolution of MAPbI3 HOIP in liquid water originates from surface defect sites that propagate laterally along the surface layer of the MAPbI3 HOIP crystal.

Microwave Dielectric Relaxation of Univalent and Bivalent Electrolyte Solutions

The microwave dielectric relaxation of aqueous solutions of univalent (KCl, NaCl, NaI) and bivalent (CaCl2, MgCl2) electrolytes at concentrations between 0.1 and 1 M at 25 °C was investigated using a vector network analyzer (0.5≤ ν ≤ 40 GHz). The spectra of these electrolyte systems are characterized by a symmetrical broadening of the main relaxation peak and were fitted using the Cole-Cole equation. In our analysis, we provide insights into the underlying physics of the relaxation events at microscopic and mesoscopic scales by using a 3D phase space trajectory that is based on the interactions of the relaxing dipole units with their surroundings and Frohlich's B function. The effect of the solutes on the H-bond network of water with increasing concentration is evident in the microwave dielectric spectra through decreasing dielectric strengths and relaxation times. It was found that the number of perturbed water molecules is higher in the case of bivalent electrolytes and appears to be proportional to the ionic radius. In our approach, the particular dependence between the broadening parameter α and the relaxation times τ reflects the rate of interactions between the elementary dipole units and their surroundings. We provide a quantitative analysis of the level of perturbation caused by the presence of ions in the hydrogen-bond network of water. It was found that the H-bonded network of water is highly perturbed in univalent systems compared to bivalent systems due to weaker bonded hydration shells. Finally, we found significant differences between the dielectric response of NaCl and NaI. The differences, originating in the counterions Cl- and I-, which are characterized by large ionic radii and consequently weaker electric fields in their vicinity, confirm that the effect of weakly hydrated ions should not be neglected in microwave dielectric spectra analysis.

Solvent Dynamics of Aqueous Halides before and after Photoionization

Electron transfer reactions can be strongly influenced by solvent dynamics. We study the photoionization of halides in water as a model system for such reactions. There are no internal nuclear degrees of freedom in the solute, allowing the dynamics of the solvent to be uniquely identified. We simulate the equilibrium solvent dynamics for Cl^-, Br^-, I^-, and their respective neutral atoms in water, comparing quantum mechanical/molecular mechanical (QM/MM) and classical molecular dynamics (MD) methods. On the basis of the obtained configurations, we calculate the extended X-ray absorption fine structure (EXAFS) spectra rigorously based on the MD snapshots and compare them in detail with other theoretical and experimental results available in the literature. We find our EXAFS spectra based on QM/MM MD simulations in good agreement with their experimental counterparts for the ions. Classical MD simulations for the ions lead to EXAFS spectra that agree equally well with the experiment when it comes to the oscillatory period of the signal, even though they differ from the QM/MM radial distribution functions extracted from the MD. The amplitude is, however, considerably overestimated. This suggests that to judge the reliability of theoretical simulation methods or to elucidate fine details of the atomistic dynamics of the solvent based on EXAFS spectra, the amplitude as well as the oscillatory period need to be considered. If simulations fail qualitatively, as does the classical MD for the aqueous neutral halogen atoms, the resulting EXAFS will also be strongly affected in both oscillatory period and amplitude. The good reliability of QM/MM-based EXAFS simulations, together with clear qualitative differences in the EXAFS spectra found between halides and their atomic counterparts, suggests that a combined theory and experimental EXAFS approach is suitable for elucidating the nonequilibrium solvent dynamics in the photoionization of halides and possibly also for electron transfer reactions in more complex systems.

Systematic Evaluation of Ion Diffusion and Water Exchange

As a fundamental property of all fluids, diffusion plays myriad roles in both science and our daily lives. Diffusive properties of many liquids including water have been extensively studied both experimentally and theoretically, while for transition metal ions, there exist significant experimental data that have not been extensively studied theoretically. Hence, high-confidence predictions for challenging systems like radioactive ions that are biohazardous cannot be reliably generated. In this work, a workflow named ISAIAH (Ion Simulation using AMBER for dIffusion Action when Hydrated) was designed to accurately simulate the diffusion coefficients of 15 monoatomic ions with charges varying from -1 to +3 in four water models. As the results indicate, good agreement with experimental values was achieved, leading us to select ²³⁹Pu⁴⁺ (for which no experimental data are available) as a candidate ion to make a theoretical prediction of its diffusion coefficient in water. Among all the force field parameter sets, the ones parametrized using an augmented 12-6-4 Lennard-Jones (LJ) potential showed lower average unsigned errors (AUE) for ions of various radii and electron configurations relative to some 12-6 LJ parameters. This observation agrees well with the fact that diffusion is affected by both the hydration free energy (HFE) and the ion-oxygen distance (IOD) between solute and solvent molecules, both of which are handled well by the 12-6-4 model.

Angstrom-scale ion channels towards single-ion selectivity

Artificial ion channels with ion permeability and selectivity comparable to their biological counterparts are highly desired for efficient separation, biosensing, and energy conversion technologies. In the past two decades, both nanoscale and sub-nanoscale ion channels have been successfully fabricated to mimic biological ion channels. Although nanoscale ion channels have achieved intelligent gating and rectification properties, they cannot realize high ion selectivity, especially single-ion selectivity. Artificial angstrom-sized ion channels with narrow pore sizes <1 nm and well-defined pore structures mimicking biological channels have accomplished high ion conductivity and single-ion selectivity. This review comprehensively summarizes the research progress in the rational design and synthesis of artificial subnanometer-sized ion channels with zero-dimensional to three-dimensional pore structures. Then we discuss cation/anion, mono-/di-valent cation, mono-/di-valent anion, and single-ion selectivities of the synthetic ion channels and highlight their potential applications in high-efficiency ion separation, energy conversion, and biological therapeutics. The gaps of single-ion selectivity between artificial and natural channels and the connections between ion selectivity and permeability of synthetic ion channels are covered. Finally, the challenges that need to be addressed in this research field and the perspective of angstrom-scale ion channels are discussed.

Molecular Dynamics Insight into the Role of Water Molecules in Ionic Liquid Mixtures of 1-Butyl-3-methylimidazolium Iodide and Ethylammonium Nitrate

Imidazolium-based ionic liquids are well known for their versatility as solvents for various applications such as dye-sensitized solar cells, fuel cells, and lithium-ion batteries; however, their complex interactions continue to be investigated to further improve upon their design. Ionic liquids (ILs) are commonly mixed with co-solvents such as water, organic solvents, or other ionic liquids to tailor their physiochemical properties. To better predict these properties and fundamentally understand the molecular interactions within the electrolyte mixtures, molecular dynamics (MD) simulations are often employed. In this study, MD simulations are performed on ternary solutions containing ionic liquids of 1-butyl-3-methylimidazolium iodide ([BMIM][I]) and ethylammonium nitrate ([EA][NO₃]) with increasing concentration of water. As previously reported, these ternary solutions displayed a wide temperature window of thermal stability and electrochemical conductivity. Utilizing MD simulations, the complex intermolecular interactions are identified, and the role of water as a co-solvent is disclosed to correlate with changes in their bulk properties. The MD results, including simulation box snapshots, radial distribution functions, and self-diffusion coefficients, reveal the formation of heterogeneous regimes with increasing water concentration, hydrogen bonding between iodide-water, iodide-[EA]⁺, and a change in IL ordering when in mixtures containing water. The simulations also display the formation of water aggregates and networks at high water concentrations, which can contribute to the thermal behavior of the respective mixtures. As the design of IL-based electrolytes grows in demand with increasing complexity, this work demonstrates the capability of MD simulations containing multiple constituents and their necessity in material development through identification of microscopic structure-property relationships.

Corrections in the CHARMM36 Parametrization of Chloride Interactions with Proteins, Lipids, and Alkali Cations, and Extension to Other Halide Anions

The nonpolarizable CHARMM force field is one of the most widely used energy functions for all-atom biomolecular simulations. Chloride is the only halide ion included in the latest version, CHARMM36m, and is used widely in simulation studies, often as an electrolyte ion but also as the biological substrate of transport proteins and enzymes. Here, we find that existing parameters systematically underestimate the interaction of Cl^- with proteins and lipids. Accordingly, when examined in solution, little to no Cl^-association can be observed with most components of the protein, including backbone, polar side chains and aromatic rings. The strength of the interaction with cationic side chains and with alkali ions is also incongruent with experimental measurements, specifically osmotic coefficients of concentrated solutions. Consistent with these findings, a 4-μs trajectory of the Cl^--specific transport protein CLC-ec1 shows irreversible Cl^- dissociation from the so-called S_cen binding site, even in a 150 mM NaCl buffer. To correct for these deficiencies, we formulate a series of pair-specific Lennard-Jones parameters that override those resulting from the conventional Lorentz-Berthelot combination rules. These parameters, referred to as NBFIX, are systematically calibrated against available experimental data as well as ab initio geometry optimizations and energy evaluations, for a wide set of binary and ternary Cl^- complexes with protein and lipid analogs and alkali cations. Analogously, we also formulate parameter sets for the other three biological halide ions, namely, fluoride, bromide, and iodide. The resulting parameters are used to calculate the potential of mean force defining the interaction of each anion and each of the protein and lipid analogues in bulk water, revealing association free energies in the range of -0.3 to -3.3 kcal/mol, with the F^- complexes being the least stable. The NBFIX corrections also preserve the Cl^- occupancy of CLC-ec1 in a second 4-μs trajectory. We posit that these optimized molecular-mechanics models provide a more realistic foundation for all-atom simulation studies of processes entailing changes in hydration, recognition, or transport of halide anions.

Recent Advances in Luminescent Recognition and Chemosensing of Iodide in Water

This Minireview covers the latest developments of chemosensors based on transition-metal receptors and organic fluorophores with specific binding sites for the luminescent detection and recognition of iodide in aqueous media and real samples. In all selected examples within the last decade (made-post 2010), the iodide sensing and recognition is probed by monitoring real-time changes of the fluorescence or phosphorescence properties of the chemosensors. This review highlights effective strategies to iodide sensing from a structural approach where the iodide recognition/sensing process, through supramolecular interactions as coordination bonds, hydrogen bonds, halogen bonds and electrostatic interactions, is transduced into an optical change easily measurable. The selective iodide sensing is an active field of research with global interest due to the importance of iodide in biological, medicinal, industrial, environmental and chemical processes.

Quadrupolar NMR Relaxation of Aqueous 127I-, 131Xe, and 133Cs+: A First-Principles Approach from Dynamics to Properties

Quadrupolar NMR relaxation rates were computed for aqueous ¹³³Cs⁺, ¹³¹Xe, and ¹²⁷I^- via Kohn-Sham (KS) density functional theory-based ab initio molecular dynamics and KS calculations of the electric field gradient (EFG) tensors along the trajectories. The resulting rates are within a factor of 1-3 of the experimental values and can be compared to available results from classical dynamics and EFGs from electrostatic models with corrections via Sternheimer antishielding factors. Relativistic effects are shown to have an enhancing effect on the magnitude of the EFGs. An analysis of the EFGs was carried out in terms of localized molecular orbitals to elucidate contributions from the solvent versus solute polarization and assess the validity of the Sternheimer approximation for these systems.

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.

New Insights into the Charge-Transfer-to-Solvent Spectrum of Aqueous Iodide: Surface versus Bulk

Liquid phase charge-transfer-to-solvent (CTTS) transitions are important, as they serve as photochemical routes to solvated electrons. In this work, broadband deep-ultraviolet electronic sum frequency generation (DUV-ESFG) and two-photon absorption (2PA) spectroscopic techniques were used to assign and compare the nature of the aqueous iodide CTTS excitations at the air/water interface and in bulk solution. In the one-photon absorption (1PA) spectrum, excitation to the 6s Rydberg-like orbital (5p → 6s) gives rise to a pair of spin-orbit split iodine states, ²P_3/2 and ²P_1/2. In the 2PA spectra, the lower-energy ²P_3/2 peak is absent and the observed 2PA peak, which is ∼0.14 eV blue-shifted relative to the upper ²P_1/2 CTTS peak seen in 1PA, arises from 5p → 6p electronic promotion. The band observed in the ESFG spectrum is attributed to mixing of excited states involving 5p → 6p and 5p → 6s promotions caused by both vibronic coupling and the external electric field generated by asymmetric interfacial solvation.

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.

Time-dependent vibrational spectral analysis of first principles trajectory of methylamine with wavelet transform

The fluctuation dynamics of amine stretching frequencies, hydrogen bonds, dangling N-D bonds, and the orientation profile of the amine group of methylamine (MA) were investigated under ambient conditions by means of dispersion-corrected density functional theory-based first principles molecular dynamics (FPMD) simulations. Along with the dynamical properties, various equilibrium properties such as radial distribution function, spatial distribution function, combined radial and angular distribution functions and hydrogen bonding were also calculated. The instantaneous stretching frequencies of amine groups were obtained by wavelet transform of the trajectory obtained from FPMD simulations. The frequency-structure correlation reveals that the amine stretching frequency is weakly correlated with the nearest nitrogen-deuterium distance. The frequency-frequency correlation function has a short time scale of around 110 fs and a longer time scale of about 1.15 ps. It was found that the short time scale originates from the underdamped motion of intact hydrogen bonds of MA pairs. However, the long time scale of the vibrational spectral diffusion of N-D modes is determined by the overall dynamics of hydrogen bonds as well as the dangling ND groups and the inertial rotation of the amine group of the molecule.

Origin of the blueshift of water molecules at interfaces of hydrophilic cyclic compounds

Water molecules at interfaces of materials exhibit enigmatic properties. A variety of spectroscopic studies have observed a high-frequency motion in these water molecules, represented by a blueshift, at both hydrophobic and hydrophilic interfaces. However, the molecular mechanism behind this blueshift has remained unclear. Using Raman spectroscopy and ab initio molecular dynamics simulations, we reveal the molecular mechanism of the blueshift of water molecules around six monosaccharide isomers. In the first hydration shell, we found weak hydrogen-bonded water molecules that cannot have a stable tetrahedral water network. In the water molecules, the vibrational state of the OH bond oriented toward the bulk solvent strongly contributes to the observed blueshift. Our work suggests that the blueshift in various solutions originates from the vibrational motions of these observed water molecules.

Structural, dynamical, and transport properties of the hydrated halides: How do At(-) bulk properties compare with those of the other halides, from F(-) to I(-)?

The properties of halides from the lightest, fluoride (F(-)), to the heaviest, astatide (At(-)), have been studied in water using a polarizable force-field approach based on molecular dynamics (MD) simulations at the 10 ns scale. The selected force-field explicitly treats the cooperativity within the halide-water hydrogen bond networks. The force-field parameters have been adjusted to ab initio data on anion/water clusters computed at the relativistic Möller-Plesset second-order perturbation theory level of theory. The anion static polarizabilities of the two heaviest halides, I(-) and At(-), were computed in the gas phase using large and diffuse atomic basis sets, and taking into account both electron correlation and spin-orbit coupling within a four-component framework. Our MD simulation results show the solvation properties of I(-) and At(-) in aqueous phase to be very close. For instance, their first hydration shells are structured and encompass 9.2 and 9.1 water molecules at about 3.70 ± 0.05 Å, respectively. These values have to be compared to the F(-), Cl(-), and Br(-) ones, i.e., 6.3, 8.4, and 9.0 water molecules at 2.74, 3.38, and 3.55 Å, respectively. Moreover our computations predict the solvation free energy of At(-) in liquid water at ambient conditions to be 68 kcal mol(-1), a value also close the I(-) one, about 70 kcal mol(-1). In all, our simulation results for I(-) are in excellent agreement with the latest neutron- and X-ray diffraction studies. Those for the At(-) ion are predictive, as no theoretical or experimental data are available to date.