Kajetan Koperwas

Scaling of volumetric data in model systems based on the Lennard-Jones potential.

The crucial problem for better understanding the nature of glass transition and related relaxation phenomena is to find proper interrelations between the molecular dynamics and thermodynamics of viscous systems. To make progress towards this goal the recently observed density scaling of viscous liquid dynamics has been very intensively and successfully studied in the past few years. However, previous attempts at related scaling of volumetric data yielded results inconsistent with those found from the density scaling of molecular dynamics. In this paper, we show that volumetric data obtained from simulations in simple molecular models based on the Lennard-Jones (LJ) potential, such as the Kob-Andersen binary LJ liquid, its repulsive inverse power-law version, and the Lewis-Wahnström o-terphenyl model, can be scaled by using the same value of the exponent, which scales dynamic quantities and is directly related to the exponent of the repulsive inverse power law that underlies short-range approximations of the LJ potential.

Activation volume in the density scaling regime: Equation of state and its test by using experimental and simulation data

In this paper, a formalism for the activation volume of glass forming materials is suggested. An isothermal equation of state for the activation volume is formulated, which is extended to a generalized equation of state that describes the activation volume as a function of temperature and pressure. Both the equations of state are very successfully validated by using experimental and simulation data collected for supercooled Kob-Andersen binary Lennard-Jones liquid and materials from various material groups such as van der Waals liquids, polymers, protic ionic liquids, and strongly hydrogen bonded liquids. Some predictions based on these equations of state for the activation volume are also very satisfactorily verified in case of each considered system, especially a kind of the activation volume scaling with the scaling exponent that also constitutes the slope of the expected linear pressure dependence of the isothermal bulk modulus for the activation volume is confirmed. The until recently unexpected negative value of the slope are explained in case of the systems that obey the thermodynamic scaling law at least to a good approximation.

Equation of state in the generalized density scaling regime studied from ambient to ultra-high pressure conditions

In this paper, based on the effective intermolecular potential with well separated density and configuration contributions and the definition of the isothermal bulk modulus, we derive two similar equations of state dedicated to describe volumetric data of supercooled liquids studied in the extremely wide pressure range related to the density range, which is extremely wide in comparison with the experimental range reached so far in pressure-volume-temperature measurements of glass-forming liquids. Both the equations comply with the generalized density scaling law of molecular dynamics versus h(ρ)/T at different densities ρ and temperatures T, where the scaling exponent can be in general only a density function γ(ρ) = d ln h/d ln ρ as recently argued by the theory of isomorphs. We successfully verify these equations of state by using data obtained from molecular dynamics simulations of the Kob-Andersen binary Lennard-Jones liquid. As a very important result, we find that the one-parameter density function h(ρ) analytically formulated in the case of this prototypical model of supercooled liquid, which implies the one-parameter density function γ(ρ), is able to scale the structural relaxation times with the value of this function parameter determined by fitting the volumetric simulation data to the equations of state. We also show that these equations of state properly describe the pressure dependences of the isothermal bulk modulus and the configurational isothermal bulk modulus in the extremely wide pressure range investigated by the computer simulations. Moreover, we discuss the possible forms of the density functions h(ρ) and γ(ρ) for real glass formers, which are suggested to be different from those valid for the model of supercooled liquid based on the Lennard-Jones intermolecular potential.

Thermodynamic consequences of the kinetic nature of the glass transition.

In this paper, we consider the glass transition as a kinetic process and establish one universal equation for the pressure coefficient of the glass transition temperature, dTg/dp, which is a thermodynamic characteristic of this process. Our findings challenge the common previous expectations concerning key characteristics of the transformation from the liquid to the glassy state, because it suggests that without employing an additional condition, met in the glass transition, derivation of the two independent equations for dTg/dp is not possible. Hence, the relation among the thermodynamic coefficients, which could be equivalent to the well-known Prigogine-Defay ratio for the process under consideration, cannot be obtained. Besides, by comparing the predictions of our universal equation for dTg/dp and Ehrenfest equations, we find the aforementioned supplementary restriction, which must be met to use the Prigogine-Defay ratio for the glass transition.

Spatially Heterogeneous Dynamics in the Density Scaling Regime: Time and Length Scales of Molecular Dynamics near the Glass Transition.

A fundamental problem of glass transition physics is to find a proper relation between length and time scales of molecular dynamics near the glass transition. Until now, this relation has been usually expected as a single variable function, for instance, as a consequence of the suggested direct relation between the structural relaxation time τ and the correlation volume defined by the maximum of the four-point correlation function χ4max. Based on high pressure data analyses, we show that it is not the case, because χ4max evaluated from its estimate based on the enthalpy fluctuations cannot be, in general, a single variable function of τ. For a wide class of real and model supercooled liquids, the molecular dynamics of which obeys a density scaling law at least to a good approximation, we argue that the important relation between the length and time scales that characterize molecular motions near the glass transition is controlled by a density factor, the exponent of which is a measure of the observed deco...Recent analyses of high pressure measurement data suggest that the degree of the dynamic heterogeneity ({\chi}_4)^max cannot be in general a single variable function of the structural relaxation time {\tau}. For a wide class of real and model supercooled liquids, the molecular dynamics of which obeys a density scaling law at least to a good approximation, we argue that the important relation between the length and time scales that characterize molecular motions near the glass transition is controlled by a density factor. If a power law density scaling is valid for both the structural relaxation times and the degrees of the dynamic heterogeneity we find that the factor is a density power, the exponent of which is a measure of the observed decoupling between {\tau} and ({\chi}_4)^max. Then, the measure can be quantified by a difference between the power law density scaling exponents, which are usually different for {\tau} and ({\chi}_4)^max.

Glass-Forming Tendency of Molecular Liquids and the Strength of the Intermolecular Attractions

When we cool down a liquid below the melting temperature, it can either crystallize or become supercooled, and then form a disordered solid called glass. Understanding what makes a liquid to crystallize readily in one case and form a stable glass in another is a fundamental problem in science and technology. Here we show that the crystallization/glass-forming tendencies of the molecular liquids might be correlated with the strength of the intermolecular attractions, as determined from the combined experimental and computer simulation studies. We use van der Waals bonded propylene carbonate and its less polar structural analog 3-methyl-cyclopentanone to show that the enhancement of the dipole-dipole forces brings about the better glass-forming ability of the sample when cooling from the melt. Our finding was rationalized by the mismatch between the optimal temperature range for the nucleation and crystal growth, as obtained for a modeled Lennard-Jones system with explicitly enhanced or weakened attractive part of the intermolecular 6-12 potential.

Effects of dynamic heterogeneity and density scaling of molecular dynamics on the relationship among thermodynamic coefficients at the glass transition.

In this paper, we define and experimentally verify thermodynamic characteristics of the liquid-glass transition, taking into account a kinetic origin of the process. Using the density scaling law and the four-point measure of the dynamic heterogeneity of molecular dynamics of glass forming liquids, we investigate contributions of enthalpy, temperature, and density fluctuations to spatially heterogeneous molecular dynamics at the liquid-glass transition, finding an equation for the pressure coefficient of the glass transition temperature, dTg/dp. This equation combined with our previous formula for dTg/dp, derived solely from the density scaling criterion, implies a relationship among thermodynamic coefficients at Tg. Since this relationship and both the equations for dTg/dp are very well validated using experimental data at Tg, they are promising alternatives to the classical Prigogine-Defay ratio and both the Ehrenfest equations in case of the liquid-glass transition.