Search

Using a combination of analytic and numerical methods, we study the polarizability of a (non-interacting) Anderson insulator in one, two, and three dimensions and demonstrate that, in a wide range of parameters, it scales proportionally to the square of the localization length, contrary to earlier claims based on the effective-medium approximation. We further analyze the effect of electron-electron interactions on the dielectric constant in quasi-1D, quasi-2D and 3D materials with large localization length, including both Coulomb repulsion and phonon-mediated attraction. The phonon-mediated attraction (in the pseudogapped state on the insulating side of the Superconductor-Insulator Transition) produces a correction to the dielectric constant, which may be detected from a linear response of a dielectric constant to an external magnetic field., Comment: 9 pages

We combine numerical diagonalization with a semi-analytical calculations to prove the existence of the intermediate non-ergodic but delocalized phase in the Anderson model on disordered hierarchical lattices. We suggest a new generalized population dynamics that is able to detect the violation of ergodicity of the delocalized states within the Abou-Chakra, Anderson and Thouless recursive scheme. This result is supplemented by statistics of random wave functions extracted from exact diagonalization of the Anderson model on ensemble of disordered Random Regular Graphs (RRG) of N sites with the connectivity K=2. By extrapolation of the results of both approaches to N->infinity we obtain the fractal dimensions D_{1}(W) and D_{2}(W) as well as the population dynamic exponent D(W) with the accuracy sufficient to claim that they are non-trivial in the broad interval of disorder strength W_{E}10^{5} reveals a singularity in D_{1,2}(W)-dependencies which provides a clear evidence for the first order transition between the two delocalized phases on RRG at W_{E}\approx 10.0. We discuss the implications of these results for quantum and classical non-integrable and many-body systems., Comment: 4 pages paper with 5 figures + Supplementary Material with 5 figures

Motivated by the problem of Many-Body Localization and the recent numerical results for the level and eigenfunction statistics on the random regular graphs, a generalization of the Rosenzweig-Porter random matrix model is suggested that possesses two localization transitions as the parameter $\gamma$ of the model varies from 0 to $\infty$. One of them is the Anderson transition from the localized to the extended states that happens at $\gamma=2$. The other one at $\gamma=1$ is the transition from the extended non-ergodic (multifractal) states to the extended ergodic states similar to the eigenstates of the Gaussian Orthogonal Ensemble. We computed the two-level spectral correlation function, the spectrum of multifractality $f(\alpha)$ and the wave function overlap which all show the transitions at $\gamma=1$ and $\gamma=2$., Comment: 8 pages, 9 figures (main text) + 7 pages, 4 figures (supplementary materials)

We employed the method of virial expansion in order to compute the retarded density correlation function (generalized diffusion propagator) in the critical random matrix ensemble in the limit of strong multifractality. We found that the long-range nature of the Hamiltonian is a common root of both multifractality and Levy flights which show up in the power-law intermediate- and long-distance behavior, respectively, of the density correlation function. We review certain models of classical random walks on fractals and show the similarity of the density correlation function in them to that for the quantum problem described by the random critical long-range Hamiltonians., Comment: 12 pages, 9 figures, v2: several typos corrected; submitted to PRB

We study level number variance in a two-dimensional random matrix model characterized by a power-law decay of the matrix elements. The amplitude of the decay is controlled by the parameter b. We find analytically that at small values of b the level number variance behaves linearly, with the compressibility chi between 0 and 1, which is typical for critical systems. For large values of b, we derive that chi=0, as one would normally expect in the metallic phase. Using numerical simulations we determine the critical value of b at which the transition between these two phases occurs., Comment: 6 pages

We consider a simple model of quantum disorder in two dimensions, characterized by a long-range site-to-site hopping. The system undergoes a metal-insulator transition -- its eigenfunctions change from being extended to being localized. We demonstrate that at the point of the transition the eigenfunctions do not become fractal. Their density moments do not scale as a power of the system size. Instead, in one of the considered limits our result suggests a power of the logarithm of the system size. In this regard, the transition differs from a similar one in the one-dimensional version of the same system, as well as from the conventional Anderson transition in more than two dimensions., Comment: 4 pages, 2 figures

