Dynamics in supercooled ionic organic liquids and mode coupling theory analysis

Optically heterodyne-detected optical Kerr effect experiments are applied to study the orientational dynamics of the supercooled ionic organic liquids N-propyl-3-methylpyridinium bis(trifluoromethylsulfonyl)imide (PMPIm) and 1-ethyl-3-methylimidazolium tosylate (EMImTOS). The orientational dynamics are complex with relaxation involving several power law decays followed by a final exponential decay. A mode coupling theory (MCT) schematic model, the Sj?gren model, was able to reproduce the PMPIm data very successfully over a wide range of times from 1 ps to hundreds of ns for all temperatures studied. Over the temperature range from room temperature down to the critical temperature Tc of 231 K, the OHD-OKE signal of PMPIm is characterized by the intermediate power law t(-1.00+/-0.04) at short times, a von Schweidler power law t(-0.51+/-0.03) at intermediate times, and a highly temperature-dependent exponential (alpha relaxation) at long times. This form of the decay is identical to the form observed previously for a large number of organic van der Waals liquids. MCT analysis indicates that the theory can explain the experimental data very well for a range of temperatures above Tc, but as might be expected, there are some deviations from the theoretical modeling at temperatures close to Tc. For EMImTOS, the orientational dynamics were studied on the ps time scale in the deeply supercooled region near its glass transition temperature. The orientational relaxation of EMImTOS clearly displays the feature associated with the boson peak at approximately 2 ps, which is the first time domain evidence of the boson peak in ionic organic liquids. Overall, all the dynamical features observed earlier for organic van der Waals liquids using the same experimental technique are also observed for organic ionic liquids.     

Diffusion trapping times and dynamic percolation in an Ising system

We address the problem of diffusion through dynamic Ising network structures using random walkers (RWs) whose net displacements are partitioned into two contributions, arising from (1) transport through neighboring "conducting" clusters and (2) self-diffusion of the site on which the RW finds itself, respectively. At finite temperatures, the conducting clusters in the network exhibit correlated dynamic behavior, making our model system different to most prior published work, which has largely been at the random percolation limit. We also present a novel heuristic scaling analysis for this system that utilizes a new scaling exponent theta(z) for representing RW trapping time as a function of "distance" from the dynamic percolation transition. Simulation results in two-dimensional networks show that when theta(z) = 2, a value found from independent physical arguments, our scaling equations appear to capture universal behavior in the system, at both the random percolation (infinite temperature) and finite temperature conditions studied. This study suggests that the model and the scaling approach given here should prove useful for studying transport in physical systems showing dynamic disorder.     

A comparison between the bond-strength-coordination number fluctuation model and the random walk model of viscosity

We have studied the temperature dependence of the viscosity of some polymeric materials by using both, the bond-strength-coordination number fluctuation model and the random walk model. The results reveal that both models show an excellent agreement with the experimental data. For the random walk model, two equations corresponding to two temperature regimes (low-T and high-T) separated by the critical temperature T _c, which is difficult to determine, are needed to describe the temperature dependence of the viscosity of a fragile system, whereas for the bond-strength-coordination number fluctuation model, a single equation with clear physical meaning describes the temperature dependence of the viscosity of both, the fragile and strong systems. We have also studied the relationship between the normalized temperature range of cooperativity and the fragility index. A theoretical expression for the relationship has been derived based on the bond-strength-coordination number fluctuation model. The comparison with the experimental data shows a good agreement, leading to the conclusion that the kinetic properties of glass forming liquids and the cooperativity of molecular relaxations are correlated.     

Rotational dynamics in supercooled water from nuclear spin relaxation and molecular simulations

