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.     

Density scaling of the transport properties of molecular and ionic liquids

Casalini and Roland [Phys. Rev. E 69, 062501 (2004); J. Non-Cryst. Solids 353, 3936 (2007)] and other authors have found that both the dielectric relaxation times and the viscosity, η, of liquids can be expressed solely as functions of the group (TV?(γ)), where T is the temperature, V is the molar volume, and γ a state-independent scaling exponent. Here we report scaling exponents γ, for the viscosities of 46 compounds, including 11 ionic liquids. A generalization of this thermodynamic scaling to other transport properties, namely, the self-diffusion coefficients for ionic and molecular liquids and the electrical conductivity for ionic liquids is examined. Scaling exponents, γ, for the electrical conductivities of six ionic liquids for which viscosity data are available, are found to be quite close to those obtained from viscosities. Using the scaling exponents obtained from viscosities it was possible to correlate molar conductivity over broad ranges of temperature and pressure. However, application of the same procedures to the self-diffusion coefficients, D, of six ionic and 13 molecular liquids leads to superpositioning of poorer quality, as the scaling yields different exponents from those obtained with viscosities and, in the case of the ionic liquids, slightly different values for the anion and the cation. This situation can be improved by using the ratio (D∕T), consistent with the Stokes-Einstein relation, yielding γ values closer to those of viscosity.     

On the density scaling of pVT data and transport properties for molecular and ionic liquids

In this work, a general equation of state (EOS) recently derived by Grzybowski et al. [Phys. Rev. E 83, 041505 (2011)] is applied to 51 molecular and ionic liquids in order to perform density scaling of pVT data employing the scaling exponent γ(EOS). It is found that the scaling is excellent in most cases examined. γ(EOS) values range from 6.1 for ammonia to 13.3 for the ionic liquid [C(4)C(1)im][BF(4)]. These γ(EOS) values are compared with results recently reported by us [E. R. López, A. S. Pensado, M. J. P. Comu?as, A. A. H. Pádua, J. Fernández, and K. R. Harris, J. Chem. Phys. 134, 144507 (2011)] for the scaling exponent γ obtained for several different transport properties, namely, the viscosity, self-diffusion coefficient, and electrical conductivity. For the majority of the compounds examined, γ(EOS) > γ, but for hexane, heptane, octane, cyclopentane, cyclohexane, CCl(4), dimethyl carbonate, m-xylene, and decalin, γ(EOS) < γ. In addition, we find that the γ(EOS) values are very much higher than those of γ for alcohols, pentaerythritol esters, and ionic liquids. For viscosities and the self-diffusion coefficient-temperature ratio, we have tested the relation linking EOS and dynamic scaling parameters, proposed by Paluch et al. [J. Phys. Chem. Lett. 1, 987-992 (2010)] and Grzybowski et al. [J. Chem. Phys. 133, 161101 (2010); Phys. Rev. E 82, 013501 (2010)], that is, γ = (γ(EOS)/φ) + γ(G), where φ is the stretching parameter of the modified Avramov relation for the density scaling of a transport property, and γ(G) is the Gru?neisen constant. This relationship is based on data for structural relaxation times near the glass transition temperature for seven molecular liquids, including glass formers, and a single ionic liquid. For all the compounds examined in our much larger database the ratio (γ(EOS)/φ) is actually higher than γ, with the only exceptions of propylene carbonate and 1-methylnaphthalene. Therefore, it seems the relation proposed by Paluch et al. applies only in certain cases, and is really not generally applicable to liquid transport properties such as viscosities, self-diffusion coefficients or electrical conductivities when examined over broad ranges of temperature and pressure.     

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.     

Experimental investigation of the spontaneous wetting of polymers and polymer blends

The dynamics of polymeric liquids and mixtures spreading on a solid surface have been investigated on completely wetting and partially wetting surfaces. Drops were formed by pushing the test liquid through a hole in the underside of the substrate, and the drop profiles were monitored as the liquid wet the surface. Silicon surfaces coated with diphenyldichlorosilane (DPDCS) and octadecyltrichlorosilane (OTS) were used as wetting and partial wetting surfaces, respectively, for the fluids we investigated. The response under complete and partial wetting conditions for a series of polypropylene glycols (PPG) with different molecular weights and the same surface tension could be collapsed onto a single curve when scaling time based on the fluid viscosity, the liquid-vapor surface tension, and the radius of a spherical drop with equivalent volume. A poly(ethylene glycol) (PEG300) and a series of poly(ethylene oxide-rand-propylene oxide) copolymers did not show the same viscosity scaling when spread on the partially wetting surface. A combined model incorporating hydrodynamic and molecular-kinetic wetting models adequately described the complete wetting results. The assumptions in the hydrodynamic model, however, were not valid under the partial wetting conditions in our work, and the molecular-kinetic model was chosen to describe our results. The friction coefficient used in the molecular-kinetic model exhibited a nonlinear dependence with viscosity for the copolymers, indicating a more complex relationship between the friction coefficient and the fluid viscosity.     

