Ion solvation and water structure in potassium halide aqueous solutions

Exploring the influence of metal cations on individual hydrogen bonds in Watson-Crick guanine-cytosine DNA base pair: An interacting quantum atoms analysis

This study delves into the nature of individual hydrogen bonds and the relationship between metal cations and hydrogen bonding in the Watson-Crick guanine-cytosine (GC) base pair and its alkali and alkaline earth cation-containing complexes (Mn+-GC). The findings reveal how metal cations affect the nature and strength of individual hydrogen bonds. The study employs interacting quantum atoms (IQA) analysis to comprehensively understand three individual hydrogen bonds within the GC base pair and its cationic derivatives. These analyses unveil the nature and strength of hydrogen bonds and serve as a valuable reference for exploring the impact of cations (and other factors) on each hydrogen bond. All the H⋯ $$ \cdots $$ D interactions (H is hydrogen and D is oxygen or nitrogen) in the GC base pair are primarily electrostatic in nature, with the charge transfer component playing a substantial role. Introducing a metal cation perturbs all H⋯ $$ \cdots $$ D interatomic interactions in the system, weakening the nearest hydrogen bond to the cation (indicated by a) and reinforcing the other (b and c) interactions. Notably, the interaction a, the strongest H⋯ $$ \cdots $$ D interaction in the GC base pair, becomes the weakest in the Mn+-GC complexes. A broader perspective on the stability of GC and Mn+-GC complexes is provided through interacting quantum fragments (IQF) analysis. This approach considers all pairwise interactions between fragments and intra-fragment components, offering a complete view of the factors that stabilize and destabilize GC and Mn+-GC complexes. The IQF analysis underscores the importance of electron sharing, with the dominant contribution arising from the inter-fragment exchange-correlation term, in shaping and sustaining GC and Mn+-GC complexes. From this point of view, alkaline and alkaline earth cations have distinct effects, with alkaline cations generally weakening inter-fragment interactions and alkaline earth cations strengthening them. In addition, IQA and IQF calculations demonstrate that the hydration of cations led to small changes in the hydrogen bonding network. Finally, the IQA interatomic energies associated with the hydrogen bonds and also inter-fragment interaction energies provide robust indicators for characterizing hydrogen bonds and complex stability, showing a strong correlation with total interaction energies.

To Pair or not to Pair? Machine-Learned Explicitly-Correlated Electronic Structure for NaCl in Water

The extent of ion pairing in solution is an important phenomenon to rationalize transport and thermodynamic properties of electrolytes. A fundamental measure of this pairing is the potential of mean force (PMF) between solvated ions. The relative stabilities of the paired and solvent shared states in the PMF and the barrier between them are highly sensitive to the underlying potential energy surface. However, direct application of accurate electronic structure methods is challenging, since long simulations are required. We develop wave function based machine learning potentials with the random phase approximation (RPA) and second order Møller-Plesset (MP2) perturbation theory for the prototypical system of Na and Cl ions in water. We show both methods in agreement, predicting the paired and solvent shared states to have similar energies (within 0.2 kcal/mol). We also provide the same benchmarks for different DFT functionals as well as insight into the PMF based on simple analyses of the interactions in the system.

Rocking-Chair Aqueous Fluoride-Ion Batteries Enabled by Hydrogen Bonding Competition

Aqueous fluoride ion batteries (FIBs) have garnered attention for their high theoretical energy density, yet they are challenged by sluggish fluorination kinetics, active material dissolution, and electrolyte instability. Here, we present a room temperature rocking-chair aqueous FIBs featuring KOAc-KF binary salt electrolytes, enabling concurrent fluorination and defluorination reactions at both cathode and anode electrodes. Experimental and theoretical results reveal that acetate ions in the electrolyte compete with fluoride ions in hydrogen bonding formation, weakening the excessively strong solvation between H2O and F- ions. This results in the suppression of detrimental HF formation and a reduced desolvation energy of F- ions, enhancing the electrochemical reaction kinetics. The bismuth-based cathode exhibits direct conversion in the optimized electrolyte, effectively suppressing the detrimental disproportionation reactions from Bi2+ intermediates. Additionally, zinc anode undergoes a typical fluorination process, forming solid KZnF3 as the electrode product, minimizing the risks of hydrogen evolution. The proposed aqueous FIBs with the optimized electrolyte demonstrate high discharge capacity, long-term cycling stability and excellent rate capabilities.

Characterization of the Coordination and Solvation Dynamics of Solvated Systems─Implications for the Analysis of Molecular Interactions in Solutions and Pure H2O

The characterization of solvation shells of atoms, ions, and molecules in solution is essential to relate solvation properties to chemical phenomena such as complex formation and reactivity. Different definitions of the first-shell coordination sphere from simulation data can lead to potentially conflicting data on the structural properties and associated ligand exchange dynamics. The definition of a solvation shell is typically based on a given threshold distance determined from the respective solute-solvent pair distribution function g(r) (i.e., GC). Alternatively, a nearest neighbor (NN) assignment based on geometric properties of the coordination complex without the need for a predetermined cutoff criterion, such as the relative angular distance (RAD) or the modified Voronoi (MV) tessellation, can be applied. In this study, the effect of different NN algorithms on the coordination number and ligand exchange dynamics evaluated for a series of monatomic ions in aqueous solution, carbon dioxide in aqueous and dichloromethane solutions, and pure liquid water has been investigated. In the case of the monatomic ions, the RAD approach is superior in achieving a well separated definition of the first solvation layer. In contrast, the MV algorithm provides a better separation of the NNs from a molecular point of view, leading to better results in the case of solvated CO2. When analyzing the coordination environment in pure water, the cutoff-based GC framework was found to be the most reliable approach. By comparison of the number of ligand exchange reactions and the associated mean ligand residence times (MRTs) with the properties of the coordination number autocorrelation functions, it is shown that although the average coordination numbers are sensitive to the different definitions of the first solvation shell, highly consistent estimates for the associated MRT of the solvated system are obtained in the majority of cases.

