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.     

Molecular dynamics study of the temperature-dependent Optical Kerr effect spectra and intermolecular dynamics of room temperature ionic liquid 1-methoxyethylpyridinium dicyanoamide

We have performed classical molecular dynamics simulations to calculate the Optical Kerr effect (OKE) spectra of 1-methoxyethylpyridinium dicyanoamide, a room-temperature ionic liquid (IL) which has been recently studied by Shirota and Castner (Shirota, H. ; Castner, E. J. Phys. Chem. A 2005, 109, 9388-9392) in comparison to its neutral isoelectronic solvent mixture. Our theoretical and computational studies show that the decay of the collective polarizability anisotropy correlation exhibits several different time scales originating from inter- and intramolecular dynamics, in good agreement with experiments. What's more, we find that the portion of the collective anisotropic polarizability relaxation due to "interaction-induced" phenomena is important at times much longer than those observed in normal solvents when these are far from their glass transition temperature. From our long (60 ns) molecular dynamics simulations, we are able to determine the appropriate time scales for orientational relaxation and interaction-induced processes occurring in the liquid. We find that the cationic contribution to the OKE signal is predominant. Because of the slow nature of relaxation processes in ILs, these calculations are very time, memory, and storage intensive. In the context of this research, we have developed a polarizable force field for this system and also theoretical methodology to generate molecular polarizabilities for arbitrarily shaped molecules and ions from corresponding atomic polarizabilities. We expect this methodology to have an important impact on the speed of molecular dynamics simulations of polarizable systems in the future.     

Hydrogen bonded structure and dynamics of liquid-vapor interface of water-ammonia mixture: an ab initio molecular dynamics study

We have carried out ab initio molecular dynamics simulations of a liquid-vapor interfacial system consisting of a mixture of water and ammonia molecules. We have made a detailed analysis of the structural and dynamical properties of the bulk and interfacial regions of the mixture. Among structural properties, we have looked at the inhomogeneous density profiles of water and ammonia molecules, hydrogen bond distributions, orientational profiles, and also vibrational frequency distributions of bulk and interfacial molecules. It is found that the interfacial molecules show preference for specific orientations so as to form water-ammonia hydrogen bonds at the interface with ammonia as the acceptor. The structure of the system is also investigated in terms of inter-atomic voids present in the system. Among the dynamical properties, we have calculated the diffusion, orientational relaxation, hydrogen bond dynamics, and vibrational spectral diffusion in bulk and interfacial regions. It is found that the diffusion and orientation relaxation of the interfacial molecules are faster than those of the bulk. However, the hydrogen bond lifetimes are longer at the interface which can be correlated with the time scales found from the decay of frequency time correlations.     

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.     

Time-resolved optical Kerr-effect investigation on CS2/polystyrene mixtures

The relaxation dynamics of carbon disulfide are investigated in mixtures with polystyrene (PS) using the time-resolved optical heterodyne-detected optical Kerr effect (OHD-OKE). The data are analyzed using both the model-dependent approach, which assumes four distinct temporal responses, and the model-independent Fourier transform approach, which generates a spectral response that can be compared with results obtained by depolarized Rayleigh scattering. A slow dynamics is observed for the OHD-OKE transient decaying exponentially with a time constant that varies from 1.68 ps for neat CS2 to 3.76 ps for the most concentrated CS2PS mixture. The increase of this time constant accompanies an increase in the viscosity of the mixture, so we can associate this component with the diffusive reorientation process of the induced polarizability anisotropy of the carbon disulfide in the mixture. The short-time nuclear response is characterized in the frequency domain by a broad band that peaks around 30 cm(-1) for neat carbon disulfide, and is associated with a complex relaxation pattern. The vibrational distribution shifts to higher frequencies when the PS concentration is increased in the mixture. This result is discussed in terms of an increase in the interaction strength between the PS phenyl rings and the carbon disulfide molecules.     

Local hydrogen bonding dynamics and collective reorganization in water: ultrafast infrared spectroscopy of HOD/D(2)O

