Abstract

The universal scaling concept is applied to the low-temperature range of any liquid states and substances located between the melting () and normal boiling () points far away from the critical region. The physical reason to develop such approach is the revealed collapse of all low-temperature isotherms onto the single universal one argued by the model of fluctuational thermodynamics (FT) proposed recently by author. The pressure reduced by the molecular parameters of the effective short-range Lennard-Jones (LJ) potential depends here only on the reduced density. To demonstrate the extraordinary predictive abilities of the developed low-temperature scaling model it has been applied to the prediction of equilibrium and transport (kinetic and dynamic viscosity, self-diffusion, and thermal conductivity) properties not only for molecular liquids but also for molten organic salts termed ionic liquids (ILs). The best argument in favor of the proposed methodology is the appropriate consistency with the scarce experiments prediction of transport coefficients for ILs on the base of universal scaling function constructed for the simplest LJ-like liquid argon. The only input data of any substance for prediction are the linear approximations of -dependent density and isobaric heat capacity taken from the standard measurements at atmospheric pressure.

1. Introduction

The concept of scaling description at low temperatures based on the chosen effective form of a pair interparticle potential is not completely novel. In particular, Hoover et al. [1] used the oversimplified (even in this specific region of a liquid state) so-called soft-sphere (ss) purely repulsive functionwith constant ss-parameters and exponent ranging from about six for metals to about twelve for rare gases . Its limitative value for a primitive hard-sphere (hs) model leads the scaling variable introduced by the above authors ,to be only -dependent. Another notable feature here is the implied PCS (principle of corresponding states) nature of such ss-scaling. Indeed one has to fix the exponent in (1) and (2) to select the group of similar substances.

Nevertheless the useful consequence of scaling transformation has been noticed by Hoover et al. [1] because any equilibrium or transport reduced properties (e.g., dynamic viscosity or self-diffusion , mass of particle) of different ss-systems with fixed exponent should coincide at the different and being expressed in terms of -variable from (2). Of course, the additional serious restriction of ss-model (common with hs-model) is an absence of gas- (vapor-) liquid transition at any value of repulsive exponent.

It was shown recently [2, 3] in the framework of FT-model [46] that the aforementioned oversimplifications of ss- and/or hs-model are not necessary if one considers the more realistic LJ-type potential with the effective -dependent parameters :The advantage of such approach in comparison with the conventional estimates of their constant LJ-counterparts , (mainly, from the experimental second-virial-coefficient’s data) follows from the exact -dependence derived by FT-model for any subcritical temperature:where , , , and is the -dependent parameter of excluded volume.

The determination of both is not trivial task for ILs [7] in which the vapor-pressure curve is inaccessible by the standard lab. vacuum (~10−3, Pa) experiment. However, just this restriction leads to the significant simplification of (4) at low temperatures as is shown in Figure 1 for the well-studied simple and complex substances.

One may note that in parallel with the well-known condition for compressibility factor measured along the -branch of CXC: its value along the -branch becomes here negligible or close to zero . It leads to the interconsistent low-temperature values of LJ-parameters:This is main assumption used in the present work to formulate the phenomenology of scaling transport coefficients and the respective equilibrium properties following from the generalized reduced form of LJ-EOS discussed below. It is interesting to compare the particular value of ss-scaling variable from (2) at hs-limit : with that from LJ-model in (5b). For the former value is irrelevant while for the latter it determines by (5a) the well-depth energy .

2. Theoretical Basis of Global Asymmetry between the Nonordered (Gas) and Condensed (Liquid) States