Going beyond Radial Hydration Models: The Hidden Structures of Chloride and Iodide Aqua Ions Revealed by the Use of Lone Pairs

A novel model of hydration for the chloride and iodide ions in water is proposed, which overcomes the limitations of conventional radial models. A new approach, based on a representation of the halide lone pairs, highlighted a subset of first shell water molecules featuring preferential strong interactions with the ion lone pairs, giving rise to tetrahedral hydration structures in both Cl- and I- aqueous solutions. By adopting a novel descriptor correlated to the halide-water interaction energy, we were able to split the conventional first solvation shell into a tight first hydration shell, composed of water molecules strongly interacting with the ions via hydrogen bonds, and a loose first shell containing molecules that are only slightly perturbed by the halide electrostatic charge. The picture emerging from our findings indicates that lone pairs play an important role in the description of systems where hydrogen bonds are the main interactions taking place in the solvation process.

Structural and dynamical properties of concentrated alkali- and alkaline-earth metal chloride aqueous solutions

Concentrated ionic aqueous electrolytes possess a diverse array of applications across various fields, particularly in the field of energy storage. Despite extensive examination, the intricate relationships and numerous physical mechanisms underpinning diverse phenomena remain incompletely understood. Molecular dynamics simulations are employed to probe the attributes of aqueous solutions containing LiCl, NaCl, KCl, MgCl2, and CaCl2, spanning various solute fractions. The primary emphasis of the simulations is on unraveling the intricate interplay between these attributes and the underlying physical mechanisms. The configurations of cation-Cl- and Cl--Cl- pairs within these solutions are disclosed. As the solute fraction increases, consistent trends manifest regardless of solute type: (i) the number of hydrogen bonds formed by the hydration water surrounding ions decreases, primarily attributed to the growing presence of counter ions in proximity to the hydration water; (ii) the hydration number of ions exhibits varying trends influenced by multiple factor; and (iii) the diffusion of ions slows down, attributed to the enhanced confinement and rebound of cations and Cl- ions from the surrounding atoms, concurrently coupled with the changes in ion vibration modes. In our analysis, we have, for the first time, clarified the reasons behind the slowing down of the diffusion of the ions with increasing solute fraction. Our research contributes to a better understanding and manipulation of the attributes of ionic aqueous solutions and may help designing high-performance electrolytes.

Novel Computational Chemistry Infrastructure for Simulating Astatide in Water: From Basis Sets to Force Fields Using Particle Swarm Optimization

Using the example of astatine, the heaviest naturally occurring halogen whose isotope At-211 has promising medical applications, we propose a new infrastructure for large-scale computational models of heavy elements with strong relativistic effects. In particular, we focus on developing an accurate force field for At- in water based on reliable relativistic density functional theory (DFT) calculations. To ensure the reliability of such calculations, we design novel basis sets for relativistic DFT, via the particle swarm optimization algorithm to optimize the coefficients of the new basis sets and the polarization-consistent basis set idea's extension to heavy elements to eliminate the basis set error from DFT calculations. The resulting basis sets enable the well-grounded evaluation of relativistic DFT against "gold-standard" CCSD(T) results. Accounting for strong relativistic effects, including spin-orbit interaction, via our redesigned infrastructure, we elucidate a noticeable dissimilarity between At- and I- in halide-water force field parameters, radial distribution functions, diffusion coefficients, and hydration energies. This work establishes the framework for the systematic development of polarization-consistent basis sets for relativistic DFT and accurate force fields for molecular dynamics simulations to be used in large-scale models of complex molecular systems with elements from the bottom of the periodic table, including actinides and even superheavy elements.

Visualization of Self-Assembly and Hydration of a β-Hairpin through Integrated Small and Wide-Angle Neutron Scattering

Fundamental understanding of the structure and assembly of nanoscale building blocks is crucial for the development of novel biomaterials with defined architectures and function. However, accessing self-consistent structural information across multiple length scales is challenging. This limits opportunities to exploit atomic scale interactions to achieve emergent macroscale properties. In this work we present an integrative small- and wide-angle neutron scattering approach coupled with computational modeling to reveal the multiscale structure of hierarchically self-assembled β hairpins in aqueous solution across 4 orders of magnitude in length scale from 0.1 Å to 300 nm. Our results demonstrate the power of this self-consistent cross-length scale approach and allows us to model both the large-scale self-assembly and small-scale hairpin hydration of the model β hairpin CLN025. Using this combination of techniques, we map the hydrophobic/hydrophilic character of this model self-assembled biomolecular surface with atomic resolution. These results have important implications for the multiscale investigation of aqueous peptides and proteins, for the prediction of ligand binding and molecular associations for drug design, and for understanding the self-assembly of peptides and proteins for functional biomaterials.

Terahertz Spectroscopy of Aqueous Solutions of Sodium Halides (NaX): Self- and Cross-Correlation Contributions of Ions and Hydration Shell Water for X- = F-, Cl-, Br-, and I. J Phys Chem B 2023;127:9323-9335. [PMID: 37871257 DOI: 10.1021/acs.jpcb.3c05228] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