Chromatography with Supercritical Fluids

Chromatography with supercritical fluids unites the features of both gas chromatography and liquid chromatography, yet retains special characteristics of its own. The diffusion coefficient and particularly the viscosity of fluid phases may approach values for low-pressure gases, while the solvent power may be similar to that of liquids. However, with supercritical fluids it is possible to control chromatographic separations very effectively by pressure programming, since the solubility increases with increasing density. Temperature programming, on the other hand, can have the opposite effect to that in gas- or liquid-chromatography since the density decreases with increasing temperature at a given pressure. Supercritical fluid chromatography is primarily of interest for the separation of higher molecular weight compounds. The efficiency of this method of separation is demonstrated on several homologous series. Thus, a styrene oligomer with nominal M_w=2200 can be resolved by a pressure and temperature program into 40 species.     

Effect of the pressure on the viscosities of ionic liquids: Experimental values for 1-ethyl-3-methylimidazolium ethylsulfate and two bis(trifluoromethyl-sulfonyl)imide salts

A new falling-body viscometer has been implemented to measure viscosity of liquids in a temperature range from (313.15 to 363.15) K at pressures up to 150 MPa. The accuracy of the viscometer was verified after comparing experimental results of squalane with previous literature data finding an average absolute deviation lower than 1.5%. With this device, we have measured viscosity values for three ionic liquids: 1-ethyl-3-methylimidazolium ethylsulfate, 1-butyl-1-methylpyrrolidinium bis(trifluoro-methylsulfonyl)imide and 1-(2-methoxyethyl)-1-methylpyrrolidinium bis(trifluoromethylsulfonyl)imide within the temperature and pressure ranges noted above. The experimental values were correlated as a function of temperature and pressure with four different equations. In addition, we have analysed the pressure–viscosity derived properties for these fluids and for other five ionic liquids using literature values.     

Transport Coefficients Of Ar—Kr Mixtures by Molecular Dynamics Computer Simulation

Abstract

Equilibrium molecular dynamics computer simulations have been used to determine the transport coefficients of model Ar—Kr mixtures, which are represented by Lennard-Jones pair potentials with Lorentz—Berthelot rules for the cross-species interactions. The component self-diffusion and mutual-diffusion coefficients are calculated from time correlation functions and mean square displacements. Time correlation functions are used to evaluate the shear and bulk viscosity, thermal conductivity and the thermal diffusion coefficient (Soret/Dufour coefficient). In the case of the thermal transport coefficients, the partial enthalpy of the two species is required at each state point to define the heat flux rigorously. We obtain this and the partial volume (and species resolved chemical potential) using particle-exchange (and particle insertion) techniques implemented in separate [NPT] simulations at the same state point.

The viscoelasticity of the fluids is characterised by the relaxation times for bulk and shear stress relaxation. The results are for dense liquids close to the triple point temperature and density. Agreement with experiment and previous simulation is particularly good for the density of the mixtures, the shear modulus, shear viscosity, shear stress relaxation time and thermal conductivity. As for the single component noble gas fluids (simulated and experiment) there is a significant qualitative difference in the temperature and, for mixtures, composition dependence of the bulk viscosity.     

Thermodynamic interpretation of the scaling of the dynamics of supercooled liquids

The recently discovered scaling law for the relaxation times, tau(T,upsilon) = I(Tupsilon(gamma)), where T is temperature and upsilon the specific volume, is derived by a revision of the entropy model of the glass transition dynamics originally proposed by Avramov [J. Non-Cryst. Solids 262, 258 (2000)]. In this modification the entropy is calculated by an alternative route. The resulting expression for the variation of the relaxation time with T and upsilon is shown to accurately fit experimental data for several glass-forming liquids and polymers over an extended range encompassing the dynamic crossover. From this analysis, which is valid for any model in which the relaxation time is a function of the entropy, we find that the scaling exponent gamma can be identified with the Gruneisen constant.     

含氰基离子液体的合成、表征及流变性质研究

摘要合成、 表征了一系列新的含氰基咪唑类离子液体. 测定了该离子液体的密度、 熔点及溶解性等物理性质, 研究了其在稳态、 瞬态和动态条件下的流变行为. 结果表明, 当剪切速率在0.1～50 s^-1范围内时, 其粘度不随剪切速率的变化而变化, 但随温度升高而降低, 粘流活化能随取代基长度变化呈现规律性变化. 对于1-丁基-3-氰乙基咪唑六氟磷酸盐离子液体, 维持剪切速率不变时, 其剪切应力和粘度均不随时间变化, 且随着温度的升高而降低; 在动态条件下, 在线性粘弹区, 复合粘度和损耗模量G″ 随温度升高而降低. 关键词     

