Source author record

Srikanth Sastry

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

39works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

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

39 published item(s)

preprint2025arXiv

Self-organization, Memory and Learning: From Driven Disordered Systems to Living Matter

Disordered systems subject to a fluctuating environment can self-organize into a complex history-dependent response, retaining a memory of the driving. In sheared amorphous solids, self-organization is established by the emergence of a persistent system of mechanical instabilities that can repeatedly be triggered by the driving, leading to a state of high mechanical reversibility. As a result of self-organization, the response of the system becomes correlated with the dynamics of its environment, which can be viewed as a sensing mechanism of the system's environment. Such phenomena emerge across a wide variety of soft matter systems, suggesting that they are generic and hence may depend very little on the underlying specifics. We review self-organization in driven amorphous solids, concluding with a discussion of what self-organization in driven disordered systems can teach us about how simple organisms sense and adapt to their changing environments.

preprint2022arXiv

Criticality and marginal stability of the shear jamming transition of frictionless soft spheres

We study numerically the critical behavior and marginal stability of the shear jamming transition for frictionless soft spheres, observed to occur over a finite range of densities, associated with isotropic jamming for densities above the minimum jamming (J-point) density. Several quantities are shown to scale near the shear jamming point in the same way as the isotropic jamming point. We compute the exponents associated with the small force distribution and the interparticle gap distribution,and show that the corresponding exponents are consistent with the marginal stability condition observed for isotropic jamming, and with predictions of the mean field theory of jamming in hard spheres.

preprint2022arXiv

Crossover in dynamics in the Kob-Andersen binary mixture glass-forming liquid

Glass-forming liquids are broadly classified as being fragile or strong, depending on the deviation from Arrhenius behavior of their relaxation times. A fragile to strong crossover is observed or inferred in liquids like water and silica, and more recently also in metallic glasses and phase change alloys, leading to the expectation that such a crossover is more widely realised among glass formers. We investigate computationally the well-studied Kob-Andersen model, accessing temperatures well below the mode coupling temperature $T_{MCT}$. We find that relaxation times exhibit a crossover in dynamics around $T_{MCT}$, and discuss whether it bears characteristics of the fragile to strong crossover. Several aspects of dynamical heterogeneity exhibit behavior mirroring the dynamical crossover, whereas thermodynamic quantities do not. In particular, the Adam-Gibbs relation describing the relation between relaxation times and configurational entropy continues to hold below the dynamical crossover, when anharmonic corrections to the vibrational entropy are included.

preprint2022arXiv

Dependence of the glass transition and jamming densities on spatial dimension

We investigate the dynamics of soft sphere liquids through computer simulations for spatial dimensions from $d =3$ to $8$, over a wide range of temperatures and densities. Employing a scaling of density-temperature dependent relaxation times, we precisely identify the density $ϕ_0$ which marks the ideal glass transition in the hard sphere limit, and a crossover from sub- to super-Arrhenius temperature dependence. The difference between $ϕ_0$ and the athermal jamming density $ϕ_J$, small in 3 and 4 dimensions, increases with dimension, with $ϕ_0 > ϕ_J$ for $d > 4$. We compare our results with recent theoretical calculations.

preprint2022arXiv

Enhancing associative memory recall in non-equilibrium materials through activity

Associative memory, a form of content-addressable memory, facilitates information storage and retrieval in many biological and physical systems. In statistical mechanics models, associative memory at equilibrium is represented through attractor basins in the free energy landscape. Here, we use the Hopfield model, a paradigmatic model to describe associate memory, to investigate the effect of non-equilibrium activity on memory retention and recall. We introduce activity into the system as gaussian-colored noise which breaks detailed balance and forces the system out of equilibrium. We observe that, under these non-equilibrium conditions, the Hopfield network has a higher storage capacity than that allowed at equilibrium. Using analytical and numerical techniques, we show that the rate of entropy production modifies the energy landscape and helps the system to access memory regions which were previously inaccessible.

preprint2022arXiv

Liquid-liquid phase transition in deeply supercooled Stillinger-Weber silicon

The existence of a phase transition between two distinct liquid phases in single-component network-forming liquids (e.g., water, silica, silicon) has elicited considerable scientific interest. The challenge, both for experiments and simulations, is that the liquid-liquid phase transition occurs under deeply supercooled conditions, where crystallization occurs very rapidly. Thus, early evidence from numerical equation of state studies was challenged, with the argument that slow spontaneous crystallization had been misinterpreted as evidence of a second liquid state. Rigorous free energy calculations have subsequently confirmed the existence of a liquid-liquid phase transition in some models of water, and exciting new experimental evidence has since supported these computational results. Similar results have so far not been found for silicon. Here, we present results from free energy calculations performed for silicon modelled with the classical, empirical Stillinger-Weber potential. Through a careful study employing state-of-the-art constrained simulation protocols and numerous checks for thermodynamic consistency, we find that there are two distinct metastable liquid states and a phase transition. Our results resolve a long-standing debate concerning the existence of a liquid-liquid transition in supercooled liquid silicon and address key questions regarding the nature of the phase transition and the associated critical point.

