Structure and dynamics of nanoconfined water and aqueous solutions.

This review is devoted to discussing recent progress on the structure, thermodynamic, reactivity, and dynamics of water and aqueous systems confined within different types of nanopores, synthetic and biological. Currently, this is a branch of water science that has attracted enormous attention of researchers from different fields interested to extend the understanding of the anomalous properties of bulk water to the nanoscopic domain. From a fundamental perspective, the interactions of water and solutes with a confining surface dramatically modify the liquid's structure and, consequently, both its thermodynamical and dynamical behaviors, breaking the validity of the classical thermodynamic and phenomenological description of the transport properties of aqueous systems. Additionally, man-made nanopores and porous materials have emerged as promising solutions to challenging problems such as water purification, biosensing, nanofluidic logic and gating, and energy storage and conversion, while aquaporin, ion channels, and nuclear pore complex nanopores regulate many biological functions such as the conduction of water, the generation of action potentials, and the storage of genetic material. In this work, the more recent experimental and molecular simulations advances in this exciting and rapidly evolving field will be reported and critically discussed.

Influence of Cholesterol on the Dynamics of Hydration in Phospholipid Bilayers.

We investigate the dynamics of interfacial waters in dipalmitoylphosphatidylcholine (DPPC) bilayers upon the addition of cholesterol, by molecular dynamics simulations. Our data reveal that the inclusion of cholesterol modifies the membrane aqueous interfacial dynamics: waters diffuse faster, their rotational decay time is shorter, and the DPPC/water hydrogen bond dynamics relaxes faster than in the pure DPPC membrane. The observed acceleration of the translational water dynamics agrees with recent experimental results, in which, by means of NMR techniques, an increment of the surface water diffusivity is measured upon the addition of cholesterol. A microscopic analysis of the lipid/water hydrogen bond network at the interfacial region suggests that the mechanism underlying the observed water mobility enhancement is given by the rupture of a fraction of interlipid water bridge hydrogen bonds connecting two different DPPC molecules, concomitant to the formation of new lipid/solvent bonds, whose dynamics is faster than that of the former. The consideration of a simple two-state model for the decay of the hydrogen bond correlation function yielded excellent results, obtaining two well-separated characteristic time scales: a slow one (∼250 ps) associated with bonds linking two DPPC molecules, and a fast one (∼15 ps), related to DPPC/solvent bonds.

Molecular Dynamics Simulations of Ibuprofen Release from pH-Gated Silica Nanochannels.

The iboprufen delivery process from cylindrical silica pores of 3 nm diameter, with polyamine chains anchored at the pore outlets, was investigated by means of massive molecular dynamics simulations. Effects from pH were introduced by considering polyamine chains with different degrees of protonation. High, low, and intermediate pH environments were investigated. The increment of the acidity of the environment leads to a significant decrease of the pore aperture, yielding an effective diameter, for the lowest pH case, that is 3.5 times smaller than the one associated with the highest pH one. Using a biased sampling procedure, Gibbs free energy profiles for the ibuprofen delivery process were obtained. The joint analysis of the corresponding profiles, time evolution of the ibuprofen position within the channel, orientation of the molecule, and instantaneous effective diameter of the gate suggest a three-step mechanism for ibuprofen delivery. A complementary analysis of the translational mobility of ibuprofen along the axial direction of the channel revealed a subdiffusive dynamics in the low and intermediate pH cases. Deviations from Brownian diffusive dynamics are discussed and compared with direct experimental results.

Equilibrium and Dynamical Characteristics of Imidazole Langmuir Monolayers on Graphite Sheets.

