Driving-induced crossover: from classical criticality to self-organized criticality

Thermal effects on fracture and the brittle-to-ductile transition

The fracture behavior of brittle and ductile materials can be strongly influenced by thermal fluctuations, especially in micro- and nanodevices as well as in rubberlike and biological materials. However, temperature effects, in particular on the brittle-to-ductile transition, still require a deeper theoretical investigation. As a step in this direction we propose a theory, based on equilibrium statistical mechanics, able to describe the temperature-dependent brittle fracture and brittle-to-ductile transition in prototypical discrete systems consisting in a lattice with breakable elements. Concerning the brittle behavior, we obtain closed form expressions for the temperature-dependent fracture stress and strain, representing a generalized Griffith criterion, ultimately describing the fracture as a genuine phase transition. With regard to the brittle-to-ductile transition, we obtain a complex critical scenario characterized by a threshold temperature between the two fracture regimes (brittle and ductile), an upper and a lower yield strength, and a critical temperature corresponding to the complete breakdown. To show the effectiveness of the proposed models in describing thermal fracture behaviors at small scales, we successfully compare our theoretical results with molecular dynamics simulations of Si and GaN nanowires.

Dislocation avalanches from strain-controlled loading: A discrete dislocation dynamics study

We study strain-controlled plastic deformation of crystalline solids via two-dimensional discrete dislocation dynamics simulations. To this end, we characterize the average stress-strain curves as well as the statistical properties of strain bursts and the related stress drops as a function of the imposed strain rate and the stiffness of the specimen-machine system. The dislocation system exhibits strain-rate sensitivity such that a larger imposed strain rate results in a higher average stress at a given strain. In the limit of small strain rate and driving spring stiffness, the sizes and durations of the dislocation avalanches are power law distributed up to a cutoff scale, and exhibit temporally asymmetric average shapes. We discuss the dependence of the results on the driving parameters and compare our results to those from previous simulations where quasistatic stress-controlled loading was used.

Non-trivial avalanches triggered by shear banding in compression of metallic glass foams

Ductile metallic glass foams (DMGFs) are a new type of structural material with a perfect combination of high strength and toughness. Owing to their disordered atomic-scale microstructures and randomly distributed macroscopic voids, the compressive deformation of DMGFs proceeds through multiple nanoscale shear bands accompanied by local fracture of cellular structures, which induces avalanche-like intermittences in stress–strain curves. In this paper, we present a statistical analysis, including distributions of avalanche size, energy dissipation, waiting times and aftershock sequence, on such a complex dynamic process, which is dominated by shear banding. After eliminating the influence of structural disorder, we demonstrate that, in contrast to the mean-field results of their brittle counterparts, scaling laws in DMGFs are characterized by different exponents. It is shown that the occurrence of non-trivial scaling behaviours is attributed to the localized plastic yielding, which effectively prevents the system from building up a long-range correlation. This accounts for the high structural stability and energy absorption performance of DMGFs. Furthermore, our results suggest that such shear banding dynamics introduce an additional characteristic time scale, which leads to a universal gamma distribution of waiting times.

Rigidity-Controlled Crossover: From Spinodal to Critical Failure

Failure in disordered solids is accompanied by intermittent fluctuations extending over a broad range of scales. The implied scaling has been previously associated with either spinodal or critical points. We use an analytically transparent mean-field model to show that both analogies are relevant near the brittle-to-ductile transition. Our study indicates that in addition to the strength of quenched disorder, an appropriately chosen global measure of rigidity (connectivity) can be also used to tune the system to criticality. By interpreting rigidity as a timelike variable we reveal an intriguing parallel between earthquake-type critical failure and Burgers turbulence.

Critical disorder and critical magnetic field of the nonequilibrium athermal random-field Ising model in thin systems

