Source author record

Smarajit Karmakar

Smarajit Karmakar 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

37works
6topics
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

37 published item(s)

preprint2026arXiv

Anomalous Dynamical Heterogeneity in Active Glasses as a Signature of Violation of Mermin-Wagner-Hohenberg Theorem

Two-dimensional (2D) systems have attracted renewed interest within the scientific community due to their anomalous dynamical behaviors, which arise from long-wavelength density fluctuations as predicted by the Mermin-Wagner-Hohenberg (MWH) theorem. In equilibrium, it is well established that continuous spontaneous symmetry breaking (SSB) in 2D is prohibited at any finite temperature ($T > 0$), resulting in the absence of true long-range positional order and establishing $d_l = 2$ as the lower critical dimension. Recent studies have demonstrated that, in active systems, the lower critical dimension can shift from $d_l = 2$ to $3$. This study examines the impact of MWH theorem violation in active systems on dynamical heterogeneity (DH). As a minimal model, glassy systems of active particles undergoing run-and-tumble (RT) motion are considered. Glass-like dynamical behavior, including anomalously enhanced DH, is observed in various biological systems such as collective cell migration, bacterial cytoplasm, and ant colonies. Furthermore, the study investigates the influence of local positional order, or medium-range crystalline order (MRCO), on DH in the presence of activity. The results indicate that the growth of DH with increasing activity differs significantly between systems with and without MRCO. These findings may have important implications, as many biological systems exhibit local structural ordering, and DH could serve as a useful indicator for quantifying the degree of ordering.

preprint2026arXiv

Enhanced Long Wavelength Mermin-Wagner Fluctuations in Active Crystals and Glasses

In two-dimensions (2D), the Mermin-Wagner-Hohenberg (MWH) fluctuation plays a significant role, giving rise to striking dimensionality effects marked by long-range density fluctuations leading to the singularities of various dynamical properties. According to the MWH theorem, a 2D equilibrium system with continuous degrees of freedom cannot achieve long-range crystalline order at non-zero temperatures. Recently, MWH fluctuations have been observed in glass-forming liquids, evidenced by the logarithmic divergence in the plateau value of mean squared displacement (MSD). Our research investigates long-wavelength fluctuations in crystalline and glassy systems influenced by non-equilibrium active noises. Active systems serve as a minimal model for understanding diverse non-equilibrium dynamics, such as those in biological systems and self-propelled colloids. We demonstrate that fluctuations from active forces can strongly couple with long-wavelength density fluctuations, altering the lower critical dimension ($d_l$) from $2$ to $3$ and leading to a novel logarithmic divergence of the MSD plateau with system size in 3D.

preprint2023arXiv

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(ω) \sim ω^4$ law changes to $D(ω) \sim ω^δ$ with $δ\approx 5$ in two dimensions and $δ\approx 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(ω) \sim ω^4$ behavior is recovered.

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.

preprint2021arXiv

Enhanced Phonon Peak in Four-point Dynamic Susceptibility in the Supercooled Active Glass-forming Liquids

Active glassy systems can be thought of as simple model systems that imitate complex biological systems. Sometimes, it becomes crucial to estimate the amount of the activity present in such biological systems, such as predicting the progression rate of the cancer cells or the healing time of the wound. In this work, we study a model active glassy system to understand a possible quantification of the degree of activity from the collective, long-range phonon response in the system. We find that the four-point dynamic susceptibility, $χ_4(t)$ at the phonon timescale, grows with increased activity. We then show how one can estimate the degree of the activity at such a small timescale by measuring the growth of $χ_4(t)$ with changing activity. A detailed finite size analysis of this measurement, shows that the peak height of $χ_4(t)$ at this phonon timescale increases strongly with increasing system size suggesting a possible existence of an intrinsic dynamic length scale that grows with increasing activity. Finally, we show that this peak height is a unique function of effective activity across all system sizes, serving as a possible parameter for characterizing the degree of activity in a system.

preprint2021arXiv