Using molecular dynamics techniques, we examine structural and dynamical characteristics of liquid-like imidazole (Im) monolayers physisorbed onto a planar graphite sheet, at T = 384 K. Our simulations reveal that molecular orientations in the saturated monolayer exhibit a bistable distribution, characterized by an inner parallel arrangement of the molecules in close contact with the substrate and a slanted alignment, in those lying in adjacent, outer locations. Compared to the results found in three-dimensional, bulk phases, the analysis of the spatial correlations between sites participating in hydrogen bonding shows a clear enhancement of the intermolecular interactions, which also leads to stronger dipolar correlations. As a result, the gross structural features of the monolayer can be cast in terms of mesoscopic domains, comprising units articulated via winding hydrogen bonds, that persist along typical time intervals of a few tens of picoseconds. On the dynamical side, a similar comparison of the characteristic decorrelation time for orientational motions shows a 4-fold increment. Contrasting, the reduction of the system dimensionality leads to a larger diffusion constant. Possible substrate-induced anisotropies in the diffusive motions are also investigated.

Solvation of Coumarin 480 within nano-confining environments: structure and dynamics.

Equilibrium and dynamical characteristics pertaining to the solvation of the fluorescent probe Coumarin 480 within different confining environments are investigated using molecular dynamics simulations. Three kinds of confining systems are examined: (i) the cetyltrimethylammonium bromide (CTAB)/isooctane/1-hexanol/water; cationic inverse micelle (IM) (ii) a CTAB/water direct micelle (DM), and (iii) a silica-surfactant nanocomposite, comprising a cylindrical silica pore (SP) containing small amounts of water and CTAB species adsorbed at the pore walls. The solvation structures in the three environments differ at a qualitative level: an exchange between bulk- and interface-like solvation states was found in the IM, whereas in the DM, the solvation states of the probe are characterized by its embedding at the interface, trapped among the surfactant heads and tails. Within the SP structure, the coumarin exhibits alternations between internal and interfacial solvation states that occur on a ∼20 ns time scale and operate via 90° rotations of its molecular plane. The solvation responses of the environment following a vertical excitation of the probe are also investigated. Solvation times resulted between 2 and 1000 times longer than those found in bulk water, with a fast-to-slow trend IM→DM→SP, which can be interpreted in terms of the solvation structures that prevail in each case.

Structure and dynamics of nonionic surfactants adsorbed at vacuum/ionic liquid interfaces.

Structural and dynamical properties related to the adsorption of nonionic surfactants at vacuum/ionic liquid interfaces were studied using molecular dynamics simulations. Specifically, the surface activity of pentaethylene glycol monododecyl ether (C12E5) was investigated at the free interface of an imidazolium-based room temperature ionic liquid (RTIL), at different surface densities. At low surface coverages, the incorporation of C12E5 does not produce meaningful changes in the vacuum/RTIL interface: the C12E5 hydrophobic tails remain entangled with those of the RTIL cation groups in the outer shell, whereas the C12E5 hydrophilic heads reside at an inner layer. At high surface coverages, the structure in the substrate-in terms of the features exhibited by the local density profiles-practically vanishes; the interface becomes wider and the surfactant molecules shift toward more external positions. Information about the local structure of the interface at high surface densities can be recovered by performing a tessellation procedure. For the sake of comparison, the surface behavior of two commonly used ionic surfactants, sodium dodecyl sulfate and dodecyl trimethyl ammonium chloride, were also studied. The modifications in the width and structure of the bare vacuum/RTIL interface due to the presence of the ionic surfactants are markedly milder than those observed for the nonionic surfactant. Moreover, the RTIL seemed to behave as a better solvent for the chloride counterions than for sodium ones; which were found to remain bound to the surfactant head groups. An analysis of the dynamics at the surface was also performed. Our results indicate that the presence of increasing amounts of nonionic surfactants leads to a gradual reduction of the mobility of the RTIL species. When ionic surfactants are adsorbed, these retardations are even more severe for the surfactant head groups, where the corresponding diffusion coefficients show reductions of practically 1 order of magnitude.

Liquid methanol confined within functionalized silica nanopores. 2. Solvation dynamics of coumarin 153.

