Source author record

Alberto Rosso

Alberto Rosso 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

53works
14topics
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

53 published item(s)

preprint2022arXiv

Avalanches and deformation in glasses and disordered systems

In this chapter, we discuss avalanches in glasses and disordered systems, and the macroscopic dynamical behavior that they mediate. We briefly review three classes of systems where avalanches are observed: depinning transition of disordered interfaces, yielding of amorphous materials, and the jamming transition. Without extensive formalism, we discuss results gleaned from theoretical approaches -- mean-field theory, scaling and exponent relations, the renormalization group, and a few results from replica theory. We focus both on the remarkably sophisticated physics of avalanches and on relatively new approaches to the macroscopic flow behavior exhibited past the depinning/yielding transition.

preprint2022arXiv

Clusters in an epidemic model with long-range dispersal

In presence of long range dispersal, epidemics spread in spatially disconnected regions known as clusters. Here, we characterize exactly their statistical properties in a solvable model, in both the supercritical (outbreak) and critical regimes. We identify two diverging length scales, corresponding to the bulk and the outskirt of the epidemic. We reveal a nontrivial critical exponent that governs the cluster number, the distribution of their sizes and of the distances between them. We also discuss applications to depinning avalanches with long range elasticity.

preprint2022arXiv

The Fate of Shear-Oscillated Amorphous Solids

The behavior of shear-oscillated amorphous materials is studied using a coarse-grained model. Samples are prepared at different degrees of annealing and then subject to athermal and quasistatic oscillatory deformations at various fixed amplitudes. The steady-state reached after several oscillations is fully determined by the initial preparation and the oscillation amplitude, as seen from stroboscopic stress and energy measurements. Under small oscillations, poorly annealed materials display shear-annealing, while ultra-stabilized materials are insensitive to them. Yet, beyond a critical oscillation amplitude, both kind of materials display a discontinuous transition to the same mixed state composed by a fluid shear-band embedded in a marginal solid. Quantitative relations between uniform shear and the steady-state reached with this protocol are established. The transient regime characterizing the growth and the motion of the shear band is also studied.

preprint2020arXiv

Active Brownian Motion in Two Dimensions

