Non‐asymptotic critical behavior of the transport properties of fluids

A scaling investigation of pattern in the spread of COVID-19: universality in real data and a predictive analytical description

We analyse the spread of COVID-19, a disease caused by a novel coronavirus, in various countries by proposing a model that exploits the scaling and other important concepts of statistical physics. Quite expectedly, for each of the considered countries, we observe that the spread at early times occurs exponentially fast. We show how the countries can be classified into groups, like universality classes in the literature of phase transitions, based on the rates of infections during late times. This method brings a new angle to the understanding of disease spread and is useful in obtaining a country-wise comparative picture of the effectiveness of lockdown-like social measures. Strong similarity, during both natural and lockdown periods, emerges in the spreads within countries having varying geographical locations, climatic conditions, population densities and economic parameters. We derive accurate mathematical forms for the corresponding scaling functions and show how the model can be used as a predictive tool, with instruction even for future waves, and, thus, as a guide for optimizing social measures and medical facilities. The model is expected to be of general relevance in the studies of epidemics.

Correlations for the Viscosity and Thermal Conductivity of Ethyl Fluoride (R161)

This paper presents new wide-ranging correlations for the viscosity and thermal conductivity of ethyl fluoride (R161) based on critically evaluated experimental data. The correlations are designed to be used with a recently published equation of state that is valid from 130 K to 450 K, at pressures up to 100 MPa. The estimated uncertainty at a 95% confidence level is 2% for the viscosity of low-density gas (pressures below 0.5 MPa), and 3% for the viscosity of the liquid over the temperature range from 243 K to 363 K at pressures up to 30 MPa. The estimated uncertainty is 3% for the thermal conductivity of the low-density gas, and 3% for the liquid over the temperature range from 234 K to 374 K at pressures up to 20 MPa. Both correlations may be used over the full range of the equation of state, but the uncertainties will be larger, especially in the critical region.

Finite-size scaling study of dynamic critical phenomena in a vapor-liquid transition

Via a combination of molecular dynamics (MD) simulations and finite-size scaling (FSS) analysis, we study dynamic critical phenomena for the vapor-liquid transition in a three dimensional Lennard-Jones system. The phase behavior of the model has been obtained via the Monte Carlo simulations. The transport properties, viz., the bulk viscosity and the thermal conductivity, are calculated via the Green-Kubo relations, by taking inputs from the MD simulations in the microcanonical ensemble. The critical singularities of these quantities are estimated via the FSS method. The results thus obtained are in nice agreement with the predictions of the dynamic renormalization group and mode-coupling theories.

Stationary and transient Soret separation in a binary mixture with a consolute critical point

The stationary and transient Soret separation in a binary mixture with a consolute critical point is studied theoretically. The mixture is placed between two parallel plates kept at different temperatures. A polymer blend is used as a model system. Analytical solutions are constructed to describe the stationary separation in a binary mixture with variable Soret coefficient. The latter strongly depends on temperature and concentration and enhances near a consolute critical point due to reduced diffusion. As a result, a large concentration gradient is observed locally, while much smaller concentration variations are found in the rest of the layer. It is shown that complete separation can be obtained by applying a small temperature difference first, waiting for the establishment of stationary state, and then increasing this difference again. In this case, the critical temperature lies between hot and cold wall temperatures, while the mixture still remains in the one-phase region. When the initial (mean) temperature or concentration are shifted away from the near-critical values, the separation decreases. The analysis of transient behavior shows that the Soret separation occurs much faster than diffusion to the homogeneous state when the initial concentration is close to the critical one. It happens due to the decrease (increase) of the local relaxation time during the Soret (Diffusion) steps. The transient times of these steps become comparable for small temperature differences or off-critical initial concentrations. An unusual (non-exponential) separation dynamics is observed when the separation starts in the off-critical domain, and then enhances greatly when the system enters into the near-critical region. It is also found that the transient time decreases with increasing the applied temperature difference.

Study of critical dynamics in fluids via molecular dynamics in canonical ensemble

With the objective of understanding the usefulness of thermostats in the study of dynamic critical phenomena in fluids, we present results for transport properties in a binary Lennard-Jones fluid that exhibits liquid-liquid phase transition. Various collective transport properties, calculated from the molecular dynamics (MD) simulations in canonical ensemble, with different thermostats, are compared with those obtained from MD simulations in microcanonical ensemble. It is observed that the Nosé-Hoover and dissipative particle dynamics thermostats are useful for the calculations of mutual diffusivity and shear viscosity. The Nosé-Hoover thermostat, however, as opposed to the latter, appears inadequate for the study of bulk viscosity.

