Source author record

Jeppe C. Dyre

Jeppe C. Dyre appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.soft cond-mat.stat-mech cond-mat.mtrl-sci cond-mat.dis-nn physics.chem-ph physics.comp-ph cond-mat.other cond-mat.mes-hall Biological Physics cond-mat physics.ed-ph physics.flu-dyn physics.geo-ph physics.ins-det physics.plasm-ph

Catalog footprint

What is connected

89works

15topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2023arXiv

Active-matter isomorphs in the size-polydisperse Ornstein-Uhlenbeck Lennard-Jones model

This paper studies size-polydisperse Lennard-Jones systems described by active Ornstein-Uhlenbeck particle dynamics. The focus is on the existence of isomorphs (curves of invariant structure and dynamics) in the model's three-dimensional phase diagram. Isomorphs are traced out from a single steady-state configuration by means of the configurational-temperature method. Good invariance of the reduced-unit radial distribution function and the mean-square displacement as a function of time is demonstrated for three uniform-distribution polydispersities, 12%, 23%, and 29%. Comparing to active-matter isomorphs generated by the analytical direct-isomorph-check method, the latter give somewhat poorer invariance of the structure, but better invariance of the dynamics. We conclude that both methods can be used to quickly get an overview of the phase diagram of polydisperse AOUP models involving a potential-energy function obeying the hidden-scale-invariance property required for isomorph theory to apply.

preprint2023arXiv

Active-parameter polydispersity in the 2d ABP Yukawa model

In both experiments and simulations the most commonly studied kind of parameter polydispersity is that of varying particles size. This paper investigates by simulations the effects of introducing polydispersity in other parameters for two-dimensional Active Brownian Particles with Yukawa pair interactions. Polydispersity is studied separately in the translational and rotational diffusion coefficients, as well as in the swim velocity $v_0$. Uniform and binary parameter distributions are considered in both the homogeneous and the motility-induced phase-separation (MIPS) phases. We find only minute changes in structure and dynamics upon the introduction of parameter polydispersity, even for situations involving 50% polydispersity. The reason for this is not clear. An exception is the case of $v_0$ polydispersity for which the average radial distribution function with changing polydispersity shows significant variations in the MIPS phase. Even in this case, however, the dynamics is only modestly affected. As a possible application of our findings, we suggest that a temporary introduction of polydispersity into a single-component active-matter model characterized by a very long equilibration time, i.e., a glass-forming active system, may be used to equilibrate the system efficiently by particle swaps.

preprint2022arXiv

Hidden scale invariance in the Gay-Berne model

This paper presents a numerical study of the Gay-Berne liquid crystal model with parameters corresponding to calamitic (rod-shaped) molecules. The focus is on the isotropic and nematic phases at temperatures above unity. There we find strong correlations between the virial and potential-energy thermal fluctuations, reflecting the hidden-scale invariance symmetry. This implies the existence of isomorphs, which are curves in the thermodynamic phase diagram of approximately invariant physics. We study numerically one isomorph in the isotropic phase and one in the nematic phase. In both cases, good invariance of the dynamics is demonstrated via data for the reduced-unit time-autocorrelation functions of the mean-square displacement, angular velocity, force, torque, and first- and second-order Legendre polynomial orientational order parameters. Deviations from isomorph invariance are observed at short times for the orientational time-autocorrelation functions, which reflects the fact that the moment of inertia is assumed to be constant and thus not isomorph invariant in reduced units. Structural isomorph invariance is demonstrated from data for the radial distribution functions of the particles and their orientations. For comparison, all quantities were also simulated along an isochore of similar temperature variation in which case invariance is not observed. We conclude that the thermodynamic phase diagram of the calamitic Gay-Berne model is essentially one-dimensional in the studied regions as predicted by isomorph theory, a fact that potentially allows for simplifications of future theories and numerical studies.

preprint2022arXiv

Isomorph invariance in the liquid and plastic-crystal phases of asymmetric-dumbbell models

We present a numerical study of the asymmetric dumbbell model consisting of ``molecules'' constructed as two different-sized Lennard-Jones spheres connected by a rigid bond. In terms of the largest (A) particle radius, we report data for the structure and dynamics of the liquid phase for the bond lengths 0.05, 0.1, 0.2, and 0.5, and analogous data for the plastic-crystal phase for the bond lengths 0.05, 0.1, 0.2, and 0.3. Structure is probed by means of the AA, AB, and BB radial distribution functions. Dynamics is probed via the A and B particle mean-square displacement as functions of time and via the rotational time-autocorrelation function. Consistent with the systems' strong virial potential-energy correlations, the structure and dynamics are found to be isomorph invariant to a good approximation in reduced units, while they generally vary considerably along isotherms of the same (20%) density variation. Even the rotational time-autocorrelation function, which due to the constant bond length is not predicted to be isomorph invariant, varies more along isotherms than along isomorphs. Our findings provide the first validation of isomorph-theory predictions for plastic crystals for which isomorph invariance, in fact, is found to apply better than in the liquid phase of asymmetric-dumbbell models.

preprint2022arXiv

Predicting nonlinear physical aging of glasses from equilibrium relaxation via the material time

The noncrystalline glassy state of matter plays a role in virtually all fields of materials science and offers complementary properties to those of the crystalline counterpart. The caveat of the glassy state is that it is out of equilibrium and therefore exhibits physical aging, i.e., material properties change over time. For half a century, the physical aging of glasses has been known to be described well by the material-time concept, although the existence of a material time has never been directly validated. We do this here by successfully predicting the aging of the molecular glass 4-vinyl-1,3-dioxolan-2-one from its linear relaxation behavior. This establishes the defining property of the material time. Via the fluctuation-dissipation theorem, our results imply that physical aging can be predicted from thermal-equilibrium fluctuation data, which is confirmed by computer simulations of a binary liquid mixture.

preprint2022arXiv

Rheological model for the alpha relaxation of glass-forming liquids and its comparison to data for DC704 and DC705

Dynamic shear-modulus data are presented for the two silicone oils DC704 and DC705 for frequencies between 1 mHz and 10 kHz at temperatures covering more than five decades of relaxation-time variation. The data are fitted to the alpha part of a phenomenological model previously shown to describe well the dynamic shear modulus of squalane, which has a large beta process [Hecksher \textit{et al.}, J. Chem. Phys. \textbf{146}, 154504 (2017)]; that model is characterized by additivity of the alpha and beta shear compliance and by a high-frequency decay of the alpha process in proportion to $ω^{-1/2}$ in which $ω$ is the angular frequency. The fits of the alpha part of this model to the DC704 and DC705 data are compared to fits by a Havriliak-Negami type model, the Barlow-Erginsav-Lamb model, and a Cole-Davidson type model. At all temperatures the best fit is obtained by the alpha part of the squalane model. This strengthens the conjecture that so-called $\sqrt{t}$-relaxation, leading to high-frequency decays proportional to $ω^{-1/2}$, is a general characteristic of the alpha relaxation of supercooled liquids [Dyre, Phys. Rev. E {\bf 74}, 021502 (2006); Nielsen \textit{et al.}, J. Chem. Phys. \textbf{130}, 154508 (2009); Pabst \textit{et al.}, J. Phys. Chem. Lett. \textbf{12}, 3685 (2021)].

preprint2022arXiv

Single-parameter aging in the weakly nonlinear limit

Physical aging deals with slow property changes over time caused by molecular rearrangements. This is relevant for non-crystalline materials like polymers and inorganic glasses, both in production and during subsequent use. The Narayanaswamy theory from 1971 describes physical aging - an inherently nonlinear phenomenon - in terms of a linear convolution integral over the so-called material time $ξ$. The resulting "Tool-Narayanaswamy (TN) formalism" is generally recognized to provide an excellent description of physical aging for small, but still highly nonlinear temperature variations. The simplest version of the TN formalism is single-parameter aging according to which the clock rate $dξ/dt$ is an exponential function of the property monitored [T. Hecksher et al., J. Chem. Phys. 142, 241103 (2015)]. For temperature jumps starting from thermal equilibrium, this leads to a first-order differential equation for property monitored, involving a system-specific function. The present paper shows analytically that the solution to this equation to first order in the temperature variation has a universal expression in terms of the zeroth-order solution, $R_0(t)$. Numerical data for a binary Lennard-Jones glass former probing the potential energy confirm that, in the weakly nonlinear limit, the theory predicts aging correctly from $R_0(t)$ (which by the fluctuation-dissipation theorem is the normalized equilibrium potential-energy time-autocorrelation function).

preprint2021arXiv

Generalized hydrodynamics of the Lennard-Jones liquid in view of hidden scale invariance

In recent years lines along which structure and dynamics are invariant to a good approximation, so-called isomorphs, have been identified in the thermodynamic phase diagrams of several model liquids and solids. This paper reports computer simulations of the transverse and longitudinal collective dynamics at different length scales along an isomorph of the Lennard-Jones system. Our findings are compared to corresponding results along an isotherm and an isochore. Confirming the theoretical prediction, the reduced-unit dynamics of the transverse momentum density is invariant to a good approximation along the isomorph at all time and length scales. Likewise, the wave-vector dependent shear-stress autocorrelation function is found to be isomorph invariant. A similar invariance is not seen along the isotherm or the isochore. Using a spatially non-local hydrodynamic model for the transverse momentum-density time-autocorrelation function, the macroscopic shear viscosity and its wave dependence are determined, demonstrating that the shear viscosity is isomorph invariant on all length scales studied. This analysis implies the existence of a novel length scale which characterizes each isomorph. The transverse sound-wave velocity, the Maxwell relaxation time, and the rigidity shear modulus are also isomorph invariant. In contrast, the reduced-unit dynamics of the mass density is not invariant at length scales longer than the inter-particle distance. By fitting to a generalized hydrodynamic model, we extract values for the wave-vector-dependent thermal diffusion coefficient, sound attenuation coefficient, and adiabatic sound velocity. The isomorph variation of these quantities in reduced units at long length scales can be eliminated by scaling with the density-scaling, a fundamental quantity in the isomorph theory framework, an empirical observation that remains to be explained theoretically.

preprint2021arXiv

Isomorphs in nanoconfined liquids

We study in this paper the possible existence of Roskilde-simple liquids and their isomorphs in a rough-wall nanoconfinement. Isomorphs are curves in the thermodynamic phase diagram along which structure and dynamics are invariant in suitable nondimensionalized units. Two model liquids using molecular dynamics computer simulations are considered: the single-component Lennard-Jones (LJ) liquid and the Kob-Andersen binary LJ mixture, both of which in the bulk phases are known to have isomorphs. Nanoconfinement is implemented by adopting a slit-pore geometry with fcc crystalline walls; this implies inhomogenous density profiles both parallel and perpendicular to the confining walls. Despite this fact and consistent with an earlier study [Ingebrigtsen et. al, Phys. Rev. Lett. 111, 235901 (2013)] we find that these nanoconfined liquids have isomorphs to a good approximation. More specifically, we show good scaling of inhomogenous density profiles, mean-square displacements, and higher-order structures probed using the topological cluster classification algorithm along the isomorphs. From this study, we conjecture that in experiments, Roskilde-simple liquids may exhibit isomorphs if confined in a suitable manner, for example with carbon nanotubes. Our study thus provides an alternative framework for understanding nanoconfined liquids.

