Molecular dynamics simulations of aqueous NaCl and KCl solutions: Effects of ion concentration on the single-particle, pair, and collective dynamical properties of ions and water molecules

Coalescence of multiple droplets induced by a constant DC electric field

In this work, the electro-coalescence process of three nanodroplets under a constant DC electric field is investigated via molecular dynamics simulations (MD), aiming to explore the electric manipulation of multiple droplets coalescence on the molecular level. The symmetrical and asymmetrical dynamic evolutions of electrocoalescence process can be observed. Our MD simulations show that there are two types of critical electric fields to induce the special dynamics. The chain configuration can be formed, when one of the critical electric field is exceeded, referred to as Ecc. On the other hand, there is another critical electric field to change the coalescence pattern from complete coalescence to partial coalescence, the so-called Ecn. Finally, we find that the use of the pulsed DC electric field can overcome the drawbacks of the constant DC electric field in the crude oil industry, and the mechanisms behind the suppressed effect of the water chain or non-coalescence are further revealed.

Anomalous Concentration Dependence of Viscosity: Hidden Role of Cross-Correlations in Aqueous Electrolyte Solutions

The viscosity of aqueous electrolyte solutions exhibits well-known composition-dependent anomalies that show certain definitive trends and universal features. The viscosity of LiCl and NaCl solutions increases with concentration in a monotonic fashion, while solutions of KCl, RbCl, and CsCl exhibit a more complex behavior. Here, the viscosity first decreases and then increases with increasing concentration, with a rather broad minimum at intermediate concentrations (ca. 1-3 m). To unearth the origin of such puzzling behavior, we carried out detailed molecular-level analyses by interrogating the exact Green-Kubo expression of viscosity in terms of the stress-stress time correlation function (SS-TCF). The total SS-TCF can be decomposed into a collection of three self- and three cross-SS-TCFs arising from the three constituent components (water, cations, and anions). Mode coupling theory (MCT) analysis for the friction on ions and the viscosity of the solution suggests the possible importance of two-particle static and time-dependent cross-correlations between water and the ions. We calculate the viscosity and other dynamical properties for all five electrolyte (LiCl, NaCl, KCl, RbCl, and CsCl) solutions over a range of concentrations, using two models of water (SPC/E and TIP4P/2005). The total viscosity derives non-negligible contributions from all of the terms. The cross-correlations are found to be surprisingly large and seen to play a hidden role in the concentration dependence. However, the importance of cross-correlations is often not discussed. Our study leads to a theoretical understanding of the microscopic origin of the observed anomalies in the composition dependence of viscosity across all five electrolytes.

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.

Toward a molecular understanding of the conductivity of lithium-ion conducting polyanion polymer electrolytes by molecular dynamics simulation

With the improved lithium-ion transference number near unity, the low conductivity of single lithium-ion conducting solid polymer electrolytes (SLIC-SPEs) still hinders their application in high-rate batteries. Though some empirical conclusions on the conducting mechanism of SLIC-SPEs have been obtained, a more comprehensive study on the quantitative relationship between the molecular structure factors and ionic conduction performance is expected. In this study, a model structure that contains adjustable main chain and anion groups in the polyethylene oxide (PEO) matrix was used to clarify the influence of molecular structural factors on ionic conductivity and electrochemical stability of SLIC-SPEs. The anionic group was further disassembled into the intermediate group and end group while the main chain structure was distinguished into different degrees of polymerization and various lengths of the spacers between anions. Therefore, a well-defined molecular structure was employed to describe its relationship with ionic conductivity. In addition, the dissociation degree of salts and mobility of ions changing with the molecular structure were also discussed to explore the fundamental causes of conductivity. It can be concluded that the anion group affects the conductivity mainly via the dissociation degree, while the main chain structure impacts the conductivity by both dissociation degree and mobility.

Molecular dynamics simulations of the dielectric constants of salt-free and salt-doped polar solvents

We develop a Stockmayer fluid model that accounts for the dielectric responses of polar solvents (water, MeOH, EtOH, acetone, 1-propanol, DMSO, and DMF) and NaCl solutions. These solvent molecules are represented by Lennard-Jones (LJ) spheres with permanent dipole moments and the ions by charged LJ spheres. The simulated dielectric constants of these liquids are comparable to experimental values, including the substantial decrease in the dielectric constant of water upon the addition of NaCl. Moreover, the simulations predict an increase in the dielectric constant when considering the influence of ion translations in addition to the orientation of permanent dipoles.

