Effect of hydrophilic walls on the hydration of sodium cations in planar nanopores

A computer simulation of the structure of Na⁺ ion hydration shells with sizes in the range of 1 to 100 molecules in a planar model nanopore 0.7 nm wide with structureless hydrophilic walls is performed using the Monte Carlo method at a temperature of 298 K. A detailed model of many-body intermolecular interactions, calibrated with reference to experimental data on the free energy and enthalpy of reactions after gaseous water molecules are added to a hydration shell, is used. It is found that perturbations produced by hydrophilic walls cause the hydration shell to decay into two components that differ in their spatial arrangement and molecular orientational order.     

Charge separation in water molecule clusters under thermal fluctuations: 1. Intermolecular interactions

A detailed model of intermolecular interactions in water molecule clusters is developed that makes it possible to describe their disintegration to ions under conditions of finite temperatures by the stochastic simulation methods. In this model, the Hamiltonian in explicit form includes Coulomb, dispersion, exchange, and polarization interactions; many-particle covalent interactions and hydrogen bonds; the interaction of induced dipoles; charge transfers from ions to molecules; and the recombination of counterion charges, as well as the effect of an ion field on the unpaired interactions of molecules. The model is consistent with experimental data on the free energy and entropy of ion hydration in water vapors and the free energy of the hydration of a recombined ion pair.     

Structure and electric properties of the hydration shell of a singly charged chloride ion in a nanopore with hydrophilic walls

The effect of hydrophilic walls on the structure of the hydration shell of a Cl^? ion is studied in terms of the model flat nanopore in contact with water vapors at room temperature by the Monte Carlo computerassisted simulations. In the field of hydrophilic walls, the hydration shell falls into two parts: the ion-enveloping part and the molecular-film spots spread over the wall surface above and under the ion. Both parts have the pronounced radial-layered structure. The three-dimensional scheme of distribution of the averaged local shell density represents a system of conical coaxial layers expanding in the direction from wall to ion. The effect of forcing out the ion from its own hydration shell is also observed for hydrophilic walls. The specific electric polarizability of the shell is strongly anisotropic. Its longitudinal component is several times larger than the transversal component and behaves nonmonotonically as the hydration shell grows, passing through the maximum. The molecular order near the walls is characterized by the preferential orientation of the molecule plane in parallel to the wall plane and the turn of symmetry axes of molecules in the direction parallel to the normal to the pore plane in the vicinity of the ion.     

Water vapor clustering in the field of Na<Superscript>+</Superscript> cation inside a nanopore with hydrophilic walls. 2. Thermodynamic properties

The bicanonical statistical ensemble method has been used to calculate at the molecular level the free energy, entropy, and work of hydration of single-charged sodium cation in a model planar nanopore with structureless hydrophilic walls. The calculations have been performed in terms of a detailed many-particle model of intermolecular interactions calibrated with respect to experimental data on the free energy and enthalpy of the initial reactions of attachment in water vapor. In contrast to chlorine anion, at initial stages of formation, the hydration shell of sodium cation has a loose chain structure, which is reflected in the character of the interaction with pore walls and the behavior of entropy. Under the conditions of weakly hydrophilic walls, the system loses its stability; however, the stability remains preserved in a pore with strongly hydrophilic walls. Hydrophilic walls stabilize the system and shift the onset of hydration toward lower vapor pressures by several orders of magnitude.     

Structural transition in the OH<Superscript>−</Superscript>(H<Subscript>2</Subscript>O)<Subscript><Emphasis Type="Italic">n</Emphasis></Subscript> cluster in water vapors

The hydration of hydroxyl ion in water vapors at temperatures corresponding to seasonal variations in natural air medium is studied by the Monte Carlo simulation in grand canonical statistical ensemble using the detailed model of intermolecular forces that takes into account many-particle covalent interactions, polarization, and charge transfer. An increase in the number of water molecules in a cluster is accompanied by a structural transition from strongly asymmetric ion environment of water molecules to the formation of enveloping shell composed of these molecules. This transition is accompanied by an abrupt increase in cluster size and qualitative changes in its structural characteristics. The displacement of ion on the surface of clusters with extremely small sizes is an entropy effect. Results of simulation are compared with data on the hydration of hydroxonium at which similar structural transition is not observed, and with data of quantum-chemical calculations.     