Thermodynamics and its correlation with dynamics in a mean-field model and pinned systems: A comparative study using two different methods of entropy calculation

Recently, some of us developed a novel model glass-forming liquid with k extra interactions with pseudo neighbours to each liquid particle over and above the regular interactions with its neighbours. Analysis of the structure and dynamics of these systems showed that with an increase in k the systems have more mean-field like properties. This work presents an extensive study of the thermodynamics of the above-mentioned model for several values of k and its correlation with the dynamics. We surprisingly find that the usual thermodynamic integration (TI) method of calculating the entropy provides unphysical results for this model. It predicts the vanishing of configurational entropy at state points at which both the collective and the single-particle dynamics of the system show complete relaxation. We then employ a new method known as the two-phase thermodynamics (2PT) method to calculate the entropy. We find that with an increase in k the difference in the entropy computed using the two methods (2PT and TI) increases. We also find that in the temperature range studied, the entropy calculated via the 2PT method satisfies the Adam-Gibbs (AG) relationship between the relaxation time and the configurational entropy, whereas the entropy calculated via the TI method shows a strong violation of the same. We then apply the 2PT method to calculate the entropy in another system where some fractions of particles are pinned randomly in their equilibrium positions. This system also shows a similar breakdown of the AG relationship as reported earlier. We show that the difference in entropy calculated via the 2PT and TI methods increases with an increase in pinning density. We also find that when the entropy is calculated using the 2PT method, the AG relationship between the dynamics and the entropy holds.

preprint2021arXiv

Universal non-Debye low-frequency vibrations in sheared amorphous solids

We study energy minimized configurations of amorphous solids with a simple shear degree of freedom. We show that the low-frequency regime of the vibrational density of states of structural glass formers is crucially sensitive to the stress-ensemble from which the configurations are sampled. In both two and three dimensions, a shear-stabilized ensemble displays a $D(ω_{\min}) \sim ω^{5}_{\min}$ regime, as opposed to the $ω^{4}_{\min}$ regime observed under unstrained conditions. We also study an ensemble of two dimensional, strained amorphous solids near a plastic event. We show that the minimum eigenvalue distribution at a strain $γ$ near the plastic event occurring at $γ_{P}$, displays a collapse when scaled by $\sqrt{γ_P - γ}$, and with the number of particles as $N^{-0.22}$. Notably, at low-frequencies, this scaled distribution displays a robust $D(ω_{\min}) \sim ω^{6}_{\min}$ power-law regime, which survives in the large $N$ limit. Finally, we probe the universal properties of this ensemble through a characterization of the second and third eigenvalues of the Hessian matrix near a plastic event.

preprint2019arXiv

Signature of Dynamical Heterogeneity in Spatial Correlations of Particle Displacement and its Temporal Evolution in Supercooled Liquids

The existence of heterogeneity in the dynamics of supercooled liquids is believed to be one of the hallmarks of the glass transition. Intense research has been carried out in the past to understand the origin of this heterogeneity in dynamics and a possible length scale associated with it. We have done extensive molecular dynamics simulations of few model glass-forming liquids in three dimensions to understand the temporal evolution of the dynamic heterogeneity and the heterogeneity length scale. We find that although the strength of the dynamic heterogeneity is maximum at a timescale close to characteristic $α$-relaxation time of the system, dynamic heterogeneity itself is well-developed at timescale as short as $β$-relaxation time and survives up to a timescale as long as few tens of $α$-relaxation time. Moreover, we discovered that temperature dependence of heterogeneity length remains the same in the whole time window although its absolute value changes over time in a non-monotonic manner.

preprint2016arXiv

Amorphization of a crystalline solid by plastic deformation