Dielectric relaxation and dielectric decrement in ionic acetamide deep eutectic solvents: Spectral decomposition and comparison with experiments

Frequency-dependent dielectric relaxation in three deep eutectic solvents (DESs), (acetamide+LiClO4/NO3/Br), was investigated in the temperature range, 329 ≤ T/K ≤ 358, via molecular dynamics simulations. Subsequently, decomposition of the real and the imaginary components of the simulated dielectric spectra was carried out to separate the rotational (dipole-dipole), translational (ion-ion), and ro-translational (dipole-ion) contributions. The dipolar contribution, as expected, was found to dominate all the frequency-dependent dielectric spectra over the entire frequency regime, while the other two components together made tiny contributions only. The translational (ion-ion) and the cross ro-translational contributions appeared in the THz regime in contrast to the viscosity-dependent dipolar relaxations that dominated the MHz-GHz frequency window. Our simulations predicted, in agreement with experiments, anion-dependent decrement of the static dielectric constant (ɛs ∼ 20 to 30) for acetamide (ɛs ∼ 66) in these ionic DESs. Simulated dipole-correlations (Kirkwood g factor) indicated significant orientational frustrations. The frustrated orientational structure was found to be associated with the anion-dependent damage of the acetamide H-bond network. Single dipole reorientation time distributions suggested slowed down acetamide rotations but did not indicate presence of any "rotationally frozen" molecule. The dielectric decrement is, therefore, largely static in origin. This provides a new insight into the ion dependence of the dielectric behavior of these ionic DESs. A good agreement between the simulated and the experimental timescales was also noticed.

A new methodology of understanding the mechanism of high shear wet granulation based on experiment and molecular dynamics stimulation

In high shear wet granulation (HSWG), the interaction mechanism between binder and powder with different sugar content is still unclear. Herein, the law and mechanism of the interaction between binder and powder were studied on the molecular level by combining experiment analysis through the Kriging model and molecular dynamics (MD) simulation. For the sticky powder with high sugar content, the ethanol in the binder played a pivotal role in dispersing water into powders, and the amount of water determined the growth of granules. In the saturating stage, the reduction of sugar content facilitates the penetration of ethanol molecules. The concentration of ethanol determines whether the mixture is blended uniformly in the merging stage. The simulation results are consistent with the actual situation and explain the competition mechanism of interaction with binder and powder. Therefore, this research offers an efficient strategy for the in-depth understanding of the HSWG process where the powder is sticky, as well as providing guidelines for the practical application of preparation for TCM granules.

Direct Correlation between Short-Range Vibrational Spectral Diffusion and Localized Ion-Cage Dynamics of Water-in-Salt Electrolytes

The molecular dynamics simulations of a "water-in-salt" electrolyte, lithium bis(trifluoromethyl sulfonyl) imide (LiNTf₂), with a varying concentration range of 3 to 20 m were performed to establish a direct connection between a dynamic property like the ion-cage lifetime with the short-range vibrational stretching frequency shift of the used probe, HOD. The properties reported here are compared to that obtained from experiments performed at the same concentrations. The time-series wavelet transform was adopted as a preferable mathematical tool for calculating the instantaneous fluctuating frequencies of the probe O-D stretch mode and the concentration-dependent vibrational stretch spectral signature based on the variable functions associated with a particular chemical bond derived from classical molecular dynamics trajectories. The decay time constants of frequency fluctuations and the lifetime of the ion cage (τ_IC) were estimated as a function of salt concentration. Herein, we emphasize the correlation between the slowest time constant (τ₃) of the decay of O-D stretch frequency fluctuations and the timescales associated with the lifetime of ion cages (τ_IC). The results exhibit that the existing relationships were also concentration-dependent. Therefore, this study highlights the connection between the ionic motions that regulate the overall system dynamics with the short-range vibrational frequency shift of the used probe, which was used similar to experiments. It also provides an understanding of the interionic interactions and the dynamical and spectral properties of the electrolytic mixtures. We establish a direct correlation between short-range frequency profile and localized ion-cage lifetime, which can fill the gap of understanding between viscosity, vibrational frequency, and ion-cage dynamics of electrolytes.

Molecular Simulation-Guided Spectroscopy of Imidazolium-Based Ionic Liquids and Effects of Methylation on Ion-Cage and -Pair Dynamics

