BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Cecile Monthus

Cecile Monthus appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
cond-mat.dis-nncond-mat.stat-mechmath-phmath.MPcond-matmath.PRcond-mat.str-el

Catalog footprint

What is connected

52works
7topics
2close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

52 published item(s)

preprint2022arXiv

Anomalous dynamical large deviations of local empirical densities and activities in the pure and in the random kinetically-constrained East Model

The East model is the simplest one-dimensional kinetically-constrained model of $N$ spins with a trivial equilibrium that displays anomalously large spatio-temporal fluctuations, with characteristic "space-time bubbles" in trajectory space, and with a discontinuity at the origin for the first derivative of the scaled cumulant generating function of the total activity. These striking dynamical properties are revisited via the large deviations at various levels for the relevant local empirical densities and activities that only involve two consecutive spins. This framework allows to characterize their anomalous rate functions and to analyze the consequences for all the time-additive observables that can be reconstructed from them, both for the pure and for the random East model. These singularities in dynamical large deviations properties disappear when the hard-constraint of the East model is replaced by the soft constraint.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Markov trajectories : Microcanonical Ensembles based on empirical observables as compared to Canonical Ensembles based on Markov generators

The Ensemble of trajectories $x(0 \leq t \leq T)$ produced by the Markov generator $M$ can be considered as 'Canonical' for the following reasons : (C1) the probability of the trajectory $x(0 \leq t \leq T)$ can be rewritten as the exponential of a linear combination of its relevant empirical time-averaged observables $E_n$, where the coefficients involving the Markov generator are their fixed conjugate parameters; (C2) the large deviations properties of these empirical observables $E_n$ for large $T$ are governed by the explicit rate function $I^{[2.5]}_M (E_.) $ at Level 2.5, while in the thermodynamic limit $T=+\infty$, they concentrate on their typical values $E_n^{typ[M]}$ determined by the Markov generator $M$. This concentration property in the thermodynamic limit $T=+\infty$ suggests to introduce the notion of the 'Microcanonical Ensemble' at Level 2.5 for stochastic trajectories $x(0 \leq t \leq T)$, where all the relevant empirical variables $E_n$ are fixed to some values $E^*_n$ and cannot fluctuate anymore for finite $T$. The goal of the present paper is to discuss its main properties : (MC1) when the long trajectory $x(0 \leq t \leq T) $ belongs the Microcanonical Ensemble with the fixed empirical observables $E_n^*$, the statistics of its subtrajectory $x(0 \leq t \leq τ) $ for $1 \ll τ\ll T $ is governed by the Canonical Ensemble associated to the Markov generator $M^*$ that would make the empirical observables $E_n^*$ typical ; (MC2) in the Microcanonical Ensemble, the central role is played by the number $Ω^{[2.5]}_T(E^*_.) $ of stochastic trajectories of duration $T$ with the given empirical observables $E^*_n$, and by the corresponding explicit Boltzmann entropy $S^{[2.5]}( E^*_. ) = [\ln Ω^{[2.5]}_T(E^*_.)]/T $. This general framework is applied to continuous-time Markov Jump processes and to discrete-time Markov chains with illustrative examples.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Microcanonical conditioning of Markov processes on time-additive observables

The recent study by B. De Bruyne, S. N. Majumdar, H. Orland and G. Schehr [arXiv:2110.07573], concerning the conditioning of the Brownian motion and of random walks on global dynamical constraints over a finite time-window $T$, is reformulated as a general framework for the 'microcanonical conditioning' of Markov processes on time-additive observables. This formalism is applied to various types of Markov processes, namely discrete-time Markov chains, continuous-time Markov jump processes and diffusion processes in arbitrary dimension. In each setting, the time-additive observable is also fully general, i.e. it can involve both the time spent in each configuration and the elementary increments of the Markov process. The various cases are illustrated via simple explicit examples. Finally, we describe the link with the 'canonical conditioning' based on the generating function of the time-additive observable for finite time $T$, while the regime of large time $T$ allows to recover the standard large deviation analysis of time-additive observables via the deformed Markov operator approach.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Inhomogeneous asymmetric exclusion processes between two reservoirs : large deviations for the local empirical observables in the Mean-Field approximation

For a given inhomogeneous exclusion processes on $N$ sites between two reservoirs, the trajectories probabilities allow to identify the relevant local empirical observables and to obtain the corresponding rate function at Level 2.5. In order to close the hierarchy of the empirical dynamics that appear in the stationarity constraints, we consider the simplest approximation, namely the Mean-Field approximation for the empirical density of two consecutive sites, in direct correspondence with the previously studied Mean-Field approximation for the steady state. For a given inhomogeneous Totally Asymmetric model (TASEP), this Mean-Field approximation yields the large deviations for the joint distribution of the empirical density profile and of the empirical current around the mean-field steady state; the further explicit contraction over the current allows to obtain the large deviations of the empirical density profile alone. For a given inhomogeneous Asymmetric model (ASEP), the local empirical observables also involve the empirical activities of the links and of the reservoirs; the further explicit contraction over these activities yields the large deviations for the joint distribution of the empirical density profile and of the empirical current. The consequences for the large deviations properties of time-additive space-local observables are also discussed in both cases.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Large deviations for metastable states of Markov processes with absorbing states with applications to population models in stable or randomly switching environment

The large deviations at Level 2.5 are applied to Markov processes with absorbing states in order to obtain the explicit extinction rate of metastable quasi-stationary states in terms of their empirical time-averaged density and of their time-averaged empirical flows over a large time-window $T$. The standard spectral problem for the slowest relaxation mode can be recovered from the full optimization of the extinction rate over all these empirical observables and the equivalence can be understood via the Doob generator of the process conditioned to survive up to time $T$. The large deviation properties of any time-additive observable of the Markov trajectory before extinction can be derived from the Level 2.5 via the decomposition of the time-additive observable in terms of the empirical density and the empirical flows. This general formalism is described for continuous-time Markov chains, with applications to population birth-death model in a stable or in a switching environment, and for diffusion processes in dimension $d$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Finite size scaling for the Many-Body-Localization Transition : finite-size-pseudo-critical points of individual eigenstates

To understand the finite-size-scaling properties of phases transitions in classical and quantum models in the presence of quenched disorder, it has proven to be fruitful to introduce the notion of a finite-size-pseudo-critical point in each disordered sample and to analyze its sample-to-sample fluctuations as a function of the size. For the Many-Body-Localization transition, where very strong eigenstate-to-eigenstate fluctuations have been numerically reported even within a given disordered sample at a given energy density [X. Yu, D. J. Luitz, B. K. Clark, arxiv:1606.01260 and V. Khemani, S. P. Lim, D. N. Sheng, D. A. Huse,arxiv:1607.05756], it seems thus useful to introduce the notion of a finite-size-pseudo-critical point for each individual eigenstate and to study its eigenstate-to-eigenstate fluctuations governed by the correlation length exponent $ν$. The scaling properties of critical eigenstates are also expected to appear much more clearly if one considers each eigenstate at its finite-size-pseudo-critical point, where it is 'truly critical', while standard averages over eigenstates and samples in the critical region actually see a mixture of states that are effectively either localized or delocalized.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Flow towards diagonalization for Many-Body-Localization models : adaptation of the Toda matrix differential flow to random quantum spin chains

The iterative methods to diagonalize matrices and many-body Hamiltonians can be reformulated as flows of Hamiltonians towards diagonalization driven by unitary transformations that preserve the spectrum. After a comparative overview of the various types of discrete flows (Jacobi, QR-algorithm) and differential flows (Toda, Wegner, White) that have been introduced in the past, we focus on the random XXZ chain with random fields in order to determine the best closed flow within a given subspace of running Hamiltonians. For the special case of the free-fermion random XX chain with random fields, the flow coincides with the Toda differential flow for tridiagonal matrices which is related to the classical integrable Toda chain and which can be seen as the continuous analog of the discrete QR-algorithm. For the random XXZ chain with random fields that displays a Many-Body-Localization transition, the present differential flow should be an interesting alternative to compare with the discrete flow that has been proposed recently to study the Many-Body-Localization properties in a model of interacting fermions (L. Rademaker and M. Ortuno, Phys. Rev. Lett. 116, 010404 (2016)).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Level repulsion exponent $β$ for Many-Body Localization Transitions and for Anderson Localization Transitions via Dyson Brownian Motion

