Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solvation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 721 records · Page 40

Relationship of the Molecular Structure and Transport Properties of Imide-Based Lithium Salts of “Acetonitrile/Water-in-Salt” Electrolytes

“Water-in-salt” (WIS) electrolytes exhibit excellent safety and electrochemical performance. However, they possess high concentrations, relatively low diffusion coefficient, and high viscosity. “Acetonitrile/water in salt” (AWIS) electrolytes can overcome the disadvantages of WIS electrolytes. Under relatively low concentrations, AWIS electrolytes show good electrochemical performance comparable to WIS electrolytes and low conductivity. Herein, we investigate the relationship between the solvation structures and the transport properties using small-angle X-ray scattering and molecular dynamics simulation. We observed two solvation behaviors of AWIS: anions dissolved in acetonitrile forming small acetonitrile/anion clusters and additional water further dissolving the small acetonitrile/anion clusters. The introduction of acetonitrile weakens the water–solute interaction and enhances the cation–anion interaction, which results in an enhanced dynamical slowdown as the concentration increases. In conclusion, this work provides molecular-level understanding of the connection between two-stage solvation structures and transport properties for imide-based lithium salt solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Correlated Ion Transport Governed by Dynamic Local Structure in High Concentration and Localized High Concentration Electrolytes

Understanding the dynamics of cluster formation and network percolation provides the mechanistic link between microscopic solvation structure and transport in concentrated electrolytes, including localized high-concentration electrolytes (LHCEs). Although recent studies have shown LHCEs to form micelle-like aggregates at specific compositions, a quantitative understanding of how solvation structures and transport properties depend on salt–solvent–diluent ratios remains limited. Here, we integrate molecular dynamics with Onsager transport analyses to chart the evolution of solvation microstructure and associated ionic transport in LiFSI/DMC/TTE, achieving good agreement with experimental conductivities across composition. We show that micelle formation arises from a dynamic instability of the cation–anion network, and that network percolation is the primary determinant of conductivity. Two compositional thresholds emerge: A critical network concentration (CNC) and a critical micelle concentration (CMC) that delineate transitions from extended percolating networks to micelle-like clusters and then to fragments. These structural transitions rationalize the nonintuitive conductivity decrease at intermediate dilution despite monotonically decreasing viscosity and provide composition-level design rules for LHCEs.