We demonstrate that plastic deformation in solids is associated with a dynamic transition that is reminiscent to the transition from a superconducting to a mixed phase in Type II superconductors. We report analytic calculations, extensive molecular dynamics and sequential umbrella sampling Monte Carlo simulations of a two dimensional triangular crystalline solid undergoing plastic deformation under strain. The solid consists of particles connected by harmonic springs. Non-affine displacement fluctuations of the solid are enhanced using an external field, causing a rich deformation behaviour. The external field leads to a mixed phase, where defect and stress-free crystallites are surrounded by a network of amorphous regions with large local internal stress --- a "stress Meissner" effect. The transition shows slow ageing dynamics caused by the presence of many competing, non-crystalline free-energy minima. Under uniform uniaxial strain, this amorphization transition gives rise to irreversible plastic deformation with the amorphous regions appearing as band-like structures. Our results may be checked by careful experiments on colloidal crystals using holographic optical tweezers.

preprint2016arXiv

Understanding Stokes-Einstein Relation in Supercooled Liquids using Random Pinning

Breakdown of Stokes-Einstein relation in supercooled liquids is believed to be one of the hallmarks of glass transition. The phenomena is studied in depth over many years to understand the microscopic mechanism without much success. Recently it was found that violation of Stokes-Einstein relation in supercooled liquids can be tuned very systematically by pinning randomly a set of particles in their equilibrium positions. This observation suggested a possible framework where breakdown of Stokes-Einstein relation in the dynamics of supercooled liquids can be studied with precise control. We have done extensive molecular dynamics simulations to understand this phenomena by analyzing the structure of appropriately defined set of dynamically slow and fast particles clusters. We have shown that the Stokes-Einstein breakdown actually become predominant once the cluster formed by the slow particles percolate the entire system size. Finally we proposed a possible close connection between fractal dimensions of these clusters and the exponents associated with the fractional Stokes-Einstein relation.

preprint2016arXiv

Understanding the Dynamics of Glass-forming Liquids with Random Pinning within the Random First Order Transition Theory

Extensive computer simulations are performed for a few model glass-forming liquids in both two and three dimensions to study their dynamics when a randomly chosen fraction of particles are frozen in their equilibrium positions. For all the studied systems, we find that the temperature-dependence of the $α$ relaxation time extracted from an overlap function related to the self part of the density autocorrelation function can be explained within the framework of the Random First Order Transition (RFOT) theory of the glass transition. We propose a scaling description to rationalize the simulation results and show that our data for the $α$ relaxation time for all temperatures and pin concentrations are consistent with this description. We find that the fragility parameter obtained from fits of the temperature dependence of the $α$ relaxation time to the Vogel-Fulcher-Tammann (VFT) form decreases by almost an order of magnitude as the pin concentration is increased from zero. Our scaling description relates the fragility parameter to the static length scale of RFOT and thus provides a physical understanding of fragility within the framework of the RFOT theory. Implications of these findings for the values of the exponents appearing in the RFOT theory are discussed.

preprint2015arXiv

Dynamics of Glass Forming Liquids with Randomly Pinned Particles

It is frequently assumed that in the limit of vanishing cooling rate, the glass transition phenomenon becomes a thermodynamic transition at a temperature $T_{K}$. However, with any finite cooling rate, the system falls out of equilibrium at temperatures near $T_g(>T_{K})$, implying that the very existence of the putative thermodynamic phase transition at $T_{K}$ can be questioned. Recent studies of systems with randomly pinned particles have hinted that the thermodynamic glass transition may be observed in simulations and experiments carried out for liquids with randomly pinned particles. This expectation is based on the results of approximate calculations that suggest that the temperature of the thermodynamic glass transition increases as the concentration of pinned particles is increased and it may be possible to equilibrate the system at temperatures near the increased transition temperature. We test the validity of this prediction through extensive molecular dynamics simulations of two model glass-forming liquids in the presence of random pinning. We fit the temperature-dependence of the structural relaxation time to the Vogel-Fulcher-Tammann form that predicts a divergence of the relaxation time at a temperature $T_{VFT}$ and identify this temperature with the thermodynamic transition temperature $T_K$. We find that $T_{VFT}$ does not show any sign of increasing with increasing concentration of pinned particles. The main effect of pinning is found to be a rapid decrease in the kinetic fragility of the system with increasing pin concentration. Implications of these observations for current theories of the glass transition are discussed.

