Incorporating solvent effects in DFT: insights from cation exchange in faujasites

Zeolites are versatile materials renowned for their extra-framework cation exchange capabilities, with applications spanning diverse fields, including nuclear waste treatment. While detailed experimental characterization offers valuable insight, density functional theory (DFT) proves particularly adept at investigating ion exchange in zeolites, owing to its atomic and electronic resolution. However, the prevalent occurrence of zeolitic ion exchange in aqueous environments poses a challenge to conventional DFT modeling, traditionally conducted in a vacuum. This study seeks to enhance zeolite modeling by systematically evaluating predictive differences across varying degrees of aqueous solvent inclusion. Specifically focusing on monovalent cation exchange in Na-X zeolites, we explore diverse modeling approaches. These range from simple dehydrated systems (representing bare reference states in vacuum) to more sophisticated models that incorporate aqueous solvent effects through explicit water molecules and/or a dielectric medium. Through comparative analysis of DFT and semi-empirical DFT approaches, along with their validation against experimental results, our findings underscore the necessity to concurrently consider explicit and implicit solvent effects for accurate prediction of zeolitic ionic exchange.

An inorganic-rich but LiF-free interphase for fast charging and long cycle life lithium metal batteries

Li metal batteries using Li metal as negative electrode and LiNi1-x-yMnxCoyO2 as positive electrode represent the next generation high-energy batteries. A major challenge facing these batteries is finding electrolytes capable of forming good interphases. Conventionally, electrolyte is fluorinated to generate anion-derived LiF-rich interphases. However, their low ionic conductivities forbid fast-charging. Here, we use CsNO3 as a dual-functional additive to form stable interphases on both electrodes. Such strategy allows the use of 1,2-dimethoxyethane as the single solvent, promising superior ion transport and fast charging. LiNi1-x-yMnxCoyO2 is protected by the nitrate-derived species. On the Li metal side, large Cs+ has weak interactions with the solvent, leading to presence of anions in the solvation sheath and an anion-derived interphase. The interphase is surprisingly dominated by cesium bis(fluorosulfonyl)imide, a component not reported before. Its presence suggests that Cs+ is doing more than just electrostatic shielding as commonly believed. The interphase is free of LiF but still promises high performance as cells with high LiNi0.8Mn0.1Co0.1O2 loading (21 mg/cm2) and low N/P ratio (~2) can be cycled at 2C (~8 mA/cm2) with above 80% capacity retention after 200 cycles. These results suggest the role of LiF and Cs-containing additives need to be revisited.

An Electrochemical Perspective of Aqueous Zinc Metal Anode

Based on the attributes of nonflammability, environmental benignity, and cost-effectiveness of aqueous electrolytes, as well as the favorable compatibility of zinc metal with them, aqueous zinc ions batteries (AZIBs) become the leading energy storage candidate to meet the requirements of safety and low cost. Yet, aqueous electrolytes, acting as a double-edged sword, also play a negative role by directly or indirectly causing various parasitic reactions at the zinc anode side. These reactions include hydrogen evolution reaction, passivation, and dendrites, resulting in poor Coulombic efficiency and short lifespan of AZIBs. A comprehensive review of aqueous electrolytes chemistry, zinc chemistry, mechanism and chemistry of parasitic reactions, and their relationship is lacking. Moreover, the understanding of strategies for suppressing parasitic reactions from an electrochemical perspective is not profound enough. In this review, firstly, the chemistry of electrolytes, zinc anodes, and parasitic reactions and their relationship in AZIBs are deeply disclosed. Subsequently, the strategies for suppressing parasitic reactions from the perspective of enhancing the inherent thermodynamic stability of electrolytes and anodes, and lowering the dynamics of parasitic reactions at Zn/electrolyte interfaces are reviewed. Lastly, the perspectives on the future development direction of aqueous electrolytes, zinc anodes, and Zn/electrolyte interfaces are presented.

Cesium activates the neurotransmitter receptor for glycine

