Quantum Diffusion in the Lowest Landau Level of Disordered Graphene

Electronic transport in the lowest Landau level of disordered graphene sheets placed in a homogeneous perpendicular magnetic field is a long-standing and cumbersome problem which defies a conclusive solution for several years. Because the modeled system lacks an intrinsic small parameter, the theoretical picture is infested with singularities and anomalies. We propose an analytical approach to the conductivity based on the analysis of the diffusive processes, and we calculate the density of states, the diffusion coefficient and the static conductivity. The obtained results are not only interesting from the purely theoretical point of view but have a practical significance as well, especially for the development of the novel high-precision calibration devices.


Introduction
Two-dimensional (2d) electronic systems in general and their transport properties in particular have been in the focus of intense research for several decades. In such systems, the effects due to quantum interference are strong and give rise to the interesting and rather unintuitive phenomena, as for instance various facets of the quantum Hall effect. Yet another effect on the transport that is supposed to be strong in 2d arises from the disorder which is always present in realistic materials. In conventional 2d electron systems, which are characterized by a parabolic and isomorphic spectrum, the presence of the disorder is widely believed to lead to the destructive interference of electronic quantum waves and consequently to the suppression of the electronic transport through the system on macroscopic scales. This phenomenon is usually called the Anderson localization of electronic wave functions, and it has received much of attention in the past [1][2][3][4]. This picture was challenged with the discovery of the unconventional behavior of electrons in the transition between Hall plateaux in quantum Hall systems. The experimental evidence from this observations points to the principal possibility for the existence of a metallic state in 2d systems under special conditions [5]. However, a real change of paradigm occured with the discovery of metallic states in graphene [6][7][8] and in a number of further low-dimensional systems, which is collectively known as the topological insulators [9][10][11][12][13][14][15][16]. A feature common to all these systems is the presence of the so-called nodes in the band structure and the linearity of the spectrum in the vicinity of these nodes. Despite being pristine 2d systems, they reveal a finite dc conductivity, which is very robust against disorder and thermal fluctuations.
The theoretical approach to the electronic transport of disordered electron gases is a rather formidable and cumbersome undertaking. Because the translational symmetry in the system is explicitly broken by the randomness, the usual methods of the theoretical analysis, which are mainly built around the duality between the position and momentum space representations and the special role of the Fourier transformation as the diagonalization tool for the quadratic Hamiltonians, no longer work. Therefore, the main idea behind every analytical approach to the macroscopic disordered systems is to reintroduce the translational invariance into the system by mapping the initial problem, which usually neglects the electron-electron and electron-phonon interactions from the outset, by a kind of an averaging operation on an effective interacting model in which the scattering of individual electrons on the randomly distributed potentials is approximated by the interaction operators expressed in terms of bilineals of second quantization operators. However, in practical terms, such an averaging procedure works well only under a weak disorder assumption, which guarantees a well-formed saddle-like shape of the free energy functional. In this case, the main effects caused by the disorder are taken into account by the summation of all contributions in partial diagrammatic channels [2][3][4][17][18][19][20][21][22][23][24][25][26][27][28][29][30].
In magnetic fields, the quantum mechanics of charge carriers with a linear spectrum specific for graphene is governed by an interplay of the intrinsic and magnetic-field-induced Berry curvatures. Several aspects of this physics remain widely untouched, though. For instance, relatively little is known about the role of disorder and its interplay with the magnetic field. The overall progress in this area has been slow because of the technical challenges, which are considerable even by the standards of the community [31][32][33][34][35]. A number of issues make the disordered electrons in the homogeneous perpendicular magnetic field look differently than the situation without a magnetic field. Due to the freedom of the gauge choice, the problem can be approached in a number of ways, which differ very much in details and in the outer appearance. The popular choice of the central gauge has the advantage that the solutions of the Schrödinger equation are states localized in the position space. Therefore, one can do computations in the position space in an exact manner.
The envisaged problem is notoriously difficult because the model lacks a small expansion parameter [36]. This inevitably leads to divergent expansion series. A powerful method developed to keep such divergences under control is the renormalization group. In the past, our understanding of the physics of disordered metals and semiconductors profited vastly from the various combinations of variational and perturbative techniques with the renormalization group, c.f. Refs. [1][2][3][25][26][27] and Refs. [37][38][39]. However, in the central gauge picture, there is no continuous variable to be sliced off by iterations in order to obtain the renormalization group equations. Of course, one can use a different gauge, which allows for a description in terms of states localized in one direction and propagating in the other. The price to pay is the loss of exactness, which is too costly to give up. In this paper, we develop a diagrammatic approach to the conductivity of the two-dimensional disordered electron gas in a strong magnetic field in a central gauge picture. While these series can still be wrapped up exactly for the single-particle propagators, as it was impressively demonstrated by Wegner in Ref. [32], additional technical issues make elusive every attempt of applying these techniques with the same success to the two-particles propagators. The available divergent series cannot be directly plugged into the Kubo formula without some not a priori obvious regularization or resummation. Hence, the usual way to approach the conductivity is via the Einstein relation and correspondingly via the notion of diffusion [40][41][42]. Because the corresponding statistical averages require normalization with respect to the vacuum fluctuations [43], this provides a tool of estimating the measurable quantities by means of some kind of analytical continuation [34,35,44,45].
To make our approach function, it relies on the information from the perturbative expansion. Therefore, we perform the exact computations of the perturbative series to the very high order. We identify the exact asymptotics of the two-particles propagator functions and approach the diffusion coefficient via the mean squared displacement using these asymptotics. It turns out that the behavior at longer time is dominated by the higher-order elements and tends toward a stationary state. On the sublaying time scales though, there is a large region with linear time dependence, which is characteristic of the diffusion. To approach this regime, we propose a self-consistent equation of motion for the mean squared displacement and extract the diffusion coefficient from there. With the obtained diffusion coefficient and density of states, we find via the Einstein relation a universal expression for the static conductivity in the lowest Landau level. All the system becomes metallic within a parameter window around the eigenvalues of the Hamiltonian of the clean system. With increasing disorder, this parameter window becomes broader. Numerically, the conductivity of disordered gapless and undoped graphene is very close to the experimentally established values.
The structure of this paper is as follows: In Section 2, we briefly discuss the main facts about the tight-binding Hamiltonian on the honeycomb lattice, its eigenvalues and eigenstates, and introduce the effective continuous model. In Section 3, we elaborate on the topological properties of the Hamiltonian and its eigenstates. In Section 4, we proceed with the consideration of the effective Hamiltonian, which describes the graphene in a strong external magnetic field and evaluates the single-particle propagator of the clean system in Section 5. In Section 6, we evaluate the Kubo-Greenwood formula for the dc conductivity of the clean gapless and chemically neutral graphene off and in an external magnetic field. In Section 7, we approach the single-particle propagator of the disordered system and discuss the Wegner's exact solution and the exact density of states. In Section 8, we give our result for the two-particles propagator and for the mean squared displacement. Finally, in Section 9, we extract the diffusion coefficient from the equation of motion of the mean squared displacement and with that the static conductivity.