Equilibrium and dynamical characteristics pertaining to the solvation of the fluorescent probe coumarin 153 in liquid methanol confined within cylindrical silica pores are investigated using molecular dynamics techniques. Three kinds of pores are examined: (i) Soft hydrophobic cavities, in which wall-solvent interactions were exclusively of the Lennard-Jones type; (ii) Hydrophilic cavities, in which unsaturated oxygen sites at the wall were transformed into hydroxyl groups; (iii) Rugged pores, in which 60% of the polar groups were transformed into bulkier and mobile trimethylsilyl moieties. Equilibrium solvation structures in the three pores differ considerably: In hydrophobic environments, the solute remains adsorbed to the pore wall, with its molecular plane mostly parallel to the interface. Upon hydroxylation, the solid interface becomes preferentially coated by methanol, leading to a bistable solvation state of the probe, with alternation of "wall-like" and "bulk-like" events. An increment in the interface roughness promotes a solvation structure characterized by the embedding of the probe within a wall domain surrounded by trimethylsilyl groups. In hydrophobic environments, the relevant dynamical modes of the probe can be cast in terms of in-the-wall rotations, whereas in hydrophilic pores, out-of-the-wall evolutions are also present. The embedding of the probe at wall domains in more rugged pores, leads to restrained angular motions, with maximum amplitudes of the order of 20°. Results of early stages of the solvation response of the environment following a vertical excitation of the probe are also presented. During the initial 30 ps, we found no evidence of modifications in the spatial localizations of the probe. The overall responses are found to be between 2 and 4.5 times slower than the one observed in the bulk, being the fastest relaxation the one associated to rugged pores whereas the slowest one corresponds to hydrophilic cavities. These features are rationalized in terms of the composition of the first solvation shells and the local dynamical inhomogeneities prevailing within the different regions of the pores.

Structure and dynamics of liquid methanol confined within functionalized silica nanopores.

Molecular dynamics simulations have been carried out to investigate the structure and dynamics of liquid methanol confined in 3.3 nm diameter cylindrical silica pores. Three cavities differing in the characteristics of the functional groups at their walls have been examined: (i) smooth hydrophobic pores in which dispersive forces prevail, (ii) hydrophilic cavities with surfaces covered by polar silanol groups, and (iii) a much more rugged pore in which 60% of the previous interfacial hydroxyl groups were replaced by the bulkier trimethylsilyl ones. Confinement promotes a considerable structure at the vicinity of the pore walls which is enhanced in the case of hydroxylated surfaces. Moreover, in the presence of the trimethylsilyl groups, the propagation of this interface-induced spatial ordering extends down to the central region of the pore. Concerning the dynamical modes, we observed an overall slowdown in both the translational and rotational motions. An analysis of these mobilities from a local perspective shows that the largest retardations operate at the vicinity of the interfaces. The gross features of the rotational dynamics were analyzed in terms of contributions arising from bulk and surface states. Compared to the bulk dynamical behavior, the characteristic timescales associated with the rotational motions show the most dramatic increments. A dynamical analysis of hydrogen bond formation and breaking processes is also included.

Confined polar mixtures within cylindrical nanocavities.

Using molecular dynamics experiments, we have extended our previous analysis of equimolar mixtures of water and acetonitrile confined between silica walls [J. Phys. Chem. B 2009, 113, 12744] to examine similar solutions trapped within carbon nanotubes and cylindrical silica pores. Two different carbon tube sizes were investigated, (8,8) tubes, with radius R(cnt) = 0.55 nm, and (16,16) ones, with R(cnt) = 1.1 nm. In the narrowest tubes, we found that the cylindrical cavity is filled exclusively by acetonitrile; as the radius of the tube reaches approximately 1 nm, water begins to get incorporated within the inner cavities. In (16,16) tubes, the analysis of global and local concentration fluctuations shows a net increment of the global acetonitrile concentration; in addition, the aprotic solvent is also the prevailing species at the vicinity of the tube walls. Mixtures confined within silica nanopores of radius approximately 1.5 nm were also investigated. Three pores, differing in the effective wall/solvent interactions, were analyzed, (i) a first class, in which dispersive forces prevail (hydrophobic cavities), (ii) a second type, where oxygen sites at the pore walls are transformed into polar silanol groups (hydrophilic cavities), and (iii) finally, an intermediate scenario, in which 60% of the OH groups are replaced by mobile trimethylsilyl groups. Within the different pores, we found clear distinctions between the solvent layers that lie in close contact with the silica substrate and those with more central locations. Dynamical modes of the confined liquid phases were investigated in terms of diffusive and rotational time correlation functions. Compared to bulk results, the characteristic time scales describing different solvent motions exhibit significant increments. In carbon nanotubes, the most prominent modifications operate in the narrower tubes, where translations and rotations become severely hindered. In silica nanopores, the manifestations of the overall retardations are more dramatic for solvent species lying at the vicinity of trimethylsilyl groups.