preprint2021arXiv

Single-parameter aging in a binary Lennard-Jones system

This paper studies physical aging by computer simulations of a 2:1 Kob-Andersen binary Lennard-Jones mixture, a system that is less prone to crystallization than the standard 4:1 composition. Starting from thermal-equilibrium states, the time evolution of the following four quantities is monitored following up and down jumps in temperature: the potential energy, the virial, the average squared force, and the Laplacian of the potential energy. Despite the fact that significantly larger temperature jumps are studied here than in previous experiments, to a good approximation all four quantities conform to the single-parameter-aging scenario derived and validated for small jumps in experiments [Hecksher et al., J. Chem. Phys. 142, 241103 (2015)]. As a further confirmation of single-parameter aging with a common material time for the different quantities monitored, their relaxing parts are found to be almost identical for all temperature jumps.

preprint2021arXiv

Time-scale coupling in hydrogen- and van der Waals-bonded liquids

The coupling behavior of time scales of structural relaxation is investigated on the basis of five different response functions for 1,2,6-hexanetriol, a hydrogen-bonded liquid with a minor secondary contribution, and 2,6,10,15,19,23-hexamethyl-tetracosane (squalane), a van der Waals bonded liquid with a prominent secondary relaxation process. Time scales of structural relaxation are derived as inverse peak frequencies for each investigated response function. For 1,2,6-hexanetriol, the time-scale indices are temperature-independent, while a decoupling of time scales is observed for squalane in accordance with literature. An alternative evaluation approach is made on the squalane data, extracting time scales from the terminal relaxation mode instead of the peak position, and in this case temperature-independent coupling is also found for squalane, despite its strong secondary relaxation contribution. Interestingly, the very same ordering of response-function-specific time scales is observed for these two liquids, which is also consistent with the observation made for simple van der Waals bonded liquids reported previously [Jakobsen \textit{et al.}, J. Chem. Phys. \textbf{136}, 081102 (2012)]. This time-scale ordering is based on the following response functions, from fast to slow dynamics: shear modulus, bulk modulus, dielectric permittivity, longitudinal thermal expansivity coefficient, and longitudinal specific heat. These findings indicate a general relation between the time scales of different response functions and, as inter-molecuar interactions apparently play a subordinate role, suggest a rather generic nature of the process of structural relaxation.

preprint2020arXiv

Solid-like mean-square displacement in glass-forming liquids

It was recently shown that the real part of the frequency-dependent fluidity for several glass-forming liquids of different chemistry conforms to the prediction of the random barrier model (RBM) devised for ac electrical conduction in disordered solids [S. P. Bierwirth \textit{et al.}, Phys. Rev. Lett. {\bf 119}, 248001 (2017)]. Inspired by these results we introduce a crystallization-resistant modification of the Kob-Andersen binary Lennard-Jones mixture for which the results of extensive graphics-processing unit (GPU)-based molecular-dynamics simulations are presented. We find that the low-temperature mean-square displacement is fitted well by the RBM prediction, which involves no shape parameters. This finding highlights the challenge of explaining why a simple model based on hopping of non-interacting particles in a fixed random energy landscape can reproduce the complex and highly cooperative dynamics of glass-forming liquids.

preprint2020arXiv

Solid-liquid coexistence of the noble elements. I. Theory illustrated by the case of argon

The noble elements constitute the simplest group of atoms. At low temperatures or high pressures they freeze into the face-centered cubic (fcc) crystal structure (except helium). We perform molecular dynamics using the recently proposed simplified ab initio atomic (SAAP) potential [Deiters and Sadus, J. Chem. Phys. 150, 134504 (2019)] . This potential is parameterized using data from accurate ab initio quantum mechanical calculations by the coupled-cluster approach on the CCSD(T) level. We compute the fcc freezing lines for Argon and find a great agreement with the experimental values. At low pressures, this agreement is further enhanced by using many-body corrections. Hidden scale invariance of the potential energy function is validated by computing lines of constant excess entropy (configurational adiabats) and shows that mean square displacement and the static structure factor are invariant. These lines (isomorphs) can be generated from simulations at a single state-point by having knowledge of the pair potential. The isomorph theory for the solid-liquid transition is used to accurately predict the shape of the freezing line in the pressure-temperature plane, the shape in the density-temperature plane, the entropy of melting and the Lindemann parameters along the melting line. We finally predict that the body-centered cubic (bcc) crystal is stable at high pressures.

preprint2020arXiv

Solid-liquid coexistence of the noble elements. II. Neon, krypton and xenon

The noble elements constitute the simplest group in the periodic table. At low temperatures or high pressures, the liquid phase solidifies into a face-centered cubic crystal structure (except helium). In the companion paper, we investigated the fcc solid-liquid coexistence of argon in the light of hidden scale invariance. Here we extend the investigation to neon, argon, krypton, and xenon. Computations are done using the SAAP potential, suggested by Deiters and Sadus [J. Chem. Phys 150, 134504 (2019)], derived from accurate ab initio calculations. The systems exhibit hidden scale invariance in the investigated part of the phase diagram, which makes it possible to predict the shape and property variations along the solid-liquid coexistence lines.

preprint2020arXiv

Structure of the Lennard-Jones liquid estimated from a single simulation

Combining the recent Piskulich-Thompson approach [Z. A. Piskulich and W. H. Thompson, {\it J. Chem. Phys.} {\bf 152}, 011102 (2020)] with isomorph theory, from a single simulation, the structure of a single-component Lennard-Jones (LJ) system is obtained at an arbitrary state point in almost the whole liquid region of the temperature-density phase diagram. The LJ system exhibits two temperature range where the van't Hoff's assumption that energetic and entropic forces are temperature independent is valid. A method to evaluate the structure at an arbitrary state point along an isochore from the knowledge of structures at two temperatures on the isochore is also discussed. We argue that, in general, the structure of any hidden scale-invariant system obeying the van't Hoff's assumption in the whole range of temperatures can be determined in the whole liquid region of the phase diagram from only a single simulation.

preprint2016arXiv

"Pseudoisomorphs" in liquids with intramolecular degrees of freedom

Computer simulations show that liquids of molecules with harmonic intramolecular bonds may have "pseudoisomorphic" lines of approximately invariant dynamics in the thermodynamic phase diagram. We demonstrate that these lines can be identified by requiring scale invariance of the inherent-structure reduced-unit low-frequency vibrational spectrum evaluated for a single equilibrium configuration. This rationalizes why excess-entropy scaling, density scaling, and isochronal superposition apply for many liquids with internal degrees of freedom.

preprint2016arXiv

A pair potential that reproduces the shape of isochrones in molecular liquids

Many liquids have curves (isomorphs) in their phase diagrams along which structure, dynamics, and some thermodynamic quantities are invariant in reduced units. A substantial part of their phase diagrams is thus effectively one dimensional. The shape of these isomorphs is described by a material-dependent function of density $h(ρ)$, which for real liquids is well approximated by a power law $ρ^γ$. However, in simulations, a power law is not adequate when density changes are large; typical models such as the Lennard-Jones liquid show that $γ(ρ) \equiv \mathrm{d} \ln h(ρ)/\mathrm{d} \ln ρ$ is a decreasing function of density. This paper presents results from computer simulations using a new pair potential that diverges at a nonzero distance and can be tuned to give a more realistic shape of $γ(ρ)$. Our results indicate that the finite size of molecules is an important factor to take into account when modeling liquids over a large density range.

preprint2016arXiv

Freezing and melting line invariants of the Lennard-Jones system

The invariance of several structural and dynamical properties of the Lennard-Jones (LJ) system along the freezing and melting lines is interpreted in terms of the isomorph theory. First the freezing/melting lines for LJ system are shown to be approximated by isomorphs. Then we show that the invariants observed along the freezing and melting isomorphs are also observed on other isomorphs in the liquid and crystalline phase. Structure is probed by the radial distribution function and the structure factor and dynamics is probed by the mean-square displacement, the intermediate scattering function, and the shear viscosity. Studying these properties by reference to the isomorph theory explains why known single-phase melting criteria holds, e.g., the Hansen-Verlet and the Lindemann criterion, and why the Andrade equation for the viscosity at freezing applies, e.g., for most liquid metals. Our conclusion is that these empirical rules and invariants can all be understood from the isomorph theory and that the invariants are not peculiar to the freezing and melting lines, but hold along all isomorphs.

preprint2016arXiv

Studies of the Lennard-Jones fluid in 2, 3, and 4 dimensions highlight the need for a liquid-state 1/d expansion

The recent theoretical prediction by Maimbourg and Kurchan [arXiv:1603.05023] that for regular pair-potential systems the virial potential-energy correlation coefficient increases towards unity as the dimension $d$ goes to infinity is investigated for the standard 12-6 Lennard-Jones fluid. This is done by computer simulations for $d=2,3,4$ going from the critical point along the critical isotherm/isochore to higher density/temperature. In all cases the virial potential-energy correlation coefficient increases significantly. For a given density and temperature relative to the critical point, with increasing number of dimension the Lennard-Jones system conforms better to the hidden-scale-invariance property characterized by high virial potential-energy correlations (a property that leads to the existence of isomorphs in the thermodynamic phase diagram, implying that it becomes effectively one-dimensional in regard to structure and dynamics). The present paper also gives the first numerical demonstration of isomorph invariance of structure and dynamics in four dimensions. Our findings emphasize the need for a universally applicable $1/d$ expansion in liquid-state theory; we conjecture that the systems known to obey hidden scale invariance in three dimensions are those for which the yet-to-be-developed $1/d$ expansion converges rapidly.

preprint2016arXiv

Thermalization Calorimetry: A simple method for investigating glass transition and crystallization of supercooled liquids

We present a simple method for fast and cheap thermal analysis on supercooled glass-forming liquids. This "Thermalization Calorimetry" technique is based on monitoring the temperature and its rate of change during heating or cooling of a sample for which the thermal power input comes from heat conduction through an insulating material, i.e., is proportional to the temperature difference between sample and surroundings. The monitored signal reflects the sample's specific heat and is sensitive to exo- and endothermic processes. The technique is useful for studying supercooled liquids and their crystallization, e.g., for locating the glass transition and melting point(s), as well as for investigating the stability against crystallization and estimating the relative change in specific heat between the solid and liquid phases at the glass transition.

preprint2015arXiv

Continuum Nanofluidics