Mohanakrishnan, Rohith Srinivaas [University of Ca↗

Early-lanthanide(III) acetonitrile–solvento adducts with iodide and noncoordinating anions

Dissolution of LnI 3 (Ln = La, Ce) in acetonitrile (MeCN) results in the highly soluble solvates LnI 3 (MeCN) 5 [Ln = La (1), Ce (2)] in good yield. The ionic complex [La(MeCN) 9 ][LaI 6 ] (4), containing a rare homoleptic La 3+ cation and anion, was also isolated as a minor product. Extending this chemistry to NdI 3 results in the consistent formation of the complex ionic structure [Nd(MeCN) 9 ] 2 [NdI 5 (MeCN)][NdI 6 ][I] (3), which contains an unprecedented pentaiodide lanthanoid anion. Also described is the synthesis, isolation, and structural characterization of several homoleptic early-lanthanide MeCN solvates with noncoordinating anions, namely, [Ln(MeCN) 9 ][AlCl 4 ] 3 [Ln = La (5), Ce (6), Nd (7)]. Notably, complex 6 is the first homoleptic cerium MeCN solvate reported to date. All reported complexes were structurally characterized by X-ray crystallography, as well as by IR spectroscopy and CHN elemental analysis. Furthermore, complexes 1–3 were also characterized by thermogravimetric analysis coupled with mass spectrometry to further elucidate their bulk composition in the solid-state.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Energetics of CO2– in Aqueous Solution

Gas-phase and aqueous solution properties of neutral and anionic clusters of CO2 with 3, 4, and 8 explicit H2O molecules are calculated at the coupled cluster (CCSD(T)) level plus a self-consistent reaction field. Anionic clusters with the radical electron density localized on the carbon of the CO2 molecule rather than localized on the H2O molecules are more favorable energetically by 10-20 kcal/mol in the gas phase (?Hgas(298 K)) and 20-30 kcal/mol in aqueous solution (?Gaq(298 K)). The most favorable structures are those with the largest number of strong hydrogen bonds between the CO2 - and the explicit H2O molecules. Adiabatic electron affinities were calculated in the gas phase and in aqueous solution for the microsolvated anion. The adiabatic electron affinity of aqueous CO2 - is predicted to be 2.35 ± 0.08 eV and is converged with as few as 3 explicit H2O molecules plus a self-consistent reaction field. The EA of aqueous CO2 is significantly greater than the aqueous solvation free energy of the electron. The vertical attachment energies to CO2 and the vertical detachment energies from CO2 - were calculated. The solvated CO2 - anion is substantially bent to 135°, which requires 1.52 eV. The large energy required for bending in combination with the vertical detachment and attachment energies shows that substantial local solvent reorganization occurs on detachment or attachment of an electron to solvated CO2. The formation of aqueous C2O4 2- from CO2 - was also explored, and dimerization is predicted to occur.

McNeill, Ashley↗

Structure of Sulfuric Acid Solutions Using Pair Distribution Function Analysis

Solvation and mesoscale ordering of sulfuric acid and other strong acid solutions leads to suppressed freezing points and strong rheological changes with varying concentration. While the solid-state structures are well-understood, studies focused on the evolving solvation structure in the solution phase have probed a limited concentration range (~1–6 M). This study applies a total scattering approach in both the wide-angle X-ray scattering (WAXS) and pair distribution function (PDF) regimes to elucidate the evolving solvation structure over its full range of acid concentration (0–18 M). The emergence of a prepeak in the WAXS regime at intermediate concentrations indicates a transition from noninteracting sulfate molecules in the dilute limit to sterically limited sulfates at concentrations near its deep eutectic point. Fits to the PDF data quantify this trend, showing a transition from octahedrally hydrated sulfates up to 6–7 M concentrations, followed by gradual dehydration, and eventually reaching a solution structure similar to that of water-in-salt electrolyte systems at high acid concentrations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Molecular Dynamics and Free Energy Calculations of Dicyclohexano-18-crown-6 Diastereoisomers with Sm 2+ , Eu 2+ , Dy 2+ , Yb 2+ , Cf 2+ , and Three Halide Salts in Tetrahydrofuran and Acetonitrile Using the AMOEBA Force Field

With the continual development of lanthanides (Ln) in current technological devices, an efficient separation process is needed that can recover greater amounts of these rare elements. Dicyclohexano-18-crown-6(DCH18C6) is a crown ether that may be a promising candidate for Ln separation, but additional research is required. As such, molecular dynamics(MD) simulations have been performed on four divalent lanthanide halide salts(Sm 2+ , Eu 2+ , Dy 2+ , and Yb 2+ ) and one divalent actinide halide salt (Cf 2+ ) bound to three diastereoisomers of DCH18C6. Dy 2+ , Yb 2+ , Cf 2+ , DCH18C6, and tetrahydrofuran (THF) solvent were parameterized for the AMOEBA polarizable force field for the first time, whereas existing parameters for Sm 2+ and Eu 2+ were utilized from our previous efforts. A coordination number (CN) of six for Ln 2+ /An 2+ –O solvated in THF indicated that the cations interacted almost entirely with the oxygens of the polyether ring. A CN of one for Ln 2+ /An 2+ -N solvated in acetonitrile for systems containing iodide suggested that theN atom of acetonitrile was competitive with I – for cation interactions. Fluctuation between five and six CNs for Dy 2+ and Yb 2+ suggested that although the cations remained in the polyether ring, the size of the ring may not be an ideal fit as these cations possess comparatively smaller ionic radii. Gibbs binding free energies of Sm 2+ in all DCH18C6 diastereoisomers solvated in THF were calculated. The binding free energy of the cis-syn-cis diastereoisomer was the most favorable, followed by cis-anti-cis, and then trans-anti-trans. Lastly, two major types of conformation were observed for each diastereoisomer that were related to the electrostatic interactions and charge density of the cations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling the Behavior of Complex Aqueous Electrolytes Using Machine Learning Interatomic Potentials: The Case of Sodium Sulfate

Understanding the structure and thermodynamics of solvated ions is essential for advancing applications in electrochemistry, water treatment, and energy storage. While ab initio molecular dynamics methods are highly accurate, they are limited by short accessible time and length scales whereas classical force fields struggle with accuracy. Herein, we explore the structure and thermodynamics of complex monovalent-divalent ion pairs using Na 2 SO 4 (aq) as a case study by applying a machine learning interatomic potential (MLIP) trained on density functional theory (DFT) data. Our MLIP-based approach reproduces key bulk properties such as density and radial distribution functions of water. We provide the hydration structure of the sodium and sulfate ions in the 0.1–2 M concentration range and the one-dimensional and two-dimensional potentials of mean force for the sodium–sulfate ion pairing at the low concentration limit (0.1 M), which are inaccessible to DFT. At low concentrations, the sulfate ion is strongly solvated, leading to the stabilization of solvent-separated ion pairs over contact ion pairs. Minimum energy pathway analysis revealed that coordinating two sodium ions with a sulfate ion is a multistep process whereby the sodium ions coordinate to the sulfate ion sequentially. Finally, we demonstrate that MLIPs allow the study of solvated ions beyond simple monovalent pairs with DFT-level accuracy in their low concentration limit (0.1 M) via statistically converged properties from ns-long simulations.

anions↗

Solvent Entropy Contributions to Catalytic Activity in Designed and Optimized Kemp Eliminases

This work reports on the role of solvation for enzymatic catalysis in two distinct, artificially designed Kemp Eliminases, KE07 and KE70, and mutated variants that were optimized by laboratory directed evolution. Using a spatially resolved analysis of hydration patterns, intermolecular vibrations, and local solvent entropies, we identify distinct classes of hydration water and follow their changes upon substrate binding and transition state formation for the designed KE07 and KE70 enzymes and their evolved variants. We observe that differences in hydration of the enzymatic systems are concentrated in the active site and undergo significant changes during substrate recruitment. For KE07, directed evolution reduces variations in the hydration of the polar catalytic center upon substrate binding, preserving strong protein-water interactions, while the evolved enzyme variant of KE70 features a more hydrophobic reaction center for which the expulsion of low-entropy water molecules upon substrate binding is substantially enhanced. While our analysis suggests a system-dependent role of solvation for the substrate binding process, we identify more subtle changes in solvation for the transition state formation, which are less affected by directed evolution.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vibrational Spectroscopy Reveals Effects of Electron Push–Pull and Solvent Polarity on Electron Delocalization in Radical Anions of Donor–Acceptor Compounds

Herein, the nature of excess electrons has been studied in donor–acceptor (D–A) compounds based on substituted triarylamines and a nitrile-functionalized fluorene, by changing the substituents on the triarylamines and also the solvent polarity. We observed that both electron push-pull capability at the distant location in the amine donor unit and solvation in solvents of varying polarity significantly affect the nitrile $ν(C≡N)$ vibrations of the fluorene acceptor unit in radical anions of these D–A compounds. Quantum calculations show that the push-pull capability translates the position of an excess electron while keeping its width relatively constant. On the other hand, solvation changes both, making an electron more compact in polar media. The current study points to the idea that solvation plays a more significant role in controlling the nature of excess electrons, while synthetic modification that influences electron push-pull capability enables further tuning.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Li-Metal Anode in Dilute Electrolyte LiFSI/TMP: Electrochemical Stability Using Ab Initio Molecular Dynamics

Ab initio molecular dynamics simulations were performed for Li + conducting electrolytes based on trimethyl phosphates (TMP) and lithium bis(fluorosulfonyl)imide (Li + FSI – ) salt in contact with a Li-metal electrode. We focused on the transient-state behavior at the electrolyte, interfacial electrolyte–Li-metal electrode, and lithium reference electrode–electrolyte–Li-metal electrode to study dynamics and activation energy barriers of the Li + ion, electrochemical and thermal stability of the interface electrode-electrolyte, and potential behavior of the Li-metal electrode, respectively. Our results show that in the most stable state, Li + ions are tetrahedrally coordinated to three TMPs and one FSI – . The inner solvation shell of a Li-ion is composed of three TMP and one FSI – in one contact ion-pair and four TMPs in a solvent-separated ion-pair. On the other hand, Li-ions transport through electrolyte cages takes place when they are coordinated with three or less molecules that could be a combination of TMPs and FSI – . The decomposition pathway of the LiFSI salt when in direct contact with the Li-metal anode starts with defluorination of FSI – , rapidly losing F– to the lithium surface, forming LiF species. The remaining FSO 2 NSO 2 –2 with the addition of 2e – from the Li-metal decomposes into SO 2 –2 and NFSO 2 –2 . SO 2 –2 deposits on the Li-surface and decomposes into Li 2 O and Li 2 S. The remaining NFSO 2 –2 defluorinates, losing F – ion to the lithium surface, resulting in LiF and the remaining NSO 2 –1 deposits on the lithium surface and decomposes in the following picoseconds, forming several binary compounds such as Li 3 N, Li 2 S and Li 2 O. In contrast, when the salt is solvated by the TMP molecules, avoiding a direct contact with the Li metal electrode, only one defluorination occurs, decomposing the FSI – into FSO 2 NSO 2 –2 and a F – . The two anions remain stable as they are solvated by the TMP molecules. Furthermore, we also analyzed the open circuit potential energy (OCPE) of the Li-metal electrode during the SEI formation. OCPE is calculated from the average local potential profile difference within the Li-metal electrode and a pristine Li-crystal reference electrode (LRE). When no SEI is formed, the Li-metal electrode has an average OCPE of +0.36 eV vs LRE. Due to the formation of a SEI, the Li-metal electrode has an average OCPE between -0.07 and -0.21 eV vs LRE. The OCPE of the Li-metal electrode decreases by ~0.42 eV when a SEI is formed.

25 ENERGY STORAGE↗

Exploration of the Detailed Conditions for Reductive Stability of Mg(TFSI) 2 in Diglyme: Implications for Multivalent Electrolytes

In this paper, we reveal the general mechanisms of partial reduction of multivalent complex cations in conditions specific for the bulk solvent and in the vicinity of the electrified metal electrode surface and disclose the factors affecting the reductive stability of electrolytes for multivalent electrochemistry. Using a combination of ab initio techniques, we clarify the relation between the reductive stability of contact-ion pairs comprising a multivalent cation and a complex anion, their solvation structures, solvent dynamics, and the electrode overpotential. We found that for ion pairs with multiple configurations of the complex anion and the Mg cation whose available orbitals are partially delocalized over the molecular complex and have antibonding character, the primary factor of the reductive stability is the shape factor of the solvation sphere of the metal cation center and the degree of the convexity of a polyhedron formed by the metal cation and its coordinating atoms. We focused specifically on the details of Mg (II) bis(trifluoromethanesulfonyl)imide in diethylene glycol dimethyl ether (Mg(TFSI) 2 )/diglyme) and its singly charged ion pair, MgTFSI + . In particular, we found that both stable (MgTFSI) + and (MgTFSI) 0 ion pairs have the same TFSI configuration but drastically different solvation structures in the bulk solution. This implies that the MgTFSI/dyglyme reductive stability is ultimately determined by the relative time scale of the solvent dynamics and electron transfer at the Mg–anode interface. In the vicinity of the anode surface, steric factors and hindered solvent dynamics may increase the reductive stability of (MgTFSI) + ion pairs at lower overpotential by reducing the metal cation coordination, in stark contrast to the reduction at high overpotential accompanied by TFSI decomposition. By examining other solute/solvent combinations, we conclude that the electrolytes with highly coordinated Mg cation centers are more prone to reductive instability due to the chemical decomposition of the anion or solvent molecules. Finally, the obtained findings disclose critical factors for stable electrolyte design and show the role of interfacial phenomena in reduction of multivalent ions.