We present an investigation into hydrogen bonding dynamics and kinetics in water using femtosecond infrared spectroscopy of the OH stretching vibration of HOD in D(2)O. Infrared vibrational echo peak shift and polarization-selective pump-probe experiments were performed with mid-IR pulses short enough to capture all relevant dynamical processes. The experiments are self-consistently analyzed with a nonlinear response function expressed in terms of three dynamical parameters for the OH stretching vibration: the frequency correlation function, the lifetime, and the second Legendre polynomial dipole reorientation correlation function. It also accounts for vibrational-relaxation-induced excitation of intermolecular motion that appears as heating. The long time, picosecond behavior is consistent with previous work, but new dynamics are revealed on the sub-200 fs time scale. The frequency correlation function is characterized by a 50 fs decay and 180 fs beat associated with underdamped intermolecular vibrations of hydrogen bonding partners prior to 1.4 ps exponential relaxation. The reorientational correlation function observes a 50 fs librational decay prior to 3 ps diffusive reorientation. Both of these correlation functions compare favorably with the predictions from classical molecular dynamics simulations. The time-dependent behavior can be separated into short and long time scales by the 340 fs correlation time for OH frequency shifts. The fast time scales arise from dynamics that are mainly local: fluctuations in hydrogen bond distances and angles within relatively fixed intermolecular configurations. On time scales longer than the correlation time, dephasing and reorientations reflect collective reorganization of the liquid structure. Since the OH transition frequency and dipole are only weakly sensitive to these collective coordinates, this is a kinetic regime which gives an effective rate for exchange of intermolecular structures.     

Hydrogen bond dynamics and water structure in glucose-water solutions by depolarized Rayleigh scattering and low-frequency Raman spectroscopy

The effect of glucose on the relaxation process of water at picosecond time scales has been investigated by depolarized Rayleigh scattering (DRS) experiments. The process is assigned to the fast hydrogen bonding dynamics of the water network. In DRS spectra this contribution can be safely separated from the slower relaxation process due to the sugar. The detected relaxation time is studied at different glucose concentrations and modeled considering bulk and hydrating water contributions. As a result, it is found that in diluted conditions the hydrogen bond lifetime of proximal water molecules becomes about three times slower than that of the bulk. The effect of the sugar on the hydrogen bond water structure is investigated by analyzing the low-frequency Raman (LFR) spectrum sensitive to intermolecular modes. The addition of glucose strongly reduces the intensity of the band at 170 cm(-1) assigned to a collective stretching mode of water molecules arranged in cooperative tetrahedral domains. These findings indicate that proximal water molecules partially lose the tetrahedral ordering typical of the bulk leading to the formation of high density environments around the sugar. Thus the glucose imposes a new local order among water molecules localized in its hydration shell in which the hydrogen bond breaking dynamics is sensitively retarded. This work provides new experimental evidences that support recent molecular dynamics simulation and thermodynamics results.     

Molecular dynamics of liquid acetone determined by depolarized Rayleigh and low-frequency Raman scattering spectroscopy

Slow to ultrafast dynamics of liquid acetone at variable temperature was investigated by depolarized Rayleigh and low-frequency Raman scattering spectroscopy, in the region 0-200 cm(-1). A detailed analysis was performed on the spectra and corresponding time responses, and a consistent view of the molecular dynamics of this dipolar solvent was obtained. The effects of temperature on the spectra were interpreted, and distinct dynamical processes identified. At very low frequencies, or long time scales, acetone dynamics is characterized by a slow diffusive reorientation obeying the Stokes-Einstein-Debye hydrodynamic theory only in the limit of subslip boundary conditions. An alternative model based on the microviscosity concept proved to be able to reproduce this correlation time and its temperature dependence. A comparative analysis of collective and single-molecule reorientational times, these latter estimated from intramolecular Raman spectra, led to an orientational correlation parameter g(2) of unity, which denotes a statistical disorder of molecular polarizability tensors. A fast local restructuring process is putatively responsible for an additional contribution at subpicosecond time scales often referred to as intermediate response in other molecular liquids. The high frequency portion of the dynamical susceptibility showed the signature of librational intermolecular motions, giving rise to an ultrafast decay of the time correlation function of polarizability anisotropy. The overall approach, which provided valuable information on dynamics, structure and molecular interactions of neat acetone, will be applied to acetone electrolytic solutions.     

