Source author record

Stefano Zapperi

Stefano Zapperi 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

59works
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

59 published item(s)

preprint2022arXiv

Predicting Creep Failure by Machine Learning -- Which Features Matter?

Spatial and temporal features are studied with respect to their predictive value for failure time prediction in subcritical failure with machine learning (ML). Data are generated from simulations of a novel, brittle random fuse model (RFM), as well as elasto-plastic finite element simulations (FEM) of a stochastic plasticity model with damage, both models considering stochastic thermally activated damage/failure processes in disordered materials. Fuse networks are generated with hierarchical and nonhierarchical architectures. Random forests - a specific ML algorithm - allow us to measure the feature importance through a feature's average error reduction. RFM simulation data are found to become more predictable with increasing system size and temperature. Increasing the load or the scatter in local materials properties has the opposite effect. Damage accumulation in these models proceeds in stochastic avalanches, and statistical signatures such as avalanche rate or magnitude have been discussed in the literature as predictors of incipient failure. However, in the present study such features proved of no measurable use to the ML models, which mostly rely on global or local strain for prediction. This suggests the strain as viable quantity to monitor in future experimental studies as it is accessible via digital image correlation.

preprint2022arXiv

Predicting the failure of two-dimensional silica glasses

Being able to predict the failure of materials based on structural information is a fundamental issue with enormous practical and industrial relevance for the monitoring of devices and components. Thanks to recent advances in deep learning, accurate failure predictions are becoming possible even for strongly disordered solids, but the sheer number of parameters used in the process renders a physical interpretation of the results impossible. Here we address this issue and use machine learning methods to predict the failure of simulated two dimensional silica glasses from their initial undeformed structure. We then exploit Gradient-weighted Class Activation Mapping (Grad-CAM) to build attention maps associated with the predictions, and we demonstrate that these maps are amenable to physical interpretation in terms of topological defects and local potential energies. We show that our predictions can be transferred to samples with different shape or size than those used in training, as well as to experimental images. Our strategy illustrates how artificial neural networks trained with numerical simulation results can provide interpretable predictions of the behavior of experimentally measured structures.

preprint2022arXiv

Size effects in micro and nanoscale materials fracture

Micro and nanoscale materials have remarkable mechanical properties, such as enhanced strength and toughness, but usually display sample-to-sample fluctuations and non-trivial size effects, a nuisance for engineering applications and an intriguing problem for science. Our understanding of size-effects in small-scale materials has progressed considerably in the past few years thanks to a growing number of experimental measurements on carbon based nanomaterials, such as graphene carbon nanotubes, and on crystalline and amorphous micro/nanopillars and micro/nanowires. At the same time, increased computational power allowed atomistic simulations to reach experimentally relevant sample sizes. From the theoretical point of view, the standard analysis and interpretation of experimental and computational data relies on traditional extreme value theories developed decades ago for macroscopic samples, with recent work extending some of the limiting assumptions of the original theories. In this review, we discuss the recent experimental and numerical literature on micro and nanoscale fracture size effects, illustrate existing theories pointing out their advantages and limitations and finally provide a tutorial for analyzing fracture data from micro and nanoscale samples. We discuss a broad spectrum of materials but provide at the same time a unifying theoretical framework that should be helpful for materials scientists working on micro and nanoscale mechanics.

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.

preprint2020arXiv

Automatic Design of Mechanical Metamaterial Actuators

Mechanical metamaterials actuators achieve pre-determined input--output operations exploiting architectural features encoded within a single 3D printed element, thus removing the need of assembling different structural components. Despite the rapid progress in the field, there is still a need for efficient strategies to optimize metamaterial design for a variety of functions. We present a computational method for the automatic design of mechanical metamaterial actuators that combines a reinforced Monte Carlo method with discrete element simulations. 3D printing of selected mechanical metamaterial actuators shows that the machine-generated structures can reach high efficiency, exceeding human-designed structures. We also show that it is possible to design efficient actuators by training a deep neural network, eliminating the need for lengthy mechanical simulations. The elementary actuators devised here can be combined to produce metamaterial machines of arbitrary complexity for countless engineering applications.

preprint2020arXiv

Protein driven lipid domain nucleation in biological membranes

Lipid rafts are heterogeneous dynamic lipid domains of the cell membranes that are involved in several biological processes, like protein and lipids specific transport and signaling. Our understanding of lipid raft formation is still limited, due to the transient and elusive nature of these domains {\it in vivo}, in contrast to the stable phase-separated domains observed in artificial membranes. Inspired by experimental findings highlighting the relevance of transmembrane proteins for lipid rafts, we investigate lipid domain nucleation by coarse-grained molecular dynamics and Ising model simulations. We find that the presence of a transmembrane protein can trigger lipid domain nucleation in a flat membrane from an otherwise mixed lipid phase. Furthermore, we study the role of the lipid domain in the diffusion of the protein showing that its mobility is hindered by the presence of the raft. The results of our coarse-grained molecular dynamics and Ising model simulations thus coherently support the important role played by transmembrane proteins in lipid domain formation and stability.

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.

preprint2020arXiv

Unjamming of active rotators