Structural dynamics in liquid water slow down dramatically in the supercooled regime. To shed further light on the origin of this super-Arrhenius temperature dependence, we report high-precision (17)O and (2)H NMR relaxation data for H(2)O and D(2)O, respectively, down to 37 K below the equilibrium freezing point. With the aid of molecular dynamics (MD) simulations, we provide a detailed analysis of the rotational motions probed by the NMR experiments. The NMR-derived rotational correlation time τ(R) is the integral of a time correlation function (TCF) that, after a subpicosecond librational decay, can be described as a sum of two exponentials. Using a coarse-graining algorithm to map the MD trajectory on a continuous-time random walk (CTRW) in angular space, we show that the slowest TCF component can be attributed to large-angle molecular jumps. The mean jump angle is ～48° at all temperatures and the waiting time distribution is non-exponential, implying dynamical heterogeneity. We have previously used an analogous CTRW model to analyze quasielastic neutron scattering data from supercooled water. Although the translational and rotational waiting times are of similar magnitude, most translational jumps are not synchronized with a rotational jump of the same molecule. The rotational waiting time has a stronger temperature dependence than the translation one, consistent with the strong increase of the experimentally derived product τ(R)?D(T) at low temperatures. The present CTRW jump model is related to, but differs in essential ways from the extended jump model proposed by Laage and co-workers. Our analysis traces the super-Arrhenius temperature dependence of τ(R) to the rotational waiting time. We present arguments against interpreting this temperature dependence in terms of mode-coupling theory or in terms of mixture models of water structure.     

Measuring diffusivity in supercooled liquid nanoscale films using inert gas permeation. II. Diffusion of Ar, Kr, Xe, and CH4 through methanol

We present an experimental technique to measure the diffusivity of supercooled liquids at temperatures near their T(g). The approach uses the permeation of inert gases through supercooled liquid overlayers as a measure of the diffusivity of the supercooled liquid itself. The desorption spectra of the probe gas are used to extract the low temperature supercooled liquid diffusivities. In the preceding companion paper, we derived equations using ideal model simulations from which the diffusivity could be extracted using the desorption peak times for isothermal or peak temperatures for temperature programmed desorption experiments. Here, we discuss the experimental conditions for which these equations are valid and demonstrate their utility using amorphous methanol with Ar, Kr, Xe, and CH(4) as probe gases. The approach offers a new method by which the diffusivities of supercooled liquids can be measured in the experimentally challenging temperature regime near the glass transition temperature.     

Moderately and strongly supercooled liquids: a temperature-derivative study of the primary relaxation time scale

The primary relaxation time scale tau(T) derived from the glass forming supercooled liquids (SCLs) is discussed within ergodic-cluster Gaussian statistics, theoretically justified near and above the glass-transformation temperature T(g). An analysis is given for the temperature-derivative data by Stickel et al. on the steepness and the curvature of tau(T). Near the mode-coupling-theory (MCT) crossover T(c), these derivatives separate by a kink and a jump, respectively, the moderately and strongly SCL states. After accounting for the kink and the jump, the steepness remains a piecewise conitnuous function, a material-independent equation for the three fundamental characteristic temperatures, T(g), T(c), and the Vogel-Fulcher-Tamman (VFT) T(0), is found. Both states are described within the heterostructured model of solidlike clusters parametrized in a self-consistent manner by a minimum set of observable parameters: the fragility index, the MCT slowing-down exponent, and the chemical excess potential of Adam and Gibbs model (AGM). Below the Arrhenius temperature, the dynamically and thermodynamically stabilized clusters emerge with a size of around of seven to nine and two to three molecules above and close to T(g) and T(c), respectively. On cooling, the main transformation of the moderately into the strongly supercooled state is due to rebuilding of the cluster structure, and is attributed to its rigidity, introduced through the cluster compressibility. It is shown that the validity of the dynamic AGM (dynamically equivalent to the standard VFT form) is limited by the strongly supercooled state (T(g) < T < T(c)) where the superrigid cooperative rearranging regions are shown to be well-chosen parametrized solidlike clusters. Extension of the basic parameter set by the observable kinetic and diffusive exponents results in prediction of a subdiffusion relaxation regime in SCLs that is distinct from that established for amorphous polymers.     

Viscosity and structure configuration properties of equilibrium and supercooled liquid cobalt

