Thermodynamic scaling of the viscosity of van der Waals, H-bonded, and ionic liquids

Viscosities eta and their temperature T and volume V dependences are reported for seven molecular liquids and polymers. In combination with literature viscosity data for five other liquids, we show that the superpositioning of relaxation times for various glass-forming materials when expressed as a function of TV(gamma), where the exponent gamma is a material constant, can be extended to the viscosity. The latter is usually measured to higher temperatures than the corresponding relaxation times, demonstrating the validity of the thermodynamic scaling throughout the supercooled and higher T regimes. The value of gamma for a given liquid principally reflects the magnitude of the intermolecular forces (e.g., steepness of the repulsive potential); thus, we find decreasing gamma in going from van der Waals fluids to ionic liquids. For some strongly H-bonded materials, such as low molecular weight polypropylene glycol and water, the superpositioning fails, due to the nontrivial change of chemical structure (degree of H bonding) with thermodynamic conditions.     

Translational diffusion in sucrose benzoate near the glass transition: probe size dependence in the breakdown of the Stokes-Einstein equation

The translational diffusion coefficient D(trans) for rubrene, 9,10-bis(phenylethynyl)anthracene (BPEA), and tetracene in the fragile molecular glass-former sucrose benzoate (SB) (Tg=337 K) was studied as a function of temperature from Tg+3 K to Tg+71 K by use of the holographic fluorescence recovery after photobleaching technique. The values of D(trans) vary by five to six orders of magnitude in this temperature range. Contrary to the predictions of the Stokes-Einstein equation, the temperature dependence of probe diffusion in SB over the temperature range of the measurements is weaker than that of T/eta, where eta is the shear viscosity. In going from the crossover temperature Tx approximately 1.2Tg to Tg, D(trans)eta/T increases by factors of 2.4+/-0.2 decades for rubrene, 3.4+/-0.2 decades for BPEA, and 3.8+/-0.4 decades for tetracene. The decoupling between probe diffusion in SB and viscosity is characterized by the scaling law D(trans) approximately T/eta(xi), with xi=0.621 for tetracene, 0.654 for BPEA, and 0.722 for rubrene. Data for probe diffusion in SB are combined with data from the literature for probe diffusion in ortho-terphenyl and alphaalphabeta-tris(naphthyl)benzene in a plot of enhancement versus the relative probe size parameter rho(m)=(m(p)m(h))(1/3), where m(p) and m(h) are, respectively, the molecular weights of the probe and host solvent. The plot clearly shows a sharp increase in enhancement of translational diffusion at rho(m) approximately 1. By applying temperature shifts, D(trans) for probe diffusion in SB and the dielectric relaxation time tau(D) can be superimposed on a single master curve based on the Williams-Landel-Ferry equation. This suggests that the dynamics of probe diffusion in SB is described by the scaling relationship D(trans) approximately 1/tau(D)(T+DeltaT), where tau(D)(T+DeltaT) is the temperature-shifted dielectric relaxation time. The results from this study are discussed within the context of dynamic heterogeneity in glass-forming liquids.     

Assessment of phenomenological models for viscosity of liquids based on nonequilibrium atomistic simulations of copper

The shear viscosity of liquid copper is studied using nonequilibrium molecular-dynamics simulations under planar shear flow conditions. We examined variation of viscosity as function of shear rate at a range of pressures (ca. 0 - 40 GPa). We analyzed these results using eight different phenomenological models and find that the observed non-Newtonian behavior is best described by the Powell-Eyring (PE) model: eta(gamma) = (eta(0)-eta(infinity))sinh(-1)(taugamma)(taugamma) + eta(infinity), where gamma is the shear rate. Here eta(0) (the zero-shear-rate viscosity) extracted from the PE fit is in excellent agreement with available experimental data. The relaxation time tau from the PE fit describes the shear response to an applied stress. This provides the framework for interpreting the shear flow phenomena in complex systems, such as liquid metal and amorphous metal alloys.     

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)相似文献   

Pressure dependence of viscosity