FT-model is described in detail elsewhere [2, 6]. For the convenience of readers its main concepts and results are discussed, in brief, below. In the field of equilibrium properties the essential feature is the well-founded rejection from the vdW-concept of a unified (i.e., common for both coexisting - and -phases) EOS. This crucial step leads to the different alternative forms of local EOS for the above phases. They are based, however, on the same mathematical form accepted for an effective pair potential. The choice here of LJ-form requires, for information, at least, the superficial comparison of FT-model’s results with the enormous number of simulated LJ-data described [8] just by the unified, purely empirical (so-called, fundamental) EOS for free energy. Despite the widespread belief of its universal applicability the arguments in favor of global indubitable asymmetry between - and -states [2] are represented below in terms of respective local thermal EOSs. Another theoretical consequence of the discussed duality is the necessity to modify significantly the description of transport coefficients (see Section 3) without the permanent appeals to the well-established PCS formalism of a diluted nonordered -phase and its properties [9, 10].

Striking asymmetry between the exact thermal EOSs of LJ-fluid may be easily revealed by the FT-model’s approach. Thermal form of EOS has the certain restrictions in the representation of caloric properties and, in particular, it needs the knowledge of so-called ideal gas heat capacities if the concept of a unified EOS is used. However, opposite to the conventional methodology of a fundamental (caloric) EOS, its exact form follows immediately for the given PCS-type of LJ-fluid (with the unspecified here integers) from the application of virial theorem [20]. One should note for that the evident formal property of underlying LJ-potential from (3). It is a sum of repulsive and attractive contributions with the different degrees of mathematical homogeneity ( and ). This observation leads straightforwardly to the equivalent forms of thermal EOS for a condensed LJ-fluid ():The important consequence of the latter is the irrelevant role of reduced along Boyle’s line (so-called Zeno-line) where the reduced density becomes constant in the whole range of -states:This result can be consistent with the known and empirically confirmed linearity of Zeno-line expressed, for example, in terms of dimensional variables for the vdW-EOS:only if the reduced LJ-density implies the inverse -dependence of cubic LJ-diameter while the vdW-coefficient is linear in .

At low temperatures where the specific conditions are accepted by (5a) and by (5b) all LJ-isochores collapse to Boyle’s line in accordance with (8). The scaling -form of reduced pressure follows immediately from (7):which collapses here at () to Boyle’s line too. Such coincidence of isobars with isochores is a sign of singular incompressible behavior of liquid in the low-temperature region. It is noteworthy that Maslov [21, 22] has revealed recently the similar phenomenon in the low-temperature, so-called soft LJ-fluid’s states located below Zeno-line in the -plane. He termed it “the jamming effect.”

The described specific features of FT-model’s LJ-liquid in the low-temperature region legalizes the simplest linear choice of input density data (Section 3) measured for any substances at atmospheric pressure and known even for ILs:where “cold” density and “hot” temperature are the substance-dependent parameters (see, for comparison, vdW-equality (9)). However, the two other input -dependencies chosen below at fixed are less motivated and linear ones only for simplicity:where is the isobaric heat capacity of liquid (not that for ideal gas as it is in an unified EOS) and is the adiabatic speed of sound in a liquid too.

To be convinced of the global -asymmetry [2, 3] existing for LJ-fluid one should consider the known virial expansion [23] for nonordered gas truncated after the second virial contribution :Its comparison to (7) leads to the physically well-founded conclusion about the dominance of in -phase but the dominance of inverse reduced temperature in -phase. Moreover, at Boyle’s condition the reduced density is irrelevant while the reduced temperature becomes constant:Such exact estimate can be quite useful at the development of non-mf (mean-field) LJ-theory [2] but in the present work has only passing interest.

3. Continual Phenomenological Model of Transport Coefficients in the Low-Temperature Liquid

The use of such standard notions [23], well-defined in the diluted -phase as the statistical mean free-path of hs-particle between two collisions and its mean transport velocity , becomes less informative due to the global -asymmetry described in Section 2. It may be rejected in favor of the respective continual scales and for -phase determined accurately by its thermodynamic equilibrium parameters , ( is the transport scale (ability) of continuum):Hoover et al. [1] have restricted the proposed scaling ss-description of transport coefficients only by the mechanical transport of mass and momentum by self-diffusion and kinematic velocity represented in terms of -variable from (2):It is natural to add here the temperature-conductivity expressed in terms of thermal conductivity:The certain similarity of scaling functions in (16b) and (17) with the known [9, 10] (inverse here) definitions of the Reynolds’ () and Pecle’s () numbers for a convective flow can be noticed taking into account the choice of characteristic thermodynamic scales by (15).

