Local order and dynamics in supercooled water: a study by IR spectroscopy and molecular dynamic simulations

Micron-sized water droplets in a cryogenic flow tube were probed by IR spectroscopy. The analysis of the IR spectra suggests that there is a relative increase of about 30% in the fraction, f(L), of low density domains in water on cooling over the temperature range between 300 and 240 K. The results derived from the experiments agree qualitatively with those of molecular dynamics (MD) simulations in terms of the increase in the f(L) values. The MD simulations show that the intensities of the mode at about 100 cm(-1) for the molecules in the low density domains are reduced in comparison to the average, while the intensities and frequencies of the librational mode at 600 cm(-1) are increased. Furthermore, the reorientations (dielectric relaxation times) in these domains are found to be somewhat slower, pointing to the fact that these low density "cages" live longer than the average local molecular environments in supercooled water.     

PEG-POSS Star Molecules Blended in Polyurethane with Flexible Hard Segments: Morphology and Dynamics

A star polymer with a polyhedral oligomeric silsesquioxanne (POSS) core and poly(ethylene glycol) (PEG) vertex groups is incorporated in a polyurethane with flexible hard segments in-situ during the polymerization process. The blends are studied in terms of morphology, molecular dynamics, and charge mobility. The methods utilized for this purpose are scanning electron and atomic force microscopies (SEM, AFM), X-ray diffraction (XRD), differential scanning calorimetry (DSC), and to a larger extent dielectric relaxation spectroscopy (DRS). It is found that POSS reduces the degree of crystallinity of the hard segments. Contrary to what was observed in a similar system with POSS pendent along the main chain, soft phase calorimetric glass transition temperature drops as a result of plasticization, and homogenization of the soft phase by the star molecules. The dynamic glass transition though, remains practically unaffected, and a hypothesis is formed to resolve the discrepancy, based on the assumption of different thermal and dielectric responses of slow and fast modes of the system. A relaxation α′, slower than the bulky segmental α and common in polyurethanes, appears here too. A detailed analysis of dielectric spectra provides some evidence that this relaxation has cooperative character. An additional relaxation g, which is not commonly observed, accompanies the Maxwell Wagner Sillars interfacial polarization process, and has dynamics similar to it. POSS is found to introduce conductivity and possibly alter its mechanism. The study points out that different architectures of incorporation of POSS in polyurethane affect its physical properties by different mechanisms.     

Molecular motions in amorphous ibuprofen as studied by broadband dielectric spectroscopy

The molecular mobility of amorphous ibuprofen has been investigated by broadband dielectric relaxation spectroscopy (DRS) covering a temperature range of more than 200 K. Four different relaxation processes, labeled as alpha, beta, gamma, and D, were detected and characterized, and a complete relaxation map was given for the first time. The gamma-process has activation energy E a = 31 kJ.mol (-1), typical for local mobility. The weak beta-relaxation, observed in the glassy state as well as in the supercooled state was identified as the genuine Johari-Goldstein process. The temperature dependence of the relaxation time of the alpha-process (dynamic glass transition) does not obey a single VFTH law. Instead two VFTH regimes are observed separated by a crossover temperature, T B = 265 K. From the low temperature VFTH regime, a T g (diel) (tau =100 s) = 226 K was estimated, and a fragility or steepness index m = 93, was calculated showing that ibuprofen is a fragile glass former. The D-process has a Debye-like relaxation function but the temperature dependence of relaxation time also follows the VFTH behavior, with a Vogel temperature and a pre-exponential factor which seem to indicate that its dynamics is governed by the alpha-process. It has similar features as the Debye-type process observed in a variety of associating liquids, related to hydrogen bonding dynamics. The strong tendency of ibuprofen to form hydrogen bonded aggregates such as dimers and trimers either cyclic or linear which seems to control in particular the molecular mobility of ibuprofen was confirmed by IR spectroscopy, electrospray ionization mass spectrometry, and MD simulations.     

Enhanced translational diffusion of rubrene in sucrose benzoate