We reanalyze the pressure dependence of viscosity of liquids of constant composition under isothermal conditions. Based exclusively on very general considerations concerning the relationship between viscosity and "free volume," we show that, at moderate values of pressure, viscosity increases, as a rule, with increasing pressure, provided the liquid is in stable or metastable (undercooled) equilibrium states. However, even if the behavior of the viscosity is governed by free volume effects, deviations from a positive pressure dependence are possible, when the liquid's thermal expansion coefficient is negative. We derive an equation that allows one to quantitatively determine the pressure dependence of viscosity, which requires, in the simplest case, only the knowledge of the temperature dependence of viscosity at constant pressure, the thermal expansion coefficient, and the isothermal compressibility of the liquid. As an example, the negative pressure dependence of water in the range of temperatures 0-4 degrees C and of several silicate liquids, such as albite, jadeite, dacite, basalts, etc., could be explained in such a way. Other glass-forming liquids initially (for moderate pressures) show a positive pressure dependence of viscosity that changes to a negative one when subjected to high (approximately GPa) isostatic pressure. A detailed analysis of water and already mentioned silicate melts at GPa pressures shows that, in addition to free volume effects, other pressure induced structural transformations may have to be accounted for in a variety of cases. By this reason, the theoretical analysis is extended (i) in order to describe the pressure dependence of viscosity for systems that are in frozen-in thermodynamic nonequilibrium states (glasses, i.e., undercooled liquids below the glass transition temperature Tg) and (ii) to systems which undergo, in addition to variations of the free volume, pressure induced changes of other structural parameters. In such cases a decrease of viscosity with increasing pressure may occur, in principle, even if the thermal expansion coefficient is positive. In this way, the present analysis grants a general tool to estimate the pressure dependence of viscosity and supposedly settles the controversy in the current literature.     

Sound attenuation, shear viscosity, and mutual diffusivity behavior in the nitroethane-cyclohexane critical mixture

The shear viscosity eta(s), mutual diffusion coefficient D, and ultrasonic attenuation spectra of the nitroethane-cyclohexane mixture of critical composition have been measured at various temperatures near the critical temperature T(c). The relaxation rate of order parameter fluctuations resulting from a combined evaluation of the eta(s) and D data follows power law behavior with the theoretical exponent and with the large amplitude Gamma(o)=(156+/-2)x10(9) s(-1). The ultrasonic spectra have been evaluated in terms of a critical contribution and a noncritical background contribution. The amplitude of the former exhibits a temperature dependence, in conformity with a temperature dependence in the adiabatic coupling constant (|g| = 0.064 near T(c) and 0.1 at T-T(c)=3 K). If the variation of the critical amplitude with T is taken into account the experimental attenuation coefficient data display a scaling function which nicely fits to the theoretical prediction from the Bhattacharjee-Ferrell dynamic scaling model [R. A. Ferrell and J. K. Bhattacharjee, Phys. Rev. A 31, 1788 (1985)].     

Crystal growth kinetics exhibit a fragility-dependent decoupling from viscosity

In this paper we establish the temperature dependence of the kinetic coefficient associated with crystal growth into the supercooled liquid for a wide range of organic and inorganic materials. We show that the kinetic coefficient for crystal growth scales with the shear viscosity eta as eta(-xi) and that the exponent depends systematically on the fragility of the liquid. The greater the fragility (i.e., deviation away from an Arrhenius temperature dependence for eta), the larger the difference 1-xi. We argue that this breakdown in scaling between the crystal growth kinetics and the viscosity is a manifestation of heterogeneous dynamics in supercooled liquids. In addition, we show that the absolute growth rate at intermediate viscosities is correlated with the entropy difference between the liquid and the crystal.     

Relationship between viscosity coefficients and volumetric properties using a scaling concept for molecular and ionic liquids

In this work, a scaling concept based on relaxation theories of the liquid state was combined with a relation previously proposed by the authors to provide a general framework describing the dependency of viscosity on pressure and temperature. Namely, the viscosity-pressure coefficient (partial differentialeta/partial differentialp)T was expressed in terms of a state-independent scaling exponent, gamma. This scaling factor was determined empirically from viscosity versus Tvgamma curves. New equations for the pressure- and temperature-viscosity coefficients were derived, which are of considerable technological interest when searching for appropriate lubricants for elastohydrodynamic lubrication. These relations can be applied over a broad range of thermodynamic conditions. The fluids considered in the present study are linear alkanes, pentaerythritol ester lubricants, polar liquids, associated fluids, and several ionic liquids, compounds selected to represent molecules of different sizes and with diverse intermolecular interactions. The values of the gamma exponent determined for the fluids analyzed in this work range from 1.45 for ethanol to 13 for n-hexane. In general, the pressure-viscosity derivative is well-reproduced with the values obtained for the scaling coefficient. Furthermore, the effects of volume and temperature on viscosity can be quantified from the ratio of the isochoric activation energy to the isobaric activation energy, Ev/Ep. The values of gamma and of the ratio Ev/Ep allow a classification of the compounds according to the effects of density and temperature on the behavior of the viscosity.     

