No Arabic abstract
Although aqueous electrolytes are among the most important solutions, the molecular simulation of their intertwined properties of chemical potentials, solubility and activity coefficients has remained a challenging problem, and has attracted considerable recent interest. In this perspectives review, we focus on the simplest case of aqueous sodium chloride at ambient conditions and discuss the two main factors that have impeded progress. The first is lack of consensus with respect to the appropriate methodology for force field (FF) development. We examine how most commonly used FFs have been developed, and emphasize the importance of distinguishing between Training Set Properties used to fit the FF parameters, and Test Set Properties, which are pure predictions of additional properties. The second is disagreement among solubility results obtained, even using identical FFs and thermodynamic conditions. Solubility calculations have been approached using both thermodynamic--based methods and direct molecular dynamics--based methods implementing coexisting solution and solid phases. Although convergence has been very recently achieved among results based on the former approach, there is as yet no general agreement with simulation results based on the latter methodology. We also propose a new method to directly calculate the electrolyte standard chemical potential in the Henry-Law ideality model. We conclude by making recommendations for calculating solubility, chemical potentials and activity coefficients, and outline a potential path for future progress.
Accurate prediction of a gas solubility in a liquid is crucial in many areas of chemistry, and a detailed understanding of the molecular mechanism of the gas solvation continues to be an active area of research. Here, we extend the idea of constant chemical potential molecular dynamics (C{mu}MD) approach to the calculation of the gas solubility in the liquid under constant gas chemical potential conditions. As a representative example, we utilize this method to calculate the isothermal solubility of carbon dioxide in water. Additionally, we provide microscopic insight into the mechanism of solvation that preferentially occurs in areas of the surface where the hydrogen network is broken.
The Seeebeck coefficients of the non-aqueous electrolytes tetrabutylammonium nitrate, tetraoctylphosphonium bromide and tetradodecylammonium nitrate in 1-octanol, 1-dodecanol and ethylene-glycol are measured in a temperature range from T=30 to T=45 C. The Seebeck coefficient is generally of the order of a few hundreds of microvolts per Kelvin for aqueous solution of inorganic ions. Here we report huge values of 7 mV/K at 0.1M concentration for tetrabutylammonium nitrate in 1-dodecanol. These striking results open the question of unexpectedly large kosmotrope or structure making effects of tetraalkylammonium ions on the structure of alcohols.
Molecular dynamics simulations of aqueous electrolytes generally rely on empirical force fields, combining dispersion interactions - described by a truncated Lennard-Jones (LJ) potential - and electrostatic interactions - described by a Coulomb potential computed with a long-range solver. Recently, force fields using rescaled ionic charges (electronic continuum correction, ECC), possibly complemented with rescaling of LJ parameters (electronic continuum correction rescaled, ECCR), have shown promising results in bulk, but their performance at interfaces has been less explored. Here we started by exploring the impact of the LJ potential truncation on the surface tension of a sodium chloride aqueous solution. We show a discrepancy between the numerical predictions for truncated LJ interactions with a large cutoff and for untruncated LJ interactions computed with a long-range solver, which can bias comparison of force field predictions with experiments. Using a long-range solver for LJ interactions, we then show that an ionic charge rescaling factor chosen to correct long-range electrostatic interactions in bulk also describes accurately image charge repulsion at the liquid-vapor interface, and that the rescaling of LJ parameters in ECCR models - aimed at capturing local ion-ion and ion-water interactions in bulk - also describes well the formation of an ionic double layer at the liquid-vapor interface. Overall, these results suggest that the molecular modeling of aqueous electrolytes at interfaces would benefit from using long-range solvers for dispersion forces, and from using ECCR models, where the charge rescaling factor should be chosen to correct long-range electrostatic interactions.
Aqueous zinc-air batteries (ZABs) are a low-cost, safe, and sustainable technology for stationary energy storage. ZABs with pH-buffered near-neutral electrolytes have the potential for longer lifetime compared to traditional alkaline ZABs due to the slower absorption of carbonates at non-alkaline pH values. However, existing near-neutral electrolytes often contain halide salts, which are corrosive and threaten the precipitation of ZnO as the dominant discharge product. This paper presents a method for designing halide-free aqueous ZAB electrolytes using thermodynamic descriptors to computationally screen components. The dynamic performance of a ZAB with one possible halide-free aqueous electrolyte based on organic salts is simulated using an advanced method of continuum modeling, and the results are validated by experiments. XRD, SEM, and EDS measurements of Zn electrodes show that ZnO is the dominant discharge product, and operando pH measurements confirm the stability of the electrolyte pH during cell cycling. Long-term full cell cycling tests are performed, and RRDE measurements elucidate the mechanism of ORR and OER. Our analysis shows that aqueous electrolytes containing organic salts could be a promising field of research for zinc-based batteries, due to their Zn$^{2+}$ chelating and pH buffering properties. We discuss the remaining challenges including the electrochemical stability of the electrolyte components.
Forced detachment of a single polymer chain, strongly-adsorbed on a solid substrate, is investigated by two complementary methods: a coarse-grained analytical dynamical model, based on the Onsager stochastic equation, and Molecular Dynamics (MD) simulations with Langevin thermostat. The suggested approach makes it possible to go beyond the limitations of the conventional Bell-Evans model. We observe a series of characteristic force spikes when the pulling force is measured against the cantilever displacement during detachment at constant velocity $v_c$ (displacement control mode) and find that the average magnitude of this force increases as $v_c$ grows. The probability distributions of the pulling force and the end-monomer distance from the surface at the moment of final detachment are investigated for different adsorption energy $epsilon$ and pulling velocity $v_c$. Our extensive MD-simulations validate and support the main theoretical findings. Moreover, the simulation reveals a novel behavior: for a strong-friction and massive cantilever the force spikes pattern is smeared out at large $v_c$. As a challenging task for experimental bio-polymers sequencing in future we suggest the fabrication of stiff, super-light, nanometer-sized AFM probe.