The monovalent cations sodium and potassium are crucial for the proper functioning of excitable cells, but, in addition, other monovalent alkali metal ions such as cesium and lithium can also affect neuronal physiology. For instance, there have been recent reports of adverse effects resulting from self-administered high concentrations of cesium in disease conditions, prompting the Food and Drug Administration (FDA) to issue an alert concerning cesium chloride. As we recently found that the monovalent cation NH4+ activates glycine receptors (GlyRs), we investigated the effects of alkali metal ions on the function of the GlyR, which belongs to one of the most widely distributed neurotransmitter receptors in the peripheral and central nervous systems. Whole-cell voltage clamp electrophysiology was performed with HEK293T cells transiently expressing different splice and RNA-edited variants of GlyR α2 and α3 homopentameric channels. By examining the influence of various milli- and sub-millimolar concentrations of lithium, sodium, potassium, and cesium on these GlyRs in comparison to its natural ligand glycine (0.1 mM), we could show that cesium activates GlyRs in a concentration- and post-transcriptional-dependent way. Additionally, we conducted atomistic molecular dynamic simulations on GlyR α3 embedded in a membrane bilayer with potassium and cesium, respectively. The simulations revealed slightly different GlyR-ion binding profiles for potassium and cesium, identifying interactions near the glycine binding pocket (potassium and cesium) and close to the RNA-edited site (cesium) in the extracellular GlyR domain. Together, these findings show that cesium acts as an agonist of GlyRs.

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.

Hydration Structure at the Calcite-Water (10.4) Interface in the Presence of Rubidium Chloride

Solid-liquid interfaces are of significant importance in a multitude of geochemical and technological fields. More specifically, the solvation structure plays a decisive role in the properties of the interfaces. Atomic force microscopy (AFM) has been used to resolve the interfacial hydration structure in the presence and absence of ions. Despite many studies investigating the calcite-water interface, the impact of ions on the hydration structure at this interface has rarely been studied. Here, we investigate the calcite-water interface at various concentrations (ranging from 0 to 5 M) of rubidium chloride (RbCl) using three-dimensional atomic force microscopy (3D AFM). We present molecularly resolved images of the hydration structure at the interface. Interestingly, the characteristic pattern of the hydration structure appears similar regardless of the RbCl concentration. The presence of the ions is detected in an indirect manner by more frequent contrast changes and slice displacements.

EXAFS spectroscopy measurements and ab initio molecular dynamics simulations reveal the hydration structure of the radium(II) ion

Radium is refocused from the viewpoint of an environmental pollutant and cancer therapy using alpha particles, where it mainly exists as a hydrated ion. We investigated the radium hydration structure and the dynamics of water molecules by extended X-ray absorption fine structure (EXAFS) spectroscopy and ab initio molecular dynamics (AIMD) simulation. The EXAFS experiment showed that the coordination number and average distance between radium ion and the oxygen atoms in the first hydration shell are 9.2 ± 1.9 and 2.87 ± 0.06 Å, respectively. They are consistent with those obtained from the AIMD simulations, 8.4 and 2.88 Å. The AIMD simulations also revealed that the water molecules in the first hydration shell of radium are less structured and more mobile than those of barium, which is an analogous element of radium. Our results indicate that radium can be more labile than barium in terms of interactions with water.

•

Extended X-ray absorption fine structure (EXAFS) measurement revealed the hydration structure of radium ion

•

Ab initio molecular dynamics (AIMD) simulation brought consistent results

•

AIMD revealed the structural and dynamic properties of the water molecules

•

The hydration structure of radium ion is more labile than that of barium ion

Kir7.1 disease mutant T153I within the inner pore affects K+ conduction