preprint2022arXiv

Mean field theory of yielding under oscillatory shear

We study a mean field elastoplastic model, embedded within a disordered landscape of local yield barriers, to shed light on the behaviour of athermal amorphous solids subject to oscillatory shear. We show that the model presents a genuine dynamical transition between an elastic and a yielded state, and qualitatively reproduces the dependence on the initial degree of annealing found in particle simulations. For initial conditions prepared below the analytically derived threshold energy, we observe a non-trivial, non-monotonic approach to the yielded state. The timescale diverges as one approaches the yielding point from above, which we identify with the fatigue limit. We finally discuss the connections to brittle yielding under uniform shear.

preprint2022arXiv

Yielding behavior of glasses under asymmetric cyclic deformation

We consider the yielding behaviour of a model glass subjected to asymmetric cyclic shear deformation, wherein the applied strain varies between 0 and a maximum value $γ_{\rm max}$, and study its dependence on the degree of annealing of the glass and system size. The yielding behaviour of well annealed glasses (unlike poorly annealed glasses) display striking differences from the symmetric case, with the emergence of an intermediate strain regime with substantial plasticity but no yielding. The observed behaviour is satisfactorily captured by a recently proposed model. For larger system sizes, the intermediate strain regime narrows, leading to a remarkable reversal of yield strain with annealing.

preprint2021arXiv

Avalanches and Structural Change in Cyclically Sheared Silica Glass

We investigate avalanches associated with plastic rearrangements and the nature of structural change in the prototypical strong glass, silica, computationally. Although qualitative aspects of yielding in silica are similar to other glasses, we find that the statistics of avalanches exhibits non-trivial behaviour. Investigating the statistics of avalanches and clusters in detail, we propose and verify a new relation between exponents characterizing the size distribution of avalanches and clusters. Across the yielding transition, anomalous structural change and densification, associated with a suppression of tetrahedral order, is observed to accompany strain localisation.

preprint2021arXiv

Correlation between plastic rearrangements and local structure in a cyclically driven glass

The correlation between local structure and the propensity for structural rearrangements has been widely investigated in glass forming liquids and glasses. In this paper we use the excess two-body entropy $S_2$ and tetrahedrality $\n_{tet}$ as the per-particle local structural order parameters to explore such correlations in a three-dimensional model glass subjected to cyclic shear deformation. We first show that for both liquid configurations and the corresponding inherent structures, local ordering increases upon lowering temperature, signaled by a decrease in the two-body entropy and an increase in tetrahedrality. When the inherent structures, or glasses, are periodically sheared athermally, they eventually reach absorbing states for small shear amplitudes, which do not change from one cycle to the next. Large strain amplitudes result in the formation of shear bands, within which particle motion is diffusive. We show that in the steady state, there is a clear difference in the local structural environment of particles that will be part of plastic rearrangements during the next shear cycle and that of particles that are immobile. In particular, particles with higher $S_2$ and lower $n_{tet}$ are more likely to go through rearrangements irrespective of the average energies of the configurations and strain amplitude. For high shear, we find very distinctive local order outside the mobile shear band region, where almost $30\%$ of the particles are involved in icosahedral clusters, contrasting strongly with the fraction of $<5\%$ found inside the shear band.

preprint2021arXiv

Temperature Protocols to Guide Selective Self-Assembly of Competing Structures

Multi-component self-assembly mixtures offer the possibility of encoding multiple target structures with the same set of interacting components. Selective retrieval of one of the stored structures has been attempted by preparing an initial state that favours the assembly of the required target, through seeding, concentration patterning or specific choices of interaction strengths. This may not be possible in an experiment where on-the-fly reconfiguration of the building blocks to switch functionality may be required. In this paper, we explore principles of inverse design of a multi-component self-assembly mixture capable of encoding two competing structures that can be selected through simple temperature protocols. We design the target structures to realise the generic situation in which one of targets has the lower nucleation barrier while the other is globally more stable. We observe that to avoid the formation of spurious or chimeric aggregates, the number of neighbouring component pairs that occur in both structures should be minimal. Our design also requires the inclusion of components that are part only of one of the target structures, but we observe that to maximize the selectivity of retrieval, the component library itself should be maximally shared by the two targets. We demonstrate that temperature protocols can be designed which lead to the formation of either one of the target structures with high selectivity. We discuss the important role played by secondary aggregation products, which we term vestigial aggregates.

