Temperature and Concentration Effects on Li+-Ion Hydration. A Molecular Dynamics Simulation Study

Study on the Superlubricity Behavior of Ions under External Electric Fields at Steel Interfaces

Realizing macroscopic superlubricity in the presence of external electric fields (EEFs) at the steel interfaces is still challenging. In this work, macroscopic superlubricity with a coefficient of friction value of approximately 0.008 was realized under EEFs with the lubrication of LiPF6-based ionic liquids at steel interfaces. The roles of cations and anions in the superlubricity realization under EEFs were studied. Based on the experimental results, the macroscopic superlubricity behavior of Li(PEG)PF6 under EEFs at steel interfaces is attributed to the strong hydration effect of Li+ cations and the complete reactions of anions that contributed to the formation of a boundary film on the appropriate surface. Moreover, the reduction in the number of iron oxides in the boundary film on the disc was beneficial for friction reduction. We also provide a calculation model to describe the relationship between the hydration effect and the optimal voltage position, at which the lowest friction might occur. Ultimately, this work proves that macroscopic superlubricity can be realized under EEFs at steel interfaces and provides a foundation for engineering applications of superlubricity in an electrical environment.

Structure and dynamics of the Li+ ion in water, methanol and acetonitrile solvents: ab initio molecular dynamics simulations

Fundamental understanding of the structure and dynamics of the Li+ ion in solution is of utmost importance in different fields of science and technology, especially in the field of ion batteries. In view of this, ab initio molecular dynamics (AIMD) simulations of the LiCl salt in water, methanol and acetonitrile were performed to elucidate structural parameters such as radial distribution function and coordination number, and dynamical properties like diffusion coefficient, limiting ion conductivity and hydrogen bond correlation function. In the present AIMD simulation, one LiCl in water is equivalent to 0.8 M, which is close to the concentration of the lithium salt used in the Li-ion battery. The first sphere of coordination number of the Li+ ion was reaffirmed to be 4. The radial distribution function for different pairs of atoms is seen to be in good agreement with the experimental results. The calculated potential of mean force indicates the stronger interaction of the Li+ ion with methanol over water followed by acetonitrile. The dynamical parameters convey quite high diffusion and limiting ionic conductivity of the Li+ ion in acetonitrile compared to that in water and methanol which has been attributed to the transport of the Li-Cl ion pair in a non-dissociated form in acetonitrile. The AIMD results were found to be in accordance with the experimental findings, i.e. the limiting ion conductivity was found to follow the order acetonitrile > methanol > water. This study shows the importance of atomistic level simulations in evaluating the structural and dynamical parameters and in implementing the results for predicting and synthesizing better next generation solvents for lithium ion batteries (LIBs).

Phase Diagram of Aqueous Solutions of LiCl: a Study of Concentration Effects on the Anomalies of Water

We perform molecular dynamics simulations in order to study thermodynamics and the structure of supercooled aqueous solutions of lithium chloride (LiCl) at concentrations c = 0.678 and 2.034 mol/kg. We model the solvent using the TIP4P/2005 potential and the ions using the Madrid-2019 force field, a force field particularly suited for studying this solution. We find that, for c = 0.678 mol/kg, the behavior of the equation of state, studied in the P-T plane, indicates the presence of a liquid-liquid phase transition, similar to what was previously found for bulk water. We estimate the position of the liquid-liquid critical point to be at T_c ≈ 174 K, P_c ≈ 1775 bar, and ρ_c ≈ 1.065 g/cm³. When the concentration is tripled to c = 2.034 mol/kg, no critical point is observed, indicating its possible disappearance at this concentration. We also study the water-water and water-ions structure in the two solutions, and we find that at the concentrations examined the effect of ions on the water-water structure is not strong, and all the features found in bulk water are preserved. We also calculate the hydration number of the Li and Cl ions, and in line with experiments, we find the value of 4 for Li⁺ and between 5.5 and 6 for Cl^-, confirming the good performances of the Madrid-2019 force field.

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.

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.

Salting-out effect promoting highly efficient ambient ammonia synthesis

