General features of the energy landscape in Lennard-Jones-like model liquids

Heterogeneous-elasticity theory of instantaneous normal modes in liquids

Since decades, the concept of vibrational density of states in glasses has been mirrored in liquids by the instantaneous-normal-mode spectrum. In glasses instantaneous configurations are believed to be situated close to minima of the potential-energy hypersurface and all eigenvalues of the associated Hessian matrix are positive. In liquids this is no longer true, and modes corresponding to both positive and negative eigenvalues exist. The instantaneous-normal-mode spectrum has been numerically investigated in the past, and it has been demonstrated to bring important information on the liquid dynamics and transport properties. A systematic deeper theoretical understanding is now needed. Heterogeneous-elasticity theory has proven to be particularly successful in explaining many details of the low-frequency excitations in glasses, ranging from the thoroughly studied boson peak, to other anomalies related to the crossover between wave-like and random-matrix-like excitations. Here we present an extension of heterogeneous-elasticity theory to the liquid state, and show that the outcome of the theory agrees well to the results of extensive molecular-dynamics simulations of a model liquid at different temperatures. We find that the spectrum of eigenvalues [Formula: see text] has a sharp maximum close to (but not at) [Formula: see text], and decreases monotonically with [Formula: see text] on both its stable and unstable side. We show that the spectral shape strongly depends on temperature, being symmetric at high temperatures and becoming rather asymmetric at low temperatures, close to the dynamical critical temperature. Most importantly, we demonstrate that the theory naturally reproduces a surprising phenomenon, a zero-energy spectral singularity with a cusp-like character developing in the vibrational spectra upon cooling. This feature, known from a few previous numerical studies, has been generally overlooked in the past due to a misleading representation of the data. We provide a thorough analysis of this issue, based on both very accurate predictions of our theory, and computational studies of model liquid systems with extended size.

Enhanced vibrational stability in glass droplets

We show through simulations of amorphous solids prepared in open-boundary conditions that they possess significantly fewer low-frequency vibrational modes compared to their periodic boundary counterparts. Specifically, using measurements of the vibrational density of states, we find that the D ( ω ) ∼ ω 4 law changes to D ( ω ) ∼ ω δ with δ ≈ 5 in two dimensions and δ ≈ 4.5 in three dimensions. Crucially, this enhanced stability is achieved when utilizing slow annealing protocols to generate solid configurations. We perform an anharmonic analysis of the minima corresponding to the lowest frequency modes in such open-boundary systems and discuss their correlation with the density of states. A study of various system sizes further reveals that small systems display a higher degree of localization in vibrations. Lastly, we confine open-boundary solids in order to introduce macroscopic stresses in the system, which are absent in the unconfined system and find that the D ( ω ) ∼ ω 4 behavior is recovered.

Modeling the instantaneous normal mode spectra of liquids as that of unstable elastic media

We present a theory, based on an unstable elasticity model with thermally fluctuating elastic constants, which explains the salient features of the instantaneous normal mode (INM) spectrum of a numerically simulated liquid over a large range of temperatures. The INM spectrum of a liquid is obtained in a molecular dynamics (MD) simulation by considering the harmonic eigenvalue spectrum of the potential energy, taken at a given instant (snapshot). Because the INM spectrum records the curvatures of the fluctuating potential landscape of the liquid, an understanding of this spectrum may pave the way toward understanding the liquid to glass transformation process.

We study the instantaneous normal mode (INM) spectrum of a simulated soft-sphere liquid at different equilibrium temperatures T. We find that the spectrum of eigenvalues ρ(λ) has a sharp maximum near (but not at) λ=0 and decreases monotonically with |λ| on both the stable and unstable sides of the spectrum. The spectral shape strongly depends on temperature. It is rather asymmetric at low temperatures (close to the dynamical critical temperature) and becomes symmetric at high temperatures. To explain these findings we present a mean-field theory for ρ(λ), which is based on a heterogeneous elasticity model, in which the local shear moduli exhibit spatial fluctuations, including negative values. We find good agreement between the simulation data and the model calculations, done with the help of the self-consistent Born approximation (SCBA), when we take the variance of the fluctuations to be proportional to the temperature T. More importantly, we find an empirical correlation of the positions of the maxima of ρ(λ) with the low-frequency exponent of the density of the vibrational modes of the glasses obtained by quenching to T=0 from the temperature T. We discuss the present findings in connection to the liquid to glass transformation and its precursor phenomena.