Active particle assemblies can exhibit a wide range of interesting dynamical phases depending on internal parameters such as density, adhesion strength or self-propulsion. Active self-rotations are rarely studied in this context, although they can be relevant for active matter systems, as we illustrate by analyzing the motion of Chlamydomonas reinhardtii algae under different experimental conditions. Inspired by this example, we simulate the dynamics of a system of interacting active disks endowed with active torques. At low packing fractions, adhesion causes the formation of small rotating clusters, resembling those observed when algae are stressed. At higher densities, the model shows a jamming to unjamming transition promoted by active torques and hindered by adhesion. Our results yield a comprehensive picture of the dynamics of active rotators, providing useful guidance to interpret experimental results in cellular systems where rotations might play a role.

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.

preprint2016arXiv

Deformation and fracture of echinoderm collagen networks

Collagen networks provide the main structural component of most tissues and represent an important ingredient for bio-mimetic materials for bio-medical applications. Here we study the mechanical properties of stiff collagen networks derived from three different echinoderms and show that they exhibit non-linear stiffening followed by brittle fracture. The disordered nature of the network leads to strong sample-to-sample fluctuations in elasticity and fracture strength. We perform numerical simulations of a three dimensional model for the deformation of a cross-linked elastic fibril network which is able to reproduce the macroscopic features of the experimental results and provide insights into the internal mechanics of stiff collagen networks. Our numerical model provides an avenue for the design of collagen membranes with tunable mechanical properties.

preprint2016arXiv

Fluctuations in protein aggregation: Design of preclinical screening for early diagnosis of neurodegenerative disease

Autocatalytic fibril nucleation has recently been proposed to be a determining factor for the spread of neurodegenerative diseases, but the same process could also be exploited to amplify minute quantities of protein aggregates in a diagnostic context. Recent advances in microfluidic technology allow analysis of protein aggregation in micron-scale samples potentially enabling such diagnostic approaches, but the theoretical foundations for the analysis and interpretation of such data are so far lacking. Here we study computationally the onset of protein aggregation in small volumes and show that the process is ruled by intrinsic fluctuations whose volume dependent distribution we also estimate theoretically. Based on these results, we develop a strategy to quantify in silico the statistical errors associated with the detection of aggregate containing samples. Our work opens a new perspective on the forecasting of protein aggregation in asymptomatic subjects.

preprint2016arXiv

Glassy features of crystal plasticity

Crystal plasticity occurs by deformation bursts due to the avalanche-like motion of dislocations. Here we perform extensive numerical simulations of a three-dimensional dislocation dynamics model under quasistatic stress-controlled loading. Our results show that avalanches are power-law distributed, and display peculiar stress and sample size dependence: The average avalanche size grows exponentially with the applied stress, and the amount of slip increases with the system size. These results suggest that intermittent deformation processes in crystalline materials exhibit an extended critical-like phase in analogy to glassy systems, instead of originating from a non-equilibrium phase transition critical point.

preprint2015arXiv

Avalanches, loading and finite size effects in 2D amorphous plasticity: results from a finite element model

Crystalline plasticity is strongly interlinked with dislocation mechanics and nowadays is relatively well understood. Concepts and physical models of plastic deformation in amorphous materials on the other hand - where the concept of linear lattice defects is not applicable - still are lagging behind. We introduce an eigenstrain-based finite element lattice model for simulations of shear band formation and strain avalanches. Our model allows us to study the influence of surfaces and finite size effects on the statistics of avalanches. We find that even with relatively complex loading conditions and open boundary conditions, critical exponents describing avalanche statistics are unchanged, which validates the use of simpler scalar lattice-based models to study these phenomena.

preprint2015arXiv

Crossover behavior in interface depinning

We study the crossover scaling behavior of the height-height correlation function in interface depinning in random media. We analyze experimental data from a fracture experiment and simulate an elastic line model with non-linear couplings and disorder. Both exhibit a crossover between two different universality classes. For the experiment, we fit a functional form to the universal crossover scaling function. For the model, we vary the system size and the strength of the non-linear term, and describe the crossover between the two universality classes with a multiparameter scaling function. Our method provides a general strategy to extract scaling properties in depinning systems exhibiting crossover phenomena.

preprint2015arXiv

Fracture size effects in nanoscale materials: the case of graphene

Nanoscale materials display enhanced strength and toughness but also larger fluctuations and more pronounced size effects with respect to their macroscopic counterparts. Here we study the system size-dependence of the failure strength distribution of a monolayer graphene sheet with a small concentration of vacancies by molecular dynamics simulations. We simulate sheets of varying size encompassing more than three decades and systematically study their deformation as a function of disorder, temperature and loading rate. We generalize the weakest-link theory of fracture size effects to rate and temperature dependent failure and find quantitative agreement with the simulations. Our numerical and theoretical results explain the crossover of the fracture strength distribution between a thermal and rate-dependent regime and a disorder-dominated regime described by extreme value theory.

preprint2015arXiv

Navigation strategies of motor proteins on decorated tracks