This paper introduces the fundamental continuum theory governing momentum transport in isotropic nanofluidic flows. The theory is an extension to the classical Navier-Stokes equation, which includes coupling between translational and rotational degrees of freedom, as well as non-local response functions that incorporates spatial correlations. The continuum theory is compared with molecular dynamics simulation data for both relaxation processes and fluid flows showing excellent agreement on the nanometer length scale. We also present practical tools to estimate when the extended theory should be used.

preprint2015arXiv

Direct tests of single-parameter aging

This paper presents accurate data for the physical aging of organic glasses just below the glass transition probed by monitoring the following quantities after temperature up and down jumps: the shear-mechanical resonance frequency (around 360 kHz), the dielectric loss at 1 Hz, the real part of the dielectric constant at 10 kHz, and the loss-peak frequency of the dielectric beta process (around 10 kHz). The setup used allows for keeping temperature constant within 100 micro Kelvin and for thermal equilibration within a few seconds after a temperature jump. The data conform to a new simplified version of the classical Tool-Narayanaswamy aging formalism, which makes it possible to calculate one relaxation curve directly from another without any fitting to analytical functions.

preprint2015arXiv

Explaining why simple liquids are quasi-universal

It has been known for a long time that many simple liquids have surprisingly similar structure as quantified, e.g., by the radial distribution function. A much more recent realization is that the dynamics are also very similar for a number of systems with quite different pair potentials. Systems with such non-trivial similarities are generally referred to as "quasi-universal". From the fact that the exponentially repulsive pair potential has strong virial potential-energy correlations in the low-temperature part of its thermodynamic phase diagram, we here show that a liquid is quasi-universal if its pair potential can be written approximately as a sum of exponential terms with numerically large prefactors. Based on evidence from the literature we moreover conjecture the converse, i.e., that quasi-universality only applies for systems with this property.

preprint2015arXiv

Hidden scale invariance of metals

Density functional theory (DFT) calculations of 58 liquid elements at their triple point show that most metals exhibit near proportionality between thermal fluctuations between virial and potential-energy in the isochoric ensemble. This demonstrates a general "hidden" scale invariance of metals making the dense part of the thermodynamic phase diagram effectively one dimensional with respect to structure and dynamics. DFT computed density scaling exponents, related to the Gr{ü}neisen parameter, are in good agreement with experimental values for 16 elements where reliable data were available. Hidden scale invariance is demonstrated in detail for magnesium by showing invariance of structure and dynamics. Computed melting curves of period three metals follow curves with invariance (isomorphs). The experimental structure factor of magnesium is predicted by assuming scale invariant inverse power-law (IPL) pair interactions. However, crystal packings of several transition metals (V, Cr, Mn, Fe, Nb, Mo, Ta, W and Hg), most post-transition metals (Ga, In, Sn, and Tl) and the metalloids Si and Ge cannot be explained by the IPL assumption. Thus, hidden scale invariance can be present even when the IPL-approximation is inadequate. The virial-energy correlation coefficient of iron and phosphorous is shown to increase at elevated pressures. Finally, we discuss how scale invariance explains the Gr{ü}neisen equation of state and a number of well-known empirical melting and freezing rules.

preprint2015arXiv

Invariants in the Yukawa system's thermodynamic phase diagram

This paper shows that several known properties of the Yukawa system can be derived from the isomorph theory, which applies to any system that has strong correlations between its virial and potential-energy equilibrium fluctuations. Such "Roskilde-simple" systems have a simplified thermodynamic phase diagram deriving from the fact that they have curves (isomorphs) along which structure and dynamics in reduced units are invariant to a good approximation. We show that the Yukawa system has strong virial potential-energy correlations and identify its isomorphs by two different methods. One method, the so-called direct isomorph check, identifies isomorphs numerically from jumps of relatively small density changes (here 10%). The second method identifies isomorphs analytically from the pair potential. The curves obtained by the two methods are close to each other; these curves are confirmed to be isomorphs by demonstrating the invariance of the radial distribution function, the static structure factor, the mean-square displacement as a function of time, and the incoherent intermediate scattering function. Since the melting line is predicted to be an isomorph, the theory provides a derivation of a known approximate analytical expression for this line in the temperature-density phase diagram. The paper's results give the first demonstration that the isomorph theory can be applied to systems like dense colloidal suspensions and strongly coupled dusty plasmas.

preprint2015arXiv

Narayanaswamy's 1971 aging theory and material time

The Bochkov-Kuzovlev nonlinear fluctuation-dissipation theorem is used to derive Narayanaswamy's phenomenological theory of physical aging, in which this highly nonlinear phenomenon is described by a linear material-time convolution integral. A characteristic property of the Narayanaswamy aging description is material-time translational invariance, which is here taken as the basic assumption of the derivation. It is shown that only one possible definition of the material time obeys this invariance, namely the square of the distance travelled from a configuration of the system far back in time. The paper concludes with suggestions for computer simulations that test for consequences of material-time translational invariance. One of these is the "unique-triangles property" according to which any three points on the system's path form a triangle such that two side lengths determine the third; this is equivalent to the well-known triangular relation for time-autocorrelation functions of aging spin glasses [L. F. Cugliandolo and J. Kurchan, J. Phys. A: Math. Gen. 27, 5749 (1994)]. The unique-triangles property implies a simple geometric interpretation of out-of-equilibrium time-autocorrelation functions, which extends to aging a previously proposed framework for such functions in equilibrium [J. C. Dyre, cond-mat/9712222].

preprint2015arXiv

Scaling of the dynamics of flexible Lennard-Jones chains. II. Effects of harmonic bonds

The previous paper [Veldhorst et al., J. Chem. Phys. 141, 054904 (2014)] demonstrated that the isomorph theory explains the scaling properties of a liquid of flexible chains consisting of ten Lennard-Jones particles connected by rigid bonds. We here investigate the same model with harmonic bonds. The introduction of harmonic bonds almost completely destroys the correlations in the equilibrium fluctuations of the potential energy and the virial. According to the isomorph theory, if these correlations are strong a system has isomorphs, curves in the phase diagram along which structure, dynamics and the excess entropy are invariant. The Lennard-Jones chain liquid with harmonic bonds does have curves in the phase diagram along which the structure and dynamics are invariant. The excess entropy is not invariant on these curves, which we refer to as "pseudoisomorphs". In particular this means that Rosenfeld's excess-entropy scaling (the dynamics being a function of excess entropy only) does not apply for the Lennard-Jones chain with harmonic bonds.

preprint2014arXiv

Cooee bitumen II: Stability of linear asphaltene nanoaggregates

Asphaltene and smaller aromatic molecules tend to form linear nanoaggregates in bitumen.Over the years bitumen undergoes chemical aging and during this process, the size of the nanoaggregate increases. This increase is associated with an increase in viscosity and brittleness of the bitumen, eventually leading to road deterioration. This paper focuses on understanding the mechanisms behind nanoaggregate size and stability. We used molecular dynamics simulations to quantify the probability of having a nanoaggregate of a given size in the stationary regime. To model this complicated behavior, we chose first to consider the simple case where only asphaltene molecules are counted in a nanoaggregate. We used a master equation approach and a related statistical mechanics model. The linear asphaltene nanoaggregates behave as a rigid linear chain. The most complicated case where all aromatic molecules are counted in a nanoaggregate is then discussed. The linear aggregates where all aromatic molecules are counted seem to behave as a flexible linear chain.

preprint2014arXiv

Density scaling and quasiuniversality of flow-event statistics for athermal plastic flows

Athermal plastic flows were simulated for the Kob-Andersen binary Lennard-Jones system and its repulsive version in which the sign of the attractive terms is changed to a plus. Properties evaluated from simulations at different densities include the distributions of energy drops, stress drops, and strain intervals between the flow events. By reference to hidden scale invariance we show that simulations at a single density in conjunction with an equilibrium-liquid simulation at the same density allows one to predict the plastic flow-event properties at other densities. We furthermore demonstrate quasiuniversality of the flow-event statistics.

preprint2014arXiv

Entropy and relaxation time

This paper discusses the possible relation between entropy and the relaxation time of liquids, in particular glass-forming systems, providing supplementing comments to the paper entitled "A brief critique of the Adam-Gibbs entropy model" by Hecksher et al. [J. Non-Cryst. Solids 325, 624-627 (2009)]. Besides the Adam-Gibbs model, we also briefly discuss Rosenfeld's excess entropy scaling and the 1964 experimental observation by Chang and Bestul of a universal excess entropy at the glass transition.

preprint2014arXiv

Estimating the density-scaling exponent of a monatomic liquid from its pair potential

This paper investigates two conjectures for calculating the density dependence of the density-scaling exponent of a single-component, pair-potential liquid with strong virial potential-energy correlations. The first conjecture gives an analytical expression for the density-scaling exponent directly in terms of the pair potential. The second conjecture is a refined version of this, which involves the most likely nearest-neighbor distance determined from the pair-correlation function. The two conjectures for the density-scaling exponent are tested by simulations of three systems, one of which is the standard Lennard-Jones liquid. While both expressions give qualitatively correct results, the second expression is more accurate.

preprint2014arXiv

Isomorph theory prediction for the dielectric loss variation along an isochrone

This paper derives a prediction for the variation of the amplitude of the dielectric loss from isomorph theory, and presents an experimental test of the prediction performed by measuring the dielectric-relaxation behavior of the van der Waals liquid 5-phenyl-4-ether (5PPE). The liquid is studied at isochronal states in the temperature range $266-333$ K and pressure range $0.1-300$ MPa, for relaxation times around $10^{-3}$ s and $10^{-4}$ s. From the isomorph statement that there is structural and dynamic invariance of isomorph states in reduced units for Roskilde simple liquids we derive four equivalent isomorph-invariant terms, one of which is used in analyzing our data. It is the frequency-dependent term $χ_{e}(f) ρ^{γ-1}$, with electric susceptibility $χ_{e}$, density $ρ$, and density-scaling factor $γ$. Due to the unique design of our experimental setup, we obtain dielectric loss data where the amplitude is reproducible $\pm 0.1 \%$. We moreover find that the empty capacitance of the capacitor cell is stable within $\pm 0.3 \%$ in our measuring range and can be assumed to be constant. Using this we predict for two isomorph states there is $-C_{2}"(f)= -C_{1}"(f) \left(ρ_1 / ρ_2 \right)^{γ-1}$ to scale the negative imaginary capacitance, where $C_{1}$ is the capacitance measurement at ambient pressure and $C_{2}$ is the predicted capacitance at elevated pressure. We visually compare the predicted and measured plots and there is good match between the two plots among the 42 pairs of isochronal states from the measurement.

preprint2014arXiv

Scaling of the dynamics of flexible Lennard-Jones chains