Critical dynamics of the binary system nitroethane/3-methylpentane: relaxation rate and scaling function

Shear viscosity and dynamic light scattering measurements as well as ultrasonic spectrometry studies of the nitroethane/3-methylpentane mixture of critical composition have been performed at various temperatures near the critical temperature, T(c). A combined evaluation of the shear viscosity and mutual diffusion coefficient data yielded the amplitude, xi(0), of the fluctuation correlation length, xi, assumed to follow power law, and the relaxation rate, Gamma, or order parameter fluctuations. The latter was found to follow power law with the theoretical universal exponent. The amplitudes xi(0) = 0.23 +/- 0.02 nm and Gamma(0) = (125 +/- 5) x 10(9) s(-1) nicely agree with literature values. Using the relaxation rates resulting from the viscosity and diffusion coefficient data, the scaling function has been calculated assuming the ultrasonic spectra to be composed of a critical part and a noncritical background contribution. The experimental scaling function fits well to the predictions of the Bhattacharjee-Ferrell dynamic scaling model with scaled half-attenuation frequency, Omega(BF)1/2= 2.1. The amplitude of the sonic spectra yields the amount |g| = 0.26 of the adiabatic coupling constant, g, in fair agreement with -0.29 from another thermodynamic relation.     

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.     

Thermodynamic scaling of diffusion in supercooled Lennard-Jones liquids

The manner in which the intermolecular potential u(r) governs structural relaxation in liquids is a long standing problem in condensed matter physics. Herein, we show, in agreement with recent experimental results, that diffusion coefficients for simulated Lennard-Jones m-6 liquids (8 < or = m < or = 36) in normal and moderately supercooled states are a unique function of the variable rhogamma/T, where rho is density and T is temperature. The scaling exponent gamma is a material specific constant whose magnitude is related to the steepness of the repulsive part of u(r), evaluated around the distance of closest approach between particles probed in the supercooled regime. Approximations of u(r) in terms of inverse power laws are also discussed.     

Critical fluctuations near the consolute point of n-pentanol-nitromethane. An ultrasonic spectrometry, dynamic light scattering, and shear viscosity study

Ultrasonic attenuation spectra, the shear viscosity, and the mutual diffusion coefficient of the n-pentanol-nitromethane mixture of critical composition have been measured at different temperatures near the critical temperature. The noncritical background contribution, proportional to frequency, to the acoustical attenuation-per-wavelength spectra has been determined and subtracted from the total attenuation to yield the critical contribution. When plotted versus the reduced frequency, with the relaxation rate of order-parameter fluctuations from the shear viscosity and diffusion coefficient measurements, the critical part in the sonic attenuation coefficient displays a scaling function which nicely fits to the data for the critical system 3-methylpentane-nitromethane and also to the empirical scaling function of the Bhattacharjee-Ferrell dynamic scaling theory. The scaled half-attenuation frequency follows from the experimental data as Omega(1/2)emp= 1.8+/-0.1. The relaxation rate of order-parameter fluctuation shows power-law behavior with the theoretically predicted universal exponent and the extraordinary high amplitude Gammao= (187+/-2) x 10(9) s(-1). The amount of the adiabatic coupling constant /g/= 0.03, as estimated from the amplitude of the critical contribution to the acoustical spectra, is unusually small.     

Sensitivity of viscosity Arrhenius parameters to polarity of liquids

Several empirical and semi-empirical equations have been proposed in the literature to estimate the liquid viscosity upon temperature. In this context, this paper aims to study the effect of polarity of liquids on the modeling of the viscosity–temperature dependence, considering particularly the Arrhenius type equations. To achieve this purpose, the solvents are classified into three groups: nonpolar, borderline polar and polar solvents. Based on adequate statistical tests, we found that there is strong evidence that the polarity of solvents affects significantly the distribution of the Arrhenius-type equation parameters and consequently the modeling of the viscosity–temperature dependence. Thus, specific estimated values of parameters for each group of liquids are proposed in this paper. In addition, the comparison of the accuracy of approximation with and without classification of liquids, using the Wilcoxon signed-rank test, shows a significant discrepancy of the borderline polar solvents. For that, we suggested in this paper new specific coefficient values of the simplified Arrhenius-type equation for better estimation accuracy. This result is important given that the accuracy in the estimation of the viscosity–temperature dependence may affect considerably the design and the optimization of several industrial processes.     