Inward-rectifier potassium channel 7.1 (Kir7.1) is present in polarized epithelium, including the RPE. A single amino acid change at position 153 in the KCNJ13 gene, a substitution of threonine to isoleucine in Kir7.1 protein, causes blindness. We hypothesized that the disease caused by this single amino acid substitution within the transmembrane protein domain could alter the translation, localization, or ion transport properties. We assessed the effects of amino acid side-chain length, arrangement, and polarity on channel structure and function. We showed that the T153I mutation yielded a full-length protein localized to the cell membrane. Whole-cell patch-clamp recordings and chord conductance analyses revealed that the T153I mutant channel had negligible K+ conductance and failed to hyperpolarize the membrane potential. However, the mutant channel exhibited enhanced inward current when Rb+ was used as a charge carrier, suggesting that an inner pore had formed, and the channel was dysfunctional. Substituting with a polar, non-polar, or short side-chain amino acid did not affect the localization of the protein. Still, it had an altered channel function due to differences in pore radius. Polar side chains (cysteine and serine) with inner pore radii comparable to wildtype exhibited normal inward K+ conductance. Short side-chains (glycine and alanine) produced a channel with wider than expected inner pore size and lacked the biophysical characteristics of the wildtype channel. Leucine substitution produced results similar to the T153I mutant channel. This study provides direct electrophysiological evidence for the structure and function of the Kir7.1 channel's narrow inner pore in regulating conductance.

Role of Counterions in the Adsorption and Micellization Behavior of 1:1 Ionic Surfactants at Fluid Interfaces─Demonstrated by the Standard Amphiphile System of Alkali Perfluoro-n-octanoates

In our latest communication, we proved experimentally that the ionic surfactant's surface excess is exclusively determined by the size of the hydrated counterion.[Lunkenheimer , Langmuir, 2017, 33, 10216-1022410.1021/acs.langmuir.7b00786]. However, at this stage of research, we were unable to decide whether this does only hold for the two or three lightest ions of lithium, sodium, and potassium, respectively. Alternatively, we could also consider the surface excess of the heavier hydrated alkali ions of potassium, rubidium, and cesium, having practically identical ion size, as being determined by the cross-sectional area of the related anionic extended chain residue. The latter assumption has represented state of art. Searching for reliable experimental results on the effect of the heavier counterions on the boundary layer, we have extended investigations to the amphiphiles' solutions of concentrations above the critical concentration of micelle formation (cmc).We provided evidence that the super-micellar solutions' equilibrium surface tension will remain constant provided the required conditions are followed. The related σ_cmc-value represents a parameter characteristic of the ionic surfactant's adsorption and micellization behavior. Evaluating the amphiphile's surface excess obtained from adsorption as a function of the related amphiphile's σ_cmc-value enables you to calculate the radius of the hydrated counterion valid in sub- and super-micellar solution likewise. The σ_cmc-value is directly proportional to the counterion's diameter concerned. Taking additionally into account the radii of naked ions known from crystal research, we succeeded in exactly discriminating the hydrated alkali ions' size from each other. There is a distinct sequence of hydration radii in absolute scale following the inequality, Li⁺ > Na⁺ > K⁺ > (NH₄)⁺ > Rb⁺ > Cs⁺. Therefore, we have to extend our model of counterion effectiveness put forward in our previous communication. It represents a general principle of the counterion effect.

Cs-RHO Goes from Worst to Best as Water Enhances Equilibrium CO2 Adsorption via Phase Change

The strong affinity of water to zeolite adsorbents has made adsorption of CO₂ from humid gas mixtures such as flue gas nearly impossible under equilibrated conditions. Here, in this manuscript, we describe a unique cooperative adsorption mechanism between H₂O and Cs⁺ cations on Cs-RHO zeolite, which actually facilitates the equilibrium adsorption of CO₂ under humid conditions. Our data demonstrate that, at a relative humidity of 5%, Cs-RHO adsorbs 3-fold higher amounts of CO₂ relative to dry conditions, at a temperature of 30 °C and CO₂ pressure of 1 bar. A comparative investigation of univalent cation-exchanged RHO zeolites with H⁺, Li⁺, Na⁺, K⁺, Rb⁺, and Cs⁺ shows an increase of equilibrium CO₂ adsorption under humid versus dry conditions to be unique to Cs-RHO. In situ powder X-ray diffraction indicates the appearance of a new phase with Im3̅m symmetry after H₂O saturation of Cs-RHO. A mixed-cation exchanged NaCs-RHO exhibits similar phase transitions after humid CO₂ adsorption; however, we found no evidence of cooperativity between Cs⁺ and Na⁺ cations in adsorption, in single-component H₂O and CO₂ adsorption. We hypothesize based on previous Rietveld refinements of CO₂ adsorption in Cs-RHO zeolite that the observed phase change is related to solvation of extra-framework Cs⁺ cations by H₂O. In the case of Cs-RHO, molecular modeling results suggest that hydration of these cations favors their migration from an original D8R position to S8R sites. We posit that this movement enables a trapdoor mechanism by which CO₂ can interact with Cs⁺ at S8R sites to access the α-cage.