The generalization of the Dyson Brownian Motion approach of random matrices to Anderson Localization (AL) models [Chalker, Lerner and Smith PRL 77, 554 (1996)] and to Many-Body Localization (MBL) Hamiltonians [Serbyn and Moore arxiv:1508.07293] is revisited to extract the level repulsion exponent $β$, where $β=1$ in the delocalized phase governed by the Wigner-Dyson statistics, $β=0$ in the localized phase governed by the Poisson statistics, and $0<β_c<1$ at the critical point. The idea is that the Gaussian disorder variables $h_i$ are promoted to Gaussian stationary processes $h_i(t)$ in order to sample the disorder stationary distribution with some time correlation $τ$. The statistics of energy levels can be then studied via Langevin and Fokker-Planck equations. For the MBL quantum spin Hamiltonian with random fields $h_i$, we obtain $β=2q^{EA}_{n,n+1}(N)/q^{EA}_{n,n}(N) $ in terms of the Edwards-Anderson matrix $q^{EA}_{nm}(N) \equiv \frac{1}{N} \sum_{i=1}^N | < ϕ_n | σ_i^z | ϕ_m> |^2 $ for the same eigenstate $m=n$ and for consecutive eigenstates $m=n+1$. For the Anderson Localization tight-binding Hamiltonian with random on-site energies $h_i$, we find $β=2 Y_{n,n+1}(N)/(Y_{n,n}(N)-Y_{n,n+1}(N)) $ in terms of the Density Correlation matrix $Y_{nm}(N) \equiv \sum_{i=1}^N | < ϕ_n | i> |^2 | <i | ϕ_m> |^2 $ for consecutive eigenstates $m=n+1$, while the diagonal element $m=n$ corresponds to the Inverse Participation Ratio $Y_{nn}(N) \equiv \sum_{i=1}^N | < ϕ_n | i> |^4 $ of the eigenstate $| ϕ_n>$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Localization transition in random Levy matrices : multifractality of eigenvectors in the localized phase and at criticality

For random Lévy matrices of size $N \times N$, where matrix elements are drawn with some heavy-tailed distribution $P(H_{ij}) \propto N^{-1} |H_{ij} |^{-1-μ}$ with $0<μ<2$ (infinite variance), there exists an extensive number of finite eigenvalues $E=O(1)$, while the maximal eigenvalue grows as $E_{max} \sim N^{\frac{1}μ}$. Here we study the localization properties of the corresponding eigenvectors via some strong disorder perturbative expansion that remains consistent within the localized phase and that yields their Inverse Participation Ratios (I.P.R.) $Y_q$ as a function of the continuous parameter $0<q<+\infty$. In the region $0<μ<1$, we find that all eigenvectors are localized but display some multifractality : the IPR are finite above some threshold $q>q_c$ but diverge in the region $0<q<q_c$ near the origin. In the region $1<μ<2$, only the sub-extensive fraction $N^{\frac{3}{2+μ}}$ of the biggest eigenvalues corresponding to the region $| E| \geq N^{\frac{(μ-1)}{μ(2+μ)}} $ remains localized, while the extensive number of other states of smaller energy are delocalized. For the extensive number of finite eigenvalues $E=O(1)$, the localization/delocalization transition thus takes place at the critical value $μ_c=1$ corresponding to Cauchy matrices : the Inverse Participation Ratios $Y_q$ of the corresponding critical eigenstates follow the Strong-Multifractality Spectrum characterized by the generalized fractal dimensions $D^{criti}(q) = \frac{1-2q}{1-q} θ(0 \leq q \leq \frac{1}{2})$, which has been found previously in various other Localization problems in spaces of effective infinite dimensionality.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Many Body Localization Transition in the strong disorder limit : entanglement entropy from the statistics of rare extensive resonances

The space of one-dimensional disordered interacting quantum models displaying a Many-Body-Localization Transition seems sufficiently rich to produce critical points with level statistics interpolating continuously between the Poisson statistics of the Localized phase and the Wigner-Dyson statistics of the Delocalized Phase. In this paper, we consider the strong disorder limit of the MBL transition, where the critical level statistics is close to the Poisson statistics. We analyse a one-dimensional quantum spin model, in order to determine the statistical properties of the rare extensive resonances that are needed to destabilize the MBL phase. At criticality, we find that the entanglement entropy can grow with an exponent $0<α< 1$ anywhere between the area law $α=0$ and the volume law $α=1$, as a function of the resonances properties, while the entanglement spectrum follows the strong multifractality statistics. In the MBL phase near criticality, we obtain the simple value $ν=1$ for the correlation length exponent. Independently of the strong disorder limit, we explain why for the Many-Body-Localization transition concerning individual eigenstates, the correlation length exponent $ν$ is not constrained by the usual Harris inequality $ν\geq 2/d$, so that there is no theoretical inconsistency with the best numerical measure $ν= 0.8 (3)$ obtained by D. J. Luitz, N. Laflorencie and F. Alet, Phys. Rev. B 91, 081103 (2015).

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Many-Body Localization : construction of the emergent local conserved operators via block real-space renormalization

A Fully Many-Body Localized (FMBL) quantum disordered system is characterized by the emergence of an extensive number of local conserved operators that prevents the relaxation towards thermal equilibrium. These local conserved operators can be seen as the building blocks of the whole set of eigenstates. In this paper, we propose to construct them explicitly via some block real-space renormalization. The principle is that each RG step diagonalizes the smallest remaining blocks and produces a conserved operator for each block. The final output for a chain of $N$ spins is a hierarchical organization of the $N$ conserved operators with $\left(\frac{\ln N}{\ln 2}\right)$ layers. The system-size nature of the conserved operators of the top layers is necessary to describe the possible long-ranged order of the excited eigenstates and the possible critical points between different FMBL phases. We discuss the similarities and the differences with the Strong Disorder RSRG-X method that generates the whole set of the $2^N$ eigenstates via a binary tree of $N$ layers. The approach is applied to the Long-Ranged Quantum Spin-Glass Ising model, where the constructed excited eigenstates are found to be exactly like ground states in another disorder realization, so that they can be either in the paramagnetic phase, in the spin-glass phase or critical.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Many-Body-Localization Transition : strong multifractality spectrum for matrix elements of local operators

For short-ranged disordered quantum models in one dimension, the Many-Body-Localization is analyzed via the adaptation to the Many-Body context [M. Serbyn, Z. Papic and D.A. Abanin, PRX 5, 041047 (2015)] of the Thouless point of view on the Anderson transition : the question is whether a local interaction between two long chains is able to reshuffle completely the eigenstates (Delocalized phase with a volume-law entanglement) or whether the hybridization between tensor states remains limited (Many-Body-Localized Phase with an area-law entanglement). The central object is thus the level of Hybridization induced by the matrix elements of local operators, as compared with the difference of diagonal energies. The multifractal analysis of these matrix elements of local operators is used to analyze the corresponding statistics of resonances. Our main conclusion is that the critical point is characterized by the Strong-Multifractality Spectrum $f(0 \leq α\leq 2)=\fracα{2}$, well known in the context of Anderson Localization in spaces of effective infinite dimensionality, where the size of the Hilbert space grows exponentially with the volume. Finally, the possibility of a delocalized non-ergodic phase near criticality is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Real-space renormalization for the finite temperature statics and dynamics of the Dyson Long-Ranged Ferromagnetic and Spin-Glass models

The finite temperature dynamics of the Dyson hierarchical classical spins models is studied via real-space renormalization rules concerning the couplings and the relaxation times. For the ferromagnetic model involving Long-Ranged coupling $J(r) \propto r^{-1-σ}$ in the region $1/2<σ<1$ where there exists a non-mean-field-like thermal Ferromagnetic-Paramagnetic transition, the RG flows are explicitly solved: the characteristic relaxation time $τ(L)$ follows the critical power-law $τ(L)\propto L^{z_c(σ)} $ at the phase transition and the activated law $\ln τ(L)\propto L^ψ $ with $ψ=1-σ$ in the ferromagnetic phase. For the Spin-Glass model involving random Long-Ranged couplings of variance $\overline{J^2(r)} \propto r^{-2σ}$ in the region $2/3<σ<1$ where there exists a non-mean-field-like thermal SpinGlass-Paramagnetic transition, the coupled RG flows of the couplings and of the relaxation times are studied numerically : the relaxation time $τ(L)$ follows some power-law $τ(L)\propto L^{z_c(σ)} $ at criticality and the activated law $\ln τ(L)\propto L^ψ $ in the Spin-Glass phase with the dynamical exponent $ψ=1-σ=θ$ coinciding with the droplet exponent governing the flow of the couplings $J(L) \propto L^θ $.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Block Renormalization for quantum Ising models in dimension $d=2$ : applications to the pure and random ferromagnet, and to the spin-glass