36 MATERIALS SCIENCE↗

Cations Enhance Hydride Transfer to Noncatalytic Metals in Concentrated Alkaline Electrolytes

Alkali metal cations are known to influence the kinetics of the hydrogen evolution reaction (HER), acting as either promoters or inhibitors to the rate-determining step depending on the metal surface, local pH, and cation concentration. Despite the importance of concentrated electrolytes for commercial electrochemical cells, the impact of cations on the HER in concentrated alkaline environments (>1 M) and in mixed cation systems remains poorly understood. Here, this study quantifies the HER kinetics at polycrystalline metal surfaces (Pt, Au, Cu, and Fe) and the equilibrium solvation environment in pure and mixed alkali metal hydroxide electrolytes at concentrations up to 3.0 M. Kinetic analyses of Au, Cu, and Fe revealed a positive cation-concentration-effect that was primarily driven by changes to the charge transfer coefficient. Multinuclear NMR spectroscopy examined the solvation of H 2 O/OH – species and the alkali cations as a function of (mixed) alkali cation concentration(s), and demonstrated rapid exchange between solvent, hydroxide, and solvated cations. Together, these findings support models where HER kinetics on noncatalytic metal surfaces in strongly alkaline conditions are primarily governed by the average polarization and polarizability of the metal/solution interface and that increasing cation activity continues to increase the transfer coefficient at metal-hydroxide concentrations up to 3.0 M. Electrolytes and additives which can outcompete weakly hydrated cations and disrupt the interfacial water structure are expected to suppress parasitic HER at electrodes for energy storage and electroplating.

