Effect of short and long range forces on the structure of water. II. Orientational ordering and the dielectric constant

Structure of Simple Dipolar Water-Like Fluids: Primitive Model and Hard Tetrahedra

Dipolar versions of two qualitatively different types of simple short range model fluids which exhibit the phenomenon of hydrogen bonding and which could thus serve as a reference in equations of state for associating fluids have been considered: the primitive model of water descending from the TIP4P model and the fluid of hard tetrahedra. The hydrogen bonding structure exhibited by the latter model results from purely repulsive interactions whereas in the first model the "hydrogen bonding interaction" is explicitly incorporated in the model. Since the water molecules bear a strong dipole moment, the effect of the added dipole-dipole interaction on the structure of the two short-range models is therefore examined considering them both in the full and screened dipole-dipole modifications. It is found that the hydrogen bonding structure in the primitive model resulting from the site-site interactions is so strong that the additional dipole-dipole interaction has only a marginal effect on its structure and contributes thus only to the internal energy. On the contrary, even only a weak dipole-dipole interaction destroys the original hydrogen bonding structure of the hard tetrahedron fluid; to preserve it, a screened dipole-dipole interaction has to be used in the equation of state development.

The long-range convergence of the energetic properties of the water monomer in bulk water at room temperature

The Interacting Quantum Atoms (IQA) energy partitioning scheme has been applied to a set of liquid water largely spherical clusters (henceforth called spheres) of up to 9 Å radius, with a maximum cluster size of 113 molecules. This constitutes half of the commonly used 216 molecules in a typical simulation box of a liquid water box, and to our knowledge is the largest analysis of this kind ever undertaken. As well as demonstrating the topological analysis of large systems, which has only recently become computationally feasible, important long range properties of liquid water are obtained. The full topological partitioning of each sphere into atomic basins is used to consider the long-range convergence of the energetic and multipolar properties of the water molecule at the centre of each sphere. It is found that the total molecular energy converges to its 9 Å value after 7 Å, which corresponds to approximately the first three solvation shells, while the molecular dipole and quadrupole moments approximately converge after 5.5 Å, which corresponds to approximately the first two solvation shells. The effect of water molecule flexibility is also considered.

Fluctuations of local electric field and dipole moments in water between metal walls

We examine the thermal fluctuations of the local electric field Ek (loc) and the dipole moment μk in liquid water at T = 298 K between metal walls in electric field applied in the perpendicular direction. We use analytic theory and molecular dynamics simulation. In this situation, there is a global electrostatic coupling between the surface charges on the walls and the polarization in the bulk. Then, the correlation function of the polarization density pz(r) along the applied field contains a homogeneous part inversely proportional to the cell volume V. Accounting for the long-range dipolar interaction, we derive the Kirkwood-Fröhlich formula for the polarization fluctuations when the specimen volume v is much smaller than V. However, for not small v/V, the homogeneous part comes into play in dielectric relations. We also calculate the distribution of Ek (loc) in applied field. As a unique feature of water, its magnitude |Ek (loc)| obeys a Gaussian distribution with a large mean value E0 ≅ 17 V/nm, which arises mainly from the surrounding hydrogen-bonded molecules. Since |μk|E0 ∼ 30kBT, μk becomes mostly parallel to Ek (loc). As a result, the orientation distributions of these two vectors nearly coincide, assuming the classical exponential form. In dynamics, the component of μk(t) parallel to Ek (loc)(t) changes on the time scale of the hydrogen bonds ∼5 ps, while its smaller perpendicular component undergoes librational motions on time scales of 0.01 ps.

Extended excluded volume: Its origin and consequences

In contrast to the common intuitive/speculative approach based on an analysis of thermodynamic or structural data of (nonpolar) fluids, the statistical mechanical approach is used to extend the excluded volume concept to all other types of fluids. The (extended) excluded volume incorporates, in addition to common nonelectrostatic interactions defining the shape and size of the molecules, also the short-range part of the repulsive interactions between the embedded Coulombic sites. In this study we show that the extended excluded volume concept correctly predicts the behavior of the partial molar volume (PMV) at infinite dilution in different solvents and, particularly, differences between nonpolar and associating solvents. The concept is then applied to estimate the PMV of methanol in water.

Water and aqueous solutions: simple non-speculative model approach

Different ways of molecular modeling of water are analyzed and their similarities and differences identified. An up-to-date summary of achievements of a general approach to common rigid site-site interaction models of molecular fluids applied to water and aqueous solutions is then presented and discussed. The method is based on considering only a short-range part of a total realistic potential (such as SPC/E or TIPxP) which determines the structure of water (and fluids in general). A simplification of the interactions at short intermolecular separations leads then to simple models, called primitive models. Quite accurate results in an analytic form for the thermodynamic properties of the models are obtained using the thermodynamic perturbation theory. It is shown that the properly constructed primitive models reproduce, qualitatively, anomalies of pure water and basic characteristics of hydrophobic hydration. The concept of an extended excluded volume, based on pseudo-hard bodies, is introduced and exemplified by considering the partial molar volume of apolar solutes. Potential future development towards a theory of water based on the primitive models as a reference with the long-range contributions added as a perturbation is discussed.

Long-distance correlations in molecular orientations of liquid water and shape-dependent hydrophobic force