There are two interconnected kinetic and field scales of velocity and which both are necessary to reduce the momentum flux of pressure:Former is namely the compressibility factor discussed in Section 2 while the latter is the generalized Euler’s number here ( is the adiabatic compressibility). One has to take into account that the spontaneously arisen in homogenous -continuum local Gaussian (small and independent) fluctuations are diffusible with the speed of sound . Thus the transport process of any type depends simultaneously on two (molecular and wave) physical mechanisms [24].

To the best of my knowledge, the contribution of wave transport by the adiabatic fluctuations of pressure in the transport properties of liquids was ignored until now [21, 22]. Its physical nature is different from the molecular-based mechanism and based on the alternative concept of continuum. Hence, two independent scales of fast and slow transport should be taken into account for any fluctuational-relaxation processes. Of course, it is hardly solvable problem to separate them and to determine accurately their relative contributions in the transport coefficients. In principle, they may be either comparable or essentially asymmetric (Section 5) for different liquids.

To make a start of consideration it was assumed below that the equal fractions of both contributions are appropriate for the thermal conductivity and the dynamic viscosity :where subscripts of isobar and of isoentrop specify, respectively, the above types of transport. Then such concept of superposition was realized in the explicit phenomenological forms by the acceptance of two adjustable (but, presumably, universal for all substances) constant coefficients and :where the generalized Prandtl’s number here is essential to estimate the transport ability of a liquid, in total.

It is possibly the most remarkable consequence of the proposed correlations here that the universal coefficients and for a variety of quite different, simple, and complex molecular and/or ionic compounds have been reasonably estimated, exclusively, from the tabular data [11] for the simplest LJ-like argon. Besides, here the chosen on the ad hoc basis correlations are based in (20)–(22) on the semiempirical PCS-functions recommended for a liquid state [9, 10]. The explicit purely empirical form of two scaling functions from (16a), (16b), and (17),is based exclusively on the argon’s data [11] too. In other words, the main aim was the verification of the low-temperature scaling concept but not the best fit of experimental data (see below) on the transport coefficients. More accurately, the above-supposed correlations (20)–(22) take a cue from their successful use in prediction of and , separately, within a group of thermodynamically similar substances. To go beyond the PCS-framework [9, 10] of molecular similarity, the more restrictive concept of scaling (which is applicable, supposedly, to the different PCS-groups) has been used. Hence, the reasonable agreement of the purely predictable coefficients below (Sections 4 and 5) (without any PCS-reference to values and at fixed ) with the experimental data , should bear testimony to general applicability of scaling methodology far away from a critical point. Thus, one may, of course, improve the chosen correlations (20)–(22) here, to fit better the universal scaling parameters , within the possibly modified equations (23a) and (23b). However, the dual molecular-wave structure of all aforementioned universal functions should be the same. The generalized Mach number () following from the ratio of in (18a) and (18b), is independent on the equilibrium parameters , and has the different degrees (linear and quadratic) in the proposed equation (20) for transport of heat and (21) for transport of momentum, respectively.

4. Problem of Input Low-Temperature Data

To avoid any subjectivity of PCS-estimates one has to use the common approximation scheme at the representation of the selected input (i.e., reference) data on each compound. Two individual coefficients of linear equations (11)-(12) determined for each respective property , , from the experimental data [1119] are reported in Table 1 for Ar, N2, C2H4, and C2H6 and the set of seven ILs. For the latter the working temperature range [290, 370 K] of a stable liquid without its upper limit (it is unknown, as a rule) was considered. The possible slight curvature of measured data often leads to the speculations about the sign of thermal expansivity [25]. For (11) it is fixed and positive, at least, at atmospheric pressure:This simple hyperbola gives, however, the good correspondence with the above experimental data. On the other hand, a variety of -data for ILs has been measured in the different ranges: (—higher, —lower), with the different levels of accuracy, and some authors [13] have used the quadratic polynomialsto obtain the best fit. The comparison of such approximation with the different variants of a linear one in (11) reported in Table 2 demonstrates the relatively steady behavior of the “cold” density for both popular ILs: [C4mim][BF4] and [C4mim][PF6].