The dynamical scaling for statistics of critical multifractal eigenstates proposed by Chalker is analytically verified for the critical random matrix ensemble in the limit of strong multifractality controlled by the small parameter $b\ll 1$. The power law behavior of the quantum return probability $P_{N}(\tau)$ as a function of the matrix size $N$ or time $\tau$ is confirmed in the limits $\tau/N\rightarrow\infty$ and $N/\tau\rightarrow\infty$, respectively, and it is shown that the exponents characterizing these power laws are equal to each other up to the order $b^{2}$. The corresponding analytical expression for the fractal dimension $d_{2}$ is found., Comment: 4 pages, 1 figure

We develop a semi-quantitative theory of electron pairing and resulting superconductivity in bulk "poor conductors" in which Fermi energy $E_F$ is located in the region of localized states not so far from the Anderson mobility edge $E_c$. We review the existing theories and experimental data and argue that a large class of disordered films is described by this model. Our theoretical analysis is based on the analytical treatment of pairing correlations, described in the basis of the exact single-particle eigenstates of the 3D Anderson model, which we combine with numerical data on eigenfunction correlations. Fractal nature of critical wavefunction's correlations is shown to be crucial for the physics of these systems. We identify three distinct phases: 'critical' superconductive state formed at $E_F=E_c$, superconducting state with a strong pseudogap, realized due to pairing of weakly localized electrons and insulating state realized at $E_F$ still deeper inside localized band. The 'critical' superconducting phase is characterized by the enhancement of the transition temperature with respect to BCS result, by the inhomogeneous spatial distribution of superconductive order parameter and local density of states. The major new feature of the pseudo-gaped state is the presence of two independent energy scales: superconducting gap $\Delta$, that is due to many-body correlations and a new "pseudogap" energy scale $\Delta_P$ which characterizes typical binding energy of localized electron pairs and leads to the insulating behavior of the resistivity as a function of temperature above superconductive $T_c$. Two gap nature of the "pseudo-gaped superconductor" is shown to lead to a number of unusual physical properties., Comment: 110 pages, 39 figures. The revised version corrects a number of typos, adds references and discussion of recent results

We develop a supersymmetric virial expansion for two point correlation functions of almost diagonal Gaussian Random Matrix Ensembles (ADRMT) of the orthogonal symmetry. These ensembles have multiple applications in physics and can be used to study universal properties of time-reversal invariant disordered systems which are either insulators or close to the Anderson localization transition. We derive a two-level contribution to the correlation functions of the generic ADRMT and apply these results to the critical (multifractal) power law banded ADRMT. Analytical results are compared with numerical ones., Comment: 13 pages, 1 figure

We consider the correlation of two single-particle probability densities $|\Psi_{E}({\bf r})|^{2}$ at coinciding points ${\bf r}$ as a function of the energy separation $\omega=|E-E'|$ for disordered tight-binding lattice models (the Anderson models) and certain random matrix ensembles. We focus on the models in the parameter range where they are close but not exactly at the Anderson localization transition. We show that even far away from the critical point the eigenfunction correlation show the remnant of multifractality which is characteristic of the critical states. By a combination of the numerical results on the Anderson model and analytical and numerical results for the relevant random matrix theories we were able to identify the Gaussian random matrix ensembles that describe the multifractal features in the metal and insulator phases. In particular those random matrix ensembles describe new phenomena of eigenfunction correlation we discovered from simulations on the Anderson model. These are the eigenfunction mutual avoiding at large energy separations and the logarithmic enhancement of eigenfunction correlations at small energy separations in the two-dimensional (2D) and the three-dimensional (3D) Anderson insulator. For both phenomena a simple and general physical picture is suggested., Comment: 16 pages, 18 figures

Some properties of $d$-dimensional disordered models with long-range random hopping amplitudes are investigated numerically at criticality. We concentrate on the correlation dimension $d_2$ (for $d=2$) and the nearest level spacing distribution $P_c(s)$ (for $d=3$) in both the weak ($b^d \gg 1$) and the strong ($b^d \ll 1$) coupling regime, where the parameter $b^{-d}$ plays the role of the coupling constant of the model. It is found that (i) the extrapolated values of $d_2$ are of the form $d_2=c_db^d$ in the strong coupling limit and $d_2=d-a_d/b^d$ in the case of weak coupling, and (ii) $P_ (s)$ has the asymptotic form $P_c(s)\sim\exp (-A_ds^{\alpha})$ for $s\gg $, with the critical exponent $\alpha=2-a_d/b^d$ for $b^d \gg 1$ and $\alpha=1+c_d b^d$ for $b^d \ll 1$. In these cases the numerical coefficients $A_d$, $a_d$ and $c_d$ depend only on the dimensionality., Comment: 9 pages, 6 .eps figures, contribution to the Festschrift for Michael Schreiber's 50th birthday