The translational diffusion of rubrene in the fragile molecular glass former, sucrose benzoate (SB) (fragility index m approximately 94), has been studied from T(g)+6 K to T(g)+71 K(T(g)=337 K) by using the technique of holographic fluorescence recovery after photobleaching. In the temperature range of the measurements, the translational relaxation functions were observed to decay exponentially, indicating that Fick's law of diffusion governs the translational motion of rubrene in sucrose benzoate. The value of the translational diffusion coefficient D(T) obtained from the 1e time of the translational relaxation function varied from 5.3 x 10(-15) cm2 s(-1) at 343 K to 5.0x10(-9) cm2 s(-1) at 408 K. The temperature dependence of D(T) for diffusion of rubrene in SB is compared with that of the viscosity and the dielectric relaxation time tau(D) of SB. The temperature dependence of D(T) is weaker than that of Teta for T<1.2T(g) but tracks the reciprocal of the dielectric relaxation time 1tau(D) for 1.05T(g)相似文献   

Chain Dynamics in Solid Polymers and Polymerizing Systems as Revealed by Broadband Dielectric Spectroscopy

Broadband Dielectric Spectroscopy (BDS) provides a direct, powerful means of studying the chain dynamics of solid polymers. The present work discusses experimental data, phenomenological theories, molecular theories and recent molecular dynamics simulations of the relaxation functions and average relaxation times for the multiple dielectric relaxations observed for amorphous polymers. Examples of current applications of BDS are given, including studies of the dynamics of ultra-thin polymer films and real-time studies of (i) the crystallization behaviour of bulk polymers and (ii) thermo-polymerizing systems.     

Hydration Dynamics of Water near an Amphiphilic Model Peptide at Low Hydration Levels: A Dielectric Relaxation Study

A dielectric relaxation study of aqueous solutions of the amphiphilic model peptide N‐acetyl‐leucine amide (NALA) at 298 K over a wide range of hydration levels is presented. The experiments range from states where water builds up several hydration layers to states where single water molecules or small water clusters are shared by several NALA molecules. The dielectric spectra reveal two modes on the 10 and 100 ps timescales. These are largely broadened with regard to the Lorentzian shape caused by simple Debye‐type relaxation, and are well described by the Kohlrausch–Williams–Watts stretched exponential function. The fast mode is assigned to water reorientation comprising bulk water as well as hydration water. Even when all water molecules are in contact with the solute, this fast component is dominant, and its mean relaxation time is retarded by less than a factor of two relative to neat water. The amplitude of the slow process is far higher than expected for the dipolar reorientation of the solute. The observations are consistent with results from molecular dynamics simulations for a similar model peptide reported in the literature. They suggest that the slow relaxation mode is mainly founded in peptide–water dipolar couplings, with some additional contribution from slowly reorienting hydration water molecules. The results are discussed with regard to the hydration dynamics of proteins and the interpretation of dielectric spectra of protein solutions.     

Impact of anisotropy on the structure and dynamics of ionic liquids: a computational study of 1-butyl-3-methyl-imidazolium trifluoroacetate

The complex ionic network of 1-butyl-3-methyl-imidazolium trifluoroacetate was simulated by means of the molecular dynamics methods over a time period of 100 ns. The influence of the anisotropy of the shape and charge distribution of both the cations and the anions on the local (molecular) and global (collective) structure and dynamics is analyzed. The distance-dependent g coefficients of the orientational probability function g(r,Omega) were found to be an excellent way to interpret local structure. Thereby, the combination and interrelation of individual g coefficients elucidate the mutual orientation. Dynamics at the molecular level is characterized by the time correlation function of the center-of-mass corrected molecular dipole moment mucm. Upon uniting the set of molecular dipoles to a single collective rotational dipole moment, MD, dynamics on a global level is studied. Decomposing into subsets of cations and anions respective self terms as well as the prominent cross term can be extracted. This decomposition also enables a detailed peak assignment in dielectric spectra.     

On the nature of the high-frequency relaxation in a molecular glass former: a joint study of glycerol by field cycling NMR, dielectric spectroscopy, and light scattering