Theoretical Estimate of the Glass Transition Line of Yukawa One-Component Plasmas

The mode coupling theory of supercooled liquids is combined with advanced closures to the integral equation theory of liquids in order to estimate the glass transition line of Yukawa one-component plasmas from the unscreened Coulomb limit up to the strong screening regime. The present predictions constitute a major improvement over the current literature predictions. The calculations confirm the validity of an existing analytical parameterization of the glass transition line. It is verified that the glass transition line is an approximate isomorphic curve and the value of the corresponding reduced excess entropy is estimated. Capitalizing on the isomorphic nature of the glass transition line, two structural vitrification indicators are identified that allow a rough estimate of the glass transition point only through simple curve metrics of the static properties of supercooled liquids. The vitrification indicators are demonstrated to be quasi-universal by an investigation of hard sphere and inverse power law supercooled liquids. The straightforward extension of the present results to bi-Yukawa systems is also discussed.

Barriers, trapping times, and overlaps between local minima in the dynamics of the disordered Ising p-spin model

We study the low-temperature out-of-equilibrium Monte Carlo dynamics of the disordered Ising p-spin Model with p=3 and a small number of spin variables. We focus on sequences of configurations that are stable against single spin flips obtained by instantaneous gradient descent from persistent ones. We analyze the statistics of energy gaps, energy barriers, and trapping times on subsequences such that the overlap between consecutive configurations does not overcome a threshold. We compare our results to the predictions of various trap models finding the best agreement with the step model when the p-spin configurations are constrained to be uncorrelated.

Single molecule probing of dynamics in supercooled polymers

Fluorescence experiments with single BODIPY molecules embedded in a poly(methyl acrylate) matrix have been performed at various temperatures in the supercooled regime. By using pulsed excitation, fluorescence lifetime and linear dichroism time trajectories were accessible at the same time. Both observables have been analyzed without data binning. While the linear dichroism solely reflects single particle dynamics, the fluorescence lifetime observable depends on the molecular environment, so that the dynamics from the polymer host surrounding a chromophore contributes to this quantity. We observe that the lifetime correlation decays slightly faster than polarization correlation, indicating the occurrence of large angular reorientations. Additionally, dichroism time trajectories have been adducted to reveal directly the geometry of rotational dynamics. We identify small but also significantly larger rotational jumps being responsible for the overall molecular reorientation.

Length scales and dynamics in the reorientational relaxation of tracers in molecular and polymeric glass formers via electron spin resonance spectroscopy

It is widely accepted that a temperature region exists above the glass transition temperature, playing a fundamental role in the physics of polymers and glass formers. In this region, several dynamic crossovers have been experimentally and numerically revealed in the past years and the onset of glassy behavior is generally located, because of cooperative and heterogeneous features exhibited by the dynamics. In this Article, the rotational dynamics of two different stiff molecular spin probes dissolved in poly(propylene glycol) has been investigated by electron spin resonance spectroscopy in a wide temperature range. Decoupling phenomena between macroscopic and microscopic transport properties have been observed in the crossover region. A comparison with previous studies carried out with the same tracers dissolved in different polymers and glass formers strongly supports the idea that the observed crossover signals the onset of spatial correlations of dynamics on the length scale probed by the tracers.

Dynamics and energy landscape in a tetrahedral network glass-former: direct comparison with models of fragile liquids

We report molecular dynamics simulations for a new model of tetrahedral network glass-former, based on short-range spherical potentials. Despite the simplicity of the forcefield employed, our model reproduces some essential physical properties of silica, an archetypal network-forming material. Structural and dynamical properties, including dynamic heterogeneities and the nature of local rearrangements, are investigated in detail and a direct comparison with models of close-packed, fragile glass-formers is performed. The outcome of this comparison is rationalized in terms of the properties of the potential energy surface, focusing on the unstable modes of the stationary points. Our results indicate that the weak degree of dynamic heterogeneity observed in network glass-formers may be attributed to an excess of localized unstable modes, associated with elementary dynamical events such as bond breaking and reformation. In contrast, the more fragile Lennard-Jones mixtures are characterized by a larger fraction of extended unstable modes, which lead to a more cooperative and heterogeneous dynamics.

Density scaling of the diffusion coefficient at various pressures in viscous liquids