For the quantum Ising chain, the self-dual block renormalization procedure of Fernandez-Pacheco [Phys. Rev. D 19, 3173 (1979)] is known to reproduce exactly the location of the zero-temperature critical point and the correlation length exponent $ν=1$. Recently, Miyazaki and Nishimori [Phys. Rev. E 87, 032154 (2013)] have proposed to study the disordered quantum Ising model in dimensions $d>1$ by applying the Fernandez-Pacheco procedure successively in each direction. To avoid the inequivalence of directions of their approach, we propose here an alternative procedure where the $d$ directions are treated on the same footing. For the pure model, this leads to the correlation length exponents $ν\simeq 0.625$ in $d=2$ (to be compared with the 3D classical Ising model exponent $ν\simeq 0.63$) and $ν\simeq 0.5018$ (to be compared with the 4D classical Ising model mean-field exponent $ν=1/2$). For the disordered model in dimension $d=2$, either ferromagnetic or spin-glass, the numerical application of the renormalization rules to samples of linear size $L=4096$ yields that the transition is governed by an Infinite Disorder Fixed Point, with the activated exponent $ψ\simeq 0.65$, the typical correlation exponent $ν_{typ} \simeq 0.44$ and the finite-size correlation exponent $ν_{FS} \simeq 1.25$. We discuss the similarities and differences with the Strong Disorder Renormalization results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dyson Hierarchical Long-Ranged Quantum Spin-Glass via real-space renormalization

We consider the Dyson hierarchical version of the quantum Spin-Glass with random Gaussian couplings characterized by the power-law decaying variance $\overline{J^2(r)} \propto r^{-2σ}$ and a uniform transverse field $h$. The ground state is studied via real-space renormalization to characterize the spinglass-paramagnetic zero temperature quantum phase transition as a function of the control parameter $h$. In the spinglass phase $h<h_c$, the typical renormalized coupling grows with the length scale $L$ as the power-law $J_L^{typ}(h) \propto Υ(h) L^θ$ with the classical droplet exponent $θ=1-σ$, where the stiffness modulus vanishes at criticality $Υ(h) \propto (h_c-h)^μ $, whereas the typical renormalized transverse field decays exponentially $ h^{typ}_L(h) \propto e^{- \frac{L}ξ}$ where the correlation length diverges at the transition $ξ\propto (h_c-h)^{-ν}$. At the critical point $h=h_c$, the typical renormalized coupling $J_L^{typ}(h_c) $ and the typical renormalized transverse field $ h^{typ}_L(h_c)$ display the same power-law behavior $L^{-z}$ with a finite dynamical exponent $z$. The RG rules are applied numerically to chains containing $L=2^{12}=4096 $ spins in order to measure these critical exponents for various values of $σ$ in the region $1/2<σ<1$.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Dyson hierarchical quantum ferromagnetic Ising chain with pure or random transverse fields

The Dyson hierarchical version of the quantum Ising chain with Long-Ranged power-law ferromagnetic couplings $J(r) \propto r^{-1-σ}$ and pure or random transverse fields is studied via real-space renormalization. For the pure case, the critical exponents are explicitly obtained as a function of the parameter $σ$, and are compared with previous results of other approaches. For the random case, the RG rules are numerically applied and the critical behaviors are compared with previous Strong Disorder Renormalization results.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Pure and Random Quantum Ising Chain : Shannon and Renyi entropies of the ground state via real space renormalization

The Shannon and the Renyi entropies of the ground state wavefunction in the pure and in the random quantum Ising chain are studied via the self-dual Fernandez-Pacheco real-space renormalization procedure. In particular, we analyze the critical behavior of the leading extensive term at the quantum phase transition : the derivative with respect to the control parameter is found to be logarithmically divergent in the pure case, and to display a cusp singularity in the random case. This cusp singularity for the random case is also derived via the Strong Disorder Renormalization approach.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Star junctions and watermelons of pure or random quantum Ising chains : finite-size properties of the energy gap at criticality

We consider $M \geq 2$ pure or random quantum Ising chains of $N$ spins when they are coupled via a single star junction at their origins or when they are coupled via two star junctions at the their two ends leading to the watermelon geometry. The energy gap is studied via a sequential self-dual real-space renormalization procedure that can be explicitly solved in terms of Kesten variables containing the initial couplings and and the initial transverse fields. In the pure case at criticality, the gap is found to decay as a power-law $Δ_M \propto N^{-z(M)} $ with the dynamical exponent $z(M)=\frac{M}{2}$ for the single star junction (the case $M=2$ corresponds to $z=1$ for a single chain with free boundary conditions) and $z(M)=M-1$ for the watermelon (the case $M=2$ corresponds to $z=1$ for a single chain with periodic boundary conditions). In the random case at criticality, the gap follows the Infinite Disorder Fixed Point scaling $\ln Δ_M = -N^ψ g$ with the same activated exponent $ψ=\frac{1}{2}$ as the single chain corresponding to $M=2$, and where $g$ is an $O(1)$ random positive variable, whose distribution depends upon the number $M$ of chains and upon the geometry (star or watermelon).

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Fractal dimension of spin glasses interfaces in dimensions $d=2$ and $d=3$ via strong disorder renormalization at zero temperature

For Gaussian Spin Glasses in low dimensions, we introduce a simple Strong Disorder renormalization procedure at zero temperature. In each disordered sample, the difference between the ground states associated to Periodic and Anti-Periodic boundary conditions defines a system-size Domain Wall. The numerical study in dimensions $d=2$ (up to sizes $2048^2$) and $d=3$ (up to sizes $128^3$) yields fractal Domain Walls of dimensions $d_s(d=2) \simeq 1.27$ and $d_s(d=3) \simeq 2.55$ respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Low-temperature dynamics of Long-Ranged Spin-Glasses : full hierarchy of relaxation times via real-space renormalization

We consider the long-ranged Ising spin-glass with random couplings decaying as a power-law of the distance, in the region of parameters where the spin-glass phase exists with a positive droplet exponent. For the Metropolis single-spin-flip dynamics near zero temperature, we construct via real-space renormalization the full hierarchy of relaxation times of the master equation for any given realization of the random couplings. We then analyze the probability distribution of dynamical barriers as a function of the spatial scale. This real-space renormalization procedure represents a simple explicit example of the droplet scaling theory, where the convergence towards local equilibrium on larger and larger scales is governed by a strong hierarchy of activated dynamical processes, with valleys within valleys.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

One-dimensional Ising spin-glass with power-law interaction : real-space renormalization at zero temperature

For the one-dimensional long-ranged Ising spin-glass with random couplings decaying with the distance $r$ as $J(r) \sim r^{-σ}$ and distributed with the Lévy symmetric stable distribution of index $1 <μ\leq 2$ (including the usual Gaussian case $μ=2$), we consider the region $σ>1/μ$ where the energy is extensive. We study two real space renormalization procedures at zero temperature, namely a simple box decimation that leads to explicit calculations, and a strong disorder decimation that can be studied numerically on large sizes. The droplet exponent governing the scaling of the renormalized couplings $J_L \propto L^{θ_μ(σ)}$ is found to be $θ_μ(σ)=\frac{2}μ-σ$ whenever the long-ranged couplings are relevant $θ_μ(σ)=\frac{2}μ-σ\geq -1$. For the statistics of the ground state energy $E_L^{GS}$ over disordered samples, we obtain that the droplet exponent $θ_μ(σ) $ governs the leading correction to extensivity of the averaged value $\overline{E_L^{GS}} \simeq L e_0 +L^{θ_μ(σ)} e_1$. The characteristic scale of the fluctuations around this average is of order $L^{\frac{1}μ}$, and the rescaled variable $u=(E_L^{GS}-\overline{E_L^{GS}})/L^{\frac{1}μ}$ is Gaussian distributed for $μ=2$, or displays the negative power-law tail in $1/(-u)^{1+μ}$ for $u \to -\infty$ in the Lévy case $1<μ<2$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Scaling of the largest dynamical barrier in the one-dimensional long-range Ising spin-glass