Fast field cycling (1)H NMR relaxometry is applied to determine the dispersion of spin-lattice relaxation time T(1)(omega) of the glass former glycerol in broad temperature (75-360 K) and frequency (10 kHz-30 MHz) ranges. The relaxation data are analyzed in terms of a susceptibility chi(")(omega) proportional, variantomegaT(1)(omega), related to the second rank (l=2) molecular orientational correlation function. Broadband dielectric spectroscopic results suggest the validity of frequency temperature superposition above the glass transition temperature T(g). This allows to combine NMR data of different temperatures into a single master curve chi(")(omegatau(alpha)) that extends over 15 decades in reduced frequency omegatau(alpha), where tau(alpha) is the structural alpha-relaxation time. This master curve is compared with the corresponding ones from dielectric spectroscopy (l=1) and depolarized light scattering (l=2). At omegatau(alpha)<1, NMR susceptibility is significantly different from both the dielectric and light scattering results. At omegatau(alpha)>1, there rather appears a difference between the susceptibilities of rank l=1 and l=2. Specifically, at omegatau(alpha)>1, where the susceptibility is dominated by the so-called excess wing, the NMR and light scattering spectra (both l=2) rather coincide with each other and are about three times more intense than the dielectric (l=1) spectrum. This is explained by assuming that the high frequency dynamics correspond to only small-angle excursions. Below T(g), dielectric and NMR susceptibility compare well and exhibit an exponential temperature dependence.     

Dielectric relaxation of aqueous solutions of hydrophilic versus amphiphilic peptides

We report on molecular dynamics simulations of the frequency-dependent dielectric relaxation spectra at room temperature for aqueous solutions of a hydrophilic peptide and an amphiphilic peptide at two concentrations. We find that only the high-concentration amphiphilic peptide solution exhibits an anomalous dielectric increment over that of pure water, while the hydrophilic peptide exhibits a significant dielectric decrement. The dielectric component analysis carried out by decomposing these peptide solutions into peptide, hydration layer, and outer layer(s) of water clearly shows the presence of a unique dipolar component with a relaxation time scale on the order of approximately 25 ps (compared to the bulk water time scale of approximately 11 ps) that originates from the interaction between the hydration layer water and the outer layer(s) of water. Results obtained from the dielectric component analysis further show the emergence of a distinct and much lower frequency relaxation process for the high-concentration amphiphilic peptide compared to the hydrophilic peptide due to strong peptide dipolar couplings to all constituents, accompanied by a slowing of the structural relaxation in all water layers, giving rise to time scales close to approximately 1 ns. We suggest that the molecular origin of the dielectric relaxation anomalies is due to frustration in the water network arising from the amphiphilic chemistry of the peptide that does not allow it to reorient on the picosecond time scale of bulk water motions. This explanation is consistent with the idea of the "slaving" of residue side chain motions to protein surface water, and furthermore offers the possibility that the anomalous dynamics observed from a number of spectroscopies arises at the interface of hydrophobic and hydrophilic domains on the protein surface.     

Solvent response and dielectric relaxation in supercooled butyronitrile

We have measured the dynamics of solvation of a triplet state probe, quinoxaline, in the glass-forming dipolar liquid butyronitrile near its glass transition temperature T(g)=95 K. The Stokes shift correlation function displays a relaxation time dispersion of considerable magnitude and the optical linewidth changes along the solvation coordinate in a nonmonotonic fashion. These features are characteristic of solvation in viscous solvents and clearly indicate heterogeneous dynamics, i.e., spatially distinct solvent response times. Using the dielectric relaxation data of viscous butyronitrile as input, a microscopic model of dipolar solvation captures the relaxation time, the relaxation dispersion, and the amplitude of the dynamical Stokes shift remarkably well.     

时域反射光谱法研究酰胺-醇混合物的介电驰豫

Using time domain reflectometry (TDR),dielectric relaxation studies were carried out on binary mixtures of amides (N-methylformamide (NMF) and N,N-dimethylformamide (DMF)) with alcohols (1-butanol,1-pentanol,1-hexanol,1-heptanol,1-octanol,and 1-decanol) for various concentrations over the frequency range from 10 MHz to 10 GHz at 303 K. The Kirkwood correlation factor and excess dielectric constant properties were determined and discussed to yield information on the molecular interactions of the systems. The relaxation time varied with the chain length of alcohols and substituted amides were noticed. The Bruggeman plot shows a deviation from linearity. This deviation was attributed to some sort of molecular interaction which may take place between the alcohols and substituted amides. The excess static permittivity and excess inverse relaxation time values varied from negative to positive for all the systems indicating that the solute-solvent interaction existed between alcohols and substituted amides for all the dynamics of the mixture.     

时域反射光谱法研究酰胺-醇混合物的介电驰豫