We investigated the terahertz (THz) absorption spectra of aqueous sodium halide solutions through molecular dynamics simulations using polarizable models of both water and ions. Specifically, we have considered aqueous solutions (∼1 M) of NaF, NaCl, NaBr, and NaI and calculated the difference THz spectrum of these solutions by subtracting the corresponding pure water contribution from the total THz spectrum of an ionic solution. The difference absorption spectrum of a given solution is then dissected into contributions from ion and ion-water correlations and also modifications of water-water correlations in the presence of the ions. The different components are further decomposed into induced dipole and permanent charge/dipole components and also into self- and cross-correlation components. The ion-water cross-correlation components are subsequently decomposed into contributions coming from different solvation shells through radially resolved calculations of such ion-water cross-correlations. Through all of these dissections, we could investigate the origin of different parts of the difference THz spectra of the sodium halide solutions studied here. It is found that while features below or around 100 cm-1 and also around 200 cm-1 arise mainly from ion and ion-water motion, that at the librational region above 600 cm-1 primarily originates from changes in water librational motion influenced by the ions. The variations of intensities of different components are also linked to the size and charge density of the anions in the solutions.

Ab Initio Molecular Dynamics Investigation of the Solvation States of Hydrated Ions in Confined Water

Ionic transport in nanoscale channels with a critical size comparable to that of ions and solutes exhibits exceptional performance in water desalination, ion separation, electrocatalysts, and supercapacitors. However, the solvation states (SSs), i.e., the hydration structures and probability distribution, of hydrated ions in nanochannels differ from those in the bulk and the perspective of continuum theory. In this work, we conduct ab initio enhanced-sampling atomistic simulations to investigate the ion-specific SSs of monovalent ions (including Li+, Na+, K+, F-, Cl-, and I-) in the graphene channel with a width of 1 nm. Our findings highlight that the SSs of those ions are primarily determined by ion-water hydration, where ion-wall interactions play a minor role. The distribution of ions in layered confined water is a result of ion-specific hydration, which arises from the synergy of entropy and enthalpy. The free energy barriers for transitions between SSs are on the order of 1kBT, allowing for modulation through applying external fields or modifying surface properties. As the ion-wall interaction strengthens, as observed in vermiculite and carbides and nitrides of transition metal channels, the probability of near-wall SSs increases. These results help to improve the performance of nanofluidic devices and provide crucial insights for developing accurate force fields of molecular simulations or advanced theoretical approaches for ion dynamics in confined channels.

Cation coordination polyhedra lead to multiple lengthscale organization in aqueous electrolytes

Understanding multiple lengthscale correlations in the pair distribution functions (PDFs) of aq. electrolytes is a persistent challenge. Here, the coordination chemistry of polyoxoanions supports an ion-network of cation-coordination polyhedra in NaNO3(aq) and NaNO2(aq) that induce long-range solution structure. Oxygen correlations associated with Na+-coordination polyhedra have two characteristics lengthscales; 3.5-5.5 Å and 5.5-7.5 Å, the latter solely associated oligomers. The PDF contraction between 5.5-7.5 Å observed in many electrolytes is attributed to the distinct O⋯O correlation found in dimers and dimer subunits within oligomers.

Observation of Solid-Liquid Phase Transitions of Brine Using the Far-Ultraviolet Charge-Transfer-to-Solvent Band

Although determining the chemical states of salts and ions is critical in numerous fields, such as elucidating biological functions and maintaining food quality, the current direct observation methods are insufficient. We propose a spectral analysis method of directly observing the phase transitions of NaCl solutions using the changes in the charge-transfer-to-solvent band and the absorption band representing the first electron transition (Ã ← X̃) of H₂O. The intensities of these bands may be observed using attenuated total reflection far-ultraviolet spectroscopy. According to the well-known phase diagram of aqueous NaCl, we observe spectral changes during freezing-thawing and may spectroscopically detect the phase transitions from liquid to mixed liquid-solid and solid phases, including eutectic crystals, in addition to their coexistence curves.

Ion/Water Network Structural Dynamics in Highly Concentrated Lithium Chloride and Lithium Bromide Solutions Probed with Ultrafast Infrared Spectroscopy

The structural dynamics of highly concentrated LiCl and LiBr aqueous solutions were observed from 1-4 to 1-16 water molecules per ion pair using ultrafast polarization-selective pump-probe (PSPP) experiments on the OD stretch of dilute HOD. At these high salt concentrations, an extended ion/water network exists with complex structural dynamics. Population decays from PSPP experiments highlight two distinct water components. From the frequency-dependent amplitudes of the decays, the spectra of hydroxyls bound to halides and to water oxygens are obtained, which are not observable in the FT-IR spectra. PSPP experiments also measure frequency-dependent water orientational relaxation. At short times, wobbling dynamics within a restricted angular cone occurs. At high concentrations, the cone angles are dependent on frequency (hydrogen bond strength), but at higher water concentrations (>10 waters per ion pair), there is no frequency dependence. The average cone angle increases as the ion concentration decreases. The slow time constant for complete HOD orientational relaxation is independent of concentration but slower in LiCl than in LiBr. Comparison to structural MD simulations of LiCl from the literature indicates that the loss of the cone angle wavelength dependence and the increase in the cone angles as the concentration decreases occur as the prevalence of large ion/water clusters gives way to contact ion pairs.

Mechanism of Cations Suppressing Proton Diffusion Kinetics for Electrocatalysis

