Spectroscopic study of water-NaCl-benzene mixtures at high temperatures and pressures

Near-infrared and ultraviolet spectra of water-NaCl-benzene mixtures have been measured in the 473-573 K and 100-400 bar range and 373-498 K and 50-300 bar range, respectively. Concentrations of water in the benzene-rich phase and benzene in the water-rich phase were estimated from integrated intensities of the absorption bands. It is found that addition of NaCl in the aqueous phase suppresses transfer of water into the benzene-rich phase, and the relative decrease in water solubility in benzene exhibits good correlation with an increase in density of the aqueous NaCl solution relative to that of neat water. The salting-out constant for the water-NaCl-benzene system, which is estimated from a relative decrease in benzene solubility in the aqueous phase by addition of sodium chloride, increases significantly with increasing temperature. It is suggested that the effect of sodium chloride on the water-benzene mutual solubilities can be explained by ion-induced electrostriction of the aqueous phase.     

Computer simulation investigation of the water-benzene interface in a broad range of thermodynamic States from ambient to supercritical conditions

The dependence of the properties of the water-benzene system on the thermodynamic conditions in a broad range of temperatures and pressures has been investigated by computer simulation methods. For this purpose, Monte Carlo simulations have been performed at 23 different thermodynamic states, ranging from ambient to supercritical conditions. The density profiles of the water and benzene molecules have been determined at each of the thermodynamic states investigated. Information on the dependence of the mutual solubility of the two components in each other as well as of the width of the interface on the temperature and pressure has been extracted from these profiles. The width of the interface has been found to increase with increasing temperature up to a certain point, where it diverges. The temperature of this divergence corresponds to the mixing of the two phases. The determination of the critical mixing temperature at various pressures allowed us to estimate the upper critical curve, separating the two-phase and one-phase liquid systems, of the phase diagram of the simulated water-benzene system. In analyzing the preferential orientation of the interfacial molecules relative to the interface, it has been found that the main orientational preference of the benzene molecules is to lie parallel with the plane of the interface, and the water molecules penetrated deepest into the benzene phase prefer to stay perpendicular to the interface, pointing by one of their O-H bonds almost straight toward the benzene phase, whereas the waters located at the aqueous side of the interface are preferentially aligned parallel with the interfacial plane. Although the strength of the observed orientational preferences decreases rapidly with increasing temperature, the preferred orientations themselves are found to be independent of the thermodynamic conditions. Remains of the orientational preferences of the molecules are found to be present up to temperatures as high as 650 K. The analysis of the relative orientation of the neighboring water-benzene pairs has revealed that the radius of the first hydration shell of the benzene molecules is independent of the thermodynamic conditions, even if the system consists of one single phase. It has been found that the nearest water neighbors of the benzene molecules are preferentially located above and below the benzene ring, whereas more distant water neighbors, belonging still to the first hydration shell, prefer to stay within the plane of the benzene molecule. In the two-phase systems the dipole vector of the nearest waters has been found to be preferentially perpendicular to the vector pointing from the center of the benzene molecule to the water O atom.     

Dynamical and structural properties of benzene in supercritical water

We have employed an anisotropic united atom model of benzene (R. O. Contreras, Ph.D. thesis, Universitat Rovira i Virgili 2002) that reproduces the quadrupolar moment of this molecule through the inclusion of seven point charges. We show that this kind of interaction is required to reproduce the solvation of these molecules in supercritical water. We have computed self-diffusion coefficient and Maxwell-Stefan coefficients as well as the shear viscosity for the mixture water-benzene at supercritical conditions. A strong density and composition dependence of these properties is observed. In addition, our simulations are in qualitative agreement with the experimental evidence that, at medium densities (0.6 g/cm(3) and 673 K), almost half of the benzene molecules have one hydrogen bond with water molecules. We also observe that these bonds are longer lived than the corresponding hydrogen bonds between water molecules. Similarly, we obtain an important reduction of the dielectric constant of the mixture with the increment of the amount of benzene molecules at medium and high densities.     

The liquid water-benzene system