The electroreduction of nitrogen to ammonia offers a promising alternative to the energy-intensive Haber-Bosch process. Unfortunately, the reaction suffers from low activity and selectivity, owing to competing hydrogen evolution and the poor accessibility of nitrogen to the electrocatalyst. Here, we report that deliberately triggering a salting-out effect in a highly concentrated electrolyte can simultaneously tackle the above challenges and achieve highly efficient ammonia synthesis. The solute ions exhibit strong affinity for the surrounding H₂O molecules, forming a hydration shell and limiting their efficacy as both proton sources and solvents. This not only effectively suppresses hydrogen evolution but also ensures considerable nitrogen flux at the reaction interface via heterogeneous nucleation of the precipitate, thus facilitating the subsequent reduction process in terms of both selectivity and activity. As expected, even when assembled with a metal-free electrocatalyst, a high Faradaic efficiency of 71 ± 1.9% is achieved with this proof-of-concept system.

Influence of Inorganic Solution Components on Lithium Carbonate Crystal Growth

Lithium-bearing brines are an increasingly attractive source of Li for extraction. One extraction mechanism is the removal of Li from the fluid phase through the precipitation of zabuyelite (Li₂CO₃). The chemistry of the brine plays an important role in this process because ions in solution can compete for the components of the Li-carbonate phase. Here we explore the effect of different brine components on the precipitation of zabuyelite using experiments and computational simulations. Crystals formed in all solutions showed morphological evidence for potential transformation from a precursor phase. Our study indicates that Ca²⁺ and SO₄ ^2- are incorporated into the precipitated zabuyelite crystals. Sulfate also interacts directly with specific surfaces on the growing crystal and is expected to form ion pairs with Li⁺ in solution. Similarly, Na⁺ appears to form ion pairs in solution with the carbonate ion, slowing nucleation of zabuyelite in the experiments. K⁺ and Cl^- may interact with the growing zabuyelite crystals but do not appear to affect zabuyelite nucleation and growth times. These experiments highlight the importance of understanding the solution chemistry on zabuyelite formation in order to predict the efficiency of extraction processes and the purity of the solids.

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

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

Spontaneous Formation of KCl Aggregates in Biomolecular Simulations: A Force Field Issue?

Realistic all-atom simulation of biological systems requires accurate modeling of both the biomolecules and their ionic environment. Recently, ion nucleation phenomena leading to the rapid growth of KCl or NaCl clusters in the vicinity of biomolecular systems have been reported. To better understand this phenomenon, molecular dynamics simulations of KCl aqueous solutions at three (1.0, 0.25, and 0.10 M) concentrations were performed. Two popular water models (TIP3P and SPC/E) and two Lennard-Jones parameter sets (AMBER and Dang) were combined to produce a total of 80 ns of molecular dynamics trajectories. Results suggest that the use of the Dang cation Lennard-Jones parameters instead of those adopted by the AMBER force-field produces a more accurate description of the ionic solution. In the later case, formation of salt aggregates is probably indicative of an artifact resulting from misbalanced force-field parameters. Because similar results were obtained with two different water parameter sets, the simulations exclude a water model dependency in the formation of anomalous ionic clusters. Overall, the results strongly suggest that for accurate modeling of ions in biomolecular systems, great care should be taken in choosing balanced ionic parameters even when using the most popular force-fields. These results invite a reexamination of older data obtained using available force-fields and a thorough check of the quality of current parameters sets by performing simulations at finite (>0.25 M) instead of minimal salt conditions.

Dynamics of Water Associated with Lithium Ions Distributed in Polyethylene Oxide

The dynamics of water in polyethylene oxide (PEO)/LiCl solution has been studied with quasielastic neutron scattering experiments and molecular dynamics (MD) simulations. Two different time scales of water diffusion representing interfacial water and bulk water dynamics have been identified. The measured diffusion coefficient of interfacial water remained 5-10 times smaller than that of bulk water, but both were slowed by approximately 50% in the presence of Li(+). Detailed analysis of MD trajectories suggests that Li(+) is favorably found at the surface of the hydration layer, and the probability to find the caged Li(+) configuration formed by the PEO is lower than for the noncaged Li(+)-PEO configuration. In both configurations, however, the slowing down of water molecules is driven by reorienting water molecules and creating water-Li(+) hydration complexes. Performing the MD simulation with different ions (Na(+) and K(+)) revealed that smaller ionic radius of the ions is a key factor in disrupting the formation of PEO cages by allowing spaces for water molecules to come in between the ion and PEO.