Shear and longitudinal viscosity of non-ionic C12E8 aqueous solutions

We present an extensive set of measurements of steady shear viscosity (eta degrees(s)), longitudinal elastic modulus (M'), and ultrasonic absorption (alpha) in the one-phase isotropic liquid region of the non-ionic surfactant C12E8 aqueous solutions. Within a given temperature interval, this phase extends along the entire surfactant concentration range that could be fully covered in the experiments. In agreement with previous studies, the overall results support the presence of two separated intervals of concentration corresponding to different structural properties. In the surfactant-rich region the temperature dependence of eta degrees(s) follows an equation characteristic of glass-like systems. The ultrasonic absorption spectra show unambiguous evidence of viscoelastic behavior that can be described by a Cole-Cole relaxation formula. In this region, when both the absorption and the frequency are scaled by the static shear viscosity (eta degrees(s)), the scaled attenuation reduces to a single universal curve for all temperatures and concentrations. In the water-rich region the behavior of eta degrees(s), M', and alpha are more complex and reflect the presence of dispersed aggregates whose size increases with temperature and concentration. At these concentrations the ultrasonic spectra are characterized by a multiple decay rate. The high-frequency tail falls in the same frequency range seen at high surfactant content and exhibits similar behaviors. This contribution is ascribed to the mixture of hydrophilic terminations and water present at the micellar interfaces that resembles the condition of a concentrated polymer solution. An additional low-frequency contribution is also observed, which is ascribed to the exchange of water molecules and/or surfactant monomers between the aggregates and the bulk solvent region.     

Pressure effects on the alpha and alpha' relaxations in polymethylphenylsiloxane

In some polymers, in addition to the usual structural alpha relaxation, a slower alpha' relaxation is observed with a non-Arrhenius temperature dependence. In order to understand better the molecular origin of this alpha' relaxation in poly(methylphenylsiloxane) (PMPS) we have studied, for the first time, the pressure dependence of its relaxation time, together with the usual temperature dependence, by means of dynamic light scattering (DLS). For the same material the alpha relaxation was also studied by means of DLS and dielectric spectroscopy (DS) in broad temperature and pressure ranges. We find that the temperature dependence of both alpha and alpha' relaxation times, at all pressures studied, can be described by a double Vogel-Fulcher-Tammann (VFT) law. The pressure dependence of the characteristic temperatures Tg (glass transition temperature) and T0 (Vogel temperature) as well as the activation volumes for both alpha and alpha' processes are very similar, indicating, that both relaxation processes originate from similar local molecular dynamics. Additionally, for both alpha and alpha' relaxations the combined temperature and pressure dependences of the relaxation times can be described using a parameter Gamma=rhon/T with the same value of the exponent n.     

Temperature dependence of the dynamic viscosity of liquid mercury

The dependence of the dynamic viscosity of mercury on temperature is calculated and expressed in terms of the cluster associate model, based on the Boltzmann distribution and with normalization at the melting point. The resulting refined equation for mercury viscosity adequately reflects this dependence over the range of the liquid state, including the critical point.     

Effect of chemical structure on the isobaric and isochoric fragility in polychlorinated biphenyls

Pressure-volume-temperature data, along with dielectric relaxation measurements, are reported for a series of polychlorinated biphenyls (PCB), differing in the number of chlorine atoms on their phenyl rings. Analysis of the results reveals that with increasing chlorine content, the relaxation times of the PCB become governed to a greater degree by density rho relative to the effect of temperature T. This result is consistent with the respective magnitudes of the scaling exponent gamma yielding superpositioning of the relaxation times measured at various temperatures and pressures, when plotted versus rho(gamma)/T. While at constant (atmospheric) pressure, fragilities for the various PCB are equivalent, the fragility at constant volume varies inversely with chlorine content. Evidently, the presence of bulkier chlorine atoms on the phenyl rings magnifies the effect which the density has on the relaxation dynamics.     

Study on laminar viscosity and zero shear viscosity of latex systems

The flowing nature and rheological properties of polymethyl methacrylate latex systems in a coaxial cylinder viscometer were studied on the basis of laminar shear flow model and rheological experimental data. The physical meaning of laminar viscosity (eta(i,j)) and zero shear viscosity (eta(0)) were described. We assumed that laminar shear flows depended on position and shear time, so microrheological parameters were the function of position and shear time. eta(i,j) was the viscosity of any shear sheet i between two neighboring laminar shear flows at time t; j was denoted as j=t/Deltat; and Deltat was the interacting time of two particles or two laminar shear flows. tau(i,j) and gamma(i,j) were shear stress and shear rate of any shear sheet i at j moment. According to Newton regulation tau(i,j)=eta(i,j)gamma(i,j), apparent viscosity eta(a) should be a statistically mean value of j shear sheets laminar viscosity at j moment, i.e., eta(a)= summation operator(i=j)eta(i,j)gamma(i,j)/ summation operator(i=j)gamma(i,j). eta(0) was defined as shear viscosity between a laminar shear flow and a still fluid surface, i.e., eta(0)=(tau(i,j)/gamma(i,j))(j-i-->0). These new ideas described above may be helpful in the study of the micromechanisms of latex particle systems and worthy of more research.     

Primary and secondary relaxations in supercooled eugenol and isoeugenol at ambient and elevated pressures: dependence on chemical microstructure

Dielectric loss spectra of two glass-forming isomers, eugenol and isoeugenol, measured at ambient and elevated pressures in the normal liquid, supercooled, and glassy states are presented. The isomeric chemical compounds studied differ only by the location of the double bond in the alkyl chain. Above the glass transition temperature T(g), the dielectric loss spectra of both isomers exhibit an excess wing on the high frequency flank of the loss peak of the alpha relaxation and an additional faster gamma process at the megahertz frequency range. By decreasing temperature below T(g) at ambient pressure or by elevating pressure above P(g), the glass transition pressure, at constant temperature, the excess wing of isoeugenol shifts to lower frequencies and is transformed into a secondary beta-loss peak, while in eugenol it becomes a shoulder. These spectral features enable the beta-relaxation time tau(beta) to be determined in the glassy state. These changes indicate that the excess wings in isoeugenol and eugenol are similar and both are secondary beta relaxations that are not resolved in the liquid state. While in both isoeugenol and eugenol the loss peak of the beta relaxation in the glassy state and the corresponding excess wing in the liquid state shifts to lower frequencies on elevating pressure, the locations of their gamma relaxation show little change with increasing pressure. The different pressure sensitivities of the excess wing and gamma relaxation are further demonstrated by the nearly perfect superposition of the alpha-loss peak together with excess wing from the data taken at ambient pressure and at elevated pressure (and higher temperature so as to have the same alpha-peak frequency), but not the gamma-loss peak in both isoeugenol and eugenol. On physical aging isoeugenol, the beta-loss peak shifts to lower frequencies, but not the gamma relaxation. Basing on these experimental facts, the faster gamma relaxation is a local intramolecular process involving a side group and the slower beta relaxation mimics the structural alpha relaxation in behavior, involves the entire molecule and satisfies the criteria for being the Johari-Goldstein beta relaxation. Analysis and interpretation of the spectra utilizing the coupling model further demonstrate that the excess wings seen in the equilibrium liquid states of these two isomers are their genuine Johari-Goldstein beta relaxation.     

PVT and oscillatory tests to analyze pressure effects on polypropylene/Rodrun LC3000 blends: Determination of the pressure dependency of the viscosity

A combined analysis of Pressure-Volume-Temperature (PVT), Dynamic Mechanical Thermal Analysis (DMTA) and oscillatory flow measurements for blends of a polypropylene (PP) with a commercial liquid crystalline polymer (Rodrun LC3000) is presented. This analysis allows the determination of the pressure-viscosity coefficient b = ∂lnη₀/∂P. This coefficient depends on the Rodrun LC3000 content, increasing with it and is of the same order of magnitude as values reported for several commercial polymers showing a similar dependence of the viscosity on pressure. The analysis of the pressure dependence of Tg (related to b) leads to the conclusion that the number of segments involved in the glass transition of PP increases with the Rodrun LC3000 content, thus demonstrating that the polymers are not totally immiscible. As far as the authors know, this is the first time that the dependence of the viscosity on the pressure has been reported for thermoplastic/liquid crystalline polymer blends.     

Ion dynamics under pressure in an ionic liquid

We investigate ion dynamics under pressure in the ionic liquid 1-butyl-1-methylpyrrolidinium bis[oxalate]borate (BMP-BOB) by conductivity relaxation measurements in the temperature range 123-300 K and varying pressures from 0.1 MPa up to 0.5 GPa. We report on the influence of pressure on the relaxation times and on the spectral shape of the conductivity relaxation process. We also analyze the pressure dependence of the glass transition temperature and find that the dynamic response under pressure in this ionic liquid shows remarkable similarities to nonionic glass formers. The main relaxation process shows temperature-pressure superposition while a secondary relaxation process, very weakly depending on pressure, is observed. The spectral shape of the main relaxation broadens with increasing pressure or decreasing temperature, but is found to be the same when the relaxation time is the same, independently of the particular pressure and temperature values.     

Solvent viscosity dependence of the folding rate of a small protein: distributed computing study

By using distributed computing techniques and a supercluster of more than 20,000 processors we simulated folding of a 20-residue Trp Cage miniprotein in atomistic detail with implicit GB/SA solvent at a variety of solvent viscosities (gamma). This allowed us to analyze the dependence of folding rates on viscosity. In particular, we focused on the low-viscosity regime (values below the viscosity of water). In accordance with Kramers' theory, we observe approximately linear dependence of the folding rate on 1/gamma for values from 1-10(-1)x that of water viscosity. However, for the regime between 10(-4)-10(-1)x that of water viscosity we observe power-law dependence of the form k approximately gamma(-1/5). These results suggest that estimating folding rates from molecular simulations run at low viscosity under the assumption of linear dependence of rate on inverse viscosity may lead to erroneous results.     

Critical behavior of 2,6-dimethylpyridine-water: measurements of specific heat, dynamic light scattering, and shear viscosity

The specific heat C(p) at constant pressure, the shear viscosity eta(s), and the mutual diffusion coefficient D of the 2,6-dimethylpyridine-water mixture of critical composition have been measured in the homogeneous phase at various temperatures near the lower critical demixing temperature T(c). The amplitude of the fluctuation correlation length xi(0)=(0.198+/-0.004) nm has been derived from a combined evaluation of the eta(s) and D data. This value is in reasonable agreement with the one obtained from the amplitude A(+)=(0.26+/-0.01) J(g K) of the critical term in the specific heat, using the two-scale-factor universality relation. Within the limits of error the relaxation rate Gamma of order parameter fluctuations follows power law with the theoretical universal exponent and with the amplitude Gamma=(25+/-1)x10(9) s(-1). No indications of interferences of the critical fluctuations with other elementary chemical reactions have been found. A noteworthy result is the agreement of the background viscosity eta(b), resulting from the treatment of eta(s) and D data, with the viscosity eta(s)(nu=0) extrapolated from high-frequency viscosity data. The latter have been measured in the frequency range of 5-130 MHz using a novel shear impedance spectrometer.     

Dielectric secondary relaxations in polypropylene glycols

Broadband dielectric measurements of polypropylene glycol of molecular weight M(w)=400 g / mol (PPG 400) were carried out at ambient pressure over the wide temperature range from 123 to 353 K. Three relaxation processes were observed. Besides the structural alpha relaxation, two secondary relaxations, beta and gamma, were found. The beta process was identified as the true Johari-Goldstein relaxation by using a criterion based on the coupling model prediction. The faster gamma relaxation, well separated from the primary process, undoubtedly exhibits the anomalous behavior near the glass transition temperature (T(g)) which is reflected in the presence of a minimum of the temperature dependence of the gamma-relaxation time. We successfully applied the minimal model [Dyre and Olsen, Phys. Rev. Lett. 91, 155703 (2003)] to describe the entire temperature dependence of the gamma-relaxation time. The asymmetric double-well potential parameters obtained by Dyre and Olsen for the secondary relaxation of tripropylene glycol at ambient pressure were modified by fitting to the minimal model at lower temperatures. Moreover, we showed that the effect of the molecular weight of polypropylene glycol on the minimal model parameters is significantly larger than that of the high pressure. Such results can be explained by the smaller degree of hydrogen bonds formed by longer chain molecules of PPG at ambient pressure than that created by shorter chains of PPG at high pressure.