preprint2015arXiv

Excess Vibrational Density of States and the Brittle to Ductile Transition in Crystalline and Amorphous Solids

The conditions which determine whether a material behaves in a brittle or ductile fashion on mechanical loading are still elusive and comprise a topic of active research among materials physicists and engineers. In this study, we present results of {\em in silico} mechanical deformation experiments from two very different model solids in two and three dimensions. The first consists of particles interacting with isotropic potentials and the other has strongly direction dependent interactions. We show that in both cases, the excess vibrational density of states is the fundamental quantity which characterises the ductility of the material. Our results can be checked using careful experiments on colloidal solids.

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

Vanishing of configurational entropy may not imply an ideal glass transition in randomly pinned liquids

Ozawa et. al [1] presented numerical results for the configurational entropy density, $s_c$, of a model glass-forming liquid in the presence of random pinning. The location of a "phase boundary" in the pin density ($c$) - temperature ($T$) plane, that separates an "ideal glass" phase from the supercooled liquid phase, is obtained by finding the points at which $s_c(T,c) \to 0$. According to the theoretical arguments by Cammarota et. al. [2], an ideal glass transition at which the $α$-relaxation time $τ_α$ diverges takes place when $s_c$ goes to zero. We have studied the dynamics of the same system using molecular dynamics simulations. We have calculated the time-dependence of the self intermediate scattering function, $F_s(k,t)$ at three state points in the $(c-T)$ plane where $s_c(T,c) \simeq 0$ according to Ref. [1]. It is clear from the plots that the relaxation time is finite [$τ_α\sim \mathcal{O}(10^6)]$ at these state points. Similar conclusions have been obtained in Ref.[3] where an overlap function was used to calculate $τ_α$ at these state points.

preprint2014arXiv

Distribution of Diffusion Constants and Stokes-Einstein Violation in supercooled liquids

It is widely believed that the breakdown of the Stokes-Einstein (SE) relation between the translational diffusivity and the shear viscosity in supercooled liquids is due to the development of dynamic heterogeneity i.e. the presence of both slow and fast moving particles in the system. In this study we \emph{directly} calculate the distribution of the diffusivity for a model system for different temperatures in the supercooled regime. We find that with decreasing temperature, the distribution evolves from Gaussian to bimodal indicating that on the time scale of the typical relaxation time, mobile (fluid like) and less mobile (solid like) particles in the system can be \emph{unambiguously} identified. We also show that less mobile particles obey the Stokes-Einstein relation even in the supercooled regime and it is the mobile particles which show strong violation of the Stokes-Einstein relation in agreement with the previous studies on different model glass forming systems. Motivated by some of the recent studies where an ideal glass transition is proposed by randomly pinning some fraction of particles, we then studied the SE breakdown as a function of random pinning concentration in our model system. We showed that degree of SE breakdown increases quite dramatically with increasing pinning concentration, thereby providing a new way to unravel the puzzles of SE violation in supercooled liquids in greater details.

preprint2014arXiv

Glass-like slow dynamics in a colloidal solid with multiple ground states

We study the phase ordering dynamics of a two dimensional model colloidal solid using molecular dynamics simulations. The colloid particles interact with each other with a Hamaker potential modified by the presence of equatorial "patches" of attractive and negative regions. The total interaction potential between two such colloids is, therefore, strongly directional and has three-fold symmetry. Working in the canonical ensemble, we determine the tentative phase diagram in the density-temperature plane which features three distinct crystalline ground states viz, a low density honeycomb solid followed by a rectangular solid at higher density, which eventually transforms to a close packed triangular structure as the density is increased further. We show that when cooled rapidly from the liquid phase along isochores, the system undergoes a transition to a "strong glass" while slow cooling gives rise to crystalline phases. We claim that geometrical frustration arising from the presence of many crystalline ground states causes glassy ordering and dynamics in this solid. Our results may be easily confirmed by suitable experiments on patchy colloids.

