Ion hydration free energies and water surface potential in water nano drops: The cluster pair approximation and the proton hydration Gibbs free energy in solution

The multi-scale polarizable pseudo-particle solvent coarse-grained approach: From NaCl salt solutions to polyelectrolyte hydration

We discuss key parameters that affect the reliability of hybrid simulations in the aqueous phase based on an efficient multi-scale coarse-grained polarizable pseudo-particle approach, denoted as pppl, to model the solvent water, whereas solutes are modeled using an all atom polarizable force field. Among those parameters, the extension of the solvent domain (SD) at the solute vicinity (domain in which each solvent particle corresponds to a single water molecule) and the magnitude of solute/solvent short range polarization damping effects are shown to be pivotal to model NaCl salty aqueous solutions and the hydration of charged systems, such as the hydrophobic polyelectrolyte polymer that we have recently investigated [Masella et al., J. Chem. Phys. 155, 114903 (2021)]. Strong short range damping is pivotal to simulate aqueous salt NaCl solutions at moderate concentration (up to 1.0M). The SD extension (as well as short range damping) has a weak effect on the polymer conformation; however, it plays a pivotal role in computing accurate polymer/solvent interaction energies. As the pppl approach is up to two orders of magnitude computationally more efficient than all atom polarizable force field methods, our results show it to be an efficient alternative route to investigate the equilibrium properties of complex charged molecular systems in extended chemical environments.

Simulating core electron binding energies of halogenated species adsorbed on ice surfaces and in solution via relativistic quantum embedding calculations

In this work, we investigate the effects of the environment on the X-ray photoelectron spectra of hydrogen chloride and chloride ions adsorbed on ice surfaces, as well as of chloride ions in water droplets. In our approach, we combine a density functional theory (DFT) description of the ice surface with that of halogen species using the recently developed relativistic core-valence separation equation of motion coupled cluster (CVS-EOM-IP-CCSD) via the frozen density embedding formalism (FDE), to determine the K and L_1,2,3 edges of chlorine. Our calculations, which incorporate temperature effects through snapshots from classical molecular dynamics simulations, are shown to reproduce experimental trends in the change of the core binding energies of Cl^- upon moving from a liquid (water droplets) to an interfacial (ice quasi-liquid layer) environment. Our simulations yield water valence band binding energies in good agreement with experiment, which vary little between the droplets and the ice surface. For halide core binding energies there is an overall trend for overestimating experimental values, though good agreement between theory and experiment is found for Cl^- in water droplets and on ice. For HCl on the other hand there are significant discrepancies between experimental and calculated core binding energies when we consider structural models that maintain the H-Cl bond more or less intact. An analysis of models that allow for pre-dissociated and dissociated structures suggests that experimentally observed chemical shifts in binding energies between Cl^- and HCl would require that H⁺ (in the form of H₃O⁺) and Cl^- are separated by roughly 4-6 Å.

Hybrid polarizable simulations of a conventional hydrophobic polyelectrolyte. Toward a theoretical tool for green science innovation

Hybrid modeling approaches based on all-atom force fields to handle a solute and coarse-grained models to account for the solvent are promising numerical tools that can be used to understand the properties of large and multi-components solutions and thus to speed up the development of new industrial products that obey the standard of green and sustainable chemistry. Here, we discuss the ability of a full polarizable hybrid approach coupled to a standard molecular dynamics scheme to model the behavior in the aqueous phase and at infinite dilution conditions of a standard hydrophobic polyelectrolyte polymer whose charge is neutralized by explicit counterions. Beyond the standard picture of a polyelectrolyte behavior governed by an interplay between opposite intra-polyelectrolyte and inter-polyelectrolyte/counterion Coulombic effects, our simulations show the key role played by both intra-solute polarization effects and long range solute/solvent electrostatics to stabilize compact globular conformations of that polyelectrolyte. Our full polarizable hybrid modeling approach is thus a new theoretical tool well suited to be used in digital strategies for accelerating innovation for green science, for instance.

Absolute ion hydration free energy scale and the surface potential of water via quantum simulation

With a goal of determining an absolute free energy scale for ion hydration, quasi-chemical theory and ab initio quantum mechanical simulations are employed to obtain an accurate value for the bulk hydration free energy of the Na⁺ ion. The free energy is partitioned into three parts: 1) the inner-shell or chemical contribution that includes direct interactions of the ion with nearby waters, 2) the packing free energy that is the work to produce a cavity of size λ in water, and 3) the long-range contribution that involves all interactions outside the inner shell. The interfacial potential contribution to the free energy resides in the long-range term. By averaging cation and anion data for that contribution, cumulant terms of all odd orders in the electrostatic potential are removed. The computed total is then the bulk hydration free energy. Comparison with the experimentally derived real hydration free energy produces an effective surface potential of water in the range -0.4 to -0.5 V. The result is consistent with a variety of experiments concerning acid-base chemistry, ion distributions near hydrophobic interfaces, and electric fields near the surface of water droplets.

Transport Models of Ammonium Nitrogen in Wastewater from Rare Earth Smelteries by Reverse Osmosis Membranes

Wastewater from rare earth smelteries contains large amounts of ammonium nitrogen (NH4+-N), which causes severe environmental problems. In this contribution, the desalination efficiency of reverse osmosis (RO) was investigated in the treatment of NH4Cl or NaCl solutions from 0.1 to 40 g/L under different operating pressures with a commercial RO membrane. Experimental results showed that when an operating pressure above 30 bar is applied to the 5 g/L NH4Cl solution, the permeate was found to meet the discharge standards of NH4+-N. Compared to NH4Cl, the permeate fluxes of NaCl solutions were higher due to the higher net driving force and lower propensity to membrane fouling. Theoretical models indicate a linear relationship between water flux and the net driving force for both NH4Cl and NaCl solutions. On the contrary, a power function between the salt flux and concentration difference correlated well with the experimental data for salt transport. The equations for water and salt transport obtained by this work would provide a facile and practical means for predicting the membrane performance in design and optimization of RO processes for the treatment of wastewater from the rare earth industry.

de Almeida M, Marques JMC, Prudente FV. A thermodynamic view on the microsolvation of ions by rare gas: application to Li+ with argon

Parallel tempering Monte Carlo calculations on the Li⁺Ar_n microsolvation clusters have shown that the two peaks appearing in the heat capacity curve as a function of temperature correspond to the melting of the second and first solvation shells.