Ultrafast anisotropy dynamics of water molecules dissolved in acetonitrile

Infrared pump-probe experiments are performed on isolated H(2)O molecules diluted in acetonitrile in the spectral region of the OH stretching vibration. The large separation between water molecules excludes intermolecular interactions, while acetonitrile as a solvent provides substantial hydrogen bonding. Intramolecular coupling between symmetric and asymmetric modes results in the anisotropy decay to the frequency-dependent values of approximately 0-0.2 with a 0.2 ps time constant. The experimental data are consistent with a theoretical model that includes intramolecular coupling, anharmonicity, and environmental fluctuations. Our results demonstrate that intramolecular processes are essential for the H(2)O stretching mode relaxation and therefore can compete with the intermolecular energy transfer in bulk water.     

Excited state hydrogen bond dynamics: coumarin 102 in acetonitrile-water binary mixtures

The time dependent change in the intermolecular response of solvent molecules following photoexcitation of Coumarin 102 (C102) has been measured in acetonitrile-water binary mixtures. Experiments were performed on mixtures of composition x(CH3CN) = 0.25, 0.50, 0.75, and 1.00. At low water concentrations (x(H2O) < or = 0.25) the solvent response is consistent with previous measurements probing dipolar solvation. With increasing water concentration (x(H2O) > or = 0.50) an additional response is found subsequent to dipolar solvation, exhibited as a rapid gain in the solvent's polarizability on a approximately 250 fs time scale. Monte Carlo simulations of the C102:binary mixture system were performed to quantify the number of hydrogen-bonding interactions between C102 and water. These simulations indicate that the probability of the C102 solute being hydrogen bound with two water molecules, both as donors at the carbonyl site, increases in a correlated fashion with the amplitude of the additional response in the measurements. We conclude that excitation of C102 simultaneously weakens and strengthens hydrogen bonding in complexes with two inequivalently bound waters.     

Rotational reorientation dynamics of oxazine 750 in polar solvents

The rotational reorientation dynamics of oxazine 750 (OX750) in the first (with pump pulse at 660 nm) and a higher excited state (with pump pulse at 400 nm) in different polar solvents have been investigated using femtosecond time-resolved stimulated emission pumping fluorescence depletion (FS TR SEP FD) spectroscopy. In both excited states, three different anisotropy decay laws have been observed for OX750 in different solvents. Only in acetone and formamide could the anisotropy decays of OX750 be described by single-exponential functions, whereas the anisotropy decays have been found to exhibit biexponential behavior in other solvents. The slower anisotropy decay observed in all of the solvents has been assigned to the overall rotational relaxation of OX750 molecules, and a quantitative analysis of this time constant has been performed using the Stokes-Einstein-Debye hydrodynamic theory and the extended charge distribution model developed by Alavi and Waldeck. In both methanol and ethanol, a faster anisotropy decay on the order of picoseconds and a slower anisotropy decay on the hundreds of picoseconds time scale are observed. The most likely explanation for the faster anisotropy involves the rotation of the transition dipole moment in the excited state of OX750 resulting from the electron transfer (ET) reaction taking place from the alcoholic solvents to the OX750 chromophore. As a possible explanation, the wobbling-in-the-cone model has been used to analyze the biexponential anisotropy decays of OX750 in dimethylformamide (DMF) and dimethyl sulfoxide (DMSO). The observed faster anisotropy decays on the hundreds of femtoseconds time scale in DMF and DMSO are ascribed to the wobbling-in-the-cone motion of the ethyl group of OX750, which is sensitive to the strength of the hydrogen bond formed between the solvent and the protonation site of OX750.     

Some structural aspects in binary aqueous mixtures of simple amides from rotational molecular motions