preprint2013arXiv

Comparison of Static Length-Scales Characterizing the Glass Transition

The dramatic dynamic slowing down associated with the glass transition is considered by many to be related to the existence of a static length scale that grows when temperature decreases. Defining, identifying and measuring such a length is a subtle and non-trivial problem. Recently, two proposals, based on very different insights regarding the relevant physics, were put forward. One approach is based on the point-to-set correlation technique and the other on the scale where the lowest eigenvalue of the Hessian matrix becomes sensitive to disorder. In this Letter we present numerical evidence that the two approaches might result in the same identical length scale. This provides further mutual support to their relevance and, at the same time, raise interesting theoretical questions, discussed in the conclusion, concerning the fundamental reason for their relationship.

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

Finite Size Scaling for the Glass Transition: the Role of a Static Length Scale

Over the last decade computer simulations have had an increasing role in shedding light on difficult statistical physical phenomena and in particular on the ubiquitous problem of the glass transition. Here in a wide variety of materials the viscosity of a super-cooled liquid increases by many orders of magnitude upon decreasing the temperature over a modest range. A natural concern in these computer simulation is the very small size of the simulated systems compared to experimental ones, raising the issue of how to assess the thermodynamic limit. Here we offer a theory for the glass transition based on finite size scaling, a method that was found very useful in the context of critical phenomena and other interesting problems. As is known, the construction of such a theory rests crucially on the existence of a growing {\em static} length scale upon decreasing the temperature. We demonstrate that the static length scale that was discovered in Ref. \cite{12KLP} fits the bill extremely well, allowing us to provide a finite size scaling theory for the $α$ relaxation time of the glass transition, including predictions for the thermodynamic limit based on simulations in small systems.

preprint2012arXiv

Random Pinning Glass Model

Glass transition where viscosity of liquids increases dramatically upon decrease of temperature without any major change in structural properties, remains one of the most challenging problems in condensed matter physics (Cavagna, 2009; Berthier and Biroli, 2011) in spite of tremendous research efforts in last decades. On the other hand disordered freezing of spins in a magnetic materials with decreasing temperature, the so-called spin glass transition, is relatively better understood (Mezard, Parisi and Virasoro, 1987; Castellani and Cavagna, 2005). Previously found similarity between some spin glass models with the structural glasses (Kirkpatrick and Thirumalai, 1987; Kirkpatrick and Wolynes, 1987; Kirkpatrick and Wolynes, 1987; Franz and Parisi, 1999; Moore and Drossel, 2002) inspired development of theories of structural glasses (Kirkpatrick, Thirumalai and Wolynes, 1989; Barrat, Franz and Parisi, 1997; Mézard and Parisi, 1999; Lubchenko and Wolynes, 2007; Biroli and Bouchaud, 2012) based on the scenario of spin glass transition. This scenario though looks very appealing is still far from being well established. One of the main differences between standard spin systems to molecular systems is the absence of quenched disorder and the presence of translational invariance: it often assumed that this difference is not relevant, but this conjecture is still far from being established. The quantities, which are well defined and characterized for spin models, are not easily calculable for molecular glasses due to the lack of quenched disorder which breaks the translational invariance in the system and the characterization of the similarity between the spin and the structural glass transition remained an elusive subject still now. In this study we introduced a model structural glass with built in quenched disorder which alleviates this main difference between the spin and molecular glasses thereby helping us to compare these two systems: the possibility of producing a good thermalization at rather low temperatures is one of the advantages of this model.

preprint2012arXiv

Relaxation Mechanisms in Glassy Dynamics: the Arrhenius and Fragile Regimes