The long-range one-dimensional Ising spin-glass with random couplings decaying as $J(r) \propto r^{-σ}$ presents a spin-glass phase $T_c(σ)>0$ for $0 \leq σ<1$ (the limit $σ=0$ corresponds to the mean-field SK-model). We use the eigenvalue method introduced in our previous work [C. Monthus and T. Garel, J. Stat. Mech. P12017 (2009)] to measure the equilibrium time $t_{eq}(N)$ at temperature $T=T_c(σ)/2$ as a function of the number $N$ of spins. We find the activated scaling $\ln t_{eq}(N) \sim N^ψ$ with the same barrier exponent $ψ\simeq 0.33$ in the whole region $0\leqσ<1$.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Zero-temperature spinglass-ferromagnetic transition : scaling analysis of the domain-wall energy

For the Ising model with Gaussian random coupling of average $J_0$ and unit variance, the zero-temperature spinglass-ferromagnetic transition as a function of the control parameter $J_0$ can be studied via the size-$L$ dependent renormalized coupling defined as the domain-wall energy $J^R(L) \equiv E_{GS}^{(AF)}(L)-E_{GS}^{(F)}(L)$ (i.e. the difference between the ground state energies corresponding to AntiFerromagnetic and and Ferromagnetic boundary conditions in one direction). We study numerically the critical exponents of this zero-temperature transition within the Migdal-Kadanoff approximation as a function of the dimension $d=2,3,4,5,6$. We then compare with the mean-field spherical model. Our main conclusion is that in low dimensions, the critical stiffness exponent $θ^c$ is clearly bigger than the spin-glass stiffness exponent $θ^{SG}$, but that they turn out to coincide in high enough dimension and in the mean-field spherical model. We also discuss the finite-size scaling properties of the averaged value and of the width of the distribution of the renormalized couplings.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Chaos properties of the one-dimensional long-range Ising spin-glass

For the long-range one-dimensional Ising spin-glass with random couplings decaying as $J(r) \propto r^{-σ}$, the scaling of the effective coupling defined as the difference between the free-energies corresponding to Periodic and Antiperiodic boundary conditions $J^R(N) \equiv F^{(P)}(N)-F^{(AP)}(N) \sim N^{θ(σ)}$ defines the droplet exponent $θ(σ)$. Here we study numerically the instability of the renormalization flow of the effective coupling $J^R(N)$ with respect to magnetic, disorder and temperature perturbations respectively, in order to extract the corresponding chaos exponents $ζ_H(σ)$, $ζ_J(σ)$ and $ζ_T(σ)$ as a function of $σ$. Our results for $ζ_T(σ) $ are interpreted in terms of the entropy exponent $θ_S(σ) \simeq 1/3$ which governs the scaling of the entropy difference $ S^{(P)}(N)-S^{(AP)}(N) \sim N^{θ_S(σ)}$. We also study the instability of the ground state configuration with respect to perturbations, as measured by the spin overlap between the unperturbed and the perturbed ground states, in order to extract the corresponding chaos exponents $ζ^{overlap}_H(σ)$ and $ζ^{overlap}_J(σ)$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamical barriers for the random ferromagnetic Ising model on the Cayley tree : traveling-wave solution of the real space renormalization flow

We consider the stochastic dynamics near zero-temperature of the random ferromagnetic Ising model on a Cayley tree of branching ratio $K$. We apply the Boundary Real Space Renormalization procedure introduced in our previous work (C. Monthus and T. Garel, J. Stat. Mech. P02037 (2013)) in order to derive the renormalization rule for dynamical barriers. We obtain that the probability distribution $P_n(B)$ of dynamical barrier for a subtree of $n$ generations converges for large $n$ towards some traveling-wave $P_n(B) \simeq P^*(B-nv) $, i.e. the width of the probability distribution remains finite around an average-value that grows linearly with the number $n$ of generations. We present numerical results for the branching ratios K=2 and K=3. We also compute the weak-disorder expansion of the velocity $v$ for K=2.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamical barriers of pure and random ferromagnetic Ising models on fractal lattices

We consider the stochastic dynamics of the pure and random ferromagnetic Ising model on the hierarchical diamond lattice of branching ratio $K$ with fractal dimension $d_f=(\ln (2K))/\ln 2$. We adapt the Real Space Renormalization procedure introduced in our previous work [C. Monthus and T. Garel, J. Stat. Mech. P02037 (2013)] to study the equilibrium time $t_{eq}(L)$ as a function of the system size $L$ near zero-temperature. For the pure Ising model, we obtain the behavior $t_{eq}(L) \sim L^α e^{β2J L^{d_s}} $ where $d_s=d_f-1$ is the interface dimension, and we compute the prefactor exponent $α$. For the random ferromagnetic Ising model, we derive the renormalization rules for dynamical barriers $B_{eq}(L) \equiv (\ln t_{eq}/β)$ near zero temperature. For the fractal dimension $d_f=2$, we obtain that the dynamical barrier scales as $ B_{eq}(L)= c L+L^{1/2} u$ where $u$ is a Gaussian random variable of non-zero-mean. While the non-random term scaling as $L$ corresponds to the energy-cost of the creation of a system-size domain-wall, the fluctuation part scaling as $L^{1/2}$ characterizes the barriers for the motion of the system-size domain-wall after its creation. This scaling corresponds to the dynamical exponent $ψ=1/2$, in agreement with the conjecture $ψ=d_s/2$ proposed in [C. Monthus and T. Garel, J. Phys. A 41, 115002 (2008)]. In particular, it is clearly different from the droplet exponent $θ\simeq 0.299$ involved in the statics of the random ferromagnet on the same lattice.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamics of Ising models near zero temperature : Real Space Renormalization Approach

We consider the stochastic dynamics of Ising ferromagnets (either pure or random) near zero temperature. The master equation satisfying detailed balance can be mapped onto a quantum Hamiltonian which has an exact zero-energy ground state representing the thermal equilibrium. The largest relaxation time $t_{eq}$ governing the convergence towards this Boltzmann equilibrium in finite-size systems is determined by the lowest non-vanishing eigenvalue $E_1=1/t_{eq}$ of the quantum Hamiltonian $H$. We introduce and study a real-space renormalization procedure for the quantum Hamiltonian associated to the single-spin-flip dynamics of Ising ferromagnets near zero temperature. We solve explicitly the renormalization flow for two cases. (i) For the one-dimensional random ferromagnetic chain with free boundary conditions, the largest relaxation time $t_{eq}$ can be expressed in terms of the set of random couplings for various choices of the dynamical transition rates. The validity of these RG results in $d=1$ is checked by comparison with another approach. (ii) For the pure Ising model on a Cayley tree of branching ratio $K$, we compute the exponential growth of $t_{eq}(N)$ with the number $N$ of generations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Typical versus averaged overlap distribution in Spin-Glasses : Evidence for the droplet scaling theory

We consider the statistical properties over disordered samples of the overlap distribution $P_{\cal J}(q)$ which plays the role of an order parameter in spin-glasses. We show that near zero temperature (i) the {\it typical} overlap distribution is exponentially small in the central region of $-1<q<1$: $ P^{typ}(q) = e^{\bar{\ln P_{\cal J}(q)}} \sim e^{- βN^θ ϕ(q)} $, where $θ$ is the droplet exponent defined here with respect to the total number $N$ of spins (in order to consider also fully connected models where the notion of length does not exist); (ii) the rescaled variable $v = - (\ln P_{\cal J}(q))/N^θ$ remains an O(1) random positive variable describing sample-to sample fluctuations; (iii) the averaged distribution $\bar{P_{\cal J}(q)} $ is non-typical and dominated by rare anomalous samples. Similar statements hold for the cumulative overlap distribution $I_{\cal J}(q_0) \equiv \int_{0}^{q_0} dq P_{\cal J}(q) $. These results are derived explicitly for the spherical mean-field model with $θ=1/3$, $ϕ(q)=1-q^2 $, and the random variable $v$ corresponds to the rescaled difference between the two largest eigenvalues of GOE random matrices. Then we compare numerically the typical and averaged overlap distributions for the long-ranged one-dimensional Ising spin-glass with random couplings decaying as $J(r) \propto r^{-σ}$ for various values of the exponent $σ$, corresponding to various droplet exponents $θ(σ)$, and for the mean-field SK-model (corresponding formally to the $σ=0$ limit of the previous model). Our conclusion is that future studies on spin-glasses should measure the {\it typical} values of the overlap distribution or of the cumulative overlap distribution to obtain clearer conclusions on the nature of the spin-glass phase.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamical Barriers in the Dyson Hierarchical model via Real Space Renormalization