Classical molecular dynamics simulations were performed to assess an atomistic interpretation of the ion-probe structural interactions in two typical ionic liquids (ILs), 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide [BMIm][NTf₂] and 1-butyl-2,3-dimethylimidazolium bis(trifluoromethylsulfonyl)imide [BDimIm][NTf₂] through computational ultrafast spectroscopy. The nitrile stretching vibrations of the thiocyanate anion, [SCN]^-, serve as the local mode of the ultrafast system dynamics within the imidazolium-based ionic liquid environment. The wavelet transform of classical trajectories determines the time-varying fluctuating frequencies and the stretch spectral signatures of SCN^- in the normalized distribution. However, computational modeling of the two-dimensional (2D) spectra from the wavelet-derived vibrational frequencies yields time evolution of the local molecular structure along with the varied time-dependent dynamics of the spectral diffusion process. We calculated the frequency-frequency correlation functions (FFCFs), time correlations associated with the ion-pair and -cage dynamics, and mean square displacements as a function of time, depicting diffusive dynamics. The calculated results based on the pair correlation functions and the distribution of atomic density suggest that the hydrogen and methylated carbon at the two-position of the imidazolium ring of [BMIm] and [BDimIm] cations, respectively, strongly interact with the probe through the N of the thiocyanate anion rather than the S atom. The center-of-mass center-of-mass (COM-COM) cation-probe radial distribution functions (RDFs) in conjunction with the site-specific structural analysis further reveal well-structured interactions of the thiocyanate ion and [BMIm]⁺ cation rather than the [BDimIm] cation. In contrast, the anion-probe COM-COM RDFs depict weak interactive associations within the vibrational probe [SCN]^- and [NTf₂]^- ions. Methylation at the two-position of the imidazolium ring predicts slower structural reorganization and breaking and reformation dynamics of the ion pairs and cages within the ionic liquid framework.

Enhanced Water Evaporation from Å-Scale Graphene Nanopores

Enhancing the kinetics of liquid-vapor transition from nanoscale confinements is an attractive strategy for developing evaporation and separation applications. The ultimate limit of confinement for evaporation is an atom thick interface hosting angstrom-scale nanopores. Herein, using a combined experimental/computational approach, we report highly enhanced water evaporation rates when angstrom sized oxygen-functionalized graphene nanopores are placed at the liquid-vapor interface. The evaporation flux increases for the smaller nanopores with an enhancement up to 35-fold with respect to the bare liquid-vapor interface. Molecular dynamics simulations reveal that oxygen-functionalized nanopores render rapid rotational and translational dynamics to the water molecules due to a reduced and short-lived water-water hydrogen bonding. The potential of mean force (PMF) reveals that the free energy barrier for water evaporation decreases in the presence of nanopores at the atomically thin interface, which further explains the enhancement in evaporation flux. These findings can enable the development of energy-efficient technologies relying on water evaporation.

Directional passive transport of nanodroplets on general axisymmetric surfaces

Rapid removal of small-sized droplets passively using fixed structures is a key challenge for various applications including anti-icing, rapid cooling, and water harvesting. In this work, we investigate the directional motion of nanodroplets on axisymmetric surfaces with curvature gradient through molecular dynamics (MD) simulations. It is found that as the shape of the axisymmetric surface is changed from a dome to a trumpet, the droplet velocity is greatly enhanced, by a factor of ∼14. Such an increase is mainly caused by the increment in the driving force. The droplet velocity changes nonlinearly as the surface wettability is varied and assumes the maximum at the contact angle of ∼75°. We derive a formula for the driving force of nanodroplets on general axisymmetric surfaces by evaluating the pressure gradient inside the droplet induced by the curvature gradient. Molecular dynamics simulations are performed to directly measure the driving force and confirm that the theoretical formula works well. By illustrating the reduced initial velocity of droplets as a function of a dimensionless number, which represents the ratio of the driving force to the retentive force due to contact angle hysteresis, we show that the onset of droplet motion on axisymmetric surfaces occurs when the dimensionless number is above a critical value. The dimensionless number reveals the effects of surface geometry, surface wettability, and droplet size on the droplet motion.

Fluorescence of KCl Aqueous Solution: A Possible Spectroscopic Signature of Nucleation

Ion pairing in water solutions alters both the water hydrogen-bond network and ion solvation, modifying the dynamics and properties of electrolyte water solutions. Here, we report an anomalous intrinsic fluorescence of KCl aqueous solution at room temperature and show that its intensity increases with the salt concentration. From the ab initio density functional theory (DFT) and time-dependent DFT modeling, we propose that the fluorescence emission could originate from the stiffening of the hydrogen bond network in the hydration shell of solvated ion-pairs that suppresses the fast nonradiative decay and allows the slower radiative channel to become a possible decay pathway. Because computations suggest that the fluorophores are the local ion-water structures present in the prenucleation phase, this band could be the signature of the incoming salt precipitation.

