BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Thomas B. Schrøder

Thomas B. Schrøder appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.softcond-mat.stat-mechphysics.chem-phcond-mat.dis-nncond-mat.mtrl-sciphysics.comp-phcond-mat.otherBiological Physicscond-mat.mes-hallphysics.plasm-ph

Catalog footprint

What is connected

41works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

41 published item(s)

preprint2022arXiv

A toy model for viscous liquid dynamics

A simple model for viscous liquid dynamics is introduced. Consider the surface of the union of hyper-spheres centered at random positions inside a hypercube with periodic boundary conditions. It is argued and demonstrated by numerical simulations that at high dimensions geodetic flows on this surface is a good model for viscous liquid dynamics. It is shown that this simple model exhibits viscous dynamics for densities above the percolation threshold in $8$, $12$ and $16$ dimensions. Thus the slowing down of the dynamics, measured by the mean-squared displacement, extends to several orders of magnitude similarly to what is observed in other models for viscous dynamics. Furthermore, the shape of the mean-squared displacement is to a very good approximation the same as for the standard model in simulations of viscous liquids: the Kob-Andersen binary Lennard Jones mixture.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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).

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access