Tight-Binding and Effective Hamiltonian of Graphene
First, we briefly review the main spectral and topological properties of the tightbinding Hamiltonian on a honeycomb lattice. In second quantization, it reads where c and d (c † and d † ) denote the annihilation (creation) operators acting on the lattice sites of each sublattice of the honeycomb lattice, respectively. The nearest neighbor positions are a 1 = a(0, −1), a 2,3 = a 2 ± √ 3, 1 , where a denotes the lattice spacing. The tightbinding Hamiltonian Equation (1) is translationally invariant and is diagonalized by a Fourier transform, giving the eigenvalues E ± = ±E = ± h 2 1 + h 2 2 and the respective eigenstates of the first-quantized Hamiltonian sin(a i · k). The eigenvalues of the tight-binding Hamiltonian vanish at nodal points at the corners of the hexagonal Brillouin zone. Each of the corners contributes with the fraction 1/3 to the total number of cones, which therefore is 2. At chemical neutrality, i.e., for Fermi energy laying precisely at nodal points, there is no extended Fermi surface, and it became common to talk about Fermi points or semimetals. Close to the Fermi points, the fermion dispersion is linear and therefore describes massless Dirac particles, cf. Figure 1. The two Dirac cones are not exactly equivalent though, but they differ by a subtle notion of chirality. The states corresponding to each of two cones can be thought of as the chiral partners of each other. The total chirality of the tight-binding Hamiltonian is therefore zero. Being interested in the physics at low energies, it is usually sufficient to use the effective low-energy Hamiltonian where ∇ ± = ∂ x ± i∂ y . To describe the 4 × 4 matrix body of the Hamiltonian, it is useful to introduce the Dirac matrices Σ ab = σ a ⊗ σ b , a, b = 0, 1, 2, 3, with σ a=1,2,3 denoting the Pauli matrices in their usual representation and σ a=0 being the two-dimensional unity matrix. The first index refers to the valley and the second refers to the sublattice degree of freedom.
With that, D ± = 1/2[Σ 01 ± iΣ 02 ] follows. The band gap ∆ 0 in pristine graphene is usually attributed to the spin-orbit coupling and has the size of roughly 10 −3 meV [11], but it can also be considered as a free parameter available for fine tuning. Finally, the chemical potential 0 is an adjustable quantity.