In Silico Drug Design of Benzothiadiazine Derivatives Interacting with Phospholipid Cell Membranes

The use of drugs derived from benzothiadiazine, a bicyclic heterocyclic benzene derivative, has become a widespread treatment for diseases such as hypertension, low blood sugar or the human immunodeficiency virus, among others. In this work we have investigated the interactions of benzothiadiazine and four of its derivatives designed in silico with model zwitterionic cell membranes formed by dioleoylphosphatidylcholine, 1,2-dioleoyl-sn-glycero-3-phosphoserine and cholesterol at the liquid–crystal phase inside aqueous potassium chloride solution. We have elucidated the local structure of benzothiadiazine by means of microsecond molecular dynamics simulations of systems including a benzothiadiazine molecule or one of its derivatives. Such derivatives were obtained by the substitution of a single hydrogen site of benzothiadiazine by two different classes of chemical groups, one of them electron-donating groups (methyl and ethyl) and another one by electron-accepting groups (fluorine and trifluoromethyl). Our data have revealed that benzothiadiazine derivatives have a strong affinity to stay at the cell membrane interface although their solvation characteristics can vary significantly—they can be fully solvated by water in short periods of time or continuously attached to specific lipid sites during intervals of 10–70 ns. Furthermore, benzothiadiazines are able to bind lipids and cholesterol chains by means of single and double hydrogen-bonds of characteristic lengths between 1.6 and 2.1 Å.

Long-range DNA-water interactions

DNA functions only in aqueous environments and adopts different conformations depending on the hydration level. The dynamics of hydration water and hydrated DNA leads to rotating and oscillating dipoles that, in turn, give rise to a strong megahertz to terahertz absorption. Investigating the impact of hydration on DNA dynamics and the spectral features of water molecules influenced by DNA, however, is extremely challenging because of the strong absorption of water in the megahertz to terahertz frequency range. In response, we have employed a high-precision megahertz to terahertz dielectric spectrometer, assisted by molecular dynamics simulations, to investigate the dynamics of water molecules within the hydration shells of DNA as well as the collective vibrational motions of hydrated DNA, which are vital to DNA conformation and functionality. Our results reveal that the dynamics of water molecules in a DNA solution is heterogeneous, exhibiting a hierarchy of four distinct relaxation times ranging from ∼8 ps to 1 ns, and the hydration structure of a DNA chain can extend to as far as ∼18 Å from its surface. The low-frequency collective vibrational modes of hydrated DNA have been identified and found to be sensitive to environmental conditions including temperature and hydration level. The results reveal critical information on hydrated DNA dynamics and DNA-water interfaces, which impact the biochemical functions and reactivity of DNA.

Structure of benzothiadiazine at zwitterionic phospholipid cell membranes

The use of drugs derived from benzothiadiazine, which is a bicyclic heterocyclic benzene derivative, has become a widespread treatment for diseases such as hypertension (treated with diuretics such as bendroflumethiazide or chlorothiazide), low blood sugar (treated with non-diuretic diazoxide), or the human immunodeficiency virus, among others. In this work, we have investigated the interactions of benzothiadiazine with the basic components of cell membranes and solvents, such as phospholipids, cholesterol, ions, and water. The analysis of the mutual microscopic interactions is of central importance to elucidate the local structure of benzothiadiazine as well as the mechanisms responsible for the access of benzothiadiazine to the interior of the cell. We have performed molecular dynamics simulations of benzothiadiazine embedded in three different model zwitterionic bilayer membranes made by dimyristoylphosphatidylcholine, dioleoylphosphatidylcholine, 1,2-dioleoyl-sn-glycero-3-phosphoserine, and cholesterol inside aqueous sodium-chloride solution in order to systematically examine microscopic interactions of benzothiadiazine with the cell membrane at liquid-crystalline phase conditions. From data obtained through radial distribution functions, hydrogen-bonding lengths, and potentials of mean force based on reversible work calculations, we have observed that benzothiadiazine has a strong affinity to stay at the cell membrane interface although it can be fully solvated by water in short periods of time. Furthermore, benzothiadiazine is able to bind lipids and cholesterol chains by means of single and double hydrogen-bonds of different characteristic lengths.