Proton transfer is crucial for electrocatalysis. Accumulating cations at electrochemical interfaces can alter the proton transfer rate and then tune electrocatalytic performance. However, the mechanism for regulating proton transfer remains ambiguous. Here, we quantify the cation effect on proton diffusion in solution by hydrogen evolution on microelectrodes, revealing the rate can be suppressed by more than 10 times. Different from the prevalent opinions that proton transport is slowed down by modified electric field, we found water structure imposes a more evident effect on kinetics. FTIR test and path integral molecular dynamics simulation indicate that proton prefers to wander within the hydration shell of cations rather than to hop rapidly along water wires. Low connectivity of water networks disrupted by cations corrupts the fast-moving path in bulk water. This study highlights the promising way for regulating proton kinetics via a modified water structure.

Molecular Influences on the Quantification of Lewis Acidity with Phosphine Oxide Probes

Gutmann-Beckett-type measurements with phosphine oxide probes can be used to estimate effective Lewis acidity with ³¹P nuclear magnetic resonance spectroscopy, but the influence of the molecular structure of a given probe on the quantification of Lewis acidity remains poorly documented in experimental work. Here, a quantitative comparison of triethyl (E), trioctyl (O), and triphenyl (P) phosphine oxides as molecular probes of Lewis acidity has been carried out via titration studies in MeCN with a test set of six mono- and divalent metal triflate salts. In comparison to E, the bulkier O displays a similar range of chemical shift values and binding affinities for the various test metal ions. Spectral linewidths and speciation properties vary for individual cation-to-probe ratios, however, confirming probe-specific properties that can impact the data quality. Importantly, P displays a consistently narrower dynamic range than both E and O, illustrating how electronic changes at phosphorus can influence the NMR response. Comparative parametrizations of the effective Lewis acidities of a broader range of metal ions, including the trivalent rare earth ions Y³⁺, Lu³⁺, and Sc³⁺ as well as the uranyl ion (UO₂²⁺), can be understood in light of these results, providing insight into the fundamental chemical processes underlying the useful approach of single-point measurements for quantification of effective Lewis acidity. Together with a study of counteranion effects reported here, these data clarify the diverse ensemble of factors that can influence the measurement of Lewis acid/base interactions.

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.

Dynamic Electric Field Alignment Determines the Water Rotational Motion around Protein

Water rotational dynamics in biomolecular solution is crucial to evaluating and controlling biomolecule stability. In this molecular dynamics simulation (MD) study on lysozyme solutions, we present how the exerted internal electric field determines water rotational dynamics. We find that the relaxation time of water rotation is equivalent to that of the reorientation of the exerted overall electric field for every single water molecule, regardless of its translation mode. Namely, water molecular rotation synchronizes with the exerted field reorientation. We also map the reorientation process of the electric field at fixed points relative to protein in the solution, which displays the local hydration dynamics commensurate with the reported time-dependent fluorescence Stokes shift (TDFSS) measurements. Comparing the spatial distribution of local field reorientation relaxation time with that of rotational relaxation time, we further suggest that water rotation dynamics are subject to the reorientation of the local overall field within the hydration layer. While outside the hydration layer, the relaxation time of the local electric field reorientation is short enough (subpicosecond) to assume the δ function, showing the electric force with randomly changing orientation is applied to each water molecule.

Alchemical Osmostat for Monte Carlo Simulation: Sampling Aqueous Electrolyte Solution in Open Systems

Molecular simulations involving electrolytes are usually performed at a fixed amount of salt ions in the simulation box, reproducing macroscopic concentration. Although this statement is valid in the bulk, the concentration of an electrolyte confined in nanoporous materials such as MOFs or zeolites is greatly affected and remains a priori unknown. The nanoporous material in equilibrium with the bulk electrolyte exchange water and ions at a given chemical potential Δμ in the semi-grand-canonical ensemble, that must be calibrated in order to determine the concentration in the nanoporous material. In this work, we propose an algorithm based on nonequilibrium candidate Monte Carlo (NCMC) moves to ultimately perform MC simulations in contact with a saline reservoir. First, we adapt the Widom insertion technique to calibrate the chemical potential by alchemically transmuting water molecules into ions by using NCMC moves. The chemical potential defines a Monte Carlo osmostat in the semi-grand-constant volume and temperature ensemble (Δμ, N, V, T) to be added in a Monte Carlo simulation where the number of ions fluctuates. In order to validate the method, we adapted the NCMC move to determine the free energy of water solvation and subsequently explore thermodynamics of electrolyte solvation at infinite dilution in water. Finally, we implemented the osmostat in MC simulations initialized with bulk water that are driven toward electrolyte solutions of similar concentration as the saline reservoir. Our results demonstrate that alchemical osmostat for MC simulation is a promising tool for use to sample electrolyte insertion in nanoporous materials.

Quantifying the Structure of Water and Hydrated Monovalent Ions by Density Functional Theory-Based Molecular Dynamics

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.

The Contribution of the Ion-Ion and Ion-Solvent Interactions in a Molecular Thermodynamic Treatment of Electrolyte Solutions

Developing molecular equations of state to treat electrolyte solutions is challenging due to the long-range nature of the Coulombic interactions. Seminal approaches commonly used are the mean spherical approximation (MSA) and the Debye-Hückel (DH) theory to account for ion-ion interactions and, often, the Born theory of solvation for ion-solvent interactions. We investigate the accuracy of the MSA and DH approaches using each to calculate the contribution of the ion-ion interactions to the chemical potential of NaCl in water, comparing these with newly computer-generated simulation data; the ion-ion contribution is isolated by selecting an appropriate primitive model with a Lennard-Jones force field to describe the solvent. A study of mixtures with different concentrations and ionic strengths reveals that the calculations from both MSA and DH theories are of similar accuracy, with the MSA approach resulting in marginally better agreement with the simulation data. We also demonstrate that the Born theory provides a good qualitative description of the contribution of the ion-solvent interactions; we employ an explicitly polar water model in these simulations. Quantitative agreement up to moderate salt concentrations and across the relevant range of temperature is achieved by adjusting the Born radius using simulation data of the free energy of solvation. We compute the radial and orientational distribution functions of the systems, thereby providing further insight on the differences observed between the theory and simulation. We thus provide rigorous benchmarks for use of the MSA, DH, and Born theories as perturbation approaches, which will be of value for improving existing models of electrolyte solutions, especially in the context of equations of state.