Topological Chern Number
The phase of the wave function plays a crucial role for the properties of the related physical system. It is associated with a topological invariant called the Chern number and is ultimately responsible for the robust macroscopic properties, such as for instance the famous universal conductance. The Chern number is defined as a contour integral [46] over the so-called Berry vector potential A(k) = −i v ± | ∇ k |v ± along any closed path in the reciprocal space. Here, |v ± denotes an eigenstate of the Hamiltonian defined in Equation (2). The Berry vector potential corresponding to the completely filled band of the full tight-binding Hamiltonian is shown in Figure 2. It appears to have the shape of a double vortex centered around the nodal points of the spectrum and demonstrates nicely the difference in chirality of the Dirac cones by whirling in opposite directions. The total Chern appears as the sum of Chern numbers from each eigenstate. Therefore, the total Chern number of the pristine graphene is zero, but this might change if a fundamental symmetry of the Hamiltonian is broken, e.g., by applying a magnetic field.

The Effective Hamiltonian in Strong Magnetic Field
In strong magnetic fields, we replace the usual derivatives by the covariant ones Here, we use the central gauge A = B/2(−y, x, 0) T , the choice which makes analytical calculations particularly convenient. Introducing complex coordinates z = x + iy,z = x − iy, and corresponding derivatives ∂ z = (∂ x − i∂ y )/2, ∂z = (∂ x + i∂ y )/2, with the properties where = 1/k is referred to as the magnetic length. The operator A annihilates the functions i.e., Aϕ n (r) = 0, for every positive integer n. The Gaussian part of Equation (6) guarantees the localization in the position space and makes it possible to carry out an integration in the position space exactly. The holomorphic part of Equation (6) contains only powers of z, and therefore, the wave function itself is manifestly chiral, which can be linked to the induced Berry curvature. The difference in the intrinsic Berry curvature discussed around Equation (4) is in the absence of the partner state with the opposite chirality, which reflects the explicit time-reversal symmetry breaking by an external magnetic field. The Hilbert space of the lowest Landau level is infinitely degenerate; i.e., n can assume every positive integer value between zero and infinity. In this notation, the Hamiltonian becomes The ground state (i.e., the eigenstate in the lowest Landau level) suffices the condition which suggests two solutions: which correspond to two valley polarization. The respective eigenvalues of the Hamiltonian for each spin projection in the lowest Landau level are found from the stationary Schrödinger equation which yields for both spectral branches (or Landau sublevels) [47] i.e., the spectrum of both Dirac electron species consists of two flat bands irrespective of the strength of the magnetic field. Moreover, in chemically neutral and gapless graphene, the spectrum in the lowest Landau level is at zero [48].

Single-Particle Propagator in the Lowest Landau Level
The advanced (+) or retarded (−) Green's function in the lowest Landau level can be calculated using the spectral representation where E s are the eigenvalues of the Hamiltonian for each spin projection in the lowest Landau level, as shown in Equation (11), and the normalization will be fixed later. The projectors P ± on the spin space are idempotent and orthogonal matrices with properties P + P − = 0, P + P + = P + , P − P − = P − . The summation over all n yields which then gives for the Green's function [32,49] Notably, the local Green's function (r = r ) is a coordinate independent constant. The propagator is normalized this way in order to satisfy the usual sum rule where the trace operator acts only on the spin space. Equation (16) gives the number of the elementary flux quanta φ 0 = h/e per unit volume. In the real-time representation, the Green's function represents a simple collection of undamped harmonic functions with the period determined by the eigenenergies of the lowest Landau level modes and the initial time is assumed to be at zero. The Green's function is totally separable on the space-time.
For the case of chemically neutral gapless graphene, the Green's function becomes particularly simple [49]: i.e., in the real-time representation, it is just a step function θ(t).