Individual Ion Activity Coefficients in Aqueous Electrolytes from Explicit-Water Molecular Dynamics Simulations

We compute individual ion activity coefficients (IIACs) in aqueous NaCl, KCl, NaF, and KF solutions from explicit-water molecular dynamics simulations. Free energy changes are obtained from insertion of single ions-accompanied by uniform neutralizing backgrounds-into solution by gradually turning on first Lennard-Jones interactions, followed by Coulombic interactions using Ewald electrostatics. Simulations are performed at multiple system sizes, and all results are extrapolated to the thermodynamic limit, thus eliminating any possible artifacts from the neutralizing backgrounds. Because of controversies associated with measurements of IIACs from electrochemical cells with ion-selective electrodes, the reported experimental data are not widely accepted; thus there remains a knowledge gap with respect to the contributions of individual ions to solution nonidealities. Our results are in good qualitative agreement with these reported measurements, though significantly larger in magnitude. In particular, the relative positioning for the activity coefficients of anions and cations matches the experimental ordering for all four systems. This work establishes a robust thermodynamic framework, without a need to invoke extra hypotheses, that sheds light on the behavior of individual ions and their contributions to nonidealities of aqueous electrolyte solutions.

Visualization of microwave near-field distribution in sodium chloride and glucose aqueous solutions by a thermo-elastic optical indicator microscope

In this study, a new optical method is presented to determine the concentrations of NaCl and glucose aqueous solutions by using a thermo-elastic optical indicator microscope. By measuring the microwave near-field distribution intensity, concentration changes of NaCl and glucose aqueous solutions were detected in the 0-100 mg/ml range, when exposed to microwave irradiation at 12 GHz frequency. Microwave near-field distribution intensity decreased as the NaCl or glucose concentration increased due to the changes of the absorption properties of aqueous solution. This method provides a novel approach for monitoring NaCl and glucose in biological liquids by using a CCD sensor capable of visualizing NaCl and glucose concentrations without scanning.

Failure of molecular dynamics to provide appropriate structures for quantum mechanical description of the aqueous chloride ion charge-transfer-to-solvent ultraviolet spectrum

The lowest band in the charge-transfer-to-solvent ultraviolet absorption spectrum of aqueous chloride ion is studied by experiment and computation. Interestingly, the experiments indicate that at concentrations up to at least 0.25 M, where calculations indicate ion pairing to be significant, there is no notable effect of ionic strength on the spectrum. The experimental spectra are fitted to aid comparison with computations. Classical molecular dynamic simulations are carried out on dilute aqueous Cl-, Na+, and NaCl, producing radial distribution functions in reasonable agreement with experiment and, for NaCl, clear evidence of ion pairing. Clusters are extracted from the simulations for quantum mechanical excited state calculations. Accurate ab initio coupled-cluster benchmark calculations on a small number of representative clusters are carried out and used to identify and validate an efficient protocol based on time-dependent density functional theory. The latter is used to carry out quantum mechanical calculations on thousands of clusters. The resulting computed spectrum is in excellent agreement with experiment for the peak position, with little influence from ion pairing, but is in qualitative disagreement on the width, being only about half as wide. It is concluded that simulation by classical molecular dynamics fails to provide an adequate variety of structures to explain the experimental CTTS spectrum of aqueous Cl-.

Interfacial and bulk properties of concentrated solutions of ammonium nitrate

We conducted molecular dynamics (MD) simulations to calculate the density and surface tension of concentrated ammonium nitrate (AN) solutions up to the solubility limit of ammonium nitrate in water, by combining the SPC/E, SPCE/F and TIP4P/2005 water models with OPLS model for ammonium and nitrate ions. This is the first time that the properties of concentrated solutions of nitrates, especially AN, have been studied by molecular dynamics. We effectively account for the polarisation effects by the electronic continuum correction (ECC), practically realised via rescaling of the ionic charges. We found that, the full-charge force field MD simulations overestimate the experimental results, as the ions experience repulsion from the interface and prefer to remain in the subsurface layer and the bulk solution. In contrast, reducing the ionic charges results in the behaviour that fits well with the experimental data. The nitrate anions display a greater propensity for the interface than the ammonium cations. We accurately predict both the density and the rise in the surface tension of concentrated solutions of AN, recommending TIP4P/2005 for water and the scaled-charge OPLS model (OPLS/ECC) for the ions in the solutions. We observe that, the adsorption of anions to the interface accompanies their depletion in the subsurface layer, which is preferentially occupied by cations, resulting in the formation of the electric double layer. We demonstrate the ion deficiency for up to 3 Å below the surface and establish the requirement to include the polarisability effects in the OPLS model for AN. While these results confirmed the findings of the previous studies for dilute solutions, they are new in the solubility limit. Concentrated solutions exhibit a strong effect of the abundance of solute on the coordination numbers of ions and on the degree of ion pairing. Surprisingly, ion pairing decreases significantly at the interface compared with the bulk. The present study identifies OPLS/ECC, along with TIP4P/2005, to yield accurate predictions of physical properties of concentrated AN, with precision required for industrial applications, such as a formulation of emulsion and fuel-oil explosives that now predominate the civilian use of AN. An application of this model will allow one to predict the surface properties of supersaturated solutions of AN which fall outside the capability of the present laboratory experiments but are important industrially.