Simulation of transport around the coexistence region of a binary fluid

We use Monte Carlo and molecular dynamics simulations to study phase behavior and transport properties in a symmetric binary fluid where particles interact via Lennard-Jones potential. Our results for the critical behavior of collective transport properties, with particular emphasis on bulk viscosity, is understood via appropriate application of finite-size scaling technique. It appears that the critical enhancements in these quantities are visible far above the critical point. This result is consistent with an earlier report from computer simulations where, however, the authors do not quantify the critical singularity.

Thermal conductivity of supercooled water

The heat capacity of supercooled water, measured down to -37°C, shows an anomalous increase as temperature decreases. The thermal diffusivity, i.e., the ratio of the thermal conductivity and the heat capacity per unit volume, shows a decrease. These anomalies may be associated with a hypothesized liquid-liquid critical point in supercooled water below the line of homogeneous nucleation. However, while the thermal conductivity is known to diverge at the vapor-liquid critical point due to critical density fluctuations, the thermal conductivity of supercooled water, calculated as the product of thermal diffusivity and heat capacity, does not show any sign of such an anomaly. We have used mode-coupling theory to investigate the possible effect of critical fluctuations on the thermal conductivity of supercooled water and found that indeed any critical thermal-conductivity enhancement would be too small to be measurable at experimentally accessible temperatures. Moreover, the behavior of thermal conductivity can be explained by the observed anomalies of the thermodynamic properties. In particular, we show that thermal conductivity should go through a minimum when temperature is decreased, as Kumar and Stanley observed in the TIP5P model of water. We discuss physical reasons for the striking difference between the behavior of thermal conductivity in water near the vapor-liquid and liquid-liquid critical points.

Dynamic structure factor of density fluctuations from direct imaging very near (both above and below) the critical point of SF(6)

Large density fluctuations were observed by illuminating a cylindrical cell filled with sulfur hexafluoride (SF(6)), very near its liquid-gas critical point (|T-T(c)|< 300 μK) and recorded using a microscope with 3 μm spatial resolution. Using a dynamic structure factor algorithm, we determined from the recorded images the structure factor (SF), which measures the spatial distribution of fluctuations at different moments, and the correlation time of fluctuations. This method authorizes local measurements in contrast to the classical scattering techniques that average fluctuations over the illuminating beam. We found that during the very early stages of phase separation the SF scales with the wave vector q according to the Lorentzian q(-2), which shows that the liquid and vapor domains are just emerging. The critical wave number, which is related to the characteristic length of fluctuations, steadily decreases over time, supporting a sustained increase in the spatial scale of the fluctuating domains. The scaled evolution of the critical wave number obeys the universal evolution for the interconnected domains at high volume fraction with an apparent power law exponent of -0.35 ± 0.02. We also determined the correlation time of the fluctuations and inferred values for thermal diffusivity coefficient very near the critical point, above and below. The values were used to pinpoint the crossing of T(c) within 13 μK.

Critical dynamics of an isothermal compressible nonideal fluid

A pure fluid at its critical point shows a dramatic slow-down in its dynamics, due to a divergence of the order-parameter susceptibility and the coefficient of heat transport. Under isothermal conditions, however, sound waves provide the only possible relaxation mechanism for order-parameter fluctuations. Here we study the critical dynamics of an isothermal, compressible nonideal fluid via scaling arguments and computer simulations of the corresponding fluctuating hydrodynamics equations. We show that, below a critical dimension of 4, the order-parameter dynamics of an isothermal fluid effectively reduces to "model A," characterized by overdamped sound waves and a divergent bulk viscosity. In contrast, the shear viscosity remains finite above two dimensions. Possible applications of the model are discussed.

Dynamics of critical fluctuations in polymer solutions

Using dynamic light scattering we have investigated the time dependence of fluctuations near the critical point of phase separation in solutions of polystyrene in cyclohexane with polymer molecular weights ranging from 196,000 to 11.4 x 10(6) g mol(-1). At the lowest polymer molecular weight the dynamic correlation function follows a single-exponential decay with a decay rate that can be represented by the mode-coupling theory of critical dynamics but with a mesoscopic viscosity that characterizes the hydrodynamic environment of the polymers in the solution. At all higher polymer molecular weights two distinct dynamic modes are observed, a slow and a fast mode, that originate from a coupling of the critical concentration fluctuations with viscoelastic relaxation of the polymer chain in solutions. This coupling causes an additional slowing down of the fluctuations on top of the well-known critical slowing down expected in the absence of a coupling between the two modes. From an analysis of the time dependence of the experimental dynamic correlation functions in terms of a theory of coupling of dynamic modes we are able to determine the viscoelastic properties of the polymers in the solution. These viscoelastic properties diverge in the theta-point limit of infinite polymer molecular weight.