Static Conductivity of the Pristine Graphene vs. the Lowest Landau Level
The static conductivity of the clean system can be evaluated from the Kubo-Greenwood formula [28,29,34,35,50]: We first evaluate this expression for the pristine graphene without a magnetic field. The Green's function of such system reads where denotes the current operator, while the second power of the position operator can be written as Therefore, the Kubo-Greenwood formula changes to Taking into account we then get to with 4 being the trace of the unity matrix. Assuming an infinitely large upper cutoff, we finally get for the conductivity a universal number which is the famous universal conductivity of graphene [6]. Remarkably, this result is also valid for the case of weakly disordered Dirac electron gas [30,[51][52][53].
For the calculation of the static conductivity in the lowest Landau level of clean gapless and undoped graphene, we employ the Green's function shown in Equation (18). Here, we can evaluate the Kubo-Greenwood formula directly in the position space where 2 is the trace of the matrix P + + P − . In addition, here is the conductivity of a universal number, but its magnitude is only a quarter of the clean graphene. It is obvious that this result is solely due to the presence of the zero mode in the spectrum of the gapless and chemically neutral graphene. A slightest doping or a smallest spectral gap would destroy this dc conductivity. Because of this fragility, we can think of the resulting Equation (27) as an anomaly in the parametric space of infinitely small thickness. In analogy to the situation without a magnetic field, we expect the widening of this line by disorder [53].

Single-Particle Propagator Renormalization Due to the Disorder
The disorder is introduced in the form of the fluctuating chemical potential v(r), which couples in the spin space to the unity matrix Σ 00 , with the white noise correlator: The averaged propagator reads To perform the disorder average perturbative, we expand the propagator in powers of v. Because of the properties of the disorder correlator Equation (28), all terms with an odd number of potentials v vanish. The series then becomes Here, the summation over repeating indices is understood. The Green's function shown in Equation (15) is spanned by the spin projectors P s . Therefore, only the disorder diagonal in the spin space is of importance. In addition to the randomly fluctuating chemical potential considered here, these might include the randomly fluctuating gap, which couples to Σ 03 , the random "chiral" chemical potential (Σ 30 ), or the random "chiral" mass (Σ 33 ). Each product of these matrices with P s projects them bar the sign back onto P s again. Therefore, the perturbative series shown in Equation (30) does not depend on a particular disorder type, and our analysis is generic and disorder type independent.
The exact Green's function of disordered electrons in the lowest Landau level was obtained by Wegner in Ref. [32] in the distinctly separable form The frequency-dependent part of the Green's function evaluated to the order g 3 by evaluation of the diagrams shown in Figure 3 reads where E 2 g = gk 2 4π . The expansion coefficients 1, 1, 5/2, 37/4... are precisely those of the Wegner's exact solution [32]. They are determined as expansion coefficients of the function in powers of b/a 2 . Following [32], we find for the frequency-dependent part of the dressed single-particle propagator with the following explicit expressions for the real and imaginary parts [33,34] They depend on the dimensionless energy . Perturbative processes contributing to the dressing of the single-particle propagator due to the disorder to order g 1 (one diagram), g 2 (three diagrams), and g 3 (fifteen diagrams). Some of the diagrams of order g 3 should be counted twice because of the degeneracy due to the mirror symmetry with respect to the imaginable vertical axis, which is accounted for by the factors 2 in front of them.
In the chosen units, the disorder-related energy E g is a dimensionless quantity, which is proportional to the ratio of two relevant lengths E g ∼ l λ / : the magnetic length ∼ 1/k and the disorder related length l λ ∼ √ g. The total density of states is correctly normalized in accordance with Equation (16). Figure 4 shows the DOS from Equation (38). For weak disorder strength, the density of states that appears has the form of two sharp peaks placed symmetrically around the energy eigenvalues in the lowest Landau level. It is plotted in units of 1 l λ being the parametric volume constructed from the two specific lengths of the model. The peaks become broader with increasing disorder strength and overlap with each other until they merge to a single structure.
For the gapless and chemically neutral graphene, both peaks overlap and form a unique structure around the zero energy Therefore, at the band center, we get