The two-level correlation function $R_{d,\beta}(s)$ of $d$-dimensional disordered models ($d=1$, 2, and 3) with long-range random-hopping amplitudes is investigated numerically at criticality. We focus on models with orthogonal ($\beta=1$) or unitary ($\beta=2$) symmetry in the strong ($b^d \ll 1$) coupling regime, where the parameter $b^{-d}$ plays the role of the coupling constant of the model. It is found that $R_{d,\beta}(s)$ is of the form $R_{d,\beta}(s)=1+\delta(s)-F_{\beta}(s^{\beta}/b^{d\beta})$, where $F_{1}(x)=\text{erfc}(a_{d,\beta} x)$ and $F_{2}(x)=\exp (-a_{d,\beta} x^2)$, with $a_{d,\beta}$ being a numerical coefficient depending on the dimensionality and the universality class. Finally, the level number variance and the spectral compressibility are also considerded., Comment: RevTex4, 7 two-column pages, 6 .eps figures, 1 table. A new section devoted to the spectral compressibility added. To be published in Phys. Rev. B

The nearest level spacing distribution $P_c(s)$ of $d$-dimensional disordered models ($d=1$ and 2) with long-range random hopping amplitudes is investigated numerically at criticality. We focus on both the weak ($b^d \gg 1$) and the strong ($b^d \ll 1$) coupling regime, where the parameter $b^{-d}$ plays the role of the coupling constant of the model. It is found that $P_c(s)$ has the asymptotic form $P_c(s)\sim\exp [-A_ds^{\alpha}]$ for $s\gg 1$, with the critical exponent $\alpha=2-a_d/b^d$ in the weak coupling limit and $\alpha=1+c_d b^d$ in the case of strong coupling., Comment: EPL style, 7 pages, 4 .eps figures, to be published in Europhys. Lett

Finite-size effects in the generalized fractal dimensions $d_q$ are investigated numerically. We concentrate on a one-dimensional disordered model with long-range random hopping amplitudes in both the strong- and the weak-coupling regime. At the macroscopic limit, a linear dependence of $d_q$ on $q$ is found in both regimes for values of $q \alt 4g^{-1}$, where $g$ is the coupling constant of the model., Comment: RevTex4, 5 two-column pages, 5 .eps figures, to be published in Phys. Rev. B

The system size dependence of the multifractal spectrum $f(\alpha)$ and its singularity strength $\alpha$ is investigated numerically. We focus on one-dimensional (1D) and 2D disordered systems with long-range random hopping amplitudes in both the strong and the weak disorder regime. At the macroscopic limit, it is shown that $f(\alpha)$ is parabolic in the weak disorder regime. In the case of strong disorder, on the other hand, $f(\alpha)$ strongly deviates from parabolicity. Within our numerical uncertainties it has been found that all corrections to the parabolic form vanish at some finite value of the coupling strength., Comment: RevTex4, 6 two-column pages, 4 .eps figures, new results added, updated references, to be published in Phys. Rev. B

