Hydration of nonelectrolytes in binary aqueous solutions

Literature data on the thermodynamic properties of binary aqueous solutions of nonelectrolytes that show negative deviations from Raoult’s law due largely to the contribution of the hydration of the solute are briefly surveyed. Attention is focused on simulating the thermodynamic properties of solutions using equations of the cluster model. It is shown that the model is based on the assumption that there exists a distribution of stoichiometric hydrates over hydration numbers. In terms of the theory of ideal associated solutions, the equations for activity coefficients, osmotic coefficients, vapor pressure, and excess thermodynamic functions (volume, Gibbs energy, enthalpy, entropy) are obtained in analytical form. Basic parameters in the equations are the hydration numbers of the nonelectrolyte (the mathematical expectation of the distribution of hydrates) and the dispersions of the distribution. It is concluded that the model equations adequately describe the thermodynamic properties of a wide range of nonelectrolytes partly or completely soluble in water.     

高温及超临界条件下 MgCl2与CaCl2水溶液中离子水化与缔合的量子化学计算

对NaCl等碱金属水溶液的研究表明,室温条件下,离子在溶液中以水合形式存在,而在高温及超临界时,阴阳离子将结合成为离子对．采用量子化学计算,研究了MgCl2与CaCl2水溶液中水化与缔合的情况．通过Gaussian98软件包计算了阳离子的水化自由能以及离子对的生成能,从而获得水合离子与离子对的热力学稳定性及其随温度、压力的变化情况．通过比较热力学稳定性,考察了两种溶液中水化与缔合的变化情况．研究结果表明,MgCl2与CaCl2水溶液中离子水化与缔合的变化趋势与碱金属溶液基本一致,但是存在一个过渡区域,该区域内离子对与水合离子共存,因此需要采用不同于碱金属溶液的方法处理MgCl2与CaCl2水溶液．     

Calculating the thermodynamic properties of aqueous solutions of alkali metal carboxylates

A modified Robinson-Stokes equation with terms that consider the formation of ionic hydrates and associates is used to describe thermodynamic properties of aqueous solutions of electrolytes. The model is used to describe data on the osmotic coefficients of aqueous solutions of alkali metal carboxylates, and to calculate the mean ionic activity coefficients of salts and excess Gibbs energies. The key contributions from ionic hydration and association to the nonideality of solutions is determined by analyzing the contributions of various factors. Relations that connect the hydration numbers of electrolytes with the parameters of the Pitzer-Mayorga equation and a modified Hückel equation are developed.     

Activity coefficients of concentrated strong and weak electrolytes by a hydration equilibrium and group contribution model

A new speciation-based group contribution model for activity coefficients is proposed to estimate the equilibrium properties of aqueous solutions containing electrolytes. The chemical part of the model accounts for the hydration equilibrium of water and ions with the formation of ion n-water complexes in a single stage process; the hydration number n and the hydration equilibrium constant K are the two independent parameters in this part. The physical part of the model is the UNIFAC group contribution model for short-range interactions. Each ion is considered as a group. Long-range interactions are accounted for by a Pitzer contribution (Debye–Hückel theory). The model is compared with experimental data at 25 °C including water activity, osmotic coefficients, activity coefficients, and pH of binary diluted and concentrated electrolyte solutions (up to 20 mol kg⁻¹ for NaOH, 16 mol kg⁻¹ for HCl, etc.).     

Intradiffusion coefficients for perchlorate ions in zinc perchlorate and zinc chloride solutions at 25°C: Comparing transport properties of zinc chloride and zinc perchlorate systems

