Search

The equipartition theorem states that in equilibrium, thermal energy is equally distributed among uncoupled degrees of freedom that appear quadratically in the system's Hamiltonian. However, for spatially coupled degrees of freedom, such as interacting particles, one may speculate that the spatial distribution of thermal energy may differ from the value predicted by equipartition, possibly quite substantially in strongly inhomogeneous or disordered systems. Here we show that for systems undergoing simple Gaussian fluctuations around an equilibrium state, the spatial distribution is universally bounded from above by 1/2k(B)T. We further show that in one-dimensional systems with short-range interactions, the thermal energy is equally partitioned even for coupled degrees of freedom in the thermodynamic limit and that in higher dimensions nontrivial spatial distributions emerge. Some implications are discussed.

Variable-amplitude oscillatory shear tests are emerging as powerful tools to investigate and quantify the nonlinear rheology of amorphous solids, complex fluids, and biological materials. Quite a few recent experimental and atomistic simulation studies demonstrated that at low shear amplitudes, an amorphous solid settles into an amplitude- and initial-conditions-dependent dissipative limit cycle, in which back-and-forth localized particle rearrangements periodically bring the system to the same state. At sufficiently large shear amplitudes, the amorphous system loses memory of the initial conditions, exhibits chaotic particle motions accompanied by diffusive behavior, and settles into a stochastic steady state. The two regimes are separated by a transition amplitude, possibly characterized by some critical-like features. Here we argue that these observations support some of the physical assumptions embodied in the nonequilibrium thermodynamic, internal-variables based, shear-transformation-zone model of amorphous viscoplasticity; most notably that "flow defects" in amorphous solids are characterized by internal states between which they can make transitions, and that structural evolution is driven by dissipation associated with plastic deformation. We present a rather extensive theoretical analysis of the thermodynamic shear-transformation-zone model for a variable-amplitude oscillatory shear protocol, highlighting its success in accounting for various experimental and simulational observations, as well as its limitations. Our results offer a continuum-level theoretical framework for interpreting the variable-amplitude oscillatory shear response of amorphous solids and may promote additional developments.

The strength and stability of frictional interfaces, ranging from tribological systems to earthquake faults, are intimately related to the underlying spatially extended dynamics. Here we provide a comprehensive theoretical account, both analytic and numeric, of spatiotemporal interfacial dynamics in a realistic rate-and-state friction model, featuring both velocity-weakening and velocity-strengthening behaviors. Slowly extending, loading-rate-dependent creep patches undergo a linear instability at a critical nucleation size, which is nearly independent of interfacial history, initial stress conditions, and velocity-strengthening friction. Nonlinear propagating rupture fronts-the outcome of instability-depend sensitively on the stress state and velocity-strengthening friction. Rupture fronts span a wide range of propagation velocities and are related to steady-state-front solutions.

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.

A wide range of glassy and disordered materials exhibit complex, nonexponential, structural relaxation (aging). We propose a simple nonlinear rate equation δ = a[1-exp(b δ)], where δ is the normalized deviation of a macroscopic variable from its equilibrium value, to describe glassy relaxation. Analysis of extensive experimental data shows that this equation quantitatively captures structural relaxation, where a and b are both temperature- and, more importantly, history-dependent parameters. This analysis explicitly demonstrates that structural relaxation cannot be accurately described by a single nonequilibrium variable. Relaxation rates extracted from the data imply the existence of cooperative rearrangements on a supermolecular scale.

We present a linearized shear-transformation-zone (STZ) theory of glassy dynamics in which the internal STZ transition rates are characterized by a broad distribution of activation barriers. For slowly aging or fully aged systems, the main features of the barrier-height distribution are determined by the effective temperature and other near-equilibrium properties of the configurational degrees of freedom. Our theory accounts for the wide range of relaxation rates observed in both metallic glasses and soft glassy materials such as colloidal suspensions. We find that the frequency-dependent loss modulus is not just a superposition of Maxwell modes. Rather, it exhibits an α peak that rises near the viscous relaxation rate and, for nearly jammed, glassy systems, extends to much higher frequencies in accord with experimental observations. We also use this theory to compute strain recovery following a period of large, persistent deformation and then abrupt unloading. We find that strain recovery is determined in part by the initial barrier-height distribution, but that true structural aging also occurs during this process and determines the system's response to subsequent perturbations. In particular, we find by comparison with experimental data that the initial deformation produces a highly disordered state with a large population of low activation barriers, and that this state relaxes quickly toward one in which the distribution is dominated by the high barriers predicted by the near-equilibrium analysis. The nonequilibrium dynamics of the barrier-height distribution is the most important of the issues raised and left unresolved in this paper.