Nuclear magnetic relaxation rates of²D and¹⁴N in binary aqueous mixtures of formamide,N-methylformamide (NMF), andN,N-dimethylformamide (DMF) are reported as a function of the mixture composition. From these intramolecular quadrupolar relaxation data separate rotational correlation times for the two components of the mixture can be determined. The relative variation of the single correlation time as a function of the composition is interpreted in terms of structural changes caused by hydrogen bonding and hydrophobic effects. The results also clearly reflect the expected characteristic variation of these effects on the rotational molecular motions in going from formamide to NMF and DMF. The maximum correlation time retardation of DMF in the aqueous mixture is compared with those of other hydrophobic solvents. A correlation between this maximum retardation and the excess enthalpy of mixing of hydrophobic solvents in aqueous solution can be established graphically.     

Hydrogen bond breaking mechanism and water reorientational dynamics in the hydration layer of lysozyme

The mechanism and the rate of hydrogen bond-breaking in the hydration layer surrounding an aqueous protein are important ingredients required to understand the various aspects of protein dynamics, its function, and stability. Here, we use computer simulation and a time correlation function technique to understand these aspects in the hydration layer of lysozyme. Water molecules in the layer are found to exhibit three distinct bond-breaking mechanisms. A large angle orientational jump of the donor water molecule is common among all of them. In the most common ( approximately 80%) bond-breaking event in the layer, the new acceptor water molecule comes from the first coordination shell (initially within 3.5 A of the donor), and the old acceptor water molecule remains within the first coordination shell, even after the bond-breaking. This is in contrast to that in bulk water, in which both of the acceptor molecules involve the second coordination shell. Additionally, the motion of the incoming and the outgoing acceptor molecules involved is not diffusive in the hydration layer, in contrast to their observed diffusive motion in the bulk. The difference in rotational dynamics between the bulk and the hydration layer water molecules is clearly manifested in the calculated time-dependent angular van Hove self-correlation function ( G(theta, t)) which has a pronounced two-peak structure in the layer, and this can be traced to the constrained translational motion in the layer. The longevity of the surrounding hydrogen bond network is found to be significantly enhanced near a hydrophilic residue.     

Femtosecond photolysis of aqueous formamide

In this work, we investigate the primary photodynamics of aqueous formamide. The formamide was photolyzed using 200 nm femtosecond pulses, and formation of products and their relaxation was followed with approximately 300 fs time resolution using probe pulses covering the range from 193 to 700 nm. Following excitation, the majority of formamide molecules (approximately 80%) converts the electronic excitation energy to vibrational excitation, which effectively is dissipated to the solvent through vibrational relaxation in just a few picoseconds. The vibrational relaxation is observed as a distinct modulation of the electronic absorption spectrum of formamide. The relaxation process is modeled by a simple one-dimensional wavepacket calculation. A smaller fraction of the excited formamide molecules dissociates to the CHO and NH2 radical pairs, of which 50% escape recombination. In addition to the electronic excitation of formamide, we also observe a small contribution from one-photon ionization of formamide and two-photon ionization and dissociation of the water solvent.     

Collective solvent flows around a protein investigated by molecular dynamics simulation

Translational motions of water molecules in various systems equilibrated at room temperature are thought to be diffusive and nondirectional. We performed molecular dynamics simulations of a protein system and showed that the water molecules collectively move around the protein. The motions of two water molecules, which were about 12 A away from each other, are correlated to each other. Such collective motions of water can be regarded as flows around the protein, and the flows exhibited various coherent patterns: fair currents, vortices, and divergent flows. The patterns were highly fluctuating: a set of patterns changed to a different set of patterns within a time scale of 10 ps. Thus, the water motions observed in a scale of length smaller than 12 A and a time scale shorter than 10 ps were nondiffusive, and the motions above these scales were diffusive, where the flows disappeared. The flows near the protein surface had an orientational propensity to be highly parallel to the protein surface, and this propensity gradually vanished with an increment of distance from the protein surface. The divergent patterns of flows, which frequently emerge during the fluctuations of flows, may temporarily cause solvent drying in the vicinity of solutes. The current simulation is supportive of a molecular interaction mechanism that the fluctuations of hydration structure induce attractive interactions between solutes.     

