Source author record

Itamar Procaccia

Itamar Procaccia 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

113works
17topics
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

113 published item(s)

preprint2022arXiv

Anomalous Elasticity in Classical Glass-formers

Amorphous solids under mechanical strains are prone to plastic responses. Recent work showed that in amorphous granular system these plastic events, that are typically quadrupolar in nature, can screen the elastic response. When the density of the quadrupoles is high, the gradients of the quadrupole field act as emergent dipole sources, leading to qualitative changes in the mechanical response, as seen for example in the displacement field. In this paper we examine the effect of screening in classical glass formers. These are made of point particles that interact via binary forces. Both inverse power law forces and Lennard-Jones interactions are examined, and it is shown that in both cases the elastic response can be strongly screened, in agreement with the novel theory. The degree of deviation from classical elasticity theory is parameterized by a proposed new measure that is shown to have a functional dependence of on the amount of energy lost to plastic responses.

preprint2022arXiv

Creep failure of amorphous solids under tensile stress

Applying constant tensile stress to a piece of amorphous solid results in a slow extension, followed by an eventual rapid mechanical collapse. This "creep" process is of paramount engineering concern, and as such was the subject of study in a variety of materials, for more than a century. Predictive theories for $τ_w$, the expected time of collapse, are lacking, mainly due to its dependence on a bewildering variety of parameters, including temperature, system size, tensile force, but also the detailed microscopic interactions between constituents. The complex dependence of the collapse time on all the parameters is discussed below, using simulations of strip of amorphous material. Different scenarios are observed for ductile and brittle materials, resulting in serious difficulties in creating an all-encompassing theory that could offer safety measures for given conditions. A central aim of this paper is to employ scaling concepts, to achieve data collapse for the probability distribution function (pdf) of $\ln{τ_w}$. The scaling ideas result in a universal function which provides a prediction of the pdf of $\ln{τ_w}$ for out-of-sample systems, from measurements at other values of these parameters. The predictive power of the scaling theory is demonstrated for both ductile and brittle systems. Finally, we present a derivation of universal scaling function for brittle materials. The ductile case appears to be due to a plastic necking instability and is left for future research.

preprint2022arXiv

Experimental and Numerical Verification of Anomalous Screening Theory in Granular Matter

The concept of mechanical screening is widely applied in solid-state systems. Examples include nucleation of defects in crystalline materials, scars and pleats in curved crystals, wrinkles in strongly confined thin sheets, and cell-rearrangements in biological tissue. Available theories of such screening usually contain a crucial ingredient, which is the existence of an ordered reference state, with respect to which screening elements nucleate to release stresses. In contradistinction, amorphous materials in which a unique reference state does not exist, nevertheless exhibits plastic events that act as screening geometric charges with significant implications on the mechanical response. In a recent paper [Phys. Rev. E 104, 024904] it was proposed that mechanical strains in amorphous solids can be either weakly or strongly screened by the formation of low or high density of plastic events. At low densities the screening effect is reminiscent of the role of dipoles in dielectrics, in only renormalizing the elastic moduli. The effect of high density screening has no immediate electrostatic analog and is expected to change qualitatively the mechanical response, as seen for example in the displacement field. On the basis of experiments and simulations, we show that in granular matter, strong screening results in significant deviation from elasticity theory. The theoretical analysis, which accounts for an emergent inherent length scale, the experimental measurements and the numerical simulations of frictional granular amorphous assemblies are in agreement with each other, and provide a strong support for the novel continuum theory.

preprint2022arXiv

Universal Density of Low Frequency States in Silica Glass at Finite Temperatures

The theoretical understanding of the low-frequency modes in amorphous solids at finite temperature is still incomplete. The study of the relevant modes is obscured by the dressing of inter-particle forces by collision-induced momentum transfer that is unavoidable at finite temperatures. Recently, it was proposed that low frequency modes of vibrations around the {\em thermally averaged} configurations deserve special attention. In simple model glasses with bare binary interactions, these included quasi-localized modes whose density of states appears to be universal, depending on the frequencies as $D(ω) \sim ω^4$, in agreement with the similar law that is obtained with bare forces at zero temperature. In this work, we report investigations of a model of silica glass at finite temperature; here the bare forces include binary and ternary interactions. Nevertheless we can establish the validity of the universal law of the density of quasi-localized modes also in this richer and more realistic model glass.

preprint2021arXiv

Density of Quasi-localized Modes in Glasses: where are the Two-Level Systems?

The existence of a constant density of two-level systems (TLS) was proposed as the basis of some intriguing universal aspects of glasses at ultra-low temperatures. Here we ask whether their existence is necessary for explaining the universal density of states quasi-localized modes (QLM) in glasses at ultra-low temperatures. A careful examination of the QLM that exist in a generic atomistic model of a glass former reveals at least two types of them, each exhibiting a different density of states, one depending on the frequency as $ω^3$ and the other as $ω^4$. The properties of the glassy energy landscape that is responsible for the two types of modes is examined here, explaining the analytic feature responsible for the creations of (at least) two families of QLM's. Although adjacent wells certainly exist in the complex energy landscape of glasses, doubt is cast on the relevance of TLS for the universal density of QLM's.

preprint2021arXiv

Fatigue and Collapse of Cyclically Bent Strip of Amorphous Solid

Fatigue caused by cyclic bending of a piece of material, resulting in its mechanical failure, is a phenomenon that had been studied for ages by engineers and physicists alike. In this Letter we study such fatigue in a strip of athermal amorphous solid. On the basis of atomistic simulations we conclude that the crucial quantity to focus on is the {\em accumulated damage}. Although this quantity exhibits large sample-to-sample fluctuations, its dependence on the loading determines the statistics of the number of cycles to failure. Thus we can provide a scaling theory for the Wöhler plots of mean number of cycles for failure as a function of the loading amplitude.

preprint2021arXiv

Scaling Theory for Wöhler plots in amorphous solids under cyclic forcing

In mechanical engineering Wöhler plots serve to measure the average number of load cycles before materials break, as a function of the maximal stress in each cycle. Although such plots are prevalent in engineering for more than 150 years, their theoretical understanding is lacking. Recently a scaling theory of Wöhler plots in the context of cyclic bending was offered [1]. Here we elaborate further on cyclic bending and extend the considerations to cyclic tensile loads on an amorphous strip of material; the scaling theory applies to both types of cyclic loading equally well. On the basis of atomistic simulations we conclude that the crucial quantities to focus on are the accumulated damage and the average damage per cycle. The dependence of these quantities on the loading determines the statistics of the number of cycles to failure. Finally we consider the probability distribution functions of the number of cycles to failure and demonstrate that the scaling theory allows prediction of these distributions at one value of the forcing amplitude from measurements and another value.

preprint2020arXiv

Fatigue and failure of a polymer chain under tension

The rupture of a polymer chain maintained at temperature $T$ under fixed tension is prototypical to a wide array of systems failing under constant external strain and random perturbations. Past research focused on analytic and numerical studies of the mean rate of collapse of such a chain. Surprisingly, an analytic calculation of the probability distribution function (PDF) of collapse rates appears to be lacking. Since rare events of rapid collapse can be important and even catastrophic, we present here a theory of this distribution, with a stress on its tail of fast rates. We show that the tail of the PDF is a power law with a {\em universal} exponent that is theoretically determined. Extensive numerics validate the offered theory. Lessons pertaining to other problems of the same type are drawn.

preprint2020arXiv

Frictional Granular Matter: Protocol Dependence of Mechanical Properties

Theoretical treatments of frictional granular matter often assume that it is legitimate to invoke classical elastic theory to describe its coarse-grained mechanical properties. Here we show, based on experiments and numerical simulations, that this is generically not the case since stress auto-correlation functions decay more slowly than the elastic Green's function. It was shown theoretically that standard elastic decay demands pressure and torque density fluctuations to be normal, with possibly one of them being hyperuniform. Generic compressed frictional assemblies exhibit however abnormal pressure fluctuations, failing to conform with the central limit theorem. The physics of this failure is linked to correlations built in the material during compression from a dilute configuration prior to jamming. By changing the protocol of compression one can observe different pressure fluctuations and stress auto-correlations decay at large scales.

preprint2020arXiv

Giant Amplification of Small Perturbations in Frictional Amorphous Solid

Catastrophic events in Nature can be often triggered by small perturbations, with "remote triggering" of earthquakes being an important example. Here we present a mechanism for the giant amplification of small perturbations that is expected to be generic in systems whose dynamics is not derivable from a Hamiltonian. We offer a general discussion of the typical instabilities involved (being oscillatory with an exponential increase of noise) and examine in detail the normal forms that determine the relevant dynamics. The high sensitivity to external perturbations is explained for systems with and without dissipation. Numerical examples are provided using the dynamics of frictional granular matter. Finally we point out the relationship of the presently discussed phenomenon to the highly topical issue of "exceptional points" in quantum models with non-Hermitian Hamiltonians.

preprint2020arXiv

Plastic Instabilities in Charged Granular Systems: Competition between Elasticity and Electrostatics