Fundamental thermodynamics and an earlier elastic solid-state point defect model [P. Varotsos and K. Alexopoulos, Phys. Rev. B 15, 4111 (1977); Phys. Rev. B 18, 2683 (1978)] are employed to formulate an analytical second-order polynomial function describing the density scaling of the diffusion coefficient in viscous liquids. The function parameters are merely determined by the scaling exponent, which is directly connected with the Grüneisen constant. Density scaling diffusion coefficient isotherms obtained at different pressures collapse on a unique master curve, in agreement with recent computer simulation results of Lennard-Jones viscous liquids [D. Coslovich and C. M. Roland, J. Phys. Chem. 112, 1329 (2008)].

Exploring the potential energy landscape of glass-forming systems: from inherent structures via metabasins to macroscopic transport

In this review a systematic analysis of the potential energy landscape (PEL) of glass-forming systems is presented. Starting from the thermodynamics, the route towards the dynamics is elucidated. A key step in this endeavor is the concept of metabasins. The relevant energy scales of the PEL can be characterized. Based on the simulation results for some glass-forming systems one can formulate a relevant model system (ideal Gaussian glass-former) which can be treated analytically. The macroscopic transport can be related to the microscopic hopping processes, using either the strong relation between energy (thermodynamics) and waiting times (dynamics) or, alternatively, the concepts of the continuous-time random walk. The relation to the geometric properties of the PEL is stressed. The emergence of length scales within the PEL approach as well as the nature of finite-size effects is discussed. Furthermore, the PEL view is compared to other approaches describing the glass transition.

Understanding fragility in supercooled Lennard-Jones mixtures. II. Potential energy surface

The connection between isobaric fragility and the properties of high-order stationary points of the potential energy surface in different supercooled Lennard-Jones mixtures was investigated. The increase of effective activation energies upon supercooling appears to be driven by the increase of average potential energy barriers measured by the energy dependence of the fraction of unstable modes. Such an increase is sharper, the more fragile the mixture. Correlations between fragility and other properties of high-order stationary points, including the vibrational density of states and the localization features of unstable modes, are also discussed.

Fickian crossover and length scales from two point functions in supercooled liquids

Particle motion of a Lennard-Jones supercooled liquid near the glass transition is studied by molecular dynamics simulations. We analyze the wave vector dependence of relaxation times in the incoherent self-scattering function and show that at least three different regimes can be identified and its scaling properties determined. The transition from one regime to another happens at characteristic length scales. The length scale associated with the onset of Fickian diffusion corresponds to the maximum size of heterogeneities in the system, and the characteristic time scale is several times larger than the alpha relaxation time. A second crossover length scale is observed, which corresponds to the typical time and length of heterogeneities, in agreement with results from four point functions. The different regimes can be traced back to the behavior of the van Hove distribution of displacements, which shows a characteristic exponential regime in the heterogeneous region before the crossover to Gaussian diffusion and should be observable in experiments. Our results show that it is possible to obtain characteristic length scales of heterogeneities through the computation of two point functions at different times.

Out-of-equilibrium dynamical fluctuations in glassy systems

In this paper we extend the earlier treatment of out-of-equilibrium mesoscopic fluctuations in glassy systems in several significant ways. First, via extensive simulations, we demonstrate that models of glassy behavior without quenched disorder display scalings of the probability of local two-time correlators that are qualitatively similar to that of models with short-ranged quenched interactions. The key ingredient for such scaling properties is shown to be the development of a criticallike dynamical correlation length, and not other microscopic details. This robust data collapse may be described in terms of a time-evolving "extreme value" distribution. We develop a theory to describe both the form and evolution of these distributions based on a effective sigma model approach.

On the mechanism of reorientational and structural relaxation in supercooled liquids: the role of border dynamics and cooperativity

Molecular dynamics simulation and analysis based upon the many-body potential energy landscape (PEL) are employed to characterize single molecule reorientation and structural relaxation, and their interrelation, in deeply supercooled liquid CS(2). The rotational mechanism changes from small-step Debye diffusion to sudden large angle reorientation (SLAR) as the temperature falls below the mode-coupling temperature T(c). The onset of SLAR is explained in terms of the PEL; it is an essential feature of low-T rotational dynamics, along with the related phenomena of dynamic heterogeneity and the bifurcation of slow and fast relaxation processes. A long trajectory in which the system is initially trapped in a low energy local minimum, and eventually escapes, is followed in detail, both on the PEL and in real space. During the trapped period, "return" dynamics occurs, always leading back to the trap. Structural relaxation is identified with irreversible escape to a new trap. These processes lead to weak and strong SLAR, respectively; strong SLAR is a clear signal of structural relaxation. Return dynamics involves small groups of two to four molecules, while a string-like structure composed of all the active groups participates in the escape. It is proposed that, rather than simple, nearly instantaneous, one-dimensional barrier crossings, relaxation involves activation of the system to the complex, multidimensional region on the borders of the basins of attraction of the minima for an extended period.