Generic glass formers exhibit at least two characteristic changes in their relaxation behavior, first to an Arrhenius-type relaxation at some characteristic temperature, and then at a lower characteristic temperature to a super-Arrhenius (fragile) behavior. We address these transitions by studying the statistics of free energy barriers for different systems at different temperatures and space dimensions. We present a clear evidence for changes in the dynamical behavior at the transition to Arrhenius and then to a super-Arrhenius behavior. A simple model is presented, based on the idea of competition between single-particle and cooperative dynamics. We argue that Arrhenius behavior can take place as long as there is enough free volume for the completion of a simple $T_1$ relaxation process. Once free volume is absent one needs a cooperative mechanism to `collect' enough free volume. We show that this model captures all the qualitative behavior observed in simulations throughout the considered temperature range.

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.

preprint2012arXiv

The nature of the β-peak in the loss modulus of amorphous solids

Glass formers exhibit, upon an oscillatory excitation, a response function whose imaginary and real parts are known as the loss and storage moduli respectively. The loss modulus typically peaks at a frequency known as the αfrequency which is associated with the main relaxation mechanism of the super-cooled liquid. In addition, the loss modulus is decorated by a smaller peak, shoulder or wing which is referred to as the β-peak. The physical origin of this secondary peak had been debated for decades, with proposed mechanisms ranging from highly localized relaxations to entirely cooperative ones. Using numerical simulations we bring an end to the debate, exposing a clear and unique cooperative mechanism for the said β-peak which is distinct from that of the $α$-peak.

preprint2011arXiv

An Athermal Brittle to Ductile Transition in Amorphous Solids

Brittle materials exhibit sharp dynamical fractures when meeting Griffith's criterion, whereas ductile materials blunt a sharp crack by plastic responses. Upon continuous pulling ductile materials exhibit a necking instability which is dominated by a plastic flow. Usually one discusses the brittle to ductile transition as a function of increasing temperature. We introduce an athermal brittle to ductile transition as a function of the cut-off length of the inter-particle potential. On the basis of extensive numerical simulations of the response to pulling the material boundaries at a constant speed we offer an explanation of the onset of ductility via the increase in the density of plastic modes as a function of the potential cutoff length. Finally we can resolve an old riddle: in experiments brittle materials can be strained under grip boundary conditions, and exhibit a dynamic crack when cut with a sufficiently long initial slot. Mysteriously, in molecular dynamics simulations it appeared that cracks refused to propagate dynamically under grip boundary conditions, and continuous pulling was necessary to achieve fracture. We argue that this mystery is removed when one understands the distinction between brittle and ductile athermal amorphous materials.

preprint2011arXiv

Direct Estimate of the Static Length-Scale Accompanying the Glass Transition

Characterizing the glass state remains elusive since its distinction from a liquid state is not obvious. Glasses are liquids whose viscosity has increased so much that they cannot flow. Accordingly there have been many attempts to define a static length-scale associated with the dramatic slowing down of supercooled liquid with decreasing temperature. Here we present a simple method to extract the desired length-scale which is highly accessible both for experiments and for numerical simulations. The fundamental idea is that low lying vibrational frequencies come in two types, those related to elastic response and those determined by plastic instabilities. The minimal observed frequency is determined by one or the other, crossing at a typical length-scale which is growing with the approach of the glass transition. This length-scale characterizes the correlated disorder in the system, where on longer length-scales the details of the disorder become irrelevant, dominated by the Debye model of elastic modes.

preprint2011arXiv

Do Athermal Amorphous Solids Exist?

We study the elastic theory of amorphous solids made of particles with finite range interactions in the thermodynamic limit. For the elastic theory to exist one requires all the elastic coefficients, linear and nonlinear, to attain a finite thermodynamic limit. We show that for such systems the existence of non-affine mechanical responses results in anomalous fluctuations of all the nonlinear coefficients of the elastic theory. While the shear modulus exists, the first nonlinear coefficient B_2 has anomalous fluctuations and the second nonlinear coefficient B_3 and all the higher order coefficients (which are non-zero by symmetry) diverge in the thermodynamic limit. These results put a question mark on the existence of elasticity (or solidity) of amorphous solids at finite strains, even at zero temperature. We discuss the physical meaning of these results and propose that in these systems elasticity can never be decoupled from plasticity: the nonlinear response must be very substantially plastic.

preprint2011arXiv

Exposing the static scale of the glass transition by random pinning

The dramatic slowing down associated with the glass transition cannot be fully understood without an associated static length that is expected to increase rapidly as the temperature is reduced. The search for such a length was long and arduous, without a universally accepted candidate at hand. Recently a natural such length $ξ_s$ was proposed, stemming from a cross-over between plastic and elastic mechanical responses of the material. In this Letter we show that supercooled liquids in which there exists random pinning sites of density $ρ_{\rm im}\sim 1/ξ_s^d$ exhibit complete jamming of all dynamics. This is a direct demonstration that the proposed length scale is indeed {\em the} static length that was long sought-after.

preprint2011arXiv

Universality of the Plastic Instability in Strained Amorphous Solids

By comparing the response to external strains in metallic glasses and in Lenard-Jones glasses we find a quantitative universality of the fundamental plastic instabilities in the athermal, quasistatic limit. Microscopically these two types of glasses are as different as one can imagine, the latter being determined by binary interactions, whereas the former by multiple interactions due to the effect of the electron gas that cannot be disregarded. In spite of this enormous difference the plastic instability is the same saddle-node bifurcation. As a result the statistics of stress and energy drops in the elasto-plastic steady state are universal, sharing the same system-size exponents.

preprint2010arXiv

Athermal Nonlinear Elastic Constants of Amorphous Solids

We derive expressions for the lowest nonlinear elastic constants of amorphous solids in athermal conditions (up to third order), in terms of the interaction potential between the constituent particles. The effect of these constants cannot be disregarded when amorphous solids undergo instabilities like plastic flow or fracture in the athermal limit; in such situations the elastic response increases enormously, bringing the system much beyond the linear regime. We demonstrate that the existing theory of thermal nonlinear elastic constants converges to our expressions in the limit of zero temperature. We motivate the calculation by discussing two examples in which these nonlinear elastic constants play a crucial role in the context of elasto-plasticity of amorphous solids. The first example is the plasticity-induced memory that is typical to amorphous solids (giving rise to the Bauschinger effect). The second example is how to predict the next plastic event from knowledge of the nonlinear elastic constants. Using the results of this paper we derive a simple differential equation for the lowest eigenvalue of the Hessian matrix in the external strain near mechanical instabilities; this equation predicts how the eigenvalue vanishes at the mechanical instability and the value of the strain where the mechanical instability takes place.

preprint2010arXiv

Plasticity-Induced Anisotropy in Amorphous Solids: the Bauschinger Effect

Amorphous solids that underwent a strain in one direction such that they responded in a plastic manner `remember' that direction also when relaxed back to a state with zero mean stress. We address the question `what is the order parameter that is responsible for this memory?' and is therefore the reason for the different subsequent responses of the material to strains in different directions. We identify such an order parameter which is readily measurable, we discuss its trajectory along the stress-strain curve, and propose that it and its probability distribution function must form a necessary component of a theory of elasto-plasticity.