Non-Ewald methods for evaluating the electrostatic interactions of charge systems: similarity and difference

In molecular simulations, it is essential to properly calculate the electrostatic interactions of particles in the physical system of interest. Here we consider a method called the non-Ewald method, which does not rely on the standard Ewald method with periodic boundary conditions, but instead relies on the cutoff-based techniques. We focus on the physicochemical and mathematical conceptual aspects of the method in order to gain a deeper understanding of the simulation methodology. In particular, we take into account the reaction field (RF) method, the isotropic periodic sum (IPS) method, and the zero-multipole summation method (ZMM). These cutoff-based methods are based on different physical ideas and are completely distinguishable in their underlying concepts. The RF and IPS methods are "additive" methods that incorporate information outside the cutoff region, via dielectric medium and isotropic boundary condition, respectively. In contrast, the ZMM is a "subtraction" method that tries to remove the artificial effects, generated near the boundary, from the cutoff sphere. Nonetheless, we find physical and/or mathematical similarities between these methods. In particular, the modified RF method can be derived by the principle of neutralization utilized in the ZMM, and we also found a direct relationship between IPS and ZMM.

The ability of trimethylamine N-oxide to resist pressure induced perturbations to water structure

Trimethylamine N-oxide (TMAO) protects organisms from the damaging effects of high pressure. At the molecular level both TMAO and pressure perturb water structure but it is not understood how they act in combination. Here, we use neutron scattering coupled with computational modelling to provide atomistic insight into the structure of water under pressure at 4 kbar in the presence and absence of TMAO. The data reveal that TMAO resists pressure-induced perturbation to water structure, particularly in retaining a clear second solvation shell, enhanced hydrogen bonding between water molecules and strong TMAO - water hydrogen bonds. We calculate an 'osmolyte protection' ratio at which pressure and TMAO-induced energy changes effectively cancel out. Remarkably this ratio translates across scales to the organism level, matching the observed concentration dependence of TMAO in the muscle tissue of organisms as a function of depth. Osmolyte protection may therefore offer a molecular mechanism for the macroscale survival of life in extreme environments.

Structure of Aqueous KNO3 Solutions by Wide-Angle X-ray Scattering and Density Functional Theory

The structure of aqueous KNO₃ solutions was studied by wide-angle X-ray scattering (WAXS) and density functional theory (DFT). The interference functions were subjected to empirical potential structure refinement (EPSR) modeling to extract the detailed hydration structure information. In aqueous KNO₃ solutions with a concentration from 0.50 to 2.72 mol·dm^-3, water molecules coordinate K⁺ with a mean coordination number (CN) from 6.6 ± 0.9 to 5.8 ± 1.2, respectively, and a K-O_W(H₂O) distance of 2.82 Å. To further describe the hydration properties of K⁺, a hydration factor (f_h) was defined based on the orientational angle between the water O-H vector and the ion-oxygen vector. The f_h value obtained for K⁺ is 0.792, 0.787, 0.766, and 0.765, and the corresponding average hydration numbers (HN) are 5.2, 5.1, 4.8, and 4.5. The reduced HN is compensated by the direct binding of oxygen atoms of NO₃^- with the average CN from 0.3 ± 0.7 to 2.6 ± 1.7, respectively, and the K-O_N(NO₃^-) distance of 2.82 Å. The average number of water molecules around NO₃^- slightly reduces from 10.5 ± 1.5 to 9.6 ± 1.7 with r_N-OW = 3.63 Å. K⁺-NO₃^- ion association is characterized by K-O_N and K-N pair correlation functions (PCFs). A K-N peak is observed at 3.81 Å as the main peak with a shoulder at 3.42 Å in g_K-N(r). This finding indicates that two occupancy sites are available for K⁺, i.e., the edge-shared bidentate (BCIP) and the corner-shared monodentate (MCIP) contact ion pairs. The structure and stability of the KNO₃(H₂O)₁₀ cluster were investigated by a DFT method and cross-checked with the results from WAXS.

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.

Structural Analysis of Molecular Materials Using the Pair Distribution Function

This is a review of atomic pair distribution function (PDF) analysis as applied to the study of molecular materials. The PDF method is a powerful approach to study short- and intermediate-range order in materials on the nanoscale. It may be obtained from total scattering measurements using X-rays, neutrons, or electrons, and it provides structural details when defects, disorder, or structural ambiguities obscure their elucidation directly in reciprocal space. While its uses in the study of inorganic crystals, glasses, and nanomaterials have been recently highlighted, significant progress has also been made in its application to molecular materials such as carbons, pharmaceuticals, polymers, liquids, coordination compounds, composites, and more. Here, an overview of applications toward a wide variety of molecular compounds (organic and inorganic) and systems with molecular components is presented. We then present pedagogical descriptions and tips for further implementation. Successful utilization of the method requires an interdisciplinary consolidation of material preparation, high quality scattering experimentation, data processing, model formulation, and attentive scrutiny of the results. It is hoped that this article will provide a useful reference to practitioners for PDF applications in a wide realm of molecular sciences, and help new practitioners to get started with this technique.

Overlapping hydration shells in salt solutions causing non-monotonic Soret coefficients with varying concentration