Cations↗

Complex Ion Dynamics in Carbonate Lithium-Ion Battery Electrolytes

Li-ion battery performance is strongly influenced by ionic conductivity, which depends on the mobility of the Li ions in solution, and is related to their solvation structure. In this work, we have performed first-principles molecular dynamics (FPMD) simulations of a LiPF6 salt solvated in different Li-ion battery organic electrolytes. We employ an analytical method using relative angles from successive time intervals to characterize complex ionic motion in multiple dimensions from our FPMD simulations. We find different characteristics of ionic motion on different time scales. We find that the Li ion exhibits a strong caging effect due to its strong solvation structure, while the counterion, PF6– undergoes more Brownian-like motion. Lastly, our results show that ionic motion can be far from purely diffusive and provide a quantitative characterization of the microscopic motion of ions over different time scales.

25 ENERGY STORAGE↗

Quantum Chemical Prediction of p K a Values of Cationic Ion-Exchange Groups in Polymer Electrolyte Membranes

The durability of anion-exchange membranes (AEMs) in alkaline fuel cells is largely determined by the stability of anion-exchange groups. Despite continuous research efforts, the commonly employed cations still have limited stability against hydroxide that can act as a strong base and nucleophile. This work is concerned with base-catalyzed degradation of organic cations initiated by proton abstraction to form reactive ylides or carbenes. We report on the performance of 24 methods combining density functional theory and electronic structure methods with implicit solvation calculations for predicting pK a values of organic cations in water and DMSO. The most accurate computational protocols are obtained using a combination of M06-2X/6-311++G** with the SMD solvation model for water (the mean absolute error of 0.4 pK a units) and B3LYP/aug-cc-pVTZ with the IEFPCM solvation model for DMSO (the mean absolute error of 1.4 pK a units). The aqueous pK a calculation protocol is cross-validated against the experimental C–H acidity constants outside the conventional range of 0–14 pK a values. In conclusion, this study rationalizes alkaline degradation of imidazolium cations with C2-alkyl substituents and provides a theoretical scale of C–H acidity for potential anion-exchange groups in AEMs.