The glass state: characterizing energy landscape features

We probed the topologies imposed on configuration hyperspace by the potential energy function-the shapes of the constant potential energy manifolds-for the glassy state of monatomic Lennard-Jones matter, by following trajectories of constant potential energy. A prominent characteristic of this model matter is well-defined regions of confinement (pockets) in configuration hyperspace. We found that there are constant potential energy hyperspace paths (tubes) between such pockets, applying even to paths linking glassy regions to crystalline regions of configuration hyperspace. Also, we found that glass and crystal pockets are interspersed. For monatomic Lennard-Jones matter at least, the transition from glass to crystal therefore does not have to involve the traversing of a potential energy barrier, as is usually assumed.

Non-Gaussian energy landscape of a simple model for strong network-forming liquids: Accurate evaluation of the configurational entropy

We present a numerical study of the statistical properties of the potential energy landscape of a simple model for strong network-forming liquids. The model is a system of spherical particles interacting through a square-well potential, with an additional constraint that limits the maximum number of bonds Nmax per particle. Extensive simulations have been carried out as a function of temperature, packing fraction, and Nmax. The dynamics of this model are characterized by Arrhenius temperature dependence of the transport coefficients and by nearly exponential relaxation of dynamic correlators, i.e., features defining strong glass-forming liquids. This model has two important features: (i) Landscape basins can be associated with bonding patterns. (ii) The configurational volume of the basin can be evaluated in a formally exact way, and numerically with an arbitrary precision. These features allow us to evaluate the number of different topologies the bonding pattern can adopt. We find that the number of fully bonded configurations, i.e., configurations in which all particles are bonded to Nmax neighbors, is extensive, suggesting that the configurational entropy of the low temperature fluid is finite. We also evaluate the energy dependence of the configurational entropy close to the fully bonded state and show that it follows a logarithmic functional form, different from the quadratic dependence characterizing fragile liquids. We suggest that the presence of a discrete energy scale, provided by the particle bonds, and the intrinsic degeneracy of fully bonded disordered networks differentiates strong from fragile behavior.

Relationship between phase transitions and topological changes in one-dimensional models

We address the question of the quantitative relationship between thermodynamic phase transitions and topological changes in the potential energy manifold analyzing two classes of one dimensional models, the Burkhardt solid-on-solid model and the Peyrard-Bishop model for DNA thermal denaturation, both in the confining and nonconfining version. These models, apparently, do not fit [M. Kastner, Phys. Rev. Lett. 93, 150601 (2004)] in the general idea that the phase transition is signaled by a topological discontinuity. We show that in both models the phase transition energy v(c) is actually noncoincident with, and always higher than, the energy v(theta) at which a topological change appears. However, applying a procedure already successfully employed in other cases as the mean field phi4 model, i.e., introducing a map M:v-->v(s) from levels of the energy hypersurface V to the level of the stationary points "visited" at temperature T, we find that M (v(c))=v(theta). This result enhances the relevance of the underlying stationary points in determining the thermodynamics of a system, and extends the validity of the topological approach to the study of phase transition to the elusive one-dimensional systems considered here.

Molecular dynamics studies of heterogeneous dynamics and dynamic crossover in supercooled atomic liquids

Supercooled liquids near the glass transition exhibit the phenomenon of heterogeneous relaxation; at any specific time, a nominally homogeneous equilibrium fluid undergoes dynamic fluctuations in its structure on a molecular distance scale with rates that are very different in different regions of the sample. Several theoretical and simulation studies have suggested a change in the nature of the dynamics of fluids as they are supercooled, leading to the concept of a dynamic crossover that is often associated with mode coupling theory. Here, we will review the use of molecular dynamics computer simulation methods to investigate heterogeneous dynamics and dynamic crossovers in models of atomic liquids.

Topology and phase transitions: from an exactly solvable model to a relation between topology and thermodynamics