Liquid water, at ambient conditions, has short-range density correlations which are well known in literature. Surprisingly, large scale molecular-dynamics simulations reveal an unusually long-distance correlation in "longitudinal" part of dipole-dipole orientational correlations. It is nonvanishing even at 75 A and falls off exponentially with a correlation length of about 24 A beyond solvation region. Numerical evidence suggests that the long-range nature of dipole-dipole correlation is due to underlying fluctuating network of hydrogen bonds in the liquid phase. This correlation is shown to give a shape dependant attraction between two hydrophobic surfaces at large distances of separation and the range of this attractive force is in agreement with experiments. In addition it is seen that quadrupolar fluctuations vanish within the first solvation peak (3 A) .

Perturbed Chain-Statistical Associating Fluid Theory Extended to Dipolar and Quadrupolar Molecular Fluids

The perturbed chain statistical associating fluid theory (PC-SAFT) is extended to polar molecular fluids, namely dipolar and quadrupolar fluids. The extension is based on the perturbation theory for polar fluids by Stell and co-workers. Appropriate expressions are proposed for dipole-dipole, quadrupole-quadrupole, and dipole-quadrupole interactions. Furthermore, induced dipole interactions are calculated explicitly in the model. The new polar PC-SAFT model is relatively complex; for this purpose, a truncated polar PC-SAFT model is proposed using only the leading term in the polynomial expansion for polar interactions. The new model is used for the calculation of thermodynamic properties of various quadrupolar pure fluids. In all cases, the agreement between experimental data and model predictions is very good.

Modeling of the Carbon Dioxide Solubility in Imidazolium-Based Ionic Liquids with the tPC-PSAFT Equation of State

In this work, an equation of state (EoS) is developed to predict accurately the phase behavior of ionic liquid + CO2 systems based on the truncated perturbed chain polar statistical associating fluid theory (tPC-PSAFT) EoS. This EoS accounts explicitly for the dipolar interactions between ionic liquid molecules, the quadrupolar interactions between CO2 molecules, and the Lewis acid-base type of association between the ionic liquid and the CO2 molecules. Physically meaningful model pure-component parameters for ionic liquids are estimated based on literature data. All experimental vapor-liquid equilibrium data are correlated with a single linearly temperature-dependent binary interaction parameter. The ability of the model to describe accurately carbon dioxide solubility in various 1-alkyl-3-methylimidazolium-based ionic liquids with different alkyl chain lengths and different anions at pressures from 0 to 100 MPa and carbon dioxide fractions from 0 to 75 mol % is demonstrated. In all cases, good agreement with experimental data is obtained.

Structural Correlations and Motifs in Liquid Water at Selected Temperatures: Ab Initio and Empirical Model Predictions

To gain further insight into liquid water's structure and explore the role of different physical forces underlying the interaction between water molecules, the radial and angular structure of water is probed as a function of temperature for a carefully selected set of theoretical models. Simulations are performed with empirical rigid, empirical polarizable, empirical flexible with classical and quantum nuclei, and ab initio models with classical nuclei at 300 and 353 K and quantum nuclei at 300 K. The predicted radial distribution functions, spatial maps, and angular distributions of the neighboring water molecules are consistent with a model of liquid water in which water molecules are tetrahedrally coordinated. In addition, three-dimensional joint distribution functions are introduced and analyzed. By comparison of the functions obtained for hexagonal ice to those of liquid water, several thermally disordered, ice-like cluster structures are identified in the fluid.

Effect of the Range of Interactions on the Properties of Fluids. 2. Structure and Phase Behavior of Acetonitrile, Hydrogen Fluoride, and Formic Acid

To complete the study on the effect of the long-range part of Coulombic interactions on properties of complex polar and associating fluids, we have investigated in detail three compounds with extreme features: acetonitrile for its unusually large dipole moment, hydrogen fluoride with very strong hydrogen bonding, and formic acid for its potential formation of different n-mers in liquid and gaseous phases. The effect of the long-range Coulombic interactions on both the structure and thermodynamics of the homogeneous phase, and on the vapor-liquid equilibria has been examined using the same decomposition of realistic potential models into a short-range part and a residual part as in the previous paper [Kettler, M.; et al. J. Phys. Chem. B 2002, 106, 7537-7546]. The present results fully confirm the previous findings that the properties of polar and associating systems are determined primarily by the short-range interactions regardless of their nature, i.e., contributions arising from the long-range interactions constitute only a small portion of the total properties, and thus that the short-range potential counterpart of full realistic models can be used as a convenient reference for a successful perturbation expansion.

Angular resolution and range of dipole–dipole correlations in water

We investigate the dipolar correlations in liquid water at angular resolution by molecular-dynamics simulations of a large periodic simulation system containing about 40,000 molecules. Because we are particularly interested in the long-range ordering, we use a simple three-point model for these molecules. The electrostatics is treated both by Ewald summation and by minimum image truncation combined with a reaction field approach. To gain insight into the angular dependence of the simulated dipolar ordering we introduce a suitable expansion of the molecular pair distribution function into a set of two-dimensional correlation functions. We show that these functions enable detailed insights into the shell structure of the dipolar ordering around a given water molecule. For these functions we derive analytical expressions in the particular case in which liquid water is conceived as a dielectric continuum. Comparisons of these continuum models with the correlation functions derived from the simulations yield the key result that liquid water behaves like a continuum dielectric beyond distances of about 15 angstroms from a given water molecule. We argue that this should be a generic property of water independent of our modeling. By comparison of the results of the two different electrostatics treatments with the continuum description we show that the boundary artifacts occurring in both methods are isotropically distributed and are locally small in the respective boundary regions.