Mean Squared Displacement of the Disordered System
The access to the diffusion goes via the mean squared displacement where r µ is the position operator and P rr (t) is the return probability density defined as where is the disorder averaged two-particles propagator. Notably, the numerator of Equation (41) appears to be essentially the Kubo-Greenwood formula shown in Equation (19). The relation between the mean squared displacement and diffusion is established via where D is the diffusion coefficient. If we would be able to determine the diffusion coefficient of the disordered system through the direct evaluation of Equation (41), then it will be possible to compute the conductivity from the Einstein relation (in this particular form adopted from [29,50]) where ρ(E) is the density of states discussed in the previous paragraph. A rigorous evaluation of the full perturbative series for the two-particles propagator G + r,0 G − 0,r g along the lines of Wegner's calculations for the single-particle propagator is principally impossible. Therefore, we need to consider the spatial averages. We evaluate both expressions from the numerator and denominator of Equation (41) perturbatively. Evaluation of all perturbative diagrams to order g 3 , shown in Figure 5 yields for the expression in the numerator of Equation (41) tr ∑ r r 2 µ P r0 (E) = 1 4 (2X s ) 4 + 343 72 (2X s ) 6 + · · · cos 2φ s + 139 72 (2X s ) 6 + · · · cos 4φ s + · · · , where According to Equations (35) and (36), X 2 s (E) and φ s (E) are dimensionless functions of the argument ν s = (E − E s )/E g . The analogous computation for the denominator of Equation (41) yields 2 (2X s ) 6 + · · · cos 2φ s + 5 2 (2X s ) 6 + · · · cos 4φ s + · · · .
A reasonable approximation for the two-particles propagator that leads beyond this partially rigorous result includes all diagrams of the so-called ladder channel. The four lowest order ladder diagrams are evaluated as [54] which suggests the following asymptotics of the two-particles propagator in the form of an infinite series: Using this expression, one can complement Equations (46) and (48) to any order.