Motor proteins display widely different stepping patterns as they move on microtubule tracks, from the deterministic linear or helical motion performed by the protein kinesin to the uncoordinated random steps made by dynein. How these different strategies produce an efficient navigation system needed to ensure correct cellular functioning is still unclear. Here, we show by numerical simulations that deterministic and random motor steps yield different outcomes when random obstacles decorate the microtubule tracks: kinesin moves faster on clean tracks but its motion is strongly hindered on decorated tracks, while dynein is slower on clean tracks but more efficient in avoiding obstacles. Further simulations indicate that dynein's advantage on decorated tracks is due to its ability to step backwards. Our results explain how different navigation strategies are employed by the cell to optimize motor driven cargo transport.

preprint2015arXiv

Osmotic stress affects functional properties of human melanoma cell lines

Understanding the role of microenvironment in cancer growth and metastasis is a key issue for cancer research. Here, we study the effect of osmotic pressure on the functional properties of primary and metastatic melanoma cell lines. In particular, we experimentally quantify individual cell motility and transmigration capability. We then perform a circular scratch assay to study how a cancer cell front invades an empty space. Our results show that primary melanoma cells are sensitive to a low osmotic pressure, while metastatic cells are less. To better understand the experimental results, we introduce and study a continuous model for the dynamics of a cell layer and a stochastic discrete model for cell proliferation and diffusion. The two models capture essential features of the experimental results and allow to make predictions for a wide range of experimentally measurable parameters.

preprint2015arXiv

Overshoot during phenotypic switching of cancer cell populations

The dynamics of tumor cell populations is hotly debated: do populations derive hierarchically from a subpopulation of cancer stem cells (CSCs), or are stochastic transitions that mutate differentiated cancer cells to CSCs important? Here we argue that regulation must also be important. We sort human melanoma cells using three distinct cancer stem cell (CSC) markers - CXCR6, CD271 and ABCG2 - and observe that the fraction of non-CSC-marked cells first overshoots to a higher level and then returns to the level of unsorted cells. This clearly indicates that the CSC population is homeostatically regulated. Combining experimental measurements with theoretical modeling and numerical simulations, we show that the population dynamics of cancer cells is associated with a complex miRNA network regulating the Wnt and PI3K pathways. Hence phenotypic switching is not stochastic, but is tightly regulated by the balance between positive and negative cells in the population. Reducing the fraction of CSCs below a threshold triggers massive phenotypic switching, suggesting that a therapeutic strategy based on CSC eradication is unlikely to succeed.

preprint2015arXiv

Role of the Number of Microtubules in Chromosome Segregation during Cell Division

Faithful segregation of genetic material during cell division requires alignment of chromosomes between two spindle poles and attachment of their kinetochores to each of the poles. Failure of these complex dynamical processes leads to chromosomal instability (CIN), a characteristic feature of several diseases including cancer. While a multitude of biological factors regulating chromosome congression and bi-orientation have been identified, it is still unclear how they are integrated so that coherent chromosome motion emerges from a large collection of random and deterministic processes. Here we address this issue by a three dimensional computational model of motor-driven chromosome congression and bi-orientation during mitosis. Our model reveals that successful cell division requires control of the total number of microtubules: if this number is too small bi-orientation fails, while if it is too large not all the chromosomes are able to congress. The optimal number of microtubules predicted by our model compares well with early observations in mammalian cell spindles. Our results shed new light on the origin of several pathological conditions related to chromosomal instability.

preprint2015arXiv

Volume changes during active shape fluctuations in cells

Cells modify their volume in response to changes in osmotic pressure but it is usually assumed that other active shape variations do not involve significant volume fluctuations. Here we report experiments demonstrating that water transport in and out of the cell is needed for the formation of blebs, commonly observed protrusions in the plasma membrane driven by cortex contraction. We develop and simulate a model of fluid mediated membrane-cortex deformations and show that a permeable membrane is necessary for bleb formation which is otherwise impaired. Taken together our experimental and theoretical results emphasize the subtle balance between hydrodynamics and elasticity in actively driven cell morphological changes.

preprint2014arXiv

Avoidance and Coalescence of Delamination Patterns

Delamination of coatings and thin films from substrates generates a fascinating variety of patterns, from circular blisters to wrinkles and labyrinth domains, in a way that is not completely understood. We report on large-scale numerical simulations of the universal problem of avoidance and coalescence of delamination wrinkles, focusing on a case study of graphene sheets on patterned substrates. By nucleating and growing wrinkles in a controlled way, we are able to characterize how their interactions, mediated by long-range stress fields, determine their formation and morphology. We also study how the interplay between geometry and stresses drives a universal transition from conformation to delamination when sheets are deposited on particle-decorated substrates. Our results are directly applicable to strain engineering of graphene and also uncover universal phenomena observed at all scales, as for example in geomembrane deposition.

preprint2014arXiv

Conformational mechanism for the stability of microtubule-kinetochore attachments