Electrostatic theory preserves charges, but allows dipolar excitations. Elasticity theory preserves dipoles, but allows quadrupolar (Eshelby like) plastic events. Charged amorphous granular systems are interesting in their own right; here we focus on their plastic instabilities and examine their mechanical response to external strain and to external electric field, to expose the competition between elasticity and electrostatics. In this paper a generic model is offered, its mechanical instabilities are examined and a theoretical analysis is presented. Plastic instabilities are discussed as saddle-node bifurcations that can be fully understood in terms of eigenvalues and eigenfunctions of the relevant Hessian matrix. This system exhibits moduli that describe how electric polarization and stress are influenced by strain and electric field. Theoretical expression for these moduli are offered and compared to the measurements in numerical simulations.

preprint2020arXiv

Robustness of Density of Low Frequency States in Amorphous Solids

Low frequency quasi-localized modes of amorphous glasses appear to exhibit universal density of states, depending on the frequencies as $D(ω) \sim ω^4$. To date various models of glass formers with short range binary interaction, and network glasses with both binary and ternary interactions, were shown to conform with this law. In this paper we examine granular amorphous solids with long-range electrostatic interactions, and find that they exhibit the same law. To rationalize this wide universality class we return to a model proposed by Gurevich, Parshin and Schober (GPS) and analyze its predictions for interaction laws with varying spatial decay, exploring this wider than expected universality class. Numerical and analytic results are provided for both the actual system with long range interaction and for the GPS model.

preprint2020arXiv

Stress Correlations in Frictional Granular Media

This paper investigates whether in frictional granular packings, like in Hamiltonian amorphous elastic solids, the stress autocorrelation matrix presents long range anisotropic contributions just as elastic Green's functions. We find that in a standard model of frictional granular packing this is not the case. We prove quite generally that mechanical balance and material isotropy constrain the stress auto-correlation matrix to be fully determined by two spatially isotropic functions: the pressure and torque auto-correlations. The pressure and torque fluctuations being respectively normal and hyperuniform force the stress autocorrelation to decay as the elastic Green's function. Since we find the torque fluctuations to be hyper-uniform, the culprit is the pressure whose fluctuations decay slower than normally as a function of the system's size. Investigating the reason for these abnormal pressure fluctuations we discover that anomalous correlations build up already during the compression of the dilute system before jamming. Once jammed these correlations remain frozen. Whether this is true for frictional matter in general or is it the consequence of the model properties is a question that must await experimental scrutiny and possible alternative models.

preprint2020arXiv

The dimension of Diffusion Limited Aggregates grown on a line

Diffusion Limited Aggregation (DLA) has served for forty years as a paradigmatic example for the creation of fractal growth patterns. In spite of thousands of references no exact result for the fractal dimension $D$ of DLA is known. In this Letter we announce an exact result for off-lattice DLA grown on a line, $D=3/2$. The result relies on representing DLA with iterated conformal maps, allowing one to prove self-affinity, a proper scaling limit and a well defined fractal dimension. Mathematical proofs of the main results are available in arXiv:2008.05792.

preprint2020arXiv

Thermodynamic Equivalence of Cyclic Shear and Deep Cooling in Glass-Formers

The extreme slowing down associated with glass formation in experiments and in simulations results in serious difficulties to prepare deeply quenched, well annealed, glassy material. Recently, methods to achieve such deep quenching were proposed, including vapor deposition on the experimental side and "Swap Monte Carlo" and oscillatory shearing on the simulation side. The relation between the resulting glasses under different protocols remains unclear. Here we show that oscillatory shear and "Swap Monte Carlo" result in thermodynamically equivalent glasses sharing the same statistical mechanics and similar mechanical responses under external strain.

preprint2020arXiv

Universal Density of Low Frequency States in Amorphous Solids at Finite Temperatures

It has been established that the low frequency quasi-localized modes of amorphous solids at zero temperature exhibit universal density of states, depending on the frequencies as $D(ω) \sim ω^4$. It remains an open question whether this universal law extends to finite temperatures. In this Letter we show that well quenched model glasses at temperatures as high as $T_g/3$ possess the same universal density of states. The only condition required is that {\em average} particle positions stabilize before thermal diffusion destroys the cage structure of the material. The universal density of quasi-localized low frequency modes refers then to vibrations around the thermally averaged configuration of the material.

preprint2020arXiv

Universal low-frequency vibrational modes in silica glasses

It was recently shown that different simple models of glass formers with binary interactions define a universality class in terms of the density of states of their quasi-localized low-frequency modes. Explicitly, once the hybridization with standard Debye (extended) modes is avoided, a number of such models exhibit a universal density of state, depending on the mode frequencies as $D(ω) \sim ω^4$. It is unknown however how wide is this universality class, and whether it also pertains to more realistic models of glass formers. To address this issue we present analysis of the quasi-localized modes in silica, a network glass which has both binary and ternary interactions. We conclude that in 3-dimensions silica exhibits the very same frequency dependence at low frequencies, suggesting that this universal form is a generic consequence of amorphous glassiness.

preprint2019arXiv

Elementary plastic events in amorphous silica

Plastic instabilities in amorphous materials are often studied using idealized models of binary mixtures that do not capture accurately molecular interactions and bonding present in real glasses. Here we study atomic scale plastic instabilities in a three dimensional molecular dynamics model of silica glass under quasi-static shear. We identify two distinct types of elementary plastic events, one is a standard quasi-localized atomic rearrangement while the second is a bond breaking event that is absent in simplified models of fragile glass formers. Our results show that both plastic events can be predicted by a drop of the lowest non-zero eigenvalue of the Hessian matrix that vanishes at a critical strain. Remarkably, we find very high correlation between the associated eigenvectors and the non-affine displacement fields accompanying the bond breaking event, predicting the locus of structural failure. Both eigenvectors and non-affine displacement fields display an Eshelby-like quadrupolar structure for both failure modes, rearrangement or bond-breaking. Our results thus clarify the nature of atomic scale plastic instabilities in silica glasses providing useful information for the development of mesoscale models of amorphous plasticity.

preprint2019arXiv

Oscillatory Instabilities in 3-Dimensional Frictional Granular Matter

The dynamics of amorphous granular matter with frictional interactions cannot be derived in general from a Hamiltonian and therefore displays oscillatory instabilities stemming from the onset of complex eigenvalues in the stability matrix. These instabilities were discovered in the context of one and two dimensional systems, while the three dimensional case was never studied in detail. Here we fill this gap by deriving and demonstrating the presence of oscillatory instabilities in a three dimensional granular packing. We study binary assemblies of spheres of two sizes interacting via classical Hertz and Mindlin force laws for the longitudinal and tangent interactions, respectively. We formulate analytically the stability matrix in 3D and observe that a couple of complex eigenvalues emerges at the onset of the instability as in the case of frictional disks in two-dimensions. The dynamics then shows oscillatory exponential growth in the Mean-Square-Displacement, followed by a catastrophic event. The generality of these results for any choice of forces that break the symplectic Hamiltonian symmetry is discussed.

preprint2019arXiv

Transition from Static to Dynamic Friction in an Array of Frictional Disks

The nature of an instability that controls the transition from static to dynamical friction is studied in the the context of an array of frictional disks that are pressed from above on a substrate. In this case the forces are all explicit and Newtonian dynamics can be employed without any phenomenological assumptions. We show that an oscillatory instability that had been discovered recently is responsible for the transition, allowing individual disks to spontaneously reach the Coulomb limit and slide with dynamic friction. The transparency of the model allows a full understanding of the phenomenon, including the speeds of the waves that travel from the trailing to the leading edge and vice versa.

preprint2016arXiv

Breakdown of Nonlinear Elasticity in Amorphous Solids at Finite Temperatures

It is known by now that amorphous solids at zero temperature do not possess a nonlinear elasticity theory: besides the shear modulus which exists, all the higher order coefficients do not exist in the thermodynamic limit. Here we show that the same phenomenon persists up to temperatures comparable to the glass transition. The zero temperature mechanism due to the prevalence of dangerous plastic modes of the Hessian matrix is replaced by anomalous stress fluctuations that lead to the divergence of the variances of the higher order elastic coefficients. The conclusion is that in amorphous solids elasticity can never be decoupled from plasticity: the nonlinear response is very substantially plastic.

preprint2016arXiv

Determining the Inter-Particle Force-Laws in Amorphous Solids from a Visual Image

We consider the problem of how to determine the force laws in an amorphous system of interacting particles. Given the positions of the centers of mass of the constituent particles we propose a new algorithm to determine the inter-particle force-laws. Having $n$ different types of constituents we determine the coefficients in the Laurent polynomials for the $n(n+1)/2$ possibly different force-laws. A visual providing the particle positions in addition to a measurement of the pressure is all that is required. The algorithm proposed includes a part that can correct for experimental errors in the positions of the particles. Such a correction of unavoidable measurement errors is expected to benefit many experiments in the field.

preprint2016arXiv

Elasticity in Amorphous Solids: Nonlinear or Piece-Wise Linear?

Quasi-static strain-controlled measurements of stress vs strain curves in macroscopic amorphous solids result in a nonlinear looking curve that ends up either in mechanical collapse or in a steady-state with fluctuations around a mean stress that remains constant with increasing strain. It is therefore very tempting to fit a nonlinear expansion of the stress in powers of the strain. We argue here that at low temperatures the meaning of such an expansion needs to be reconsidered. We point out the enormous difference between quenched and annealed averages of the stress vs. strain curves, and propose that a useful description of the mechanical response is given by a stress (or strain) dependent shear modulus for which a theoretical evaluation exists. The elastic response is piece-wise linear rather than nonlinear.

preprint2016arXiv