The isomorph theory provides an explanation for the so-called power law density scaling which has been observed in many molecular and polymeric glass formers, both experimentally and in simulations. Power law density scaling (relaxation times and transport coefficients being functions of $ρ^{γ_S}/T$, where $ρ$ is density, $T$ is temperature, and $γ_S$ is a material specific scaling exponent) is an approximation to a more general scaling predicted by the isomorph theory. Furthermore, the isomorph theory provides an explanation for Rosenfeld scaling (relaxation times and transport coefficients being functions of excess entropy) which has been observed in simulations of both molecular and polymeric systems. Doing molecular dynamics simulations of flexible Lennard-Jones chains (LJC) with rigid bonds, we here provide the first detailed test of the isomorph theory applied to flexible chain molecules. We confirm the existence of isomorphs, which are curves in the phase diagram along which the dynamics is invariant in the appropriate reduced units. This holds not only for the relaxation times but also for the full time dependence of the dynamics, including chain specific dynamics such as the end-to-end vector autocorrelation function and the relaxation of the Rouse modes. As predicted by the isomorph theory, jumps between different state points on the same isomorph happen instantaneously without any slow relaxation. Since the LJC is a simple coarse-grained model for alkanes and polymers, our results provide a possible explanation for why power-law density scaling is observed experimentally in alkanes and many polymeric systems. The theory provides an independent method of determining the scaling exponent, which is usually treated as a empirical scaling parameter.

preprint2014arXiv

Simplicity of condensed matter at its core: Generic definition of a Roskilde-simple system

The theory of isomorphs is reformulated by defining Roskilde-simple systems (those with isomorphs) by the property that the order of the potential energies of configurations at one density is maintained when these are scaled uniformly to a different density. Isomorphs remain curves in the thermodynamic phase diagram along which structure, dynamics, and excess entropy are invariant, implying that the phase diagram is effectively one-dimensional with respect to many reduced-unit properties. In contrast to the original formulation of the isomorph theory, however, the density-scaling exponent is not exclusively a function of density and the isochoric heat capacity is not an exact isomorph invariant. A prediction is given for the latter quantity's variation along the isomorphs. Molecular dynamics simulations of the Lennard-Jones and Lennard-Jones Gaussian systems validate the new approach.

preprint2014arXiv

The variation of the dynamic susceptibility along an isochrone

Koperwas {\it et al.} showed in a recent paper, Phys. Rev. Lett. {\bf 111}, 125701 (2013), that the dynamic susceptibility $χ_4$ as estimated by dielectric measurements for certain glass-forming liquids decreases substantially with increasing pressure along a curve of constant relaxation time. This observation is at odds with other measures of dynamics being invariant and seems to pose a problem for theories of glass formation. We show that this variation is in fact consistent with predictions for liquids with hidden scale invariance: measures of dynamics at constant volume are invariant along isochrones, called isomorphs in such liquids, but contributions to fluctuations from long-wavelength fluctuations can vary. This is related to the known non-invariance of the isothermal bulk modulus. Considering the version of $χ_4$ defined for the NVT ensemble, data from simulations of a binary Lennard-Jones liquid show in fact a slight increase with increasing density. This is a true departure from the formal invariance expected for this quantity.

preprint2013arXiv

Aging of CKN: Modulus versus conductivity analysis

It was recently reported that the electrical modulus peaks narrows upon annealing of the ionic system CKN [Paluch et al., Phys. Rev. Lett. 110, 015702 (2013)], which was interpreted as providing evidence of dynamic heterogeneity of this glass-forming liquid. An analysis of the same data in terms of the ac conductivity shows no shape changes, however. We discuss the relation between both findings and show further that the ac conductivity conforms to the prediction of the random barrier model (RBM) at all times during the annealing.

preprint2013arXiv

COOEE bitumen. Chemical aging

We study chemical aging in "COOEE bitumen" using molecular dynamic simulations. This model bitumen is composed of four realistic molecule types: saturated hydrocarbon, resinous oil, resin, and asphaltene. The aging reaction is modelled by the chemical reaction: "2 resins $\rightarrow$ 1 asphaltene". Molecular dynamic simulations of four bitumen compositions, obtained by a repeated application of the aging reaction, are performed. The stress autocorrelation function, the fluid structure, the rotational dynamics of the plane aromatic molecules, and the diffusivity of each molecule, are determined for the four different compositions. The aging reaction causes a significant dynamics slowdown, which is correlated to the aggregation of asphaltene molecules in larger and dynamically slower nanoaggregates. Finally, a detailed description of the role of each molecule types in the aggregation and aging processes is given.

preprint2013arXiv

Four-component united-atom model of bitumen

We propose a four-component molecular model of bitumen. The model includes realistic chemical constituents and introduces a coarse-graining level that suppresses the highest frequency modes. Molecular dynamics simulations of the model are being carried out using Graphic-Processor-Units based software in time spans in order of microseconds, and this enables the study of slow relaxation processes characterizing bitumen. This paper focuses on the high-temperature dynamics as expressed through the mean-square displacement, the stress autocorrelation function, and rotational relaxation. The diffusivity of the individual molecules changes little as a function of temperature and reveals distinct dynamical time scales as a result of the different constituents in the system. Different time scales are also observed for the rotational relaxation. The stress autocorrelation function features a slow non-exponential decay for all temperatures studied. From the stress autocorrelation function, the shear viscosity and shear modulus are evaluated at the highest temperature, showing a viscous response at frequencies below 100 MHz. The model predictions of viscosity and diffusivities are compared to experimental data, giving reasonable agreement. The model shows that the asphaltene, resin and resinous oil tend to form nano-aggregates. The characteristic dynamical relaxation time of these aggregates is different from the homogeneously distributed parts of the system, leading to strong dynamical heterogeneity.

preprint2013arXiv

Isomorph invariance of Couette shear flows simulated by the SLLOD equations of motion

Non-equilibrium molecular dynamics simulations were performed to study the thermodynamic, structural, and dynamical properties of the single-component Lennard-Jones and the Kob-Andersen binary Lennard-Jones liquids. Both systems are known to be strongly correlating, i.e., have strong correlations between equilibrium thermal fluctuations of virial and potential energy. Such systems have good isomorphs, i.e., curves in the thermodynamic phase diagram along which structural, dynamical, and some thermodynamic quantities are invariant when expressed in reduced units. The SLLOD equations of motion were used to simulate Couette shear flows of the two systems. We show analytically that these equations are isomorph invariant provided the reduced strain rate is fixed along the isomorph. Since isomorph invariance is generally only approximate, a range of shear rates were simulated to test for the predicted invariance, covering both the linear and non-linear regimes. For both systems, when represented in reduced units the radial distribution function and the intermediate scattering function collapse for state points that are isomorphic. The strain-rate dependence of the viscosity, which exhibits shear thinning, is also invariant along an isomorph. Our results extend the isomorph concept to the state-state non-equilibrium situation of a shear flow, in which the phase diagram is three dimensional because the shear rate defines the third dimension.

preprint2013arXiv

Isomorphs, hidden scale invariance, and quasiuniversality

This paper first establishes an approximate scaling property of the potential-energy function of a classical liquid with good isomorphs (a Roskilde-simple liquid). This "pseudohomogeneous" property makes explicit that - and in which sense - such a system has a hidden scale invariance. The second part gives a potential-energy formulation of the quasiuniversality of monatomic Roskilde-simple liquids, which was recently rationalized in terms of the existence of a quasiuniversal single-parameter family of reduced-coordinate constant-potential-energy hypersurfaces [J. C. Dyre, Phys. Rev. E 87, 022106 (2013)]. The new formulation involves a quasiuniversal reduced-coordinate potential-energy function. A few consequences of this are discussed.

preprint2013arXiv

NVU perspective on simple liquids' quasiuniversality

The last half century of research into the structure, dynamics, and thermodynamics of simple liquids has revealed a number of approximate universalities. This paper argues that simple liquids' reduced-coordinate constant-potential-energy hypersurfaces constitute a quasiuniversal family of compact Riemannian manifolds parameterized by a single number, from which follows these liquids' quasiuniversalities.

preprint2013arXiv

Predicting how nanoconfinement changes the relaxation time of a supercooled liquid

The properties of nanoconfined fluids can be strikingly different from those of bulk liquids. A basic unanswered question is whether the equilibrium and dynamic consequences of confinement are related to each other in a simple way. We study this question by simulation of a liquid comprising asymmetric dumbbell-shaped molecules, which can be deeply supercooled without crystallizing. We find that the dimensionless structural relaxation times - spanning six decades as a function of temperature, density, and degree of confinement - collapse when plotted versus excess entropy. The data also collapse when plotted versus excess isochoric heat capacity, a behaviour that follows from the existence of isomorphs in the bulk and confined states.

preprint2013arXiv

Statistical mechanics of Roskilde liquids: Configurational adiabats, specific heat contours and density dependence of the scaling exponent

We derive exact results for the rate of change of thermodynamic quantities, in particular the configurational specific heat at constant volume, $C_V$, along configurational adiabats (curves of constant excess entropy $S_{\textrm{ex}}$). Such curves are designated isomorphs for so-called Roskilde liquids, in view of the invariance of various structural and dynamical quantities along them. Their slope in a double logarithmic representation of the density-temperature phase diagram, $γ$ can be interpreted as one third of an effective inverse power-law potential exponent. We show that in liquids where $γ$ increases (decreases) with density, the contours of $C_V$ have smaller (larger) slope than configurational adiabats. We clarify also the connection between $γ$ and the pair potential. A fluctuation formula for the slope of the $C_V$-contours is derived. The theoretical results are supported with data from computer simulations of two systems, the Lennard-Jones fluid and the Girifalco fluid. The sign of $dγ/dρ$ is thus a third key parameter in characterizing Roskilde liquids, after $γ$ and the virial-potential energy correlation coefficient $R$. To go beyond isomorph theory we compare invariance of a dynamical quantity, the self-diffusion coefficient along adiabats and $C_V$-contours, finding it more invariant along adiabats.

preprint2013arXiv

The impact range for smooth wall-liquid interactions in nanoconfined liquids

Bulk and nanoconfined liquids have initially very different physics; for instance, nanoconfined liquids show stratification and position-dependent relaxation processes. A number of similarities between bulk and nanoconfined liquids have nevertheless been reported in computer simulations during the last decade. Inspired by these observations, we present results from molecular dynamics computer simulations of three nanoconfined liquids (i.e., single-component Lennard-Jones (LJ) liquid, Kob-Andersen binary LJ mixture, and an asymmetric dumbbell model) demonstrating also a microscopic similarity between bulk and nanoconfined liquids. The results show that the interaction range for the wall-liquid and liquid-liquid interactions of the nanoconfined liquid are identical to the bulk liquid as long as the liquid remains "Roskilde simple" in nanoconfinement, i.e., the liquid has strong correlations between virial and potential energy equilibrium fluctuations in the NVT ensemble.

preprint2013arXiv

The Rosenfeld-Tarazona expression for liquids' specific heat: A numerical investigation of eighteen systems