Regulating the stability of microtubule(MT)-kinetochore attachments is fundamental to avoiding mitotic errors and ensure proper chromosome segregation during cell division. While biochemical factors involved in this process have been identified, its mechanics still needs to be better understood. Here we introduce and simulate a mechanical model of MT-kinetochore interactions in which the stability of the attachment is ruled by the geometrical conformations of curling MT-protofilaments entangled in kinetochore fibrils. The model allows us to reproduce with good accuracy in vitro experimental measurements of the detachment times of yeast kinetochores from MTs under external pulling forces. Numerical simulations suggest that geometrical features of MT-protofilaments may play an important role in the switch between stable and unstable attachments.

preprint2014arXiv

Fracture strength: Stress concentration, extreme value statistics and the fate of the Weibull distribution

The fracture strength distribution of materials is often described in terms of the Weibull law which can be derived by using extreme value statistics if elastic interactions are ignored. Here, we consider explicitly the interplay between elasticity and disorder and test the asymptotic validity of the Weibull distribution through numerical simulations of the two-dimensional random fuse model. Even when the local fracture strength follows the Weibull distribution, the global failure distribution is dictated by stress enhancement at the tip of the cracks and sometimes deviates from the Weibull law. Only in the case of a pre-existing power law distribution of crack widths do we find that the failure strength is Weibull distributed. Contrary to conventional assumptions, even in this case, the Weibull exponent can not be simply inferred from the exponent of the initial crack width distribution. Our results thus raise some concerns on the applicability of the Weibull distribution in most practical cases.

preprint2014arXiv

Mechanical Properties of Growing Melanocytic Nevi and the Progression to Melanoma

Melanocytic nevi are benign proliferations that sometimes turn into malignant melanoma in a way that is still unclear from the biochemical and genetic point of view. Diagnostic and prognostic tools are then mostly based on dermoscopic examination and morphological analysis of histological tissues. To investigate the role of mechanics and geometry in the morpholgical dynamics of melanocytic nevi, we study a computation model for cell proliferation in a layered non-linear elastic tissue. Numerical simulations suggest that the morphology of the nevus is correlated to the initial location of the proliferating cell starting the growth process and to the mechanical properties of the tissue. Our results also support that melanocytes are subject to compressive stresses that fluctuate widely in the nevus and depend on the growth stage. Numerical simulations of cells in the epidermis releasing matrix metalloproteinases display an accelerated invasion of the dermis by destroying the basal membrane. Moreover, we suggest experimentally that osmotic stress and collagen inhibit growth in primary melanoma cells while the effect is much weaker in metastatic cells. Knowing that morphological features of nevi might also reflect geometry and mechanics rather than malignancy could be relevant for diagnostic purposes

preprint2014arXiv

Scalar model for frictional precursors dynamics

Recent experiments indicate that frictional sliding occurs by the nucleation of detachment fronts at the contact interface that may appear well before the onset of global sliding. This intriguing precursory activity is not accounted for by traditional friction theories but is extremely important for friction dominated geophysical phenomena such as earthquakes, landslides or avalanches. Here we simulate the onset of slip of a three dimensional elastic body resting on a surface and show that experimentally observed frictional precursors depend in a complex non-universal way on the sample geometry and the loading conditions. Our model satisfies Archard's law and Amontons' first and second laws, reproducing with remarkable precision the real contact area dynamics, the precursors' envelope dynamics prior to the transition to sliding, and the normal and shear internal stress distributions close to the slider-substrate interface. Moreover, it allows to assess which experimental features can be attributed to the elastic equilibrium, and which are attributed to the out-of-equilibrium dynamics, suggesting that precursory activity is an intrinsically quasi-static physical process. A direct calculation of the evolution of the Coulomb stress before and during precursors nucleation shows large variations across the sample, which helps to explain why earthquake forecasting methods based only on accumulated slip and Coulomb stress monitoring are often ineffective.

preprint2013arXiv

Avalanche localization and crossover scaling in amorphous plasticity

We perform large scale simulations of a two dimensional lattice model for amorphous plasticity with random local yield stresses and long-range quadrupolar elastic interactions. We show that as the external stress increases towards the yielding phase transition, the scaling behavior of the avalanches crosses over from mean-field theory to a different universality class. This behavior is associated with strain localization, which significantly depends on the short-range properties of the interaction kernel.

preprint2013arXiv

Avalanches in 2D Dislocation Systems: Plastic Yielding is not Depinning

We study the properties of strain bursts (dislocation avalanches) occurring in two-dimensional discrete dislocation dynamics models under quasistatic stress-controlled loading. Contrary to previous suggestions, the avalanche statistics differs fundamentally from predictions obtained for the depinning of elastic manifolds in quenched random media. Instead, we find an exponent τ=1 of the power-law distribution of slip or released energy, with a cut-off that increases exponentially with the applied stress and diverges with system size at all stresses. These observations demonstrate that the avalanche dynamics of 2D dislocation systems is scale-free at every applied stress and, therefore, can not be envisaged in terms of critical behavior associated with a depinning transition.

preprint2013arXiv

Barkhausen instabilities from labyrinthine magnetic domains

Experimental investigations of the scaling behavior of Barkhausen avalanches in out-of-plane ferromagnetic films yield widely different results for the values of the critical exponents despite similar labyrinthine domain structures, suggesting that universality may not hold for this class of materials. Analyzing a phase field model for magnetic reversal, we show that avalanche scaling is bounded by characteristic lengthscales arising from the competition between dipolar forces and exchange interactions. We compare our results with the experiments and find a good qualitative and quantitative agreement, reconciling apparent contradictions. Finally, we make some prediction, amenable to experimental verification, on the dependence of the avalanches behavior from the film thickness and disorder.