Hydration breaking and chemical ordering in a levitated NaCl solution droplet beyond the metastable zone width limit: evidence for the early stage of two-step nucleation

For over two decades, NaCl nucleation from a supersaturated aqueous solution has been predicted to occur via a two-step nucleation (TSN) mechanism, i.e., two sequential events, the formation of locally dense liquid regions followed by structural ordering. However, the formation of dense liquid regions in the very early stage of TSN has never been experimentally observed. By using a state-of-the-art technique, a combination of electrostatic levitation (ESL) and in situ synchrotron X-ray and Raman scatterings, we find experimental evidence that indicates the formation of dense liquid regions in NaCl bulk solution at an unprecedentedly high level of supersaturation (S = 2.31). As supersaturation increases, evolution of ion clusters leads to chemical ordering, but no topological ordering, which is a precursor for forming the dense disordered regions of ion clusters at the early stage of TSN. Moreover, as the ion clusters proceed to evolve under highly supersaturated conditions, we observe the breakage of the water hydration structure indicating the stability limit of the dense liquid regions, and thus leading to nucleation. The evolution of solute clusters and breakage of hydration in highly supersaturated NaCl bulk solution will provide new insights into the detailed mechanism of TSN for many other aqueous solutions.

This work provides evidence for two-step nucleation in highly supersaturated bulk NaCl solution, using electrostatic levitation combined with Raman/X-ray scatterings.

Solvation structure and dynamics of the dimethylammonium cation diluted in liquid water: A molecular dynamics approach

Classical molecular dynamics simulation techniques were employed to investigate the local solvation structure and related dynamics of the dimethylammonium cation diluted in liquid water at ambient conditions. The translational and orientational order around the dimethylammonium cation was investigated in terms of the corresponding radial and angular distribution functions. The results obtained revealed that the first solvation shell of the dimethylammonium consists mainly of two and, less frequently, three water molecules. The two nearest water neighbors form hydrogen bonds with the ammonium hydrogen atoms of the cation, whereas the third neighbor interacts with the methyl hydrogen atoms as well. The distribution of the trigonal order parameter exhibits a bimodal behavior, signifying the existence of local orientational heterogeneities in the solvation shell of the dimethylammonium cation. The calculated continuous and intermittent residence and hydrogen bond lifetimes for the cation-water pairs have also been found to be longer in comparison with the water-water ones. The very similar self-diffusion coefficients of the dimethylammonium cation and the water molecules in the bulk dilute solution indicate that the translational motions of the cation are mainly controlled by the translational mobility of the surrounding water molecules.

Role of Viscosity in Deviations from the Nernst-Einstein Relation

Deviations from the Nernst-Einstein relation are commonly attributed to ion-ion correlation and ion pairing. Despite the fact that these deviations can be quantified by either experimental measurements or molecular dynamics simulations, there is no rule of thumb to tell the extent of deviations. Here, we show that deviations from the Nernst-Einstein relation are proportional to the inverse viscosity by exploring the finite-size effect on transport properties under periodic boundary conditions. This conclusion is in accord with the established experimental results of ionic liquids.

Salt parameterization can drastically affect the results from classical atomistic simulations of water desalination by MoS2 nanopores