25 ENERGY STORAGE↗

Probing Conformational Evolution and Associated Dynamics of Mg(N(SO 2 CF 3 ) 2 ) 2 ·Dimethoxyethane Adduct Using Solid-State 19 F and 1 H NMR

Bis(trifluoromethanesulfonimide) or TFSI is widely used as a counter anion in electrolyte design due to its structural flexibility and chemical stability. Here we studied the conformational variations and associated dynamics of TFSI in adduct of Mg(TFSI) 2 with dimethoxyethane (DME), a solvate crystalline material using solid-state 1 H and 19 F NMR. TFSI molecular motion in this solvate structure falls within the timescale of the 19 F NMR experiment, yielding spectroscopic signatures for unique TFSI conformers under the coordination environment of Mg 2+ cation. Within the temperature range of -5 to 82 °C, we observe nine distinct TFSI sites in both crystalline and disordered regions using 19 F NMR, reflecting complexity of structural and dynamics of TFS1 anions within solvate structure. The four distinguishable sites in the disordered region for the two CF 3 groups of the same TFSI molecule are identified using chemical shift analysis. The exchange rate constants from site to site are calculated through variable temperature 19 F NMR and two-dimensional (2D) exchange spectroscopy (EXSY) experiments, along with respective activation enthalpies using Eyring's formulation. The flip rate of CF 3 around the S-C bond is estimated as ~ 15 s -1 at 8 °C with ΔH ≠ ~ 22 kJ/mol, but the rotation of the entire TFSI is 4.8 s -1 at 8 °C with a significantly greater ΔH ≠ = 98 ± 10 kJ/mol. Furthermore, the slow conversion of trans to cis conformers at a lower temperature (T ≤ 1 °C) in the crystalline region is monitored, with a conversion rate of ~ 2 x 10 -5 s -1 at -5 °C. Density functional theory (DFT)-based calculations were performed to support further the assignment of experimental chemical shifts, and the activation energy E a = 21.1 kJ/mol obtained for the cis to trans transition is consistent with experimental values. The combined set of 19 F and 1 H under both one-dimensional (1D) and 2D NMR methods demonstrated here can be further used for examining electrode-electrolyte interfaces to probe the motions of various constituents that can enable detailed studies of interfacial processes and dynamics. Ultimately, such studies will aid in the design and discovery of interfacial constructs in which directed defect chemistry, chemical moiety distribution, and nanostructure are employed to drive efficient charge transport.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Condensed Phase Water Molecular Multipole Moments from Deep Neural Network Models Trained on Ab Initio Simulation Data