In the present study of the nonequilibrium athermal random-field Ising model we focus on the behavior of the critical disorder R_{c}(l) and the critical magnetic field H_{c}(l) under different boundary conditions when the system thickness l varies. We propose expressions for R_{c}(l) and H_{c}(l) as well as for the effective critical disorder R_{c}^{eff}(l,L) and effective critical magnetic field H_{c}^{eff}(l,L) playing the role of the effective critical parameters for the L×L×l lattices of finite lateral size L. We support these expressions by the scaling collapses of the magnetization and susceptibility curves obtained in extensive simulations. The collapses are achieved with the two-dimensional (2D) exponents for l below some characteristic value, providing thus a numerical evidence that the thin systems exhibit a 2D-like criticality which should be relevant for the experimental analyses of thin ferromagnetic samples.

Acoustic Emissions during Structural Changes in Shape Memory Alloys

Structural changes (martensitic transformation, rearrangements of martensitic variants) in shape memory alloys have an intermittent character that is accompanied by the emission of different (thermal, acoustic, and magnetic) noises, which are fingerprints of the driven criticality, resulting in a damped power-law behaviour. We will illustrate what kinds of important information can be obtained on the structural changes in shape memory alloys. It was established that the power exponents of distributions of acoustic emission (AE) parameters (energy, amplitude, etc.), belonging to martensitic transformations, show quite a universal character and depend only on the symmetry of the martensite. However, we have shown that the asymmetry of the transformation (the exponents are different for the forward and reverse transformations) results in as large differences as those due to the martensite symmetry. We will also demonstrate how the recently introduced AE clustering method can help to identify the different contributions responsible for the asymmetry. The usefulness of the investigations of time correlations between the subsequent events and correlations between acoustic and magnetic noise events in ferromagnetic shape memory alloys will be demonstrated too. Finally, examples of acoustic and magnetic emissions during variant rearrangements (superplastic or superelastic behaviour) in the martensitic state will be described.

Fast generation of ultrastable computer glasses by minimization of an augmented potential energy

We present a model and protocol that enable the generation of extremely stable computer glasses at minimal computational cost. The protocol consists of an instantaneous quench in an augmented potential energy landscape, with particle radii as additional degrees of freedom. We demonstrate how our glasses' mechanical stability, which is readily tunable in our approach, is reflected in both microscopic and macroscopic observables. Our observations indicate that the stability of our computer glasses is at least comparable to that of computer glasses generated by the celebrated Swap Monte Carlo algorithm. Strikingly, some key properties support even qualitatively enhanced stability in our scheme: the density of quasilocalized excitations displays a gap in our most stable computer glasses, whose magnitude scales with the polydispersity of the particles. We explain this observation, which is consistent with the lack of plasticity we observe at small stress. It also suggests that these glasses are depleted from two-level systems, similarly to experimental vapor-deposited ultrastable glasses.

Random critical point separates brittle and ductile yielding transitions in amorphous materials

We combine an analytically solvable mean-field elasto-plastic model with molecular dynamics simulations of a generic glass former to demonstrate that, depending on their preparation protocol, amorphous materials can yield in two qualitatively distinct ways. We show that well-annealed systems yield in a discontinuous brittle way, as metallic and molecular glasses do. Yielding corresponds in this case to a first-order nonequilibrium phase transition. As the degree of annealing decreases, the first-order character becomes weaker and the transition terminates in a second-order critical point in the universality class of an Ising model in a random field. For even more poorly annealed systems, yielding becomes a smooth crossover, representative of the ductile rheological behavior generically observed in foams, emulsions, and colloidal glasses. Our results show that the variety of yielding behaviors found in amorphous materials does not necessarily result from the diversity of particle interactions or microscopic dynamics but is instead unified by carefully considering the role of the initial stability of the system.

Bifurcation analysis and phase diagram of a spin-string model with buckled states