Ultrafast vibrational dynamics of SCN(-) and N3(-) in polar solvents studied by nonlinear infrared spectroscopy

In this paper, we report on our investigation into the vibrational dynamics of the antisymmetric stretching modes of SCN(-) and N(3)(-) in several polar solvents. We used an infrared (IR) pump-probe method to study orientational relaxation processes. In two aprotic solvents (N,N-dimethylformamide (DMF) and dimethyl sulfoxide (DMSO)), the anisotropy decay shows a bimodal feature, whereas in other solvents the anisotropy decay can be fitted well by a single exponential function. We consider that the relative contribution of fast-decaying components is smaller in the other solvents than in DMF and DMSO. We discuss the possible origins of the different anisotropy decay behavior in different solvents. From the three-pulse IR photon echo measurements for SCN(-) and N(3)(-), we found that the time-correlation functions (TCFs) of vibrational frequency fluctuations decay on two different time scales, one of which is less than 100 fs and the other is approximately 3-6 ps. In aprotic solvents, the fast-decaying components of the TCFs on a <100 fs time scale play an important role in the vibrational frequency fluctuation, although the contribution of collective solvent reorganization in aprotic solvents was clearly observed to have small amplitudes. On the other hand, we found that the amplitude of components that decay in a few picoseconds and/or the constant offset of the TCF in protic solvents is relatively large compared with that in aprotic solvents. With the formation and dissociation of hydrogen bonds between ion solute and solvent molecules, the spectra of different solvated species are exchanged with each other and merged into one band. We considered that this exchange may be an origin of slow-decaying components of the TCFs and that the decay of the TCFs corresponds to the time scales of the exchange for protic solvents such as formamide. The mechanism of vibrational frequency fluctuations for the antisymmetric stretching modes of SCN(-) and N(3)(-) is discussed in terms of the difference between protic and aprotic solvents.     

Water Dynamics in Nafion Fuel Cell Membranes: the Effects of Confinement and Structural Changes on the Hydrogen Bond Network

The complex environments experienced by water molecules in the hydrophilic channels of Nafion membranes are studied by ultrafast infrared pump-probe spectroscopy. A wavelength dependent study of the vibrational lifetime of the O-D stretch of dilute HOD in H(2)O confined in Nafion membranes provides evidence of two distinct ensembles of water molecules. While only two ensembles are present at each level of membrane hydration studied, the characteristics of the two ensembles change as the water content of the membrane changes. Time dependent anisotropy measurements show that the orientational motions of water molecules in Nafion membranes are significantly slower than in bulk water and that lower hydration levels result in slower orientational relaxation. Initial wavelength dependent results for the anisotropy show no clear variation in the time scale for orientational motion across a broad range of frequencies. The anisotropy decay is analyzed using a model based on restricted orientational diffusion within a hydrogen bond configuration followed by total reorientation through jump diffusion.     

Structure and dynamics of the solvation of bovine pancreatic trypsin inhibitor in explicit water: a comparative study of the effects of solvent and protein polarizability

To isolate the effects of the inclusion of polarizability in the force field model on the structure and dynamics of the solvating water in differing electrostatic environments of proteins, we present the results of molecular dynamics simulations of the bovine pancreatic trypsin inhibitor (BPTI) in water with force fields that explicitly include polarization for both the protein and the water. We use three model potentials for water and two model potentials for the protein. Two of the water models and one of the protein models are polarizable. A total of six systems were simulated representing all combinations of these polarizable and nonpolarizable protein and water force fields. We find that all six systems behave in a similar manner in regions of the protein that are weakly electrostatic (either hydrophobic or weakly hydrophilic). However, in the vicinity of regions of the protein with relatively strong electrostatic fields (near positively or negatively charged residues), we observe that the water structure and dynamics are dependent on both the model of the protein and the model of the water. We find that a large part of the dynamical dependence can be described by small changes in the local environments of each region that limit the local density of non-hydrogen-bonded waters, precisely the water molecules that facilitate the dynamical relaxation of the water-water hydrogen bonds. We introduce a simple method for rescaling for this effect. When this is done, we are able to effectively isolate the influence of polarizability on the dynamics. We find that the solvating water's relaxation is most affected when both the protein and the water models are polarizable. However, when only one model (or neither) is polarizable, the relaxation is similar regardless of the models used.     

