Researcher profile

Jeppe C. Dyre

Jeppe C. Dyre contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
5topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

15 published item(s)

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.