The system size dependence of the fluctuations in generalized inverse participation ratios (IPR's) $I_{\alpha}(q)$ at criticality is investigated numerically. The variances of the IPR logarithms are found to be scale-invariant at the macroscopic limit. The finite size corrections to the variances decay algebraically with nontrivial exponents, which depend on the Hamiltonian symmetry and the dimensionality. The large-$q$ dependence of the asymptotic values of the variances behaves as $q^2$ according to theoretical estimates. These results ensure the self-averaging of the corresponding generalized dimensions., Comment: RevTex4, 5 pages, 4 .eps figures, to be published in Phys. Rev. B

We investigate numerically the inverse participation ratio, $P_2$, of the 3D Anderson model and of the power-law random banded matrix (PRBM) model at criticality. We found that the variance of $\ln P_2$ scales with system size $L$ as $\sigma^2(L)=\sigma^2(\infty)-A L^{-D_2/2d}$, being $D_2$ the correlation dimension and $d$ the system dimension. Therefore the concept of a correlation dimension is well defined in the two models considered. The 3D Anderson transition and the PRBM transition for $b=0.3$ (see the text for the definition of $b$) are fairly similar with respect to all critical magnitudes studied., Comment: RevTex, 5 pages, 4 eps figures, to be published in Phys. Rev. Lett

We investigate numerically the power-law random matrix ensembles. Wavefunctions are fractal up to a characteristic length whose logarithm diverges asymmetrically with different exponents, 1 in the localized phase and 0.5 in the extended phase. The characteristic length is so anomalously large that for macroscopic samples there exists a finite critical region, in which this length is larger than the system size. The Green's functions decrease with distance as a power law with an exponent related to the correlation dimension., Comment: RevTex, 4 pages, 4 eps figures. Final version to be published in Phys. Rev. Lett

We perform a numerical simulation of energy relaxation in three-dimensional electron glasses in the strongly localized regime at finite temperatures. We consider systems with no interactions, with long-range Coulomb interactions and with short-range interactions, obtaining a power law relaxation with an exponent of 0.15, which is independent of the parameters of the problem and of the type of interaction. At very long times, we always find an exponential regime whose characteristic time strongly depends on temperature, system size, interaction type and localization radius. We extrapolate the longest relaxation time to macroscopic sizes and, for interacting samples, obtain values much larger than the measuring time. We finally study the number of electrons participating in the relaxation processes of very low energy configurations., Comment: 6 eps figures. To be published in Phys. Rev. B

The effect of Coulomb and short-range interactions on the spectral properties of two-dimensional disordered systems with two spinless fermions is investigated by numerical scaling techniques. The size independent universality of the critical nearest level-spacing distribution $P(s)$ allows one to find a delocalization transition at a critical disorder $W_{\rm c}$ for any non-zero value of the interaction strength. At the critical point the spacings distribution has a small-$s$ behavior $P_c(s)\propto s$, and a Poisson-like decay at large spacings., Comment: 4 two-column pages, 3 eps figures, RevTeX, new results added

We show by a numerical procedure that a short-range interaction $u$ induces extended two-particle states in a two-dimensional random potential. Our procedure treats the interaction as a perturbation and solve Dyson's equation exactly in the subspace of doubly occupied sites. We consider long bars of several widths and extract the macroscopic localization and correlation lengths by an scaling analysis of the renormalized decay length of the bars. For $u=1$, the critical disorder found is $W_{\rm c}=9.3\pm 0.2$, and the critical exponent $\nu=2.4\pm 0.5$. For two non-interacting particles we do not find any transition and the localization length is roughly half the one-particle value, as expected., Comment: 4 two-column pages, 4 eps figures, Revtex, to be published in Europhys. Lett

The elastic mean free path of carriers in a recently introduced model of quantum chaotic billiards in two and three dimensions is calculated. The model incorporates surface roughness at a microscopic scale by randomly choosing the atomic levels at the surface sites between -W/2 and W/2. Surface roughness yields a mean free path l that decreases as L/W^2 as W increases, L being the linear size of the system. But this diminution ceases when the surface layer begins to decouple from the bulk for large enough values of W, leaving more or less unperturbed states on the bulk. Consequently, the mean free path shows a minimum of about L/2 for W of the order of the band width. Energy fluctuations reflect the behavior of the mean free path. At small energy scales, strong level correlations manifest themselves by small values of the number of levels variance Sigma^2(E) that are close to Random Matrix Theory (RMT) in all cases. At larger energy scales, fluctuations are below the logarithmic behavior of RMT for l > L, and above RMT value when l < L., Comment: 8 twocolumn pages, seven figures, revtex and epsf macros. To be published in Physical Review B.

Quantum chaotic dynamics is obtained for a tight-binding model in which the energies of the atomic levels at the boundary sites are chosen at random. Results for the square lattice indicate that the energy spectrum shows a complex behavior with regions that obey the Wigner-Dyson statistics and localized and quasi-ideal states distributed according to Poisson statistics. Although the averaged spatial extension of the eigenstates in the present model scales with the size of the system as in the Gaussian Orthogonal Ensemble, the fluctuations are much larger., Comment: 4 twocolumnn pages in revtex style, 4 postscript figures, to be published in PRL, send comments to jav@boston.icmm.csic.es