The inflection point in the pressure dependence of viscosity under high pressure: a comprehensive study of the temperature and pressure dependence of the viscosity of propylene carbonate

The pressure dependence of the prototypical glass-former propylene carbonate has been investigated over a broad range of temperature and pressure that were inaccessible in previous investigations using dielectric spectroscopy. We find that the viscosity measurements validate the scaling relation, eta(T,V)=J(TV gamma), with a scaling parameter gamma close to that found from dielectric relaxation measurements. In the pressure dependence of the viscosity, we observe an inflection point in the log eta versus P response, similar to that found previously for other materials. However, this inflection has never been observed in dielectric relaxation measurements. Using the scaling property above, it is possible to determine the behavior of the dielectric relaxation time in this otherwise inaccessible experimental range and compare it with the viscosity measurements. We find that the behaviors of eta and tau are very similar, and a very good agreement between the function phi P calculated for these two quantities is found. Starting from the validity of the scaling properties, we show that the inflection point in the pressure dependence of the viscosity can be attributed to the convolution of the pressure dependences of the compressibility kappa T and the apparent activation energy at constant volume EV.     

Thermodynamics of viscous flow and elasticity of glass forming liquids in the glass transition range

The elastic moduli of glasses from different chemical systems, including oxide, chalcogenide, oxynitride, and metallic, were investigated through the glass transition (T(g)), typically from 0.4 to 1.3 T(g). These data were used to interpret the temperature sensitivity of the shear viscosity coefficient obtained on the same materials. The relevant Gibbs free activation energy was estimated from the apparent heat of flow by means of the temperature dependence of the shear elastic modulus. The activation entropy associated with the viscous flow was also derived and was found to correlate with the fragile versus strong character of the glass forming liquids. Finally, the physicochemistry of the flow process was described on the basis of the glass network de-structuration which shows up through the temperature dependence of Poisson's ratio, and an expression for the shear viscosity coefficient is proposed which is chiefly based on the high temperature elastic behavior.     

On the density scaling of liquid dynamics

Superpositioning of relaxation data as a function of the product variable TV(γ), where T is temperature, V the specific volume, and γ a material constant, is an experimental fact demonstrated for approximately 100 liquids and polymers. Such scaling behavior would result from the intermolecular potential having the form of an inverse power law (IPL), suggesting that an IPL is a good approximation for certain relaxation properties over the relevant range of intermolecular distances. However, the derivation of the scaling property of an IPL liquid is based on reduced quantities, for example, the reduced relaxation time equal to T(1∕2V - 1∕3) times the actual relaxation time. The difference between scaling using reduced rather than unreduced units is negligible in the supercooled regime; however, at higher temperature the difference can be substantial, accounting for the purported breakdown of the scaling and giving rise to different values of the scaling exponent. Only the γ obtained using reduced quantities can be sensibly related to the intermolecular potential.     

Viscosity of liquid Cu–Sn alloys

We investigated of the kinematic viscosity of liquid Cu–Sn alloys upon heating and subsequent cooling by the method of the oscillating cylinder. For the liquids alloys Cu75Sn25, Cu50Sn50, Cu48Sn52, Cu32Sn68, and Cu17Sn83, the temperature dependencies of the viscosity upon heating deviate from the Arrhenius relation. The temperature dependencies of viscosity show the Arrhenius-like behaviour upon cooling for all investigated alloys. A discrepancy between the temperature dependencies of viscosity obtained upon heating and cooling arised. We built the concentration dependences of the kinematic viscosity of liquid Cu–Sn alloys upon cooling. The increase of the values of viscosity and activation energy of viscous flow in the concentration range corresponding to the existence of intermetallic compounds Cu₃Sn in the solid state was observed. These results were qualitatively interpreted using the concept of microheterogeneities of liquid alloys.     

Novel halogen-free chelated orthoborate-phosphonium ionic liquids: synthesis and tribophysical properties

We report on the synthesis, characterisation, and physical and tribological properties of halogen-free ionic liquids based on various chelated orthoborate anions with different phosphonium cations, both without halogen atoms in their structure. Important physical properties of the ILs including glass transition temperatures, density, viscosity and ionic conductivity were measured and are reported here. All of these new halogen-free orthoborate ionic liquids (hf-BILs) are hydrophobic and hydrolytically stable liquids at room temperature. As lubricants, these hf-BILs exhibit considerably better antiwear and friction reducing properties under boundary lubrication conditions for steel-aluminium contacts as compared with fully formulated (15W-50 grade) engine oil. Being halogen free these hf-BILs offer a more environmentally benign alternative to ILs being currently developed for lubricant applications.