Water vapor clustering in the field of Na<Superscript>+</Superscript> cation inside a nanopore with hydrophilic walls. 1. Spatial organization

Computer simulation has been employed to study the structure of a hydration shell of a Na⁺ ion under the conditions of a planar nanopore with structureless hydrophilic walls at 298 K. Intermolecular interactions have been described in terms of a detailed model calibrated with respect to experimental data on the free energy and enthalpy of the initial reactions of vapor molecule attachment to the ion. In the field of hydrophilic walls, the hydration shell is disrupted into an enveloping part and that spread over the surface of the walls. At the final stage of hydration, states with asymmetric distribution of molecules on opposite walls survive and the phenomenon of ion displacement out of its shell is stably reproduced. The orientational molecular order in the system strongly depends on the degree of wall hydrophilicity. The hydration shell of a sodium ion is less stable with respect to disturbances generated by the field of hydrophilic walls than the shell of a chlorine ion is.     

Water vapor nucleation on ion pairs under the conditions of a planar nanopore

Computer simulation has been employed to study the effect of a confined space of a planar model pore with structureless hydrophobic walls on the hydration of Na⁺Cl^– ion pairs in water vapor at room temperature. A detailed many-body model of intermolecular interactions has been used. The model has been calibrated relative to experimental data on the free energy and enthalpy of the initial reactions of water molecule attachment to ions and the results of quantum-chemical calculations of the geometry and energy of Na⁺Cl^– (H₂O)_N clusters in stable configurations, as well as spectroscopic data on Na⁺Cl^– dimer vibration frequencies. The free energy and work of hydration, as well as the adsorption curve, have been calculated from the first principles by the bicanonical statistical ensemble method. The dependence of hydration shell size on interionic distance has been calculated by the method of compensation potential. The transition between the states of a contact (CIP) and a solvent-separated ion pair (SSIP) has been reproduced under the conditions of a nanopore. The influence of the pore increases with the hydration shell size and leads to the stabilization of the SSIP states, which are only conditionally stable in bulk water vapor.     

Structure of the hydration shell of the Na+ ion in a planar nanopore with hydrophobic walls

The effect of steric hindrances in extremely narrow planar pores on the structure of the hydration shell of the single-charged sodium cation in water vapors at room temperature was studied by computer simulation. The deficiency of empty space for the motion in the slit-like pore was shown to slightly affect the radial distribution of molecules around the ion. The integrated (over the directions) numbers of ion-oxygen atom bonds of molecules in the ion’s hydration shell did not change despite the change in the shape of the hydration cluster from three- to two-dimensional. It was concluded that the changes in the positions of molecules relative to the ion were mainly reduced to azimuthal displacements; as a result, the local bulk density of molecules in the pore was higher than at the same distances outside the pore for the same total number of molecules. The distribution of molecules over layers inside the pore demonstrates the effect of molecules spread over the walls. The effect of ion displacement from its own hydration shell found earlier for the free chloride ion is steadily reproduced under the pore conditions. An alternative explanation to this effect was proposed that does not suggest high ion polarizability.     

Hydration of Cl<Superscript>–</Superscript> ion in a planar nanopore with hydrophilic walls. 2. Thermodynamic stability

The Monte Carlo bicanonical statistical ensemble method has been employed to calculate the dependences of the Gibbs free energy, formation work, and entropy on the size of a hydration shell grown from water vapor on single-charged chlorine anion in a model planar nanopore with hydrophilic structureless walls at 298 K. A refined model comprising many-particle polarization interactions and calibrated with respect to experimental data on the free energy and enthalpy of the initial reactions of attachment of water molecules to the ion has been used. It has been found that a weak hydrophilicity of pore walls leads to destabilization of the hydration shell, while a strong one, on the contrary, causes its stabilization. The physical reason for the instability in the field of hydrophilic walls qualitatively differs from that under the conditions of hydration in bulk water vapor.     

On the molecular origins of volumetric data