We analyze a one-dimensional spin-string model, in which string oscillators are linearly coupled to their two nearest neighbors and to Ising spins representing internal degrees of freedom. String-spin coupling induces a long-range ferromagnetic interaction among spins that competes with a spin-spin antiferromagnetic coupling. As a consequence, the complex phase diagram of the system exhibits different flat rippled and buckled states, with first or second order transition lines between states. This complexity translates to the two-dimensional version of the model, whose numerical solution has been recently used to explain qualitatively the rippled to buckled transition observed in scanning tunneling microscopy experiments with suspended graphene sheets. Here we describe in detail the phase diagram of the simpler one-dimensional model and phase stability using bifurcation theory. This gives additional insight into the physical mechanisms underlying the different phases and the behavior observed in experiments.

Avalanche criticality in thermal-driven martensitic transitions: the asymmetry of the forward and reverse transitions in shape-memory materials

Martensitic transitions take place intermittently as a sequence of avalanches which are accompanied by the emission of acoustic waves. The study of this acoustic emission (AE) reveals the scale-free nature of the avalanches. In a number of shape memory materials undergoing a martensitic transition it has been found that, in spite of relatively low hysteresis, the dynamics of forward and reverse transitions are different, which may explain the fact that the AE activity is different in both forward and reverse transitions. The asymmetry could be a consequence of the fact that, while nucleation is required for the transition from the parent to martensitic phase to take place, reverse transition occurs by fast shrinkage of martensitic domains. We have analysed in detail the distribution of avalanches in cooling and heating runs in Fe-Pd and Cu-Zn-Al shape-memory alloys. In the former, the martensitic transition is weakly first order while it shows a significant first order character in the latter. We have found that in Fe-Pd the distributions are power law for the forward and reverse transitions characterized by the same critical exponents. For Cu-Zn-Al the distribution of avalanches is critical in forward transitions but exponentially damped in the reverse transition. It is suggested that this different behaviour could originate from the different dynamic mechanisms in forward and reverse transitions.

Avalanches and force drops in displacement-driven compression of porous glasses

Similarities between force-driven compression experiments of porous materials and earthquakes have been recently proposed. In this paper, we measure the acoustic emission during displacement-driven compression of a porous glass. The energy of acoustic-emission events shows that the failure process exhibits avalanche scale-invariance and therefore follows the Gutenberg-Richter law. The resulting exponents do not exhibit significant differences with respect the force-driven case. Furthermore, the force exhibits an avalanche-type behavior for which the force drops are power-law distributed and correlated with the acoustic emission events.

Layer-layer competition in multiplex complex networks

The coexistence of multiple types of interactions within social, technological and biological networks has moved the focus of the physics of complex systems towards a multiplex description of the interactions between their constituents. This novel approach has unveiled that the multiplex nature of complex systems has strong influence in the emergence of collective states and their critical properties. Here we address an important issue that is intrinsic to the coexistence of multiple means of interactions within a network: their competition. To this aim, we study a two-layer multiplex in which the activity of users can be localized in each of the layers or shared between them, favouring that neighbouring nodes within a layer focus their activity on the same layer. This framework mimics the coexistence and competition of multiple communication channels, in a way that the prevalence of a particular communication platform emerges as a result of the localization of user activity in one single interaction layer. Our results indicate that there is a transition from localization (use of a preferred layer) to delocalization (combined usage of both layers) and that the prevalence of a particular layer (in the localized state) depends on the structural properties.

Effects of local and global network connectivity on synergistic epidemics

Epidemics in networks can be affected by cooperation in transmission of infection and also connectivity between nodes. An interplay between these two properties and their influence on epidemic spread are addressed in the paper. A particular type of cooperative effects (called synergy effects) is considered, where the transmission rate between a pair of nodes depends on the number of infected neighbors. The connectivity effects are studied by constructing networks of different topology, starting with lattices with only local connectivity and then with networks that have both local and global connectivity obtained by random bond-rewiring to nodes within a certain distance. The susceptible-infected-removed epidemics were found to exhibit several interesting effects: (i) for epidemics with strong constructive synergy spreading in networks with high local connectivity, the bond rewiring has a negative role in epidemic spread, i.e., it reduces invasion probability; (ii) in contrast, for epidemics with destructive or weak constructive synergy spreading on networks of arbitrary local connectivity, rewiring helps epidemics to spread; (iii) and, finally, rewiring always enhances the spread of epidemics, independent of synergy, if the local connectivity is low.