We investigate the accuracy of the expression of Rosenfeld and Tarazona (RT) for the excess isochoric heat capacity, C_V^{ex} \propto T^{-2/5}, for eighteen model liquids. Previous investigations have reported no unifying features of breakdown for the RT expression. Here liquids with different stoichiometric composition, molecular topology, chemical interactions, degree of undercooling, and environment are investigated. We find that the RT expression is a better approximation for liquids with strong correlations between equilibrium fluctuations of virial and potential energy, i.e., Roskilde simple liquids [Ingebrigtsen et al., Phys. Rev. X 2, 011011 (2012)]. This observation holds even for molecular liquids under severe nanoscale confinement, the physics of which is completely different from the original RT bulk hard-sphere fluid arguments. The density dependence of the specific heat is predicted from the isomorph theory for Roskilde simple liquids, which in combination with the RT expression provides a complete description of the specific heat's density and temperature dependence.

preprint2013arXiv

Two measures of isochronal superposition

A liquid obeys isochronal superposition if its dynamics is invariant along the isochrones in the thermodynamic phase diagram (the curves of constant relaxation time). This paper introduces two quantitative measures of isochronal superposition. The measures are used to test the following six liquids for isochronal superposition: 1,2,6 hexanetriol, glycerol, polyphenyl ether, diethyl phthalate, tetramethyl tetraphenyl trisiloxane, and dibutyl phthalate. The latter four van der Waals liquids obey isochronal superposition to a higher degree than the two hydrogen-bonded liquids. This is a predic- tion of the isomorph theory, and it confirms findings by other groups.

preprint2012arXiv

"Cooling by heating" - demonstrating the significance of the longitudinal specific heat

Heating a solid sphere at the surface induces mechanical stresses inside the sphere. If a finite amount of heat is supplied, the stresses gradually disappear as temperature becomes homogeneous throughout the sphere. We show that before this happens, there is a temporary lowering of pressure and density in the interior of the sphere, inducing a transient lowering of the temperature here. For ordinary solids this effect is small because c_p is almost equal to c_V. For fluent liquids the effect is negligible because their dynamic shear modulus vanishes. For a liquid at its glass transition, however, the effect is generally considerably larger than in solids. This paper presents analytical solutions of the relevant coupled thermoviscoelastic equations. In general, there is a difference between the isobaric specific heat, c_p, measured at constant isotropic pressure and the longitudinal specific heat, c_l, pertaining to mechanical boundary conditions that confine the associated expansion to be longitudinal. In the exact treatment of heat propagation the heat diffusion constant contains c_l rather than c_p. We show that the key parameter controlling the magnitude of the "cooling-by-heating" effect is the relative difference between these two specific heats. For a typical glass-forming liquid, when temperature at the surface is increased by 1 K, a lowering of the temperature in the sphere center of order 5 mK is expected if the experiment is performed at the glass transition. The cooling-by-heating effect is confirmed by measurements on a 19 mm diameter glucose sphere at the glass transition.

preprint2012arXiv

Do the repulsive and attractive pair forces play separate roles for the physics of liquids?

According to standard liquid-state theory repulsive and attractive pair forces play distinct roles for the physics of liquids. This paradigm is put into perspective here by demonstrating a continuous series of pair potentials that have virtually the same structure and dynamics, although only some of them have attractive forces of significance. Our findings reflect the fact that the motion of a given particle is determined by the total force on it, whereas the quantity usually discussed in liquid-state theory is the individual pair force.

preprint2012arXiv

Identical temperature dependence of the time scales of several linear-response functions of two glass-forming liquids

The frequency-dependent dielectric constant, shear and adiabatic bulk moduli, longitudinal thermal expansion coefficient, and longitudinal specific heat have been measured for two van der Waals glass-forming liquids, tetramethyl-tetraphenyl-trisiloxane (DC704) and 5-polyphenyl-4-ether. Within the experimental uncertainties the loss-peak frequencies of the measured response functions have identical temperature dependence over a range of temperatures, for which the Maxwell relaxation time varies more than nine orders of magnitude. The time scales are ordered from fastest to slowest as follows: Shear modulus, adiabatic bulk modulus, dielectric constant, longitudinal thermal expansion coefficient, longitudinal specific heat. The ordering is discussed in light of the recent conjecture that van der Waals liquids are strongly correlating, i.e., approximate single-parameter liquids.

preprint2012arXiv

Measuring the dynamic thermal expansivity of molecular liquids near the glass transition

Based on previous works on polymers by Bauer et al. [Phys, Rev. B (2000)], this paper describes a capacitative method for measuring the dynamical expansion coefficient of a viscous liquid. Data are presented for the glass-forming liquid tetramethyl tetraphenyl trisiloxane (DC704) in the ultraviscous regime. Compared to the method of Bauer et al. the dynamical range has been extended by making time-domain experiments and by making very small and fast temperature steps. The modelling of the experiment presented in this paper includes the situation where the capacitor is not full because the liquid contracts when cooling from room temperature down to around the glass-transition temperature, which is relevant when measuring on a molecular liquid rather than polymer.

preprint2012arXiv

NVU dynamics. III. Simulating molecules at constant potential energy

This is the final paper in a series that introduces geodesic molecular dynamics at constant potential energy. This dynamics is entitled NVU dynamics in analogy to standard energy-conserving Newtonian NVE dynamics. In the first two papers [Ingebrigtsen et al., J. Chem. Phys. 135, 104101 (2011); ibid, 104102 (2011)], a numerical algorithm for simulating geodesic motion of atomic systems was developed and tested against standard algorithms. The conclusion was that the NVU algorithm has the same desirable properties as the Verlet algorithm for Newtonian NVE dynamics, i.e., it is time-reversible and symplectic. Additionally, it was concluded that NVU dynamics becomes equivalent to NVE dynamics in the thermodynamic limit. In this paper, the NVU algorithm for atomic systems is extended to be able to simulate geodesic motion of molecules at constant potential energy. We derive an algorithm for simulating rigid bonds and test this algorithm on three different systems: an asymmetric dumbbell model, Lewis-Wahnstrom OTP, and rigid SPC/E water. The rigid bonds introduce additional constraints beyond that of constant potential energy for atomic systems. The rigid-bond NVU algorithm conserves potential energy, bond lengths, and step length for indefinitely long runs. The quantities probed in simulations give results identical to those of Nose-Hoover NVT dynamics. Since Nose-Hoover NVT dynamics is known to give results equivalent to those of NVE dynamics, the latter results show that NVU dynamics becomes equivalent to NVE dynamics in the thermodynamic limit also for molecular systems.

preprint2012arXiv

The instantaneous shear modulus in the shoving model

We point out that the instantaneous shear modulus of the shoving model for the non-Arrhenius temperature dependence of viscous liquids' relaxation time is the experimentally accessible high-frequency plateau modulus, not the idealized instantaneous affine shear modulus that cannot be measured. Data for a large selection of metallic glasses are compared to three different versions of the shoving model. The original shear-modulus based version shows a slight correlation to the Poisson ratio, which is eliminated by the energy-landscape formulation of the model in which the bulk modulus plays a minor role.

preprint2012arXiv

Thermodynamics of condensed matter with strong pressure-energy correlations

We show that for any liquid or solid with strong correlation between its $NVT$ virial and potential-energy equilibrium fluctuations, the temperature is a product of a function of excess entropy per particle and a function of density, $T=f(s)h(ρ)$. This implies that 1) the system's isomorphs (curves in the phase diagram of invariant structure and dynamics) are described by $h(ρ)/T={\rm Const.}$, 2) the density-scaling exponent is a function of density only, 3) a Gr{ü}neisen-type equation of state applies for the configurational degrees of freedom. For strongly correlating atomic systems one has $h(ρ)=\sum_nC_nρ^{n/3}$ in which the only non-zero terms are those appearing in the pair potential expanded as $v(r)=\sum_n v_n r^{-n}$. Molecular dynamics simulations of Lennard-Jones type systems confirm the theory.

preprint2012arXiv

What constitutes a simple liquid?

Simple liquids are traditionally defined as many-body systems of classical particles interacting via radially symmetric pair potentials. We suggest that a simple liquid should be defined instead by the property of having strong correlation between virial and potential energy equilibrium fluctuations in the NVT ensemble. There is considerable overlap between the two definitions, but also some notable differences. For instance, in the new definition simplicity is not a property of the intermolecular potential only because a liquid is usually only strongly correlating in part of its phase diagram. Moreover, according to the new definition not all simple liquids are atomic (i.e., with radially symmetric pair potentials) and not all atomic liquids are simple. The main part of the paper motivates the new definition of liquid simplicity by presenting evidence that a liquid is strongly correlating if and only if its intermolecular interactions may be ignored beyond the first coordination shell (FCS). This is demonstrated by NVT simulations of structure and dynamics of 15 atomic and molecular model liquids with a shifted-forces cutoff placed at the first minimum of the radial distribution function. No proof is given that the chemical characterization follows from the strong correlation property, but it is shown to be consistent with the existence of isomorphs in strongly correlating liquids' phase diagram. Finally, we note that the FCS characterization of simple liquids calls into question the basis for standard perturbation theory, according to which the repulsive and attractive forces play fundamentally different roles for the physics of liquids.

preprint2011arXiv

Beyond power-law density scaling: Theory, Simulation, and Experiment

Supercooled liquids are characterized by relaxation times that increase dramatically by cooling or compression. Many liquids have been shown to obey power-law density scaling, according to which the relaxation time is a function of density to some power over temperature. We show that power-law density scaling breaks down for larger density variations than usually studied. This is demonstrated by simulations of the Kob-Andersen binary Lennard-Jones mixture and two molecular models, as well as by experimental results for two van der Waals liquids. A more general form of density scaling is derived, which is consistent with results for all the systems studied. An analytical expression for the scaling function for liquids of particles interacting via generalized Lennard-Jones potentials is derived and shown to agree very well with simulations. This effectively reduces the problem of understanding the viscous slowing down from being a quest for a function of two variables to a search for a single-variable function.

preprint2011arXiv

Do all liquids become strongly correlating at high pressure?

We present molecular dynamics simulations studying the influence of pressure on the correlation between the constant-volume thermal equilibrium fluctuations of virial W and potential energy U, focusing on liquids that are not strongly correlating at low pressure, i.e., do not have a WU correlation coefficient above 0.9. The systems studied are the two hydrogen-bonded liquids GROMOS methanol and TIT5P water, the ionic liquid defined by a united-atom model of the 1-butyl-3-methyl-imidazolium nitrate and, for reference, the standard single-component Lennard-Jones liquid. The simulations were performed for pressures varying from 0 GPa to 10 GPa. For all systems studied we find that the virial / potential energy correlation increases with increasing pressure. This suggests that if crystallization is avoided, all liquids become strongly correlating at sufficiently high pressure.

preprint2011arXiv

Isomorphs in model molecular liquids