The shear viscosity of liquid cobalt at the pressure p = 1.5 bar and at the temperatures corresponding to equilibrium liquid and supercooled liquid states is measured experimentally and evaluated by means of molecular dynamics simulations. Further, the shear viscosity is also calculated within the microscopic theoretical model. Comparison of our experimental, simulation and theoretical results with other available data allows one to examine the issue about the correct temperature dependence of the shear viscosity of liquid cobalt. A strong correlation between the viscosity and the configuration entropy of liquid cobalt over the considered temperature range is found, which can be taken into account by the Rosenfeld’s model.     

Structural relaxation and rheological response of a driven amorphous system

The interplay between the structural relaxation and the rheological response of a simple amorphous system {a 80:20 binary Lennard-Jones mixture [W. Kob and H. C. Andersen, Phys. Rev. Lett. 73, 1376 (1994)]} is studied via molecular dynamics simulations. In the quiescent state, the model is well known for its sluggish dynamics and a two step relaxation of correlation functions at low temperatures. An ideal glass transition temperature of Tc=0.435 has been identified in the previous studies via the analysis of the system's dynamics in the framework of the mode coupling theory of the glass transition [W. Kob and H. C. Andersen, Phys. Rev. E 51, 4626 (1995)]. In the present work, we focus on the question whether a signature of this ideal glass transition can also be found in the case where the system's dynamics is driven by a shear motion. Indeed, the following distinction in the structural relaxation is found: In the supercooled state, the structural relaxation is dominated by the shear at relatively high shear rates gamma, whereas at sufficiently low gamma the (shear-independent) equilibrium relaxation is recovered. In contrast to this, the structural relaxation of a glass is always driven by shear. This distinct behavior of the correlation functions is also reflected in the rheological response. In the supercooled state, the shear viscosity eta decreases with increasing shear rate (shear thinning) at high shear rates, but then converges toward a constant as the gamma is decreased below a (temperature-dependent) threshold value. Below Tc, on the other hand, the shear viscosity grows as eta proportional, etax 1/gamma, suggesting a divergence at gamma=0. Thus, within the accessible observation time window, a transition toward a nonergodic state seems to occur in the driven glass as the driving force approaches zero. As to the flow curves (stress versus shear rate), a plateau forms at low shear rates in the glassy phase. A consequence of this stress plateau for Poiseuille-type flows is demonstrated.     

粘度法研究胶团的球-棒转变

本文基于粘度测量求得了溴化十四烷基吡啶胶团水溶液的特性粘度随盐量(0.5—2.0 mol L~(-1)NaBr), 温度(20°—40 ℃)和活性剂浓度(3-25 g L~(-1))的变化。利用Missel等提出的热力学模型,由此求算了支配球-棒转变的平衡常数K和多分散胶团的重均聚集数, 当盐浓度为0.5和1.0 mol L~(-1)时, 由粘度法在不同温度下得到的K值以及在30 ℃下随活性剂浓度的变化, 与光散射的相应结果相当一致。此外, 在不同温度下求得的参数K在一定的胶团浓度范围内保持为常量。这符合模型的要求。盐量增至2 mol L~(-1)时, K明显地具有浓度依赖性,并且两法给出的K和值都相差较大。上述模型不再适用。对于符合球-棒转变的胶团体系, 从特性粘度具有浓度依赖性出发的粘度方法是研究胶团长大的简单而可靠的方法, 能够提供很多重要的信息。     

DFT based Monte Carlo simulations of poly(9,9-dialkylfluorene-2,7-diyl) polymers in solution

The dilute solution properties of poly(9,9-dihexylfluorene-2,7-diyl) (PDHF) were studied by coupled SEC/light scattering and MALDI-TOF over a large molecular weight (MW) span ranging from PDHF oligomers (1-8-mer) to high MW polymer. The results were compared with Monte Carlo simulations based on realistic PDHF models obtained from X-ray data and density functional theory (DFT) calculations and with a DFT based Kratky-Porod-Benoit-Doty (KPBD) worm-like chain. The simulations called "selective random walk" (SRW) and the corresponding "selective self-avoiding random walk" (SSAW) explicitly take into account the rotationally labile bonds between the fluorene units in that four distinct torsion angles (+/-37.5 and +/-143 degrees) between the units are chosen randomly. The simulations better account than the KPBD model for the experimental data obtained by us and others for various poly(9,9-dialkylfluorene-2,7-diyl) polymers but still give somewhat larger values for the radii of gyration and hydrodynamic volumes. The torsion angle selectivity of the SRW and SSAW simulations predict long chain sections punctuated by sudden sharp loops.     