The 500 MHz NMR spectra of water-benzene solution near saturation at 303.15, 323.15, and 343.15 K indicate that there is a proton-proton exchange between the water and benzene molecules. In the solution water appears to be present as a dimer attached to the benzene pi cloud on one side of each of the two (initially degenerate) fundamental energy levels, as predicted by the Jahn-Teller effect. This view is reinforced by the fact that one of its hydrogen atoms hovers above one of the carbon atoms and the other three are spread upward around the C6 axis of the benzene molecule. It is also supported by the calculated NMR spectra. Both effects are responsible for the change in the NMR spectra of the water molecules from a single line into four AB signals.     

Near infrared study of water-benzene mixtures at high temperatures and pressures

Near-infrared absorption of water-benzene mixtures has been measured at temperatures and pressures in the ranges of 473-673 K and 100-400 bar, respectively. Concentrations of water and benzene in the water-rich phase of the mixtures were obtained from the integrated absorption intensities of the OH stretching overtone transition of water and the CH stretching overtone transition of benzene, respectively. Using these concentrations, the densities of the water-rich phase were estimated and compared with the average densities before mixing, which were calculated from literature densities of neat water and neat benzene. It is found that anomalously large volume expansion on the mixing occurs in the region enclosed by an extended line of the three-phase equilibrium curve and the one-phase critical curve of the mixtures, and the gas-liquid equilibrium curve of water. Furthermore, magnitude of the relative volume change increases with decreasing molar fraction of benzene in the present experimental range. It is suggested that dissolving a small amount of benzene in water induces a change in the fluid density from a liquidlike condition to a gaslike condition in the vicinity of the critical region.     

Volumetric behavior of water–methanol mixtures in the vicinity of the critical region

Densities of water–methanol mixtures at 573 and 588 K and at pressures in the 100–200 bar range have been measured with a vibrating-tube densimeter. Temperature and pressure dependence of the excess molar volumes together with the previous results was discussed. A large negative-to-positive sigmoidal change of the excess molar volumes as a function of methanol mole fraction was interpreted on the basis of an estimated critical locus of the mixtures. The volumetric behavior of the mixtures was compared with that of the previously reported water–benzene mixtures by estimating the relative volume change on mixing. A large negative volume change at the lower methanol concentrations is in sharp contrast to the large positive change for the water–benzene mixtures. This contrast may be attributable to characteristic features of aqueous solutions of hydrophilic and hydrophobic substances in the vicinity of the critical region. The behavior of the water–methanol mixtures at the lower methanol mole fractions was discussed in terms of the local solute–solvent structure by estimating radial distribution functions and self-diffusion coefficients from molecular dynamics calculations.     

Spectroscopic study of mutual solubilities of water and benzene at high temperatures and pressures

Near-infrared and ultraviolet absorption of water-benzene mixtures has been measured at temperatures and pressures in the ranges of 323-673 K and 50-400 bar, respectively. Concentrations of water and benzene in both the water-rich phase and the benzene-rich phase of the mixtures were obtained from absorption intensities of near-infrared bands of water and benzene and ultraviolet bands of benzene. Mutual solubilities in molar fractions increase remarkably with increasing temperature at pressures in the two-liquid-phase coexistence region, and are consistent with previously reported values. It proves that the solubility of benzene in water is an order of magnitude smaller than that of water in benzene throughout the two-phase region. In addition, it is found that effect of pressure on the solubilities is opposite between water in benzene and benzene in water. These solubility properties are discussed on the basis of a cavity-based solvation model. It is suggested that the asymmetry in the mutual solubility and the opposite direction of the pressure effect are caused by difference in molecular size and difference in thermal compressibility, respectively, between water and benzene.     

Structure of liquid chlorobenzene in the temperature range 293-363 K

The classical molecular dynamics method is used to simulate the structure of liquid chlorobenzene in the temperature range 293-363 K. Theradial angular distribution functions for the distances between the benzene molecular planes and the angle between them; the radial distribution functions for the distances between chlorine atoms; self-diffusion coefficients, local dipole moments, and permittivities are calculated. In the entire temperature range molecules are joined into agglomerates due to contacts between chlorine atoms (halogen aggregation) and to the specific interactions of benzene rings, which causes mainly parallel and perpendicular orientations of the rings in the first coordination sphere. The molecules in the agglomerates are mostly organized in the 1D motif: chains of chlorine atoms and stacks of benzene rings. With increasing temperature, the agglomerate structure is reorganized, with the most visible changes occurring in temperature ranges 293-298 K, 313-323 K, and 343-353 K.     