Coaxial cross-diffusion through carbon nanotubes.

We present results from nonequilibrium molecular dynamics experiments describing the relaxation of local concentrations at two reservoirs, initially filled with water (W) and acetonitrile (ACN), as they become connected through a membrane composed of (16,16) carbon nanotubes. Within the hydrophobic nanotube cavities, the equilibrium concentrations contrast sharply to those observed at the reservoirs, with a clear enhancement of ACN, in detriment of W. From the dynamical side, the relaxation involves three well-differentiated stages; the first one corresponds to the equilibration of individual concentrations within the nanotubes. An intermediate interval with Fickian characteristics follows, during which the overall transport can be cast in terms of coaxial opposite fluxes, with a central water domain segregated from an external ACN shell, in close contact with the tube walls. We also found evidence of a third, much slower, mechanism to reach equilibration, which involves structural modifications of tightly bound solvation shells, in close contact with the nanotube rims.

Anomalous dynamics of hydration water in carbohydrate solutions.

Molecular dynamics simulations have been carried out to investigate the dynamics of fructose aqueous solutions up to 70 wt % concentration. We find that the hydrogen (H)-bonded network of fructose molecules extends with increasing sugar content and forms a structurally heterogeneous system around and above 45 wt % concentration, characterized as a percolated-like solute domain permeated by patchy regions of solvent. The presence of such aggregates in concentrated solutions promotes the slowing down of water translational, reorientational, and H-bonding dynamics, typical of many biomolecular environments. Analysis of the effects of the topological and energetic disorder of the sugar aggregates on vicinal water dynamics, similar to that recently carried out for the hydration layer of proteins by Pizzitutti et al. (J. Phys. Chem. B 2007, 111, 7584), reveals many similarities between the dynamical anomaly of the hydration layers of both systems. Like a protein surface, topological and energetic disorders of the sugar aggregates both contribute to the translational diffusion anomaly. However, unlike in the vicinity of a protein surface, the rotational relaxation is also hindered by the topological disorder created by the intertwined, percolating sugar clusters in concentrated solutions.

Polar mixtures under nanoconfinement.

We present results from molecular dynamics simulations describing structural and dynamical characteristics of equimolar mixtures of water and acetonitrile, confined between two silica walls separated at interplate distances of d=0.6, 1, and 1.5 nm. Two different environments were investigated: a first one where wall-solvent dispersion forces prevail (hydrophobic confinement) and a second one in which the terminal O atoms at the silica surface are transformed into silanol groups (hydrophilic confinement). For the former case, we found that, at the shortest interplate distance examined, the confined region is devoid of water molecules. At an interplate distance of the order of 1 nm, water moves into the confined region, although, in all cases, there is a clear enhancement of the local concentration of acetonitrile in detriment of that of water. Within hydrophilic environments, we found clear distinctions between a layer of bound water lying in close contact with the silica substrates and a minority of confined water that occupies the inner liquid slab. The bound aqueous layer is fully coordinated to the silanol groups and exhibits minimal hydrogen bonding with the second solvation layer, which exclusively includes acetonitrile molecules. Dynamical characteristics of the solvent mixture are analyzed in terms of diffusive and rotational motions in both environments. Compared to bulk mixtures, we found significant retardations in all dynamical modes, with those ascribed to water molecules bound to the hydrophilic plates being the most dramatic.

Encapsulation of small ionic molecules within alpha-cyclodextrins.