The influence of discussed approximations on the “hot” temperature from (11) is more pronounced due to the rather large differences in the negative coefficient of (26). The coefficients , reported in Table 3 and calculated without any fit from - and -points provide the maximum deviations in the range of .

Since the ultimate goal of developed here methodology is the scaling prediction of transport coefficients, the search for the best fit of input data is not the task of current work.

Most likely such fit cannot affect essentially the predicted values of self-diffusivity, thermal conductivity, and viscosity of liquids if the chosen scaling form of their description is inherently consistent. Nevertheless, to determine the dependences of input data on the scaling variable from (5b) one needs the reliable estimate of excluded low-temperature volume (see Table 3).

Its value for any ILs can be reasonably estimated [6, 7] just from the long-range extrapolated “cold” volume . The choice of this work emphasizes that the “cold” volume for the variety of molecular and ionic liquids [7] should be always slightly less than the respective mf-estimate of following from the unified vdW-EOS (the mf-estimate is impossible for ILs):where coefficient 1.1 has [6, 7] certain physical meaning.

Indeed, the important, nearly linear correlations of for the wide groups of compounds have been revealed as is shown in Figure 2 and represented in Table 3 for 39 compounds.

There are three (apparent to the naked eye) groups of thermodynamically similar, high-molecular-mass substances: (1) heavy normal -alkanes CnH2(n+1) with ; (2) “friable” ILs (mainly, pyrrolidinium- and ammonium-based ones which can be extrapolated to the low-molecular-mass values for Ar, N2, Li, and H2O as well as above -alkanes to the low-molecular-mass C2H6 and C2H4); (3) “compact” ILs (only imidazolium-based ones). The terms of friable and compact ILs used here are related exclusively to the interrelation between the “cold” volumes at the equal large molecular masses g/mol. The average percent deviations of recommended values (Table 3) from their linear correlations shown in Figure 2 are as follows: for 1st group: 0.22%, for 2nd group: 3.96%, and for 3rd group: 2.36%.

All input linear functions of (11)-(12) expressed in terms of ,can be now directly used for predictions by (20)–(22). In spite of the dramatic differences between values from Table 3 the proposed scaling methodology has to be applicable to all different PCS-types of compounds (Section 5).

5. Comparison of Predicted Transport Coefficients with Experimental Ones

5.1. Thermal Conductivity

The spectacular confirmation of universality, which demonstrates the proposed scaling approach, is represented for thermal conductivity in Figure 3. In contrast with the standard PCS approach [9, 10] any input values of this quantity itself taken at the certain fixed temperature have not been used. Besides, only low-temperature data for argon [11] measured within the quite narrow range of temperature [84, 88 K] are necessary to predict, for example, thermal conductivity of ILs and other compounds in the much wider ranges of temperature shown in Figure 3.

Unfortunately, the experimental data on thermal conductivity (and viscosity) are often scarce and scattered (white symbols) even for the simple molecular liquids. All measured values for ILs (black symbols) have been taken from the work of Nieto de Castro et al. [25]. I have not found in literature the necessary input data on speed of sound for [pmmim][Tf2N] and used those (Table 2) for the structurally close [bmim][Tf2N] with the same molecular mass and common anion. The predicted data in Figure 3 are reasonably compatible with the available literature data on and the maximum deviations are in the range from −6 to 9%.

It is also important and interesting to calculate the relative percentage contribution of molecular-based mechanism in the overall transport of heat and momentum by equations following from (20) and (21):where coefficients and are adjusted to the argon data. Such calculations are represented in Table 4.