Photoelectron spectroscopy and ab initio calculations of Li(H₂O)n(-) and Cs(H₂O)n(-) (n = 1-6) clusters

The Li(H2O)n(-) and Cs(H2O)n(-) (n = 0-6) clusters were studied using anion photoelectron spectroscopy combined with ab initio calculations. It was found that Li tends to be surrounded by water molecules with no water-water H-bonds being formed in the first hydration shell; while Cs sticks on the surface of water-water H-bonds network. The Li atom in its anionic or neutral state is surrounded by four water molecules through Li-O interactions within the first hydration shell; while the case of Cs is different. For the anionic Cs(H2O)n(-) clusters, two types of structures, namely H-end and O-end structures, were identified, with nearly degenerate energies. For the neutral Cs(H2O)n clusters, only O-end structures exist and the first hydration shell of the Cs atom has four water molecules. The different hydration nature of Li and Cs atoms can be ascribed to the delicate balance between the alkali metal-water interactions and the water-water interactions as well as the effect of excess electron.

Aqueous solutions: state of the art in ab initio molecular dynamics

The simulation of liquids by ab initio molecular dynamics (AIMD) has been a subject of intense activity over the last two decades. The significant increase in computational resources as well as the development of new and efficient algorithms has elevated this method to the status of a standard quantum mechanical tool that is used by both experimentalists and theoreticians. As AIMD computes the electronic structure from first principles, it is free of ad hoc parametrizations and has thus been applied to a large variety of physical and chemical problems. In particular, AIMD has provided microscopic insight into the structural and dynamical properties of aqueous solutions which are often challenging to probe experimentally. In this review, after a brief theoretical description of the Born-Oppenheimer and Car-Parrinello molecular dynamics formalisms, we show how AIMD has enhanced our understanding of the properties of liquid water and its constituent ions: the proton and the hydroxide ion. Thereafter, a broad overview of the application of AIMD to other aqueous systems, such as solvated organic molecules and inorganic ions, is presented. We also briefly describe the latest theoretical developments made in AIMD, such as methods for enhanced sampling and the inclusion of nuclear quantum effects.

Temperature effects on the structure and dynamics of the Jahn-Teller distorted Cr(2+) ion in aqueous solution: a hybrid QM/MM molecular dynamics simulation

The influences of an elevated temperature on the structure and dynamics of the Jahn-Teller distorted [Cr(H(2)O)(6)](2+) complex have been studied using an ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics simulation, showing that the increased temperature affects the lifetime distortions of the hydrated Cr(2+) ion by decreasing the inversion time to 0.5-2 ps.

The role of Li + , Na + , and K + in the ligand binding inside the human acetylcholinesterase gorge

Alkali cations can affect the catalytic efficiency of enzymes. This is particularly true when dealing with enzymes whose substrate bears a formal positive charge. Computational and biochemical approaches have been combined to shed light on the atomic aspects of the role of Li(+), Na(+), and K(+) on human acetylcholinesterase (hAChE) ligand binding. In this respect, molecular dynamics simulations and our recently developed metadynamics method were applied to study the entrance of the three cations in the gorge of hAChE, and their effect on the dynamical motion of a ligand (tetramethylammonium) from the bulk of the solvent into the deep narrow enzyme gorge. Furthermore, in order to support the theoretical results, K(M) and k(cat) for the acetylcholine hydrolysis in the presence of the three cations were evaluated by using an approach based on the Ellman's method. The combination of computational and biochemical experiments clearly showed that Li(+), Na(+), and K(+) may influence the ligand binding at the hAChE gorge.

First-Principles Molecular Dynamics Evaluation of Thermal Effects on the NMR1JLi,C Spin–Spin Coupling