preprint2021arXiv

Yielding transition of a two dimensional glass former under athermal cyclic shear deformation

We study numerically the yielding transition of a two dimensional model glass subjected to athermal quasi-static cyclic shear deformation, with the aim of investigating the effect on the yielding behaviour of the degree of annealing, which in turn depends on the preparation protocol. We find two distinct regimes of annealing separated by a threshold energy. Poorly annealed glasses progressively evolve towards the threshold energy as the strain amplitude is increased towards the yielding value. Well annealed glasses with initial energies below the threshold energy exhibit stable behaviour, with negligible change in energy with increasing strain amplitude, till they yield. Discontinuities in energy and stress at yielding increase with the degree of annealing, consistently with recent results found in three dimensions. We observe significant structural change with strain amplitude that closely mirrors the changes in energy and stresses. We investigate groups of particles that are involved in plastic rearrangements. We analyse the distributions of avalanche sizes, of clusters of connected rearranging particles, and related quantities, employing finite size scaling analysis. We verify previously investigated relations between exponents characterising these distributions, and a newly proposed relation between exponents describing avalanche and cluster size distributions.

preprint2020arXiv

Is friction essential for dilatancy and shear jamming in granular matter?

Granular packings display the remarkable phenomenon of dilatancy [1], wherein their volume increases upon shear deformation. Conventional wisdom and previous results suggest that dilatancy, as also the related phenomenon of shear-induced jamming, requires frictional interactions [2, 3]. Here, we investigate the occurrence of dilatancy and shear jamming in frictionless packings. We show that the existence of isotropic jamming densities ϕj above the minimal density, the J-point density ϕJ [4, 5], leads both to the emergence of shear-induced jamming and dilatancy. Packings at ϕJ form a significant threshold state into which systems evolve in the limit of vanishing pressure under constant pressure shear, irrespective of the initial jamming density ϕj. While packings for different ϕj display equivalent scaling properties under compression [6], they exhibit striking differences in rheological behaviour under shear. The yield stress under constant volume shear increases discontinuously with density when ϕj > ϕJ, contrary to the continuous behavior in generic packings that jam at ϕJ [4, 7].

preprint2020arXiv

Timescale divergence at the shear jamming transition

We find that in simulations of quasi-statically sheared frictional disks, the shear jamming transition can be characterized by an abrupt jump in the number of force bearing contacts between particles. This mechanical coordination number increases discontinuously from $Z = 0$ to $Z \gtrsim d +1$ at a critical shear value $γ_c$, as opposed to a smooth increase in the number of geometric contacts. This is accompanied by a diverging timescale $τ^*$ that characterizes the time required by the system to attain force balance when subjected to a perturbation. As the global shear $γ$ approaches the critical value $γ_c$ from below, one observes the divergence of the time taken to relax to a state where all the inter-particle contacts have uniformly zero force. Above $γ_{c}$, the system settles into a state characterized by finite forces between particles, with the timescale also increasing as $γ\to γ_{c}^{+}$. By using two different protocols to generate force balanced configurations, we show that this timescale divergence is a robust feature that accompanies the shear jamming transition.

preprint2019arXiv

Information density, structure and entropy in equilibrium and non-equilibrium systems

During a spontaneous change, a macroscopic physical system will evolve towards a macro-state with more realizations. This observation is at the basis of the Statistical Mechanical version of the Second Law of Thermodynamics, and it provides an interpretation of entropy in terms of probabilities. However, we cannot rely on the statistical-mechanical expressions for entropy in systems that are far from equilibrium. In this paper, we compare various extensions of the definition of entropy, which have been proposed for non-equilibrium systems. It has recently been proposed that measures of information density may serve to quantify entropy in both equilibrium and nonequilibrium systems. We propose a new "bit-wise" method to measure the information density for off lattice systems. This method does not rely on coarse-graining of the particle coordinates. We then compare different estimates of the system entropy, based on information density and on the structural properties of the system, and check if the various entropies are mutually consistent and, importantly, whether they can detect non-trivial ordering phenomena. We find that, except for simple (one-dimensional) cases, the different methods yield answers that are at best qualitatively similar, and often not even that, although in several cases, different entropy estimates do detect ordering phenomena qualitatively. Our entropy estimates based on bit-wise data compression contain no adjustable scaling factor, and show large quantitative differences with the thermodynamic entropy obtained from equilibrium simulations. Hence, our results suggest that, at present, there is not yet a single, structure-based entropy definition that has general validity for equilibrium and non equilibrium systems.