preprint2010arXiv

Predicting plastic flow events in athermal shear-strained amorphous solids

We propose a method to predict the value of the external strain where a generic amorphous solid will fail by a plastic response (i.e. an irreversible deformation), solely on the basis of measurements of the nonlinear elastic moduli. While usually considered fundamentally different, with the elastic properties describing reversible phenomena and plastic failure epitomizing irreversible behavior, we show that the knowledge of some nonlinear elastic moduli is enough to predict where plasticity sets in.

preprint2010arXiv

Statistical Physics of Elasto-Plastic Steady States in Amorphous Solids: Finite Temperatures and Strain Rates

The effect of finite temperature $T$ and finite strain rate $\dotγ$ on the statistical physics of plastic deformations in amorphous solids made of $N$ particles is investigated. We recognize three regimes of temperature where the statistics are qualitatively different. In the first regime the temperature is very low, $T<T_{\rm cross}(N)$, and the strain is quasi-static. In this regime the elasto-plastic steady state exhibits highly correlated plastic events whose statistics are characterized by anomalous exponents. In the second regime $T_{\rm cross}(N)<T<T_{\rm max}(\dotγ)$ the system-size dependence of the stress fluctuations becomes normal, but the variance depends on the strain rate. The physical mechanism of the cross-over is different for increasing temperature and increasing strain rate, since the plastic events are still dominated by the mechanical instabilities (seen as an eigenvalue of the Hessian matrix going to zero), and the effect of temperature is only to facilitate the transition. A third regime occurs above the second cross-over temperature $T_{\rm max}(\dotγ)$ where stress fluctuations become dominated by thermal noise. Throughout the paper we demonstrate that scaling concepts are highly relevant for the problem at hand, and finally we present a scaling theory that is able to collapse the data for all the values of temperatures and strain rates, providing us with a high degree of predictability.