Solubilities in aqueous nitrate solutions that appear to reverse the law of mass action

Non-ideal aqueous electrolyte solutions have been studied since the start of the application of thermodynamics to chemistry in the late 19th century. The present study examines some of the most extreme non-ideal behavior ever observed: solubilities of alkali and NH₄⁺ nitrate salts in water that appear to behave the opposite of how the Law of Mass Action would predict. A literature review discovered that the solubilities of NH₄NO₃ and many alkali nitrate salts increases when another nitrate-bearing electrolyte is added to solution. These occurrences were in concentrated solutions with insufficient water to provide all ions their preferred hydration number without sharing waters between ions. This water deficit results in the formation of contact ion-pairs as well as larger ion-clusters. These ion-clusters may be favored when there is more than one type of monovalent cation present.

Hydration of Heavy Alkaline-Earth Cations Studied by Molecular Dynamics Simulations and X-ray Absorption Spectroscopy

The physicochemical properties of the three heaviest alkaline-earth cations, Sr²⁺, Ba²⁺, and Ra²⁺ in water have been studied by means of classical molecular dynamics (MD) simulations. A specific set of cation-water intermolecular potentials based on ab initio potential energy surfaces has been built on the basis of the hydrated ion concept. The polarizable and flexible model of water MCDHO2 was adopted. The theoretical-experimental comparison of structural, dynamical, energetic, and spectroscopical properties of Sr²⁺ and Ba²⁺ aqueous solutions is satisfactory, which supports the methodology developed. This good behavior allows a reasonable reliability for the predicted Ra²⁺ physicochemical data not experimentally determined yet. Simulated extended X-ray absorption fine-structure (EXAFS) and X-ray absorption near-edge spectroscopy spectra have been computed from the snapshots of the MD simulations and compared with the experimental information available for Sr²⁺ and Ba²⁺. For the Ra²⁺ case, the Ra L₃-edge EXAFS spectrum is proposed. Structural and dynamical properties of the aqua ions for the three cations have been obtained and analyzed. Along the [M(H₂O)_n]^m+ series, the M-O distance for the first-hydration shell is 2.57, 2.81, and 2.93 Å for Sr²⁺, Ba²⁺, and Ra²⁺, respectively. The hydration number also increases when one is going down along the group: 8.1, 9.4, and 9.8 for Sr²⁺, Ba²⁺, and Ra²⁺, respectively. Whereas [Sr(H₂O)₈]²⁺ is a typical aqua ion with a well-defined structure, the Ba²⁺ and Ra²⁺ hydration provides a picture exhibiting an average between the ennea- and the deca-hydration. These results show a similar chemical behavior of Ba²⁺ and Ra²⁺ aqueous solutions and support experimental studies on the removal of Ra-226 of aquifers by different techniques, where Ra²⁺ is replaced by Ba²⁺. A comparison of the heavy alkaline ions, Rb⁺ and Cs⁺, with the heavy alkaline-earth ions is made.

Nonselective cation permeation in an AMPA-type glutamate receptor

