---
abstract: 'In quantum systems, signatures of multifractality are rare. They have been found only in the multiscaling of eigenfunctions at critical points. Here we demonstrate multifractality in the magnetic-field-induced universal conductance fluctuations of the conductance in a quantum condensed-matter system, namely, high-mobility single-layer graphene field-effect transistors. This multifractality decreases as the temperature increases or as doping moves the system away from the Dirac point. Our measurements and analysis present evidence for an incipient Anderson-localization near the Dirac point as the most plausible cause for this multifractality. Our experiments suggest that multifractality in the scaling behaviour of local eigenfunctions are reflected in macroscopic transport coefficients. We conjecture that an incipient Anderson-localization transition may be the origin of this multifractality. It is possible that multifractality is ubiquitous in transport properties of low-dimensional systems. Indeed, our work suggests that we should look for multifractality in transport in other low-dimensional quantum condensed-matter systems.'
author:
- Kazi Rafsanjani Amin
- Samriddhi Sankar Ray
- Nairita Pal
- Rahul Pandit
- Aveek Bid
title: Exotic Multifractal Conductance Fluctuations in Graphene 
---

Introduction
============

One of the unsolved problems in single-layer graphene (SLG) is the nature of the electronic wave-function near the charge-neutrality (Dirac) point. In principle, the charge-carrier density of SLG should be continuously tunable, down to zero, leading to the largely unexplored regime of extremely weak interactions in a low-carrier-density system. The interaction parameter $r_{\mathrm{s}}$, which parametrizes the ratio of the average inter-electron Coulomb interaction energy to the Fermi energy, turns out to be independent of charge-carrier density in the case of SLG where $r_{\mathrm{s}} = e^2/(\kappa\hbar v_{\mathrm{F}})$. Here, $\kappa$ is the dielectric constant of the surrounding medium and $v_{\mathrm{F}}$ is the Fermi velocity. In the case of SLG on an SiO$_2$ substrate, $r_{\mathrm{s}}$ = 0.8. Thus SLG is a very weakly interacting system, when compared to other conventional two-dimensional systems such as GaAs/AlGaAs and Si inversion layers, where the values of $r_{\mathrm{s}}$ are typically much higher. A naïve application of the scaling theory to such a system in this regime would predict Anderson localization— a disorder-driven quantum phase transition, leading to a complete localization of the charge-carriers, and thence, to an insulator [@RevModPhys.80.1355; @PhysRevB.50.7526; @PhysRevB.74.041403; @PhysRevLett.97.236801; @Katsnelson2006]. Indeed, theoretical calculations for graphene indicate that intervalley scattering can lead to changes in the local and averaged electronic density of states (DOS) with the creation of localized states [@PhysRevLett.96.036801; @PhysRevLett.116.126804]. For example, graphene-terminated SiC (0001) surfaces undergo an Anderson-localization transition upon dosing with small amounts of atomic hydrogen [@PhysRevLett.103.056404]. This intervalley-scattering-induced localization is most effective near the Dirac point, where the screening of the impurity scatterers is negligible [@PhysRevLett.116.126804]. Experiments, however, find the appearance of a minimum conductance value at the Dirac point with $\sigma_{\textrm{min}}\approx 2e^2/\pi\hbar$, instead of a diverging resistance which is the hallmark of a truly localized state; notable exceptions being carrier localization [@Jung2011] and an Anderson-localization transition in bi-layer graphene heterostructures [@Ponomarenko2011].

By studying the scaling behaviour of the universal conductance fluctuations (UCF), we look for signs of charge-carrier localization near the Dirac point in ultra-high mobility SLG and uncover, as a result, an exotic multifractal behaviour in the UCF. Multifractality, characterized by an infinite number of scaling exponents, is ubiquitous in classical systems. Since the pioneering work of Mandelbrot [@mandelbrot19827he], the detection and analysis of multifractal scaling in such systems have enhanced our understanding of several complex phenomena, e.g., the dynamics of the human heart-beat [@Ivanov1999], the form of critical wave-functions at the Anderson-localization transition [@RevModPhys.80.1355], the time series of the Sun’s magnetic field [@lawrence1993multifractal], in medical-signal analysis (for instance, in pattern recognition, texture analysis and segmentation) [@6602370], fully-developed turbulence and in a variety of chaotic systems [@frisch1995turbulence; @0305-4470-17-18-021]. In condensed-matter systems, signatures of multifractality are usually sought in the scaling of eigenfunctions at critical points [@0305-4470-19-8-004; @JANSSEN19981; @Mirlin2000259; @PhysRevB.51.663; @PhysRevB.91.085427; @PhysRevLett.111.066601; @PhysRevB.87.125144; @2053-1583-1-1-011009]. Despite compelling theoretical predictions [@Facchini2007266; @PhysRevB.79.205120; @PhysRevLett.54.1718; @PhysRevB.43.3601; @PhysRevLett.72.3582; @PhysRevLett.81.4696], there are no reports of the successful observation of multifractality in transport coefficients.

Simple fractal conductance fluctuations, on the other hand, have been observed in several condensed-matter systems [@PhysRevLett.77.3885; @PhysRevLett.87.036802; @PhysRevLett.80.1948; @PhysRevB.54.10841; @PhysRevLett.77.3885; @PhysRevLett.87.036802; @1742-6596-109-1-012035]. These arise through semi-classical electron-wave interference processes whenever a system has mixed chaotic-regular dynamics [@MEISS1986387; @PhysRevLett.59.2503; @PhysRevB.54.10841; @PhysRevLett.77.3885]. In all these systems, a primary prerequisite for the observation of fractal transport is that the electron dynamics be amenable to semi-classical analysis - the fractal nature of conductance fluctuations seen in these systems disappears as the system is driven deep into the quantum limit. For such a mixed-phase semi-classical system, the graph of conductance ($G$) versus an externally applied magnetic field ($B$) has the same statistical properties as a Gaussian random process with increments of mean zero and variance $(\Delta B)^\gamma$. These processes are known as fractional Brownian motion and have the property that their graph is a fractal of dimension $D_{\mathrm{F}} = 2 - \frac{\gamma}{2}$  [@PhysRevLett.77.3885; @PhysRevB.54.10841], with $1\leq D_{\mathrm{F}}\leq 2$.