Water scarcity is a reality in our world, and scenarios predicted by leading scientists in this area indicate that it will worsen in the next decades. However, new technologies based on low-cost seawater desalination can prevent the worst scenarios, providing fresh water for humanity. With this goal, membranes based on nanoporous materials have been suggested in recent years. One of the materials suggested is MoS₂, and classical Molecular Dynamics (MD) simulation is one of the most powerful tools to explore these nanomaterials. However, distinct force fields employed in MD simulations are parameterized based on distinct experimental quantities. In this paper, we compare two models of salt that were built based on distinct properties of water-salt mixtures. One model fits the hydration free energy and lattice properties, and the second fits the crystal density and the density and the dielectric constant of water and salt mixtures. To compare the models, MD simulations for salty water flow through nanopores of two sizes were used - one pore big enough to accommodate hydrated ions, and one smaller in which the ion has to dehydrate to enter - and two rigid water models from the TIP4P family - TIP4P/2005 and TIP4P/ε. Our results indicate that the water permeability and salt rejection by the membrane are more influenced by the salt model than by the water model, especially for the narrow pore. In fact, completely distinct mechanisms were observed, and they are related to the characteristics employed in the ion model parameterization. The results show that not only can the water model influence the outcomes, but the ion model plays a crucial role when the pore is small enough.

Decelerated Liquid Dynamics Induced by Component-Dependent Supercooling in Hydrogen and Deuterium Quantum Mixtures

Isotopic mixtures of p-H₂ and o-D₂ molecules have been an attractive binary system because they include two kinds of purely isotopic molecules which possess the same electronic potential but the twice different mass inducing differently pronounced nuclear quantum effects (NQEs). Accessing details of structures and dynamics in such quantum mixtures combining complex molecular dynamics with NQEs of different strengths remains a challenging problem. Taking advantage of the nonempirical molecular dynamics method which describes p-H₂ and o-D₂ molecules, we found that the liquid dynamics slows down at a specific mixing ratio, which can be connected to the observed anomalous slowdown of crystallization in the quantum mixtures. We attributed the decelerated dynamics to the component-dependent supercooling of p-H₂ taking place in the mixtures, demonstrating that there is an optimal mixing ratio to hinder crystallization. The obtained physical insights will help in experimentally controlling and achieving unknown quantum mixtures including superfluid.

Surface Energies and Structure of Salt-Brine Interfaces

Permeability of salt formations is controlled by the equilibrium between the salt-brine and salt-salt interfaces described by the dihedral angle, which can change with the composition of the intergranular brine. Here, classical molecular dynamics (MD) simulations were used to investigate the structure and properties of the salt-brine interface to provide insight into the stability of salt systems. Mixed NaCl-KCl brines were investigated to explore differences in ion size on the surface energy and interface structure. Nonlinearity was noted in the salt-brine surface energy with increasing KCl concentration, and the addition of 10% KCl increased surface energies by 2-3 times (5.0 M systems). Size differences in Na⁺ and K⁺ ions altered the packing of dissolved ions and water molecules at the interface, impacting the surface energy. Additionally, ions at the interface had lower numbers of coordinating water molecules than those in the bulk and increased hydration for ions in systems with 100% NaCl or 100% KCl brines. Ultimately, small changes in brine composition away from pure NaCl altered the structure of the salt-brine interface, impacting the dihedral angle and the predicted equilibrium permeability of salt formations.

Shared hydrogen bonds: water in aluminated faujasite

Water confined in faujasite, a zeolite, with aluminium content, exhibits properties different from those of bulk water as well as water confined in siliceous faujasite. The RDF between oxygen of water (OW) and oxygen of aluminium (OAl) shows a prominent first peak near to 2.9 Å similar to any oxygen-oxygen RDF seen in bulk water and unlike water confined in siliceous faujasite. Further, HW-OAl shows a peak near 1.9 Å suggesting hydrogen bonding between hydrogen of water and OAl. The water satisfies the hydrogen bond criteria with both O1Al and O2Al indicating that it is participating in a shared hydrogen bond. The hydrogen bond exchange between such a water forming a shared hydrogen bond to OAl and another water molecule H2Ob is investigated through the changes in the distances and appropriate angles. The O-Al-O angle of the zeolite increases by about 7 degrees on the formation of the shared hydrogen bond. The jump dynamics of the shared hydrogen bond when the two bonds break simultaneously has been obtained and this is reported. This jump reorientation dynamics is different compared to normal hydrogen bonding reported by Laage and Hynes: it has a short lifetime, around 50-100 fs computed from SHB(t). The intermittent and continuous hydrogen bond correlation functions are also reported.

Temperature effects on the ionic conductivity in concentrated alkaline electrolyte solutions