Anomalies of Liquid Water

Water at low temperatures (T < 300 K) reveals a multitude of unusual physico-chemical properties, which originate from the strong directivity of the hydrogen bonds. This directivity forces the molecules into a tetrahedral arrangement of first neighbors and thus leads to a poor packing efficiency. All anomalies of the liquid become more pronounced in the metastable range below the melting pressure curve. In the following article, methods for the investigation of the supercooled range are presented. A few anomalous static and dynamic properties are described and compared with the properties of normal liquids. The experimental results are discussed in the context of computer simulations and the more recent theories on cold water. It can be concluded that the anomalies of water, for instance the density maximum of the liquid phase and the initial decrease in viscosity with pressure, are confined to the pT-range: T < 300 K, p ? 200 MPa. At these pressures, lowering the temperature leads to an unexpected behavior of almost all physical properties. They seem to approach a phase transition at T ? (T_M ? 50 K), the nature of which has not yet been fully characterized. This phenomenon has hitherto been observed only in supercooled water.—In the future it will be possible to use supercooled aqueous solutions for kinetic studies, thus expanding the available dynamic range for the investigation of aqueous solutions considerably.     

Simplified crossover droplet model for adsorption of pure fluids in slit pores

We present a generalized crossover (GC) model for the excess adsorption of pure fluids at a flat solid-liquid interface, which reproduces scaling behavior of the excess adsorption in the critical region and is reduced to the classical, van der Waals-type analytical model far away from the bulk critical point. In developing this model, we used the density-functional theory (DFT) approach for the order parameter profile calculations with a generalized corresponding states model for the local free-energy density. The GC DFT model well represents the available experimental adsorption data for Kr/graphite, C2H4/graphite, C3H8/graphite, CO2/silica, and SF6/graphite systems in the entire density range 0 < rho < or = 3rhoc and temperatures up to 1.7Tc. In the critical region 0.5 rhoc < r < or = 1.5rhoc and T < or = 1.15Tc, the GC DFT model is consistent with the predictions of the asymptotic renormalization-group crossover model for the critical adsorption in a semi-infinite system developed earlier. For the excess adsorption on the critical isochore, both theories predict a scaling-law behavior Gamma proportional tau(-nu+beta), but fail to reproduce a "critical depletion" of the excess adsorption along the critical isochore of the SF6/graphite system near Tc. We show that an anomalous decrease of adsorption observed in this system at tau = T/Tc - 1 < 10(-2) can be explained by finite-size effect and develop a simplified crossover droplet (SCD) model for the excess adsorption in a slit pore. With the effective size of the pore of L = 50 nm, the SCD model reproduces all available experimental data for SF6/graphite, including the critical isochore data where tau-->0, within experimental accuracy. At L > xib (where xib is a bulk correlation length) the SCD model is transformed into the GC DFT model for semi-infinite systems. Application of the SCD model to the excess adsorption of carbon dioxide on the silica gel is also discussed.     

An Alternate Solution of Fluorescence Recovery Kinetics after Spot-Bleaching for Measuring Diffusion Coefficients. 2. Diffusion of Fluorescein in Aqueous Sucrose Solutions

The traditional analysis of the fluorescence recovery kinetics after spot bleaching yields expressions for the diffusion coefficient of the probe that are not suitable for linear fittings. In a previous work we developed an improved recovery function that is a better alternative for data analysis. To illustrate its application to real cases and compare it with the previous data treatment, we measured the time response of fluorescein in aqueous sucrose solutions, covering the unsaturated and the supercooled region, where decoupling between diffusion and viscosity is observed. The results are compared with the mobility of different types of solutes in aqueous sucrose solutions and are discussed in terms of the classical hydrodynamic model.     

Heat capacity of tetrahydrofuran clathrate hydrate and of its components, and the clathrate formation from supercooled melt