The recently developed weakly nonlinear theory of dynamic fracture predicts 1/r corrections to the standard asymptotic linear elastic 1/(square root of r) displacement gradients, where r is measured from the tip of a tensile crack. We show that the 1/r singularity does not automatically conform with the notion of autonomy (autonomy means that any crack tip nonlinear solution is uniquely determined by the surrounding linear elastic 1/(square root of r) fields) and that it does not automatically satisfy the resultant Newton's equation in the crack parallel direction. We show that these two properties are interrelated and that by requiring that the resultant Newton's equation is satisfied, autonomy of the 1/r singular solution is retained. We further show that the resultant linear momentum carried by the 1/r singular fields vanishes identically. Our results, which reveal the physical and mathematical nature of the solution, are in favorable agreement with recent near-tip measurements.

We use the internal-variable, effective-temperature thermodynamics developed in two preceding papers to reformulate the shear-transformation-zone (STZ) theory of amorphous plasticity. As required by the preceding analysis, we make explicit approximations for the energy and entropy of the STZ internal degrees of freedom. We then show that the second law of thermodynamics constrains the STZ transition rates to have an Eyring form as a function of the effective temperature. Finally, we derive an equation of motion for the effective temperature for the case of STZ dynamics.

We develop a theory of the effective disorder temperature in glass-forming materials driven away from thermodynamic equilibrium by external forces. Our basic premise is that the slow configurational degrees of freedom of such materials are weakly coupled to the fast kinetic-vibrational degrees of freedom and therefore that these two subsystems can be described by different temperatures during deformation. We use results from the preceding paper on the nonequilibrium thermodynamics of systems with internal degrees of freedom to derive an equation of motion for the effective temperature and to learn how this temperature couples to the dynamics of the system as a whole.

This is the first of three papers devoted to the nonequilibrium thermodynamics of amorphous materials. Our focus here is on the role of internal degrees of freedom in determining the dynamics of such systems. For illustrative purposes, we study a solid whose internal degrees of freedom are vacancies that govern irreversible volume changes. Using this model, we compare a thermodynamic theory based on the Clausius-Duhem inequality to a statistical analysis based directly on the law of increase of entropy. The statistical theory is used first to derive the Clausius-Duhem inequality. We then use the theory to go beyond those results and obtain detailed equations of motion, including a rate factor that is enhanced by deformation-induced noisy fluctuations. The statistical analysis points to the need for understanding how both energy and entropy are shared by the vacancies and their environments.

The dynamics of rapid brittle cracks is commonly studied in the framework of linear elastic fracture mechanics where nonlinearities are neglected. However, recent experimental and theoretical work demonstrated explicitly the importance of elastic nonlinearities in fracture dynamics. We study two simple one-dimensional models of fracture in order to gain insights about the role of elastic nonlinearities and the implications of their exclusion in the common linear elastic approximation. In one model we consider the decohesion of a nonlinear elastic membrane from a substrate. In a second model we follow the philosophy of linear elastic fracture mechanics and study a linearized version of the nonlinear model. By analyzing the steady state solutions, the velocity-load relations and the response to perturbations of the two models we show that the linear approximation fails at finite crack tip velocities. We highlight certain features of the breakdown of the linear theory and discuss possible implications of our results to higher dimensional systems.

The existence of a finite threshold, the yield stress, for the onset of plastic yielding is a universal feature of plasticity. This jamming-unjamming transition is naturally accounted for by the dynamics of a bistable internal state field. We show, within the athermal shear transformation zone theory of amorphous plasticity, that the transition is accompanied by the propagation of plastic fronts. We further show that the mean-field theory cannot select the velocity of these fronts, and go beyond the mean-field description to include fluctuations and correlations effects, resulting in additional nonlocal terms in the equations. Finally, we demonstrate that these terms, with an associated intrinsic length scale, provide a velocity selection mechanism for the plastic fronts.

