Despite its importance to experiments, numerical simulations, and the development of theoretical models, self-averaging in many-body quantum systems out of equilibrium remains underinvestigated. Usually, in the chaotic regime, self-averaging is taken for granted. The numerical and analytical results presented here force us to rethink these expectations. They demonstrate that self-averaging properties depend on the quantity and also on the time scale considered. We show analytically that the survival probability in chaotic systems is not self-averaging at any time scale, even when evolved under full random matrices. We also analyze the participation ratio, R\'enyi entropies, the spin autocorrelation function from experiments with cold atoms, and the connected spin-spin correlation function from experiments with ion traps. We find that self-averaging holds at short times for the quantities that are local in space, while at long times, self-averaging applies for quantities that are local in time. Various behaviors are revealed at intermediate time scales.
The self-averaging behavior of interacting many-body quantum systems has been mostly studied at equilibrium. The present paper addresses what happens out of equilibrium, as the increase of the strength of on-site disorder takes the system to the localized phase. We consider two local and two nonlocal quantities of great experimental and theoretical interest. In the delocalized phase, self-averaging depends on the observable and on the timescale, but the picture simplifies substantially when localization is reached. In the localized phase, the local observables become self-averaging at all times while the nonlocal quantities are throughout non-self-averaging. These behaviors are explained and scaling analysis is provided using the $\ensuremath{\ell}$-bit model and a toy model.
We investigate how the dynamical fluctuations of many-body quantum systems out of equilibrium can be mitigated when they are opened to a dephasing environment. We consider the survival probability (spectral form factor with a filter) evolving under different kinds of random matrices and under a spin-1/2 model with weak and strong disorder. In isolated many-body quantum systems, the survival probability is non-self-averaging at any timescale, that is, the relative variance of its fluctuations does not decrease with system size. By opening the system, we find that the fluctuations are always reduced, but self-averaging can only be ensured away from critical points. Self-averaging is achieved for the long-time dynamics of full random matrices, power-law banded random matrices deep in the delocalized phase, and the Rosenzweig-Porter ensemble in all phases except at the localization transition point. For the spin model, the survival probability becomes self-averaging only in the chaotic regime provided the initial states are in the middle of the spectrum. Overall, a strongly non-self-averaging survival probability in open systems is an indicator of criticality.
An observable of a disordered system is self-averaging when its properties do not depend on the specific realization considered. Lack of self-averaging, on the other hand, implies that sample to sample fluctuations persist no matter how large the system is. The latter scenario is often found in the vicinity of critical points, such as at the metal-insulator transition of interacting many-body quantum systems. Much attention has been devoted to these systems at equilibrium, but little is known about their self-averaging behavior out of equilibrium, which is the subject of this work. We consider two local and two non-local quantities in real space that are of great experimental and theoretical interest. In the metallic phase, we show that their self-averaging behavior is highly dependent on the observable itself and on the time scale, but the picture simplifies substantially as we approach localization. In this phase, the local quantities are self-averaging at any time, while the non-local ones are non-self-averaging at all time scales.
We consider the sample to sample fluctuations that occur in the value of a thermodynamic quantity
P in an ensemble of finite systems with quenched disorder, at equilibrium. The variance of
P,
VP, which characterizes these fluctuations is calculated as a function of the systems' linear size
l, focusing on the behavior at the critical point. The specific model considered is the bond-disordered Ashkin-Teller model on a square lattice [Phys. Rev. 64, 178 (1943)]. Using extensive Monte Carlo simulations, several bond-disordered Ashkin-Teller models were examined, including the bond-disordered Ising model and the bond-disordered four-state Potts model. It was found that far from criticality all thermodynamic quantities which were examined (energy, magnetization, specific heat, susceptibility) are strongly self-averaging, that is ${V}_{P}\ensuremath{\sim}{l}^{\ensuremath{-}d}$ (where
d=2 is the dimension). At criticality though, the results indicate that the magnetization
M and the susceptibility $\ensuremath{\chi}$ are non-self-averaging, i.e., $\frac{{V}_{\ensuremath{\chi}}}{{\ensuremath{\chi}}^{2}}, \frac{{V}_{M}}{{M}^{2}}\ensuremath{\nrightarrow}0$. The energy
E at criticality is clearly weakly self-averaging, that is ${V}_{E}\ensuremath{\sim}{l}^{\ensuremath{-}{y}_{v}}$ with
0<yv<d. Less conclusively, and possibly only as a transient behavior, the specific heat too is found to be weakly self-averaging. A phenomenological theory of finite size scaling for disordered systems is developed, based on physical considerations similar to those leading to the Harris criterion. Its main prediction is that when the specific heat exponent $\ensuremath{\alpha}<0$ ($\ensuremath{\alpha}$ of the disordered model) then, for a quantity
P which scales as ${l}^{\ensuremath{\rho}}$ at criticality, its variance
VP will scale asymptotically as ${l}^{2\ensuremath{\rho}+\frac{\ensuremath{\alpha}}{\ensuremath{\nu}}}$. The theory is not applicable in the asymptotic limit ($l\ensuremath{\rightarrow}\ensuremath{\infty}$) to the bond-disordered Ashkin-Teller model where $\frac{\ensuremath{\alpha}}{\ensuremath{\nu}}={0}_{+}$. Nonetheless in the accessible range of lattice sizes we found very good agreement between the theory and the data for ${V}_{\ensuremath{\chi}}$ and
VE. The theory may also be compatible with the data for the variance of the magnetization
VM and the variance of the specific heat
VC, but evidence for this is less convincing.
We simulated site dilute Ising models in
d=3 dimensions for several lattice sizes
L. For each
L singular thermodynamic quantities
X were measured at criticality and their distributions
P(X) were determined for ensembles of several thousand random samples. For $L\ensuremath{\rightarrow}\ensuremath{\infty}$ the relative width of
P(X) tends to a universal constant: there is no self-averaging. The width of the distribution of the sample
(i) dependent pseudocritical temperatures
Tc(i,L) scales as $\ensuremath{\delta}{T}_{c}(L)\ensuremath{\sim}{L}^{\ensuremath{-}1/\ensuremath{\nu}}$ and not as $\mathrm{\ensuremath{\sim}}{L}^{\ensuremath{-}d/2}$. The sample dependence of
Xi(T,L) enters dominantly, but not exclusively, via
Tc(i,L).