The Dyson hierarchical one-dimensional Ising model of parameter $σ>0$ contains long-ranged ferromagnetic couplings decaying as $1/r^{1+σ}$ in terms of the distance $r$. We study the stochastic dynamics near zero-temperature via the Real Space Renormalization introduced in our previous work (C. Monthus and T. Garel, arxiv:1212.0643) in order to compute explicitly the equilibrium time $t_{eq}(L)$ as a function of the system size $L$. For $σ<1$ where the static critical temperature for the ferromagnetic transition is finite $T_c>0$, we obtain that dynamical barriers grow as the power-law: $\ln t_{eq}(L) \simeq β(\frac{4 J_0}{3(2^{1-σ}-1)}) L^{1-σ}$. For $σ=1$ where the static critical temperature vanishes $T_c=0$, we obtain that dynamical barriers grow logarithmically as : $\ln t_{eq}(L) \simeq [β(\frac{4 J_0}{3 \ln 2}) -1] \ln L $. We also compute finite contributions to the dynamical barriers that can depend on the choice of transition rates satisfying detailed balance.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Random Transverse Field Ising model in $d=2$ : analysis via Boundary Strong Disorder Renormalization

To avoid the complicated topology of surviving clusters induced by standard Strong Disorder RG in dimension $d>1$, we introduce a modified procedure called 'Boundary Strong Disorder RG' where the order of decimations is chosen a priori. We apply numerically this modified procedure to the Random Transverse Field Ising model in dimension $d=2$. We find that the location of the critical point, the activated exponent $ψ\simeq 0.5$ of the Infinite Disorder scaling, and the finite-size correlation exponent $ν_{FS} \simeq 1.3$ are compatible with the values obtained previously by standard Strong Disorder RG.Our conclusion is thus that Strong Disorder RG is very robust with respect to changes in the order of decimations. In addition, we analyze in more details the RG flows within the two phases to show explicitly the presence of various correlation length exponents : we measure the typical correlation exponent $ν_{typ} \simeq 0.64$ in the disordered phase (this value is very close to the correlation exponent $ν^Q_{pure}(d=2) \simeq 0.63$ of the {\it pure} two-dimensional quantum Ising Model), and the typical exponent $ν_h \simeq 1$ within the ordered phase. These values satisfy the relations between critical exponents imposed by the expected finite-size scaling properties at Infinite Disorder critical points. Within the disordered phase, we also measure the fluctuation exponent $ω\simeq 0.35$ which is compatible with the Directed Polymer exponent $ω_{DP}(1+1)=1/3$ in $(1+1)$ dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Random Transverse Field Ising Model in dimension $d>1$ : scaling analysis in the disordered phase from the Directed Polymer model

For the quantum Ising model with ferromagnetic random couplings $J_{i,j}>0$ and random transverse fields $h_i>0$ at zero temperature in finite dimensions $d>1$, we consider the lowest-order contributions in perturbation theory in $(J_{i,j}/h_i)$ to obtain some information on the statistics of various observables in the disordered phase. We find that the two-point correlation scales as : $\ln C(r) \sim - \frac{r}{ξ_{typ}} +r^ω u$, where $ξ_{typ} $ is the typical correlation length, $u$ is a random variable, and $ω$ coincides with the droplet exponent $ω_{DP}(D=d-1)$ of the Directed Polymer with $D=(d-1)$ transverse directions. Our main conclusions are (i) whenever $ω>0$, the quantum model is governed by an Infinite-Disorder fixed point : there are two distinct correlation length exponents related by $ν_{typ}=(1-ω)ν_{av}$ ; the distribution of the local susceptibility $χ_{loc}$ presents the power-law tail $P(χ_{loc}) \sim 1/χ_{loc}^{1+μ}$ where $μ$ vanishes as $ξ_{av}^{-ω} $, so that the averaged local susceptibility diverges in a finite neighborhood $0<μ<1$ before criticality (Griffiths phase) ; the dynamical exponent $z$ diverges near criticality as $z=d/μ\sim ξ_{av}^ω$ (ii) in dimensions $d \leq 3$, any infinitesimal disorder flows towards this Infinite-Disorder fixed point with $ω(d)>0$ (for instance $ω(d=2)=1/3$ and $ω(d=3) \sim 0.24$) (iii) in finite dimensions $d > 3$, a finite disorder strength is necessary to flow towards the Infinite-Disorder fixed point with $ω(d)>0$ (for instance $ω(d=4) \simeq 0.19$), whereas a Finite-Disorder fixed point remains possible for a small enough disorder strength. For the Cayley tree of effective dimension $d=\infty$ where $ω=0$, we discuss the similarities and differences with the case of finite dimensions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Random Transverse Field Ising model on the Cayley Tree : analysis via Boundary Strong Disorder Renormalization

Strong Disorder Renormalization for the Random Transverse Field Ising model leads to a complicated topology of surviving clusters as soon as $d>1$. Even if one starts from a Cayley tree, the network of surviving renormalized clusters will contain loops, so that no analytical solution can been obtained. Here we introduce a modified procedure called 'Boundary Strong Disorder Renormalization' that preserves the tree structure, so that one can write simple recursions with respect to the number of generations. We first show that this modified procedure allows to recover exactly most of the critical exponents for the one-dimensional chain. After this important check, we study the RG equations for the quantum Ising model on a Cayley tree with a uniform ferromagnetic coupling $J$ and random transverse fields with support $[h_{min},h_{max}]$. We find the following picture (i) for $J>h_{max}$, only bonds are decimated, so that the whole tree is a quantum ferromagnetic cluster (ii) for $J<h_{min}$, only sites are decimated, so that no quantum ferromagnetic cluster is formed, and the ferromagnetic coupling to the boundary coincides with the partition function of a Directed Polymer model in a random medium (iii) for $h_{min}<J<h_{max}$, both sites and bonds can be decimated : the quantum ferromagnetic clusters can either remain finite (the physics is then similar to (ii), with a quantitative mapping to a modified Directed Polymer model) or an infinite quantum ferromagnetic cluster appears. We find that the quantum transition can be of two types : (a) either the quantum transition takes place in the region where quantum ferromagnetic clusters remain finite, and the singularity of the ferromagnetic coupling to the boundary involves the typical correlation length exponent $ν_{typ}=1$ (b) or the quantum transition takes place at the point where an extensive quantum ferromagnetic cluster appears.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Strong Disorder RG principles within a fixed cell-size real space renormalization : application to the Random Transverse Field Ising model on various fractal lattices

Strong Disorder Renormalization is an energy-based renormalization that leads to a complicated renormalized topology for the surviving clusters as soon as $d>1$. In this paper, we propose to include Strong Disorder Renormalization ideas within the more traditional fixed cell-size real space RG framework. We first consider the one-dimensional chain as a test for this fixed cell-size procedure: we find that all exactly known critical exponents are reproduced correctly, except for the magnetic exponent $β$ (because it is related to more subtle persistence properties of the full RG flow). We then apply numerically this fixed cell-size procedure to two types of renormalizable fractal lattices (i) the Sierpinski gasket of fractal dimension $D=\ln 3/\ln 2$, where there is no underlying classical ferromagnetic transition, so that the RG flow in the ordered phase is similar to what happens in $d=1$ (ii) a hierarchical diamond lattice of fractal dimension $D=4/3$, where there is an underlying classical ferromagnetic transition, so that the RG flow in the ordered phase is similar to what happens on hypercubic lattices of dimension $d>1$. In both cases, we find that the transition is governed by an Infinite Disorder Fixed Point : besides the measure of the activated exponent $ψ$, we analyze the RG flow of various observables in the disordered and ordered phases, in order to extract the 'typical' correlation length exponents of these two phases which are different from the finite-size correlation length exponent.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A critical Dyson hierarchical model for the Anderson localization transition

A Dyson hierarchical model for Anderson localization, containing non-random hierarchical hoppings and random on-site energies, has been studied in the mathematical literature since its introduction by Bovier [J. Stat. Phys. 59, 745 (1990)], with the conclusion that this model is always in the localized phase. Here we show that if one introduces alternating signs in the hoppings along the hierarchy (instead of choosing all hoppings of the same sign), it is possible to reach an Anderson localization critical point presenting multifractal eigenfunctions and intermediate spectral statistics. The advantage of this model is that one can write exact renormalization equations for some observables. In particular, we obtain that the renormalized on-site energies have the Cauchy distributions for exact fixed points. Another output of this renormalization analysis is that the typical exponent of critical eigenfunctions is always $α_{typ}=2$, independently of the disorder strength. We present numerical results concerning the whole multifractal spectrum $f(α)$ and the compressibility $χ$ of the level statistics, both for the box and the Cauchy distributions of the random on-site energies. We discuss the similarities and differences with the ensemble of ultrametric random matrices introduced recently by Fyodorov, Ossipov and Rodriguez [J. Stat. Mech. L12001 (2009)].

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Anderson localization on the Cayley tree : multifractal statistics of the transmission at criticality and off criticality