Recently, the existence and properties of unbounded cavity modes, resulting in extensive plastic deformation failure of two-dimensional sheets of amorphous media, were discussed in the context of the athermal shear-transformation-zones (STZ) theory. These modes pertain to perfect circular symmetry of the cavity and the stress conditions. In this paper we study the shape stability of the expanding circular cavity against perturbations, in both the unbounded and the bounded growth regimes (for the latter the unperturbed theory predicts no catastrophic failure). Since the unperturbed reference state is time dependent, the linear stability theory cannot be cast into standard time-independent eigenvalue analysis. The main results of our study are as follows: (i) sufficiently small perturbations are stable; (ii) larger perturbations within the formal linear decomposition may lead to an instability; this dependence on the magnitude of the perturbations in the linear analysis is a result of the nonstationarity of the growth; and (iii) the stability of the circular cavity is particularly sensitive to perturbations in the effective disorder temperature; in this context we highlight the role of the rate sensitivity of the limiting value of this effective temperature. Finally we point to the consequences of the form of the stress dependence of the rate of STZ transitions. The present analysis indicates the importance of nonlinear effects that were not taken into account yet. Furthermore, the analysis suggests that details of the constitutive relations appearing in the theory can be constrained by the modes of macroscopic failure in these amorphous systems.

Recent developments in the theory of amorphous plasticity point to the central role played by the concept of an effective disorder temperature T_{eff} . An athermal dynamics for T_{eff} is proposed in the framework of a deformation theory and discussed in light of recent steady-state simulations by Haxton and Liu [T. K. Haxton and A. J. Liu, Phys. Rev. Lett. 99, 195701 (2007)]. The structure of the resulting theory, its parameters and transient dynamics are discussed and compared to available data.

The understanding of dynamic failure in amorphous materials via the propagation of free boundaries like cracks and voids must go beyond elasticity theory, since plasticity intervenes in a crucial and poorly understood manner near the moving free boundary. We focus on failure via a cavitation instability in a radially symmetric stressed material and set up the free boundary dynamics taking both elasticity and viscoplasticity into account using the recently proposed athermal shear transformation zone theory. We demonstrate that this theory predicts the existence (in amorphous systems) of fast cavitation modes accompanied by extensive plastic deformations and discuss the revealed physics.

We develop an athermal shear-transformation-zone (STZ) theory of plastic deformation in spatially inhomogeneous, amorphous solids. Our ultimate goal is to describe the dynamics of the boundaries of voids or cracks in such systems when they are subjected to remote, time-dependent tractions. The theory is illustrated here for the case of a circular hole in an infinite two-dimensional plate, a highly symmetric situation that allows us to solve much of the problem analytically. In spite of its special symmetry, this example contains many general features of systems in which stress is concentrated near free boundaries and deforms them irreversibly. We depart from conventional treatments of such problems in two ways. First, the STZ analysis allows us to keep track of spatially heterogeneous, internal state variables such as the effective disorder temperature, which determines plastic response to subsequent loading. Second, we subject the system to stress pulses of finite duration, and therefore are able to observe elastoplastic response during both loading and unloading. We compute the final deformations and residual stresses produced by these stress pulses. Looking toward more general applications of these results, we examine the possibility of constructing a boundary-layer theory that might be useful in less symmetric situations.

The long-range elastic model, which is believed to describe the evolution of a self-affine rough crack front, is analyzed to linear and nonlinear orders. It is shown that the nonlinear terms, while important in changing the front dynamics, do not change the scaling exponent which characterizes the roughness of the front. The scaling exponent thus predicted by the model is much smaller than the one observed experimentally. The inevitable conclusion is that the gap between the results of experiments and the model that is supposed to describe them is too large and some new physics has to be invoked for another model.

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.

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.

We address the role of the nature of material disorder in determining the roughness of cracks, which grow by damage nucleation and coalescence ahead of the crack tip. We highlight the role of quenched and annealed disorders in relation to the length scales d and xic associated with the disorder and the damage nucleation, respectively. In two related models, one with quenched disorder in which d approximately xic, the other with annealed disorder in which d<

The apparent similarity of microbranching instabilities in different brittle materials gave rise to a widely held belief that many aspects of the postinstability physics were universal. We propose that the physics determining the typical length and time scales characterizing the postinstability patterns differ greatly from material to material. We offer a scaling theory connecting the pattern characteristics to material properties (like molecular weight) in brittle plastics like PMMA, and stress the fundamental differences with patterns in glass which are crucially influenced by three-dimensional dynamics. In both cases the present ab initio theoretical models are still too far from reality, disregarding some fundamental physics of the phenomena.