Static and dynamic critical behavior of a symmetrical binary fluid: A computer simulation

A symmetrical binary, A+B Lennard-Jones mixture is studied by a combination of semi-grand-canonical Monte Carlo (SGMC) and molecular dynamics (MD) methods near a liquid-liquid critical temperature T(c). Choosing equal chemical potentials for the two species, the SGMC switches identities (A-->B-->A) to generate well-equilibrated configurations of the system on the coexistence curve for T<T(c) and at the critical concentration, x(c)=12, for T>T(c). A finite-size scaling analysis of the concentration susceptibility above T(c) and of the order parameter below T(c) is performed, varying the number of particles from N=400 to 12 800. The data are fully compatible with the expected critical exponents of the three-dimensional Ising universality class. The equilibrium configurations from the SGMC runs are used as initial states for microcanonical MD runs, from which transport coefficients are extracted. Self-diffusion coefficients are obtained from the Einstein relation, while the interdiffusion coefficient and the shear viscosity are estimated from Green-Kubo expressions. As expected, the self-diffusion constant does not display a detectable critical anomaly. With appropriate finite-size scaling analysis, we show that the simulation data for the shear viscosity and the mutual diffusion constant are quite consistent both with the theoretically predicted behavior, including the critical exponents and amplitudes, and with the most accurate experimental evidence.

Critical dynamics in a binary fluid: simulations and finite-size scaling

We report comprehensive simulations of the critical dynamics of a symmetric binary Lennard-Jones mixture near its consolute point. The self-diffusion coefficient exhibits no detectable anomaly. The data for the shear viscosity and the mutual-diffusion coefficient are fully consistent with the asymptotic power laws and amplitudes predicted by renormalization-group and mode-coupling theories provided finite-size effects and the background contribution to the relevant Onsager coefficient are suitably accounted for. This resolves a controversy raised by recent molecular simulations.

Critical fluctuations of the micellar triethylene glycol monoheptyl ether-water system

Using the equal volume criterion and also the pseudospinodal conception the critical demixing point of the triethylene glycol monoheptyl ether/water system (C7E3H2O) has been determined as Ycrit=0.1 and Tcrit=296.46 K (Y, mass fraction of surfactant). From density measurements the critical micelle concentration (cmc) followed as Ycmc=0.007 at 288.15 K and Ycmc=0.0066 at 298.15 K. The (static) shear viscosity etas and the mutual diffusion coefficient D of the C7E3H2O mixture of critical composition have been evaluated to yield their singular and background parts. From a combined treatment of both quantities the relaxation rate Gamma of order parameter fluctuations has been derived. Gamma follows power law with universal critical exponent and amplitude Gamma0=3.1 x 10(9) s(-1). Broadband ultrasonic spectra of C7E3H2O mixtures exhibit a noncritical relaxation, reflecting the monomer exchange between micelles and the suspending phase, and a critical term due to concentration fluctuations. The former is subject to a relaxation time distribution that broadens when approaching the critical temperature. The latter can be well represented with the aid of the dynamic scaling model by Bhattacharjee and Ferrell (BF) [Phys. Rev. A. 31, 1788 (1985)]. The half-attenuation frequency in the scaling function of the latter model is noticeably smaller (Omega12 (BF) approximately 1) than the theoretically predicted value Omega12 (BF)=2.1. This result has been taken as an indication of a coupling between the fluctuations in the local concentration and the kinetics of micelle formation, in correspondence with the idea of a fluctuation controlled monomer exchange [T. Telgmann and U. Kaatze, Langmuir 18, 3068 (2002)].

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.

Droplet motion with phase change in a temperature gradient

We examine the droplet motion in one-component fluids in a small temperature gradient by solving linearized hydrodynamic equations supplemented with appropriate surface boundary conditions. We show that the velocity field and the temperature around the droplet are strongly influenced by first-order phase transition taking place at the interface. Latent heat released or absorbed at the interface drastically changes the hydrodynamic flow around the droplet. As a result, the temperature becomes almost homogeneous inside the droplet and the Marangoni effect arising from the surface tension gradient is much suppressed. The droplet velocity is also much decelerated.