Mechanical Yield in Amorphous Solids: a First-Order Phase Transition

Amorphous solids yield at a critical value of the strain (in strain controlled experiments); for larger strains the average stress can no longer increase - the system displays an elasto-plastic steady state. A long standing riddle in the materials community is what is the difference between the microscopic states of the material before and after yield. Explanations in the literature are material specific, but the universality of the phenomenon begs a universal answer. We argue here that there is no fundamental difference in the states of matter before and after yield, but the yield is a bona-fide first order phase transition between a highly restricted set of possible configurations residing in a small region of phase space to a vastly rich set of configurations which include many marginally stable ones. To show this we employ an order parameter of universal applicability, independent of the microscopic interactions, that is successful in quantifying the transition in an unambiguous manner.

preprint2016arXiv

Microalloying and the mechanical properties of amorphous solids

The mechanical properties of amorphous solids like metallic glasses can be dramatically changed by adding small concentrations (as low as 0.1\%) of foreign elements. The glass-forming-ability, the ductility, the yield stress and the elastic moduli can all be greatly effected. This paper presents theoretical considerations with the aim of explaining the magnitude of these changes in light of the small concentrations involved. The theory is built around the experimental evidence that the microalloying elements organize around them a neighborhood that differs from both the crystalline and the glassy phases of the material in the absence of the additional elements. These regions act as {\em isotropic} defects that in unstressed systems dress the shear moduli. When strained, these defects interact with the incipient plastic responses which are quadrupolar in nature. It will be shown that this interaction interferes with the creation of system-spanning shear bands and increases the yield strain. We offer experimentally testable estimates of the lengths of nano-shear bands in the presence of the additional elements.

preprint2016arXiv

Scaling Theory of the Mechanical Properties of Amorphous Nano-Films

Numerical Simulations are employed to create amorphous nano-films of a chosen thickness on a crystalline substrate which induces strain on the film. The films are grown by a vapor deposition technique which was recently developed to create very stable glassy films. Using the exact relations between the Hessian matrix and the shear and bulk moduli we explore the mechanical properties of the nano-films as a function of the density of the substrate and the film thickness. The existence of the substrate dominates the mechanical properties of the combined substrate-film system. Scaling concepts are then employed to achieve data collapse in a wide range of densities and film thicknesses.

preprint2016arXiv

Statistics of Plastic Events in Post-Yield Strain-Controlled Amorphous Solids

Amorphous solids yield in strain-controlled protocols at a critical value of the strain. For larger strains the stress and energy display a generic complex serrated signal with elastic segments punctuated by sharp energy and stress plastic drops having a wide range of magnitudes. Here we provide a theory of the scaling properties of such serrated signals taking into account the system-size dependence. We show that the statistics are not homogeneous - they separate sharply to a regime of `small' and `large' drops, each endowed with its own scaling properties. A scaling theory is first derived solely by data analysis, showing a somewhat complex picture. But after considering the physical interpretation one discovers that the scaling behavior and the scaling exponents are in fact very simple and universal.

preprint2016arXiv

The sandpile revisited: Computer assisted determination of constitutive relations and the breaking of scaling

We revisit the problem of the stress distribution in a frictional sandpile under gravity, equipped with a new numerical model of granular assemblies with both normal and tangential (frictional) inter-granular forces. Numerical simulations allow a determination of the spatial dependence of all the components of the stress field as a function of systems size, the coefficient of static friction and the frictional interaction with the bottom surface. Our study clearly demonstrates that interaction with the bottom surface plays a crucial role in the formation of a pressure dip under the apex of a granular pile. Basic to the theory of sandpiles are assumptions about the form of scaling solutions and constitutive relations for cohesive-less hard grains for which no typical scale is available. We find that these constitutive relations must be modified; moreover for smaller friction coefficients and smaller piles these scaling assumptions break down in the bulk of the sandpile due to the presence of length scales that must be carefully identified. After identifying the crucial scale we provide a predictive theory to when scaling solutions are expected to break down. At the bottom of the pile the scaling assumption always breaks, due to the different interactions with the bottom surface. The consequences for measurable quantities like the pressure distribution and shear stress at the bottom of the pile are discussed. For example one can have a transition from no dip in the base-pressure to a dip at the center of the pile as a function of the system size.

preprint2016arXiv

What determines static force chains in stressed granular media?

The determination of the normal and transverse (frictional) inter-particle forces within a granular medium is a long standing, taunting, and yet unresolved problem. We present a new formalism which employs the knowledge of the external forces and the orientations of contacts between particles (of any given sizes), to compute all the inter-particle forces. Having solved this problem we exemplify the efficacy of the formalism showing that the force chains in such systems are determined by an expansion in the eigenfunctions of a newly defined operator.

preprint2015arXiv

Consequences of Disorder on the Stability of Amorphous Solids

Highly acurate numerical simulations are employed to highlight the subtle but important differences in the mechanical stability of perfect crystalline solids versus amorphous solids. We stress the difference between strain values at which the shear modulus vanishes and strain values at which a plastic instability insues. The temperature dependence of the yield strain is computed for the two types of solids, showing different scaling laws: $γ_{_{Y}}\simeq γ_{_{Y}}^{0}-C_1 T^{1/3}$ for crystals versus $γ_{_{Y}}\simeq γ_{_{Y}}^{0}- C_2 T^{2/3}$ for amorphous solids.

preprint2015arXiv

Counter-flow Induced Decoupling in Super-Fluid Turbulence