We use a statistical thermodynamic approach and a simple thermodynamic model of hydration to examine the molecular origins of the volumetric properties of solutes. In this model, solute-solvent interactions are treated as a binding reaction. The free energy of hydration of the noninteracting solute species coincides with the free energy of cavity formation, while the free energy of solute-solvent interactions is given by the binding polynomial. By differentiating the relationship for the free energy of hydration with respect to temperature and pressure, one obtains the complete set of equations describing the thermodynamic profile of hydration, including enthalpy, entropy, volume, compressibility, expansibility, and so forth. The model enables one to rigorously define in thermodynamic terms the hydration number and the related concept of hydration shell, which are both widely used as operational definitions in experimental studies. Hydration number, nh, is the effective number of water molecules solvating the solute and represents the derivative of the free energy of hydration with respect to the logarithm of water activity. One traditional way of studying hydration relies on the use of volumetric measurements. However, microscopic interpretation of macroscopic volumetric data is complicated and currently relies on empirical models that are not backed by theory. We use our derived model to link the microscopic determinants of the volumetric properties of a solute and its statistical thermodynamic parameters. In this treatment, the partial molar volume, V degrees, of a solute depends on the cavity volume, hydration number, and the properties of waters of hydration. In contrast, the partial molar isothermal compressibility, K degrees T, and expansibility, E degrees, observables, in addition to the intrinsic compressibility and expansibility of the cavity enclosing the solute, hydration number, and the properties of waters of hydration, contain previously unappreciated relaxation terms that originate from pressure- and temperature-induced perturbation of the equilibrium between the solvated solute species. If significant, the relaxation terms may bring about a new level of nonadditivity to compressibility and expansibility group contributions that goes beyond the overlap of the hydration shells of adjacent groups. We apply our theoretical results to numerical analyses of the volume and compressibility responses to changes in the distribution of solvated species of polar compounds.     

Mapping Hydration Water around Alcohol Chains by THz Calorimetry

THz spectroscopy was used to probe changes that occur in the dynamics of the hydrogen bond network upon solvation of alcohol chains. The THz spectra can be decomposed into the spectrum of bulk water, tetrahedral hydration water, and more disordered (or interstitial) hydration water. The tetrahedrally ordered hydration water exhibits a band at 195 cm⁻¹ and is localized around the hydrophobic moiety of the alcohol. The interstitial component yields a band at 164 cm⁻¹ which is associated with hydration water in the first hydration shell. These temperature‐dependent changes in the low‐frequency spectrum of solvated alcohol chains can be correlated with changes of heat capacity, entropy, and free energy upon solvation. Surprisingly, not the tetrahedrally ordered component but the interstitial hydration water is found to be mainly responsible for the temperature‐dependent change in ΔC_p and ΔG. The solute‐specific offset in free energy is attributed to void formation and scales linearly with the chain length.     

Thermodynamic characteristics of the hydrate shell of a Na+ ion in a plane nanopore with hydrophobic walls

The chemical potential, free energy, and work of hydration of a single-charged sodium cation are calculated using the Monte Carlo method for a bicanonic statistic ensemble at the molecular level at 298 K in plane model nanopores 0.5 and 0.7 nm wide. It is shown that the nanopores have a stabilizing effect on the hydrate shell of an ion. It is concluded that the crisis of stability that occurs outside a pore is transformed into an abrupt acceleration of growth with the conservation of a stable equilibrium with vapor under the conditions of plane nanopores. It is established that the mechanism of the threshold acceleration of growth inside a pore is associated with an ion being displaced from its own hydrate shell.     

Nucleation of water vapor on Na<Superscript>+</Superscript>Cl<Superscript>−</Superscript> ion pairs: Computer simulation

The Monte Carlo method is used to calculate the free energy, entropy, and work of water cluster formation in the field of Na⁺Cl⁻ ion pairs. A detailed model is used that allows for polarization and covalent many-particle interactions, as well as the effects of ion charge reversal. The model is matched to the experimental data on the free energy of ion hydration and the results of the quantum-chemical calculations of stable configurations. The hydration leads to the cleavage of an ion pair in a molecular cluster after approximately ten water molecules are captured. As vapor molecules are added, the stable interion distance monotonically elongates. The low free energy barrier separating the dissociated and nondissociated states of the ion pair in an equilibrium cluster does not hinders the reversible spontaneous transitions between the states, which are responsible for strong fluctuations and the instability of the system. Unlike hydroxonium-containing ion pairs, the formation of long-lived metastable states of hydrated Na⁺Cl⁻ pairs is impossible.     