We develop an intuitive picture that overlapping hydration shells in salt solutions cause non-monotonic Soret coefficients with varying concentration.

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.

Computational and solubility equilibrium experimental insight into Ca2+-fluoride complexation and their dissociation behaviors in aqueous solutions: implication for the association constant measured using fluoride ion selective electrodes

Although the Ca²⁺-F^- association is of great importance for aqueous environments and industrial systems containing F^-, as well as for defluorination processes, many details of the association solvation structures and behavior remain unclear. Herein, a combination of classical/ab initio molecular dynamics simulations and density functional theory calculations was used to investigate the structure and hydration of CaF_x^2-x (x = 1, 2) and the association/dissociation behavior of Ca²⁺-F^- in aqueous CaF₂ solutions. The primary shell of Ca²⁺ is found to be very flexible in the association of Ca²⁺-F^-, with coordination numbers dynamically oscillating in the range of 6-9, with 6 and 7 being the most favorable. The calculations show that for CaF(H₂O)₁₄⁺, the contact ion pair (CIP) is more favorable and occurs with no energy barrier, whereas the formation of CaF₂(aq.) must overcome a ∼3.6 kJ mol^-1 energy barrier; moreover, the CIP and solvent shared ion pair (SSIP) dynamically coexist for CaF₂(H₂O)₁₄ in aqueous CaF₂ solutions. Calculations for the dissociation process of CaF(H₂O)₆⁺ show a dramatic energy increase going from SSIP to free Ca²⁺ and F^-, ascribed to the surprisingly long-range electrostatic attraction between Ca²⁺ and F^- rather than to special F⋯H interactions. The energy increase results in the estimated association constant of CaF⁺ being larger than that previously measured using fluoride ion selective electrodes. This is attributed to the fact that the latter value might correspond to the ligand reaction of free Ca²⁺ and F^- to form the Ca²⁺-F^- SSIP. The combination of these results with CaF_2(s) solubility measurements suggests that the higher-order Ca²⁺-F^- complexes are absent in aqueous CaF₂ solutions.

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.

Renormalized site density functional theory for models of ion hydration

The development of accurate statistical mechanics models of molecular liquid systems is a problem of great practical and fundamental importance. Site-density functional theory (SDFT) is one of the promising directions in this area, but its success hinges upon the ability to efficiently reconcile the co-existence of two distinct intra- and inter-molecular interaction regimes in a molecular liquid. The renormalized formulation of SDFT (RSDFT), which we have recently developed, resolves this problem by introducing an additional potential field variable that decouples two interaction scales and maps the molecular liquid problem onto the effective simple liquid mixture. This work provides a critical assessment of RSDFT for the hydrated ion system-a problem that historically has always been one of the most difficult cases for SDFT applications. Using a two-site model of water, we perform a comprehensive analysis of hydrated alkali metal and halogen ions, including both structural and free energy based characteristics. The results indicate that RSDFT provides a significant improvement over conventional three-dimensional reference interaction site model implementations and may prove useful in coarse grained simulations based on two-site solvent models.

The Development of the Innovative Synthesis Methodology of Albumin Nanoparticles Supported by Their Physicochemical, Cytotoxic and Hemolytic Evaluation

Many studies are being performed to develop effective carriers for controlled cytostatic delivery wherein albumin is a promising material due to its tendency to accumulate near cancer cells. The novelty of this work involves the development of the synthesis methodology of albumin nanoparticles and their biological and physicochemical evaluation. Albumin particles were obtained via the salt-induced precipitation and K3PO4 was used as a salting-out agent. Various concentrations of protein and salting-out agent solutions were mixed using a burette or a syringe system. It was proved that the size of the particles depended on the concentrations of the reagents and the methodology applied. As a result of a process performed using a burette and 2 M K3PO4, albumin spheres having a size 5-25 nm were obtained. The size of nanospheres and their spherical shape was confirmed via TEM analysis. The use of a syringe system led to preparation of particles of large polydispersity. The highest albumin concentration allowing for synthesis of homogeneous particles was 2 g/L. The presence of albumin in spheres was confirmed via the FT-IR technique and UV-Vis spectroscopy. All samples showed no cytotoxicity towards normal human dermal fibroblasts and no hemolytic properties against human erythrocytes (the hemolysis did not exceed 2.5%).

A unified computational view of DNA duplex, triplex, quadruplex and their donor-acceptor interactions

DNA can assume various structures as a result of interactions at atomic and molecular levels (e.g., hydrogen bonds, π–π stacking interactions, and electrostatic potentials), so understanding of the consequences of these interactions could guide development of ways to produce elaborate programmable DNA for applications in bio- and nanotechnology. We conducted advanced ab initio calculations to investigate nucleobase model structures by componentizing their donor-acceptor interactions. By unifying computational conditions, we compared the independent interactions of DNA duplexes, triplexes, and quadruplexes, which led us to evaluate a stability trend among Watson–Crick and Hoogsteen base pairing, stacking, and even ion binding. For a realistic solution-like environment, the influence of water molecules was carefully considered, and the potassium-ion preference of G-quadruplex was first analyzed at an ab initio level by considering both base-base and ion-water interactions. We devised new structure factors including hydrogen bond length, glycosidic vector angle, and twist angle, which were highly effective for comparison between computationally-predicted and experimentally-determined structures; we clarified the function of phosphate backbone during nucleobase ordering. The simulated tendency of net interaction energies agreed well with that of real world, and this agreement validates the potential of ab initio study to guide programming of complicated DNA constructs.

Core and Valence Level Photoelectron Spectroscopy of Nanosolvated KCl