Dynamics of fluids near the consolute critical point: A molecular-dynamics study of the Widom–Rowlinson mixture

Molecular-dynamics simulations are presented for the dynamic behavior of the Widom-Rowlinson mixture [B. Widom, and J. S. Rowlinson, J. Chem. Phys. 52, 1670 (1970)] at its critical point. This model consists of two components where like species do not interact and unlike species interact via a hard-core potential. Critical exponents are obtained from a finite-size scaling analysis. The self-diffusion coefficient shows no anomalous behavior near the critical point. The shear viscosity and thermal conductivity show no divergent behavior for the system sizes considered, although there is a significant critical enhancement. The mutual diffusion coefficient, D(AB), vanishes as D(AB) approximately xi(-1.26 +/- 0.08), where xi is the correlation length. This is different from the renormalization-group (D(AB) approximately xi(-1.065)) mode coupling theory (D(AB) approximately xi(-1)) predictions. The theories and simulations can be reconciled if we assume that logarithmic corrections to scaling are important.

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)].

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.

Shear viscosity studies above and below the critical consolute point in a nitrobenzene-decane mixture

The shear viscosity has been studied in a nitrobenzene-decane critical mixture above the critical consolute temperature T(C), in the homogeneous phase, and below T(C), in coexisting phases. The form of background viscosity for coexisting phases has been postulated. The same value of the critical exponent straight phi has been obtained in the lower (L), upper (U), and homogeneous (H) phases. The pretransitional amplitudes (A(L,U)) in coexisting phases are approximately the same, whereas A(L,U)/A(H) approximately 0.965. In the homogeneous phase the possibility of the appearance of the quasinematic, field-induced structure of critical fluctuations has been discussed.

Critical light scattering in liquids

We compare theoretical results for the characteristic frequency of the Rayleigh peak calculated in one-loop order within the field theoretical method of the renormalization group theory with experiments and other theoretical results. Our expressions describe the nonasymptotic crossover in temperature, density, and wave vector. In addition we discuss the frequency dependent shear viscosity evaluated within the same model and compare our theoretical results with recent experiments in microgravity.

Frequency-dependent viscosity of xenon near the critical point

We used a novel, overdamped oscillator aboard the Space Shuttle to measure the viscosity eta of xenon near its critical density rho(c) and temperature Tc. In microgravity, useful data were obtained within 0.1 mK of Tc, corresponding to a reduced temperature t=(T-Tc)/Tc=3 x 10(-7). Because they avoid the detrimental effects of gravity at temperatures two decades closer to T(c) than the best ground measurements, the data directly reveal the expected power-law behavior eta proportional, variant t(-nuz(eta)). Here nu is the correlation length exponent, and our result for the viscosity exponent is z(eta)=0.0690+/-0.0006. (All uncertainties are one standard uncertainty.) Our value for z(eta) depends only weakly on the form of the viscosity crossover function, and it agrees with the value 0.067+/-0.002 obtained from a recent two-loop perturbation expansion [H. Hao, R.A. Ferrell, and J.K. Bhattacharjee, (unpublished)]. The measurements spanned the frequency range 2 Hz< or = f < or =12 Hz and revealed viscoelasticity when t < or = 10(-5), further from Tc than predicted. The viscoelasticity's frequency dependence scales as Aftau, where tau is the fluctuation-decay time. The fitted value of the viscoelastic time-scale parameter A is 2.0+/-0.3 times the result of a one-loop perturbation calculation. Near Tc, the xenon's calculated time constant for thermal diffusion exceeded days. Nevertheless, the viscosity results were independent of the xenon's temperature history, indicating that the density was kept near rho(c) by judicious choices of the temperature versus time program. Deliberately bad choices led to large density inhomogeneities. At t>10(-5), the xenon approached equilibrium much faster than expected, suggesting that convection driven by microgravity and by electric fields slowly stirred the sample.

Frequency-dependent shear viscosity, sound velocity, and sound attenuation near the critical point in liquids. III. The shear viscosity

We compare theoretical results for the shear viscosity calculated in one-loop order within the field-theoretical method of the renormalization-group theory with experiments. Our expressions describe the nonasymptotic crossover in both temperature and density, and allow us to consider effects of finite gravitation and finite frequency at which the experiments are performed. In doing so we treat the critical exponent x(eta) of the shear viscosity as an independent parameter, keeping the one-loop value of the Kawasaki amplitude fixed. Within our model we also consider the temperature and density dependence of the thermal diffusion including gravitational effects.