Results from molecular dynamics experiments pertaining to the encapsulation of ClO4- within the hydrophobic cavity of an aqueous alpha-cyclodextrin (alpha-CD) are presented. Using a biased sampling procedure, we constructed the Gibbs free energy profile associated with the complexation process. The profile presents a global minimum at the vicinity of the primary hydroxyl groups, where the ion remains tightly coordinated to four water molecules via hydrogen bonds. Our estimate for the global free energy of encapsulation yields DeltaGenc approximately -2.5 kBT. The decomposition of the average forces acting on the trapped ion reveals that the encapsulation is controlled by Coulomb interactions between the ion and OH groups in the CD, with a much smaller contribution from the solvent molecules. Changes in the previous results, arising from the partial methylation of the host CD and modifications in the charge distribution of the guest molecule are also discussed. The global picture that emerges from our results suggests that the stability of the ClO4- encapsulation involves not only the individual ion but also its first solvation shell.

Molecular dynamics simulations of AOT-water/formamide reverse micelles: structural and dynamical properties.

We present results from molecular dynamics simulations performed on reverse micelles immersed in cyclohexane. Three different inner polar phases are considered: water (W), formamide (FM), and an equimolar mixture of the two solvents. In all cases, the surfactant was sodium bis(2-ethylhexyl) sulfosuccinate (usually known as AOT). The initial radii of the micelles were R approximately 15 A, while the corresponding polar solvent-to-surfactant molar ratios were intermediate between w(0)=4.3 for FM and w(0)=7 for W. The resulting overall shapes of the micelles resemble distorted ellipsoids, with average eccentricities of the order of approximately 0.75. Moreover, the pattern of the surfactant layer separating the inner pool from the non-polar phase looks highly irregular, with a roughness characterized by length scales comparable to the micelle radii. Solvent dipole orientation polarization along radial directions exhibit steady growths as one moves from central positions toward head group locations. Local density correlations within the micelles indicate preferential solvation of sodium ionic species by water, in contrast to the behavior found in bulk equimolar mixtures. Still, a sizable fraction of approximately 90% of Na(+) remains associated with the head groups. Compared to bulk results, the translational and rotational modes of the confined solvents exhibit important retardations, most notably those operated in rotational motions where the characteristic time scales may be up to 50 times larger. Modifications of the intramolecular connectivity expressed in terms of the average number of hydrogen bonds and their lifetimes are also discussed.

Intermolecular polarizability dynamics of aqueous formamide liquid mixtures studied by molecular dynamics simulations.

A molecular dynamics simulation study is presented for the relaxation of the polarizability anisotropy in liquid mixtures of formamide and water, using a dipolar induction scheme that involves the intrinsic polarizability and first hyperpolarizability tensors of the molecules, and the dipole-quadrupole polarizability of water species. The long time diffusive decay of the collective polarizability anisotropy correlations exhibits a substantial slowing down as the formamide mole fraction increases in the mixture. The diffusive times for the polarizability relaxation obtained from the authors' simulations are in good agreement with optical Kerr effect experimental data, and they are found to correlate nearly linearly with the estimated mean lifetimes of the hydrogen bonds within the mixture, suggesting that the relaxation of the hydrogen bond network is responsible to some extent for the collective relaxation of the polarizability anisotropy of the mixture. The short time behavior of the polarizability anisotropy relaxation was investigated by computing the nuclear response function, R(t), which is very rapidly dominated by the formamide contribution as it is added to water, due to the much larger polarizability anisotropy of formamide molecules compared to that of water. Several contributions to the Raman spectrum were also analyzed as a function of composition, and the dynamical origin of the different bands was determined.

Computational study of structural and dynamical properties of formamide-water mixtures.