The solvation of alkali and halide ions in the aqueous environment has been a subject of intense experimental and theoretical research with multidisciplinary interests; yet, a comprehensive molecular-level understanding has still not been obtained. In recent years, electron spectroscopy has been increasingly applied to study the electronic and structural properties of aqueous ions with implications, especially in atmospheric chemistry. In this work, we report core and valence level (Cl 2p, Cl 3p, and K 3p) photoelectron spectra of the common alkali halide, KCl, doped in gas-phase water clusters in the size range of a few hundred water molecules. The results indicate that the electronic structure of these nanosolutions shows a distinct character from that observed at the liquid-vapor interface in liquid microjets and ambient pressure setups. Insights are provided into the unique solvation properties of ions in a nanoaqueous environment, emerging properties of bulk electrolyte solutions with growing cluster size, and sensitivity of the electronic structure to varying solvation configurations.

Role of water model on ion dissociation at ambient conditions

We study ion pair dissociation in water at ambient conditions using a combination of classical and ab initio approaches. The goal of this study is to disentangle the sources of discrepancy observed in computed potentials of mean force. In particular, we aim to understand why some models favor the stability of solvent-separated ion pairs vs contact ion pairs. We found that some observed differences can be explained by non-converged simulation parameters. However, we also unveil that for some models, small changes in the solution density can have significant effects on modifying the equilibrium balance between the two configurations. We conclude that the thermodynamic stability of contact and solvent-separated ion pairs is very sensitive to the dielectric properties of the underlying simulation model. In general, classical models are very robust in providing a similar estimation of the contact ion pair stability, while this is much more variable in density functional theory-based models. The barrier to transition from the solvent-separated to contact ion pair is fundamentally dependent on the balance between electrostatic potential energy and entropy. This reflects the importance of water intra- and inter-molecular polarizability in obtaining an accurate description of the screened ion-ion interactions.

Binary structure and dynamics of the hydrogen bonds in the hydration shells of ions

Ion-specific effects of cations (Li+, Na+, K+, Mg2+, Ca2+) and anions (F-, Cl-) on the hydrogen bond structure and dynamics of the coordination waters in the hydration shells have been studied using molecular dynamics simulations. Our simulations indicate that the hydrogen bonds between the first and second hydration shell waters show binary structural and dynamic properties. The hydrogen bond with a first shell water as the donor (HD) is strengthened, while those with a first shell water as the acceptor (HA) are weakened. For a hydrated anion, this binary effect reverses, but is less significant. This ion-specific binary effect correlates with the size and the valence of the ion, and is more significant for the strong kosmotropic ions of high charge density.

Resolving the Equal Number Density Puzzle: Molecular Picture from Simulations of LiCl(aq) and NaCl(aq)

The change in number densities of aqueous solutions of alkali chlorides should be qualitatively predictable. Typically, as cations get larger, the number density of the solution decreases. However, aqueous solutions of lithium and sodium chloride exhibit at ambient conditions practically identical number densities at equal molalities despite different ionic sizes. Here, we provide an atomistic interpretation of this experimentally observed anomalous behavior using molecular dynamics simulations. The obtained results show that the rigidity of the Li⁺ first and second solvation shells and the associated compromised hydrogen bonding result in practically equal average water densities in the local hydration regions for Li⁺ and Na⁺ despite different sizes of the cations. In addition, in more distant regions from the cations, the water densities of these two solutions also coincide. These findings thus provide an atomistic interpretation for matching number densities of LiCl and NaCl solutions. In contrast, the number density differences between NaCl and KCl solutions as well as between LiCl and KCl solutions behave in a regular fashion with lower number densities of solutions observed for larger cations.

Protein Nanofibrils and Their Hydrogel Formation with Metal Ions

Protein nanofibrils (PNFs) have been prepared by whey protein fibrillation at low pH and in the presence of different metal ions. The effect of the metal ions was systematically studied both in terms of PNF suspension gelation behavior and fibrillation kinetics. A high valence state and a small ionic radius (e.g., Sn4+) of the metal ion resulted in the formation of hydrogels already at a metal ion concentration of 30 mM, whereas an intermediate valence state and larger ionic radius (Co2+, Ni2+, Al3+) resulted in the hydrogel formation occurring at 60 mM. A concentration of 120 mM of Na+ was needed to form a PNF hydrogel, while lower concentrations showed liquid behaviors similar to the reference PNF solution where no metal ions had been introduced. The hydrogel mechanics were investigated at steady-state conditions after 24 h of incubation/gelation, revealing that more acidic (smaller and more charged) metal ions induced ca. 2 orders of magnitude higher storage modulus as compared to the less acidic metal ions (with smaller charge and larger radius) for the same concentration of metal ions. The viscoelastic nature of the hydrogels was attributed to the ability of the metal ions to coordinate water molecules in the vicinity of the PNFs. The presence of metal ions in the solutions during the growth of the PNFs typically resulted in curved fibrils, whereas an upper limit of the concentration existed when oxides/hydroxides were formed, and the hydrogels lost their gel properties due to phase separation. Thioflavin T (ThT) fluorescence was used to determine the rate of the fibrillation to form 50% of the total PNFs (t1/2), which decreased from 2.3 to ca. 0.5 h depending on the specific metal ions added.

Salt-Induced Liquid-Liquid Phase Separation: Combined Experimental and Theoretical Investigation of Water-Acetonitrile-Salt Mixtures