Isomorphs are curves in the phase diagram along which a number of static and dynamic quantities are invariant in reduced units. A liquid has good isomorphs if and only if it is strongly correlating, i.e., the equilibrium virial/potential energy fluctuations are more than 90% correlated in the NVT ensemble. This paper generalizes isomorphs to liquids composed of rigid molecules and study the isomorphs of two systems of small rigid molecules, the asymmetric dumbbell model and the Lewis-Wahnstrom OTP model. In particular, for both systems we find that the isochoric heat capacity, the excess entropy, the reduced molecular center-of-mass self part of the intermediate scattering function, the reduced molecular center-of-mass radial distribution function to a good approximation are invariant along an isomorph. In agreement with theory, we also find that an instantaneous change of temperature and density from an equilibrated state point to another isomorphic state point leads to no relaxation. The isomorphs of the Lewis-Wahnstrom OTP model were found to be more approximative than those of the asymmetric dumbbell model, which is consistent with the OTP model being less strongly correlating. For both models we find "master isomorphs", i.e., isomorphs have identical shape in the virial/potential energy phase diagram.

preprint2011arXiv

Isomorphs in the phase diagram of a model liquid without inverse power law repulsion

It is demonstrated by molecular dynamics simulations that liquids interacting via the Buckingham potential are strongly correlating, i.e., have regions of their phase diagram where constant-volume equilibrium fluctuations in the virial and potential energy are strongly correlated. A binary Buckingham liquid is cooled to a viscous phase and shown to have isomorphs, which are curves in the phase diagram along which structure and dynamics in appropriate units are invariant to a good approximation. To test this, the radial distribution function, and both the incoherent and coherent intermediate scattering function are calculated. The results are shown to reflect a hidden scale invariance; despite its exponential repulsion the Buckingham potential is well approximated by an inverse power-law plus a linear term in the region of the first peak of the radial distribution function. As a consequence the dynamics of the viscous Buckingham liquid is mimicked by a corresponding model with purely repulsive inverse-power-law interactions. The results presented here closely resemble earlier results for Lennard-Jones type liquids, demonstrating that the existence of strong correlations and isomorphs does not depend critically on the mathematical form of the repulsion being an inverse power law.

preprint2011arXiv

NVU dynamics. I. Geodesic motion on the constant-potential-energy hypersurface

An algorithm is derived for computer simulation of geodesics on the constant potential-energy hypersurface of a system of N classical particles. First, a basic time-reversible geodesic algorithm is derived by discretizing the geodesic stationarity condition and implementing the constant potential energy constraint via standard Lagrangian multipliers. The basic NVU algorithm is tested by single-precision computer simulations of the Lennard-Jones liquid. Excellent numerical stability is obtained if the force cutoff is smoothed and the two initial configurations have identical potential energy within machine precision. Nevertheless, just as for NVE algorithms, stabilizers are needed for very long runs in order to compensate for the accumulation of numerical errors that eventually lead to "entropic drift" of the potential energy towards higher values. A modification of the basic NVU algorithm is introduced that ensures potential-energy and step-length conservation; center-of-mass drift is also eliminated. Analytical arguments confirmed by simulations demonstrate that the modified NVU algorithm is absolutely stable. Finally, simulations show that the NVU algorithm and the standard leap-frog NVE algorithm have identical radial distribution functions for the Lennard-Jones liquid.

preprint2011arXiv

NVU dynamics. II. Comparing to four other dynamics

In the companion paper [Ingebrigtsen et al., arXiv:1012.3447] an algorithm was developed for tracing out a geodesic curve on the constant-potential-energy hypersurface. Here simulations of this NVU dynamics are compared to results for four other dynamics, both deterministic and stochastic. First, NVU dynamics is compared to the standard energy-conserving Newtonian NVE dynamics by simulations of the Kob-Andersen binary Lennard-Jones liquid, its WCA version (i.e., with cut-off's at the pair potential minima), and the Gaussian Lennard-Jones liquid. We find identical results for all quantities probed: radial distribution functions, incoherent intermediate scattering functions, and mean-square displacement as function of time. Arguments are then presented for the equivalence of NVU and NVE dynamics in the thermodynamic limit; in particular to leading order in 1/N these two dynamics give identical time-autocorrelation functions. In the final section NVU dynamics is compared to Monte Carlo dynamics, to a diffusive dynamics of small-step random walks on the constant-potential-energy hypersurface, and to Nose-Hoover NVT dynamics. If time is scaled for the two stochastic dynamics to make their single-particle diffusion constants identical to those of NVE dynamics, the simulations show that all five dynamics are equivalent at low temperatures except at short times.

preprint2011arXiv

Pressure-energy correlations in liquids. V. Isomorphs in generalized Lennard-Jones systems

This series of papers is devoted to identifying and explaining the properties of strongly correlating liquids, i.e., liquids with more than 90% correlation between their virial W and potential energy U fluctuations in the NVT ensemble. Paper IV [N. Gnan et al., J. Chem. Phys. v131, 234504 (2009)] showed that strongly correlating liquids have "isomorphs", which are curves in the phase diagram along which structure, dynamics, and some thermodynamic properties are invariant in reduced units. In the present paper, using the fact that reduced-unit radial distribution functions are isomorph invariant, we derive an expression for the shapes of isomorphs in the WU phase diagram of generalized Lennard-Jones systems of one or more types of particles. The isomorph shape depends only on the Lennard-Jones exponents; thus all isomorphs of standard Lennard-Jones systems (with exponents 12 and 6) can be scaled onto to a single curve. Two applications are given. One is testing the prediction that the solid-liquid coexistence curve follows an isomorph by comparing to recent simulations by Ahmed and Sadus [J. Chem. Phys. v131, 174504 (2009)]. Excellent agreement is found on the liquid side of the coexistence, whereas the agreement is worse on the solid side. A second application is the derivation of an approximate equation of state for generalized Lennard-Jones systems by combining the isomorph theory with the Rosenfeld-Tarazona expression for the temperature dependence of potential energy on isochores. It is shown that the new equation of state agrees well with simulations.

preprint2011arXiv

Role of the first coordination shell in determining the equilibrium structure and dynamics of simple liquids

The traditional view that the physical properties of a simple liquid are determined primarily by its repulsive forces was recently challenged by Berthier and Tarjus, who showed that in some cases ignoring the attractions leads to large errors in the dynamics [L. Berthier and G. Tarjus, Phys. Rev. Lett. 103, 170601 (2009); J. Chem. Phys. 134, 214503 (2011)]. We present simulations of the standard Lennard-Jones liquid at several condensed-fluid state points, including a fairly low density state and a very high density state, as well as simulations of the Kob-Andersen binary Lennard-Jones mixture at several temperatures. By varying the range of the forces, results for the thermodynamics, dynamics, and structure show that the determining factor for getting the correct statics and dynamics is not whether or not the attractive forces {\it per se} are included in the simulations. What matters is whether or not interactions are included from all particles within the first coordination shell (FCS) - the attractive forces can thus be ignored, but only at extremely high densities. The recognition of the importance of a local shell in condensed fluids goes back to van der Waals; our results confirm this idea and thereby the basic picture of the old hole- and cell theories for simple condensed fluids.

preprint2011arXiv

Shifted forces in molecular dynamics

Simulations involving the Lennard-Jones potential usually employ a cut-off at $r=2.5σ$. This paper investigates the possibility of reducing the cut-off. Two different cut-off implementations are compared, the standard shifted potential cut-off and the less commonly used shifted forces cut-off. The first has correct forces below the cut-off, whereas the shifted forces cut-off modifies Newton's equations at all distances. The latter is nevertheless superior; we find that for most purposes realistic simulations may be obtained using a shifted-forces cut-off at $r=1.5σ$, even though the pair force is here 30 times larger than at $r=2.5σ$.

preprint2011arXiv

Simplistic Coulomb forces in molecular dynamics: Comparing the Wolf and shifted-force approximations

This paper compares the Wolf method to the shifted forces (SF) method for efficient computer simulation of isotropic systems interacting via Coulomb forces, taking results from the Ewald summation method as representing the true behavior. We find that for the Hansen-McDonald molten salt model the SF approximation overall reproduces the structural and dynamical properties as accurately as does the Wolf method. It is shown that the optimal Wolf damping parameter depends on the property in focus, and that neither the potential energy nor the radial distribution function are useful measures for the convergence of theWolf method to the Ewald summation method. The SF approximation is also tested for the SPC/Fw model of liquid water at room temperature, showing good agreement with both the Wolf and the particle mesh Ewald methods; this confirms previous findings [Fennell & Gezelter, J. Chem. Phys. {\bf 124}, 234104 (2006)]. Beside its conceptualsimplicity the SF approximation implies a speed-up of a factor 2 to 3 compared to the Wolf method (which is in turn much faster than the Ewald method).

preprint2010arXiv

A repulsive reference potential reproducing the dynamics of a liquid with attractions

A well-known result of liquid state theory is that the structure of dense fluids is mainly determined by repulsive forces. The WCA potential, which cuts intermolecular potentials at their minima, is therefore often used as a reference. However, this reference gives quite wrong results for the viscous dynamics of the Kob-Andersen binary Lennard-Jones liquid [Berthier and Tarjus, Phys. Rev. Lett. 103, 170601 (2009)]. We show that repulsive inverse-power law potentials provide a useful reference for this liquid by reproducing its structure, dynamics, and isochoric heat capacity.

preprint2010arXiv

Aging effects manifested in the potential energy landscape of a model glass former

We present molecular dynamics simulations of a model glass-forming liquid (the binary Kob-Anderson Lennard-Jones model) and consider the distributions of inherent energies and metabasins during aging. In addition to the typical protocol of performing a temperature jump from a high temperature to a low destination temperature, we consider the temporal evolution of the distributions after an 'up-jump', i.e. from a low to a high temperature. In this case the distribution of megabasin energies exhibits a transient two-peak structure. Our results can qualitatively be rationalized in terms of a trap model with a Gaussian distribution of trap energies. The analysis is performed for different system sizes. A detailed comparison with the trap model is possible only for a small system because of major averging effects for larger systems.

preprint2010arXiv

Beta relaxation in the shear mechanics of equilibrium viscous liquids: Phenomenology and network modeling of the alpha-beta merging region

The phenomenology of the beta relaxation process in the shear-mechanical response of glass-forming liquids is summarized and compared to that of the dielectric beta process. Furthermore, we discuss how to model the observations by means of standard viscoelastic modeling elements. Necessary physical requirements to such a model are outlined, and it is argued that physically relevant models must be additive in the shear compliance of the alpha and beta parts. A model based on these considerations is proposed and fitted to data for Polyisobutylene 680.

preprint2010arXiv

Physical aging studied by a device allowing for rapid thermal equilibration

Aging to the equilibrium liquid state of organic glasses is studied. The glasses were prepared by cooling the liquid to temperatures just below the glass transition. Aging following a temperature jump was studied by measuring the dielectric loss at a fixed frequency using a microregulator in which temperature is controlled by means of a Peltier element. Compared to conventional equipment the new device adds almost two orders of magnitude to the span of observable aging times. Data for the following five glass-forming liquids are presented: Dibutyl phthalate, diethyl phthalate, 2,3-epoxy propyl-phenyl-ether, 5-polyphenyl-ether, and triphenyl phosphite. The aging data were analyzed using the Tool-Narayanaswamy formalism. The following features are found for all five liquids: 1) Each liquid has an "internal clock", a fact that is established by showing that the aging of the structure is controlled by the same material time that controls the dielectric properties. 2) There are no so-called expansion gaps between the long-time limits of the relaxation rates following up and down jumps to the same temperature. 3) At long times the structural relaxation is not stretched, but a simple exponential decay. 4) For small temperature steps the rate of the long-time exponential structural relaxation is identical to that of the long-time decay of the dipole autocorrelation function.