Alkaline electrolyte solutions are important components in rechargeable batteries and alkaline fuel cells. As the ionic conductivity is thought to be a limiting factor in the performance of these devices, which are often operated at elevated temperatures, its temperature dependence is of significant interest. Here we use NaOH as a prototypical example of alkaline electrolytes, and for this system we have carried out reactive molecular dynamics simulations with an experimentally verified high-dimensional neural network potential derived from density-functional theory calculations. It is found that in concentrated NaOH solutions elevated temperatures enhance both the contributions of proton transfer to the ionic conductivity and deviations from the Nernst-Einstein relation. These findings are expected to be of practical relevance for electrochemical devices based on alkaline electrolyte solutions.

Dielectric Decrement for Aqueous NaCl Solutions: Effect of Ionic Charge Scaling in Nonpolarizable Water Force Fields

We investigate the dielectric constant and the dielectric decrement of aqueous NaCl solutions by means of molecular dynamic simulations. We thereby compare the performance of four different force fields and focus on disentangling the origin of the dielectric decrement and the influence of scaled ionic charges, as often used in nonpolarizable force fields to account for the missing dynamic polarizability in the shielding of electrostatic ion interactions. Three of the force fields showed excessive contact ion pair formation, which correlates with a reduced dielectric decrement. In spite of the fact that the scaling of charges only weakly influenced the average polarization of water molecules around an ion, the rescaling of ionic charges did influence the dielectric decrement, and a close-to-linear relation of the slope of the dielectric constant as a function of concentration with the ionic charge was found.

Molecular Dynamics Simulations of the Electrocoalescence Behaviors of Two Unequally Sized Conducting Droplets

The electrocoalescence of droplets plays a crucial role in various fields. However, studies on the effects of droplet radius on the electrocoalescence behaviors of droplets have not been conducted until now. In this work, the electrocoalescence behaviors of two unequally sized conducting droplets are investigated via molecular dynamics (MD) simulations. The influences of electric field strength and droplet radius on the electrocoalescence behaviors of two unequally sized droplets are investigated. When the electric field strength increases, the contact cone angle between the droplets increases, and the two droplets are more likely to partially coalesce and bounce. When the radius of the smaller droplet between the two droplets increases at the same electric field strength, the contact cone angle, daughter droplet size, and ions in the daughter droplet increase, whereas the critical electric field strength ( E_d) for generating the daughter droplet decreases. Furthermore, the daughter droplet is ejected from the smaller droplet when the two droplets have different radii.

Fluid release pressure for nanochannels: the Young-Laplace equation using the effective contact angle

Releasing fluids from nanochannels is quite challenging, yet crucial for the application of nanofluidic systems, e.g. drug delivery and nanoprinting. Previous work suggests that the pressure required to activate the releasing is enormously high (50 to above 300 MPa), while its underlying mechanism still remains unclear. In this work, through molecular dynamics simulations, we have identified a critical tilt angle of the hydrophilic nanochannel, below which spontaneous release of water is achieved. A significant increase in the contact angle is observed during the fluid releasing process due to the transition from the solid fluid contact to the fluid vapor contact. Such transition in nanoscale channels can significantly raise the release pressure by at most ∼30 MPa depending on the channel height and surface property, which makes the classical Young-Laplace equation underestimate the release pressure. By incorporating the derived formula for the largest effective contact angle, a modified Young-Laplace equation is developed, which predicts the release pressure well for both hydrophobic and hydrophilic channels down to the nanoscale. Furthermore, it is discovered that for nanoscale channels, the decreased rate of the normalized release pressure as a function of the contact angle becomes fast when the surface energy of the channel grows strong. The fast decreased rate is mainly caused by the adsorption of water molecules at the exit when the surface becomes highly hydrophilic.

Correlations from Ion Pairing and the Nernst-Einstein Equation

We present a new approximation to ionic conductivity well suited to dynamical atomic-scale simulations, based on the Nernst-Einstein equation. In our approximation, ionic aggregates constitute the elementary charge carriers, and are considered as noninteracting species. This approach conveniently captures the dominant effect of ion-ion correlations on conductivity, short range interactions in the form of clustering. In addition to providing better estimates to the conductivity at a lower computational cost than exact approaches, this new method allows us to understand the physical mechanisms driving ion conduction in concentrated electrolytes. As an example, we consider Li^{+} conduction in poly(ethylene oxide), a standard solid-state polymer electrolyte. Using our newly developed approach, we are able to reproduce recent experimental results reporting negative cation transference numbers at high salt concentrations, and to confirm that this effect can be caused by a large population of negatively charged clusters involving cations.