Equation of Motion for the Mean Squared Displacement
The equation of motion for the mean squared displacement has the form of the secondorder ordinary differential equation where N refers to the order of the perturbative expansion. The expression for the integral I N can be found in [54]. For our purposes, it is only important that it decreases with increasing N. The large-time asymptotics (t t c , t c being some crossover time far in the past, which can be chosen zero) is given by the solution of the self-consistent equation where we introduced the scattering time as In the limit t → ∞, the mean squared displacement approaches its upper bound which for N → ∞ lies in the infinity and is therefore never reached. Hence, it can be only approached from below, which requires C to be negative. If I N is small, then τ is large, and the regime with linear time dependence should be broad. The diffusion coefficient is Formally, |C| should follow from the initial condition at t = t c , but for this, we need to know r 2 µ (t c ) > , which lies far in the past and is therefore forgotten. In order to be a physical quantity, we demand for D an invariance with respect to N. This is similar to the version of the renormalization group typically used in the high-energy physics. This implies from where then follows Even though this constrain might appear not entirely transparent, it has a natural analogy in the case of disordered electron gas without a magnetic field. Here, the diffusion coefficient is determined from the self-consistent Born approximation and appears unchanged in the partial series, e.g., cooperon or diffuson [22,23]. A comparison with Equation (44) suggests this constant to be 2. Then, the physical diffusion coefficient becomes i.e., it is proportional to the parametric volume of the model. Inserting the density of states from Equation (38) and the diffusion coefficient Equation (61) into the Einstein relation Equation (45) yields the conductivity. The system is conducting within a parametric window located around each of the Landau sublevels. The width of the conducting window is determined by the parameters of the microscopic model and by the disorder. The transition g → 0 is smooth, and the conductivity degenerates to two sharp peaks at the Landau sublevels. With increasing disorder, the peaks become broader and merge at some point to an amorphous structure. Simultaneously, the amplitude becomes smaller, signaling the suppression of the conductivity in the strong disorder limit.
We are now in the position to compute the conductivity at an arbitrary Landau sublevel defined in Equation (11). The corresponding density of the states is shown in Figure 4. For weak disorder, the contribution from the other sublevel is negligible, and we get for the density of states ρ sl (0) ≈ 1 Using the units of e 2 /h instead of e 2 /h adds an extra factor of 2π, i.e., Numerically, it is close to the dc conductivity in clean graphene off the magnetic field e 2 πh ≈ 0.32 e 2 h evaluated in Equation (26), which is valid even for the weakly disordered systems [30,[51][52][53]. For us, the most interesting limit is the case of the neutral and gapless graphene, for which a rich empirical knowledge is available. Our estimation would give for this zero mode Respected experimental studies of Ref. [55,56] determine the room-temperature longitudinal resistivity at the band center in the lowest Landau level as roughly 35 kΩ and 42 kΩ, respectively, which corresponds to (with h/e 2 ≈ 25, 812 Ω) to which is surprisingly close to our estimation. The comparison is justified, since the disorder can be considered as an effective temperature, cf. [57] and references therein.

Discussions
The diffusion of electrons in random environments confined to the lowest Landau level in two spatial dimensions is a long-standing and conceptually challenging problem of quantum statistical mechanics. Without disorder, the quantum mechanical description of the problem is simply that of the harmonic oscillator with a discrete, though highly degenerate spectrum, comprising of the so-called Landau levels. In the strong magnetic field, the gap between the lowest and the first Landau levels is very large and only the lowest Landau level is relevant. In this regime, the electrons are distributed between the stationary Landau orbits in the position space and should stay there forever, thus forbidding any transport across the sample. This is due to the disorder that the electrons can move from one Landau orbit to another, producing an observable current. While several analytical approaches have been developed in the past for systems of disorder electrons off a magnetic field, a meaningful formulation of the problem in the strong magnetic field is conceptually difficult, because the problem lacks a small parameter, and therefore, the perturbative expansions in powers of the disorder potential diverge. At the single-particle level, the problem has been solved by Wegner [32], who found an exact expression for the single-particle propagator. The Wegner propagator does not reveal any singularities and describes a state of the matter without pronounced resonances and consequently without clearly defined quasiparticles. However, the difficulties aggregate by far if one goes beyond the single-particle picture and considers processes involving two or more particles.
Our main intention is to calculate the static conductivity of disordered graphene in a strong magnetic field. In the lowest Landau level, the spectrum of the gapless and chemically neutral graphene Hamiltonian has zero energy, which are the consequence of its band topology and responsible for the exceptional transport properties under normal conditions. The spectral gap or fluctuations of the Fermi energy due to the hopping between second-nearest neighbors on the honeycomb lattice split this zero mode in two sublevels. Surprisingly, the static conductivity of the clean system evaluated from the Kubo-Greenwood formula gives a conductivity within infinitely thin parametric windows around the zero mode. One would expect that the disorder broadens this window to considerable sizes. However, because the perturbative series for the two-particles propagator diverges, a naive use of the Kubo formula fails. Therefore, we approach the conductivity via the Einstein relation, which requires the knowledge of the density of states and of the diffusion coefficient. While the former is known from the Wegner's solution, the latter is not. To deal with such divergences, we develop an analytical approach based on a self-consistent equation for the mean square displacement, which allows one to directly extract the diffusion coefficient and static conductivity.
Following the line of Wegner's exact considerations, we determine the general expression of the density of states of graphene. With a gap, the density of states of weakly disordered graphene represents two sharp peaks centered around each of the sublevels. For the case of gapless and chemically neutral graphene, both peaks coalesce to a single one with twice the height. The diffusion coefficient is extracted from the time evolution of the mean squared displacement. The latter tends toward a stationary state, which would reestablish the situation we observe in the clean system with all electrons distributed between stationary orbits. However, our findings suggest an infinitely large time needed for the system to arrive in this state. At the intermediary time scales, the mean squared displacement behaves lineary in time from which the diffusion coefficient is extracted. The combination of the density of states and diffusion coefficient, known as the Einstein relation, gives a universal, i.e., disorder independent value for the static conductivity. At the band center of the lowest Landau level, we find for the conductivity a universal value ∼ 0.72 e 2 /h, which is surprisingly close to the established results for the conductivity of the disordered Dirac electrons. In the subsequent work, we intend to extend our analysis to higher Landau levels and to address the Hall conductivity with the aim of arriving at an effective description of a kind of the Chern-Simons theories .
The quantum Hall effect has long become the standard tool for high-precision measurements and adjustments [58]. With our clearly laid out prediction for the dc conductivity in the lowest Landau level of graphene, we have provided a benchmark for prospective graphene-based metrological standardization devices.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.