Ionic solvation phenomena in liquids involve intense interactions in the inner solvation shell. For interactions beyond the first shell, the ion–solvent interaction energies result from the sum of many smaller-magnitude contributions that can still include polarization effects. Deep neural network (DNN) methods have recently found wide application in developing efficient molecular models that maintain near-quantum accuracy. In this work, we extend the DeePMD-kit code to produce accurate molecular multipole moments in the bulk and near interfaces. The new method is validated by comparing the DNN moments with those generated by ab initio simulations. The moments are used to compute the electrostatic potential at the center of a molecular-sized hydrophobic cavity in water. The results show that the fields produced by the DNN models are in quantitative agreement with the AIMD-derived values. These efficient methods will open the door to more accurate solvation models for large solutes such as proteins.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Polarity of Co-solvents Regulates the Charge Storage Mechanisms in Supercapacitors with Concentrated Electrolytes

Developing better energy storage devices depends on comprehending the underlying mechanisms involved in charge storage. With the continuous conception of new electrolytes, this task becomes progressively more urgent and complex. An example is the utilization of co-solvated concentrated solutions. While these show promising electrochemical responses, their dynamic properties (especially under confinement) and their relationships with performance are not fully understood. Here, we combined modified step potential electrochemical spectroscopy and quasielastic neutron scattering to investigate systems composed of activated mesoporous carbon (AMC) and concentrated solutions of lithium bis(trifluoromethanesulfonyl)imide in acetonitrile co-solvated with either toluene or acetone. Further, we report that acetone does not impair surface-controlled mechanisms, contrary to the case with toluene, which competes with charged species to populate the AMC’s pores without contributing to charge storage. In turn, toluene promotes a greater overall capacitance owing to Faradaic processes, which may be related to changes in the solvation structures under confinement.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Local and Global Effects of Dissolved Sodium Chloride on the Structure of Water

Determining how the structure of water is modified by the presence of salts is instrumental to understanding the solvation of biomolecules and, in general, the role played by salts in biochemical processes. However, the extent of hydrogen bonding disruption induced by salts remains controversial. We performed extensive first-principles simulations of solutions of a simple salt (NaCl) and found that, while the cation does not significantly change the structure of water beyond the first solvation shell, the anion has a further reaching effect, modifying the hydrogen-bond network even outside its second solvation shell. We found that a distinctive fingerprint of hydrogen bonding modification is the change in polarizability of water molecules. Molecular dipole moments are instead insensitive probes of long-range modifications induced by Na + and Cl - ions. Though noticeable, the long-range effect of Cl - is expected to be too weak to affect solubility of large biomolecules.

Gaiduk, Alex P.↗