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

Ulf R. Pedersen

Ulf R. Pedersen 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-mechcond-mat.mtrl-scicond-mat.dis-nnphysics.chem-phBiological Physicsphysics.comp-phphysics.geo-ph

Catalog footprint

What is connected

18works
8topics
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

18 published item(s)

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.

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

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Computing the crystal growth rate by the interface pinning method

An essential parameter for crystal growth is the kinetic coefficient given by the proportionality between super-cooling and average growth velocity. Here we show that this coefficient can be computed in a single equilibrium simulation using the interface pinning method where two-phase configurations are stabilized by adding an spring-like bias field coupling to an order-parameter that discriminates between the two phases. Crystal growth is a Smoluchowski process and the crystal growth rate can therefore be computed from the terminal exponential relaxation of the order parameter. The approach is investigated in detail for the Lennard-Jones model. We find that the kinetic coefficient scales as the inverse square-root of temperature along the high temperature part of the melting line. The practical usability of the method is demonstrated by computing the kinetic coefficient of the elements Na, Mg, Al and Si from first principles. It is briefly discussed how a generalized version of the method is an alternative to forward flux sampling methods for computing rates along trajectories of rare events.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Computing Gibbs free energy differences by interface pinning

We propose an approach for computing the Gibbs free energy difference between phases of a material. The method is based on the determination of the average force acting on interfaces that separate the two phases of interest. This force, which depends on the Gibbs free energy difference between the phases, is computed by applying an external harmonic field that couples to a parameter which specifies the two phases. Validated first for the Lennard-Jones model, we demonstrate the flexibility, efficiency and practical applicability of this approach by computing the melting temperatures of sodium, magnesium, aluminum and silicon at ambient pressure using density functional theory. Excellent agreement with experiment is found for all four elements, except for silicon, for which the melting temperature is, in agreement with previous simulations, seriously underestimated.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Direct calculation of the solid-liquid Gibbs free energy difference in a single equilibrium simulation

Computing phase diagrams of model systems is an essential part of computational condensed matter physics. In this paper we discuss in detail the interface pinning (IP) method for calculation of the Gibbs free energy difference between a solid and a liquid. This is done in a single equilibrium simulation by applying a harmonic field that biases the system towards two-phase configurations. The Gibbs free energy difference between the phases is determined from the average force that the applied field exerts on the system. As a test system we study the Lennard-Jones model. It is shown that the coexistence line can be computed efficiently to a high precision when the IP method is combined with the Newton-Raphson method for finding roots. Statistical and systematic errors are investigated. Advantages and drawbacks of the IP method are discussed. The high pressure part of the temperature-density coexistence region is outlined by isomorphs.

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
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
preprint2010arXiv

Crystallization of the Lewis-Wahnström ortho-terphenyl model

Crystallization is observed during long molecular dynamics simulations of bent trimers, a molecular model proposed by Lewis and Wahnström for ortho-terphenyl. In the crystal, the three spheres that make up the rigid molecule sit near sites of a body centered cubic lattice, the trimer bond angle being almost optimal for this structure. The crystal exhibits orientational disorder with the molecules aligned randomly along the three Cartesian axis (an example of cubatic orientational order). The rotational and translational mobilities exhibit only modest decreases on crystallization, by factors of 10 and 3 respectively. The rotational relaxation does change from Debye-like in the liquid to large angle jumps in the crystal. We consider the origin of the superior glass forming ability of the the trimer over the monatomic liquid.

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

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

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