Using time domain reflectometry （TDR）, dielectric relaxation studies were carded out on binary mixtures of amides （N-methylformamide （NMF） and N,N-dimethylformamide （DMF）） with alcohols （1-butanol, 1-pentanol, 1- hexanol, 1-heptanol, 1-octanol, and 1-decanol） for various concentrations over the frequency range from 10 MHz to 10 GHz at 303 K. The Kirkwood correlation factor and excess dielectric constant properties were determined and discussed to yield information on the molecular interactions of the systems. The relaxation time varied with the chain length of alcohols and substituted amides were noticed. The Bruggeman plot shows a deviation from linearity. This deviation was attributed to some sort of molecular interaction which may take place between the alcohols and substituted amides. The excess static permittivity and excess inverse relaxation time values varied from negative to positive for all the systems indicating that the solute-solvent interaction existed between alcohols and substituted amides for all the dynamics of the mixture.     

Solvation dynamics and electric field relaxation in an imidazolium-PF6 ionic liquid: from room temperature to the glass transition

Time-resolved phosphorescence spectra and anisotropy of quinoxaline were measured in an ionic liquid, 1-butyl-3-methylimidazolium hexafluorophosphate (BMIM-HFP), in its supercooled state near the glass-transition temperature. The solvation dynamics results are compared with the rotational motion of the probe and with the dielectric behavior of the neat ionic liquid. The dynamics in the viscous state are highly dispersive and show a super-Arrhenius temperature dependence, as typical for glass-forming materials. Combined with room-temperature results, solvation dynamics is observed to follow the structural relaxation times in terms of eta/T for more than 10 decades, from subnanoseconds at room temperature to seconds near the glass-transition temperature T(g). The dielectric modulus relaxation follows this trend only for temperatures T > 1.2T(g) and departs significantly from eta/T in the 1.1T(g) > T > T(g) range. This deviation is reminiscent of the enhanced translational diffusion or fractional Stokes-Einstein behavior observed in many fragile supercooled liquids. Because the electric field relaxation in BMIM-HFP includes dc conductivity, this correlation function involves translational motion and thus displays the effect of enhanced diffusivity. A microscopic model is required for rationalizing the decoupling of solvation dynamics from the longitudinal time scales and the limitation of this effect to the viscous regime with T < 1.2T(g).     

Solute rotational dynamics at the water liquid/vapor interface

The rotational dynamics of a number of diatomic molecules adsorbed at different locations at the interface between water and its own vapors are studied using classical molecular dynamics computer simulations. Both equilibrium orientational and energy correlations and nonequilibrium orientational and energy relaxation correlations are calculated. By varying the dipole moment of the molecule and its location, and by comparing the results with those in bulk water, the effects of dielectric and mechanical frictions on reorientation dynamics and on rotational energy relaxation can be studied. It is shown that for nonpolar and weekly polar solutes, the equilibrium orientational relaxation is much slower in the bulk than at the interface. As the solute becomes more polar, the rotation slows down and the surface and bulk dynamics become similar. The energy relaxation (both equilibrium and nonequilibrium) has the opposite trend with the solute dipole (larger dipoles relax faster), but here again the bulk and surface results converge as the solute dipole is increased. It is shown that these behaviors correlate with the peak value of the solvent-solute radial distribution function, which demonstrates the importance of the first hydration shell structure in determining the rotational dynamics and dependence of these dynamics on the solute dipole and location.     

Correlation of nonexponentiality with dynamic heterogeneity from four-point dynamic susceptibility χ4(t) and its approximation χT(t)

Various properties of vitrifying liquids are correlated with the dispersity of the dynamics, the latter reflected in the magnitude of the nonexponentiality parameter, β(K), describing the distribution of relaxation times. These properties include the mean relaxation time, τ(α), the fragility, and the dynamic crossover. The correlations with β(K) are observed in both experimental data and the results from molecular dynamics simulations on Lennard-Jones (LJ) type systems. Another, rather obvious property to correlate with β(K) is the dynamic heterogeneity, which can be quantified from the number of molecules, N(c), dynamically correlated over a time span τ(α). For a given LJ system, N(c) can be rigorously calculated and we find that it does indeed correlate with β(K) over a range of thermodynamic conditions. However, the analysis of experimental data for a broad range of real materials, wherein an approximation is required to obtain N(c), reveals the absence of any relationship between N(c) and β(K) among different materials.     

Decomposing total IR spectra of aqueous systems into solute and solvent contributions: a computational approach using maximally localized Wannier orbitals