Translational and rotational energy content of benzene molecules IR-desorbed from an in vacuo liquid surface

Benzene molecules were desorbed from an in vacuo aqueous liquid beam by direct irradiation of the beam with an IR laser tuned to the 2.85 μm absorption band of water. Spectroscopic interrogation of the desorbed benzene molecules was performed via 1 + 1 Resonance-Enhanced Multi-photon Ionisation (REMPI). Rotational contour analyses of the 6 vibronic transition of benzene were performed to determine the rotational temperature of those molecules ejected during the desorption event. At the peak of the desorption plume density, the rotational temperatures were found to be up to ～100 K lower than that recorded for molecules spontaneously evaporating from the liquid surface. At longer IR-UV laser delay times the benzene rotational temperatures are found to return to those observed following spontaneous evaporation. No evidence of IR desorbed neutral or cationic benzene-containing clusters was observed. However, ionic clusters were observed to be formed after REMPI of the benzene monomer. Analysis of the benzene intensity and that of post-REMPI formed clusters as a function of IR-UV delay shows that number density and local translational temperature vary along the desorption plume.     

Analysis of the pressure effect on the local composition in a water-alkanol mixture using Kirkwood-Buff integrals

Kirkwood-Buff integrals are calculated from the thermodynamic data for binary mixtures of water with methanol, ethanol, 1-propanol, and 2-propanol at a temperature of 298.15 K in the pressure range from atmospheric to 100 MPa. The values of local compositions Δn _ij are calculated which characterize the excess (or deficit) of molecules i around the central molecule j. It is found that the pressure affects destructively the homoassociation in all mixtures studied. In a series MeOH < EtOH < 2-PrOH < 1-PrOH an excess of molecules around the similar type molecules increases in the local environment and the pressure effect on the local composition is enhanced.     

Solvation free energies and hydration structure of N-methyl-p-nitroaniline

Solvation Gibbs energies of N-methyl-p-nitroaniline (MNA) in water and 1-octanol are calculated using the expanded ensemble molecular dynamics method with a force field taken from the literature. The accuracy of the free energy calculations is verified with the experimental Gibbs free energy data and found to reproduce the experimental 1-octanol∕water partition coefficient to within ±0.1 in log unit. To investigate the hydration structure around N-methyl-p-nitroaniline, an independent NVT molecular dynamics simulation was performed at ambient conditions. The local organization of water molecules around the solute MNA molecule was investigated using the radial distribution function (RDF), the coordination number, and the extent of hydrogen bonding. The spatial distribution functions (SDFs) show that the water molecules are distributed above and below the nitrogen atoms parallel to the plane of aromatic ring for both the methylamino and nitro functional groups. It is found that these groups have a significant effect on the hydration of MNA with water molecules forming two weak hydrogen bonds with both the methylamino and nitro groups. The hydration structures around the functional groups in MNA in water are different from those that have been found for methylamine, nitrobenzene, and benzene in aqueous solutions, and these differences together with weak hydrogen bonds explain the lower solubility of MNA in water. The RDFs together with SDFs provide a tool for the understanding the hydration of MNA (and other molecules) and therefore their solubility.     

Dynamics of benzene guest inside a self-assembled cylindrical capsule: a combined solid-state 2H NMR and molecular dynamics simulation study