In contrast to finite dimensions where disordered systems display multifractal statistics only at criticality, the tree geometry induces multifractal statistics for disordered systems also off criticality. For the Anderson tight-binding localization model defined on a tree of branching ratio K=2 with $N$ generations, we consider the Miller-Derrida scattering geometry [J. Stat. Phys. 75, 357 (1994)], where an incoming wire is attached to the root of the tree, and where $K^{N}$ outcoming wires are attached to the leaves of the tree. In terms of the $K^{N}$ transmission amplitudes $t_j$, the total Landauer transmission is $T \equiv \sum_j | t_j |^2$, so that each channel $j$ is characterized by the weight $w_j=| t_j |^2/T$. We numerically measure the typical multifractal singularity spectrum $f(α)$ of these weights as a function of the disorder strength $W$ and we obtain the following conclusions for its left-termination point $α_+(W)$. In the delocalized phase $W<W_c$, $α_+(W)$ is strictly positive $α_+(W)>0$ and is associated with a moment index $q_+(W)>1$. At criticality, it vanishes $α_+(W_c)=0$ and is associated with the moment index $q_+(W_c)=1$. In the localized phase $W>W_c$, $α_+(W)=0$ is associated with some moment index $q_+(W)<1$. We discuss the similarities with the exact results concerning the multifractal properties of the Directed Polymer on the Cayley tree.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Non-equilibrium steady states : maximization of the Shannon entropy associated to the distribution of dynamical trajectories in the presence of constraints

Filyokov and Karpov [Inzhenerno-Fizicheskii Zhurnal 13, 624 (1967)] have proposed a theory of non-equilibrium steady states in direct analogy with the theory of equilibrium states : the principle is to maximize the Shannon entropy associated to the probability distribution of dynamical trajectories in the presence of constraints, including the macroscopic current of interest, via the method of Lagrange multipliers. This maximization leads directly to generalized Gibbs distribution for the probability distribution of dynamical trajectories, and to some fluctuation relation of the integrated current. The simplest stochastic dynamics where these ideas can be applied are discrete-time Markov chains, defined by transition probabilities $W_{i \to j}$ between configurations $i$ and $j$ : instead of choosing the dynamical rules $W_{i \to j} $ a priori, one determines the transition probabilities and the associate stationary state that maximize the entropy of dynamical trajectories with the other physical constraints that one wishes to impose. We give a self-contained and unified presentation of this type of approach, both for discrete-time Markov Chains and for continuous-time Master Equations. The obtained results are in full agreement with the Bayesian approach introduced by Evans [Phys. Rev. Lett. 92, 150601 (2004)] under the name 'Non-equilibrium Counterpart to detailed balance', and with the 'invariant quantities' derived by Baule and Evans [Phys. Rev. Lett. 101, 240601 (2008)], but provide a slightly different perspective via the formulation in terms of an eigenvalue problem.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Random elastic networks : strong disorder renormalization approach

For arbitrary networks of random masses connected by random springs, we define a general strong disorder real-space renormalization (RG) approach that generalizes the procedures introduced previously by Hastings [Phys. Rev. Lett. 90, 148702 (2003)] and by Amir, Oreg and Imry [Phys. Rev. Lett. 105, 070601 (2010)] respectively. The principle is to eliminate iteratively the elementary oscillating mode of highest frequency associated with either a mass or a spring constant. To explain the accuracy of the strong disorder RG rules, we compare with the Aoki RG rules that are exact at fixed frequency.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Random field Ising model : statistical properties of low-energy excitations and of equilibrium avalanches

With respect to usual thermal ferromagnetic transitions, the zero-temperature finite-disorder critical point of the Random-field Ising model (RFIM) has the peculiarity to involve some 'droplet' exponent $θ$ that enters the generalized hyperscaling relation $2-α= ν(d-θ)$. In the present paper, to better understand the meaning of this droplet exponent $θ$ beyond its role in the thermodynamics, we discuss the statistics of low-energy excitations generated by an imposed single spin-flip with respect to the ground state, as well as the statistics of equilibrium avalanches i.e. the magnetization jumps that occur in the sequence of ground-states as a function of the external magnetic field. The droplet scaling theory predicts that the distribution $dl/l^{1+θ}$ of the linear-size $l$ of low-energy excitations transforms into the distribution $ds/s^{1+{θ/d_f}}$ for the size $s$ (number of spins) of excitations of fractal dimension $d_f$ ($s \sim l^{d_f}$). In the non-mean-field region $d<d_c$, droplets are compact $d_f=d$, whereas in the mean-field region $d>d_c$, droplets have a fractal dimension $d_f=2 θ$ leading to the well-known mean-field result $ds/s^{3/2}$. Zero-field equilibrium avalanches are expected to display the same distribution $ds/s^{1+{θ/d_f}}$. We also discuss the statistics of equilibrium avalanches integrated over the external field and finite-size behaviors. These expectations are checked numerically for the Dyson hierarchical version of the RFIM, where the droplet exponent $θ(σ)$ can be varied as a function of the effective long-range interaction $J(r) \sim 1/r^{d+σ}$ in $d=1$.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Random Transverse Field Ising Model in dimension $d=2,3$ : Infinite Disorder scaling via a non-linear transfer approach

The 'Cavity-Mean-Field' approximation developed for the Random Transverse Field Ising Model on the Cayley tree [L. Ioffe and M. Mézard, PRL 105, 037001 (2010)] has been found to reproduce the known exact result for the surface magnetization in $d=1$ [O. Dimitrova and M. Mézard, J. Stat. Mech. (2011) P01020]. In the present paper, we propose to extend these ideas in finite dimensions $d>1$ via a non-linear transfer approach for the surface magnetization. In the disordered phase, the linearization of the transfer equations correspond to the transfer matrix for a Directed Polymer in a random medium of transverse dimension $D=d-1$, in agreement with the leading order perturbative scaling analysis [C. Monthus and T. Garel, arxiv:1110.3145]. We present numerical results of the non-linear transfer approach in dimensions $d=2$ and $d=3$. In both cases, we find that the critical point is governed by Infinite Disorder scaling. In particular exactly at criticality, the one-point surface magnetization scales as $\ln m_L^{surf} \simeq - L^{ω_c} v$, where $ω_c(d)$ coincides with the droplet exponent $ω_{DP}(D=d-1)$ of the corresponding Directed Polymer model, with $ω_c(d=2)=1/3$ and $ω_c(d=3) \simeq 0.24$. The distribution $P(v)$ of the positive random variable $v$ of order O(1) presents a power-law singularity near the origin $P(v) \propto v^a$ with $a(d=2,3)>0$ so that all moments of the surface magnetization are governed by the same power-law decay $\bar{(m_L^{surf})^k} \propto L^{- x_S}$ with $x_S=ω_c (1+a)$ independently of the order $k$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Anderson localization of phonons in dimension $d=1,2,3$ : finite-size properties of the Inverse Participation Ratios of eigenstates

We study by exact diagonalization the localization properties of phonons in mass-disordered harmonic crystals of dimension $d=1,2,3$. We focus on the behavior of the typical Inverse Participation Ratio $Y_2(ω,L)$ as a function of the frequency $ω$ and of the linear length $L$ of the disordered samples. In dimensions $d=1$ and $d=2$, we find that the low-frequency part $ω\to 0$ of the spectrum satisfies the following finite-size scaling $L Y_2(ω,L)=F_{d=1}(L^{1/2} ω)$ in dimension $d=1$ and $L^2 Y_2(ω,L)=F_{d=2}((\ln L)^{1/2} ω)$ in dimension $d=2$, with the following conclusions (i) an eigenstate of any fixed frequency $ω$ becomes localized in the limit $L \to +\infty$ (ii) a given disordered sample of size $L^d$ contains a number $N_{deloc}(L)$ of delocalized states growing as $N_{deloc}(L)\sim L^{1/2}$ in $d=1$ and as $N_{deloc}(L)\sim L^2/(\ln L)$ in $d=2$. In dimension $d=3$, we find a localization-delocalization transition at some finite critical frequency $ω_c(W)>0$ (that depends on the disorder strength $W$). Our data are compatible with the finite-size scaling $L^{D(2)} Y_2(ω,L)=F_{d=3}(L^{1/ν} (ω-ω_c))$ with the values $D(2) \simeq 1.3$ and $ν\simeq 1.57$ corresponding to the universality class of the localization transition for the Anderson tight-binding electronic model in dimension $d=3$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Anderson localization transition with long-ranged hoppings : analysis of the strong multifractality regime in terms of weighted Levy sums