We report a thermodynamic study of the formation of tetrahydrofuran clathrate hydrate by explosive crystallization of water-deficient, near stoichiometric, and water-rich solutions, as well as of the heat capacity, C(p), of (i) supercooled tetrahydrofuran-H2O solutions and of the clathrate hydrate, (ii) tetrathydrofuran (THF) liquid, and (iii) supercooled water and the ice formed on its explosive crystallization. In explosive freezing of supercooled solutions at a temperature below 257 K, THF clathrate hydrate formed first. The nucleation temperature depends on the cooling rate, and excess water freezes on further cooling. The clathrate hydrate melts reversibly at 277 K and C(p) increases by 770 J/mol K on melting. The enthalpy of melting is 99.5 kJ/mol and entropy is 358 J/mol K. Molar C(p) of the empty host lattice is less than that of the ice, which is inconsistent with the known lower phonon frequency of H2O in the clathrate lattice. Analysis shows that C(p) of THF and ice are not additive in the clathrate. C(p) of the supercooled THF-H2O solutions is the same as that of water at 247 K, but less at lower temperatures and more at higher temperatures. The difference tends to become constant at 283 K. The results are discussed in terms of the hydrogen-bonding changes between THF and H2O.     

Compatibility studies on solutions of polymer blends by viscometric and ultrasonic techniques

Detailed viscometric and ultrasonic velocity studies have been conducted on solutions of blends of poly(methyl methacrylate) with poly(vinyl acetate), poly(vinyl chloride) with poly(vinyl acetate) and poly(methyl methacrylate) with polystyrene over an extended range of concentrations and temperatures in toluene, chlorobenzene and toluene respectively. The plots of both absolute viscosity and ultrasonic velocity vs composition deviate from linearity according to the degree of compatibility of polymer blends, at all concentrations and temperatures. The curves for compatible systems are linear. These investigations offer an entirely new approach to the study of the compatibility of polymer blends.     

Viscosity Investigations on the Binary Systems of (1 ChCl:2 Ethylene Glycol) DES and Methanol or Ethanol

In this study, the viscosity behavior of two mixtures of Ethaline (1 ChCl:2 ethylene glycol) with either methanol or ethanol were investigated over the temperature range of 283.15–333.15 K at atmospheric pressure. The measured viscosities of neat Ethaline, methanol, and ethanol showed reliable agreement with the corresponding reported literature values. The mixture viscosities were modeled by an Arrhenius-like model to determine the behavior of viscosity with respect to temperature. The data were also modeled by the four well-known mixture viscosity models of Grunberg–Nissan, Jouyban–Acree, McAllister, and Preferential Solvation. All of the model results were reliable, with the Jouyban–Acree and Preferential Solvation models showing the most accurate agreement with the experimental measurements. The Jones–Dole viscosity model was also investigated for the measured viscosities, and by analyzing the results of this model, strong interactions among Ethaline and the alcohol molecules were proposed for both systems. As a final analysis, viscosity deviations of the investigated systems were calculated to study the deviations of the viscosity behaviors with respect to ideal behavior. Both systems showed negative viscosity deviations at all of the investigated temperatures, with the negative values tending towards zero, and hence more ideal behavior, with increasing temperatures. Moreover, in order to correlate the calculated viscosity deviations, the Redlich–Kister model was successfully used for both systems and at each investigated temperature.     

Diffusion-viscosity decoupling in supercooled aqueous trehalose solutions

The diffusional mobility of disodium fluorescein has been measured in supercooled aqueous solutions of trehalose, a widely used cryoprotectant disaccharide. The results were analyzed on the basis of the classical continuum hydrodynamic theory (Stokes-Einstein relationship) and compared with results for the diffusion and electrical conductivity of other ionic and nonionic solutes in trehalose and sucrose aqueous solutions. Disodium fluorescein obeys the classical model over a restricted range of inverse reduced temperatures, T g/ T, scaled by the glass transition temperature. Decoupling in neutral solutes takes place at higher values of T g/ T, while in ionic solutes it occurs all over the range of T g/ T studied, as observed for the water mobility in supercooled sugar solutions.     