Fast excitatory synaptic transmission in the central nervous system relies on the AMPA-type glutamate receptor (AMPAR). This receptor incorporates a nonselective cation channel, which is opened by the binding of glutamate. Although the open pore structure has recently became available from cryo-electron microscopy (Cryo-EM), the molecular mechanisms governing cation permeability in AMPA receptors are not understood. Here, we combined microsecond molecular dynamic (MD) simulations on a putative open-state structure of GluA2 with electrophysiology on cloned channels to elucidate ion permeation mechanisms. Na⁺, K⁺, and Cs⁺ permeated at physiological rates, consistent with a structure that represents a true open state. A single major ion binding site for Na⁺ and K⁺ in the pore represents the simplest selectivity filter (SF) structure for any tetrameric cation channel of known structure. The minimal SF comprised only Q586 and Q587, and other residues on the cytoplasmic side formed a water-filled cavity with a cone shape that lacked major interactions with ions. We observed that Cl^- readily enters the upper pore, explaining anion permeation in the RNA-edited (Q586R) form of GluA2. A permissive architecture of the SF accommodated different alkali metals in distinct solvation states to allow rapid, nonselective cation permeation and copermeation by water. Simulations suggested Cs⁺ uses two equally populated ion binding sites in the filter, and we confirmed with electrophysiology of GluA2 that Cs⁺ is slightly more permeant than Na⁺, consistent with serial binding sites preferentially driving selectivity.

Confinement effects on the solvation structure of solvated alkaline metal cations in a single-digit 1T-MoS2 nanochannel: A first-principles study

Confinement plays an important role in determining ion transport in porous materials, which, in turn, may influence the performance of many energy storage and desalination devices. In this work, we combined density functional theory (DFT) with an implicit solvation model and ab initio molecular dynamics (AIMD) to investigate the effects of nanoconfinement on several solvated alkaline metal cations in a single-digit 1T-MoS₂ nanochannel. Our DFT calculations with a solvation model indicated that cations with stronger hydration energy introduce a higher number of co-intercalated water molecules into the channel, consistent with early experimental observation obtained for MXene (2D transition metal carbide) channels. The predicted optimal water numbers for the cations were then used for AIMD simulations that explicitly include the effects of the solvent. When compared with the cations in bulk solution, our simulations showed that the hydration structure and coordination number (CN) of the solvated cations confined in the MoS₂ channel can be significantly altered. We found that larger cations with weaker hydration energy (K⁺, Rb⁺, and Cs⁺) exhibited a distinctive CN decrease under confinement, while smaller cations (Li⁺ and Na⁺) retained a similar hydration shell as in the bulk solution. More specifically, the hydration shell of large cations (K⁺, Rb⁺, and Cs⁺) in MoS₂ showed similar features of the coordination angle to the bulk, which suggests the partially broken hydration shell with no geometry change under confinement. Our simulations provided insights into the change of the hydration structure of alkaline metal cations under confinement, which may have important implications on their transport in the 1T-MoS₂ channel.

Perspectives of fast magic-angle spinning 87 Rb NMR of organic solids at high magnetic fields

We report solid-state ⁸⁷ Rb NMR spectra from two Rb-ionophore complexes obtained with fast magic-angle spinning (MAS) (up to 60 kHz) at 21.1 T. These Rb-ionophore complexes containing macrocycles such as benzo-15-crown-5 and cryptand [2.2.2] are typical of organic Rb salts that exhibit very large ⁸⁷ Rb quadrupole coupling constants (close to 20 MHz). We have also obtained static ⁸⁷ Rb NMR spectra for these two compounds and determined both ⁸⁷ Rb quadrupole coupling and chemical shift tensors. The experimental ⁸⁷ Rb NMR tensor parameters are compared with those obtained by quantum chemical computations. Our results demonstrate that the combination of fast MAS (60 kHz or higher) and a high magnetic field (21.1 T or higher) is sufficient to produce high-quality solid-state ⁸⁷ Rb NMR spectra for organic Rb solids at the natural abundance level. We anticipate that, with additional ⁸⁷ Rb isotope enrichment (up to 99%), the sensitivity of solid-state ⁸⁷ Rb NMR will be 400 times higher than ³⁹ K NMR, which makes the former an attractive surrogate probe for studying K⁺ ion binding in biological systems.