A molecular dynamics simulation study of structural and dynamical properties in liquid mixtures of formamide and water is presented. Site-site radial pair distribution functions, local mole fractions, pair energy distributions, and tetrahedral orientational order are the quantities analyzed to investigate the local structure in the simulated mixtures, along with a review of the intermolecular structure in terms of the distribution of hydrogen bonds. Our results indicate that there is a substitution of formamide molecules by water in the hydrogen bonds and a formation of a common hydrogen bond network. By analyzing the extent of tetrahedral order in the liquid as a function of composition, it is observed that whereas the tetrahedral network of liquid water is progressively lost by increasing the formamide concentration, the water structure within the first coordination shell is preserved and somewhat enhanced. The hydrogen-bond mean lifetimes were estimated by performing a time integration of the autocorrelation functions of bond occupation numbers. The lifetimes associated with hydrogen bonds between water, formamide, and interspecies pairs are found to increase with increasing formamide concentration. The lifetimes of the water hydrogen bonds show the largest variations, supporting the picture of an enhancement of the water structure among the nearest neighbors within the first coordination shell. We have used two different force field models for water, SPC/E [J. C. Berendsen et al., J. Phys. Chem. 91, 6269 (1987)] and TIP4P/2005 [J. L. F. Abascal and C. Vega, J. Chem. Phys. 123, 234505 (2005)]. Our results for structural and dynamical properties yield very small differences between those models, the TIP4P/2005 predicting a slightly more structured liquid and, consequently, exhibiting a slightly slower translational and librational dynamics.

Molecular dynamics study of polarizability anisotropy relaxation in aromatic liquids and its connection with local structure.

The collective polarizability anisotropy dynamics in a set of three aromatic liquids, benzene (Bz), hexafluorobenzene (HFB), and 1,3,5-trifluorobenzene (TFB), has been studied by molecular dynamics simulation. These liquids have very similar shapes, but different electrostatic interactions due to opposite polarities of C-H and C-F bonds, giving rise to different local intermolecular structures in the liquid phase. We have investigated how these structural arrangements affect polarizability anisotropy dynamics observed in optical Kerr-effect (OKE) spectroscopy. We have modeled the interaction-induced polarizability with the first-order dipole-induced dipole approximation, with the molecular polarizability distributed over the carbon sites. Local contributions to the librational OKE spectrum were computed separately for molecules participating in parallel or perpendicular relative orientations within the first coordination shell. We found that the relative locations of parallel and perpendicular librational bands of the OKE spectra are closely related to the corresponding pair energy distributions of the closest four neighbors of a given molecule, corresponding to a model of a harmonic oscillator in a cage of nearest neighbors. This model predicts higher librational frequencies for more attractive intermolecular interactions, which in all three liquids correspond to parallel local arrangements. On the diffusive orientational time scale, all three liquids exhibit slower relaxation of molecules in parallel arrangements, although the difference in relaxation rates is substantial only in TFB, which has the strongest tendency toward parallel stacking. The analysis of the collective polarizability relaxation was performed using two different approaches, the projection scheme (J. Chem. Phys. 1980, 72, 2801) and the theory developed by Steele (Mol. Phys. 1987, 61, 1031) for the second time derivatives applied to collective time correlations. Both approaches allow the decomposition of the OKE response into contributions from orientational relaxation and other dynamical processes. We find that they lead to different predictions on how the response depends on collective reorientation and processes arising from fluctuations in the interaction-induced polarizability. We discuss the reasons for these differences and the advantages and disadvantages of the two analysis schemes.

Effects of molecular association on polarizability relaxation in liquid mixtures of benzene and hexafluorobenzene.