For Anderson tight-binding models in dimension $d$ with random on-site energies $ε_{\vec r}$ and critical long-ranged hoppings decaying typically as $V^{typ}(r) \sim V/r^d$, we show that the strong multifractality regime corresponding to small $V$ can be studied via the standard perturbation theory for eigenvectors in quantum mechanics. The Inverse Participation Ratios $Y_q(L)$, which are the order parameters of Anderson transitions, can be written in terms of weighted Lévy sums of broadly distributed variables (as a consequence of the presence of on-site random energies in the denominators of the perturbation theory). We compute at leading order the typical and disorder-averaged multifractal spectra $τ_{typ}(q)$ and $τ_{av}(q)$ as a function of $q$. For $q<1/2$, we obtain the non-vanishing limiting spectrum $τ_{typ}(q)=τ_{av}(q)=d(2q-1)$ as $V \to 0^+$. For $q>1/2$, this method yields the same disorder-averaged spectrum $τ_{av}(q)$ of order $O(V)$ as obtained previously via the Levitov renormalization method by Mirlin and Evers [Phys. Rev. B 62, 7920 (2000)]. In addition, it allows to compute explicitly the typical spectrum, also of order $O(V)$, but with a different $q$-dependence $τ_{typ}(q) \ne τ_{av}(q)$ for all $q>q_c=1/2$. As a consequence, we find that the corresponding singularity spectra $f_{typ}(α)$ and $f_{av}(α)$ differ even in the positive region $f>0$, and vanish at different values $α_+^{typ} > α_+^{av}$, in contrast to the standard picture. We also obtain that the saddle value $α_{typ}(q)$ of the Legendre transform reaches the termination point $α_+^{typ}$ where $f_{typ}(α_+^{typ})=0 $ only in the limit $q \to +\infty$.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Many-body localization transition in a lattice model of interacting fermions: statistics of renormalized hoppings in configuration space

We consider the one-dimensional lattice model of interacting fermions with disorder studied previously by Oganesyan and Huse [Phys. Rev. B 75, 155111 (2007)]. To characterize a possible many-body localization transition as a function of the disorder strength $W$, we use an exact renormalization procedure in configuration space that generalizes the Aoki real-space RG procedure for Anderson localization one-particle models [H. Aoki, J. Phys. C13, 3369 (1980)]. We focus on the statistical properties of the renormalized hopping $V_L$ between two configurations separated by a distance $L$ in configuration space (distance being defined as the minimal number of elementary moves to go from one configuration to the other). Our numerical results point towards the existence of a many-body localization transition at a finite disorder strength $W_c$. In the localized phase $W>W_c$, the typical renormalized hopping $V_L^{typ} \equiv e^{\bar{\ln V_L}}$ decays exponentially in $L$ as $ (\ln V_L^{typ}) \simeq - \frac{L}{ξ_{loc}}$ and the localization length diverges as $ξ_{loc}(W) \sim (W-W_c)^{-ν_{loc}}$ with a critical exponent of order $ν_{loc} \simeq 0.5$. In the delocalized phase $W<W_c$, the renormalized hopping remains a finite random variable as $L \to \infty$, and the typical asymptotic value $V_{\infty}^{typ} \equiv e^{\bar{\ln V_{\infty}}}$ presents an essential singularity $(\ln V_{\infty}^{typ}) \sim - (W_c-W)^{-κ}$ with an exponent of order $κ\sim 1.4$. Finally, we show that this analysis in configuration space is compatible with the localization properties of the simplest two-point correlation function in real space.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Matching between typical fluctuations and large deviations in disordered systems : application to the statistics of the ground state energy in the SK spin-glass model

For the statistics of global observables in disordered systems, we discuss the matching between typical fluctuations and large deviations. We focus on the statistics of the ground state energy $E_0$ in two types of disordered models : (i) for the directed polymer of length $N$ in a two-dimensional medium, where many exact results exist (ii) for the Sherrington-Kirkpatrick spin-glass model of $N$ spins, where various possibilities have been proposed. Here we stress that, besides the behavior of the disorder-average $E_0^{av}(N)$ and of the standard deviation $ ΔE_0(N) \sim N^{ω_f}$ that defines the fluctuation exponent $ω_f$, it is very instructive to study the full probability distribution $Π(u)$ of the rescaled variable $u= \frac{E_0(N)-E_0^{av}(N)}{ΔE_0(N)}$ : (a) numerically, the convergence towards $Π(u)$ is usually very rapid, so that data on rather small sizes but with high statistics allow to measure the two tails exponents $η_{\pm}$ defined as $\ln Π(u \to \pm \infty) \sim - | u |^{η_{\pm}}$. In the generic case $1< η_{\pm} < +\infty$, this leads to explicit non-trivial terms in the asymptotic behaviors of the moments $\bar{Z_N^n}$ of the partition function when the combination $[| n | N^{ω_f}]$ becomes large (b) simple rare events arguments can usually be found to obtain explicit relations between $η_{\pm}$ and $ω_f$. These rare events usually correspond to 'anomalous' large deviation properties of the generalized form $R(w_{\pm} = \frac{E_0(N)-E_0^{av}(N)}{N^{κ_{\pm}}}) \sim e^{- N^{ρ_{\pm}} {\cal R}_{\pm}(w_{\pm})}$ (the 'usual' large deviations formalism corresponds to $κ_{\pm}=1=ρ_{\pm}$).

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Random cascade models of multifractality : real-space renormalization and travelling-waves

Random multifractals occur in particular at critical points of disordered systems. For Anderson localization transitions, Mirlin and Evers [PRB 62,7920 (2000)] have proposed the following scenario (a) the Inverse Participation Ratios (I.P.R.) $Y_q(L)$ display the following fluctuations between the disordered samples of linear size $L$ : with respect to the typical value $Y^{typ}_q(L) = e^{\bar{\ln Y_q(L)}} \sim L^{- τ_{typ}(q)} $ that involve the typical multifractal spectrum $τ_{typ}(q)$, the rescaled variable $y=Y_q(L)/Y^{typ}_q(L) $ is distributed with a scale-invariant distribution presenting the power-law tail $1/y^{1+β_q}$, so that the disorder-averaged I.P.R. $\bar{Y_q(L)} \sim L^{- τ_{av}(q)} $ have multifractal exponents $τ_{av}(q) $ that differ from the typical ones $τ_{typ}(q)$ whenever $β_q<1$; (b) the tail exponents $β_q$ and the multifractal exponents are related by the relation $β_q τ_{typ}(q)=τ_{av}(q β_q)$. Here we show that this scenario can be understood by considering the real-space renormalization equations satisfied by the I.P.R. For the simplest multifractals described by random cascades, these renormalization equations are formally similar to the recursion relations for disordered models defined on Cayley trees and they admit travelling-wave solutions for the variable $(\ln Y_q)$ in the effective time $t_{eff}=\ln L$ : the exponent $τ_{typ}(q)$ represents the velocity, whereas the tail exponent $β_q$ represents the usual exponential decay of the travelling-wave tail. In addition, we obtain that the relation (b) above can be obtained as a self-consistency condition from the self-similarity of the multifractal spectrum at all scales. Our conclusion is thus that the Mirlin-Evers scenario should apply to other types of random critical points, and even to random multifractals occurring in other fields.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Random walk in a two-dimensional self-affine random potential : properties of the anomalous diffusion phase at small external force