An XAFS study of Cs adsorption by the precipitation bands of Mn-Fe-based Prussian blue analogues spontaneously formed in agarose gel

The adsorption of Cs⁺ ions by the precipitation bands of a Mn-Fe based Prussian blue analogue (Mn-Fe PBA) that form spontaneously in agarose gel was investigated by X-ray absorption fine structure spectroscopy coupled with scanning electron microscopy (SEM) and X-ray fluorescence (XRF) distribution analysis. Two gel samples were prepared by contacting a gel containing 0.05 M [Fe(CN)₆]^3- and 2.3 mass% agarose with a 0.50 M MnSO₄ solution, into one of which a 0.10 M CsCl solution was introduced. The SEM images and the XRF intensity distributions reveal that Mn-Fe PBA forms cubic crystallites (approx. 3 × 3 × 3 μm in size) in the gels that trap Cs⁺ ions with considerably high affinity. Cs L₃-edge and Mn K-edge X-ray absorption near-edge structure (XANES) spectra, which were analyzed with the aid of FEFF simulations, strongly suggest that Cs adsorption occurs at relatively large defect sites close to the sub-cube faces in the PBA. This suggestion is supported by Cs L₃-edge extended X-ray absorption fine structure spectroscopy, which suggested that the first and second coordination shells around the Cs⁺ ions are at a Cs-O distance of 0.35 ± 0.02 nm and a Cs-N distance of 0.43 ± 0.01 nm, respectively, with coordination numbers of 1.5 ± 0.5 and 3.0 ± 0.5. The Mn K-edge XANES data also suggest that H₂O molecules, which initially occupy many cubic centers in the Mn-Fe PBAs, are mostly displaced during Cs adsorption. These findings provide valuable insight toward fully understanding Cs adsorption by Mn-Fe PBA.

Ion-Specific and pH-Dependent Hydration of Mica-Electrolyte Interfaces

Hydration forces play a crucial role in a wide range of phenomena in physics, chemistry, and biology. Here, we study the hydration of mica surfaces in contact with various alkali chloride solutions over a wide range of concentrations and pH values. Using atomic force microscopy and molecular dynamics simulations, we demonstrate that hydration forces consist of a superposition of a monotonically decaying and an oscillatory part, each with a unique dependence on the specific type of cation. The monotonic hydration force gradually decreases in strength with decreasing bulk hydration energy, leading to a transition from an overall repulsive (Li⁺, Na⁺) to an attractive (Rb⁺, Cs⁺) force. The oscillatory part, in contrast, displays a binary character, being hardly affected by the presence of strongly hydrated cations (Li⁺, Na⁺), but it becomes completely suppressed in the presence of weakly hydrated cations (Rb⁺, Cs⁺), in agreement with a less pronounced water structure in simulations. For both aspects, K⁺ plays an intermediate role, and decreasing pH follows the trend of increasing Rb⁺ and Cs⁺ concentrations.

Many-Body Effects Determine the Local Hydration Structure of Cs+ in Solution

A systematic analysis of the hydration structure of Cs⁺ ions in solution is derived from simulations carried out using a series of molecular models built upon a hierarchy of approximate representations of many-body effects in ion-water interactions. It is found that a pairwise-additive model, commonly used in biomolecular simulations, provides poor agreement with experimental X-ray spectra, indicating an incorrect description of the underlying hydration structure. Although the agreement with experiment improves in simulations with a polarizable model, the predicted hydration structure is found to lack the correct sequence of water shells. Progressive inclusion of explicit many-body effects in the representation of Cs⁺-water interactions as well as accounting for nuclear quantum effects is shown to be necessary for quantitatively reproducing the experimental X-ray spectra. Besides emphasizing the importance of many-body effects, these results suggest that molecular models rigorously derived from many-body expansions hold promise for realistic simulations of aqueous solutions.