preprint2010arXiv

Predicting The Effective Temperaure of a Glass

We explain the findings by Di Leonardo et al. [Phys. Rev. Lett. 84, 6054 (2000)] that the effective temperature of a Lennard-Jones glass depends only on the final value of the density in the volume and/or temperature jump that produces the glass phase. This is not only a property of the Lennard-Jones liquid, but a feature of all strongly correlating liquids. For such liquids data from a single quench simulation provides enough information to predict the effective temperature of any glass produced by jumping from an equilibrium state. This prediction is validated by simulations of the Kob-Andersen binary Lennard-Jones liquid and shown not to apply for the non-strongly correlating monatomic Lennard-Jones Gaussian liquid.

preprint2010arXiv

Strongly correlating liquids and their isomorphs

This paper summarizes the properties of strongly correlating liquids, i.e., liquids with strong correlations between virial and potential energy equilibrium fluctuations at constant volume. We proceed to focus on the experimental predictions for strongly correlating glass-forming liquids. These predictions include i) density scaling, ii) isochronal superposition, iii) that there is a single function from which all frequency-dependent viscoelastic response functions may be calculated, iv) that strongly correlating liquids are approximately single-parameter liquids with close to unity Prigogine-Defay ratio, and v) that the fictive temperature initially decreases for an isobaric temperature up jump. The "isomorph filter", which allows one to test for universality of theories for the non-Arrhenius temperature dependence of the relaxation time, is also briefly discussed.

preprint2010arXiv

Unique Dynamic Correlation Length in Supercooled Liquids

We present a direct comparison of the number of dynamically correlated molecules in the shear-mechanical and dielectric relaxations of the following seven supercooled organic liquids: triphenylethylene, tetramethyl-tetraphenyl-trisiloxane, polyphenyl ether, perhydrosqualene, polybutadiene, decahydroisoquinoline, and tripropylene glycol. For each liquid we observe that the numbers of dynamically correlated molecules in the shear and in the dielectric relaxation are proportional. We show that this proportionality can be explained by the constancy of the decoupling index of the shear and dielectric relaxation times in conjunction with time-temperature superposition. Moreover the value of this proportionality constant is related to the difference in stretching of the shear and dielectric response functions. The most significant deviations from unity of this constant are found in a liquid with strong hydrogen bonds and in a polymer.

preprint2009arXiv

"Isomorphs" in liquid state diagrams

A liquid is termed strongly correlating if its virial and potential energy thermal equilibrium fluctuations in the NVT ensemble are more than 90% correlated [Phys. Rev. Lett. 100, 015701 (2008)]. The fluctuations of a strongly correlating liquid are well approximated by those of an inverse power-law intermolecular potential. Building on this fact we here define "isomorphic lines" in the state diagram of a strongly correlating liquid. It is shown from computer simulations of the Kob-Andersen binary Lennard-Jones liquid that no aging is associated with jumps between two isomorphic points. Isomorphic state points have the same excess entropy, the same reduced average relaxation time, the same (reduced) dynamics, and the same scaled radial distribution functions. Finally we calculate the equation for isomorphs in the virial / potential energy diagram for Lennard-Jones type liquids and show that all such 12-6 liquids have the same isomorphs; these may be scaled into a master isomorph.

preprint2009arXiv

A brief critique of the Adam-Gibbs entropy model

This paper critically reviews the entropy model proposed by Adam and Gibbs in 1965 for explaining the dramatic temperature dependence of glass-forming liquids' average relaxation time, one of the most influential models during the last three decades. We discuss the Adam-Gibbs model's theoretical bases as well as the reported experimental model confirmations; in the process of doing this a number of problems with the model are identified.

preprint2009arXiv

Generalized fluctuation-dissipation relation and effective temperature in off-equilibrium colloids

The fluctuation-dissipation relation (FDR), a fundamental result of equilibrium statistical physics, ceases to be valid when a system is taken out of the equilibrium. A generalization of FDR has been theoretically proposed for out-of-equilibrium systems: the kinetic temperature entering FDR is substituted by a time-scale dependent effective temperature. We combine the measurements of the correlation function of the rotational dynamics of colloidal particles obtained via dynamic light scattering with those of the birefringence response to study the generalized FDR in an off-equilibrium clay suspension undergoing aging. We i) find that FDR is strongly violated in the early stage of the aging process and is gradually recovered as the aging time increases and, ii), we determine the aging time evolution of the effective temperature, giving support to the proposed generalization scheme.

preprint2009arXiv

Pressure-energy correlations in liquids. III. Statistical mechanics and thermodynamics of liquids with hidden scale invariance

In this third paper of the series, which started with [N. P. Bailey et al., J. Chem. Phys. 129, 184507 and 184508 (2008)], we continue the development of the theoretical understanding of strongly correlating liquids - those whose instantaneous potential energy and virial are strongly correlated in their thermal equilibrium fluctuations at constant volume. The existence of such liquids was detailed in previous work which identified them, based on computer simulations, as a large class of liquids, including van der Waals liquids but not, e.g., hydrogen-bonded liquids. We here discuss the following: (1) The scaling properties of inverse power-law and extended inverse power-law potentials (the latter include a linear term which "hides" the approximate scale invariance); (2) results from computer simulations of molecular models concerning out-of-equilibrium conditions; (3) ensemble dependence of the virial / potential energy correlation coefficient; (4) connection to the Gruneisen parameter; (5) interpretation of strong correlations in terms of the energy-bond formalism.

preprint2009arXiv

Pressure-energy correlations in liquids. IV. "Isomorphs" in liquid state diagrams

This paper is the fourth in a series devoted to identifying and explaining the properties of strongly correlating liquids, i.e., liquids where virial and potential energy correlate better than 90% in their thermal equilibrium fluctuations in the NVT ensemble. For such liquids we here introduce the concept of "isomorphic" curves in the state diagram. A number of thermodynamic, static, and dynamic isomorph invariants are identified. These include the excess entropy, the isochoric specific heat, reduced-unit static and dynamic correlation functions, as well as reduced-unit transport coefficients. The dynamic invariants apply for both Newtonian and Brownian dynamics. It is shown that after a jump between isomorphic state points the system is instantaneously in thermal equilibrium; consequences of this for generic aging experiments are discussed. Selected isomorph predictions are validated by computer simulations of the Kob-Andersen binary Lennard-Jones mixture, which is a strongly correlating liquid. The final section of the paper relates the isomorph concept to phenomenological melting rules, Rosenfeld's excess entropy scaling, Young and Andersen's approximate scaling principle, and the two-order parameter maps of Debenedetti and coworkers. This section also shows how the existence of isomorphs implies an "isomorph filter" for theories for the non-Arrhenius temperature dependence of viscous liquids' relaxation time, as well as explains isochronal superposition for strongly correlating viscous liquids.

preprint2009arXiv

The Geometry of Slow Structural Fluctuations in a Supercooled Binary Alloy

The liquid structure of a glass-forming binary alloy is studied using molecular dynamics simulations. The analysis combines common neighbour analysis with the geometrical approach of Frank and Kasper to establish that the supercooled liquid contains extended clusters characterised by the same short range order as the crystal. Fluctuations in these clusters exhibit strong correlations with fluctuations in the inherent structure energy. The steep increase in the heat capacity on cooling is, thus, directly coupled to the growing fluctuations of the Frank-Kasper clusters. The relaxation of particles in the clusters dominates the slow tail of the self-intermediate scattering function.

preprint2008arXiv

Ac hopping conduction at extreme disorder takes place on the percolating cluster

Simulations of the random barrier model show that ac currents at extreme disorder are carried almost entirely by the percolating cluster slightly above threshold; thus contradicting traditional theories contributions from isolated low-activation-energy clusters are negligible. The effective medium approximation in conjunction with the Alexander-Orbach conjecture leads to an excellent analytical fit to the universal ac conductivity with no nontrivial fitting parameters.

preprint2008arXiv

Crystallization of the Kob-Andersen binary Lennard-Jones liquid

At density 1.2 the Kob-Andersen binary Lennard-Jones liquid partly crystallizes in the temperature interval [0.39, 0.44] after typically 10-100 microseconds (Argon units). The crystallization is initiated by a phase separation where the large (A) particles cluster in a volume void of B particles. We investigate a modification of the Kob-Andersen system where the attraction between particles of the same type is removed, thus disfavoring phase separation. We have not been able to crystallize this new system.

preprint2008arXiv

Density scaling as a property of strongly correlating viscous liquids

We address a recent conjecture according to which the relaxation time $τ$ of a viscous liquid obeys density scaling ($τ=F(ρ^γ/T)$ where $ρ$ is density) if the liquid is ``strongly correlating,'' i.e., has almost 100% correlation between equilibrium virial and potential-energy fluctuations [Pedersen {\it et al.}, PRL {\bf 100}, 011201 (2008)]. Computer simulations of two model liquids - an asymmetric dumbbell model and the Lewis-Wahnström OTP model - confirm the conjecture and demonstrate that the scaling exponent $γ$ can be accurately predicted from equilibrium fluctuations.

preprint2008arXiv

Supercooled Liquid Dynamics Studied via Shear-Mechanical Spectroscopy

We report dynamical shear-modulus measurements for five glass-forming liquids (pentaphenyl trimethyl trisiloxane, diethyl phthalate, dibutyl phthalate, 1,2-propanediol, and m-touluidine). The shear-mechanical spectra are obtained by the piezoelectric shear-modulus gauge (PSG) method. This technique allows one to measure the shear modulus ($10^{5} -10^{10}$ Pa) of the liquid within a frequency range from 1 mHz to 10 kHz. We analyze the frequency-dependent response functions to investigate whether time-temperature superposition (TTS) is obeyed. We also study the shear-modulus loss-peak position and its high-frequency part. It has been suggested that when TTS applies, the high-frequency side of the imaginary part of the dielectric response decreases like a power law of the frequency with an exponent -1/2. This conjecture is analyzed on the basis of the shear mechanical data. We find that TTS is obeyed for pentaphenyl trimethyl trisiloxane and in 1,2-propanediol while in the remaining liquids evidence of a mechanical $β$ process is found. Although the the high-frequency power law behavior $ω^{-α}$ of the shear-loss may approach a limiting value of $α=0.5$ when lowering the temperature, we find that the exponent lies systematically above this value (around 0.4). For the two liquids without beta relaxation (pentaphenyl trimethyl trisiloxane and 1,2-propanediol) we also test the shoving model prediction, according to which the the relaxation-time activation energy is proportional to the instantaneous shear modulus. We find that the data are well described by this model.