The Energy Barrier to Recombination in Hydrated Plasma

An abnormally high potential barrier that separates the H₃O⁺ and Cl^? ions in a cluster of water molecules was revealed. The profile of the barrier was calculated by computer simulation. The calculation was based on a detailed model of intermolecular interactions developed on the basis of experimental data on the free energy and entropy for the addition reactions of water molecules in the vapor phase to the hydration shell of ions in conjunction with the results of quantum-chemical calculations.     

Thermodynamics of thiocyanate complexes of trivalent actinides and lanthanides

The thermodynamic parameters ΔF, ΔH and ΔS of the complexes of Cm(III), C(III) and Tm(III) with the SCN^? ion have been determined at 30°C in ammonium ion medium of unit ionic strength by the temperature variation method. It has been concluded that both the thiocyanate complexes of trivalent actinides and lanthanides are predominantly inner-sphere type. The higher stability of the second complexes of trivalent actinides is reflected either in the enthalpy or the entropy change depending on the degree of hydration of the trivalent actinide ions. The implications of the greater free energy change for PuSCN²⁺ as compared with other trivalent actinide or lanthanide first thiocyanate complexes are discussed.     

Monte Carlo bicanonical ensemble simulation for sodium cation hydration free energy in liquid water

A series of bicanonical ensemble Monte Carlo (BC MC) simulations has been performed to calculate Na⁺ hydration Gibbs energy in aqueous solution. The hydration Gibbs energy of Na⁺ ion in aqueous solution is the difference between formation free energies of Na⁺ (H₂O)_n and (H₂O)_n clusters at n → α. The convergence of the hydration free energy to bulk water value is fast, and the results at n = 60 turned out to be in good agreement with experimental ones and those calculated using free energy perturbation method [1]. The ion–water interaction has been described by Aqvist's pair potential [1] and SPC model [2] has been used for water–water interactions. The behaviour of the absolute Gibbs energy, the entropy, the internal energy of the clusters and the development of hydration shells’ structure with the increase of the number of water molecules are discussed.     

Structure and electric properties of sodium ion hydrate shell in nanopore with hydrophilic walls

The method of molecular–level computer simulation at the temperature of 298 K was used to study the fundamental regularities of formation of electric properties of the hydrate shell of the Na⁺ cation in a planar model nanopore with hydrophilic structureless walls in contact with water vapors. Electric polarizability changes nonmonotonously: as consistent with the changes in the molecular structure of the system. Hydration within the pore occurs in several stages, from formation of chain structures, microdrop compaction and ejection of the ion from its own hydrate shell to encapsulation and absorption of the ion by the solvent preceding formation of nanoelectrolyte. Despite the significant differences in the energy of retaining hydrate shells for Na⁺ and Cl^─ ions, polarizabilities of the two systems are close and behave similarly under variation of conditions. Strong spatial anisotropy of the polarizability tensor of the ion–hydrate complex is due to the effect of the nanopore walls on multiparticle spatial correlations in the system.     

Phenomenon of the ousting of a monatomic ion from its hydration shell in flat nanopores

The structure and stability of hydrate shells of singly charged sodium and chlorine ions are studied by computer simulations under the conditions of nanoscopic flat pores with the use of the previously proposed detailed force field model containing polarization interactions, transferring charge effects as well as manybody interactions of covalent type. It is found that the effect of ousting a monatomic ion from its hydration shell, which has previously been observed by independent authors in bulk vapor, is also reproduced persistently in nanoscopic pores. Whereas the ousting of the ion from its hydration shell in bulk vapor is accompanied by the loss of thermodynamic stability of the system and at sufficiently high vapor pressure causes avalanche-like condensation, under the conditions of a nanoscopic pore the thermodynamic stability is retained. The obtained data show that the ousting of the ion from its hydration shell is a universal phenomenon covering the majority, if not all, of monatomic and, possibly, some of molecular ions.