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.     

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.     

Hydration of Cl<Superscript>–</Superscript> ion in a planar nanopore with hydrophilic walls. 1. Molecular structure

Computer simulation has been employed to obtain equilibrium molecular configurations, as well as spatial and angular distributions of water molecules, under the action of the field of a single-charged chlorine anion in a model planar nanopore with structureless walls at room temperature. A detailed many-body model of intermolecular interactions calibrated in accordance with experimental data relative to the free energy of hydration in water vapor has been used. The effect of the hydrophilicity of the walls on the ion hydration shell consists in its disintegration into two parts, i.e., molecules retained exclusively due to the interactions with the ion and those adsorbed on the walls. In the regime of strong interactions with the walls, two relatively stable states arise with asymmetric distribution of molecules between opposite walls. The existence of the two metastable states destabilizes the position of ions inside a pore and is expected to accelerate their adsorption on the walls.     

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.     

Water vapor clustering in the field of a chlorine anion occurring in a planar nanopore with structureless walls

The Monte Carlo bicanonical statistical ensemble method has been employed to calculate the free energy, entropy, and work of Cl^? ion hydration in model planar pores 0.5 and 0.7 nm wide at 298 and 400 K. A detailed model of many-body interactions with the ion has been used, the model being matched to experimental data with respect to the free energy and enthalpy of attachment reaction in water vapor. Under the conditions of a restricted volume, the equilibrium size of a hydration shell substantially decreases, with the effect becoming stronger in the range of moderate and large sizes. In moderately supersaturated vapors, under the conditions of a nanopore, the ion loses its hydration shell as the temperature is decreased. In supersaturated vapors, the hydration shell formed on the ion is thermodynamically stable, while the stability crisis shifts to the region of larger sizes. The enhancement of the thermodynamic stability in the pore results from a rise in the chemical potential of molecules due to the deficiency of closet neighbors and a reduction in the entropy under the conditions of the restricted volume. As the temperature is elevated, the effect of ion displacement out of its hydration shell is leveled. The regularities derived in terms of the estimation model based on the capillary approximation are in qualitative agreement with the results of computer simulation.     

Ion pairs in aqueous electrolyte microdrops under conditions of a flat nanopore

The molecular mechanisms of aqueous solvent penetration into a flat nanopore with hydrophobic structureless walls containing a Na⁺Cl^? ion pair with nonfixed distance between ions is studied by computer simulations. A detailed many-body polycenter model of intermolecular interactions calibrated with respect to experimental data for the free energy of attachment of water vapor molecules and quantum-chemical calculations in clusters is used. The ion pair hydration results in its decomposition. Drawing the molecules into the gap between ions makes easier penetration of solvent and filling of the nanopore with electrolyte. The ion-pair dissociation is accompanied by dramatic changes in the chemical potential of molecules and electric properties of the whole system. The thermodynamic characteristics of decomposition are stable as regards variations in the pore width. The post-decomposition electric polarizability demonstrates strong anisotropy associated with the nanopore flatness.     

Charge separation in Na<Superscript>+</Superscript>Cl<Superscript>-</Superscript>(H<Subscript>2</Subscript>O)<Subscript><Emphasis Type="Italic">n</Emphasis></Subscript> clusters in water vapors. 1. Intermolecular interactions

Charge separation in “soft” nanoparticles composed of water molecules, as well as sodium and chlorine ions, is studied by computer simulation. The detailed model of intermolecular interactions that includes, in addition to Coulomb, exchange, and dispersion forces, many-particle polarization and covalent interactions, as well as the effect due to the transfer of excess ion charges and influence of ion field on molecular interactions, is constructed. Model potentials are calibrated using experimental data on the free energy and enthalpy of the addition of vapor molecules to the hydration shells of ions, as well as the data of quantum-chemical calculations for stable cluster configurations and the vibration frequency of interionic bonds. The allowance for many-particle interactions makes it possible to improve the agreement between experimental and quantum-chemical data by more than an order of magnitude. The disregard for many-particle interactions leads to the significant overestimation of cluster stability.     

Speciation of the curium(III) ion in aqueous solution: a combined study by quantum chemistry and molecular dynamics simulation

The structures of aquo complexes of the curium(III) ion have been systematically studied using quantum chemical and molecular dynamics (MD) methods. The first hydration shell of the Cm3+ ion has been calculated using density functional theory (DFT), with and without inclusion of the conductor-like polarizable continuum medium (CPCM) model of solvation. The calculated results indicate that the primary hydration number of Cm3+ is nine, with a Cm-O bond distance of 2.47-2.48 A. The calculated bond distances and the hydration number are in excellent agreement with available experimental data. The inclusion of a complete second hydration shell of Cm3+ has been investigated using both DFT and MD methods. The presence of the second hydration shell has significant effects on the primary coordination sphere, suggesting that the explicit inclusion of second-shell effects is important for understanding the nature of the first shell. The calculated results indicate that 21 water molecules can be coordinated in the second hydration shell of the Cm3+ ion. MD simulations within the hydrated-ion model suggest that the second-shell water molecules exchange with the bulk solvent with a lifetime of 161 ps.     

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.     

Charge separation in Na<Superscript>+</Superscript>Cl<Superscript>-</Superscript>(H<Subscript>2</Subscript>O)<Subscript><Emphasis Type="Italic">n</Emphasis></Subscript> clusters in water vapors. 2. Free energy