Glasslike behavior in aqueous electrolyte solutions

When salts are added to water, generally the viscosity increases, suggesting that the ions increase the strength of the water's hydrogen-bond network. However, infrared pump-probe measurements on electrolyte solutions have found that ions have no influence on the rotational dynamics of water molecules, implying no enhancement or breakdown of the hydrogen-bond network. Here, we report optical Kerr effect and dielectric relaxation spectroscopic measurements, which have enabled us to separate the effects of rotational and transitional motions of the water molecules. These data show that electrolyte solutions behave like a supercooled liquid approaching a glass transition in which rotational and translational molecular motions are decoupled. It is now possible to understand previously conflicting viscosity data, nuclear magnetic resonance relaxation, and ultrafast infrared spectroscopy in a single unified picture.     

Uptake measurements of ethanol on ice surfaces and on supercooled aqueous solutions doped with nitric acid between 213 and 243 K

Uptake of ethanol either on pure frozen ice surfaces or supercooled solutions doped with HNO3 (0.63 and 2.49 wt %) has been investigated using a coated wall flow tube coupled to a mass spectrometric detection. The experiments were conducted over the temperature range of 213-243 K. Uptake of ethanol on these surfaces was always found to be totally reversible whatever were the experimental conditions. The number of ethanol molecules adsorbed per surface unit was conventionally plotted as a function of ethanol concentration in the gas phase and subsequently analyzed using Langmuir's model. The amount of ethanol molecules taken up on nitric acid doped-ice surfaces was found to increase largely with increasing nitric acid concentrations. For example at 223 K, and for an ethanol gas-phase concentration of 1x10(13) molecules cm3, the number of adsorbed molecules are (in units of molecules cm-2): approximately 1.3x10(14) on pure ice; approximately 1.4x10(15) on ice doped with HNO3 0.63 wt %; approximately 7.5x10(15) on ice doped with HNO3, 2.49 wt %, i.e. 60 times larger than on pure ice. Since, according to the shape of the isotherms, the adsorption did not proceed beyond monolayer coverage, the enormous increase of ethanol uptake was explained by considering its dissolution in either a supercooled liquid layer (T<230 K) or a liquid solution (T>230 K). The formation of both was indeed favored by the presence of the HNO3. Our experimental results suggest that the amount of ethanol dissolved in such supercooled solutions follows Henry's law and that the Henry's law constants at low temperatures, i.e., 223-243 K, can be estimated by extrapolation from higher temperatures. Such supercooled solutions which exist in the troposphere either in deep convective clouds or in mixed clouds for temperature above 233 K, might be responsible for the scavenging of large amounts of soluble species, such as nitric and sulfuric acids, oxygenated VOCs including alcohols, carboxylic acids, and formaldehyde.     

Insight into the structural and transport properties of methyl and benzyl triphenyl phosphonium based deep eutectic solvents using molecular dynamics simulations

Molecular dynamics simulation of twelve phosphonium based DESs based on methyl triphenyl phosphonium bromide and benzyl triphenyl phosphonium chloride as hydrogen bond acceptors (HBA) with glycerol, ethylene glycol, triethylene glycol and trifluoroacetamide as hydrogen bond donors (HBDs) are carried out with the Generalized Amber Force Field between 25 °C and 95 °C at 1 atm pressure. Isobaric-isothermal ensemble is used for predicting the structural properties, namely, radial distribution function and coordination number. Further, time resolved trajectory data is correlated to predict the transport properties such as, self-diffusivity, binary-diffusivity, ionic conductivity and viscosity. The temperature dependence of predicted density is observed with 4.96 × 10⁻² absolute average deviations from experimental dataset. Further, the self-diffusivities are predicted using Einstein's approach and ionic conductivity is calculated from Nernst-Einstein equation as well as Green-Kubo formulation. Further, the viscosity is also predicted using the Green-Kubo formulation. The predicted ionic conductivities and viscosities are fitted in the Vogel-Fulcher-Tammann (VFT) equation to calculate the VFT parameters for the respective DESs.