preprint2007arXiv

Approximate square-root-time relaxation in glass-forming liquids

We present data for the dielectric relaxation of 43 glass-forming organic liquids, showing that the primary (alpha) relaxation is often close to square-root-time relaxation. The better an inverse power-law description of the high-frequency loss applies, the more accurately is square-root-time relaxation obeyed. These findings suggest that square-root-time relaxation is generic to the alpha process, once a common view, but since long believed to be incorrect. Only liquids with very large dielectric losses deviate from this picture by having consistently narrower loss peaks. As a further challenge to the prevailing opinion, we find that liquids with accurate square-root-time relaxation cover a wide range of fragilities.

preprint2007arXiv

Conventional methods fail to measure c_p(omega) of glass-forming liquids

The specific heat is frequency dependent in highly viscous liquids. By solving the full one-dimensional thermo-viscoelastic problem analytically it is shown that, because of thermal expansion and the fact that mechanical stresses relax on the same time scale as the enthalpy relaxes, the plane thermal-wave method does not measure the isobaric frequency-dependent specific heat c_p(omega). This method rather measures a "longitudinal" frequency-dependent specific heat, a quantity defined and detailed here that is in-between c_p(omega) and c_v(omega). This result means that no wide-frequency measurements of c_p(omega) on liquids approaching the calorimetric glass transition exist. We briefly discuss consequences for experiment.

preprint2007arXiv

Volume-energy correlations in the slow degrees of freedom of computer-simulated phospholipid membranes

Constant-pressure molecular-dynamics simulations of phospholipid membranes in the fluid phase reveal strong correlations between equilibrium fluctuations of volume and energy on the nanosecond time-scale. The existence of strong volume-energy correlations was previously deduced indirectly by Heimburg from experiments focusing on the phase transition between the fluid and the ordered gel phases. The correlations, which are reported here for three different membranes (DMPC, DMPS-Na, and DMPSH), have volume-energy correlation coefficients ranging from 0.81 to 0.89. The DMPC membrane was studied at two temperatures showing that the correlation coefficient increases as the phase transition is approached.

preprint2006arXiv

An energy landscape model for glass forming liquids in three dimensions

We present a three-dimensional lattice-gas model with trivial thermodynamics, but nontrivial dynamics. The model is characterized by each particle having its own random energy landscape. The equilibrium dynamics of the model were investigated by continuous time Monte Carlo simulations at two different densities at several temperatures. At high densities and low temperatures the model captures the important characteristics of viscous liquid dynamics. We thus observe non-exponential relaxation in the self part of the density auto-correlation function, and fragility plots of the self-diffusion constant and relaxation times show non-Arrhenius behavior.

preprint2006arXiv

Elastic models for the non-Arrhenius viscosity of glass-forming liquids

This paper first reviews the shoving model for the non-Arrhenius viscosity of viscous liquids. According to this model the main contribution to the activation energy of a flow event is the energy needed for molecules to shove aside the surrounding, an energy which is proportional to the instantaneous shear modulus of the liquid. Data are presented supporting the model. It is shown that the fractional Debye-Stokes-Einstein relation, that quantitatively expresses the frequently observed decoupling of, e.g., conductivity from viscous flow, may be understood within the model. The paper goes on to review several related explanations for the non-Arrhenius viscosity. Most of these are also "elastic models," i.e., they express the viscosity activation energy in terms of short-time elastic properties of the liquid. Finally, two new arguments for elastic models are given, a general solid-state defect argument and an Occam's razor type argument.

preprint2001arXiv

Hopping models and ac universality

Some general relations for hopping models are established. We proceed to discuss the universality of the ac conductivity which arises in the extreme disorder limit of the random barrier model. It is shown that the relevant dimension entering into the diffusion cluster approximation (DCA) is the harmonic (fracton) dimension of the diffusion cluster. The temperature scaling of the dimensionless frequency entering into the DCA is discussed. Finally, some open questions about ac universality are mentioned.

preprint1999arXiv

Crossover to Potential Energy Landscape Dominated Dynamics in a Model Glass-forming Liquid

An equilibrated model glass-forming liquid is studied by mapping successive configurations produced by molecular dynamics simulation onto a time series of inherent structures (local minima in the potential energy). Using this ``inherent dynamics'' approach we find direct numerical evidence for the long held view that below a crossover temperature, $T_x$, the liquid's dynamics can be separated into (i) vibrations around inherent structures and (ii) transitions between inherent structures (M. Goldstein, J. Chem. Phys. {\bf 51}, 3728 (1969)), i.e., the dynamics become ``dominated'' by the potential energy landscape. In agreement with previous proposals, we find that $T_x$ is within the vicinity of the mode-coupling critical temperature $T_c$. We further find that at the lowest temperature simulated (close to $T_x$), transitions between inherent structures involve cooperative, string like rearrangements of groups of particles moving distances substantially smaller than the average interparticle distance.

preprint1997arXiv

Why are glass-forming liquids non-Arrhenius?

A major mystery of glass-forming liquids is the non-Arrhenius temperature-dependence of the average relaxation time. This paper briefly reviews the classical phenomenological models for this phenomenon - the free-volume model and the entropy model - and critiques against these models. We then discuss a recent model [Dyre, Olsen, and Christensen, Phys. Rev. B 53, 2171 (1996)] according to which the activation energy for the average relaxation time is determined by the work done in shoving aside the surrounding liquid to create space needed for a flow event. In this model the non-Arrhenius temperature-dependence is a consequence of the fact that the instantaneous (infinite-frequency) shear modulus increases upon cooling.

Jeppe C. Dyre

What is connected

Connect this record

See the researcher in context

Building this map preview

89 published item(s)

Active-matter isomorphs in the size-polydisperse Ornstein-Uhlenbeck Lennard-Jones model

Active-parameter polydispersity in the 2d ABP Yukawa model

Hidden scale invariance in the Gay-Berne model

Isomorph invariance in the liquid and plastic-crystal phases of asymmetric-dumbbell models

Predicting nonlinear physical aging of glasses from equilibrium relaxation via the material time

Rheological model for the alpha relaxation of glass-forming liquids and its comparison to data for DC704 and DC705

Single-parameter aging in the weakly nonlinear limit

Generalized hydrodynamics of the Lennard-Jones liquid in view of hidden scale invariance

Isomorphs in nanoconfined liquids

Single-parameter aging in a binary Lennard-Jones system

Time-scale coupling in hydrogen- and van der Waals-bonded liquids

Solid-like mean-square displacement in glass-forming liquids

Solid-liquid coexistence of the noble elements. I. Theory illustrated by the case of argon

Solid-liquid coexistence of the noble elements. II. Neon, krypton and xenon

Structure of the Lennard-Jones liquid estimated from a single simulation

"Pseudoisomorphs" in liquids with intramolecular degrees of freedom

A pair potential that reproduces the shape of isochrones in molecular liquids

Freezing and melting line invariants of the Lennard-Jones system

Studies of the Lennard-Jones fluid in 2, 3, and 4 dimensions highlight the need for a liquid-state 1/d expansion

Thermalization Calorimetry: A simple method for investigating glass transition and crystallization of supercooled liquids

Continuum Nanofluidics

Direct tests of single-parameter aging

Explaining why simple liquids are quasi-universal

Hidden scale invariance of metals

Invariants in the Yukawa system's thermodynamic phase diagram

Narayanaswamy's 1971 aging theory and material time

Scaling of the dynamics of flexible Lennard-Jones chains. II. Effects of harmonic bonds

Cooee bitumen II: Stability of linear asphaltene nanoaggregates

Density scaling and quasiuniversality of flow-event statistics for athermal plastic flows

Entropy and relaxation time

Estimating the density-scaling exponent of a monatomic liquid from its pair potential

Isomorph theory prediction for the dielectric loss variation along an isochrone

Scaling of the dynamics of flexible Lennard-Jones chains

Simplicity of condensed matter at its core: Generic definition of a Roskilde-simple system

The variation of the dynamic susceptibility along an isochrone

Aging of CKN: Modulus versus conductivity analysis

COOEE bitumen. Chemical aging

Four-component united-atom model of bitumen

Isomorph invariance of Couette shear flows simulated by the SLLOD equations of motion

Isomorphs, hidden scale invariance, and quasiuniversality

NVU perspective on simple liquids' quasiuniversality

Predicting how nanoconfinement changes the relaxation time of a supercooled liquid

Statistical mechanics of Roskilde liquids: Configurational adiabats, specific heat contours and density dependence of the scaling exponent

The impact range for smooth wall-liquid interactions in nanoconfined liquids

The Rosenfeld-Tarazona expression for liquids' specific heat: A numerical investigation of eighteen systems

Two measures of isochronal superposition

"Cooling by heating" - demonstrating the significance of the longitudinal specific heat

Do the repulsive and attractive pair forces play separate roles for the physics of liquids?

Identical temperature dependence of the time scales of several linear-response functions of two glass-forming liquids

Measuring the dynamic thermal expansivity of molecular liquids near the glass transition

NVU dynamics. III. Simulating molecules at constant potential energy

The instantaneous shear modulus in the shoving model

Thermodynamics of condensed matter with strong pressure-energy correlations

What constitutes a simple liquid?

Beyond power-law density scaling: Theory, Simulation, and Experiment

Do all liquids become strongly correlating at high pressure?

Isomorphs in model molecular liquids

Isomorphs in the phase diagram of a model liquid without inverse power law repulsion

NVU dynamics. I. Geodesic motion on the constant-potential-energy hypersurface

NVU dynamics. II. Comparing to four other dynamics

Pressure-energy correlations in liquids. V. Isomorphs in generalized Lennard-Jones systems

Role of the first coordination shell in determining the equilibrium structure and dynamics of simple liquids

Shifted forces in molecular dynamics

Simplistic Coulomb forces in molecular dynamics: Comparing the Wolf and shifted-force approximations

A repulsive reference potential reproducing the dynamics of a liquid with attractions

Aging effects manifested in the potential energy landscape of a model glass former

Beta relaxation in the shear mechanics of equilibrium viscous liquids: Phenomenology and network modeling of the alpha-beta merging region

Physical aging studied by a device allowing for rapid thermal equilibration

Predicting The Effective Temperaure of a Glass

Strongly correlating liquids and their isomorphs

Unique Dynamic Correlation Length in Supercooled Liquids

"Isomorphs" in liquid state diagrams

A brief critique of the Adam-Gibbs entropy model

Generalized fluctuation-dissipation relation and effective temperature in off-equilibrium colloids