One may note that the contribution of molecular mechanism gives roughly one-third in total for simple molecular substances (Ar, N2, C2H4) but becomes more than two-thirds for ILs ([bmim][BF4], [bmim][PF6], [pmmim][Tf2N]). The reverse interrelation of contributions is observable for the transport of momentum (see below). It is more than two-thirds for simple molecular liquids but becomes relatively small for ILs especially near their melting temperatures. Moreover, the contribution of molecular mechanism in for all pyridinium-based (friable) ILs (not represented in Table 4) is less than 3% (i.e., negligible) even for the highest temperature of the measured range. Hence, the natural conclusion is that the transport of momentum in such ILs has the dominantly wave mechanism.

5.2. Viscosity

Taking into account the results of Table 4 it is especially interesting to discuss the scaling prediction of , namely, for the pyridinium-based ILs (Figures 5 and 6) and to compare it with the results for simple molecular liquids (Figure 4).

Again the spectacular predictive abilities of the scaling approach as well as its universality are evident. All measured -data represented by black symbols in Figures 5 and 6 for the pairs of ILs with the common anions ([EtSO4] and [MeSO4], resp.) have been taken from the work of Gonzále et al. [19]. This work contains also the necessary input data on density and speed of sound . Again the reasonable (<30%) overall level of accuracy was demonstrated, though for some ILs ([EtMePyr][MeSO4]) and ([EtMePyr][EtSO4]) the deviations from experiment [19] were large and systematic while for the others ([BuEtPyr][EtSO4] and [BuMePyr][MeSO4]) were close to the experimental uncertainties.

At the same time, our attempt to predict the viscosity for two compact imidazolium-based ILs ([bmim][BF4] and [bmim][PF6]) with the common cation was not so successful. The use of input data on density and speed of sound reported in the work [13] predicts by (21) the significantly underestimated quantities of (located between the limits 0,006 and 0,001 Pa·s). Such predicted values are larger than those for molecular liquids (located (see Figure 4) for C2H4 and C2H6 between the limits of 0,0015 and 0,0002 Pa·s and for Ar and N2 between the limits of 0,0003 and 0,00015 Pa·s) but much lower than ones reported by different authors for [bmim][PF6] (e.g., 0,173; 0,2179; 0,272; 0,450 Pa·s [2628]). Of course, the accuracy of the above scattered experimental data must be improved because here the presumable sharp jump between the average viscosities of friable and compact ILs may be the illuminating factor at the choice of IL for practical usage. High sensitivity of results to the input value has been confirmed, namely, for the viscosity (21), while the thermal conductivity (20) remains relatively insensitive to its value.

5.3. Prandtl’s Number

The remarkable predictive abilities of the scaling approach are again confirmed by Figure 7(a) for molecular liquids and, partially, for the imidazolium-based ILs by Figure 7(b). Taking into account the relative failure of proposed methodology at the prediction of viscosity for [bmim][BF4] and [bmim][PF6], it is interesting to test its ratio at the respective, reasonably predicted value for them as is shown in Figure 3. The result of such test is rather unwonted. The scaling approach predicts the absolutely realistic values for both [bmim][BF4] and [bmim][PF6] (compatible with that (experimental) for [pmmim][Tf2N]) shown in Figure 7(b). The only reason of such observation follows from the definition of -number by (22). To compensate the presumably underestimated dynamic viscosity one needs the significantly overestimated values for the prediction of realistic . This trend is confirmed, partially, in the measurements of isobaric heat capacity of ILs by different authors [29]. The comparison of for molecular liquids (Figure 7(a)) with that for ILs leads to the simple practical conclusion. The latter is preferable, in total, as the thermal storage media while the imidazolium-based ILs can be, additionally, useful as the potential heat-transfer liquids too.

6. Conclusions

The predictive abilities and universality of the proposed scaling model of low-temperature transport coefficients are obvious. This model is testable by the MD-simulations based on the thermodynamically predicted LJ-parameters and here. The derived exact form of EOS for LJ-liquid leads to the potential possibility to construct the crossover formalism connecting critical () and normal boiling () regions by the FT-model’s concepts and results.

Conflict of Interests

The author declares that there is no conflict of interests regarding the publication of this paper.