preprint2019arXiv

Unified phase diagram of reversible-irreversible, jamming and yielding transitions in cyclically sheared soft sphere packings

Self-organization, and transitions from reversible to irreversible behaviour, of interacting particle assemblies driven by externally imposed stresses or deformation is of interest in comprehending diverse phenomena in soft matter. They have been investigated in a wide range of systems, such as colloidal suspensions, glasses, and granular matter. In different density and driving regimes, such behaviour is related to yielding of amorphous solids, jamming, and memory formation, \emph{etc.} How these phenomena are related to each other has not, however, been much studied. In order to obtain a unified view of the different regimes of behaviour, and transitions between them, we investigate computationally the response of soft sphere assemblies to athermal cyclic shear deformation over a wide range of densities and amplitudes of shear deformation. Cyclic shear deformation induces transitions from reversible to irreversible behaviour in both unjammed and jammed soft sphere packings. Well above isotropic jamming density ($\bf{ϕ_J}$), this transition corresponds to yielding. In the vicinity of the jamming point, up to a higher density limit we designate ${\bf ϕ_J^{cyc}}$, an unjammed phase emerges between a localised, \emph{absorbing} phase, and a diffusive, {\emph irreversible} phase. The emergence of the unjammed phase signals the shifting of the jamming point to higher densities as a result of annealing, and opens a window where shear jamming becomes possible for frictionless packings. Below $\bf{ϕ_J}$, two distinct localised states, termed point and loop reversibile, are observed. We characterise in detail the different regimes and transitions between them, and obtain a unified density-shear amplitude phase diagram.

preprint2016arXiv

Effect of total and pair configurational entropy in determining dynamics of supercooled liquids over a range of densities

In this paper, we present a study of supercooled liquids interacting with the Lennard Jones (LJ) potential and the corresponding purely repulsive (Weeks-Chandler-Andersen or WCA) potential, over a range of densities and temperatures, in order to understand the origin of their different dynamics in spite of their structures being similar. Using the configurational entropy as the thermodynamic marker via the Adam Gibbs (AG) relation, we show that the difference in the dynamics of these two systems at low temperatures can be explained from thermodynamics. At higher densities both the thermodynamical and dynamical difference between these model systems decrease, which is quantitatively demonstrated in this paper by calculating different parameters. The study also reveals the origin of the difference in pair entropy despite the similarity in the structure. Although the maximum difference in structure is obtained in the partial radial distribution function (rdf) of the B type of particles, the rdf of AA pairs and AB pairs give rise to the differences in the entropy and dynamics. This work supports the observation made in an earlier study (Phys. Rev. Lett.,\textbf {113}, 225701, 2014) and shows that they are generic in nature, independent of density.

preprint2015arXiv

Disentangling the role of structure and friction in shear jamming

Amorphous packings of spheres have been intensely investigated in order to understand the mechanical and flow behaviour of dense granular matter, and to explore universal aspects of the transition from fluid to structurally arrested or jammed states. Considerable attention has recently been focussed on anisotropic packings of frictional grains generated by shear deformation leading to shear jamming, which occurs below the jamming density for isotropic packings of frictionless grains. With the aim of disentangling the role of shear deformation induced structures and friction in generating shear jamming, we study sheared assemblies of frictionless spheres computationally, over a wide range of densities, extending far below the isotropic jamming point. We demonstrate the emergence of a variety of geometric features characteristic of jammed packings with the increase of shear strain. The average contact number and the distributions of contact forces suggest the presence of a threshold density, well below the isotropic jamming point, above which a qualitative change occurs in the jamming behaviour of sheared configurations. We show that above this threshold density, friction stabilizes the sheared configurations we generate. Our results thus reveal the emergence of geometric features characteristic of jammed states as a result of shear deformation alone, while friction is instrumental in stabilising packings over a range of densities below the isotropic jamming point.

preprint2015arXiv

Free volume under shear

Using an athermal quasistatic simulation protocol, we study the distribution of free volumes in sheared hard-particle packings close to, but below, the random-close packing threshold. We show that under shear, and independent of volume fraction, the free volumes develop features similar to close-packed systems -- particles self-organize in a manner as to mimick the isotropically jammed state. We compare athermally sheared packings with thermalized packings and show that thermalization leads to an erasure of these structural features. The temporal evolution, in particular the opening-up and the closing of free-volume patches is associated with the single-particle dynamics, showing a crossover from ballistic to diffusive behavior.

preprint2015arXiv