The theoretical principles underpinning the calculation of infrared spectra for condensed-phase systems in the context of ab initio molecular dynamics have been recently developed in literature. At present, most ab initio molecular dynamics calculations are restricted to relatively small systems and short simulation times. In this paper we devise a method that allows well-converged results for infrared spectra from ab initio molecular dynamics simulations using small systems and short trajectories characteristic of simulations typically performed in practice. We demonstrate the utility of our approach by computing the imaginary part of the dielectric constant epsilon"(omega) for H2O and D2O in solid and liquid phases and show that it compares well with experimental data. We further demonstrate that maximally localized Wannier orbitals can be used to separate the individual contributions of different molecular species to the linear spectrum of complex systems. The new spectral decomposition method is shown to be useful in present-day ab initio molecular dynamics calculations to compute the magnitude of the "continuous absorption" generated by excess protons in aqueous solutions with good accuracy even when other species present in the solutions absorb strongly in the same frequency window.     

OH-stretch vibrational relaxation of HOD in liquid to supercritical D2O

The population relaxation of the OH-stretching vibration of HOD diluted in D2O is studied by time-resolved infrared (IR) pump-probe spectroscopy for temperatures of up to 700 K in the density range 12 1 OH stretching transition with a 200 fs laser pulse centered at approximately 3500 cm(-1). Above 400 K these spectra show no indication of spectral diffusion after pump-probe delays of 0.3 ps. Over nearly the entire density range and for sufficiently high temperatures (T > 360 K), the vibrational relaxation rate constant, kr, is strictly proportional to the dielectric constant, epsilon, of water. Together with existing molecular dynamics simulations, this result suggests a simple linear dependence of kr on the number of hydrogen-bonded D2O molecules. It is shown that, for a given temperature, an isolated binary collision model is able to adequately describe the density dependence of vibrational energy relaxation even in hydrogen-bonded fluids. However, dynamic hydrogen bond breakage and formation is a source of spectral diffusion and affects the nature of the measured kr. For sufficiently high temperatures when spectral diffusion is much faster than energy transfer, the experimentally observed decays correspond to ensemble averaged population relaxation rates. In contrast, when spectral diffusion and vibrational relaxation occur on similar time scales, as is the case for ambient conditions, deviations from the linear kr(epsilon) relation occur because the long time decay of the v = 1 population is biased to slower relaxing HOD molecules that are only weakly connected to the hydrogen bond network.     

Dielectric behavior of some small ketones as ideal polar molecules

The dielectric behaviors of some small symmetric ketone molecules, including acetone, 3-pentanone, cyclopentanone, 4-heptanone, and cyclohexanone, were investigated as a function of temperature (T) over a wide frequency range from 50 MHz (3.14 × 10(8) s(-1), in angular frequency) to 3 THz (1.88 × 10(13) s(-1)). The temperature dependencies of the rotational diffusion times (τ(r)) determined using (17)O NMR spin-lattice relaxation time (T(1)) measurements and viscosities of the ketones were also examined. The obtained temperature dependencies of the parameters for the ketones were compared with those of ideal polar molecules, which obey the Stokes-Einstein-Debye (SED) relationship without the formation of intermolecular dimeric associations and without orientational correlations between dipoles (molecular axes), that is, free rotation. Kirkwood correlation factors (g(K)) of only acetone and 3-pentanone were close to unity over a wide temperature range, whereas those of other ketones were obviously less than unity. These results revealed that no correlations exist between the rotational motions of dipoles in acetone and 3-pentanone, as expected in ideal polar molecules. However, other ketones exhibited orientational correlations in their dipoles because of dipole-dipole interactions via antiparallel configurations. Furthermore, because acetone and 3-pentanone satisfied the SED relationship and because their microscopic dielectric relaxation times (τ(μ)), which were calculated from the determined dielectric relaxation times (τ(D)) via the relationship τ(μ) = τ(D)g(K)(-1), were identical to 3τ(r) and were proportional to Vη(k(B)T)(-1) over the wide temperature range examined, where V, k(B), and η represent the effective molecular volume, Boltzmann's constant, and the viscosity of the liquid molecules, respectively, these two ketone molecules behave as ideal polar molecules. In addition, other ketones not significantly larger than acetone and 3-pentanone in molecular size likely form dimeric intermolecular associations with antiparallel cyclic configurations, which demonstrates the g(K) values less than unity.