preprint2010arXiv

Statistical Physics of the Yielding Transition in Amorphous Solids

The art of making structural, polymeric and metallic glasses is rapidly developing with many applications. A limitation to their use is their mechanical stability: under increasing external strain all amorphous solids respond elastically to small strains but have a finite yield stress which cannot be exceeded without effecting a plastic response which typically leads to mechanical failure. Understanding this is crucial for assessing the risk of failure of glassy materials under mechanical loads. Here we show that the statistics of the energy barriers ΔE that need to be surmounted changes from a probability distribution function (pdf) that goes smoothly to zero to a pdf which is finite at ΔE=0. This fundamental change implies a dramatic transition in the mechanical stability properties with respect to external strain. We derive exact results for the scaling exponents that characterize the magnitudes of average energy and stress drops in plastic events as a function of system size.

preprint2010arXiv

What Determines the Yield Stress in Amorphous Solids?

A crucially important material parameter for all amorphous solids is the yield stress, which is the value of the stress for which the material yields to plastic flow when it is strained quasi-statically at zero temperature. It is difficult in laboratory experiments to determine what parameters of the inter-particle potential effect the value of the yield stress. Here we use the versatility of numerical simulations to study the dependence of the yield stress on the parameters of the inter-particle potential. We find a very simple dependence on the fundamental scales which characterize the repulsive and attractive parts of the potential respectively, and offer a scaling theory that collapses the data for widely different potentials and in different space dimensions.

preprint2009arXiv

Size of Plastic Events in Strained Amorphous Solids at Finite Temperatures

We address the system-size dependence of typical plastic flow events when an amorphous solid is put under a fixed external strain rate at a finite temperature. For system sizes that are accessible to numerical simulations at reasonable strain rates and at low temperatures the magnitude of plastic events grows with the system size. We explain however that this must be a finite size effect; for larger systems there exist two cross-over length-scales $ξ_1$ and $ξ_2$, the first determined by the elastic time-scale and the second by the thermal energy-scale. For system of linear size $L$ larger than any of these scales the magnitude of plastic events must be bounded. For systems of size $L\gg ξ$ there must exist $(L/ξ)^d$ uncorrelated plastic events which occur simultaneously. We present a scaling theory that culminates with the dependence of the cross-over scales on temperature and strain rate. Finally we relate these findings to the temperature and size dependence of the stress fluctuations. We comment on the importance of these considerations for theories of elasto-plasticity.

preprint2005arXiv

Equilibrium glassy phase in a polydisperse hard sphere system

The phase diagram of a polydisperse hard sphere system is examined by numerical minimization of a discretized form of the Ramakrishnan-Yussouff free energy functional. Crystalline and glassy local minima of the free energy are located and the phase diagram in the density-polydispersity plane is mapped out by comparing the free energies of different local minima. The crystalline phase disappears and the glass becomes the equilibrium phase beyond a "terminal" value of the polydispersity. A crystal to glass transition is also observed as the density is increased at high polydispersity. The phase diagram obtained in our study is qualitatively similar to that of hard spheres in a quenched random potential.