The reorientational dynamics of benzene-d(6) molecules hosted into the cavity of a cavitand-based, self-assembled capsule was investigated by Molecular Dynamics (MD) simulations and temperature-dependent solid-state (2)H NMR spectroscopy. MD simulations were preliminarily performed to assess the motional models of the guest molecules inside the capsules. An in-plane fast reorientation of the benzene guest around the C(6) symmetry axis (B1 motion), characterized by correlation times of the order of picoseconds, was predicted with an activation barrier ( approximately 8 kJ/mol) very similar to that found for neat benzene in the liquid state. An out-of-plane reorientation corresponding to a nutation of the C(6) symmetry axis in a cone angle of 39 degrees (B2 motion, 373 K) with an activation barrier ( approximately 39 kJ/mol) definitely larger than that of liquid benzene was also anticipated. In the temperature range 293-373 K correlation times of the order of a nanosecond have been calculated and a transition from fast to slow regime in the (2)H NMR scale has been predicted between 293 and 173 K. (2)H NMR spectroscopic analysis, carried out in the temperature range 173-373 K on the solid capsules containing the perdeuterated guest (two benzene molecules/capsule), confirmed the occurrence of the B1 and B2 motions found in slow exchange in the (2)H NMR time scale. Line shape simulation of the (2)H NMR spectral lines permitted defining a cone angle value of 39 degrees at 373 K and 35 degrees at 173 K for the nutation axis. The T(1) values measured for the (2)H nuclei of the encapsulated aromatic guest gave correlation times and energetic barrier for the in-plane motion B1 in fine agreement with theoretical calculation. The experimental correlation time for B2 as well as the corresponding energetic barrier are in the same range found for B1. A molecular mechanism for the encapsulated guest accounting for the B1 and B2 motions was also provided.     

Detailed insight into the hydrogen bonding interactions in acetone-methanol mixtures. A molecular dynamics simulation and Voronoi polyhedra analysis study

Voronoi polyhedra (VP) analysis of mixtures of acetone and methanol is reported on the basis of molecular dynamics computer simulations, performed at 300 K and 1 bar. The composition of the systems investigated covers the entire range from neat acetone to neat methanol. Distribution of the volume, reciprocal volume and asphericity parameter of the VP as well as that of the area of the individual VP faces and of the radius of the empty voids located between the molecules are calculated. To investigate the tendency of the like molecules to self-associate the analyses are repeated by disregarding one of the two components. The self-aggregates of the disregarded component thus turn into large empty voids, which are easily detectable in VP analysis. The obtained results reveal that both molecules show self-association, but this behavior is considerably stronger among the acetone than among the methanol molecules. The strongest self-association of the acetone and methanol molecules is found in their mole fraction ranges of 02-0.5 and 0.5-0.6, respectively. The caging effect around the methanol molecules is found to be stronger than around acetones. Finally, the local environment of the acetone molecules turns out to be more spherical than that of the methanols, not only in the respective neat liquids, but also in their mixtures.     

Vapor-phase chemical equilibrium for the hydrogenation of benzene to cyclohexane from reaction-ensemble molecular simulation

The reaction-ensemble Monte-Carlo (REMC) molecular simulation method was used to study the vapor-phase chemical equilibrium for the reaction of hydrogenation of benzene to cyclohexane. A one-center modified Buckingham exponential-6 (1CMBE6) effective pair potential model (that had already been used to predict thermodynamic properties and liquid–liquid equilibria of helium+hydrogen mixtures) was used for hydrogen. Six-center modified Buckingham exponential-6 (6CMBE6) effective pair potential models (that had already been used to reproduce the saturated liquid and vapor densities, vapor pressures, second virial coefficients, and critical parameters of the six-membered ring molecules), were used for benzene and cyclohexane. No binary adjustable parameters were needed to compute the unlike-pair Buckingham exponential-6 interactions in the ternary system. Simulation results were obtained for the effect of some operating variables such as temperature (from 500 to 650 K), pressure (from 1 to 30 bar), and hydrogen to benzene feed mole ratio (from 1.5:1 to 6:1) on the reaction conversion, molar composition, and mass density of the ternary system at equilibrium. These results were found to be in excellent agreement with calculations using the predictive Soave–Redlich–Kwong (PSRK) group contribution equation of state.     

Molecular dynamics simulation of liquid mixtures of benzene with chlorobenzene

A molecular dynamics method is used to simulate liquid mixtures of benzene and chlorobenzene at different concentrations. Radial angular distribution functions (RADFs) for distances between the benzene ring planes and the angle between them were calculated to analyze the structure of pure components and mixtures. In chlorobenzene, the highest RADF maximum at a distance between the mass centers of the benzene rings of about 4 Å corresponds to the stacked configurations of molecules, and at 5–7 Å the number of stacked contacts are much less than that at 4 Å and is comparable with the orthogonal ones. In liquid benzene, the number of stacked and orthogonal configurations is approximately equal in a range from 4 Å to 7 Å. RADF for benzene reveals extended regions of correlation, which gives evidence of the occurrence of agglomerates bound by specific interactions between the benzene rings. These agglomerates are not characteristic of chlorobenzene, but the presence of maxima on the radial distribution function for the distances between chlorine atoms indicates chlorine aggregation. The effect of halogen aggregation on the structure of benzene-chlorobenzene mixtures is considered. The obtained results are compared with the data on molecular light scattering.     