Orientational dynamics of water confined on a nanometer length scale in reverse micelles

The time-resolved orientational anisotropies of the OD hydroxyl stretch of dilute HOD in H(2)O confined on a nanometer length scale in sodium bis(2-ethylhexyl) sulfosuccinate (AOT) reverse micelles are studied using ultrafast infrared polarization and spectrally resolved pump-probe spectroscopy, and the results are compared to the same experiments on bulk water. The orientational anisotropy data for three water nanopool sizes (4.0, 2.4, and 1.7 nm) can be fitted well with biexponential decays. The biexponential decays are analyzed using a wobbling-in-a-cone model that involves fast orientational diffusion within a cone followed by slower, full orientational relaxation. The data provide the cone angles, the diffusion constants for motion within the cones, and the final diffusion constants as a function of the nanopool size. The two processes can be interpreted as a local angular fluctuation of the OD and a global hydrogen bond network rearrangement process. The trend in the relative amplitudes of the long and short exponential decays suggest an increasing rigidity as the nanopool size decreases. The trend in the long decay constants indicates a longer hydrogen bond network rearrangement time with decreasing reverse micelle size. The anisotropy measurements for the reverse micelles studied extrapolate to approximately 0.33 rather than the ideal value of 0.4, suggesting the presence of an initial inertial component in the anisotropy decay that is too fast to resolve. The very fast decay component is consistent with initial inertial orientational motion that is seen in published molecular-dynamics simulations of water in AOT reverse micelles. The angle over which the inertial orientational motion occurs is determined. The results are in semiquantitative agreement with the molecular-dynamics simulations.     

Rotational dynamics of benzene and water in an ionic liquid explored via molecular dynamics simulations and NMR T1 measurements

The rotational dynamics of benzene and water in the ionic liquid (IL) 1-butyl-3-methylimidazolium chloride are studied using molecular dynamics (MD) simulation and NMR T(1) measurements. MD trajectories based on an effective potential are used to calculate the (2)H NMR relaxation time, T(1) via Fourier transform of the relevant rotational time correlation function, C(2R)(t). To compensate for the lack of polarization in the standard fixed-charge modeling of the IL, an effective ionic charge, which is smaller than the elementary charge is employed. The simulation results are in closest agreement with NMR experiments with respect to the temperature and Larmor frequency dependencies of T(1) when an effective charge of ±0.5e is used for the anion and the cation, respectively. The computed C(2R)(t) of both solutes shows a bi-modal nature, comprised of an initial non-diffusive ps relaxation plus a long-time ns tail extending to the diffusive regime. Due to the latter component, the solute dynamics is not under the motional narrowing condition with respect to the prevalent Larmor frequency. It is shown that the diffusive tail of the C(2R)(t) is most important to understand frequency and temperature dependencies of T(1) in ILs. On the other hand, the effect of the initial ps relaxation is an increase of T(1) by a constant factor. This is equivalent to an "effective" reduction of the quadrupolar coupling constant (QCC). Thus, in the NMR T(1) analysis, the rotational time correlation function can be modeled analytically in the form of aexp (-t/τ) (Lipari-Szabo model), where the constant a, the Lipari-Szabo factor, contains the integrated contribution of the short-time relaxation and τ represents the relaxation time of the exponential (diffusive) tail. The Debye model is a special case of the Lipari-Szabo model with a = 1, and turns out to be inappropriate to represent benzene and water dynamics in ILs since a is as small as 0.1. The use of the Debye model would result in an underestimation of the QCC by a factor of 2-3 as a compensation for the neglect of the Lipari-Szabo factor.