Memory Effects in Schematic Models of Glasses Subjected to Oscillatory Deformation

We consider two schematic models of glasses subjected to oscillatory shear deformation, motivated by the observations, in computer simulations of a model glass, of a nonequilibrium transition from a localized to a diffusive regime as the shear amplitude is increased, and of persistent memory effects in the localized regime. The first of these schematic models is the NK model, a spin model with disordered multi-spin interactions previously studied as a model for sheared amorphous solids. The second model, a transition matrix model, is an abstract formulation of the manner in which occupancy of local energy minima evolves under oscillatory deformation cycles. In both of these models, we find a behavior similar to that of an atomic model glass studied earlier. We discuss possible further extensions of the approaches outlined.

preprint2015arXiv

Short-time $β$-relaxation in glass-forming liquids is cooperative in nature

Temporal relaxation of density fluctuations in supercooled liquids near the glass transition occurs in multiple steps. The short-time $β$-relaxation is generally attributed to spatially local processes involving the rattling motion of a particle in the transient cage formed by its neighbors. Using molecular dynamics simulations for three model glass-forming liquids, we show that the $β$-relaxation is actually cooperative in nature. Using finite-size scaling analysis, we extract a growing length-scale associated with $β$-relaxation from the observed dependence of the $β$-relaxation time on the system size. Remarkably, the temperature dependence of this length scale is found to be the same as that of the length scale that describes the spatial heterogeneity of local dynamics in the long-time $α$-relaxation regime. These results show that the conventional interpretation of $β$-relaxation as a local process is too simplified and provide a clear connection between short-time dynamics and long-time structural relaxation in glass-forming liquids.

preprint2015arXiv

Stochastic Approach to Plasticity and Yield in Amorphous Solids

We focus on the probability distribution function (pdf) $P(Δγ; γ)$ where $Δγ$ are the {\em measured} strain intervals between plastic events in an athermal strained amorphous solids, and $γ$ measures the accumulated strain. The tail of this distribution as $Δγ\to 0$ (in the thermodynamic limit) scales like $Δγ^η$. The exponent $η$ is related via scaling relations to the tail of the pdf of the eigenvalues of the {\em plastic modes} of the Hessian matrix $P(λ)$ which scales like $λ^θ$, $η=(θ-1)/2$. The numerical values of $η$ or $θ$ can be determined easily in the unstrained material and in the yielded state of plastic flow. Special care is called for in the determination of these exponents between these states as $γ$ increases. Determining the $γ$ dependence of the pdf $P(Δγ; γ)$ can shed important light on plasticity and yield. We conclude that the pdf's of both $Δγ$ and $λ$ are not continuous functions of $γ$. In slowly quenched amorphous solids they undergo two discontinuous transitions, first at $γ=0^+$ and then at the yield point $γ=γ_{_{\rm Y}}$ to plastic flow. In quickly quenched amorphous solids the second transition is smeared out due to the non existing stress peak before yield. The nature of these transitions and scaling relations with the system size dependence of $\langle Δγ\rangle$ are discussed.

preprint2015arXiv

Unraveling the success and failure of mode coupling theory from consideration of entropy

We analyze the dynamics of model supercooled liquids in a temperature regime where predictions of mode coupling theory (MCT) are known to be valid qualitatively. In this regime, the Adam-Gibbs (AG) relation, based on an activation picture of dynamics also describes the dynamics satisfactorily, and we explore the mutual consistency and interrelation of these descriptions. Although entropy and dynamics are related via phenomenological theories, the connection between MCT and entropy has not been argued for. In this work we explore this connection and provide a microscopic derivation of the phenomenological Rosenfeld theory. At low temperatures the overlap between MCT power law regime and AG relation implies that the AG relation predicts an avoided divergence at $T_c$, the origin of which is traced back to the vanishing of pair configurational entropy, which we find occurs at the same temperature. We also show that the residual multiparticle entropy plays an important role in describing the relaxation time

preprint2014arXiv

Free volume distribution of nearly jammed hard sphere packings

We calculate the free volume distributions of nearly jammed packings of monodisperse and bidisperse hard sphere configurations. These distributions differ qualitatively from those of the fluid, displaying a power law tail at large free volumes, which constitutes a distinct signature of nearly jammed configurations, persisting for moderate degrees of decompression. We reproduce and explain the observed distribution by considering the pair correlation function within the first coordination shell for jammed hard sphere configurations. We analyze features of the equation of state near jamming, and discuss the significance of observed asphericities of the free volumes to the equation of state.

preprint2014arXiv

Solvability-Based Comparison of Failure Detectors