The elsewhere surmized topological origin of phase transitions is given here important evidence through the analytic study of an exactly solvable model for which both topology of submanifolds of configuration space and thermodynamics are worked out. The model is a mean-field one with a k-body interaction. It undergoes a second-order phase transition for k=2 and a first-order one for k >2 . This opens a perspective for the understanding of the deep origin of first and second-order phase transitions, respectively. In particular, a remarkable theoretical result consists of a mathematical characterization of first-order transitions. Moreover, we show that a "reduced" configuration space can be defined in terms of collective variables, such that the correspondence between phase transitions and topology changes becomes one-to-one, for this model. Finally, an unusual relationship is worked out between the microscopic description of a classical N -body system and its macroscopic thermodynamic behavior. This consists of a functional dependence of thermodynamic entropy upon the Morse indexes of the critical points (saddles) of the constant energy hypersurfaces of the microscopic 2N-dimensional phase space. Thus phase space (and configuration space) topology is directly related to thermodynamics.

Dynamic phase diagram of the number partitioning problem

We study the dynamic phase diagram of a spin model associated with the number partitioning problem, as a function of temperature and of the fraction K/N of spins allowed to flip simultaneously. The case K=1 reproduces the activated behavior of Bouchaud's trap model, whereas the opposite limit K=N can be mapped onto the entropic trap model proposed by Barrat and Mézard. In the intermediate case 1<<K<<N , the dynamics corresponds to a modified version of the Barrat and Mézard model, which includes a slow (rather than instantaneous) decorrelation at each step. A transition from an activated regime to an entropic one is observed at temperature T(g) /2 in agreement with recent work on this model. Ergodicity breaking occurs for T< T(g) /2 in the thermodynamic limit, if K/N-->0 . In this temperature range, the model exhibits a nontrivial fluctuation-dissipation relation leading for K<<N to a single effective temperature equal to T(g) /2 . These results give insights into the relevance and limitations of the picture proposed by simple trap models.

Saddle index properties, singular topology, and its relation to thermodynamic singularities for a phi4 mean-field model

We investigate the potential energy surface of a phi(4) model with infinite range interactions. All stationary points can be uniquely characterized by three real numbers alpha(+) , alpha(0) , alpha(-) with alpha(+) + alpha(0) + alpha(-) =1 , provided that the interaction strength mu is smaller than a critical value. The saddle index n(s) is equal to alpha(0) and its distribution function has a maximum at n(max)(s) =1/3 . The density p (e) of stationary points with energy per particle e , as well as the Euler characteristic chi (e) , are singular at a critical energy e(c) (mu) , if the external field H is zero. However, e(c) (mu) not equal upsilon(c) (mu) , where upsilon(c) (mu) is the mean potential energy per particle at the thermodynamic phase transition point T(c) . This proves that previous claims that the topological and thermodynamic transition points coincide is not valid, in general. Both types of singularities disappear for H not equal 0 . The average saddle index n (s) as function of e decreases monotonically with e and vanishes at the ground state energy, only. In contrast, the saddle index n(s) as function of the average energy e ( n(s) ) is given by n(s) ( e ) =1+4 e (for H=0 ) that vanishes at e =-1/4> upsilon(0) , the ground state energy.

Saddles in the energy landscape: extensivity and thermodynamic formalism

We formally extend the energy landscape approach for the thermodynamics of liquids to account for saddle points. By considering the extensive nature of macroscopic potential energies, we derive the scaling behavior of saddles with system size, as well as several approximations for the properties of low-order saddles (i.e., those with only a few unstable directions). We then cast the canonical partition function in a saddle-explicit form and develop, for the first time, a rigorous energy landscape approach capable of reproducing trends observed in simulations, in particular, the temperature dependence of the energy and fractional order of sampled saddles.

Real space origin of temperature crossovers in supercooled liquids

We show that the various crossovers between dynamical regimes observed in experiments and simulations of supercooled liquids can be explained in simple terms from the existence and statistical properties of dynamical heterogeneities. We confirm that dynamic heterogeneity is responsible for the slowing down of glass formers at temperatures well above the dynamic singularity Tc predicted by mode-coupling theory. Our results imply that activated processes govern the long-time dynamics even in the temperature regime where they are neglected by mode-coupling theory. We show that alternative interpretations based on topographic properties of the potential energy landscape are inefficient ways of describing simple physical features which are naturally accounted for within our approach. We show in particular that the reported links between mode coupling and landscape singularities do not exist.