Slow crack propagation in ductile, and in certain brittle materials, appears to take place via the nucleation of voids ahead of the crack tip due to plastic yields, followed by the coalescence of these voids. Postmortem analysis of the resulting fracture surfaces of ductile and brittle materials on the microm-mm and the nm scales, respectively, reveals self-affine cracks with anomalous scaling exponent zeta approximately = 0.8 in 3 dimensions and zeta approximately = 0.65 in 2 dimensions. In this paper we present an analytic theory based on the method of iterated conformal maps aimed at modelling the void formation and the fracture growth, culminating in estimates of the roughening exponents in 2 dimensions. In the simplest realization of the model we allow one void ahead of the crack, and address the robustness of the roughening exponent. Next we develop the theory further, to include two voids ahead of the crack. This development necessitates generalizing the method of iterated conformal maps to include doubly connected regions (maps from the annulus rather than the unit circle). While mathematically and numerically feasible, we find that the employment of the stress field as computed from elasticity theory becomes questionable when more than one void is explicitly inserted into the material. Thus further progress in this line of research calls for improved treatment of the plastic dynamics.

We propose a theoretical model for branching instabilities in 2-dimensional fracture, offering predictions for when crack branching occurs, how multiple cracks develop, and what is the geometry of multiple branches. The model is based on equations of motion for crack tips which depend only on the time dependent stress intensity factors. The latter are obtained by invoking an approximate relation between static and dynamic stress intensity factors, together with an essentially exact calculation of the static ones. The results of this model are in qualitative agreement with a number of experiments in the literature.

We address the interaction of fast moving cracks in stressed materials with microcracks on their way, considering it as one possible mechanism for fluctuations in the velocity of the main crack (irrespective whether the microcracks are existing material defects or they form during the crack evolution). We analyze carefully the dynamics (in two space dimensions) of one macrocrack and one microcrack, and demonstrate that their interaction results in a large and rapid velocity fluctuation, in qualitative correspondence with typical velocity fluctuations observed in experiments. In developing the theory of the dynamical interaction we invoke an approximation that affords a reduction in mathematical complexity to a simple set of ordinary differential equations for the positions of the crack tips; we propose that this kind of approximation has a range of usefulness that exceeds the present context.

The calculation of the stress field around an arbitrarily shaped crack in an infinite two-dimensional elastic medium is a mathematically daunting problem. With the exception of few exactly soluble crack shapes the available results are based on either perturbative approaches or on combinations of analytic and numerical techniques. We present here a general solution of this problem for any arbitrary crack. Along the way we develop a method to compute the conformal map from the exterior of a circle to the exterior of a line of arbitrary shape, offering it as a superior alternative to the classical Schwartz-Cristoffel transformation. Our calculation results in an accurate estimate of the full stress field and in particular of the stress intensity factors K(I) and K(II) and the T-stress which are essential in the theory of fracture.

We address the theory of quasistatic crack propagation in a strip of glass that is pulled from a hot oven towards a cold bath. This problem had been carefully studied in a number of experiments that offer a wealth of data to challenge the theory. We improve upon previous theoretical treatments in a number of ways. First, we offer a technical improvement of the discussion of the instability towards the creation of a straight crack. This improvement consists in employing Padé approximants to solve the relevant Wiener-Hopf factorization problem that is associated with this transition. Next we improve the discussion of the onset of oscillatory instability towards an undulating crack. We offer a way of considering the problem as a sum of solutions of a finite strip without a crack and an infinite medium with a crack. This allows us to present a closed form solution of the stress intensity factors in the vicinity of the oscillatory instability. Most importantly we develop a dynamical description of the actual trajectory in the regime of oscillatory crack. This theory is based on the dynamical law for crack propagation proposed by Hodgdon and Sethna. We show that this dynamical law results in a solution of the actual crack trajectory in post-critical conditions; we can compute from first principles the critical value of the control parameters, and the characteristics of the solution such as the wavelength of the oscillations. We present detailed comparison with experimental measurements without any free parameters. The comparison appears quite excellent. Finally we show that the dynamical law can be translated to an equation for the amplitude of the oscillatory crack; this equation predicts correctly the scaling exponents observed in experiments.