Car-Parrinello (CP) molecular dynamics were applied to sample conformations of various models of organolithium aggregates which are chosen to estimate (1)J(Li,C) NMR coupling constants. The results show that the deviations from the values computed using static (optimized) geometries are small provided no large-amplitude motions occur within the timescale of the simulations. In the case of the vinyllithium dimer, for which rotation of the vinyl chain is observed, this approach allows analysis of the various contributions to the experimentally measured constants. For the trisolvated methyllithium monomer, partial decoordination of solvating dimethyl ether is observed and results in a significant shift of (1)J(Li,C). All these results highlight that a varied physicochemical machinery is hidden behind general empirical formulas, such as the Bauer-Winchester-Schleyer rule used experimentally.

X-ray and neutron scattering studies of the hydration structure of alkali ions in concentrated aqueous solutions

The presence of ions in water provides a rich and varied environment in which many natural processes occur with important consequences in biology, geology and chemistry. This article will focus on the structural properties of ions in water and it will be shown how the 'difference' methods of neutron diffraction with isotopic substitution (NDIS) and anomalous X-ray diffraction (AXD) can be used to obtain direct information regarding the radial pair distribution functions of many cations and anions in solution. This information can subsequently be used to calculate coordination numbers and to determine ion-water conformation in great detail. As well as enabling comparisons to be made amongst ions in particular groups in the periodic table, such information can also be contrasted with results provided by molecular dynamics (MD) simulation techniques. To illustrate the power of these 'difference' methods, reference will be made to the alkali group of ions, all of which have been successfully investigated by the above methods, with the exception of the radioactive element francium. Additional comments will be made on how NDIS measurements are currently being combined with MD simulations to determine the structure around complex ions and molecules, many of which are common in biological systems.

Neutron diffraction studies of electrolytes in null water: a direct determination of the first hydration zone of ions

A method of neutron diffraction is described which enables the first hydration zone of small cations to be investigated at atomic resolution. It is shown that the cation structures of aqueous electrolyte solutions dissolved in a 'null' mixture of water (H(2)O) and heavy water (D(2)O), can be calculated directly from the neutron scattering patterns. The hitherto unresolved structure around Na(+) is used to illustrate the power of this method, the accuracy of which is discussed formally with reference to standard nickel chloride solutions. Possible applications to a variety of other systems and at different thermodynamic states are proposed.

Temperature effects on the structural and dynamical properties of the Zn(II)-water complex in aqueous solution: a QM/MM molecular dynamics study

An ab initio quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) simulation at double-zeta restricted Hartree-Fock (RHF) level was performed at an elevated temperature of 363 K (90 degrees C) to study the temperature effects on the structural and dynamical properties of a Zn(II)-water complex in aqueous solution. The first hydration shell, consisting of 6 water molecules at a mean Zn-O distance of 2.16 A, was found to remain stable also at 90 degrees C with respect to exchange processes. The flexible second shell contains, in average, approximately 27 water ligands. To fully characterize the hydration structure, several other parameters such as radial and angular distribution functions (RDF and ADF) and tilt- and theta-angle distributions were evaluated and compared to data obtained at 298 K (25 degrees C). Temperature effects on the dynamics of the Zn(II)-water complex were studied in terms of water reorientations, mean ligand residence times (MRTs), and number of ligand exchange processes. To get further insight into the solute dynamics, additional data, in particular, librational and vibrational motions of water ligands and Zn-O stretching frequencies, were calculated. The second shell is considerably influenced by the elevated temperature, as the ligands' mean residence time is shortened to 4 ps from the value of 10.5 observed at room temperature. The values of the QM/MM MD simulation were also compared to the results of a classical molecular dynamics (CMD) simulation with two- plus three-body potential performed at 90 degrees C, revealing that an accurate description of the second shell and the dynamics of the Zn(II) hydrate needs the inclusion of quantum mechanics in the description.

Structure and Dynamics of Br- Ion in Liquid Methanol

A Car-Parrinello molecular dynamics simulation has been performed on a solution of Br- in liquid methanol analyzing with particular attention charge transfer and polarization effects. The first solvation shell has been characterized in terms of H-bonds, and it has been found that the high polarization of the bromide gives rise to a stable solvent cage. The differences in the coordination number with the chloride can be ascribed to the ionic radius and to the stronger perturbations brought by the solvent to the bromide ion.