Salt-induced liquid-liquid phase separation in liquid mixtures is a common phenomenon in nature and in various applications, such as in separation and extraction of chemicals. Here, we present results of a systematic investigation of the phase behaviors in water-acetonitrile-salt mixtures using a combination of experiment and theory. We obtain complete ternary phase diagrams for nine representative salts in water-acetonitrile mixtures by cloud point and component analysis. We construct a thermodynamic free energy model by accounting for the nonideal mixing of the liquids, ion hydration, electrostatic interactions, and Born energy. Our theory yields phase diagrams in good agreement with the experimental data. By comparing the contributions due to the electrostatic interaction, Born energy, and hydration, we find that hydration is the main driving force for the liquid-liquid separation and is a major contributor to the specific ion effects. Our theory highlights the important role of entropy in the hydration driving force. We discuss the implications of our findings in the context of salting-out assisted liquid-liquid extraction and make suggestions for selecting salt ions to optimize the separation performance.

A Ternary Hybrid-Cation Room-Temperature Liquid Metal Battery and Interfacial Selection Mechanism Study

The dendrite-free sodium-potassium (Na-K) liquid alloy composed of two alkali metals is one of the ideal alternatives for Li metal as an anode material while maintaining large capacity, low potential, and high abundance. However, Na- or K-ion batteries have limited cathode materials that can deliver stably large capacity. Combining advantages of both, a hybrid-cation liquid metal battery is designed for a Li-ion-insertion-based cathode to deliver stable high capacity using a Na-K liquid anode to avoid dendrites. The mechanical property of the Na-K alloy is confirmed by simulation and experimental characterization, which leads to stable cycling performance. The charge carrier selection principle in this ternary hybrid-cation system is investigated, showing consistency with the proposed interfacial layer formation and ion distribution mechanism for the electrochemical process as well as the good stability. With Li ions contributing stable cycling as the cathode charge carrier, the K ion working as charge carrier on the anode, and Na as the medium to liquefy K metal, such a ternary hybrid battery system not only inherits the rich battery chemistry of Li-insertion cathodes but also broadens the understanding of alkali metal alloys and hybrid-ion battery chemistry.

Deciphering the Solvent Effect for the Solvation Structure of Ca2+ in Polar Molecular Liquids

Although the crystal structures for many inorganic compounds are readily available, researchers are still working hard to understand the relations between the structures and chemical properties of solutions because most of the chemical reactions take place in solutions. A huge amount of effort has been put toward modeling the ion solvation structure from the perspectives of both experiments and theories. In this study, the solvation structures of Ca²⁺ ions in aqueous and alcoholic solutions at different concentrations were carefully evaluated by Ca K-edge X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) analyses. Density functional theory (DFT) calculations were also performed to correlate the results with the experimental data and then further extended to other similar systems. It was found that the number of coordinating solvent molecules decreases with increasing Ca²⁺ concentration and increasing solvent molecule sizes. From the EXAFS data, it was observed that the first solvation shell of Ca²⁺ splits into two Ca-O distances in a methanol solution and the counter ion Cl^- might also be within the first shell at high concentrations. For the first time, the effects of solvents with different polarities and sizes on the ion solvation environment were systematically evaluated.

Vegetation canopy effects on total and dissolved Cl, Br, F and I concentrations in soil and their fate along the hydrological flow path

Although halogens are omnipresent in the environment, detailed understanding of processes involving chlorine (Cl), bromine (Br), fluorine (F) and iodine (I) in the terrestrial halogen cycle is still sparse. Our objectives were to (i) assess vertical depth profiles of total and water-extractable inorganic halogen concentrations (Cl_tot, Br_tot, F_tot, I_tot) in solid soil, (ii) test the effect of a tree canopy, and (iii) follow the fate of dissolved inorganic halogens along the hydrological flow path. More than 200 soil samples and ecosystem solutions (rainwater, soil solution, adit and creek water) collected in the Schwarzwald, SW Germany, were analyzed by combustion ion chromatography and ion chromatography for total and inorganic halogen concentrations. We found decreasing Cl_tot concentrations with increasing soil depth which were indicative of biological chlorination of organic matter and nutrient uplift, both associated with Cl accumulation in upper soil horizons. Vertical patterns of total Br, F and I were contrary to Cl_tot concentrations and were related significantly (positively) to pedogenic oxides, revealing their dependence on abiotic processes. The presence of a canopy at our study site resulted in significantly higher halogen concentrations in throughfall compared to rainfall and higher Br_tot concentrations in the organic layer. We attribute this difference to leaching from leaves and needles and wash-off of dry deposition. There were hardly any differences in halogen concentrations along the hydrological flow path except for significantly higher inorganic I concentrations in soil solution compared to rainfall due to equilibrium reactions between the soil solution and the solid soil phase. Highest inorganic F concentrations of up to 0.2 mg L^-1 were detected in creek water samples and may originate from the weathering of fluorite-bearing veins. Our study indicates halogen-specific processes underlying Cl, Br, I and F cycling in ecosystems.

A sensitive and specific nanosensor for monitoring extracellular potassium levels in the brain

Extracellular potassium concentration affects the membrane potential of neurons, and, thus, neuronal activity. Indeed, alterations of potassium levels can be related to neurological disorders, such as epilepsy and Alzheimer's disease, and, therefore, selectively detecting extracellular potassium would allow the monitoring of disease. However, currently available optical reporters are not capable of detecting small changes in potassium, in particular, in freely moving animals. Furthermore, they are susceptible to interference from sodium ions. Here, we report a highly sensitive and specific potassium nanosensor that can monitor potassium changes in the brain of freely moving mice undergoing epileptic seizures. An optical potassium indicator is embedded in mesoporous silica nanoparticles, which are shielded by an ultrathin layer of a potassium-permeable membrane, which prevents diffusion of other cations and allows the specific capturing of potassium ions. The shielded nanosensor enables the spatial mapping of potassium ion release in the hippocampus of freely moving mice.