We consider the random walk of a particle in a two-dimensional self-affine random potential of Hurst exponent $H=1/2$ in the presence of an external force $F$. We present numerical results on the statistics of first-passage times that satisfy closed backward master equations. We find that there exists a zero-velocity phase in a finite region of the external force $0<F<F_c$, where the dynamics follows the anomalous diffusion law $ x(t) \sim ξ(F) \ t^{μ(F)} $. The anomalous exponent $0<μ(F)<1$ and the correlation length $ξ(F)$ vary continuously with $F$. In the limit of vanishing force $F \to 0$, we measure the following power-laws : the anomalous exponent vanishes as $μ(F) \propto F^a$ with $a \simeq 0.6$ (instead of $a=1$ in dimension $d=1$), and the correlation length diverges as $ξ(F) \propto F^{-ν}$ with $ν\simeq 1.29$ (instead of $ν=2$ in dimension $d=1$). Our main conclusion is thus that the dynamics renormalizes onto an effective directed trap model, where the traps are characterized by a typical length $ξ(F)$ along the direction of the force, and by a typical barrier $1/μ(F)$. The fact that these traps are 'smaller' in linear size and in depth than in dimension $d=1$, means that the particle uses the transverse direction to find lower barriers.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Random walk in two-dimensional self-affine random potentials : strong disorder renormalization approach

We consider the continuous-time random walk of a particle in a two-dimensional self-affine quenched random potential of Hurst exponent $H>0$. The corresponding master equation is studied via the strong disorder renormalization procedure introduced in Ref. [C. Monthus and T. Garel, J. Phys. A: Math. Theor. 41 (2008) 255002]. We present numerical results on the statistics of the equilibrium time $t_{eq}$ over the disordered samples of a given size $L \times L$ for $10 \leq L \leq 80$. We find an 'Infinite disorder fixed point', where the equilibrium barrier $Γ_{eq} \equiv \ln t_{eq}$ scales as $Γ_{eq}=L^H u $ where $u$ is a random variable of order O(1). This corresponds to a logarithmically-slow diffusion $ | \vec r(t) - \vec r(0) | \sim (\ln t)^{1/H}$ for the position $\vec r(t)$ of the particle.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Eigenvalue method to compute the largest relaxation time of disordered systems

We consider the dynamics of finite-size disordered systems as defined by a master equation satisfying detailed balance. The master equation can be mapped onto a Schrödinger equation in configuration space, where the quantum Hamiltonian $H$ has the generic form of an Anderson localization tight-binding model. The largest relaxation time $t_{eq}$ governing the convergence towards Boltzmann equilibrium is determined by the lowest non-vanishing eigenvalue $E_1=1/t_{eq}$ of $H$ (the lowest eigenvalue being $E_0=0$). So the relaxation time $t_{eq}$ can be computed {\it without simulating the dynamics} by any eigenvalue method able to compute the first excited energy $E_1$. Here we use the 'conjugate gradient' method to determine $E_1$ in each disordered sample and present numerical results on the statistics of the relaxation time $t_{eq}$ over the disordered samples of a given size for two models : (i) for the random walk in a self-affine potential of Hurst exponent $H$ on a two-dimensional square of size $L \times L$, we find the activated scaling $\ln t_{eq}(L) \sim L^ψ$ with $ψ=H$ as expected; (ii) for the dynamics of the Sherrington-Kirkpatrick spin-glass model of $N$ spins, we find the growth $\ln t_{eq}(N) \sim N^ψ$ with $ψ=1/3$ in agreement with most previous Monte-Carlo measures. In addition, we find that the rescaled distribution of $(\ln t_{eq})$ decays as $e^{- u^η}$ for large $u$ with a tail exponent of order $η\simeq 1.36$. We give a rare-event interpretation of this value, that points towards a sample-to-sample fluctuation exponent of order $ψ_{width} \simeq 0.26$ for the barrier.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Statistics of first-passage times in disordered systems using backward master equations and their exact renormalization rules

We consider the non-equilibrium dynamics of disordered systems as defined by a master equation involving transition rates between configurations (detailed balance is not assumed). To compute the important dynamical time scales in finite-size systems without simulating the actual time evolution which can be extremely slow, we propose to focus on first-passage times that satisfy 'backward master equations'. Upon the iterative elimination of configurations, we obtain the exact renormalization rules that can be followed numerically. To test this approach, we study the statistics of some first-passage times for two disordered models : (i) for the random walk in a two-dimensional self-affine random potential of Hurst exponent $H$, we focus on the first exit time from a square of size $L \times L$ if one starts at the square center. (ii) for the dynamics of the ferromagnetic Sherrington-Kirkpatrick model of $N$ spins, we consider the first passage time $t_f$ to zero-magnetization when starting from a fully magnetized configuration. Besides the expected linear growth of the averaged barrier $\bar{\ln t_{f}} \sim N$, we find that the rescaled distribution of the barrier $(\ln t_{f})$ decays as $e^{- u^η}$ for large $u$ with a tail exponent of order $η\simeq 1.72$. This value can be simply interpreted in terms of rare events if the sample-to-sample fluctuation exponent for the barrier is $ψ_{width}=1/3$.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Disorder-dominated phases of random systems : relations between tails exponents and scaling exponents

We consider various random models (directed polymer, ferromagnetic random Potts model, Ising spin-glasses) in their disorder-dominated phases, where the free-energy cost $F(L)$ of an excitation of length $L$ present fluctuations that grow as a power-law $ΔF(L) \sim L^ω$ with the so-called droplet exponent $ω>0$. We study the tails of the probability distribution $Π(x)$ of the rescaled free-energy cost $x= \frac{F_L- \bar{F_L}}{L^ω}$, which are governed by two exponents $(η_-,η_+)$ defined by $\ln Π(x \to \pm \infty) \sim - | x |^{η_{\pm}}$. The aim of this paper is to establish simple relations between these tail exponents $(η_-,η_+) $ and the droplet exponent $ω$. We first prove these relations for disordered models on diamond hierarchical lattices where exact renormalizations exist for the probability distribution $Π(x)$. We then interpret these relations via an analysis of the measure of the rare disorder configurations governing the tails. Our conclusion is that these relations, when expressed in terms of the dimensions of the bulk and of the excitation surface are actually valid for general lattices.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Localization properties of the anomalous diffusion phase $x ~ t^μ$ in the directed trap model and in the Sinai diffusion with bias

We study the anomalous diffusion phase $x ~ t^μ$ with $0<μ<1$ which exists both in the Sinai diffusion at small bias, and in the related directed trap model presenting a large distribution of trapping time $p(τ) \sim 1/τ^{1+μ}$. Our starting point is the Real Space Renormalization method in which the whole thermal packet is considered to be in the same renormalized valley at large time : this assumption is exact only in the limit $μ\to 0$ and corresponds to the Golosov localization. For finite $μ$, we thus generalize the usual RSRG method to allow for the spreading of the thermal packet over many renormalized valleys. Our construction allows to compute exact series expansions in $μ$ of all observables : at order $μ^n$, it is sufficient to consider a spreading of the thermal packet onto at most $(1+n)$ traps in each sample, and to average with the appropriate measure over the samples. For the directed trap model, we show explicitly up to order $μ^2$ how to recover the diffusion front, the thermal width, and the localization parameter $Y_2$. We moreover compute the localization parameters $Y_k$ for arbitrary $k$, the correlation function of two particles, and the generating function of thermal cumulants. We then explain how these results apply to the Sinai diffusion with bias, by deriving the quantitative mapping between the large-scale renormalized descriptions of the two models.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

On the localization of random heteropolymers at the interface between two selective solvents

To study the localization of random heteropolymers at an interface separating two selective solvents within the model of Garel, Huse, Leibler and Orland, Europhys. Lett. {\bf 8} 9 (1989), we propose an approach based on a disorder-dependent real space renormalization procedure. This approach allows to recover that a chain with a symmetric distribution in hydrophobic/hydrophilic components is localized at any temperature in the thermodynamic limit, whereas a dissymmetric distribution in hydrophobic/hydrophilic components leads to a delocalization phase transition. It yields in addition explicit expressions for the thermodynamic quantities as well as a very detailed description of the statistical properties of the behaviors of the heteropolymers in the high temperature limit. For the case of a small dissymmetry in hydrophobic/hydrophilic components, the renormalization approach yields explicit predictions for the delocalization transition temperature and for the critical behaviors of various quantities : in particular, the free energy presents an essential singularity at the transition, the typical length of blobs in the preferred solvent diverges with an essential singularity, whereas the typical length of blobs in the other solvent diverges algebraically. Finite-size properties are also characterized in details for both cases. In particular, we give the probability distribution of the delocalization temperature for the ensemble of chains of finite (large) length $L$. Finally, we discuss the non-equilibrium dynamics at temperature $T$ starting from a zero-temperature initial condition.

0 citations0 reviewsNo code linkedOpen access