In this work we have studied the relaxation dynamics of the many-body polarizability anisotropy in liquid mixtures of benzene (Bz) and hexafluorobenzene (Hf) at room temperature by femtosecond optical heterodyne-detected Raman-induced Kerr effect spectroscopy (OHD-RIKES) experiments and molecular dynamics (MD) simulations. The computed polarizability response arising from intermolecular interactions was included using the first-order dipole-induced-dipole model with the molecular polarizability distributed over the carbon sites of each molecule. We found good qualitative agreement between experiments and simulations in the features exhibited by the nuclear response function R(t) for pure liquids and mixtures. The long-time diffusive decay of R(t) was observed to vary substantially with composition, slowing down noticeably with dilution of each of the species as compared with that in the corresponding pure liquids. MD simulation shows that the effect on R(t) is due to the formation of strong and localized intermolecular association between Bz and Hf species that hinder the rotational diffusive dynamics. The formation of these Bz-Hf complexes in the liquid mixtures also modifies the rotational diffusive dynamics of the component species in such a way that cannot be explained solely in terms of a viscosity effect. Even though the computed orientational diffusive relaxation times associated with Bz and Hf are larger by a factor of approximately 2 than those from experiments, we found similar trends in experiments and simulations for these characteristic times as a function of composition. Namely, the collective and single-molecule orientational correlation times associated with Bz are observed to grow monotonically with the dilution of Bz, while those corresponding to Hf species exhibit a maximum at the equimolar composition. We attribute the quantitative discrepancy between experiments and simulations to the use of the Williams potential, which seems to overestimate the intermolecular interactions and thus predicts not only a slower translational dynamics but also a slower rotational diffusion dynamics than in real fluids.

Polarizability response in polar solvents: molecular-dynamics simulations of acetonitrile and chloroform.

The relaxation of the many-body polarizability in liquid acetonitrile and chloroform at room temperature was studied by molecular-dynamics simulations. The collective polarizability induced by intermolecular interactions was included using first- and all-orders dipole-induced-dipole models and calculated considering both molecule-centered and distributed site polarizabilities. The anisotropic response was analyzed using a separation scheme that allows a decomposition of the total response in terms of orientational and collision-induced effects. We found the method effective in approximately separating the contributions of these relaxation mechanisms, although the orientational-collision-induced interference makes a non-negligible contribution to the total response. In both liquids the main contribution to the anisotropic response is due to orientational dynamics, but intermolecular collision-induced (or translational) effects are important, especially at short times. We found that higher-order interaction-induced effects were essentially negligible for both liquids. Larger differences were found between the center-center and site-site models, with the latter showing faster polarizability relaxation and better agreement with experiment. Isotropic and anisotropic spectra were computed from the corresponding time correlation functions. The lowest-frequency contributions are largely suppressed in the isotropic spectra and their overall shape is similar to the purely collision-induced contribution to the anisotropic spectra, but with an amplitude which is smaller by a factor of approximately 5 in acetonitrile and approximately 3 in chloroform.

Investigation of benzene-hexafluorobenzene dynamics in liquid binary mixtures.

The structure and microscopic dynamics of liquid mixtures of benzene and hexafluorobenzene at room temperature and several compositions have been studied by molecular-dynamics simulations. In this implementation we have rescaled the intermolecular H-F cross potential parameters obtained from the Lorentz-Berthelot combining rules, in order to avoid the substantial overestimation of the energy of mixing predicted by the model when the usual rules are employed. We found that a reduction in the strength of cross H-F interactions by 50% relative to the geometric mean is required in order to get a good agreement with experiments. Radial-angular pair-correlation functions between like and unlike species have been computed and analyzed, by comparing them with the correlations in the corresponding neat liquids. We have also studied the microscopic intermolecular momentum transfer, by computing the time correlation function between the initial velocity of a central molecule and later velocities of neighboring molecules. Structural and dynamical information extracted from the mentioned functions seem to be consistent with the picture of relatively long-lived benzene-hexafluorobenzene (Bz-Hf) complexes present in the mixtures, which would be responsible for the considerable perturbation of the structure in the first shell of like species, and would be moving within the liquid in a parallel face-to-face configuration. Using the tools developed originally to estimate hydrogen-bond lifetimes in liquids, we have computed the lifetimes of the Bz-Hf complexes as a function of the mixture composition, by two different methods: the direct time-averaging scheme and from the autocorrelation function of bond occupation numbers. The obtained lifetimes are strongly dependent on the scheme chosen to compute the characteristic times. We have obtained for the Bz-Hf dimer in solution, at room temperature, lifetimes in the range of 30-40 ps from averaging schemes and around 60-120 ps from autocorrelation function methods. In the latter case, the longest times correspond to the equimolar mixture.