Intradiffusion coefficients for³⁶ClO ₄ ^– have been measured in solutions of zinc perchlorate of concentration 0.1 to 3 mol dm^–3 at 25°C by the diaphragm cell technique. In addition, intradiffusion coefficients for perchlorate ions in zinc chloride solutions have been measured over a concentration range at 25°C. The results confirm previous work on the effect of complexation on diffusion in zinc chloride solutions above a salt concentration of 0.1M. The present data, together with literature data for diffusion coefficients of the other species present in the zinc perchlorate electrolyte system, have enabled a simple analysis of the hydration around the zinc ions to be carried out. This indicates that the water diffusion data are consistent with the zinc ions having an effective hydration sphere of 11 (±2) water molecules. This is in keeping with values obtained for other simple divalent electrolytes using the same model. The model is extended here to allow analysis of water diffusion in zinc chloride solutions taking into account the presence of complexed chloro-zinc species. The experimental data are consistent with the effective hydration of the chloro-zinc complexes being independent of the number of chloride ligands and equal to 18±3 over a concentration range of 0 tol mol-dm^–3. This postulate is discussed in terms of its consequences on the water ligand dynamics for the complex equilibria.     

Diffusion and Ion Association in Concentrated Solutions of Aqueous Lithium, Sodium, and Potassium Sulfates

Taylor dispersion is used to measure mutual diffusion coefficients for aqueous Li₂SO₄ solutions at concentrations from 0.09 to 2.62 mol-dm^-3 at 25°C. The Li₂SO₄ results and previously reported diffusion coefficients for aqueous Na₂SO₄ and K₂SO₄ are compared with predictions made by treating the limiting electrolyte diffusion coefficients as reference values and applying corrections for nonideal solution behavior, ionic hydration, and viscosity changes as the concentration is raised. Good agreement is obtained if the M⁺ + SO ₄ ^2- ? MSO ₄ ^- (M = Li, Na, K) association equilibria are included in the analysis. Extents of formation of the MSO ₄ ^- ion pairs are evaluated by fitting Pitzer's mixed electrolyte equations for aqueous M⁺–MSO ₄ ^- –SO ₄ ^2- ions to osmotic coefficient data. Diffusion coefficients for hypothetical solutions of the completely dissociated M₂SO₄ electrolytes are calculated to illustrate the effects of ion association on diffusion. Association of the M⁺ and SO ₄ ^2- ions increases the overall mobility and thermodynamic driving forces for their diffusion.     

Structure of the nearest surrounding of the Na<Superscript>+</Superscript> ion in aqueous solutions of its salts

Published data obtained by various research methods on structural characteristics of sodium ion hydration in aqueous solutions of its salts and authors, X-ray diffraction data have been generalized. Structural parameters of the nearest surrounding of Na⁺ ion, such as its coordination number, interparticle distances, and types of ion association, have been discussed. It has been noted that the coordination number of the cation changes from four to six upon dilution of the solutions.     

过饱和五硼酸钠溶液结构

用快速X射线衍射法测量了298和323K时过饱和五硼酸钠溶液的时间空间平均结构,得到了溶液径向分布函数.通过模型设计及理论计算获得了B-B、B-O、O-O、Na-O、Na-B原子对相互作用的理论偏径向分布函数,并讨论了浓度和温度对五硼酸钠溶液结构的影响.在五硼酸钠过饱和溶液中六水合Na+形成八面体结构,其平均配位数随浓度和温度变化不大,作用距离随温度升高及浓度减小而减小;给出了溶液中主要硼酸盐离子B3O3(OH)4-、B5O6(OH)4-和B(OH)3的水合结构,较高温度及较高浓度有利于更高聚合的多聚硼酸根离子的形成;浓度对硼酸盐离子的第一水合层的水合数影响较大,在较浓的五硼酸钠过饱和溶液中,五硼酸根离子的一个端氧单齿配位到Na+上形成离子对,Na-B特征距离为0.328nm.     