Failure detectors are oracles that have been introduced to provide processes in asynchronous systems with information about faults. This information can then be used to solve problems otherwise unsolvable in asynchronous systems. A natural question is on the "minimum amount of information" a failure detector has to provide for a given problem. This question is classically addressed using a relation that states that a failure detector D is stronger (that is, provides "more, or better, information") than a failure detector D' if D can be used to implement D'. It has recently been shown that this classic implementability relation has some drawbacks. To overcome this, different relations have been defined, one of which states that a failure detector D is stronger than D' if D can solve all the time-free problems solvable by D'. In this paper we compare the implementability-based hierarchy of failure detectors to the hierarchy based on solvability. This is done by introducing a new proof technique for establishing the solvability relation. We apply this technique to known failure detectors from the literature and demonstrate significant differences between the hierarchies.

preprint2014arXiv

The role of structure and entropy in determining differences in dynamics for glass formers with different interaction potentials

We present a study of two model liquids with different interaction potentials, exhibiting similar structure but significantly different dynamics at low temperatures. By evaluating the configurational entropy, we show that the differences in the dynamics of these systems can be understood in terms of their thermodynamic differences. Analyzing their structure, we demonstrate that differences in pair correlation functions between the two systems, through their contribution to the entropy, dominate the differences in their dynamics, and indeed overestimate the differences. Including the contribution of higher order structural correlations to the entropy leads to smaller estimates for the relaxation times, as well as smaller differences between the two studied systems.

preprint2013arXiv

Density-temperature scaling of the fragility in a model glass-former

Dynamical quantities such as the diffusion coefficient and relaxation times for some glass-formers may depend on density and temperature through a specific combination, rather than independently, allowing the representation of data over ranges of density and temperature as a function of a single scaling variable. Such a scaling, referred to as density - temperature (DT) scaling, is exact for liquids with inverse power law (IPL) interactions but has also been found to be approximately valid in many non-IPL liquids. We have analyzed the consequences of DT scaling on the density dependence of the fragility in a model glass-former. We find the density dependence of kinetic fragility to be weak, and show that it can be understood in terms of DT scaling and deviations of DT scaling at low densities. We also show that the Adam-Gibbs relation exhibits DT scaling and the scaling exponent computed from the density dependence of the activation free energy in the Adam-Gibbs relation, is consistent with the exponent values obtained by other means.

preprint2013arXiv

Oscillatory athermal quasi-static deformation of a model glass

We report computer simulations of oscillatory athermal quasi-static shear deformation of dense amorphous samples of a three dimensional model glass former. A dynamical transition is observed as the amplitude of the deformation is varied: for large values of the amplitude the system exhibits diffusive behavior and loss of memory of the initial conditions, whereas localization is observed for small amplitudes. Our results suggest that the same kind of transition found in driven colloidal systems is present in the case of amorphous solids (e.g. metallic glasses). The onset of the transition is shown to be related to the onset of energy dissipation. Shear banding is observed for large system sizes, without, however, affecting qualitative aspects of the transition.

preprint2013arXiv

The Relationship of Dynamical Heterogeneity to the Adam-Gibbs and Random First-Order Transition Theories of Glass Formation

We carefully examine common measures of dynamical heterogeneity for a model polymer melt and test how these scales compare with those hypothesized by the Adam and Gibbs (AG) and random first-order transition (RFOT) theories of relaxation in glass-forming liquids. To this end, we first analyze clusters of highly mobile particles, the string-like collective motion of these mobile particles, and clusters of relative low mobility. We show that the time scale of the high-mobility clusters and strings is associated with a diffusive time scale, while the low-mobility particles' time scale relates to a structural relaxation time. We find a high degree of similarity in the geometrical structure of these particle clusters. In particular, we show that the fractal dimensions of these clusters are consistent with those of swollen branched polymers or branched polymers with screened excluded-volume interactions. In contrast, the fractal dimension of the strings crosses over from that of self-avoiding walks for small strings, to simple random walks for longer, {more strongly interacting, strings, corresponding to flexible polymers with screened excluded-volume interactions. We examine the appropriateness of identifying the size scales of either mobile particle clusters or strings with the size of cooperatively rearranging regions (CRR) in the AG and RFOT theories. We find that the string size appears to be the most consistent measure of CRR for both the AG and RFOT models. Identifying strings or clusters with the "mosaic" length of the RFOT model relaxes the conventional assumption that the "entropic droplets" are compact. We also confirm the validity of the entropy formulation of the AG theory, constraining the exponent values of the RFOT theory. This constraint, together with the analysis of size scales, enables us to estimate the characteristic exponents of RFOT.

preprint2012arXiv