preprint2013arXiv

Damage accumulation in quasi-brittle fracture

The strength of quasi-brittle materials depends on the ensemble of defects inside the sample and on the way damage accumulates before failure. Using large scale numerical simulations of the random fuse model, we investigate the evolution of the microcrack distribution that is directly related to the strength distribution and its size effects. We show that the broadening of the distribution tail originates from the dominating microcracks in each sample and is related to a tendency of crack coalescence that increases with system size. We study how the observed behavior depends on the disorder present in the sample.

preprint2013arXiv

Dislocation mutual interactions mediated by mobile impurities and the conditions for plastic instabilities

Matallic alloys, such as Al or Cu, or mild steel, display plastic instabilities in a well defined range of temperatures and deformation rates, a phenomenon known as the Portevin-Le Chatelelier (PLC) effect. The stick-slip behavior, or serration, typical of this effect is due to the discontinuous motion of dislocations as they interact with solute atoms. Here we study a simple model of interacting dislocations and show how the classical Einstein fluctuation-dissipation relation can be used to define the temperature in a range of model parameters and to construct a phase diagram of serration that can be compared to experimental results. Furthermore, performing analytical calculations and numerically integrating the equations of motion, we clarify the crucial role played by dislocation mutual interactions in serration.

preprint2013arXiv

Growth and form of melanoma cell colonies

We study the statistical properties of melanoma cell colonies grown in vitro by analyzing the results of crystal violet assays at different concentrations of initial plated cells and for different growth times. The distribution of colony sizes is described well by a continuous time branching process. To characterize the shape fluctuations of the colonies, we compute the distribution of eccentricities. The experimental results are compared with numerical results for models of random division of elastic cells, showing that experimental results are best reproduced by restricting cell division to the outer rim of the colony. Our results serve to illustrate the wealth of information that can be extracted by a standard experimental method such as the crystal violet assay.

preprint2013arXiv

Quasi-periodic events in crystal plasticity and the self-organized avalanche oscillator

When external stresses in a system - physical, social or virtual - are relieved through impulsive events, it is natural to focus on the attributes of these avalanches. However, during the quiescent periods in between, stresses may be relieved through competing processes, such as slowly flowing water between earthquakes or thermally activated dislocation flow between plastic bursts. Such unassuming, smooth responses can have dramatic effects on the avalanche properties. Here we report a thorough experimental investigation of slowly compressed Ni microcrystals, covering three orders of magnitude in nominal strain rate, that exhibits unconventional quasi-periodic avalanche bursts and higher critical exponents as the strain rate is decreased. Our analytic and computational study naturally extends dislocation avalanche modeling to incorporate dislocation relaxations and reveals the emergence of the self-organized avalanche oscillator, a novel critical state exhibiting oscillatory approaches toward a depinning critical point. We demonstrate that the predictions of our theory are faithfully exhibited in our experiments.

preprint2013arXiv

Role of the sample thickness in planar crack propagation

We study the effect of the sample thickness in planar crack front propagation in a disordered elastic medium using the random fuse model. We employ different loading conditions and we test their stability with respect to crack growth. We show that the thickness induces characteristic lengths in the stress enhancement factor in front of the crack and in the stress transfer function parallel to the crack. This is reflected by a thickness-dependent crossover scale in the crack front morphology that goes from from multi-scaling to self-affine with exponents in agreement with line depinning models and experiments. Finally, we compute the distribution of crack avalanches which is shown to depend on the thickness and the loading mode.

preprint2013arXiv

Senescent Cells in Growing Tumors: Population Dynamics and Cancer Stem Cells

Tumors are defined by their intense proliferation, but sometimes cancer cells turn senescent and stop replicating. In the stochastic cancer model in which all cells are tumorigenic, senescence is seen as the result of random mutations, suggesting that it could represent a barrier to tumor growth. In the hierarchical cancer model a subset of the cells, the cancer stem cells, divide indefinitely while other cells eventually turn senescent. Here we formulate cancer growth in mathematical terms and obtain predictions for the evolution of senescence. We perform experiments in human melanoma cells which are compatible with the hierarchical model and show that senescence is a reversible process controlled by survivin. We conclude that enhancing senescence is unlikely to provide a useful therapeutic strategy to fight cancer, unless the cancer stem cells are specifically targeted

preprint2013arXiv

Size effects in dislocation depinning models for plastic yield

Typically, the plastic yield stress of a sample is determined from a stress-strain curve by defining a yield strain and reading off the stress required to attain it. However, it is not a priori clear that yield strengths of microscale samples measured this way should display the correct finite size scaling. Here we study plastic yield as a depinning transition of a 1+1 dimensional interface, and consider how finite size effects depend on the choice of yield strain, as well as the presence of hardening and the strength of elastic coupling. Our results indicate that in sufficiently large systems, the choice of yield strain is unimportant, but in smaller systems one must take care to avoid spurious effects.