Thermodynamic studies of ionic interactions in aqueous solutions of imidazolium-based ionic liquids [Emim][Br] and [Bmim][Cl

Experimental measurements of density at different temperatures ranging from 293.15 to 313.15 K, the speed of sound and osmotic coefficients at 298.15 K for aqueous solution of 1-ethyl-3-methylimidazolium bromide ([Emim][Br]), and osmotic coefficients at 298.15 K for aqueous solutions of 1-butyl-3-methylimidazolium chloride ([Bmim][Cl]) in the dilute concentration region are taken. The data are used to obtain compressibilities, expansivity, apparent and limiting molar properties, internal pressure, activity, and activity coefficients for [Emim][Br] in aqueous solutions. Experimental activity coefficient data are compared with that obtained from Debye-Hückel and Pitzer models. The activity data are further used to obtain the hydration number and the osmotic second virial coefficients of ionic liquids. Partial molar entropies of [Bmim][Cl] are also obtained using the free-energy and enthalpy data. The distance of the closest approach of ions is estimated using the activity data for ILs in aqueous solutions and is compared with that of X-ray data analysis in the solid phase. The measured data show that the concentration dependence for aqueous solutions of [Emim][Br] can be accounted for in terms of the hydrophobic hydration of ions and that this IL exhibits Coulombic interactions as well as hydrophobic hydration for both the cations and anions. The small hydration numbers for the studied ILs indicate that the low charge density of cations and their hydrophobic nature is responsible for the formation of the water-structure-enforced ion pairs.     

Interfacial compositions of cationic and mixed non-ionic micelles by chemical trapping: a new method for characterizing the properties of amphiphilic aggregates

A specially synthesized arenediazonium ion bound to amphiphilic aggregates decomposes spontaneously via rate determining loss of N₂ to give a highly reactive, unselective, aryl cation intermediate. This intermediate is trapped competitively by weakly basic nucleophiles in the interfacial region of aggregates such as micelles and other association colloids. Product yields, analyzed by HPLC with UV detection, are used to estimate, simultaneously, the interfacial concentrations of a number of different nucleophiles, including water, that are commonly found at the surfaces of biomembranes and in many commercial products. Two applications of the method are discussed. First, we show that the interfacial concentrations of X⁻ (X=Br, Cl) increase steadily with increasing cetyltrimethylammonium halide (CTAX) and tetramethylammonium halide (TMAX) concentrations and that the interfacial concentrations of these counterions increase continuously with their aqueous phase concentrations at a constant degree of micelle ionization. Interfacial Br⁻ and Cl⁻ concentrations also show marked increases at their respective sphere-to-rod transitions. This steady increase in interfacial counterion concentration with increasing aqueous counterion concentration contradicts a basic assumption of the pseudophase ion exchange (PIE) model of chemical reactivity in aggregates, i.e. that the total concentrations of ions at aggregate interfaces is constant and independent of the amphiphile and salt concentrations. The consequences for the PIE model are discussed. Second, the chemical trapping reaction is used to estimate: (a) distributions of terminal OH groups of non-ionic amphiphiles in mixed non-ionic micelles composed of amphiphiles with different lengths of oligoethylene oxide chains and (b) hydration numbers of the inner layers of interfacial region next to the hydrocarbon core in these mixed micelles. Terminal OH groups distributions are well fitted by a radial one-dimensional random walk model. The average hydration number for the inner layers at 40°C is about 3, in agreement with estimates from NMR water (D₂O) self-diffusion measurements and with the hydration number of 3 for aqueous solutions of polyethylene oxide. The results suggest that the hydration states of the ethylene oxide (EO) units near the micellar core are near their minimum value. Recent and potential applications of the chemical trapping method are briefly discussed.     

Transfer of Liquid Cathode Components to the Gas Phase and Their Effect on the Parameters of the Atmospheric Pressure DC Discharge

The transfer of solvent and components of dissolved substances from aqueous solutions to the gas phase under the action of atmospheric pressure air glow discharge was experimentally studied. solutions of NaCl, KCl, CuCl₂, MgCl₂, CaCl₂, SrCl₂, BaCl₂, NaNO₃, KNO₃, Ba(NO₃)₂, Na₂SO₄, K₂SO₄, CuSO₄ with concentrations of 0.1–0.5 mol/L were used as the cathodes of DC discharge at a current of 10–70 mA. The influence of the solution composition on cathode voltage drop and the electric field strength in plasma has been shown. Plasma emission spectra showed the appearance of metal atoms in the plasma requires threshold discharge current or threshold power input to the liquid cathode by ion bombardment. The threshold power values depend on the mass of hydrated cations and their concentration in solution. The efficiency of the transfer processes was characterized by transfer coefficients—the number of particles transferred from the liquid to the gas phase per one ion bombarding the cathode. Dependences of the transfer coefficients on the power dissipated in the cathode region and on the hydration energy of the cations were obtained. Experimental data on the rate of condensate accumulation in the special trap were used to estimate the concentrations of water molecules in the plasma.     

Thermodynamics of aqueous solutions of the alkali metal sulfates

The available thermodynamic properties for aqueous solutions of each of the alkali metal sulfates have been combined and analyzed within the framework of the ion interaction model at temperatures up to 225°C. It was necessary to set ₁ equal to 1.4kg^1/2-mol^–1/2 in order to obtain a satisfactory fit. The temperature dependence of the ion interaction parameters was given the functional form used by Rogers and Pitzer⁽¹⁾ in their study of Na₂SO₄(aq). With few exceptions, it was possible to reproduce the available thermodynamic data for aqueous solutions of the alkali metal to within the estimated experimental error. Thermodynamic results for Na₂SO₄(aq) appear to be adequate in this temperature range, but enthalpy and heat capacity data for the other alkali metal sulfate solutions are conspicuously lacking. Activity coefficients of these electrolytes decreased to less than 0.1 at moderate molalities at the higher temperatures, and their order changed with increasing temperature; two results which could be due to a combination of hydration and association effects.Research sponsored by the Division of Chemical Sciences, Office of Basic Energy Sciences of the U.S. Department of Energy under contract DE-AC05-840R21400 with the Martin Marietta Energy Systems, Inc.     

A theoretical study of the hydration of Rb+ by Monte Carlo simulations with refined ab initio-based model potentials

An infinitely diluted aqueous solution of Rb⁺ was studied using ab initio-based model potentials in classical Monte Carlo simulations to describe its structural and thermodynamic features. An existing flexible and polarizable model [Saint-Martin et al. in J Chem Phys 113(24) 10899, 2000] was used for water–water interactions, and the parameters of the Rb⁺–water potential were fitted to reproduce the polarizability of the cation and a sample of ab initio pair interaction energies. It was necessary to calibrate the basis set to be employed as a reference, which resulted in a new determination of the complete basis set (CBS) limit energy of the optimal Rb⁺–OH₂ configuration. Good agreement was found for the values produced by the model with ab initio calculations of three- and four-body nonadditive contributions to the energy, as well as with ab initio and experimental data for the energies, the enthalpies and the geometric parameters of Rb⁺(H₂O) _n clusters, with n = 1,  2,…, 8. Thus validated, the potential was used for simulations of the aqueous solution with three versions of the MCDHO water model; this allowed to assess the relative importance of including flexibility and polarizability in the molecular model. In agreement with experimental data, the Rb⁺–O radial distribution function (RDF) showed three maxima, and hence three hydration shells. The average coordination number was found to be 6.9, with a broad distribution from 4 to 12. The dipole moment of the water molecules in the first hydration shell was tilted to 55° with respect to the ion’s electric field and had a lower value than the average in bulk water; this latter value was recovered at the second shell. The use of the nonpolarizable version of the MCDHO water model resulted in an enhanced alignment to the ion’s electric field, not only in the first, but also in the second hydration shell. The hydration enthalpy was determined from the numerical simulation, taking into account corrections to the interfacial potential and to the spurious effects due to the periodicity imposed by the Ewald sums; the resulting value lied within the range of the various different experimental data. An analysis of the interaction energies between the ion and the water molecules in the different hydration shells and the bulk showed the same partition of the hydration enthalpy as for K⁺. The reason for this similarity is that at distances longer than 3 Å, the ion–water interaction is dominated by the charge-(enhanced) dipole term. Thus, it was concluded that starting at K⁺, the hydration properties of the heavier alkali metal cations should be very similar.     

A neutron inelastic scattering investigation of the diffusion kinetics of H2O molecules and hydration complexes in concentrated ionic solutions

Intermolecular frequencies of H₂O's and the diffusion kinetics have been investigated by neutron inelastic scattering for concentrated ionic solutions containing small and/or multiply charged cations (e.g., Cr⁺³, Mg⁺², Ca⁺², and Li⁺¹). As higher concentrations are approached such that the majority of H₂O's are in hydration layers, their exchange time can exceed the neutron interaction time. Then the diffusion kinetics depart functionally from activated reorientations of individual H₂O's characteristic of lower concentrations and evolve to continuous diffusion processes of hydration complexes characterized by small self-diffusion coefficients. The general features of the observed evolution in the functionality of the diffusion kinetics are found to be functionally consistent with an approximate model which includes contributions from the delayed diffusional exchange of individual H₂O's as well as the continuous diffusion of hydrated ions. At a given concentration, the temperature interval over which this evolution in functionality occurs increases both with increasing strength of the primary cation-H₂O coordination and with anion basicity. Further, as the temperature decreases, frequencies of defined cation-water hydration complexes gradually sharpen in a continuous manner, showing no abrupt variations at glass transitions. Anions of increasing basicity decrease the self-diffusion coefficients of the ion-water complexes and perturbed frequencies characteristic of cation-water hydration complexes. Such anion effects, at high concentrations, correspond to an increasing degree of time-average indirect or direct ion pairing with increasing anion basicity. This results, in turn, both in a distortion or partial disruption of the cation hydration sheaths and in a degree of coupling and/or bridging between anions and hydrated cations so as to increase the effective masses and friction coefficients associated with their diffusional motions.     

Modeling aqueous electrolyte solutions. Part 2. Weak electrolytes

In this work the ePC-SAFT model is applied to weak electrolytes, such as weak acids or salts that do form ion pairs. Considering an association/dissociation equilibrium accounts for the fact that the electrolytes are not fully dissociated. Applying this approach, modeling the mean ionic activity coefficients (MIAC) as well as the water activity coefficients (WAC) is in very good agreement with experimental data for the aqueous HF system as well as for solutions of cadmium halides or alkali acetates. Experimental MIACs of ZnBr₂ and ZnI₂ reveal the formation of more than one complex in aqueous solutions. Implementing a simultaneous two-step ion-pairing mechanism also allows the modeling of the MIAC of these zinc salts in water.     

Ion association in lithium metaborate solution: a Raman and ab initio insight

Ion association and hydration clusters in aqueous lithium borate solution are extremely important to understand some extraordinary properties of lithium borates. In the present work, polyborate distribution in aqueous LiBO₂ solution was investigated through Raman and thermodynamics equilibrium analysis. Geometry and stability of hydrated clusters LiB(OH)₄(H₂O)_n up to n = 8 were calculated at the B3LYP/aug-cc-pVDZ level. Three different types of ion association, namely, contact ion pairs (CIP), solvent-shared ion pairs (SIP) and solvent separated ion pairs (SSIP) were obtained; characteristics of all of these stable configurations were determined, and the most stable hydrated clusters were chosen. Then the mechanisms of ion aggregation and crystal nuclei formation in the LiB(OH)₄ solution were proposed. The tight four-hydrated sphere of Li⁺ makes it difficult for the dehydrated form of its first hydration sphere to from a CIP, which is the passible reason that lithium borate always has a large super-saturation degree.     

Activities of the components of glycerol-water binary solutions at 298.15 K

Based on the assumption that there exists a distribution of hydrates over the hydration number, expressions for the hydration number, activity coefficients of the components, and excess Gibbs energy were obtained. It was demonstrated that the van Laar formula for the activity coefficients corresponds to the Poisson distribution of hydrates. It was established that, at a constant ratio of the variance to the mathematical expectation of the distribution of hydrate, the model’s equations adequately describe the available experimental data on the vapor pressure and water activity for the glycerol-water system over the entire concentration range.