Charge separation in “soft” nanoparticles composed of water molecules, as well as sodium and chlorine ions, is studied by computer simulation. The detailed model of intermolecular interactions that includes, in addition to Coulomb, exchange, and dispersion forces, many-particle polarization and covalent interactions, as well as the effect due to the transfer of excess ion charges and influence of ion field on molecular interactions, is constructed. Model potentials are calibrated using experimental data on the free energy and enthalpy of the addition of vapor molecules to the hydration shells of ions, as well as the data of quantum-chemical calculations for stable cluster configurations and the vibration frequency of interionic bonds. The allowance for many-particle interactions makes it possible to improve the agreement between experimental and quantum-chemical data by more than an order of magnitude. The disregard for many-particle interactions leads to the significant overestimation of cluster stability.     

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.     

Mean force potential of interaction between Na<Superscript>+</Superscript> and Cl<Superscript>−</Superscript> ions in planar nanopores in contact with water under pressure

The mean force potential (MFP) of interaction between counterions Na⁺ and Cl^? in a planar nanopore with structureless hydrophobic walls is calculated via computer simulation under the condition that the nanopore is in contact with water at an external pressure that exceeds the saturation pressure but remains insufficient to fill the nanopore with water. For a nanopore with a liquid phase, the MFP dependence on the interionic distance indicates the dissociation of an ion pair into two hydrated ions in a nanopore that is not completely filled with water. Fluctuations in the number of water molecules drawn into the interionic space decisively influence the dissociation. The attraction between counterions, averaged over thermal fluctuations, depends largely on the pore width and grows as the shielding of the ions’ electric field by water molecules in a narrow pore diminishes. The contributions from energy and entropy to the free energy of hydration are analyzed.     

Nature of hydration shells of a polyoxy‐anion with a large cationic centre: The case of iodate ion in water

The structural nature of the solvation shells of an iodate ion, which is known to be a polyoxy‐anion with a large cationic centre, is investigated by means of Born–Oppenheimer molecular dynamics (BOMD) simulations using BLYP and the dispersion corrected BLYP‐D3 functionals. The iodate ion is found to have two distinct solvation regions around the positively charged iodine (iodine solvation shell or ISS) and the negatively charged oxygens (oxygen solvation shell or OSS). We have looked at the spatial, orientational, and hydrogen bond distributions of water in the two solvation regions. It is found that the water orientational profile in the ISS is typical of a cation hydration shell. The hydrogen bonded structure of water in the OSS is found to be very similar to that of the bulk water structure. Thus, the iodate ion essentially behaves like a positively charged iodine ion in water as if there is no anionic part. This explains why the cationic character of the iodate ion was prominently seen in earlier studies. The arrangement of water molecules in the two solvation shells and in the intervening regions around the iodate ion is further resolved by looking at structural cross‐correlations. The electronic properties of the solvation shells are also looked at by calculating the solute–solvent orbital overlap and dipole moments of the solute and solvation shell water. We have also performed BOMD simulations of iodate ion‐water clusters at experimentally relevant conditions. The simulation results are found to be in agreement with experimental results. © 2018 Wiley Periodicals, Inc.     

High Temperature Stability of Hydrated Ion Pairs Na<Superscript>+</Superscript>Cl<Superscript>–</Superscript>(H<Subscript>2</Subscript>O)<Subscript><Emphasis Type="Italic">N</Emphasis></Subscript> under Conditions of a Flat Nanopore

The high-temperature stability of hydrated ion pairs under conditions of a nanoscopic flat pore with hydrophobic structureless walls is studied by computer simulations. The limited space of the nanopore stimulates dissociation of the contact ion pair (CIP) with its transition to the state of the solvent-separated ion pair (SSIP); moreover, the ion pair demonstrates a high degree of stability on heating. The inverse temperature effect where the heating renders a moderate consolidating effect on the state of a hydrated contact ion pair is observed: when heated to the electrolyte boiling point, the free energy barrier that separates the CIP and SSIP states shifts by 2 molecules towards the larger hydration shells. On the pressure scale, the boundary between CIP and SSIP states shifts at the same rate as the saturating pressure with the increase in the temperature.     

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.     

Car-Parrinello molecular dynamics simulation of Fe 3+ (aq)

The optimized geometry and energetic properties of Fe(D2O)n 3+ clusters, with n = 4 and 6, have been studied with density-functional theory calculations and the BLYP functional, and the hydration of a single Fe 3+ ion in a periodic box with 32 water molecules at room temperature has been studied with Car-Parrinello molecular dynamics and the same functional. We have compared the results from the CPMD simulation with classical MD simulations, using a flexible SPC-based water model and the same number of water molecules, to evaluate the relative strengths and weaknesses of the two MD methods. The classical MD simulations and the CPMD simulations both give Fe-water distances in good agreement with experiment, but for the intramolecular vibrations, the classical MD yields considerably better absolute frequencies and ion-induced frequency shifts. On the other hand, the CPMD method performs considerably better than the classical MD in describing the intramolecular geometry of the water molecule in the first hydration shell and the average first shell...second shell hydrogen-bond distance. Differences between the two methods are also found with respect to the second-shell water orientations. The effect of the small box size (32 vs 512 water molecules) was evaluated by comparing results from classical simulations using different box sizes; non-negligible effects are found for the ion-water distance and the tilt angles of the water molecules in the second hydration shell and for the O-D stretching vibrational frequencies of the water molecules in the first hydration shell.