Slip statistics of dislocation avalanches under different loading modes

Slowly compressed microcrystals deform via intermittent slip events, observed as displacement jumps or stress drops. Experiments often use one of two loading modes: an increasing applied stress (stress driven, soft), or a constant strain rate (strain driven, hard). In this work we experimentally test the influence of the deformation loading conditions on the scaling behavior of slip events. It is found that these common deformation modes strongly affect time series properties, but not the scaling behavior of the slip statistics when analyzed with a mean-field model. With increasing plastic strain, the slip events are found to be smaller and more frequent when strain driven, and the slip-size distributions obtained for both drives collapse onto the same scaling function with the same exponents. The experimental results agree with the predictions of the used mean-field model, linking the slip behavior under different loading modes.

Effects of variable-state neighborhoods for spreading synergystic processes on lattices

A theoretical framework for the description of susceptible-infected-removed (SIR) spreading processes with synergistic transmission of infection on a lattice is developed. The model incorporates explicitly the effects of time-dependence of the state of the hosts in the neighborhood of transmission events. Exact solution of the model shows that time-dependence of the state of nearest neighbors of recipient hosts is a key factor for synergistic spreading processes. It is demonstrated that the higher the connectivity of a lattice, the more prominent is the effect of synergy on spread.

Crackling noise during failure of alumina under compression: the effect of porosity

We study acoustic emission avalanches during the process of failure of porous alumina samples (Al2O3) under compression. Specimens with different porosities ranging from 30% to 59% have been synthesized from a mixture of fine-grained alumina and graphite. The compressive strength as well as the characteristics of the acoustic activity have been determined. The statistical analysis of the recorded acoustic emission pulses reveals, for all porosities, a broad distribution of energies with a fat tail, compatible with the existence of an underlying critical point. In the region of 35%-55% porosity, the energy distributions of the acoustic emission signals are compatible with a power-law behaviour over two decades in energy with an exponent ϵ = 1.8 ± 0.1.

High junction and twin boundary densities in driven dynamical systems

A novel mechanism for the generation of device materials with very high domain boundary densities is described: we shear the sample in a computer experiment and achieve higher twin densities than in rapid quench. These domain patterns are very stable. Elastically soft materials (image with 6.4$ \times $10(5) atoms) has greater twin densities than hard materials, even for nano-crystals.

Analysis of power-law exponents by maximum-likelihood maps

Maximum-likelihood exponent maps have been studied as a technique to increase the understanding and improve the fit of power-law exponents to experimental and numerical simulation data, especially when they exhibit both upper and lower cutoffs. The use of the technique is tested by analyzing seismological data, acoustic emission data, and avalanches in numerical simulations of the three-dimensional random field Ising model. In the different examples we discuss the nature of the deviations observed in the exponent maps and some relevant conclusions are drawn for the physics behind each phenomenon.

Ripples in a string coupled to Glauber spins

Each oscillator in a linear chain (a string) interacts with a local Ising spin in contact with a thermal bath. These spins evolve according to Glauber dynamics. Below a critical temperature, there appears an equilibrium, time-independent, rippled state in the string that is accompanied by a nonzero spin polarization. On the other hand, the system is shown to form "metastable," nonequilibrium long-lived ripples in the string for slow spin relaxation. The system vibrates rapidly about these quasistationary states, which can be described as snapshots of a coarse-grained stroboscopic map. For moderate observation times, ripples are observed irrespective of the final thermodynamically stable state (rippled or not). Interestingly, the system can be considered as a "minimal" model to understand rippling in clamped graphene sheets.