Breakdown of the Stokes-Einstein relation in two, three and four dimensions

The breakdown of the Stokes-Einstein (SE) relation between diffusivity and viscosity at low temperatures is considered to be one of the hallmarks of glassy dynamics in liquids. Theoretical analyses relate this breakdown with the presence of heterogeneous dynamics, and by extension, with the fragility of glass formers. We perform an investigation of the breakdown of the SE relation in 2, 3 and 4 dimensions, in order to understand these interrelations. Results from simulations of model glass formers show that the degree of the breakdown of the SE relation decreases with increasing spatial dimensionality. The breakdown itself can be rationalized via the difference between the activation free energies for diffusivity and viscosity (or relaxation times) in the Adam-Gibbs relation in three and four dimensions. The behavior in two dimensions also can be understood in terms of a generalized Adam-Gibbs relation that is observed in previous work. We calculate various measures of heterogeneity of dynamics and find that the degree of the SE breakdown and measures of heterogeneity of dynamics are generally well correlated but with some exceptions. The two dimensional systems we study show deviations from the pattern of behavior of the three and four dimensional systems both at high and low temperatures. The fragility of the studied liquids is found to increase with spatial dimensionality, contrary to the expectation based on the association of fragility with heterogeneous dynamics.

preprint2012arXiv

Dependence of the fragility of a glass former on the softness of interparticle interactions

We study the influence of the softness of the interparticle interactions on the fragility of a glass former, by considering three model binary mixture glass formers. The interaction potential between particles is a modified Lennard-Jones type potential, with the repulsive part of the potential varying with an inverse power $q$ of the interparticle distance, and the attractive part varying with an inverse power $p$. We consider the combinations (12,11) (model I), (12,6) (model II) and (8,5) (model III) for (q,p) such that the interaction potential becomes softer from model I to III. We evaluate the kinetic fragilities from the temperature variation of diffusion coefficients and relaxation times, and a thermodynamic fragility from the temperature variation of the configuration entropy. We find that the kinetic fragility increases with increasing softness of the potential, consistent with previous results for these model systems, but at variance with the thermodynamic fragility, which decreases with increasing softness of the interactions, as well as expectations from earlier results. We rationalize our results by considering the full form of the Adam-Gibbs relation, which requires, in addition to the temperature dependence of the configuration entropy, knowledge of the high temperature activation energies ino rder to determine fragility. We show that consideration of the scaling of the high temperature activation energy with the liquid density, analyzed in recent studies, provides a partial rationalization of the observed behavior.

preprint2012arXiv

The Adam-Gibbs relation for glass-forming liquids in 2, 3 and 4 dimensions

The Adam-Gibbs relation between relaxation times and the configurational entropy has been tested extensively for glass formers using experimental data and computer simulation results. Although the form of the relation contains no dependence on the spatial dimensionality in the original formulation, subsequent derivations of the Adam-Gibbs relation allow for such a possibility. We test the Adam-Gibbs relation in 2, 3, and 4 spatial dimensions using computer simulations of model glass formers. We find that the relation is valid in 3 and 4 dimensions. But in 2 dimensions, the relation does not hold, and interestingly, no single alternate relation describes the results for the different model systems we study.

preprint2011arXiv

Computer Simulation Study of the Phase Behavior and Structural Relaxation in a Gel-Former Modeled by Three Body Interactions

We report a computer simulation study of a model gel-former obtained by modifying the three-body interactions of the Stillinger-Weber potential for silicon. This modification reduces the average coordination number and consequently shifts the liquid-gas phase coexistence curve to low densities, thus facilitating the formation of gels without phase separation. At low temperatures and densities, the structure of the system is characterized by the presence of long linear chains interconnected by a small number of three coordinated junctions at random locations. At small wave-vectors the static structure factor shows a non-monotonic dependence on temperature, a behavior which is due to the competition between the percolation transition of the particles and the stiffening of the formed chains. We compare in detail the relaxation dynamics of the system as obtained from molecular dynamics with the one obtained from Monte Carlo dynamics. We find that the bond correlation function displays stretched exponential behavior at moderately low temperatures and densities, but exponential relaxation at low temperatures. The bond lifetime shows an Arrhenius behavior, independent of the microscopic dynamics. For the molecular dynamics at low temperatures, the mean squared displacement and the (coherent and incoherent) intermediate scattering function display at intermediate times a dynamics with ballistic character and we show that this leads to compressed exponential relaxation. For the Monte Carlo dynamics we find always an exponential or stretched exponential relaxation. Thus we conclude that the compressed exponential relaxation observed in experiments is due to the out-of-equilibrium dynamics.

preprint2011arXiv