A semi-classical description, however, breaks down in the case of materials where the charge-carriers obey a Dirac dispersion relation, e.g., in graphene and topological insulators. In this report, we ask: can we find signatures of multifractality in a quantum system through transport measurements, specifically through the statistical properties of the graph of conductance fluctuations versus an external parameter like the magnetic field? We address this question by studying in detail the statistics of conductance fluctuations in high-mobility SLG-FET devices as a function of the perpendicular magnetic field over a wide range of temperature and doping levels. We report, for the first time, the occurrence of multifractality in the UCF in high mobility graphene devices deep in the quantum limit. Our measurements and analysis suggest at an incipient Anderson localization transition in graphene near the Dirac point.

Results
=======

![image](fig1_v2){width=".7\textwidth"}

Measurement of Universal Conductance Fluctuations. {#measurement-of-universal-conductance-fluctuations. .unnumbered}
--------------------------------------------------

SLG-FET devices with mobilities in the range $[20,000 - 30,000]$ cm$^2$V$^{-1}$s$^{-1}$ were fabricated on SiO$_2$ substrates by mechanical exfoliation from natural graphite, followed by conventional, electron-beam lithography [@novoselov2005two]\[\[fig:device\] (a-b)\]. We begin with the dependence of the resistance ($R$) on the gate voltage ($V_{\mathrm{G}}$), for the device G28M6. In \[fig:device\](c) we show plots of $R$ versus $\Delta V_{\mathrm{G}} = V_{\mathrm{G}} -V_{\mathrm{D}}$, where $V_{\mathrm{D}}$ is the Dirac point, measured at different temperatures ($T$). The high mobility and the position of the charge-neutrality point very close to $V_{\mathrm{G}}=0$ V attests to the high quality of the devices. We measured the magnetoconductance ($G$) as a function of the magnetic field $(\textbf{B}=(0,0,B) ) $, applied perpendicular to the plane of the device , in the range $-0.2$ T $\le$ $B$ $\le$ $0.2$ T. The presence of UCF was confirmed by the appearance of reproducible, non-periodic, but magnetic-field-symmetric, oscillations in $G$. The measurements were performed on multiple devices, over a wide range of $V_{\mathrm{G}}$ and $T$ . We find our UCF data to be in excellent agreement with previous studies of magnetoresistance and conductance fluctuations in single-layer graphene [@0953-8984-22-20-205301; @Berger1191; @Bohra2012; @PhysRevLett.100.056802; @Horsell20091041]. \[fig:ucf\_T\](a) shows illustrative plots of $G(B)$ with the Fermi energy ($E_{\mathrm{F}}$) maintained very close to the Dirac point \[$\Delta V_{\mathrm{G}}  ~\simeq 0$\] for the device G28M6. The data for other devices are similar and are shown in the Supplementary Figure 1 \[see Supplementary Note 1\]. At low values of $|B|$, near the minimum of $G$ at $B=0$, weak-localization corrections are visible \[\[fig:ucf\_T\](a)\]. As we move away from $B=0$, the conductance fluctuations become prominent. The amplitudes of the UCF peaks, and the values of the charge-carrier phase-coherence length $L_\upphi$, obtained from the variance of the UCF \[\[fig:ucf\_T\](b)\], decrease with increasing temperature because of thermal dephasing \[see Supplementary Note 2\]. The temperature dependence of the intervalley-scattering length and the intravalley-scattering length, extracted from weak-localization measurements at $T=$20 mK and $\Delta V_{\mathrm{G}}=$0.2 V, are shown in the Supplementary Figure 4 \[see Supplementary Note 3\]. We find them to be in excellent agreement with previous studies of localization in single-layer graphene [@PhysRevLett.98.176805; @PhysRevLett.100.056802]. We also observe an increase in $L_\upphi$ with increasing $|\Delta V_{\mathrm{G}}|=|V_{\mathrm{G}}-V_{\mathrm{D}}|$ \[\[fig:ucf\_T\](c)\], which we attribute to the increase in the screening of impurities by charge carriers [@0953-8984-22-20-205301]. We note an apparant saturation of $L_\upphi$ below a temperature of $\simeq100$ mK. The saturation of the phase-coherence length ($L_\upphi$) with decreasing temperature is an issue that has been at the forefront of research in several other semiconducting materials including doped Si [@PhysRevLett.91.066604; @PhysRevB.55.R13452; @PhysRevB.68.085413]. There are many effects, e.g. the presence of magnetic impurities [@PhysRevB.68.085413] or finite size, which can lead to a saturation of $L_\upphi$. We note here that a decoupling of the electron and lattice temperature leading to a saturation of the $L_\upphi$ is also possible. However, the data shown in \[fig:device\](c) and \[fig:ucf\_T\](a) show a continuous evolution of both the conductance and conductance fluctuations down to 20 mK. This rules out any saturation of the electron temperature down to 20 mK, and hence, the observed saturation of $L_\upphi$ below a certian temperature is not an experimental artifact \[see Supplementary Note 3\]. In a separate set of measurements, we obtain the magnitude of the UCF, at a given magnetic field, by sweeping over $V_{\mathrm{G}}$ and calculating the rms value of the fluctuations. We found that this quantity decreased sharply with increasing magnetic field \[see Supplementary Note 4, Supplementary Figure 5\] in conformity with theoretical predictions [@Liu2016].

![image](fig2_v2){width="\textwidth"}

![image](fig3){width=".75\textwidth"}

Analysis of Fractal scaling of the UCF. {#analysis-of-fractal-scaling-of-the-ucf. .unnumbered}
---------------------------------------

UCF represents quantum correction to Drude conductivity arising from the interference of electronic wavefunctions; and it is fingerprint of the disorder configuration in the conducting channel. Besides information about the phase-coherence of the charge carriers, it can also provide crucial insights into the electron dynamics and distribution of eigenstates through a scaling-dimension analysis of the magnetoconductance traces [@PhysRevB.54.10841; @PhysRevLett.77.3885; @PhysRevLett.87.036802]. We first compute the simple fractal dimensions $D_{\mathrm{F}}$ of the UCF curves via the Ketzmeric-variance method \[Supplementary Note 5\]. \[fig:lphi\_df\](a) shows plots of $D_{\mathrm{F}}$ versus $T$ for the device G28M6. At very low $T$ and small $|\Delta V_{\mathrm{G}}|$, we find $1<D_{\mathrm{F}}<2$. With increasing temperature, $D_{\mathrm{F}}$ $\to$1 monotonically. In this high-$T$ regime, the thermal-diffusion length $L_{\mathrm{T}} = \sqrt{\hbar D/k_{\mathrm{B}}T} \ll L_\upphi$, with $D$ the thermal-diffusion coefficient of the charge carriers; therefore, quasiparticle phase-decoherence, induced by inelastic thermal scattering, suppresses quantum interference. For large $|\Delta V_{\mathrm{G}}|$, the magnitude of the UCF is comparable to, or smaller than the background electrical noise, so $D_{\mathrm{F}}$$\to$$2$, the value for Gaussian white noise. In \[fig:lphi\_df\](b), we plot $D_{\mathrm{F}}$ versus $L_\upphi$ for two different devices: G28M6 and G30M4; remarkably, all the data points from these two devices cluster in the vicinity of a curve, with $D_{\mathrm{F}} \propto \ln (L_\upphi)$. In the limits $L_\upphi \ll L$, the UCF is non-fractal ($D_{\mathrm{F}}\simeq 1$), where as for large $L_\upphi$, the UCF is a fractal, so $1 < D_{\mathrm{F}} < 2$.

![image](fig4){width=".8\textwidth"}

Analysis of Multifractal scaling of the UCF. {#analysis-of-multifractal-scaling-of-the-ucf. .unnumbered}
--------------------------------------------

We build upon the predictions of multifractal scaling of conductance fluctuations in quantum systems [@Facchini2007266] by carrying out a multifractal detrended fluctuation analysis \[Supplementary Note 6\] of our UCF. The multifractality can be represented in the following two ways: (1) By the generalized Hurst exponent $h(q)$, defined using the order-$q$ moment of the UCF as $\langle {\rm rms}[\Delta G(\Delta B)]^q \rangle^{1/q} \sim [\Delta B]^{h(q)}$, and (2) by the multifractal spectrum $f(\alpha)$, obtained from the Legendre transform of $h(q)$. For a monofractal function, $h(q)$ has a single, $q$-independent value. For each one of our UCF plots, we obtain $h(q)$ in the range $-4 \le q \le 4$. In \[fig:mf\](b), we give an illustrative plot of $h(q)$ versus $q$ that we obtain from our magnetoconductance data at 20 mK \[\[fig:mf\](a)\]; $h(q)$ goes smoothly from $\simeq 1.9$, at $q=-4$, to $\simeq 0.85$, at $q=4$; the corresponding $f(\alpha)$ spectrum is plotted in \[fig:mf\](c). The singularity spectra have a definite maximum value of 1 (which is the dimension of the support graph). The width of the multifractal spectrum is defined as $\Delta \alpha  \equiv $ $h(q)_{\rm max}$ -$h(q)_{\rm min}$. The wide range of $h(q)$, or, equivalently, the wide spectrum ($\Delta \alpha$=1.05) quantifies the multifractality of the UCF. This is the first observation of multifractality of a conductance in any quantum-condensed-matter system and is the central result of our work.

We note that there are two distinct properties of UCF that can give rise to its multifractal behaviour [@PhysRevE.82.011136; @Zhou2009]: (i) a fat-tailed, non-Gaussian distribution of the UCF differences (as a function of $\delta B$) or (ii) long-range, in the magnetic field $B$, correlations of the fluctuations of $\delta g(B)$. We have verified that the distribution of our measured $\delta g(B)$ is not log-normal. Having ruled out (i), we now give a convenient test for (ii): we check for multifractality in a data set obtained from a random shuffling \[see Supplementary Note 7\] of the original sequence of $\delta g(B)$. If the long-range correlations (ii) exist, then signatures of multifractality must be absent in the reshuffled data. Indeed, in our experimental data, we observe a near-complete suppression of multifractality in the shuffled data set with $h_{\rm shuf}(q) \simeq 0.5$ for all values of $q$ and $\Delta \alpha = $0.05 \[Supplmentary Figure 9\]. Hence it is reasonable to infer that the multifractality in our UCF can be traced back to long-ranged correlations, that are otherwise difficult to measure.

\[fig:mf\](d) shows plots of $f(\alpha)$ for different values of $T$ and $\Delta V_G=0.2$V. The symmetry of $f(\alpha)$ about $\alpha_0$ \[\[fig:mf\](c)\] reflects the distribution of fluctuations, about the mean of the UCF. The large-fluctuation (small-fluctuation) segments contribute predominantly to the $q>0$ ($q<0$) part of $h(q)$. The $q>0$ ($q<0$) part of $h(q)$ maps onto the $\alpha < \alpha_0$ ( $\alpha > \alpha_0$) region of $f(\alpha)$, which is more-or-less symmetric about $\alpha_0$ at low $T$ \[\[fig:mf\](d)\]. As we increase $T$, this symmetry is lost. The magnitude of the skewness $\langle[\delta\alpha]^3\rangle/\langle[\delta\alpha]^2\rangle^{3/2}$, where $\delta\alpha = (\alpha - \alpha_0)$, increases with $T$. Therefore, as $T$ increases, large-amplitude conductance fluctuations become rarer than small-amplitude fluctuations. This is consistent with our plots of the UCF \[\[fig:ucf\_T\](a)\].

Discussion
==========

A naïve characterization of the fractal property of a curve, say by the measurement of one fractal dimension, does not rule out multifractality of this curve, which requires the calculation of an infinite number of dimensions [@PhysRevLett.50.346] (related to $h(q)$). One dimension suffices for monofractal scaling (as, e.g., in the scaling of velocity structure functions in the inverse-cascade region of forced, two-dimensional fluid turbulence [@Kellay2002; @Boffetta2012]). Our measurements of the UCF in single-layer graphene show that it is multifractal only if (i) the temperature is low and (ii) $\Delta V_{\mathrm{G}} \simeq 0$. If either one of these conditions is not met, the plot of $\delta g$ versus $B$ is a monofractal \[see Supplementary Note 8\]; at sufficiently large $T$, $D_{\mathrm{F}} \to$1 and the plots are non-fractal.

What can be the possible origin of our multifractal UCF? We list three potential causes: (1) scarring of wave functions (e.g., because of classically-chaotic billiards [@PhysRevLett.53.1515]); (2) quasi-periodicity in the Hamiltonian induced by a magnetic field [@0370-1298-68-10-305; @PhysRevB.14.2239; @PhysRevLett.50.1873] and its analogue for graphene [@0953-8984-22-46-465305]; (3) Anderson-localization-induced multifractality [@RevModPhys.80.1355]. We critically examine each one of these possibilities below and conclude that our results are most compatible with the last of these mechanisms.

While describing a quantum system, whose classical analogue is chaotic, one encounters scarred wavefunctions, whose intensity is enhanced along unstable, periodic orbits of the classical system. This non-uniform distribution of intensity of wavefunctions results from quasiparticle interference. Quantum scars can lead to pointer states [@PhysRevD.24.1516] with long trapping times, and, consequently, to large conductance fluctuations. Relativistic quantum scars have been predicted theoretically in geometrically confined graphene stadia, which exhibit classical chaos [@PhysRevLett.103.054101]. However, our devices are not shaped like billiards that are classically chaotic; and the charge carriers are not in the ballistic regime. Therefore, quantum scars cannot be the underlying cause for our multifractal UCF.

Fractal energy spectra can arise in tight-binding problems with an external magnetic field, which can be mapped onto Schrödinger problems with quasiperiodic potentials [@0370-1298-68-10-305; @PhysRevB.14.2239; @PhysRevLett.50.1873; @PhysRevB.29.1394]. It has been argued [@0953-8984-22-46-465305], therefore, that a fractal conductance can also arise, via Hofstadter-butterfly-type spectra, in Dirac systems at sufficiently high magnetic fields $B_{\mathrm{H}} \simeq \phi_0/A_0$ ($\simeq 10^5$ T for our samples), where $\phi_0 = 2 \pi\hbar/e$ is the flux quantum and $A_0$ is the unit-cell area. Our measurements use very-low-magnitude magnetic fields ($\lesssim 0.2$ T); this rules out a quasi-periodicity-induced multifractal UCF.

![image](fig5){width=".65\textwidth"}

The most compelling explanation of the multifractal UCF we observe in our SLG samples is an incipient Anderson localization near the charge-neutrality point. Multifractality of the local amplitudes of critical eigenstates near Anderson localization has been studied, theoretically, in several quantum-condensed-matter systems [@Mirlin2000259; @PhysRevLett.111.066601; @PhysRevB.87.125144; @2053-1583-1-1-011009]. The multifractality of the eigenstates near the critical point directly affects the two-particle correlation function through the generalized diffusion coefficient [@ANDP:ANDP2065080803; @0295-5075-4-12-003], which, in turn, affects the local current fluctuations in the system via the Kubo formula. It is not obvious that this must be reflected in the (macroscopic) conductance, or its moments; however, it is plausible that near the critical point, the UCF may inherit multifractal behaviour from its counterpart in the eigenfunctions [@PhysRevLett.87.014101]. Indeed, there are several theoretical predictions of multifractality in transport coefficients including conductance jumps near the percolation threshold in random-resistor networks [@PhysRevLett.54.1718; @PhysRevB.43.3601], conductance fluctuations in quantum-Hall transitions [@PhysRevLett.72.3582], and the temperature dependence of the peak height of the conductance at the Anderson-localization transition [@PhysRevLett.81.4696].

Spectroscopic studies on single-layer on-substrate graphene devices have revealed that the local potential fluctuations in this system are strongest when $E_F$ is close to the Dirac point [@PhysRevLett.116.126804; @Jung2011]. This leads to electronic states that are quasi-localized  [@PhysRevLett.104.096804; @PhysRevLett.99.236801; @RevModPhys.83.407]. Such quasi-localized states have a high inverse participation ratio [@PhysRevLett.96.036801] that can lead to the multifractality seen in our experiments. If this is true, then the multifractality in the UCF should be largest near the Dirac point; and then fall off on either side of it. In \[fig:width\] we show the dependence of $\Delta\alpha$ on $T$ and $\Delta V_{\mathrm{G}}$ for the device G28M6 \[data obtained for other devices are qualitatively similar\]. We observe that $\Delta\alpha$ is indeed largest near $\Delta V_{\mathrm{G}}=0$ and at low $T$, where the conductance of the device is of the order of $e^2/h$, and it sharply decreases as either $T$ or the magnitude of $\Delta V_{\mathrm{G}}$ increases. Similarly, as $T$ is increased, thermal scattering increases quasiparticle dephasing, and eventually at high $T$, when $L_\upphi \sim L_{\mathrm{T}} \ll L$, quantum-interference effects are masked. From our observation that a large multifractality arises only when quantum-interference-induced charge-carrier localization is significant, we propose that an incipient Anderson localization near the Dirac point is the most plausible origin of multifractal UCF in SLG.

This interpretation of the multifractality of the UCF in single-layer graphene devices is based on previous theoretical predictions and analysis. We summarize our argument below:

Conductance fluctuations, as a function of the magnetic field, have been shown to have a fractional fractal dimensions in some simple, one-dimensional, quantum systems, e.g., the kicked rotor [@PhysRevLett.87.014101; @PhysRevLett.84.63]. In these systems, such a fractional fractal dimension arises if one of the following conditions holds: (1) The PDF $P(t)$, of the charge-carrier survival time $t$, has a power-law form $P(t) \propto t^{-\gamma}$ at large $t$ (as opposed to an exponential decay); (2) the energy correlation $C(\Delta E)$ of elements of the S-matrix exhibit power laws (i.e., $C(\Delta E) \propto (\Delta E)^{-\gamma}$) [@PhysRevB.54.10841; @PhysRevLett.84.63; @PhysRevLett.77.3885; @PhysRevLett.68.3491]. Such survival probabilities are related to the conductance [@PhysRevLett.84.63; @PhysRevLett.87.014101]; and their multifractal behaviour has been explored [@Facchini2007266].

At the Anderson-localization transition, it is known that both the probability density function describing the diffusion of a wave-packet and the two-particle correlation function decay algebraically with a fractional power [@ANDP:ANDP2065080803; @0295-5075-4-12-003; @nakayama2013fractal]. Hence, as in the case of the simple quantum systems mentioned above, we may expect multifractal fluctuations of the conductance at the Anderson-localization transition. However, to the best of our knowledge, there are no exact, analytical results that yield a one-to-one correspondence between the multifractality of a critical wavefunction and the multifractality of the magnetoconductance. Thus, we propose that the multifractality of the critical wavefunctions at the Anderson localization is the most plausible cause of the multifractality of the UCF we have observed.

Conclusions
===========

In conclusion, we have uncovered and quantified the multifractal structure of mesoscopic conductance fluctuations in single-layer graphene devices. We speculate that our results are indicative of an incipient Anderson localization in this system. In particular, we quantify the multifractality of transport in a quantum condensed-matter system. There may well be multifractality in transport properties in systems other than graphene, and that multifractality is not unique to graphene. Our work provides a natural framework for studying the multifractality of such transport properties.

Acknowledgements
================

We acknowledge discussions with A. D. Mirlin and S. Bhattacharjee. A.B. acknowledges financial support from Nanomission, DST, Govt. of India project SR/NM/NS-35/2012; SERB, DST, Govt. of India and Indo-French Centre for the Promotion of Advanced Recearch (CEFIPRA). K.R.A. thanks CSIR, MHRD, Govt. of India for financial support. S.S.R. acknowledge financial support from the AIRBUS Group Corporate Foundation Chair in Mathematics of Complex Systems established in ICTS-TIFR, the support from DST, Govt. of India, project ECR/2015/000361, and the Indo-French Center for Applied Mathematics (IFCAM). R.P. acknowledges DST, Govt. of India for support. N.P. thanks UGC, Govt. of India for financial support.

Author contributions
====================

K.R.A and A.B. designed the experiment, fabricated the devices and carried out the measurements. K.R.A., A.B., N.P. and S.S.R. carried out the data analysis. All authors discussed the results and wrote the paper.

References {#references .unnumbered}
==========

[70]{}ifxundefined \[1\][ ifx[\#1]{} ]{}ifnum \[1\][ \#1firstoftwo secondoftwo ]{}ifx \[1\][ \#1firstoftwo secondoftwo ]{}““\#1””@noop \[0\][secondoftwo]{}sanitize@url \[0\][‘\
12‘\$12 ‘&12‘\#12‘12‘\_12‘%12]{}@startlink\[1\]@endlink\[0\]@bib@innerbibempty [****,  ()](\doibase 10.1103/RevModPhys.80.1355) [****,  ()](\doibase 10.1103/PhysRevB.50.7526) [****, ()](\doibase 10.1103/PhysRevB.74.041403) [****,  ()](\doibase 10.1103/PhysRevLett.97.236801) [****,  ()](\doibase 10.1038/nphys384) [****,  ()](\doibase 10.1103/PhysRevLett.96.036801) [****,  ()](\doibase 10.1103/PhysRevLett.116.126804) [****,  ()](\doibase
  10.1103/PhysRevLett.103.056404) [****,  ()](\doibase 10.1038/nphys1866) [****,  ()](\doibase
  10.1038/nphys2114) @noop [**]{} (, ) [****,  ()](http://dx.doi.org/10.1038/20924) @noop [****,  ()]{} in [**](\doibase 10.1109/WoSSPA.2013.6602370) () pp. @noop [**]{} (, ) [****,  ()](http://stacks.iop.org/0305-4470/17/i=18/a=021) [****,  ()](http://stacks.iop.org/0305-4470/19/i=8/a=004) [****, ()](\doibase
  http://dx.doi.org/10.1016/S0370-1573(97)00050-1) [****,  ()](\doibase http://dx.doi.org/10.1016/S0370-1573(99)00091-5) [****,  ()](\doibase 10.1103/PhysRevB.51.663) [****,  ()](\doibase 10.1103/PhysRevB.91.085427) [****,  ()](\doibase 10.1103/PhysRevLett.111.066601) [****,  ()](\doibase 10.1103/PhysRevB.87.125144) [****,  ()](http://stacks.iop.org/2053-1583/1/i=1/a=011009) [****,  ()](\doibase http://dx.doi.org/10.1016/j.physa.2006.10.012) [****, ()](\doibase 10.1103/PhysRevB.79.205120) [****,  ()](\doibase 10.1103/PhysRevLett.54.1718) [****,  ()](\doibase 10.1103/PhysRevB.43.3601) [****,  ()](\doibase 10.1103/PhysRevLett.72.3582) [****,  ()](\doibase 10.1103/PhysRevLett.81.4696) [****,  ()](\doibase
  10.1103/PhysRevLett.77.3885) [****,  ()](\doibase
  10.1103/PhysRevLett.87.036802) [****,  ()](\doibase
  10.1103/PhysRevLett.80.1948) [****,  ()](\doibase 10.1103/PhysRevB.54.10841) [****,  ()](http://stacks.iop.org/1742-6596/109/i=1/a=012035) [****,  ()](\doibase http://dx.doi.org/10.1016/0167-2789(86)90041-2) [****,  ()](\doibase 10.1103/PhysRevLett.59.2503) [****,  ()](\doibase
  10.1038/nature04233) [****,  ()](http://stacks.iop.org/0953-8984/22/i=20/a=205301) [****,  ()](\doibase
  10.1126/science.1125925),  [****,  ()](\doibase 10.1063/1.4748167),  [****,  ()](\doibase 10.1103/PhysRevLett.100.056802) [****, ()](\doibase
  http://dx.doi.org/10.1016/j.ssc.2009.02.058),  [****,  ()](\doibase 10.1103/PhysRevLett.98.176805) [****,  ()](\doibase 10.1103/PhysRevLett.91.066604) [****, ()](\doibase 10.1103/PhysRevB.55.R13452) [****,  ()](\doibase
  10.1103/PhysRevB.68.085413) [****,  ()](http://stacks.iop.org/0953-8984/28/i=13/a=135302) [****,  ()](\doibase 10.1103/PhysRevE.82.011136) [****,  ()](http://stacks.iop.org/0295-5075/88/i=2/a=28004) [****,  ()](\doibase 10.1103/PhysRevLett.50.346) [****,  ()](http://stacks.iop.org/0034-4885/65/i=5/a=204) [****,  ()](\doibase 10.1146/annurev-fluid-120710-101240), [****,  ()](\doibase 10.1103/PhysRevLett.53.1515) [****,  ()](http://stacks.iop.org/0370-1298/68/i=10/a=305) [****,  ()](\doibase 10.1103/PhysRevB.14.2239) [****,  ()](\doibase
  10.1103/PhysRevLett.50.1873) [****,  ()](http://stacks.iop.org/0953-8984/22/i=46/a=465305) [****, ()](\doibase 10.1103/PhysRevD.24.1516) [****,  ()](\doibase
  10.1103/PhysRevLett.103.054101) [****, ()](\doibase 10.1103/PhysRevB.29.1394) [****,  ()](\doibase 10.1002/andp.2065080803) [****,  ()](http://stacks.iop.org/0295-5075/4/i=12/a=003) [****,  ()](\doibase 10.1103/PhysRevLett.87.014101) [****, ()](\doibase 10.1103/PhysRevLett.104.096804) [****,  ()](\doibase 10.1103/PhysRevLett.99.236801) [****,  ()](\doibase 10.1103/RevModPhys.83.407) [****,  ()](\doibase 10.1103/PhysRevLett.84.63) [****,  ()](\doibase
  10.1103/PhysRevLett.68.3491) [**](https://books.google.co.in/books?id=D2r7CAAAQBAJ), Springer Series in Solid-State Sciences (, )

UCF data
========

![image](fig_S0_ucf){width="85.00000%"}

\[fig:df\_ucfplot\] shows plots of the magnetoconductance at different temperatures ($T$), for our device G30M4. The data have been vertically shifted for clarity. The measurements were carried out at a fixed gate voltage $\Delta V_{\mathrm{G}} = -0.7$ V. Each of the magnetoconductance plots show aperiodic oscillations, which are symmetric in magnetic field ($B$). The amplitude of the oscillations gets suppressed with increasing $T$, however, all the oscillation features (peaks and dips in the magnetoconductance) appear at exact same value of $B$. These symmetric and reproducible oscillations are the signatures of UCF.

We first symmetrize each of the the magnetoconductance $$G_{\rm sym}(B) = \frac{G(B) + G(-B)}{2},$$ to remove the anti-symmetric noise contribution from our data. This symmetric component, henceforth referred to as $G(B)$, is used for all our subsequent analysis.

Phase Coherence Length
======================

We estimate the phase-coherence length ($L_\upphi$) by two standard methods, other than obtaining from the UCF variance [@SBeenakker19911], which provid, qualitatively, similar result \[$L_\upphi$ dependence of $D_{\mathrm{F}}$, shown in Fig. 3(b) of the manuscript\]. First, we extract $L_\upphi$ by fitting the weak-localization peak with the Hikami-Larkin-Nagaoka equation [@SPhysRevLett.100.056802], $$\Delta \sigma = \frac{e^2}{\pi h} \times \bigg[ F \bigg(\frac{8 \pi B}{\Phi_0 L_\upphi^{-2}} \bigg) - F \bigg( \frac{8 \pi B}{\Phi_0 \big[L_\upphi^{-2} + 2L_{\mathrm{i}}^{-2} \big] } \bigg) -2F \bigg( \frac{8 \pi B}{\Phi_0 \big[ L_\upphi^{-2} + L_{\mathrm{i}}^{-2} + L_\ast^{-2} \big] } \bigg)  \bigg]
\label{eqn:wl}$$ where $F(z) = \ln z + \psi(0.5 + z^{-1})$, and $\psi(x)$ is the digamma function, $L_{\mathrm{i}}$ and $L_\ast$ are the characteristic length scale for inter-valley scattering and intra-valley scattering, respectively. The first term in Eqn. \[eqn:wl\] gives the weak-localization correction, and the other terms are attributed to weak anti-localization.

![image](fig_S_wlfit){width="55.00000%"}

![image](fig_S_autocorr){width="52.50000%"}

       **Method**       UCF variance   HLN fit   autocorrelation of UCF
  -------------------- -------------- --------- ------------------------
   $L_\upphi$($\mu$m)       1.20        1.27              0.96

  :  Comparison of $L_\upphi$ extracted via different methods\[tbl:lphi\]

\[fig:wlfit\] shows plots of a typical $\Delta\sigma=\sigma(B)-\sigma(0)$ versus $B$, and the fit to the data with Eqn. \[eqn:wl\]. The extracted value of $L_\upphi$ from the fit is $1.27~\mu$m (\[tbl:lphi\]).

$L_\upphi$ was also estimated from the auto-correlation ($C(\Delta B)$) of the UCF [@SBeenakker19911; @SHorsell20091041]. The autocorrelation of the magnetoconductance is defined as $$C(\Delta B) = \int dB G(B) G(B+\Delta B) .$$

The correlation magnetic field $B_{\mathrm{C}}$ is then obtained from the autocorrelation $C(\Delta B)$ using the relation, $$C(B_{\mathrm{C}}) = \frac{1}{2} C(0).$$ $L_\upphi$, is then calculated from $B_{\mathrm{C}}$ using the relation $$L_\upphi = \sqrt{h/eB_{\mathrm{C}}}.$$

\[fig:autocorr\] shows plot of normalized autocorrelation $C(\Delta B)/C(0)$ versus $\Delta B$ for the same UCF data \[for which the fit with the HLN equation is shown in \[fig:wlfit\]\].

A comparison of values of $L_\upphi$ extracted from a UCF data, measured at $T=0.2$ mK, using the above-mentioned different methods are shown in \[tbl:lphi\].

Effect of Temperature on UCF
============================

Plot of the rms value of the magnetoconductance fluctuation ($\delta G_{\rm rms}$) versus $T$ is shown in  \[fig:lphi\_vg\](a). The UCF corresponding to this data were measured at $\Delta V_{\mathrm{G}} = $0.2 V \[UCF series shown in Fig. 2(a) of main text\]. We observe a strong suppression of the magnetoconductance fluctuations with increasing temperature which is directly attributed to charge carrier dephasing induced by thermal scattering [@S0953-8984-22-20-205301; @SBerger1191; @SBohra2012; @SPhysRevLett.100.056802; @SHorsell20091041].

![image](Fig_S_lphi_vg_2){width="75.00000%"}

\[fig:lphi\_vg\](b) shows plots of $L_\upphi$ (in blue filled circles), $L_{\mathrm{i}}$ (in yellow filled squares) and $L_\ast$ (in red asterisk) versus $T$, extracted from fitting the weak localization correction in the magnetoresistance data with the HLN equation (Eqn. \[eqn:wl\]), at $T=20$ mK and $\Delta V_{\mathrm{G}}=0.2$ V. We find that the goodness of the fit is strongly influenced by the presence of UCF, as was observed by other groups as well [@SPhysRevLett.98.176805; @SPhysRevLett.100.056802].

We observe some evidence of a labelling off of $L_\upphi(T)$ around 100 mK. This is an issue that has been at the forefront of research in several other semiconducting materials including doped Si [@SPhysRevLett.91.066604; @SPhysRevB.55.R13452; @SPhysRevB.68.085413]. Saturation of $L_\upphi(T)$ have been observed in graphene at a much higher temperatures ($\sim 4$ K) and had been argued to arise from finite-size effects [@SPhysRevLett.100.056802]. The presence of disorder also enhances the saturation of $L_\upphi$. This is because the presence of disorder leads to the formation of electron-hole puddles, which modify the conducting paths and decrease the effective size of the device [@PhysRevLett.100.056802]. This, in turn, results in a size-limited saturation of $L_\upphi$. Earlier experiments have also shown that, scattering from extremely low concentration of magnetic impurity, even at a level undetectable by other means, can lead to the widely observed saturation of $L_\upphi$ (or equivalently, the dephasing time $\tau_\upphi$ ) in weakly disordered systems [@SPhysRevB.68.085413]. Our device has a higher mobility and a longer channel length than in Ref [@SPhysRevLett.100.056802]. Thus we expect the effects of disorder and small channel length, which typically leads to a saturation of $L_\upphi$, to be much weaker in our device. Hence, it is reasonable to observe a stronger temperature dependence of $L_\upphi$ down to lower temperatures, and a probable saturation only below 0.1 K.

 UCF at different magnetic fields
=================================

We have also measured UCF at a fixed magnetic field, by varying gate-voltage. We compare the rms value of the magnetoconductance fluctuations ($\delta G_{\rm rms}$) measured in two different ways: (1) at a constant magnetic field by varying gate voltage, and (2) at a constant gate voltage (with $\Delta V_{\mathrm{G}} \sim$ 0) by varying the magnetic field. We observe that the magnitude of the UCF measured during the gate sweep is $\sim$ 1.6 times larger than the UCF measured during the field sweep, which indicates non-ergodicity of UCF [@SLiu2016].

![image](Fig_S_var_g_b){width="75.00000%"}

Furthermore, we have carried out UCF measurements by sweeping the gate voltage at different values of the magnetic field in the range of 0 T to 0.2 T. In \[fig:rms\_b\], we plot the rms value of conductance fluctuations, $\delta G_{\rm rms}$ [versus $B$]{}. We find that $\delta G_{\rm rms}$ decreases with increasing $B$. These observations are in excellent agreement with predictions of behaviour of UCF in the presence of long-range component of random scattering potential [@SLiu2016].

Calculation of fractal dimensions
=================================

![image](fig_S1_scaling){width="70.00000%"}

We calculate the fractal dimensions, from our UCF data, in the following way.

We begin by defining and calculating $\delta G(B)$ = $G(B) - \langle G(B) \rangle$, where the angular brackets define the average of over $B$. To calculate the fractal dimension, we begin [@SPhysRevB.54.10841] by dividing our UCF data into $N_s$ overlapping segments of size $\Delta B$. We then define $\Delta G(\Delta B) \equiv \delta G(B) - \delta G(B+\Delta B)$ from which we calculate $\langle (\Delta G(\Delta B))^2 \rangle $, where we average over all the $N_s$ segments. Thence we evaluate $\langle (\Delta G(\Delta B))^2 \rangle $ for different values of $\Delta B$, and study the scaling. Typically, $\langle (\Delta G(\Delta B))^2 \rangle ~\sim ~(\Delta B)^\gamma$. From fits to log-log plots of $\langle (\Delta G(\Delta B))^2 \rangle$ versus $\Delta B$ (see \[fig:df\_scaling\] for a representative example), we obtain from our data the value of $\gamma$, which yields the fractal dimension $D_{\mathrm{F}}$ via $D_{\mathrm{F}} = 2 - \gamma /2$; for the data shown in \[fig:df\_scaling\], this yields $D_{\mathrm{F}} \simeq 1.5$.

Calculating multifractal exponents
==================================

![image](figure_S2){width="100.00000%"}

We compute the multifractal singularity spectrum form our UCF by using the method of Multifractal Detrended Fluctuation Analysis (MFDFA) [@Skantelhardt2012fractal], which we describe below.

We begin by a [*local detrending*]{} of our data. To do this, we choose a segment of the magnetic field $\Delta B$ and then divide the UCF into different overlapping segments. We then fit the different segments of our conductance data ${G}$ with a polynomial $P_m$ and remove the local trend: $$g(i) = G(i) - P_m(i).$$ It is important to keep the same order $m$ of the polynomial throughout this analysis and for all the segments. Alternatively, we can divide the conductance data ${G}$, with $N$ data points, into $N_s$ segments: $N_s = int(N/s)$, with $s$ points in each segment. In our analysis, the detrending is done with a polynomial of order 1.

Next, we evaluate the root-mean-square (rms) fluctuations and the order-$q$ moment of these fluctuation via $$g_{\rm rms} = \left (\frac{1}{s} \sum_{i=1}^s g(i)^2\right )^{1/2}.$$

We then define $$F_q(s) = \left (\frac{1}{N_s} \sum_{j=i}^{N_s} g_{\rm rms}(j)^q\right )^{1/q}. 
\label{eqn:fqs}$$ In \[fig:fq\](a-c) we show representative (log-log) plots of $F_q$ versus $s$ for three different temperatures and for different values of $q$. From such plots it is possible to calculate the generalized order-$q$ scaling exponent through the scaling relation $$F_q(s) \sim s^{H(q)}. \label{eqn:Hq}$$ We calculate the order-$q$ scaling exponent $H(q)$ for different values of $q$ and at different temperatures. In \[fig:fq\](d) we show a representative plot of $H(q)$ versus $q$ from data at the three temperatures shown in \[fig:fq\](a-c).

![image](Fig_S3_exponents){width=".75\textwidth"}

We now calculate the singularity spectrum as follows: We first calculate $$\tau(q) =  H(q) q - 1.$$

From this, we calculate the generalized Hurst exponent $$h(q) = \frac{d \tau(q)}{ dq}.  
\label{eqn_hq}$$ In \[fig:mf\_exp\](a) we show plot of $\tau (q)$, for a representative UCF. By using this, and making use of the relation between $\tau(q)$ and $h(q)$ shown above, we obtain plots of $h(q)$ such as the one shown in \[fig:mf\_exp\](b).

$\alpha$ and $f(\alpha)$ follow from $\tau(q)$ and $h(q)$: $$\alpha = \frac{d \tau(q)}{dq} = h(q);$$ $$f(\alpha) = \alpha q -  \tau(q).$$

In \[fig:mf\_exp\](c) we show a plot of the multifractal singularity spectrum $f(\alpha)$, which shows a clear maximum at $f(\alpha_0) = 1$. It is important to keep in mind that the position of this maxima $\alpha_0 \equiv h(0)$ depends on both the gate voltage $V_{\mathrm{G}}$ and the temperature $T$.

The error in $\tau(q)$ was estimated from the fitting error of the plot of $F_q(s)$ versus $s$ \[Eqn. \[eqn:Hq\]\]. The error in $h(q)$ was obtained from error in $H(q)$ by using Eqn. \[eqn\_hq\]. Absolute errors in $\tau(q)$ and $h(q)$ are represented using error-bars in \[fig:mf\_exp\](a) and (b), respectively. Finally, the error in the spectral width $\Delta \alpha$ is obtained by adding the error in $h(q)$, obtained at the two end values of q, *i.e.* $q=-4$ and $q=4$.

Multifractality of a randomly shuffled series
=============================================

We generate a randomly shuffled sequence following the prescription mentioned below:

\(i) We consider a data sequence $\{ g(i) \}$, $i=1,2,...,N$.

\(ii) We then generate a sequence $\{ I(i) \}$ of length $N$, which is a random permutation of the positive integers from $1$ to $N$.

\(iii) We then generate a new sequence $\{ g_{\rm rs}(i) \}$, where $g_{\rm rs}(i) = g(I(i))$.

\(iv) We repeat this shuffling for several times, to obtain the final randomly shuffled sequence.

\[fig:mf\_shuffle\] shows plots the generalized Hurst exponent $h(q)$ and the singularity spectrum $f(\alpha)$ of a UCF series (measured at $T$=20 mK and $\Delta V_{\mathrm{G}} = 0.2~$V) and its randomly shuffled counterpart. The spectrum for the original data series is shown in blue continuous line, where as that for the randomly shuffled series is shown in red dashed line. For the original data series, $h(q)$ has a large range of values and the singularity spectrum $f(\alpha)$ is wide ($\Delta \alpha= 1.05$ ) showing that the data series is multifractal. On the other hand, for the shuffled series, $h(q)$ is almost $q$-independent and $f(\alpha)$ is extremely narrow ($\Delta \alpha=0.05$), indicating a clear suppression of multifractality.

![image](mf_shuffled_data){width=".75\textwidth"}

Monofractality versus Multifractality
=====================================

![image](Fig_rebuttal_multifractal){width="65.00000%"}

A näive characterization of the fractal property of a curve, say by the measurement of one fractal dimension, does not rule out multifractality of this curve, which requires the calculation of an infinite number of dimensions [@SPhysRevLett.50.346] (related to our $h(q)$). For monofractal scaling, a single dimension suffices (as, *e.g.*, in the scaling of velocity structure functions in the inverse-cascade region of forced, two-dimensional fluid turbulence [@SKellay2002; @SBoffetta2012]).

We note here, that for the single-layer graphene devices studied by us, even at the lowest temperature, we find $f(\alpha)$ depends on $\Delta V_{\mathrm{G}}$. All the UCFs measured by us are not multifractal; they are multifractal only in a narrow range of gate-voltages around the Dirac point. We show in \[fig:notallfractalmultifractal2\] an example of a UCF, that is monofractal, but not multifractal. A UCF measured near the Dirac point ($\Delta V_{\mathrm{G}} =0.4$ V) \[\[fig:notallfractalmultifractal2\](a)\] is multifractal, as seen in \[fig:notallfractalmultifractal2\](c-d) (black solid lines). On the other hand, UCF measured away from the Dirac point ($\Delta V_{\mathrm{G}} =3.4$ V) \[\[fig:notallfractalmultifractal2\](b)\], despite having a fractal dimension $D_{\mathrm{F}} = 1.41$, is barely multifractal, as seen from the *almost* $q$ independent $h(q)$ and an extremely narrow singularity spectrum \[\[fig:notallfractalmultifractal2\](c-d) (blue dashed lines)\].

References {#references-1 .unnumbered}
==========

[16]{}ifxundefined \[1\][ ifx[\#1]{} ]{}ifnum \[1\][ \#1firstoftwo secondoftwo ]{}ifx \[1\][ \#1firstoftwo secondoftwo ]{}““\#1””@noop \[0\][secondoftwo]{}sanitize@url \[0\][‘\
12‘\$12 ‘&12‘\#12‘12‘\_12‘%12]{}@startlink\[1\]@endlink\[0\]@bib@innerbibempty in [**](\doibase
  http://dx.doi.org/10.1016/S0081-1947(08)60091-0), , Vol. ,  (, ) pp. [****,  ()](\doibase 10.1103/PhysRevLett.100.056802) [****, ()](\doibase
  http://dx.doi.org/10.1016/j.ssc.2009.02.058),  [****,  ()](http://stacks.iop.org/0953-8984/22/i=20/a=205301) [****,  ()](\doibase
  10.1126/science.1125925),  [****,  ()](\doibase 10.1063/1.4748167),  [****,  ()](\doibase 10.1103/PhysRevLett.98.176805) [****,  ()](\doibase 10.1103/PhysRevLett.91.066604) [****, ()](\doibase 10.1103/PhysRevB.55.R13452) [****,  ()](\doibase
  10.1103/PhysRevB.68.085413) [****,  ()](http://stacks.iop.org/0953-8984/28/i=13/a=135302) [****,  ()](\doibase 10.1103/PhysRevB.54.10841) in @noop [**]{} (, ) pp.  [****,  ()](\doibase 10.1103/PhysRevLett.50.346) [****,  ()](http://stacks.iop.org/0034-4885/65/i=5/a=204) [****,  ()](\doibase 10.1146/annurev-fluid-120710-101240),