In mechanically driven superfluid turbulence the mean velocities of the normal- and superfluid components are known to coincide: $\mathbf U_{\text{n}} =\mathbf U_{\text{s}}$. Numerous laboratory, numerical and analytical studies showed that under these conditions the mutual friction between the normal- and superfluid velocity components couples also their fluctuations: $\mathbf u'_{\text{n}}(\mathbf r,t) \approx \mathbf u'_{\text{s}}(\mathbf r,t)$ almost at all scales. In this paper we show that this is not the case in thermally driven superfluid turbulence; here the counterflow velocity $\mathbf U_{\text{ns}}\equiv \mathbf U_{\text{n}} -\mathbf U_{\text{s}}\ne 0$. We suggest a simple analytic model for the cross correlation function $\left\langle \mathbf u'_{\text{n}}(\mathbf r,t) \cdot \mathbf u'_{\text{s}}(\mathbf r',t)\right \rangle$ and its dependence on $U_{\text{ns}}$. We demonstrate that $\mathbf u'_{\text{n}}(\mathbf r,t)$ and $ \mathbf u'_{\text{s}}(\mathbf r,t)$ are decoupled almost in the entire range of separations $|\mathbf r-\mathbf r'|$ between the energy containing scale and intervortex distance.

preprint2015arXiv

Energy and Vorticity Spectra in Turbulent Superfluid $^4$He from $T=0$ to $T_λ$

We discuss the energy and vorticity spectra of turbulent superfluid $^4$He in all the temperature range from $T=0$ up to the phase transition "$λ$ point", $T_λ\simeq 2.17\,$K. Contrary to classical developed turbulence in which there are only two typical scales, i.e. the energy injection $L$ and the dissipation scales $η$, here the quantization of vorticity introduces two additional scales, i.e the vortex core radius $a_0$ and the mean vortex spacing $\ell$. We present these spectra for the super- and normal-fluid components in the entire range of scales from $L$ to $a_0$ including the cross-over scale $\ell$ where the hydrodynamic eddy-cascade is replaced by the cascade of Kelvin waves on individual vortices. At this scale a bottleneck accumulation of the energy was found earlier at $T=0$. We show that even very small mutual friction dramatically suppresses the bottleneck effect due to the dissipation of the Kelvin waves. Using our results for the spectra we estimate the Vinen "effective viscosity" $ν'$ in the entire temperature range and show agreement with numerous experimental observation for $ν'(T)$.

preprint2015arXiv

Modeling Barkhausen Noise in Magnetic Glasses with Dipole-Dipole Interactions

Long-ranged dipole-dipole interactions in magnetic glasses give rise to magnetic domains having labyrinthine patterns. Barkhausen Noise is then expected to result from the movement of domain boundaries which is supposed to be modeled by the motion of elastic membranes with random pinning. We propose an atomistic model of such magnetic glasses in which we measure the Barkhausen Noise which indeed results from the movement of domain boundaries. Nevertheless the statistics of the Barkhausen Noise is found in striking disagreement with the expectations in the literature. In fact we find exponential statistics without any power law, stressing the fact that Barkhausen Noise can belong to very different universality classes. In this glassy system the essence of the phenomenon is the ability of spin-carrying particles to move and minimize the energy without any spin flip. A theory is offered in excellent agreement with the measured data without any free parameter.

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

The Anatomy of Plastic Events in Magnetic Amorphous Solids

Plastic events in amorphous solids can be much more than just "shear transformation zones" when the positional degrees of freedom are coupled non-trivially to other degrees of freedom. Here we consider magnetic amorphous solids where mechanical and magnetic degrees of freedom interact, leading to rather complex plastic events whose nature must be disentangled. In this paper we uncover the anatomy of the various contributions to some typical plastic events. These plastic events are seen as Barkhausen Noise or other "serrated noises". Using theoretical considerations we explain the observed statistics of the various contributions to the considered plastic events. The richness of contributions and their different characteristics imply that in general the statistics of these "serrated noises" cannot be universal, but rather highly dependent on the state of the system and on its microscopic interactions.

preprint2015arXiv

The yield of Amorphous Solids Under Stress Control at Low Temperatures

The yield of amorphous solids like metallic glasses under external stress was discussed asserting that it is related to the glass transition by increasing temperature, or that it can be understood using statistical theories of various sorts. Here we study the approach to stress-controlled yield and argue that neither assertions can be supported, at least at low temperatures. The yield of amorphous solids at low temperatures is a highly structured phenomenon, characterized by a specific series of mechanical instabilities, and having no similarity at all to fluidization by increased temperature, real or fictive. The series of instabilities followed by stress controlled yield at low but finite temperature protocols can be predicted by analyzing athermal quasi-static strain controlled protocols, making the latter highly relevant for the deep understanding of the mechanical properties of amorphous solids.

preprint2014arXiv

Elasticity and Plasticity in Stiff and Flexible Oligomeric Glasses

In this paper we focus on the mechanical properties of oligomeric glasses (waxes), employing a microscopic model that provides, via numerical simulations, information about the shear modulus of such materials, the failure mechanism via plastic instabilities and about the geometric responses of the oligomers themselves to a mechanical load. We present a microscopic theory that explains the numerically observed phenomena, including an exact theory of the shear modulus and of the plastic instabilities, both local and system spanning. In addition we present a model to explain the geometric changes in the oligomeric chains under increasing strains.

preprint2014arXiv

Evolution of Neutron-Initiated Micro-Big-Bang in superfluid He 3B

A nuclear capture reaction of a single neutron by ultra-cold superfluid $^3$He results in a rapid overheating followed by the expansion and subsequent cooling of the hot subregion, in a certain analogy with the Big Bang of the early Universe. It was shown in a Grenoble experiment that a significant part of the energy released during the nuclear reaction was not converted into heat even after several seconds. It was thought that the missing energy was stored in a tangle of quantized vortex lines. This explanation, however, contradicts the expected lifetime of a bulk vortex tangle, $10^{-5}-10^{-4}\,$s, which is much shorter than the observed time delay of seconds. In this Letter we propose a scenario that resolves the contradiction: the vortex tangle, created by the hot spot, emits isolated vortex loops that take with them a significant part of the tangle's energy. These loops quickly reach the container walls. The dilute ensemble of vortex loops attached to the walls can survive for a long time, while the remaining bulk vortex tangle decays quickly.

preprint2014arXiv

Mechanical Properties and Plasticity of a Model Glass Loaded Under Stress Control

Much of the progress achieved in understanding plasticity and failure in amorphous solids had been achieved using experiments and simulations in which the materials were loaded using strain control. There is paucity of results under stress control. Here we present a new method that was carefully geared to allow loading under stress control either at $T=0$ or at any other temperature, using Monte-Carlo techniques. The method is applied to a model perfect crystalline solid, to a crystalline solid contaminated with topological defects, and to a generic glass. The highest yield stress belongs to the crystal, the lowest to the crystal with a few defects, with the glass in between. Although the glass is more disordered than the crystal with a few defects, it yields stress is much higher than that of the latter. We explain this fact by considering the actual microscopic interactions that are typical to glass forming materials, pointing out the reasons for the higher cohesive nature of the glass. The main conclusion of this paper is that the instabilities encountered in stress-control condition are the identical saddle-node bifurcation seen in strain-control. Accordingly one can use the latter condition to infer about the former. Finally we discuss temperature effects and comment on the time needed to see a stress controlled material failure.

preprint2014arXiv

Shear Transformation Zones: State Determined or Protocol Dependent?

The concept of a Shear Transformation Zone (STZ) refers to a region in an amorphous solid that undergoes a plastic event when the material is put under an external mechanical load. An important question that had accompanied the development of the theory of plasticity in amorphous solids for many years now is whether an STZ is a {\em region} existing in the material (which can be predicted by analyzing the unloaded material), or is it an {\em event} that depends on the loading protocol (i.e., the event cannot be predicted without following the protocol itself). In this Letter we present strong evidence that the latter is the case. Infinitesimal changes of protocol result in macroscopically big jumps in the positions of plastic events, meaning that these can never be predicted from considering the unloaded material.

preprint2014arXiv

Statistical Physics of Pure Barkhausen Noise

We discuss a model metallic glass in which Barkhausen Noise can be studied in exquisite detail, free of thermal effects and of the rate of ramping of the magnetic field. The mechanism of the jumps in magnetic moment that cause the Barkhausen Noise can be fully understood as consecutive instabilities where an eigenvalue of the Hessian matrix hits zero, leading to a magnetization jump $Δm$ which is simultaneous with a stress and energy changes $Δσ$ and $ΔU$ respectively. Contrary to common belief we find no "movements of magnetic domain boundaries" across pinning sites, no fractal domains, no self-organized criticality and no exact scaling behaviour. We present a careful numerical analysis of the statistical properties of the phenomenon, and show that with every care taken this analysis is tricky, and easily misleading. Without a guiding theory it is almost impossible to get the right answer for the statistics of Barkhausen Noise. We therefore present an analytic theory, showing that the probability distribution function (pdf) of Barkhausen Noise is not a power law times an exponential cutoff.

preprint2014arXiv

The Glass Transition in Fluids with Magnetic Interactions

We study the glass transition in fluids where particles are endowed with spins, such that magnetic and positional degrees of freedom are coupled. Novel results for slowing down in the spin time-correlation functions are described, and the effects of magnetic fields on the glass transition are studied. Aging effects in such systems and the corresponding data collapse are presented and discussed.

preprint2013arXiv

Analytic solution of the dynamics of quantum vortex reconnection

Experimental and simulational studies of the dynamics of vortex reconnections in quantum fluids showedthat the distance $d$ between the reconnecting vortices is close to a universal time dependence $d=D[κ|t_0-t|]^α$ with $α$ fluctuating around 1/2 and $κ=h/m$ is the quantum of circulation. Dimensional analysis, based on the assumption that the quantum of circulation $κ=h/m$ is the only relevant parameter in the problem, predicts $α=1/2$. The theoretical calculation of the dimensionless coefficient $D$ in this formula remained an open problem. In this Letter we present an analytic calculation of $D$ in terms of the given geometry of the reconnecting vortices. We start from the numerically observed generic geometry on the way to vortex reconnection and demonstrate that the dynamics is well described by a self-similar analytic solution which provides the wanted information.

preprint2013arXiv

Atomistic Simulations of Magnetic Amorphous Solids: Magnetostriction, Barkhausen noise and novel singularities

We present results of atomistic simulations of a new model of a magnetic amorphous solid subjected to external mechanical strains and magnetic fields. The model employed offers new perspectives on important effects like Barkhausen noise and magnetostriction. It is shown that the plastic response in such systems exhibit singularities characterized by unexpected exponents requiring careful theoretical reasoning. The spatial structure of the plastic events requires a new coarse grained elasto-magnetic theory which is provided here.

preprint2013arXiv

Atomistic theory of the shear band direction in amorphous solids

One of the major theoretical riddles in shear banding instabilities is the angle that the shear band chooses spontaneously with respect to the principal stress axis. Here we employ our recent atomistic theory to compute analytically the angle in terms of the characteristics of the Eshelby inclusion that models faithfully the eigenfunction of the Hessian matrix that goes soft at the plastic instability. We show that loading protocols that do not conserve volume result in shear bands at angles different from 45$^o$ to the strain axis; only when the external strains preserve volume like in pure shear, the shear bands align precisely at 45$^o$ to the strain axis. We compute an analytic formula for the angle of the shear band in terms of the characteristics of the loading protocol; quantitative agreement with computer simulations is demonstrated.

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.

preprint2013arXiv

Fokker-Planck equation with memory: the crossover from ballistic to diffusive processes in many-particle systems and incompressible media

The unified description of diffusion processes that cross over from a ballistic behavior at short times to normal or anomalous diffusion (sub- or superdiffusion) at longer times is constructed on the basis of a non-Markovian generalization of the Fokker-Planck equation. The necessary non-Markovian kinetic coefficients are determined by the observable quantities (mean- and mean square displacements). Solutions of the non-Markovian equation describing diffusive processes in the physical space are obtained. For long times, these solutions agree with the predictions of the continuous random walk theory; they are, however, much superior at shorter times when the effect of the ballistic behavior is crucial.

preprint2013arXiv

Micro-alloying and the Toughness of Glasses: Modeling with Pinned Particles

The usefulness of glasses, and particularly of metallic glasses, in technological applications is often limited by their toughness, which is defined as the area under the stress vs. strain curve before plastic yielding. Recently toughness was found to increase significantly by the addition of small concentrations of foreign atoms that act as pinning centers. We model this phenomenon at zero temperature and quasi-static straining with randomly positioned particles that participate in the elastic deformation but are pinned in the non-affine return to mechanical equilibrium. We find a very strong effect on toughness via the increase of both the shear modulus and the yield stress as a function of the density of pinned particles. Understanding the results calls for analyzing separately the elastic, or "Born term" and the contributions of the "excess modes" that result from glassy disorder. Finally we present a scaling theory that collapses the data on one universal curve as a function of rescaled variables.

preprint2013arXiv

On the Effect of Micro-alloying on the Mechanical Properties of Metallic Glasses

"Micro-alloying", referring to the addition of small concentration of a foreign metal to a given metallic glass, was used extensively in recent years to attempt to improve the mechanical properties of the latter. The results are haphazard and nonsystematic. In this paper we provide a microscopic theory of the effect of micro-alloying, exposing the delicate consequences of this procedure and the large parameter space which needs to be controlled. In particular we consider two very similar models which exhibit opposite trends for the change of the shear modulus, and explain the origins of the difference as displayed in the different microscopic structure and properties.

preprint2013arXiv

Percolating Plastic Failure as a Mechanism for Shear Softening in Amorphous Solids

``Shear softening" refers to the observed reduction in shear modulus when the stress on an amorphous solid is increased beyond the initial linear region. Careful numerical quasi-static simulations reveal an intimate relation between plastic failure and shear softening. The attaintment of the steady-state value of the shear modulus associated with plastic flow is identified with a percolation of the regions that underwent a plastic event. We present an elementary ``two-state" model that interpolates between failed and virgin regions and provides a simple and effective characterization of the shear softening.

preprint2013arXiv

Plasticity-Induced Magnetization in Amorphous Magnetic Solids

Amorphous magnetic solids, like metallic glasses, exhibit a novel effect: the growth of magnetic order as a function of mechanical strain under athermal conditions in the presence of a magnetic field. The magnetic moment increases in steps whenever there is a plastic event. Thus plasticity induces the magnetic ordering, acting as the effective noise driving the system towards equilibrium. We present results of atomistic simulations of this effect in a model of a magnetic amorphous solid subjected to pure shear and a magnetic field. To elucidate the dependence on external strain and magnetic field we offer a mean-field theory that provides an adequate qualitative understanding of the observed phenomenon.

preprint2013arXiv

Shear localization in 3-Dimensional Amorphous Solids

In this paper we extend the recent theory of shear-localization in 2-dimensional amorphous solids to 3-D. In 2-D the fundamental instability of shear-localization is related to the appearance of a line of displacement quadrupoles, that makes an angle of 45 degrees with the principal stress axis. In 3-D the fundamental plastic instability is also explained by the formation of a lattice of anisotropic elastic inclusions. In the case of pure external shear stress, we demonstrate that this is a 2-dimensional triangular lattice of similar elementary events. It is shown that this lattice is arranged on a plane, that, similarly to the 2-D case, makes an angle of 45 degrees with respect to the principal stress axis. This solution is energetically favorable only if the external strain exceeds a yield-strain value, which is determined by the strain parameters of the elementary events and the Poisson ratio. The predictions of the theory are compared to numerical simulations and very good agreement is observed.

preprint2013arXiv

Structure of quantum vortex tangle in He-4 counterflow turbulence

The main goal of this paper is to present a comprehensive characterization of well developed vortex tangles in a turbulent counterflow in quantum fluids (with a laminar normal fluid component). We analyze extensive numerical simulations using the vortex filament method, solving the full Biot-Savart equations for the vortex dynamics in a wide range of temperatures and counter-flow velocities. In addition to a detailed analysis of traditional characteristics such as vortex line density, anisotropic and curvature parameters of the vortex tangle, we stress other dynamical and statistical characteristics which are either much less studied or even unstudied. The latter include reconnection rates, mean mutual friction forces, drift velocities and the probability distribution functions of various tangle parameters: the loop length, the line curvature, the mean curvature of loops with a given length, etc. During these studies we compare the three main reconnection procedures which are widely used in the literature, and identify which properties are strongly affected by the choice of the reconnection criteria and which of them are practically insensitive to the reconnection procedure. The conclusion is that the vortex filament method in the framework of the Biot-Savart equation sufficiently robust and well suited for the description of the steady state vortex tangle in a quantum counterflow. The Local-Induction Approximation to this equation may be successfully used to analytically establish relationships between mean characteristics of the stochastic vortex tangle.

preprint2013arXiv

The cooling rate dependence of the shear modulus of amorphous solids

Rapidly cooling a liquid may result in a glass transition, creating an amorphous solid whose shear and bulk moduli are finite. Even when done with constant density, these resulting moduli depend strongly on the rate of cooling. Understanding this phenomenon calls for analyzing separately the "Born term" that exists also in perfectly ordered materials and the contributions of the "excess modes" that result from glassy disorder. We show that the Born term is very insensitive to the cooling rate, and all the variation in the shear modulus is due to the excess modes. We argue that this approach provides a quantitative understanding of the cooling rate dependence of a basic linear response coefficient, i.e. the shear modulus.

preprint2013arXiv

The Yield-Strain and Shear-Band Direction in Amorphous Solids Under General Loading

It is well known experimentally that well-quenched amorphous solids exhibit a plastic instability in the form of a catastrophic shear localization at a well defined value of the external strain. The instability may develop to a shear-band that in some cases is followed by a fracture. It is also known that the values of the yield-strain (and yield-stress), as well as the direction of the shear band with respect to the principal stress axis, vary considerably with variations in the external loading conditions. In this paper we present a microscopic theory of these phenomena for 2-dimensional athermal amorphous solids that are strained quasi-statically. We present analytic formulae for the yield-strains for different loading conditions, and well as for the angles of the shear bands. We explain that the external loading conditions determine the eigenvalues of the quadrupolar Eshelby inclusions which model the non-affine displacement field. These inclusions model elementary plastic events and determine both the yield-strain and the direction of the shear-band. We show that the angles of the shear bands with respect to the principal stress axis are limited theoretically between $30^o$ and $60^o$. Available experimental data conform to this prediction.

preprint2012arXiv

Analytical modeling for the heat transfer in sheared flows of nanofluids

We developed a model for the enhancement of the heat flux by spherical and elongated nano- particles in sheared laminar flows of nano-fluids. Besides the heat flux carried by the nanoparticles the model accounts for the contribution of their rotation to the heat flux inside and outside the particles. The rotation of the nanoparticles has a twofold effect, it induces a fluid advection around the particle and it strongly influences the statistical distribution of particle orientations. These dynamical effects, which were not included in existing thermal models, are responsible for changing the thermal properties of flowing fluids as compared to quiescent fluids. The proposed model is strongly supported by extensive numerical simulations, demonstrating a potential increase of the heat flux far beyond the Maxwell-Garnet limit for the spherical nanoparticles. The road ahead which should lead towards robust predictive models of heat flux enhancement is discussed.

preprint2012arXiv

Derivation of the Johnson-Samwer $T^{(2/3)}$ Temperature Dependence of the Yield Strain in Metallic Glasses

Metallic Glasses are prone to fail mechanically via a shear-banding instability. In a remarkable paper Johnson and Samwer demonstrated that this failure enjoys a high degree of universality in the sense that a large group of metallic glasses appears to possess a yield-strain that decreases with temperature following a $-T^{2/3}$ law up to logarithmic corrections. In this Letter we offer a theoretical derivation of this law. We show that our formula fits very well simulational data on typical amorphous solids.

preprint2012arXiv

Elastic Moduli in Nano-Size Samples of Amorphous Solids: System Size Dependence

This Letter is motivated by some recent experiments on pan-cake shaped nano-samples of metallic glass that indicate a decline in the measured shear modulus upon decreasing the sample radius. Similar measurements on crystalline samples of the same dimensions showed a much more modest change. In this Letter we offer a theory of this phenomenon; we argue that such results are generically expected for any amorphous solid, with the main effect being related to the increased contribution of surfaces with respect to bulk when the samples get smaller. We employ exact relations between the shear modulus and the eigenvalues of the system's Hessian matrix to explore the role of surface modes in affecting the elastic moduli.

preprint2012arXiv

Enhancement of intermittency in superfluid turbulence

We consider the intermittent behavior of superfluid turbulence in $^4$He. Due to the similarity in the nonlinear structure of the two-fluid model of superfluidity and the Euler and Navier-Stokes equations one expects the scaling exponents of the structure functions to be the same as in classical turbulence for temperatures close to the superfluid transition $T_λ$ and also for $T\ll T_λ$. This is not the case when mutual friction becomes important. Using shell model simulations, we propose that for an intermediate regime of temperatures, such that the density of normal and superfluid components are comparable to each other, there exists a range of scales in which the effective exponents indicate stronger intermittency. We offer a bridge relation between these effective and the classical scaling exponents. Since this effect occurs at accessible temperatures and Reynolds numbers, we propose that experiments should be conducted to further assess the validity and implications of this prediction.

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

Microscopic Mechanism of Shear Bands in Amorphous Solids

The fundamental instability responsible for the shear localization which results in shear bands in amorphous solids remains unknown despite enormous amount of research, both experimental and theoretical. As this is the main mechanism for the failure of metallic glasses, understanding the instability is invaluable in finding how to stabilize such materials against the tendency to shear localize. In this Letter we explain the mechanism for shear localization under shear, which is the appearance of highly correlated lines of Eshelby-like quadrupolar singularities which organize the non-affine plastic flow of the amorphous solid into a shear band. We prove analytically that such highly correlated solutions in which $\C N$ quadrupoles are aligned with equal orientations are minimum energy states when the strain is high enough. The line lies at 45 degrees to the compressive stress.

preprint2012arXiv

Propagation Mechanism of Brittle Cracks

We employ a recently developed model that allows the study of two-dimensional brittle crack propagation under fixed grip boundary conditions. The crack development highlights the importance of voids which appear ahead of the crack as observed in experiments on the nano-scale. The appearance of these voids is responsible for roughening the crack path on small scales, in agreement with theoretical expectations. With increasing speed of propagation one observes the branching instabilities that were reported in experiments. The simulations allow understanding the phenomena by analyzing the elastic stress field that accompanies the crack dynamics.

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

Temperature suppression of Kelvin-wave turbulence in superfluids

Kelvin waves propagating on quantum vortices play a crucial role in the phenomenology of energy dissipation of superfluid turbulence. Previous theoretical studies have consistently focused on the zero-temperature limit of the statistical physics of Kelvin-wave turbulence. In this letter, we go beyond this athermal limit by introducing a small but finite temperature in the form of non-zero mutual friction dissipative force; A situation regularly encountered in actual experiments of superfluid turbulence. In this case we show that there exists a new typical length-scale separating a quasi-inertial range of Kelvin wave turbulence from a far dissipation range. The letter culminates with analytical predictions for the energy spectrum of the Kelvin-wave turbulence in both of these regimes.

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.

preprint2012arXiv

The Plastic Response of Magnetoelastic Amorphous Solids

We address the cross effects between mechanical strains and magnetic fields on the plastic response of magnetoelastic amorphous solids. It is well known that plasticity in non-magnetic amorphous solids under external strain $γ$ is dominated by the co-dimension 1 saddle-node bifurcation in which an eigenvalue of the Hessian matrix vanishes at $γ_P$ like $\sqrt{γ_P-γ}$. This square-root singularity determines much of the statistical physics of elasto-plasticity, and in particular that of the stress-strain curves under athermal-quasistatic conditions. In this Letter we discuss the much richer physics that can be expected in magnetic amorphous solids. Firstly, magnetic amorphous solids exhibit co-dimension 2 plastic instabilities, when an external strain and an external magnetic field are applied simultaneously. Secondly, the phase diagrams promise a rich array of new effects that have been barely studied; this opens up a novel and extremely rich research program for magnetoplastic materials.

preprint2012arXiv

The Yield-Strain in Shear Banding Amorphous Solids

In recent research it was found that the fundamental shear-localizing instability of amorphous solids under external strain, which eventually results in a shear band and failure, consists of a highly correlated array of Eshelby quadrupoles all having the same orientation and some density $ρ$. In this paper we calculate analytically the energy $E(ρ,γ)$ associated with such highly correlated structures as a function of the density $ρ$ and the external strain $γ$. We show that for strains smaller than a characteristic strain $γ_Y$ the total strain energy initially increases as the quadrupole density increases, but that for strains larger than $γ_Y$ the energy monotonically decreases with quadrupole density. We identify $γ_Y$ as the yield strain. Its value, derived from values of the qudrupole strength based on the atomistic model, agrees with that from the computed stress-strain curves and broadly with experimental results.

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

Energy Spectra of Superfluid Turbulence in $^3$He

In superfluid $^3$He turbulence is carried predominantly by the superfluid component. To explore the statistical properties of this quantum turbulence and its differences from the classical counterpart we adopt the time-honored approach of shell models. Using this approach we provide numerical simulations of a Sabra-shell model that allows us to uncover the nature of the energy spectrum in the relevant hydrodynamic regimes. These results are in qualitative agreement with analytical expressions for the superfluid turbulent energy spectra that were found using a differential approximation for the energy flux.

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

Flame front propagation I: The Geometry of Developing Flame Fronts: Analysis with Pole Decomposition

The roughening of expanding flame fronts by the accretion of cusp-like singularities is a fascinating example of the interplay between instability, noise and nonlinear dynamics that is reminiscent of self-fractalization in Laplacian growth patterns. The nonlinear integro-differential equation that describes the dynamics of expanding flame fronts is amenable to analytic investigations using pole decomposition. This powerful technique allows the development of a satisfactory understanding of the qualitative and some quantitative aspects of the complex geometry that develops in expanding flame fronts.

preprint2011arXiv

Flame front propagation II: Random Noise and Pole-Dynamics in Unstable Front Propagation

The problem of flame propagation is studied as an example of unstable fronts that wrinkle on many scales is studied. The analytic tool of pole expansion in the complex plane is emloyed to address the interaction of the unstable growth process with random initial conditions and perturbations. We argue that the effect of random noise is immense and that it can never be neglected in sufficiently large systems. We present simulations that lead to scaling laws for the velocity and acceleration of the front as a function of the system size and the level of noise, and analytic arguments that explain these results in terms of the noisy pole dynamics.

preprint2011arXiv

Flame front propagation III: Random Noise and Pole-Dynamics in Unstable Front Propagation (new version)

The problem of flame propagation is studied as an example of unstable fronts that wrinkle on many scales. The analytic tool of pole expansion in the complex plane is employed to address the interaction of the unstable growth process with random initial conditions and perturbations. We argue that the effect of random noise is immense and that it can never be neglected in sufficiently large systems. We present simulations that lead to scaling laws for the velocity and acceleration of the front as a function of the system size and the level of noise, and analytic arguments that explain these results in terms of the noisy pole dynamics.terms of the dynamics of singularities in the complex plane, yielding detailed understanding of the physics of the eigenfunctions and eigenvalues.

preprint2011arXiv

Flame front propagation V: Stability Analysis of Flame Fronts: Dynamical Systems Approach in the Complex Plane

We consider flame front propagation in channel geometries. The steady state solution in this problem is space dependent, and therefore the linear stability analysis is described by a partial integro-differential equation with a space dependent coefficient. Accordingly it involves complicated eigenfunctions. We show that the analysis can be performed to required detail using a finite order dynamical system in terms of the dynamics of singularities in the complex plane, yielding detailed understanding of the physics of the eigenfunctions and eigenvalues.

preprint2011arXiv

Flame front propagation VI: Dynamics and Wrinkling of Radially Propagating Fronts Inferred from Scaling Laws in Channel Geometries

Flame Propagation is used as a prototypical example of expanding fronts that wrinkle without limit in radial geometries but reach a simple shape in channel geometry. We show that the relevant scaling laws that govern the radial growth can be inferred once the simpler channel geometry is understood in detail. In radial geometries (in contrast to channel geometries) the effect of external noise is crucial in accelerating and wrinkling the fronts. Nevertheless, once the interrelations between system size, velocity of propagation and noise level are understood in channel geometry, the scaling laws for radial growth follow.

preprint2011arXiv

Kelvin-wave turbulence in superfluids

We study the statistical and dynamical behavior of turbulent Kelvin waves propagating on quantized vortices in superfluids, and address the controversy concerning the energy spectrum that is associated with these excitations. Finding the correct energy spectrum is important because Kelvin waves play a major role in the dissipation of energy in superfluid turbulence at near-zero temperatures. In this paper, we show analytically that the solution proposed in Ref. \cite{10LN} enjoys existence, uniqueness and regularity of the pre-factor. Furthermore, we present numerical results of the dynamical equation that describes to leading order the non-local regime of the Kelvin wave dynamics. We compare our findings with the analytical results from the proposed local and non-local theories for Kelvin wave dynamics and show an agreement with the non-local predictions. Accordingly, the spectrum proposed in Ref. \cite{10LN} should be used in future theories of quantum turbulence. Finally, for weaker wave forcing we observe an intermittent behavior of the wave spectrum with a fluctuating dissipative scale, which we interpreted as a finite-size effect characteristic to mesoscopic wave turbulence.

preprint2011arXiv

Statistical Mechanics of Glass Formation in Molecular Liquids with OTP as an Example

We extend our statistical mechanical theory of the glass transition from examples consisting of point particles to molecular liquids with internal degrees of freedom. As before, the fundamental assertion is that super-cooled liquids are ergodic, although becoming very viscous at lower temperatures, and are therefore describable in principle by statistical mechanics. The theory is based on analyzing the local neighborhoods of each molecule, and a statistical mechanical weight is assigned to every possible local organization. This results in an approximate theory that is in very good agreement with simulations regarding both thermodynamical and dynamical properties.

preprint2011arXiv

Statistical Mechanics of the Anomalous Behavior of Tetrahedral Liquids

Tetrahedral liquids such as water and silica-melt show unusual thermodynamic behavior such as a density maximum and an increase in specific-heat when cooled to low temperatures. There is a debate in the literature whether these phenomena stem from a phase transition into a low-density and high-density liquid phases, which occur in the supercooled regime. Here we consider a model of tetrahedral liquids for which we construct a volume-constrained statistical mechanical theory which quantifies the local structure of the liquid. We compare the theory to molecular dynamics simulations and show that the theory can rationalize the simulations semi-quantitatively. We show that the anomalous density and specific heat behavior arise naturally from this theory without exhibiting a liquid-liquid phase-transition. We explain that this theory may or may not have a phase transition, depending on the volume and temperature dependence of the energy and entropy which are sensitive to small changes in the parameters of the model.

preprint2011arXiv

Turbulence in non-integer dimensions by fractal Fourier decimation

Fractal decimation reduces the effective dimensionality of a flow by keeping only a (randomly chosen) set of Fourier modes whose number in a ball of radius $k$ is proportional to $k^D$ for large $k$. At the critical dimension D=4/3 there is an equilibrium Gibbs state with a $k^{-5/3}$ spectrum, as in [V. L'vov {\it et al.}, Phys. Rev. Lett. {\bf 89}, 064501 (2002)]. Spectral simulations of fractally decimated two-dimensional turbulence show that the inverse cascade persists below D=2 with a rapidly rising Kolmogorov constant, likely to diverge as $(D-4/3)^{-2/3}$.

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

Density Scaling of Avalanche Statistics in Amorphous Solids

Stress vs. strain fluctuations in athermal amorphous solids are an example of `crackling noise' of the type studied extensively in the context of elastic membranes moving through random potentials. Contrary to the latter, we do not have a stochastic equation whose statistics agree with the measured ones. On the other hand we show in this Letter that the statistics of the former exhibit 'density scaling' when the interparticle potential can be well approximated by a power law. The distributions of sizes of dissipative events for a wide range of densities and system sizes follow a single scaling function. We find that both the system-size scaling of energy drops and the entire strain interval statistics are invariant to changes in density. We use this to determine accurately the exponents in the scaling laws, establishing that the present crackling noise belongs to a different universality class.

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

Stress Intensity Factor of Mode III Cracks in Thin Sheets

The stress field at the tip of a crack of a thin plate of elastic material that is broken due to a mode III shear tearing has a universal form with a non-universal amplitude, known as the stress intensity factor, which depends on the crack length and the boundary conditions. We present in this paper exact analytic results for this stress intensity factor, thus enriching the small number of exact results that can be obtained within Linear Elastic Fracture Mechanics (LEFM).

preprint2010arXiv

The dynamics of cracks in torn thin sheets

Motivated by recent experiments, we present a study of the dynamics of cracks in thin sheets. While the equations of elasticity for thin plates are well known, there remains the question of path selection for a propagating crack. We invoke a generalization of the principle of local symmetry to provide a criterion for path selection and demonstrate qualitative agreement with the experimental findings. The nature of the singularity at the crack tip is studied with and without the interference of nonlinear terms.

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

Drying Patterns: Sensitivity to Residual Stresses

Volume alteration in solid materials is a common cause of material failure. Here we investigate the crack formation in thin elastic layers attached to a substrate. We show that small variations in the volume contraction and substrate restraint can produce widely different crack patterns ranging from spirals to complex hierarchical networks. The networks are formed when there is no prevailing gradient in material contraction whereas spirals are formed in the presence of a radial gradient in the contraction of a thin elastic layer.

preprint2009arXiv

Finite-Dimensional Turbulence of Planetary Waves

Finite-dimensional wave turbulence refers to the chaotic dynamics of interacting wave `clusters' consisting of finite number of connected wave triads with exact three-wave resonances. We examine this phenomenon using the example of atmospheric planetary (Rossby) waves. It is shown that the dynamics of the clusters is determined by the types of connections between neighboring triads within a cluster; these correspond to substantially different scenarios of energy flux between different triads. All the possible cases of the energy cascade termination are classified. Free and forced chaotic dynamics in the clusters are investigated: due to the huge fluctuations of the energy exchange between resonant triads these two types of evolution have a lot in common. It is confirmed that finite-dimensional wave turbulence in finite wave systems is fundamentally different from kinetic wave turbulence in infinite systems; the latter is described by wave kinetic equations that account for interactions with overlapping quasi-resonances of finite amplitude waves. The present results are directly applicable to finite-dimensional wave turbulence in any wave system in finite domains with 3-mode interactions as encountered in hydrodynamics, astronomy, plasma physics, chemistry, medicine, etc.

preprint2009arXiv

Locality and Non-locality in Elasto-plastic Responses of Amorphous Solids

A number of current theories of plasticity in amorphous solids assume at their basis that plastic deformations are spatially localized. We present in this paper a series of numerical experiments to test the degree of locality of plastic deformation. These experiments increase in terms of the stringency of the removal of elastic contributions to the observed elasto-plastic deformations. It is concluded that for all our simulational protocols the plastic deformations are not localized, and their scaling is sub-extensive. We offer a number of measures of the magnitude of the plastic deformation, all of which display sub-extensive scaling characterized by non-trivial exponents. We provide some evidence that the scaling exponents governing the sub-extensive scaling laws are non-universal, depending on the degree of disorder and on the parameters of the systems. Nevertheless understanding what determines these exponents should shed considerable light on the physics of amorphous solids.

preprint2009arXiv

Predictive Statistical Mechanics for Glass Forming Systems

Using two extremely different models of glass formers in two and three dimensions we demonstrate how to encode the subtle changes in the geometric rearrangement of particles during the scenario of the glass transition. We construct a statistical mechanical description that is able to explain and predict the geometric rearrangement, the temperature dependent thermodynamic functions and the $α$-relaxation time within the measured temperature range and beyond. The theory is based on an up-scaling to proper variables (quasi-species) which is validated using a simple criterion. Once constructed, the theory provides an accurate predictive tool for quantities like the specific heat or the entropy at temperatures that cannot be reached by measurements. In addition, the theory identifies a rapidly increasing typical length scale $ξ$ as the temperature decreases. This growing spatial length scale determines the $α$-relaxation time as $τ_α\sim \exp(μξ/T)$ where $μ$ is a typical chemical potential per unit length.

preprint2009arXiv

Relations between the Material Mechanical Parameters and the Inter-particle Potential in Amorphous Solids

The shear-modulus and yield-stress of amorphous solids are important material parameters, with the former determining the rate of increase of stress under external strain and the latter being the stress value at which the material flows in a plastic manner. It is therefore important to understand how these parameters can be related to the inter-particle potential. Here a scaling theory is presented such that given the inter-particle potential, the dependence of the yield stress and the shear modulus on the density of the solid can be predicted in the athermal limit. It is explained when such prediction is possible at all densities and when it is only applicable at high densities. These results open up exciting possibilities for designing in principle new materials with desirable mechanical properties.

preprint2009arXiv

Scaling Theory for Steady State Plastic Flows in Amorphous Solids

Strongly correlated amorphous solids are a class of glass-formers whose inter-particle potential admits an approximate inverse power-law form in a relevant range of inter-particle distances. We study the steady-state plastic flow of such systems, firstly in the athermal, quasi-static limit, and secondly at finite temperatures and strain rates. In all cases we demonstrate the usefulness of scaling concepts to reduce the data to universal scaling functions where the scaling exponents are determined a-priori from the inter-particle potential. In particular we show that the steady plastic flow at finite temperatures with efficient heat extraction is uniquely characterized by two scaled variables; equivalently, the steady state displays an equation of state that relates one scaled variable to the other two. We discuss the range of applicability of the scaling theory, and the connection to density scaling in supercooled liquid dynamics. We explain that the description of transient states calls for additional state variables whose identity is still far from obvious.

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.

preprint2009arXiv

Stochastic Processes Crossing from Ballistic to Fractional Diffusion with Memory: Exact Results

We address the now classical problem of a diffusion process that crosses over from a ballistic behavior at short times to a fractional diffusion (sub- or super-diffusion) at longer times. Using the standard non-Markovian diffusion equation we demonstrate how to choose the memory kernel to exactly respect the two different asymptotics of the diffusion process. Having done so we solve for the probability distribution function (pdf) as a continuous function which evolves inside a ballistically expanding domain. This general solution agrees for long times with the pdf obtained within the continuous random walk approach but it is much superior to this solution at shorter times where the effect of the ballistic regime is crucial.

preprint2009arXiv

The Effective Temperature in Elasto-Plasticity of Amorphous Solids

An effective temperature $T_{\rm eff}$ which differs from the bath temperature is believed to play an essential role in the theory of elasto-plasticity of amorphous solids. The definition of a measurable $T_{\rm eff}$ in the literature on sheared solids suffers however from being connected to a fluctuation-dissipation theorem which is correct only in equilibrium. Here we introduce a natural definition of $T_{\rm eff}$ based on measurable structural features without recourse to any questionable assumption. The value of $T_{\rm eff}$ is connected, using theory and scaling concepts, to the flow stress and the mean energy that characterize the elasto-plastic flow.

preprint2007arXiv

Athermal Shear-Transformation-Zone Theory of Amorphous Plastic Deformation II: Analysis of Simulated Amorphous Silicon

In the preceding paper, we developed an athermal shear-transformation-zone (STZ) theory of amorphous plasticity. Here we use this theory in an analysis of numerical simulations of plasticity in amorphous silicon by Demkowicz and Argon (DA). In addition to bulk mechanical properties, those authors observed internal features of their deforming system that challenge our theory in important ways. We propose a quasithermodynamic interpretation of their observations in which the effective disorder temperature, generated by mechanical deformation well below the glass temperature, governs the behavior of other state variables that fall in and out of equilibrium with it. Our analysis points to a limitation of either the step-strain procedure used by DA in their simulations, or the STZ theory in its ability to describe rapid transients in stress-strain curves, or perhaps to both. Once we allow for this limitation, we are able to bring our theoretical predictions into accurate agreement with the simulations.

preprint2006arXiv

Athermal Shear-Transformation-Zone Theory of Amorphous Plastic Deformation I: Basic Principles

We develop an athermal version of the shear-transformation-zone (STZ) theory of amorphous plasticity in materials where thermal activation of irreversible molecular rearrangements is negligible or nonexistent. In many respects, this theory has broader applicability and yet is simpler than its thermal predecessors. For example, it needs no special effort to assure consistency with the laws of thermodynamics, and the interpretation of yielding as an exchange of dynamic stability between jammed and flowing states is clearer than before. The athermal theory presented here incorporates an explicit distribution of STZ transition thresholds. Although this theory contains no conventional thermal fluctuations, the concept of an effective temperature is essential for understanding how the STZ density is related to the state of disorder of the system.

preprint2006arXiv

Direct Identification of the Glass Transition: Growing Length Scale and the Onset of Plasticity

Understanding the mechanical properties of glasses remains elusive since the glass transition itself is not fully understood, even in well studied examples of glass formers in two dimensions. In this context we demonstrate here: (i) a direct evidence for a diverging length scale at the glass transition (ii) an identification of the glass transition with the disappearance of fluid-like regions and (iii) the appearance in the glass state of fluid-like regions when mechanical strain is applied. These fluid-like regions are associated with the onset of plasticity in the amorphous solid. The relaxation times which diverge upon the approach to the glass transition are related quantitatively.

preprint2006arXiv

Oscillatory Instability in Two-Dimensional Dynamic Fracture

The stability of a rapid dynamic crack in a two dimensional infinite strip is studied in the framework of Linear Elasticity Fracture Mechanics supplemented with a modified principle of local symmetry. It is predicted that a single crack becomes unstable by a finite wavelength oscillatory mode at a velocity $v_c$, $0.8c_R<v_c<0.85c_R$, where $c_R$ is the Rayleigh wave speed. The relevance of this theoretical calculation to the oscillatory instability reported in the companion experimental Letter is discussed.

preprint2005arXiv

Fracture Surfaces as Multiscaling Graphs

Fracture paths in quasi-two-dimenisonal (2D) media (e.g thin layers of materials, paper) are analyzed as self-affine graphs $h(x)$ of height $h$ as a function of length $x$. We show that these are multiscaling, in the sense that $n^{th}$ order moments of the height fluctuations across any distance $\ell$ scale with a characteristic exponent that depends nonlinearly on the order of the moment. Having demonstrated this, one rules out a widely held conjecture that fracture in 2D belongs to the universality class of directed polymers in random media. In fact, 2D fracture does not belong to any of the known kinetic roughening models. The presence of multiscaling offers a stringent test for any theoretical model; we show that a recently introduced model of quasi-static fracture passes this test.

preprint2004arXiv

Saturation of Turbulent Drag Reduction in Dilute Polymer Solutions

Drag reduction by polymers in turbulent wall-bounded flows exhibits universal and non-universal aspects. The universal maximal mean velocity profile was explained in a recent theory. The saturation of this profile and the crossover back to the Newtonian plug are non-universal, depending on Reynolds number Re, concentration of polymer $c_p$ and the degree of polymerization $N_p$. We explain the mechanism of saturation stemming from the finiteness of extensibility of the polymers, predict its dependence on $c_p$ and $N$ in the limit of small $c_p$ and large Re, and present the excellent comparison of our predictions to experiments on drag reduction by DNA.

preprint2001arXiv

Retardation of the onset of yurbulence by minor viscosity contrasts

Motivated by the large effect of turbulent drag reduction by minute concentrations of polymers we study the effects of minor viscosity contrasts on the stability of hydrodynamic flows. The key player is a localized region where the energy of fluctuations is produced by interactions with the mean flow (the "critical layer"). We show that a layer of weakly space-dependent viscosity placed near the critical layer can have very large stabilizing effect on hydrodynamic fluctuations, retarding significantly the onset of turbulence. The effect is not due to a modified dissipation (as is assumed in theories of drag reduction), but due to reduced energy intake from the mean flow to the fluctuations. We propose that similar physics act in turbulent drag reduction.

preprint1998arXiv

Continued Fraction Representation of Temporal Multi Scaling in Turbulence

It was shown recently that the anomalous scaling of simultaneous correlation functions in turbulence is intimately related to the breaking of temporal scale invariance, which is equivalent to the appearance of infinitely many times scales in the time dependence of time-correlation functions. In this paper we derive a continued fraction representation of turbulent time correlation functions which is exact and in which the multiplicity of time scales is explicit. We demonstrate that this form yields precisely the same scaling laws for time derivatives and time integrals as the "multi-fractal" representation that was used before. Truncating the continued fraction representation yields the "best" estimates of time correlation functions if the given information is limited to the scaling exponents of the simultaneous correlation functions up to a certain, finite order. It is worth noting that the derivation of a continued fraction representation obtained here for an operator which is not Hermitian or anti-Hermitian may be of independent interest.

preprint1997arXiv

Computing the Scaling Exponents in Fluid Turbulence from First Principles: Demonstration of Multi-scaling

This manuscript is a draft of work in progress, meant for network distribution only. It will be updated to a formal preprint when the numerical calculations will be accomplished. In this draft we develop a consistent closure procedure for the calculation of the scaling exponents $ζ_n$ of the $n$th order correlation functions in fully developed hydrodynamic turbulence, starting from first principles. The closure procedure is constructed to respect the fundamental rescaling symmetry of the Euler equation. The starting point of the procedure is an infinite hierarchy of coupled equations that are obeyed identically with respect to scaling for any set of scaling exponents $ζ_n$. This hierarchy was discussed in detail in a recent publication [V.S. L'vov and I. Procaccia, Phys. Rev. E, submitted, chao-dyn/9707015]. The scaling exponents in this set of equations cannot be found from power counting. In this draft we discuss in detail low order non-trivial closures of this infinite set of equations, and prove that these closures lead to the determination of the scaling exponents from solvability conditions. The equations under consideration after this closure are nonlinear integro-differential equations, reflecting the nonlinearity of the original Navier-Stokes equations. Nevertheless they have a very special structure such that the determination of the scaling exponents requires a procedure that is very similar to the solution of linear homogeneous equations, in which amplitudes are determined by fitting to the boundary conditions in the space of scales. The re-normalization scale that i necessary for any anomalous scaling appears at this point. The Holder inequalities on the scaling exponents select the renormalizaiton scale as the outer scale of turbulence $L$.

preprint1997arXiv

Computing the Scaling Exponents in Fluid Turbulence from First Principles: the Formal Setup

We propose a scheme for the calculation from the NS equations of the scaling exponents $ζ_n$ of the $n$th order correlators in fully developed hydrodynamic turbulence. The scheme is nonperturbative and constructed to respect the fundamental rescaling symmetry of the Euler equation. It constitutes an infinite hierarchy of coupled equations that are obeyed identically with respect to scaling for any set of $ζ_n$. As a consequence $ζ_n$ are determined by solvability conditions and not from power counting. It is arguedthat in order to achieve such a formulation one must recognize that the many-point space-time correlators are not scale invariant in time arguments. The assumption of full scale invariance leads unavoidably to K41 exponents. It is argued the determination of of $ζ_n$ requires equations for infinitely many renormalized objects. One ca n however proceed in controlled successive approximations by truncations of the infinite hierarchy of equations. Clues as to how to truncate without reintroducing power counting can be obtained from renormalized perturbation theory. To this aim we show that the fully resummed perturbation theory is equivalent in its contents to the exact hierarchy of equations obeyed by the $n$th order correlators and Green's function. In light of this important result we can safely use finite resummations to construct successive closures. This paper presents the conceptual and technical details of the scheme. The analysis of the high-order closure procedures which do not destroy the rescaling symmetry and the actual calculations for truncated models will be presented in a forthcoming paper in collaboration with V. Belinicher.

preprint1996arXiv

Hydrodynamic Turbulence Has Infinitely Many Anomalous Dynamical Exponents

On the basis of the Navier-Stokes equations we develop the statistical theory of many space-time correlation functions of velocity differences. Their time dependence is {\em not} scale invariant: $n$-order correlations functions exhibit $n-1$ distinct decorrelation times that are characterized by $n-1$ anomalous dynamical scaling exponents. We derive exact scaling relations that bridge all these dynamical exponents to the static anomalous exponents $ζ_n$ of the standard structure functions.

preprint1996arXiv

Non-Perturbative Zero Modes in the Kraichnan Model for Turbulent Advection

The anomalous scaling behavior of the n-th order correlation functions ${\cal F}_n$ of the Kraichnan model of turbulent passive scalar advection is believed to be dominated by the homogeneous solutions (zero-modes) of the Kraichnan equation $\hat{\cal B}_n {\cal F}_n=0$. Previous analysis found zero-modes in perturbation theory in a small parameter. We present non-perturbative analysis of the simplest (non-trivial) case of n=3 and compare the results with the perturbative predictions.

preprint1995arXiv

Exact Resummations in the Theory of Hydrodynamic Turbulence: 0. Line-Resummed Diagrammatic Perturbation Approach

The lectures presented by one of us (IP) at the Les Houches summer school dealt with the scaling properties of high Reynolds number turbulence in fluid flows. The results presented are available in the literature and there is no real need to reproduce them here. Quite on the contrary, some of the basic tools of the field and theoretical techniques are not available in a pedagogical format, and it seems worthwhile to present them here for the benefit of the interested student. We begin with a detailed exposition of the naive perturbation theory for the ensemble averages of hydrodynamic observables (the mean velocity, the response functions and the correlation functions). The effective expansion parameter in such a theory is the Reynolds number (Re); one needs therefore to perform infinite resummations to change the effective expansion parameter. We present in detail the Dyson-Wyld line resummation which allows one to dress the propagators, and to change the effective expansion parameter from Re to O(1). Next we develop the ``dressed vertex" representation of the diagrammatic series. Lastly we discuss in full detail the path-integral formulation of the statistical theory of turbulence, and show that it is equivalent order by order to the Dyson-Wyld theory. On the basis of the material presented here one can proceed smoothly to read the recent developments in this field.