Liquid-liquid critical point in supercooled silicon

A novel liquid-liquid phase transition has been proposed and investigated in a wide variety of pure substances recently, including water, silica and silicon. From computer simulations using the Stillinger-Weber classical empirical potential, Sastry and Angell [1] demonstrated a first order liquid-liquid transition in supercooled silicon, subsequently supported by experimental and simulation studies. Here, we report evidence for a liquid-liquid critical end point at negative pressures, from computer simulations using the SW potential. Compressibilities exhibit a growing maximum upon lowering temperature below 1500 K and isotherms exhibit density discontinuities below 1120 K, at negative pressure. Below 1120 K, isotherms obtained from constant volume-temperature simulations exhibit non-monotonic, van der Waals-like behavior signaling a first order transition. We identify Tc ~ 1120 +/- 12 K, Pc -0.60 +/- 0.15 GPa as the critical temperature and pressure for the liquid-liquid critical point. The structure of the liquid changes dramatically upon decreasing the temperature and pressure. Diffusivities vary over 4 orders of magnitude, and exhibit anomalous pressure dependence near the critical point. A strong relationship between local geometry quantified by the coordination number, and diffusivity, is seen, suggesting that atomic mobility in both low and high density liquids can usefully be analyzed in terms of defects in the tetrahedral network structure. We have constructed the phase diagram of supercooled silicon. We identify the lines of compressibility, density extrema (maxima and minima) and the spinodal which reveal the interconnection between thermodynamic anomalies and the phase behaviour of the system as suggested in previous works [2-9]

preprint2009arXiv

Jamming transitions in amorphous packings of frictionless spheres occur over a continuous range of volume fractions

We numerically produce fully amorphous assemblies of frictionless spheres in three dimensions and study the jamming transition these packings undergo at large volume fractions. We specify four protocols yielding a critical value for the jamming volume fraction which is sharply defined in the limit of large system size, but is different for each protocol. Thus, we directly establish the existence of a continuous range of volume fraction where nonequilibrium jamming transitions occur. However, these jamming transitions share the same critical behaviour. Our results suggest that, even in the absence of partial crystalline ordering, a unique location of a random close packing does not exist, and that volume fraction alone is not sufficient to describe the properties of jammed states.

preprint2009arXiv

Structural Relaxation of a Gel Modeled by Three Body Interactions

We report a molecular dynamics simulation study of a model gel whose interaction potential is obtained by modifying the three body Stillinger-Weber model potential for silicon. The modification reduces the average coordination number, and suppresses the liquid-gas phase coexistence curve. The low density, low temperature equilibrium gel that can thus form exhibits interesting dynamical behavior, including compressed exponential relaxation of density correlations. We show that motion responsible for such relaxation has ballistic character, and arises from the motion of chain segments in the gel without the restructuring of the gel network.

preprint2008arXiv

On the relaxation dynamics of glass-forming systems: Insights from computer simulations

We discuss the relaxation dynamics of a simple lattice gas model for glass-forming systems and show that with increasing density of particles this dynamics slows down very quickly. By monitoring the trajectory of tagged particles we find that their motion is very heterogeneous in space and time, leading to regions in space in which there is a fast dynamics and others in which it is slow. We determine how the geometric properties of these quickly relaxing regions depend on density and time. Motivated by this heterogeneous hopping dynamics, we use a simple model, a variant of a continuous time random walk, to characterize the relaxation dynamics. In particular we find from this model that for large displacements the self part of the van Hove function shows an exponential tail, in agreement with recent findings from experiments and simulations of glass-forming systems.

preprint2001arXiv

Spectral Statistics of Instantaneous Normal Modes in Liquids and Random Matrices

We study the statistical properties of eigenvalues of the Hessian matrix ${\cal H}$ (matrix of second derivatives of the potential energy) for a classical atomic liquid, and compare these properties with predictions for random matrix models (RMM). The eigenvalue spectra (the Instantaneous Normal Mode or INM spectra) are evaluated numerically for configurations generated by molecular dynamics simulations. We find that distribution of spacings between nearest neighbor eigenvalues, s, obeys quite well the Wigner prediction $s exp(-s^2)$, with the agreement being better for higher densities at fixed temperature. The deviations display a correlation with the number of localized eigenstates (normal modes) in the liquid; there are fewer localized states at higher densities which we quantify by calculating the participation ratios of the normal modes. We confirm this observation by calculating the spacing distribution for parts of the INM spectra with high participation ratios, obtaining greater conformity with the Wigner form. We also calculate the spectral rigidity and find a substantial dependence on the density of the liquid.

preprint1999arXiv

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

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