preprint2013arXiv

Universality classes and crossover scaling of Barkhausen noise in thin films

We study the dynamics of head-to-head domain walls separating in-plane domains in a disordered ferromagnetic thin film. The competition between the domain wall surface tension and dipolar interactions induces a crossover between a rough domain wall phase at short length-scales and a large-scale phase where the walls display a zigzag morphology. The two phases are characterized by different critical exponents for Barkhausen avalanche dynamics that are in quantitative agreement with experimental measurements on MnAs thin films.

preprint2012arXiv

Critical behavior and size dependence in dislocation plasticity

Microcrystals deform very differently from their macroscopic counterparts, displaying a size-dependent yield stress and intermittent plastic strain bursts. Here we demonstrate that size effects and strain bursts may be just two different signatures of fundamental collective phenomena that govern the dynamics of dislocation systems on all scales. We confirm this conjecture by demonstrating agreement between the results of 2D and 3D dislocation dynamics simulations and experimental results on compressed micropillars. Our results indicate that many peculiarities of small-sample plasticity may arise from the finite-size scaling properties of collective dislocation dynamics, rather than from changes in the physical mechanisms governing dislocation motion.

preprint2012arXiv

From damage percolation to crack nucleation through finite size criticality

We present a unified theory of fracture in disordered brittle media that reconciles apparently conflicting results reported in the literature. Our renormalization group based approach yields a phase diagram in which the percolation fixed point, expected for infinite disorder, is unstable for finite disorder and flows to a zero-disorder nucleation-type fixed point, thus showing that fracture has mixed first order and continuous character. In a region of intermediate disorder and finite system sizes, we predict a crossover with mean-field avalanche scaling. We discuss intriguing connections to other phenomena where critical scaling is only observed in finite size systems and disappears in the thermodynamic limit.

preprint2011arXiv

Avalanche Spatial Structure and Multivariable Scaling Functions; Sizes, Heights, Widths, and Views through Windows

We introduce a systematic method for extracting multivariable universal scaling functions and critical exponents from data. We exemplify our insights by analyzing simulations of avalanches in an interface using simulations from a driven quenched Kardar-Parisi-Zhang (qKPZ) equation. We fully characterize the spatial structure of these avalanches- we report universal scaling functions for size, height and width distributions, and also local front heights. Furthermore, we resolve a problem that arises in many imaging experiments of crackling noise and avalanche dynamics, where the observed distributions are strongly distorted by a limited field of view. Through artificially windowed data, we show these distributions and their multivariable scaling functions may be written in terms of two control parameters, the window size and the characteristic length scale of the dynamics. For the entire system and the windowed distributions we develop accurate parameterizations for the universal scaling functions, including corrections to scaling and systematic error bars, facilitated by a novel software environment SloppyScaling.

preprint2011arXiv

Dielectric breakdown and avalanches at non-equilibrium metal-insulator transitions

Motivated by recent experiments on the finite temperature Mott transition in VO2 films, we propose a classical coarse-grained dielectric breakdown model where each degree of freedom represents a nanograin which transitions from insulator to metal with increasing temperature and voltage at random thresholds due to quenched disorder. We describe the properties of the resulting non-equilibrium metal-insulator transition and explain the universal characteristics of the resistance jump distribution. We predict that by tuning voltage, another critical point is approached, which separates a phase of "bolt"-like avalanches from percolation-like ones.

preprint2011arXiv

Dynamic instability of a growing adsorbed polymorphic filament

The intermittent transition between slow growth and rapid shrinkage in polymeric assemblies is termed dynamic instability, a feature observed in a variety of biochemically distinct assemblies including microtubules, actin and their bacterial analogs. The existence of this labile phase of a polymer has many functional consequences in cytoskeletal dynamics, and its repeated appearance suggests that it is relatively easy to evolve. Here, we consider the minimal ingredients for the existence of dynamic instability by considering a single polymorphic filament that grows by binding to a substrate, undergoes a conformation change, and may unbind as a consequence of the residual strains induced by this change. We identify two parameters that control the phase space of possibilities for the filament: a structural mechanical parameter that characterizes the ratio of the bond strengths along the filament to those with the substrate (or equivalently the ratio of longitudinal to lateral interactions in an assembly), and a kinetic parameter that characterizes the ratio of time scales for growth and conformation change. In the deterministic limit, these parameters serve to demarcate a region of uninterrupted growth from that of collapse. However, in the presence of disorder in either the structural or the kinetic parameter the growth and collapse phases can coexist where the filament can grow slowly, shrink rapidly, and transition between these phases, thus exhibiting dynamic instability. We exhibit the window for the existence of dynamic instability in a phase diagram that allows us to quantify the evolvability of this labile phase.

preprint2011arXiv

Fracture strength of disordered media: Universality, interactions and tail asymptotics