We study the dynamics of a single active Brownian particle (ABP) in two spatial dimensions. The ABP has an intrinsic time scale $D_R^{-1}$ set by the rotational diffusion constant $D_R$. We show that, at short-times $t \ll D_R^{-1}$, the presence of `activness' results in a strongly anisotropic and non-diffusive dynamics in the $(xy)$ plane. We compute exactly the marginal distributions of the $x$ and $y$ position coordinates along with the radial distribution, which are all shown to be non-Brownian. In addition, we show that, at early times, the ABP has anomalous first-passage properties, characterized by non-Brownian exponents.

preprint2020arXiv

Current fluctuations in non-interacting run-and-tumble particles in one-dimension

We present a general framework to study the distribution of the flux through the origin up to time $t$, in a non-interacting one-dimensional system of particles with a step initial condition with a fixed density $ρ$ of particles to the left of the origin. We focus principally on two cases: (i) when the particles undergo diffusive dynamics (passive case) and (ii) run-and-tumble dynamics for each particle (active case). In analogy with disordered systems, we consider the flux distribution both for the annealed and the quenched initial conditions, for the passive and active particles. In the annealed case, we show that, for arbitrary particle dynamics, the flux distribution is a Poissonian with a mean $μ(t)$ that we compute exactly in terms of the Green's function of the single particle dynamics. For the quenched case, we show that, for the run-and-tumble dynamics, the quenched flux distribution takes an anomalous large deviation form at large times $P_{\rm qu}(Q,t) \sim \exp\left[-ρ\, v_0\, γ\, t^2 ψ_{\rm RTP}\left(\frac{Q}{ρv_0\,t} \right) \right]$, where $γ$ is the rate of tumbling and $v_0$ is the ballistic speed between two successive tumblings. In this paper, we compute the rate function $ψ_{\rm RTP}(q)$ and show that it is nontrivial. Our method also gives access to the probability of the rare event that, at time $t$, there is no particle to the right of the origin. For diffusive and run-and-tumble dynamics, we find that this probability decays with time as a stretched exponential, $\sim \exp(-c\, \sqrt{t})$ where the constant $c$ can be computed exactly. We verify our results for these large deviations by using an importance sampling Monte-Carlo method.

preprint2020arXiv

Spatial clustering of depinning avalanches in presence of long-range interactions

Disordered elastic interfaces display avalanche dynamics at the depinning transition. For short-range interactions, avalanches correspond to compact reorganizations of the interface well described by the depinning theory. For long-range elasticity, an avalanche is a collection of spatially disconnected clusters. In this paper we determine the scaling properties of the clusters and relate them to the roughness exponent of the interface. The key observation of our analysis is the identification of a Bienaym{é}-Galton-Watson process describing the statistics of the cluster number. Our work has a concrete importance for experimental applications where the cluster statistics is a key probe of avalanche dynamics.

preprint2020arXiv

Stochastic growth in time dependent environments

We study the Kardar-Parisi-Zhang (KPZ) growth equation in one dimension with a noise variance $c(t)$ depending on time. We find that for $c(t)\propto t^{-α}$ there is a transition at $α=1/2$. When $α>1/2$, the solution saturates at large times towards a non-universal limiting distribution. When $α<1/2$ the fluctuation field is governed by scaling exponents depending on $α$ and the limiting statistics are similar to the case when $c(t)$ is constant. We investigate this problem using different methods: (1) Elementary changes of variables mapping the time dependent case to variants of the KPZ equation with constant variance of the noise but in a deformed potential (2) An exactly solvable discretization, the log-gamma polymer model (3) Numerical simulations.

preprint2020arXiv

The influence of the brittle-ductile transition zone on aftershock and foreshock occurrence

Aftershock occurrence is characterized by scaling behaviors with quite universal exponents. At the same time, deviations from universality have been proposed as a tool to discriminate aftershocks from foreshocks. Here we show that the change in rheological behavior of the crust, from velocity weakening to velocity strengthening, represents a viable mechanism to explain statistical features of both aftershocks and foreshocks. More precisely, we present a model of the seismic fault described as a velocity weakening elastic layer coupled to a velocity strengthening visco-elastic layer. We show that the statistical properties of aftershocks in instrumental catalogs are recovered at a quantitative level, quite independently of the value of model parameters. We also find that large earthquakes are often anticipated by a preparatory phase characterized by the occurrence of foreshocks. Their magnitude distribution is significantly flatter than the aftershock one, in agreement with recent results for forecasting tools based on foreshocks.

preprint2020arXiv

Universal scaling of the velocity field in crack front propagation

The propagation of a crack front in disordered materials is jerky and characterized by bursts of activity, called avalanches. These phenomena are the manifestation of an out-of-equilibrium phase transition originated by the disorder. As a result avalanches display universal scalings which are however difficult to characterize in experiments at finite drive. Here we show that the correlation functions of the velocity field along the front allow to extract the critical exponents of the transition and to identify the universality class of the system. We employ these correlations to characterize the universal behavior of the transition in simulations and in an experiment of crack propagation. This analysis is robust, efficient and can be extended to all systems displaying avalanche dynamics.

preprint2019arXiv

Exact stationary state of a run-and-tumble particle with three internal states in a harmonic trap

We study the motion of a one-dimensional run-and-tumble particle with three discrete internal states in the presence of a harmonic trap of stiffness $μ.$ The three internal states, corresponding to positive, negative and zero velocities respectively, evolve following a jump process with rate $γ$. We compute the stationary position distribution exactly for arbitrary values of $μ$ and $γ$ which turns out to have a finite support on the real line. We show that the distribution undergoes a shape-transition as $β=γ/μ$ is changed. For $β<1,$ the distribution has a double-concave shape and shows algebraic divergences with an exponent $(β-1)$ both at the origin and at the boundaries. For $β>1,$ the position distribution becomes convex, vanishing at the boundaries and with a single, finite, peak at the origin. We also show that for the special case $β=1,$ the distribution shows a logarithmic divergence near the origin while saturating to a constant value at the boundaries.

preprint2016arXiv

Nuclear Magnetic Resonance studies of DNP-ready trehalose obtained by solid state mechanochemical amorphization

$^1$H nuclear spin-lattice relaxation and Dynamic Nuclear Polarization (DNP) have been studied in amorphous samples of trehalose sugar doped with TEMPO radicals by means of mechanical milling, in the 1.6 K $÷$ 4.2 K temperature range. The radical concentration was varied between 0.34 and 0.81 $\%$. The highest polarization of 15 \% at 1.6 K, observed in the sample with concentration $0.50 \%$, is of the same order of magnitude of that reported in standard frozen solutions with TEMPO. The temperature and concentration dependence of the spin-lattice relaxation rate $1/T_{\text{1}}$, dominated by the coupling with the electron spins, were found to follow power laws with an exponent close to $3$ in all samples. The observed proportionality between $1/T_{\text{1}}$ and the polarization rate $1/T_{\text{pol}}$, with a coefficient related to the electron polarization, is consistent with the presence of Thermal Mixing (TM) and a good contact between the nuclear and the electron spins. At high electron concentration additional relaxation channels causing a decrease in the nuclear polarization must be considered. These results provide further support for a more extensive use of amorphous DNP-ready samples, obtained by means of comilling, in dissolution DNP experiments and possibly for $\textit{in vivo}$ metabolic imaging.

preprint2016arXiv

Occupation time statistics of the random acceleration model

The random acceleration model is one of the simplest non-Markovian stochastic systems and has been widely studied in connection with applications in physics and mathematics. However, the occupation time and related properties are non-trivial and not yet completely understood. In this paper we consider the occupation time $T_+$ of the one-dimensional random acceleration model on the positive half-axis. We calculate the first two moments of $T_+$ analytically and also study the statistics of $T_+$ with Monte Carlo simulations. One goal of our work was to ascertain whether the occupation time $T_+$ and the time $T_m$ at which the maximum of the process is attained are statistically equivalent. For regular Brownian motion the distributions of $T_+$ and $T_m$ coincide and are given by Lévy's arcsine law. We show that for randomly accelerated motion the distributions of $T_+$ and $T_m$ are quite similar but not identical. This conclusion follows from the exact results for the moments of the distributions and is also consistent with our Monte Carlo simulations.

preprint2016arXiv

Scaling Description of Non-Local Rheology

Non-locality is crucial to understand the plastic flow of an amorphous material, and has been successfully described by the fluidity, along with a cooperativity length scale ξ. We demonstrate, by applying the scaling hypothesis to the yielding transition, that non-local effects in non-uniform stress configurations can be explained within the framework of critical phenomena. From the scaling description, scaling relations between different exponents are derived, and collapses of strain rate profiles are made both in shear driven and pressure driven flow. We find that the cooperative length in non-local flow is governed by the same correlation length in finite dimensional homogeneous flow, excluding the mean field exponents. We also show that non-locality also affects the finite size scaling of the yield stress, especially the large finite size effects observed in pressure driven flow. Our theoretical results are nicely verified by the elasto-plastic model, and experimental data.

preprint2016arXiv

Thermalization and many-body localization in systems under dynamic nuclear polarization

We study the role of dipolar interactions in the standard protocol used to achieve dynamic nuclear polarization (DNP). In the so-called spin-temperature regime, where the interactions establish an effective thermodynamic behavior in the out-of-equilibrium stationary state, we provide numerical predictions for the level of hyperpolarization. We show that nuclear spins equilibrate to the effective spin-temperature established among the electron spins of radicals, as expected from the quantum theory of thermalization. Moreover, we present an analytical technique to estimate the spin temperature, and thus, the nuclear hyperpolarization in the steady state, as a function of interaction strength and quenched disorder. This reproduces both our numerical data and experimental results. Our central finding is that the nuclear hyperpolarization increases steadily upon reducing the interaction strength (by diluting the radical density). Interestingly, the highest polarization is reached at a point where the establishment of a spin temperature is just about to break down due to the incipient many-body localization transition in the electron spin system.

preprint2015arXiv

Conformal invariance of loop ensembles under Kardar-Parisi-Zhang dynamics

We study scaling properties of the honeycomb fully packed loop ensemble associated with a lozenge tiling model of rough surface, when the latter is driven out of equilibrium by Kardar-Parisi-Zhang (KPZ) type dynamics. We show numerically that conformal invariance and signatures of critical percolation appear in the stationary KPZ state. In terms of the two-component Coulomb gas description of the Edwards-Wilkinson stationary state, our finding is understood as the invariance of one component under the effect of the non-linear KPZ term. On the other hand, we show a breaking of conformal invariance when the level lines of the other component are considered.

preprint2015arXiv

Criticality in the approach to failure in amorphous solids

Failure of amorphous solids is fundamental to various phenomena, including landslides and earthquakes. Recent experiments indicate that highly plastic regions form elongated structures that are especially apparent near the maximal shear stress $Σ_{\max}$ where failure occurs. This observation suggested that $Σ_{\max}$ acts as a critical point where the length scale of those structures diverges, possibly causing macroscopic transient shear bands. Here we argue instead that the entire solid phase ($Σ<Σ_{\max}$) is critical, that plasticity always involves system-spanning events, and that their magnitude diverges at $Σ_{\max}$ independently of the presence of shear bands. We relate the statistics and fractal properties of these rearrangements to an exponent $θ$ that captures the stability of the material, which is observed to vary continuously with stress, and we confirm our predictions in elastoplastic models.

preprint2015arXiv

Dynamic nuclear polarization and the paradox of Quantum Thermalization

Dynamic Nuclear Polarization (DNP) is to date the most effective technique to increase the nuclear polarization up to a factor $100,000$ opening disruptive perspectives for medical applications. In DNP, the nuclear spins are driven to an - out of equilibrium - hyperpolarized state by microwave saturation of the electron spins in interaction with them. Here we show that the electron dipolar interactions compete with the local magnetic fields resulting in two distinct dynamical phases: for strong interactions the electron spins equilibrate to an extremely low effective temperature that boosts DNP efficiency. For weak interaction this spin temperature is not defined and the polarization profile has an 'hole burning' shape characteristic of the non interacting case. The study of the many-body eigenstates reveals that these two phases are intimately related to the problem of thermalization in closed quantum systems where breaking of ergodicity is expected varying the strength of the interactions.

preprint2015arXiv

Extreme value statistics of 2d Gaussian Free Field: effect of finite domains

We study minima statistics of the 2d Gaussian Free Field on circles in the unit disk with Dirichlet boundary condition. Free energy distributions of the associated Random Energy models are exactly calculated in the high temperature phase, and shown to satisfy the duality property, which enables us to predict the minima distribution by assuming the freezing scenario. Numerical tests are provided. Related questions concerning the GFF on a sphere are also considered.

preprint2015arXiv

Frictional dynamics of viscoelastic solids driven on a rough surface

We study the effect of viscoelastic dynamics on the frictional properties of a (mean field) spring-block system pulled on a rough surface by an external drive. When the drive moves at constant velocity V, two dynamical regimes are observed: at fast driving, above a critical threshold Vc, the system slides at the drive velocity and displays a friction force with velocity weakening. Below Vc the steady sliding becomes unstable and a stick-slip regime sets in. In the slide-hold-slide driving protocol, a peak of the friction force appears after the hold time and its amplitude increases with the hold duration. These observations are consistent with the frictional force encoded phenomenologically in the rate-and-state equations. Our model gives a microscopical basis for such macroscopic description.

preprint2015arXiv

Role of the glassy dynamics and thermal mixing in the dynamic nuclear polarization and relaxation mechanisms of pyruvic acid

The temperature dependence of $^1$H and $^{13}$C nuclear spin-lattice relaxation rate $1/T_1$ has been studied in the 1.6 K - 4.2 K temperature range in pure pyruvic acid and in pyruvic acid containing trityl radicals at a concentration of 15 mM. The temperature dependence of $1/T_1$ is found to follow a quadratic power law for both nuclei in the two samples. Remarkably the same temperature dependence is displayed also by the electron spin-lattice relaxation rate $1/T_{1e}$ in the sample containing radicals. These results are explained by considering the effect of the structural dynamics on the relaxation rates in pyruvic acid. Dynamic nuclear polarization experiments show that below 4 K the $^{13}$C build up rate scales with $1/T_{\text{1e}}$, in analogy to $^{13}$C $1/T_1$ and consistently with a thermal mixing scenario where all the electrons are collectively involved in the dynamic nuclear polarization process and the nuclear spin reservoir is in good thermal contact with the electron spin system.

preprint2015arXiv

The critical catastrophe revisited

The neutron population in a prototype model of nuclear reactor can be described in terms of a collection of particles confined in a box and undergoing three key random mechanisms: diffusion, reproduction due to fissions, and death due to absorption events. When the reactor is operated at the critical point, and fissions are exactly compensated by absorptions, the whole neutron population might in principle go to extinction because of the wild fluctuations induced by births and deaths. This phenomenon, which has been named critical catastrophe, is nonetheless never observed in practice: feedback mechanisms acting on the total population, such as human intervention, have a stabilizing effect. In this work, we revisit the critical catastrophe by investigating the spatial behaviour of the fluctuations in a confined geometry. When the system is free to evolve, the neutrons may display a wild patchiness (clustering). On the contrary, imposing a population control on the total population acts also against the local fluctuations, and may thus inhibit the spatial clustering. The effectiveness of population control in quenching spatial fluctuations will be shown to depend on the competition between the mixing time of the neutrons (i.e., the average time taken for a particle to explore the finite viable space) and the extinction time.

preprint2015arXiv

Viscoelastic Effects in Avalanche Dynamics: A Key to Earthquake Statistics

In many complex systems a continuous input of energy over time can be suddenly relaxed in the form of avalanches. Conventional avalanche models disregard the possibility of internal dynamical effects in the inter-avalanche periods, and thus miss basic features observed in some real systems. We address this issue by studying a model with viscoelastic relaxation, showing how coherent oscillations of the stress field can emerge spontaneously. Remarkably, these oscillations generate avalanche patterns that are similar to those observed in seismic phenomena.

preprint2014arXiv

Clustering of branching Brownian motions in confined geometries

We study the evolution of a collection of individuals subject to Brownian diffusion, reproduction and disappearance. In particular, we focus on the case where the individuals are initially prepared at equilibrium within a confined geometry. Such systems are widespread in physics and biology and apply for instance to the study of neutron populations in nuclear reactors and the dynamics of bacterial colonies, only to name a few. The fluctuations affecting the number of individuals in space and time may lead to a strong patchiness, with particles clustered together. We show that the analysis of this peculiar behaviour can be rather easily carried out by resorting to a backward formalism based on the Green's function, which allows the key physical observables, namely, the particle concentration and the pair correlation function, to be explicitly derived.

preprint2014arXiv

Diffusion in periodic, correlated random forcing landscapes

We study the dynamics of a Brownian particle in a strongly correlated quenched random potential defined as a periodically-extended (with period $L$) finite trajectory of a fractional Brownian motion with arbitrary Hurst exponent $H \in (0,1)$. While the periodicity ensures that the ultimate long-time behavior is diffusive, the generalised Sinai potential considered here leads to a strong logarithmic confinement of particle trajectories at intermediate times. These two competing trends lead to dynamical frustration and result in a rich statistical behavior of the diffusion coefficient $D_L$: Although one has the typical value $D^{\rm typ}_L \sim \exp(-βL^H)$, we show via an exact analytical approach that the positive moments ($k>0$) scale like $\langle D^k_L \rangle \sim \exp{[-c' (k βL^{H})^{1/(1+H)}]}$, and the negative ones as $\langle D^{-k}_L \rangle \sim \exp(a' (k βL^{H})^2)$, $c'$ and $a'$ being numerical constants and $β$ the inverse temperature. These results demonstrate that $D_L$ is strongly non-self-averaging. We further show that the probability distribution of $D_L$ has a log-normal left tail and a highly singular, one-sided log-stable right tail reminiscent of a Lifshitz singularity.

preprint2014arXiv

Relevance of electron spin dissipative processes to dynamic nuclear polarization via thermal mixing

The available theoretical approaches aiming at describing Dynamic Nuclear spin Polarization (DNP) in solutions containing molecules of biomedical interest and paramagnetic centers are not able to model the behaviour observed upon varying the concentration of trityl radicals or the polarization enhancement caused by moderate addition of gadolinium complexes. In this manuscript, we first show experimentally that the nuclear steady state polarization reached in solutions of pyruvic acid with 15 mM trityl radicals is substantially independent from the average internuclear distance. This evidences a leading role of electron (over nuclear) spin relaxation processes in determining the ultimate performances of DNP. Accordingly, we have devised a variant of the Thermal Mixing model for inhomogenously broadened electron resonance lines which includes a relaxation term describing the exchange of magnetic anisotropy energy of the electron spin system with the lattice. Thanks to this additional term, the dependence of the nuclear polarization on the electron concentration can be properly accounted for. Moreover, the model predicts a strong increase of the final polarization on shortening the electron spin-lattice relaxation time, providing a possible explanation for the effect of gadolinium doping.

preprint2014arXiv

Revisiting Directed Polymers with heavy-tailed disorder

In this mostly numerical study, we revisit the statistical properties of the ground state of a directed polymer in a $d=1+1$ "hilly" disorder landscape, i.e. when the quenched disorder has power-law tails. When disorder is Gaussian, the polymer minimizes its total energy through a collective optimization, where the energy of each visited site only weakly contributes to the total. Conversely, a hilly landscape forces the polymer to distort and explore a larger portion of space to reach some particularly deep energy sites. As soon as the fifth moment of the disorder diverges, this mechanism radically changes the standard "KPZ" scaling behaviour of the directed polymer, and new exponents prevail. After confirming again that the Flory argument accurately predicts these exponent in the tail-dominated phase, we investigate several other statistical features of the ground state that shed light on this unusual transition and on the accuracy of the Flory argument. We underline the theoretical challenge posed by this situation, which paradoxically becomes even more acute above the upper critical dimension.

preprint2014arXiv

Scaling description of the yielding transition in soft amorphous solids at zero temperature

Yield stress materials flow if a sufficiently large shear stress is ap- plied. Although such materials are ubiquitous and relevant for indus- try, there is no accepted microscopic description of how they yield, even in the simplest situations where temperature is negligible and where flow inhomogeneities such as shear bands or fractures are ab- sent. Here we propose a scaling description of the yielding transition in amorphous solids made of soft particles at zero temperature. Our description makes a connection between the Herschel-Bulkley expo- nent characterizing the singularity of the flow curve near the yield stress Σc, the extension and duration of the avalanches of plasticity observed at threshold, and the density P(x) of soft spots, or shear transformation zones, as a function of the stress increment x be- yond which they yield. We argue that the critical exponents of the yielding transition can be expressed in terms of three independent exponents θ, df and z, characterizing respectively the density of soft spots, the fractal dimension of the avalanches, and their duration. Our description shares some similarity with the depinning transition that occurs when an elastic manifold is driven through a random potential, but also presents some striking differences. We test our arguments in an elasto-plastic model, an automaton model similar to those used in depinning, but with a different interaction kernel, and find satisfying agreement with our predictions both in two and three dimensions.

preprint2014arXiv

Strong pinning of propagation fronts in adverse flow

Reaction fronts evolving in a porous medium exhibit a rich dynamical behaviour. In presence of an adverse flow, experiments show that the front slows down and eventually gets pinned, displaying a particular sawtooth shape. Extensive numerical simulations of the hydrodynamic equations confirm the experimental observations. Here we propose a stylized model, predicting two possible outcomes of the experiments for large adverse flow: either the front develops a sawtooth shape, or it acquires a complicated structure with islands and overhangs. A simple criterion allows to distinguish between the two scenarios and its validity is reproduced by direct hydrodynamical simulations. Our model gives a better understanding of the transition and is relevant in a variety of domains, when the pinning regime is strong and only relies on a small number of sites.

preprint2013arXiv

Anomalous fluctuations of currents in Sinai-type random chains with strongly correlated disorder

We study properties of a random walk in a generalized Sinai model, in which a quenched random potential is a trajectory of a fractional Brownian motion with arbitrary Hurst parameter H, 0< H <1, so that the random force field displays strong spatial correlations. In this case, the disorder-average mean-square displacement grows in proportion to log^{2/H}(n), n being time. We prove that moments of arbitrary order k of the steady-state current J_L through a finite segment of length L of such a chain decay as L^{-(1-H)}, independently of k, which suggests that despite a logarithmic confinement the average current is much higher than its Fickian counterpart in homogeneous systems. Our results reveal a paradoxical behavior such that, for fixed n and L, the mean square displacement decreases when one varies H from 0 to 1, while the average current increases. This counter-intuitive behavior is explained via an analysis of representative realizations of disorder.

preprint2013arXiv

Asymmetric Lévy flights in the presence of absorbing boundaries

We consider a one dimensional asymmetric random walk whose jumps are identical, independent and drawn from a distribution ϕ(η) displaying asymmetric power law tails (i.e. ϕ(η) \sim c/η^{α+1} for large positive jumps and ϕ(η) \sim c/(γ|η|^{α+1}) for large negative jumps, with 0 < α< 2). In absence of boundaries and after a large number of steps n, the probability density function (PDF) of the walker position, x_n, converges to an asymmetric Lévy stable law of stability index αand skewness parameter β=(γ-1)/(γ+1). In particular the right tail of this PDF decays as c n/x_n^{1+α}. Much less is known when the walker is confined, or partially confined, in a region of the space. In this paper we first study the case of a walker constrained to move on the positive semi-axis and absorbed once it changes sign. In this case, the persistence exponent θ_+, which characterizes the algebraic large time decay of the survival probability, can be computed exactly and we show that the tail of the PDF of the walker position decays as c \, n/[(1-θ_+) \, x_n^{1+α}]. This last result can be generalized in higher dimensions such as a planar Lévy walker confined in a wedge with absorbing walls. Our results are corroborated by precise numerical simulations.

preprint2013arXiv

Dynamics of a tagged monomer: Effects of elastic pinning and harmonic absorption

We study the dynamics of a tagged monomer of a Rouse polymer for different initial configurations. In the case of free evolution, the monomer displays subdiffusive behavior with strong memory of the initial state. In presence of either elastic pinning or harmonic absorption, we show that the steady state is independent of the initial condition which however strongly affects the transient regime, resulting in non-monotonous behavior and power-law relaxation with varying exponents.

preprint2013arXiv

Effect of disorder geometry on the critical force in disordered elastic systems

We address the effect of disorder geometry on the critical force in disordered elastic systems. We focus on the model system of a long-range elastic line driven in a random landscape. In the collective pinning regime, we compute the critical force perturbatively. Not only our expression for the critical force confirms previous results on its scaling with respect to the microscopic disorder parameters, it also provides its precise dependence on the disorder geometry (represented by the disorder two-point correlation function). Our results are successfully compared to the results of numerical simulations for random field and random bond disorders.

preprint2013arXiv

From microstructural features to effective toughness in disordered brittle solids

The relevant parameters at the microstructure scale that govern the macroscopic toughness of disordered brittle materials are investigated theoretically. We focus on planar crack propagation and describe the front evolution as the propagation of a long-range elastic line within a plane with random distribution of toughness. Our study reveals two regimes: in the collective pinning regime, the macroscopic toughness can be expressed as a function of a few parameters only, namely the average and the standard deviation of the local toughness distribution and the correlation lengths of the heterogeneous toughness field; in the individual pinning regime, the passage from micro to macroscale is more subtle and the full distribution of local toughness is required to be predictive. Beyond the failure of brittle solids, our findings illustrate the complex filtering process of microscale quantities towards the larger scales into play in a broad range of systems governed by the propagation of an elastic interface in a disordered medium.

preprint2013arXiv

On the density of shear transformation zones in amorphous solids

We study the stability of amorphous solids, focusing on the distribution P(x) of the local stress increase x that would lead to an instability. We argue that this distribution is singular P(x)x^θ, where the exponent θ is non-zero if the elastic interaction between rearranging regions is non-monotonic, and increases with the interaction range. For a class of finite dimensional models we show that stability implies a lower bound on θ, which is found to lie near saturation. For quadrupolar interactions these models yield θ ~ 0.6 for d=2 and θ~ 0.4 in d=3 where d is the spatial dimension, accurately capturing previously unresolved observations in atomistic models, both in quasi-static flow and after a fast quench.

preprint2013arXiv

Sampling fractional Brownian motion in presence of absorption: a Markov Chain method

We study fractional Brownian motion (fBm) characterized by the Hurst exponent H. Using a Monte Carlo sampling technique, we are able to numerically generate fBm processes with an absorbing boundary at the origin at discrete times for a large number of 10^7 time steps even for small values like H=1/4. The results are compatible with previous analytical results that the distribution of (rescaled) endpoints y follow a power law P(y) y^ϕwith ϕ=(1-H)/H, even for small values of H. Furthermore, for the case H=0.5 we also study analytically the finite-length corrections to the first order, namely a plateau of P(y) for y->0 which decreases with increasing process length. These corrections are compatible with the numerical results.

preprint2013arXiv

Spatial extent of an outbreak in animal epidemics

Characterizing the spatial extent of epidemics at the outbreak stage is key to controlling the evolution of the disease. At the outbreak, the number of infected individuals is typically small, so that fluctuations around their average are important: then, it is commonly assumed that the susceptible-infected-recovered (SIR) mechanism can be described by a stochastic birth-death process of Galton-Watson type. The displacements of the infected individuals can be modelled by resorting to Brownian motion, which is applicable when long-range movements and complex network interactions can be safely neglected, as in case of animal epidemics. In this context, the spatial extent of an epidemic can be assessed by computing the convex hull enclosing the infected individuals at a given time. We derive the exact evolution equations for the mean perimeter and the mean area of the convex hull, and compare them with Monte Carlo simulations.

preprint2012arXiv

Counting function fluctuations and extreme value threshold in multifractal patterns: the case study of an ideal $1/f$ noise

To understand the sample-to-sample fluctuations in disorder-generated multifractal patterns we investigate analytically as well as numerically the statistics of high values of the simplest model - the ideal periodic $1/f$ Gaussian noise. By employing the thermodynamic formalism we predict the characteristic scale and the precise scaling form of the distribution of number of points above a given level. We demonstrate that the powerlaw forward tail of the probability density, with exponent controlled by the level, results in an important difference between the mean and the typical values of the counting function. This can be further used to determine the typical threshold $x_m$ of extreme values in the pattern which turns out to be given by $x_m^{(typ)}=2-c\ln{\ln{M}}/\ln{M}$ with $c=3/2$. Such observation provides a rather compelling explanation of the mechanism behind universality of $c$. Revealed mechanisms are conjectured to retain their qualitative validity for a broad class of disorder-generated multifractal fields. In particular, we predict that the typical value of the maximum $p_{max}$ of intensity is to be given by $-\ln{p_{max}} = α_{-}\ln{M} + \frac{3}{2f'(α_{-})}\ln{\ln{M}} + O(1)$, where $f(α)$ is the corresponding singularity spectrum vanishing at $α=α_{-}>0$. For the $1/f$ noise we also derive exact as well as well-controlled approximate formulas for the mean and the variance of the counting function without recourse to the thermodynamic formalism.

preprint2012arXiv

Lévy flights on the half line

We study the probability distribution function (pdf) of the position of a Lévy flight of index 0<α<2 in presence of an absorbing wall at the origin. The solution of the associated fractional Fokker-Planck equation can be constructed using a perturbation scheme around the Brownian solution (corresponding to α= 2), as an expansion in ε= 2 - α. We obtain an explicit analytical solution, exact at the first order in ε, which allows us to conjecture the precise asymptotic behavior of this pdf, including the first subleading corrections, for any α. Careful numerical simulations, as well as an exact computation for α= 1, confirm our conjecture.

preprint2012arXiv

Short time growth of a KPZ interface with flat initial conditions

The short time behavior of the 1+1 dimensional KPZ growth equation with a flat initial condition is obtained from the exact expressions of the moments of the partition function of a directed polymer with one endpoint free and the other fixed. From these expressions, the short time expansions of the lowest cumulants of the KPZ height field are exactly derived. The results for these two classes of cumulants are checked in high precision lattice numerical simulations. The short time limit considered here is relevant for the study of the interface growth in the large diffusivity/weak noise limit, and describes the universal crossover between the Edwards-Wilkinson and KPZ universality classes for an initially flat interface.

preprint2012arXiv

Tuning spreading and avalanche-size exponents in directed percolation with modified activation probabilities

We consider the directed percolation process as a prototype of systems displaying a nonequilibrium phase transition into an absorbing state. The model is in a critical state when the activation probability is adjusted at some precise value p_c. Criticality is lost as soon as the probability to activate sites at the first attempt, p1, is changed. We show here that criticality can be restored by "compensating" the change in p1 by an appropriate change of the second time activation probability p2 in the opposite direction. At compensation, we observe that the bulk exponents of the process coincide with those of the normal directed percolation process. However, the spreading exponents are changed, and take values that depend continuously on the pair (p1, p2). We interpret this situation by acknowledging that the model with modified initial probabilities has an infinite number of absorbing states.

preprint2011arXiv

Perturbation Theory for Fractional Brownian Motion in Presence of Absorbing Boundaries

Fractional Brownian motion is a Gaussian process x(t) with zero mean and two-time correlations <x(t)x(s)> ~ t^{2H} + s^{2H} - |t-s|^{2H}, where H, with 0<H<1 is called the Hurst exponent. For H = 1/2, x(t) is a Brownian motion, while for H unequal 1/2, x(t) is a non-Markovian process. Here we study x(t) in presence of an absorbing boundary at the origin and focus on the probability density P(x,t) for the process to arrive at x at time t, starting near the origin at time 0, given that it has never crossed the origin. It has a scaling form P(x,t) ~ R(x/t^H)/t^H. Our objective is to compute the scaling function R(y), which up to now was only known for the Markov case H=1/2. We develop a systematic perturbation theory around this limit, setting H = 1/2 + epsilon, to calculate the scaling function R(y) to first order in epsilon. We find that R(y) behaves as R(y) ~ y^phi as y -> 0 (near the absorbing boundary), while R(y) ~ y^gamma exp(-y^2/2) as y -> oo, with phi = 1 - 4 epsilon + O(epsilon^2) and gamma = 1 - 2 epsilon + O(epsilon^2). Our epsilon-expansion result confirms the scaling relation phi = (1-H)/H proposed in Ref. [28]. We verify our findings via numerical simulations for H = 2/3. The tools developed here are versatile, powerful, and adaptable to different situations.

preprint2011arXiv

Shock statistics in higher-dimensional Burgers turbulence

We conjecture the exact shock statistics in the inviscid decaying Burgers equation in D>1 dimensions, with a special class of correlated initial velocities, which reduce to Brownian for D=1. The prediction is based on a field-theory argument, and receives support from our numerical calculations. We find that, along any given direction, shocks sizes and locations are uncorrelated.

preprint2011arXiv

The convex hull for a random acceleration process in two dimensions

We compute exactly the mean perimeter <L(T)> and the mean area <A(T)> of the convex hull of a random acceleration process of duration T in two dimensions. We use an exact mapping that relates, via Cauchy's formulae, the computation of the perimeter and the area of the convex hull of an arbitrary two dimensional stochastic process [x(t); y(t)] to the computation of the extreme value statistics of the associated one dimensional component process x(t). The latter can be computed exactly for the one dimensional random acceleration process even though the process in non-Markovian. Physically, our results are relevant in describing theaverage shape of a semi-flexible ideal polymer chain in two dimensions.

preprint2010arXiv

Free-energy distribution of the directed polymer at high temperature

We study the directed polymer of length $t$ in a random potential with fixed endpoints in dimension 1+1 in the continuum and on the square lattice, by analytical and numerical methods. The universal regime of high temperature $T$ is described, upon scaling 'time' $t \sim T^5/κ$ and space $x = T^3/κ$ (with $κ=T$ for the discrete model) by a continuum model with $δ$-function disorder correlation. Using the Bethe Ansatz solution for the attractive boson problem, we obtain all positive integer moments of the partition function. The lowest cumulants of the free energy are predicted at small time and found in agreement with numerics. We then obtain the exact expression at any time for the generating function of the free energy distribution, in terms of a Fredholm determinant. At large time we find that it crosses over to the Tracy Widom distribution (TW) which describes the fixed $T$ infinite $t$ limit. The exact free energy distribution is obtained for any time and compared with very recent results on growth and exclusion models.

preprint2010arXiv

Freezing Transition in Decaying Burgers Turbulence and Random Matrix Dualities

We reveal a phase transition with decreasing viscosity $ν$ at ν=ν_c>0 in one-dimensional decaying Burgers turbulence with a power-law correlated random profile of Gaussian-distributed initial velocities <v(x,0)v(x',0)>\sim|x-x'|^{-2}. The low-viscosity phase exhibits non-Gaussian one-point probability density of velocities, continuously dependent on ν, reflecting a spontaneous one step replica symmetry breaking (RSB) in the associated statistical mechanics problem. We obtain the low orders cumulants analytically. Our results, which are checked numerically, are based on combining insights in the mechanism of the freezing transition in random logarithmic potentials with an extension of duality relations discovered recently in Random Matrix Theory. They are essentially non mean-field in nature as also demonstrated by the shock size distribution computed numerically and different from the short range correlated Kida model, itself well described by a mean field one step RSB ansatz. We also provide some insights for the finite viscosity behaviour of velocities in the latter model.

preprint2010arXiv

Maximum of N Independent Brownian Walkers till the First Exit From the Half Space

We consider the one-dimensional target search process that involves an immobile target located at the origin and $N$ searchers performing independent Brownian motions starting at the initial positions $\vec x = (x_1,x_2,..., x_N)$ all on the positive half space. The process stops when the target is first found by one of the searchers. We compute the probability distribution of the maximum distance $m$ visited by the searchers till the stopping time and show that it has a power law tail: $P_N(m|\vec x)\sim B_N (x_1x_2... x_N)/m^{N+1}$ for large $m$. Thus all moments of $m$ up to the order $(N-1)$ are finite, while the higher moments diverge. The prefactor $B_N$ increases with $N$ faster than exponentially. Our solution gives the exit probability of a set of $N$ particles from a box $[0,L]$ through the left boundary. Incidentally, it also provides an exact solution of the Laplace's equation in an $N$-dimensional hypercube with some prescribed boundary conditions. The analytical results are in excellent agreement with Monte Carlo simulations.

preprint2010arXiv

The universal high temperature regime of pinned elastic objects

We study the high temperature regime within the glass phase of an elastic object with short ranged disorder. In that regime we argue that the scaling functions of any observable are described by a continuum model with a $δ$-correlated disorder and that they are universal up to only two parameters that can be explicitly computed. This is shown numerically on the roughness of directed polymer models and by dimensional and functional renormalization group arguments. We discuss experimental consequences such as non-monotonous behaviour with temperature.

preprint2009arXiv

Avalanche-size distribution at the depinning transition: A numerical test of the theory

We calculate numerically the sizes S of jumps (avalanches) between successively pinned configurations of an elastic line (d=1) or interface (d=2), pulled by a spring of (small) strength m^2 in a random-field landscape. We obtain strong evidence that the size distribution, away from the small-scale cutoff, takes the form P(S) = p(S/S_m) <S>/S_m^2, where S_m:=<S^2>/(2<S>), proportional to m^(-d-zeta), is the scale of avalanches, and zeta the roughness exponent at the depinning transition. Measurement of the scaling function f(s) := s^tau p(s) is compared with the predictions from a recent Functional RG (FRG) calculation, both at mean-field and one-loop level. The avalanche-size exponent tau is found in good agreement with the conjecture tau = 2- 2/(d+zeta), recently confirmed to one loop via the FRG. The function f(s) exhibits a shoulder and a stretched exponential decay at large s, with ln f(s) proportional to - s^delta, and delta approximately 7/6 in d=1. The function f(s), universal ratios of moments, and the generating function <exp(lambda s)> are found in excellent agreement with the one-loop FRG predictions. The distribution of local avalanche sizes, i.e. of the jumps of a subspace of the manifold of dimension d', is also computed and compared to our FRG predictions, and to the conjecture tau' =2- 2/(d'+zeta).

preprint2009arXiv

Hitting probability for anomalous diffusion processes

We present the universal features of the hitting probability $Q(x,L)$, the probability that a generic stochastic process starting at $x$ and evolving in a box $[0,L]$ hits the upper boundary $L$ before hitting the lower boundary at 0. For a generic self-affine process (describing, for instance, the polymer translocation through a nanopore) we show that $Q(x,L)=Q(x/L)$ and the scaling function $Q(z)\sim z^ϕ$ as $z\to 0$ with $ϕ=θ/H$ where $H$ and $θ$ are respectively the Hurst exponent and the persistence exponent of the process. This result is verified in several exact calculations including when the process represents the position of a particle diffusing in a disordered potential. We also provide numerical supports for our analytical results.

preprint2009arXiv

Statistical Mechanics of Logarithmic REM: Duality, Freezing and Extreme Value Statistics of $1/f$ Noises generated by Gaussian Free Fields

We compute the distribution of the partition functions for a class of one-dimensional Random Energy Models (REM) with logarithmically correlated random potential, above and at the glass transition temperature. The random potential sequences represent various versions of the 1/f noise generated by sampling the two-dimensional Gaussian Free Field (2dGFF) along various planar curves. Our method extends the recent analysis of Fyodorov Bouchaud from the circular case to an interval and is based on an analytical continuation of the Selberg integral. In particular, we unveil a {\it duality relation} satisfied by the suitable generating function of free energy cumulants in the high-temperature phase. It reinforces the freezing scenario hypothesis for that generating function, from which we derive the distribution of extrema for the 2dGFF on the $[0,1]$ interval. We provide numerical checks of the circular and the interval case and discuss universality and various extensions. Relevance to the distribution of length of a segment in Liouville quantum gravity is noted.

preprint2009arXiv

The longest excursion of fractional Brownian motion : numerical evidence of non-Markovian effects

We study, using exact numerical simulations, the statistics of the longest excursion l_{\max}(t) up to time t for the fractional Brownian motion with Hurst exponent 0<H<1. We show that in the large t limit, < l_{\max}(t) > \propto Q_\infty t where Q_\infty \equiv Q_\infty(H) depends continuously on H, and in a non trivial way. These results are compared with exact analytical results obtained recently for a renewal process with an associated persistence exponent θ= 1-H. This comparison shows that Q_\infty(H) carries the clear signature of non-Markovian effects for H\neq 1/2. The pre-asymptotic behavior of < l_{\max}(t)> is also discussed.

preprint2005arXiv

Geometry of Gaussian signals

We consider Gaussian signals, i.e. random functions $u(t)$ ($t/L \in [0,1]$) with independent Gaussian Fourier modes of variance $\sim 1/q^α$, and compute their statistical properties in small windows $[x, x+δ]$. We determine moments of the probability distribution of the mean square width of $u(t)$ in powers of the window size $δ$. We show that the moments, in the small-window limit $δ\ll 1$, become universal, whereas they strongly depend on the boundary conditions of $u(t)$ for larger $δ$. For $α> 3$, the probability distribution is computed in the small-window limit and shown to be independent of $α$.

preprint2004arXiv

X-ray spectrum of a pinned charge density wave

We calculate the X-ray diffraction spectrum produced by a pinned charge density wave (CDW). The signature of the presence of a CDW consists of two satellite peaks, asymmetric as a consequence of disorder. The shape and the intensity of these peaks are determined in the case of a collective weak pinning using the variational method. We predict divergent asymmetric peaks, revealing the presence of a Bragg glass phase. We deal also with the long range Coulomb interactions, concluding that both peak divergence and anisotropy are enhanced. Finally we discuss how to detect experimentally the Bragg glass phase in the view of the role played by the finite resolution of measurements.