The interface between water and a hydrophobic gas

Classical molecular dynamics simulations have been performed to investigate the interface between liquid water and methane gas under methane hydrate forming conditions. The local environments of the water molecules were studied using order parameters which distinguish between liquid water, ice and methane hydrate phases. Bulk water and water/air interfaces were also studied to allow comparisons to be made between water molecules in the different environments and to determine the effects of the different methane densities studied. Good agreement between experimental and calculated surface tensions is obtained if long range corrections are included. The water surface is found to have a structure which is very similar to that of bulk water, but more tetrahedral, and more clathrate-like than ice-like. In these simulations the concentration of methane in water at the interface is shown to be appropriate for clathrates at higher gas densities (pressures). The orientation of water molecules around methane molecules in the interfacial region appears to depend only weakly on pressure and one of the difficulties in forming hydrate is the availability of water molecules tangential to the hydrate cage. At the interface, the water structure is more disordered than in the bulk water region with increased occurrence compared with the bulk of those angles and orientations found in the clathrate structure.     

A molecular dynamics study of nitric oxide in water: diffusion and structure

We present molecular dynamics simulations of the diffusion coefficients and structure of water-nitric oxide mixtures at ambient (298 K) and in vivo (310 K) conditions. A two-site rigid-body molecular model with partial charges and a Lennard-Jones potential on both sites is proposed for nitric oxide and used in conjunction with the extended simple point-charge model for liquid water in our simulations. The diffusion coefficients obtained from the simulations are in good agreement with experimental data. The results from intermolecular partial pair functions show that under these thermodynamic conditions, the existence of nitric oxide in liquid water has little impact on the structure of water and the tendency to form H bonds between water molecules. We also find that it is unlikely that H bonds form between the hydrogen atoms in water and either the nitrogen or the oxygen atom on the nitric oxide at the temperatures and densities examined in this study. This study suggests that in low concentrations nitric oxide molecules exist as free molecules in liquid water rather than forming complexes with water molecules.     

Partial molar volume of n-alcohols at infinite dilution in water calculated by means of scaled particle theory

The partial molar volume of n-alcohols at infinite dilution in water is smaller than the molar volume in the neat liquid phase. It is shown that the formula for the partial molar volume at infinite dilution obtained from the scaled particle theory equation of state for binary hard sphere mixtures is able to reproduce in a satisfactory manner the experimental data over a large temperature range. This finding implies that the packing effects play the fundamental role in determining the partial molar volume at infinite dilution in water also for solutes, such as n-alcohols, forming H bonds with water molecules. Since the packing effects in water are largely related to the small size of its molecules, the latter feature is the ultimate cause of the decrease in partial molar volume associated with the hydrophobic effect.     

The Kirkwood-Buff theory of solutions and the local composition of liquid mixtures

The present paper is devoted to the local composition of liquid mixtures calculated in the framework of the Kirkwood-Buff theory of solutions. A new method is suggested to calculate the excess (or deficit) number of various molecules around a selected (central) molecule in binary and multicomponent liquid mixtures in terms of measurable macroscopic thermodynamic quantities, such as the derivatives of the chemical potentials with respect to concentrations, the isothermal compressibility, and the partial molar volumes. This method accounts for an inaccessible volume due to the presence of a central molecule and is applied to binary and ternary mixtures. For the ideal binary mixture it is shown that because of the difference in the volumes of the pure components there is an excess (or deficit) number of different molecules around a central molecule. The excess (or deficit) becomes zero when the components of the ideal binary mixture have the same volume. The new method is also applied to methanol + water and 2-propanol + water mixtures. In the case of the 2-propanol + water mixture, the new method, in contrast to the other ones, indicates that clusters dominated by 2-propanol disappear at high alcohol mole fractions, in agreement with experimental observations. Finally, it is shown that the application of the new procedure to the ternary mixture water/protein/cosolvent at infinite dilution of the protein led to almost the same results as the methods involving a reference state.