We study the asymptotic properties of fracture strength distributions of disordered elastic media by a combination of renormalization group, extreme value theory, and numerical simulation. We investigate the validity of the `weakest-link hypothesis' in the presence of realistic long-ranged interactions in the random fuse model. Numerical simulations indicate that the fracture strength is well described by the Duxbury-Leath-Beale (DLB) distribution which is shown to flow asymptotically to the Gumbel distribution. We explore the relation between the extreme value distributions and the DLB type asymptotic distributions, and show that the universal extreme value forms may not be appropriate to describe the non-universal low-strength tail.

preprint2010arXiv

Avalanches and clusters in planar crack front propagation

We study avalanches in a model for a planar crack propagating in a disordered medium. Due to long-range interactions, avalanches are formed by a set of spatially disconnected local clusters, the sizes of which are distributed according to a power law with an exponent $τ_{a}=1.5$. We derive a scaling relation $τ_a=2τ-1$ between the local cluster exponent $τ_a$ and the global avalanche exponent $τ$. For length scales longer than a cross-over length proportional to the Larkin length, the aspect ratio of the local clusters scales with the roughness exponent of the line model. Our analysis provides an explanation for experimental results on planar crack avalanches in Plexiglas plates, but the results are applicable also to other systems with long-range interactions.

preprint2010arXiv

Beyond power laws: Universality in the average avalanche shape

We report the measurement of multivariable scaling functions for the temporal average shape of Barkhausen noise avalanches, and show that they are consistent with the predictions of simple mean-field theories. We bypass the confounding factors of time-retarded interactions (eddy currents) by measuring thin permal- loy films, and bypass thresholding effects and amplifier distortions by applying Wiener deconvolution. We find experimental shapes that are approximately symmetric, and track the evolution of the scaling function. We solve a mean- field theory for the magnetization dynamics and calculate the form of the scaling function in the presence of a demagnetizing field and a finite field ramp-rate, yielding quantitative agreement with the experiment.

preprint2010arXiv

From Brittle to Ductile Fracture in Disordered Materials

We introduce a lattice model able to describe damage and yielding in heterogeneous materials ranging from brittle to ductile ones. Ductile fracture surfaces, obtained when the system breaks once the strain is completely localized, are shown to correspond to minimum energy surfaces. The similarity of the resulting fracture paths to the limits of brittle fracture or minimum energy surfaces is quantified. The model exhibits a smooth transition from brittleness to ductility. The dynamics of yielding exhibits avalanches with a power-law distribution.

preprint2010arXiv

Roughness and multiscaling of planar crack fronts

We consider numerically the roughness of a planar crack front within the long-range elastic string model, with a tunable disorder correlation length $ξ$. The problem is shown to have two important length scales, $ξ$ and the Larkin length $L_c$. Multiscaling of the crack front is observed for scales below $ξ$, provided that the disorder is strong enough. The asymptotic scaling with a roughness exponent $ζ\approx 0.39$ is recovered for scales larger than both $ξ$ and $L_c$. If $L_c > ξ$, these regimes are separated by a third regime characterized by the Larkin exponent $ζ_L \approx 0.5$. We discuss the experimental implications of our results.

preprint2009arXiv

Effect of dipolar interactions for domain wall dynamics in magnetic thin films

We study the effect of long range dipolar forces on the dynamics and morphology of domain walls in magnetic thin films by numerical simulations of the spin-1 random field Ising model. By studying the size distribution of avalanches of domain wall motion arising as a response to quasistatic external driving, we observe a cross-over from the case dominated by short range interactions to another universality class where the long range dipolar forces become important. This crossover is accompanied with a change of the domain wall morphology from a rough wall to walls with zigzag structure.

preprint2009arXiv

Modeling domain wall dynamics in thin magnetic strips with disorder

We present a line-based model of transverse domain walls in thin magnetic strips, to study the effect of bulk disorder on the domain wall dynamics within the thermally activated creep regime. The creep velocity is found to exhibit a non-linear dependence on both applied magnetic fields and electric currents, characterized by similar creep exponents for both forms of the external drive. We discuss briefly the significance of the inherently stochastic thermally activated domain wall motion from the point of view of spintronics applications, where it generally is essential to be able to control the domain wall displacement in a deterministic manner.

preprint2009arXiv

Optimization and plasticity in disordered media

We study the plastic yielding of disordered media using the perfectly plastic random fuse model. The yield surfaces are shown to be different from those obtained minimizing the sum of the local yield thresholds, i.e. the so-called minimum 'energy' surfaces. As a result, the global yield stress is lower than expected from naive optimization and the difference persists as the sample size increases. At variance with minimum energy surfaces, height-height fluctuations of yield surfaces exhibit multiscaling. We provide a theoretical argument that explains how this behavior arises from the very different nature of the optimization problem in both cases.

preprint2009arXiv

Size effects in statistical fracture

We review statistical theories and numerical methods employed to consider the sample size dependence of the failure strength distribution of disordered materials. We first overview the analytical predictions of extreme value statistics and fiber bundle models and discuss their limitations. Next, we review energetic and geometric approaches to fracture size effects for specimens with a flaw. Finally, we overview the numerical simulations of lattice models and compare with theoretical models.

preprint2009arXiv

Suppression of friction by mechanical vibrations

Mechanical vibrations are known to affect frictional sliding and the associated stick-slip patterns causing sometimes a drastic reduction of the friction force. This issue is relevant for applications in nanotribology and to understand earthquake triggering by small dynamic perturbations . We study the dynamics of repulsive particles confined between a horizontally driven top plate and a vertically oscillating bottom plate. Our numerical results show a suppression of the high dissipative stick-slip regime in a well defined range of frequencies that depends on the vibrating amplitude, the normal applied load, the system inertia and the damping constant. We propose a theoretical explanation of the numerical results and derive a phase diagram indicating the region of parameter space where friction is suppressed. Our results allow to define better strategies for the mechanical control of friction.

preprint2005arXiv

Signature of effective mass in crackling noise asymmetry

Crackling noise is a common feature in many dynamic systems [1-9], the most familiar instance of which is the sound made by a sheet of paper when crumpled into a ball. Although seemingly random, this noise contains fundamental information about the properties of the system in which it occurs. One potential source of such information lies in the asymmetric shape of noise pulses emitted by a diverse range of noisy systems [8-12], but the cause of this asymmetry has lacked explanation [1]. Here we show that the leftward asymmetry observed in the Barkhausen effect [2] - the noise generated by the jerky motion of domain walls as they interact with impurities in a soft magnet - is a direct consequence of a magnetic domain wall's negative effective mass. As well as providing a means of determining domain wall effective mass from a magnet's Barkhausen noise our work suggests an inertial explanation for the origin of avalanche asymmetries in crackling noise phenomena more generally.

preprint2003arXiv

Tearing transition and plastic flow in superconducting thin films

A new class of artificial atoms, such as synthetic nanocrystals or vortices in superconductors, naturally self-assemble into ordered arrays. This property makes them applicable to the design of novel solids, and devices whose properties often depend on the response of such assemblies to the action of external forces. Here we study the transport properties of a vortex array in the Corbino disk geometry by numerical simulations. In response to an injected current in the superconductor, the global resistance associated to vortex motion exhibits sharp jumps at two threshold current values. The first corresponds to a tearing transition from rigid rotation to plastic flow, due to the reiterative nucleation around the disk center of neutral dislocation pairs that unbind and glide across the entire disk. After the second jump, we observe a smoother plastic phase proceeding from the coherent glide of a larger number of dislocations arranged into radial grain boundaries.

preprint2000arXiv

Absorbing-state phase transitions in fixed-energy sandpiles

We study sandpile models as closed systems, with conserved energy density $ζ$ playing the role of an external parameter. The critical energy density, $ζ_c$, marks a nonequilibrium phase transition between active and absorbing states. Several fixed-energy sandpiles are studied in extensive simulations of stationary and transient properties, as well as the dynamics of roughening in an interface-height representation. Our primary goal is to identify the universality classes of such models, in hopes of assessing the validity of two recently proposed approaches to sandpiles: a phenomenological continuum Langevin description with absorbing states, and a mapping to driven interface dynamics in random media. Our results strongly suggest that there are at least three distinct universality classes for sandpiles.

preprint2000arXiv

Paths to Self-Organized Criticality

We present a pedagogical introduction to self-organized criticality (SOC), unraveling its connections with nonequilibrium phase transitions. There are several paths from a conventional critical point to SOC. They begin with an absorbing-state phase transition (directed percolation is a familiar example), and impose supervision or driving on the system; two commonly used methods are extremal dynamics, and driving at a rate approaching zero. We illustrate this in sandpiles, where SOC is a consequence of slow driving in a system exhibiting an absorbing-state phase transition with a conserved density. Other paths to SOC, in driven interfaces, the Bak-Sneppen model, and self-organized directed percolation, are also examined. We review the status of experimental realizations of SOC in light of these observations.

preprint2000arXiv

Sandpiles and absorbing-state phase transitions: recent results and open problems

We review some recent results on the relations between sandpiles and a class of absorbing state phase transitions. We use the concept of fixed energy sandpiles (FES), in which external driving and dissipation are absent. FES are shown to exhibit an absorbing state transition with critical properties coinciding with those of the corresponding sandpile model. We propose a set of Langevin equations capturing the relevant features of this transition. These equations characterize the universality class of systems with an infinite number of absorbing states and a static conserved field coupled to the order parameter. Different models in this class are identified, and strong evidence is presented showing that the Manna sandpile, as well as some other stochastic sandpiles, belong in this universality class. Finally some open problems and questions are discussed.

preprint1996arXiv

Sandpile Model on Sierpinski Gasket Fractal

We investigate the sandpile model on the two--dimensional Sierpinski gasket fractal. We find that the model displays novel critical behavior, and we analyze the distribution functions of avalanche sizes, lifetimes and topplings and calculate the associated critical exponents $τ= 1.51 \pm 0.04$, $α= 1.63 \pm 0.04$ and $μ= 1.36 \pm 0.04$. The avalanche size distribution shows power law behavior modulated by logarithmic oscillations which can be related to the discrete scale invariance of the underlying lattice. Such a distribution can be formally described by introducing a complex scaling exponent $τ^{*} \equiv τ+ i δ$, where the real part $τ$ corresponds to the power law and the imaginary part $δ$ is related to the period of the logarithmic oscillations.