Prerpints.org logo
Preprint
Article

Sub-Sharvin Conductance and Incoherent Shot-Noise in Graphene Disks at Magnetic Field

Altmetrics

Downloads

33

Views

14

Comments

0

A peer-reviewed article of this preprint also exists.

This version is not peer-reviewed

Submitted:

22 April 2024

Posted:

23 April 2024

You are already at the latest version

Alerts
Abstract
Highly-doped graphene samples show the conductance reduced and the shot-noise power enhanced compared to standard ballistic systems in two-dimensional electron gas. These features can be understood within a model assuming incoherent scattering of Dirac electrons between two interfaces separating the sample and the leads. Here we find, by adopting the above-mentioned model for the edge-free (Corbino) geometry and by means of the computer simulation of quantum transport, that another graphene-specific feature should be observable when the current flow through a doped disk is blocked by high magnetic field. In case the conductance drops to zero, the Fano factor approaches the value of $F\approx{}0.56$, with a very weak dependence on the disk radii ratio. The role of finite source-drain voltages and the system behavior upon tuning the electrostatic potential barrier from a rectangular to parabolic shape are also discussed.
Keywords: 
Subject: Physical Sciences  -   Condensed Matter Physics

1. Introduction

Although electronic properties of matter are governed by the rules of quantum mechanics [1], it is very unlikely to find that any measurable characteristic of a macroscopic system is determined solely by the universal constants of nature, such as the elementary charge (e) or the Planck constant (h). In the last century, two notable exceptions arrived with the phenomena of superconductivity [2], namely, the quantization of magnetic flux piercing the superconducting circuit, being the multiplicity of the flux quantum Φ 0 = h / ( 2 e ) [3,4], and the ac Josephson effect, with the universal frequency-to-voltage ratio given by 2 e / h = 1 / Φ 0 [5]. Later, with the advent of semiconducting heterostructures [6], came the quantum Hall effect [7,8,9,10,11,12] and the conductance quantization [13], bringing us with the conductance quantum g 0 = s e 2 / h (with the degeneracy s = 1 , 2, or 4). Further development of nanosystems led to the observation of Aharonov-Bohm effect manifesting itself by magnetoconductance oscillations with the period 2 Φ 0 = h / e [14], as well as the universal conductance fluctuations [15,16,17,18], characterized by a variance ∝ β − 1 ( s e 2 / h ) 2 , with an additional symmetry-dependent prefactor ( β = 1 , 2, or 4). Related, but slightly different issue concerns the Wiedemann-Franz (WF) law defining the Lorentz number, L 0 = π 2 3 ( k B / e ) 2 (with the Boltzmann constant k B ) [2], as the proportionality coefficient between electronic part of the thermal conductivity and electrical conductivity multiplied by absolute temperature. Although the WF law is followed, with a few-percent accuracy, in various condensed-matter systems, it has never been shown to have metrological accuracy [19,20,21,22,23,24].
Some new ’magic numbers’ similar to the mentioned above have arrived with the discovery of graphene, an atomically-thin form of carbon [11,12]. For undoped graphene samples, charge transport is dominated by transport via evanescent modes [25], resulting in the universal dc conductivity 4 e 2 / ( π h ) accompanied by the sub-Poissonian shot noise, with a Fano factor F = 1 / 3 [26,27,28,29,30,31]. For high frequencies, ac conductivity is given by π e 2 / ( 2 h ) , leading to the quantized visible light opacity π α (with α ≈ 1 / 137 . 036 being the fine-structure constant) [32,33,34]. A possible new universal value is predicted for the maximum absolute thermopower, which approaches ≈ k B / e near the charge neutrality point, for both monolayer and gapless bilayer graphene [39,40,41,42,43].
Away from the charge-neutrality point, ballistic graphene samples show the sub-Sharvin charge transport [35,36], characterized by the conductance reduced by a factor of π / 4 compared to standard Sharvin contacts in two-dimensional electron gas (2DEG) [37,38]. What is more, the shot noise is enhanced (comparing to 2DEG) up to F ≈ 1 / 8 far from the charge-neutrality [30,31]. Detailed dependence of the above-mentioned factors on a sample geometry was recently discussed in analytical terms [36], on the example of edge-free (Corbino) setup, characterized by the inner radius R i and the outer radius R o (see Figure 1). It is further found in Refs. [35,36] that the ballistic values of the conductance and Fano factor are gradually restored when the potential barrier, defining a sample area in the effective Dirac-Weyl Hamiltonian, evolves from a rectangular toward a parabolic shape.
Here, we focus on the Corbino geometry, which is often considered when discussing fundamental aspects of graphene [25,44,45,46,47,48,49,50,51,52,53,54]. In this geometry, charge transport at high magnetic fields is unaffected by edge states, allowing one to probe the bulk transport properties [48,49,50,51,52,53]. Recently, we have shown numerically that thermoelectric properties in such a situation are determined by the energy interval separating consecutive Landau levels rather then by the transport gap (being the energy interval, for which the cyclotron diameter 2 r c < R o − R i ) [54]. In this paper, we address a question how the shot-noise behaves when the tunneling conductance regime is entered by increasing magnetic field at a fixed doping (or decreasing the doping at a fixed field)? Going beyond the linear-response regime, we find that the threshold voltage U on , defined as a source-drain voltage difference that activates the current at minimal doping, is accompanied by quasi-universal (i.e., weakly-dependent on the radii ratio R o / R i ) value of F ≈ 0 . 56 . The robustness of the effect is also analyzed when smoothing the electrostatic potential barrier.
The paper is organized as follows. In Sec. Section 2 we briefly present the effective Dirac Hamiltonian and the numerical approach applied in remaining parts of the paper. In Sec. Section 3, we derive an approximation for the transmission through a doped Corbino disk at non-zero magnetic field and subsequent formulas for charge-transfer characteristics: the conductance and the Fano factor. Our numerical results, for both the rectangular and smooth potential barriers, are presented in Sec. Section 4. The conclusions are given in Sec. Section 5.

2. Model and Methods

2.1. Dirac Equation for the Disk Geometry

Our analysis of the device shown schematically in Figure 1 starts from effective wave equation for Dirac fermions in graphene, near the K valley,
v F ( p + e A ) · σ + V ( r ) Ψ = E Ψ ,
where the Fermi velocity is given by v F = 3 t 0 a / ( 2 ℏ ) , with t 0 = 2 . 7 eV the nearest-neighbor hopping integral and a = 0 . 246 nm the lattice parameter, p = − i ℏ ( ∂ x , ∂ y ) is the in-plane momentum operator, we choose the symmetric gauge A = B 2 ( − y , x ) corresponding to the perpendicular, uniform magnetic field B = ( 0 , 0 , B ) , and σ = ( σ x , σ y ) , where σ j are the Pauli matrices [55]. The electrostatic potential energy in Eq. (1), V ( r ) , is given by
V ( r ) = − V 0 × 2 m | r − R av | m | R o − R i | m if | r − R av | ⩽ R o − R i 2 , 1 if | r − R av | > R o − R i 2 ,
where we have defined R av = ( R i + R o ) / 2 . In particular, the limit of m → ∞ corresponds to the rectangular barrier (with a cylindrical symmetry); any finite m ⩾ 2 defines a smooth potential barrier, interpolating between the parabolic ( m = 2 ) and rectangular shape. In principle, barrier smoothing can be regarded as a feature of a self consistent solution originating from the diffusion of carriers; we expect this feature to strongly depend on the experimental details, with graphene-on-hBN devices [48] showing rectangular, rather then smooth, profiles.
Symmetry of the problem allows one to look for the wave function in the form
Ψ j ( r , φ ) = e i ( j − 1 / 2 ) φ χ a χ b e i φ ,
where j = ± 1 / 2 , ± 3 / 2 , ⋯ is the total angular-momentum quantum number, the components χ a = χ a ( r ) , χ b = χ b ( r ) , and we have introduced the polar coordinates ( r , φ ) . Substituting the above into Eq. (1) bring us to the system of ordinary differential equations
χ a ′ = j − 1 / 2 r + e B r 2 ℏ χ a + i E − V ( r ) ℏ v F χ b ,
χ b ′ = i E − V ( r ) ℏ v F χ a − j + 1 / 2 r + e B r 2 ℏ χ b ,
where primes denote derivatives with respect to r.

2.2. Analytic Solutions

For the disk area, R i < r < R o , Eqs. (4), (5) typically need to be integrated numerically; key details of the procedure are presented in Appendix A. Here we focus on the special case of rectangular barrier ( m = ∞ ), for which some analytic solutions were reported [44,45,56].
In particular, in the absence of magnetic field ( B = 0 ), the spinors χ j = ( χ a , χ b ) T corresponding to different j-s can be written as linear combinations [44]
χ j ( disk ) = A j H j − 1 / 2 ( 2 ) ( k r ) i η H j + 1 / 2 ( 2 ) ( k r ) + B j H j − 1 / 2 ( 1 ) ( k r ) i η H j + 1 / 2 ( 1 ) ( k r ) ,
where H ν ( 1 ) ( ρ ) [ H ν ( 2 ) ( ρ ) ] is the Hankel function of the first [second] kind, k = | E | / ( ℏ v F ) , the doping sign η = s g n E = ± 1 (with η = + 1 indicating electron doping and η = − 1 indicating hole doping), and A j , B j , are arbitrary complex coefficients. For B > 0 , Eq. (6) is replaced by [45,56]
χ j ( disk ) = A j ξ j ↑ ( 1 ) i η z j , 1 ξ j ↓ ( 1 ) + B j ξ j ↑ ( 2 ) i η z j , 2 ξ j ↓ ( 2 ) ,
where z j , 1 = [ 2 ( j + s j ) ] − 2 s j , z j , 2 = 2 ( β / k 2 ) s j + 1 / 2 (with s j ≡ 1 2 s g n j , β = e B / ( 2 ℏ ) ), and
ξ j s ( ν ) = e − β r 2 / 2 ( k r ) | l s | M ( α j s , γ j s , β r 2 ) , ν = 1 , U ( α j s , γ j s , β r 2 ) , ν = 2 ,
with l s = j ∓ 1 2 for s = ↑ , ↓ , α j s = 1 4 [ 2 ( l − s + | l s | + 1 ) − k 2 / β ] , and γ j s = | l s | + 1 . M ( a , b , z ) and U ( a , b , z ) are the confluent hypergeometric functions [57].
For the leads, r < R i or r > R o , the electrostatic potential energy is constant, V ( r ) = − V 0 . We further assume B = 0 and E > − V 0 (electron doping) in the leads, allowing one to adapt the wave function given by Eq. (6); i.e., for the inner lead, r < R i ,
χ j ( inner ) = H j − 1 / 2 ( 2 ) ( K r ) i H j + 1 / 2 ( 2 ) ( K r ) + r j H j − 1 / 2 ( 1 ) ( K r ) i H j + 1 / 2 ( 1 ) ( K r ) ,
and for the outer lead, r > R o ,
χ j ( outer ) = t j H j − 1 / 2 ( 2 ) ( K r ) i H j + 1 / 2 ( 2 ) ( K r ) ,
where K = | E + V 0 | / ( ℏ v F ) and we have introduced the reflection and transmission coefficient. The first spinor in each of Eqs. (9) and (10) represents the incoming (i.e., propagating from r = 0 ) wave, the second spinor in Eq. (9) represents the outgoing (propagating from r = ∞ ) wave.

2.3. Mode-Matching Method

Since the current-density operator following from Eq. (1), j = e v F σ does not involve differentiation, the mode-matching conditions for r = R i and r = R o reduce to the equalities for spinor components, namely
χ j ( inner ) ( R o ) = χ j ( disk ) ( R o ) and χ j ( disk ) ( R i ) = χ j ( outer ) ( R i ) .
The resulting formula for transmission probability for j-th mode becomes particularly simple upon taking the limit of heavily doped leads, U 0 → ∞ . In particular, for B = 0 , substituting Eq. (6) into the above gives [58]
T j = | t j | 2 = 16 π 2 k 2 R i R o 1 D j ( + ) 2 + D j ( − ) 2 ,
where
D j ( ± ) = I m H j − 1 / 2 ( 1 ) ( k R i ) H j ∓ 1 / 2 ( 2 ) ( k R o ) ± H j + 1 / 2 ( 1 ) ( k R i ) H j ± 1 / 2 ( 2 ) ( k R o ) .
Analogously, for B > 0 one finds, using Eqs. (7) and (8),
T j = | t j | 2 = 16 ( k 2 / β ) | 2 j − 1 | k 2 R i R o ( X j 2 + Y j 2 ) Γ ( γ j ↑ ) Γ ( α j ↑ ) 2 ,
where Γ ( z ) is the Euler Gamma function, and
X j = w j ↑ ↑ − + z j , 1 z j , 2 w j ↓ ↓ − , Y j = z j , 2 w j ↑ ↓ + − z j , 1 w j ↓ ↑ + , w j s s ′ ± = ξ j s ( 1 ) ( R i ) ξ j s ′ ( 2 ) ( R o ) ± ξ j s ( 1 ) ( R o ) ξ j s ′ ( 2 ) ( R i ) .
For B < 0 , one gets T j ( B ) = T − j ( − B ) .
Details of numerical mode-matching, applicable for smooth potentials, are given in Appendix A.

2.4. Landauer-Büttiker Formalism

In case the nanoscopic system is connected to external reservoirs, characterized by the electrochemical potentials μ and μ + e U eff (for simplicity, the two reservoirs are considered; for more general discussion see Ref. [59]), the conductance of the system is related to the transmission probabilities for normal modes ( T j -s) via
G ( U eff ) = 〈 I 〉 U eff = g 0 U eff ∫ μ μ + e U eff d ϵ ∑ j T j ( ϵ ) ,
where 〈 I 〉 denotes the average electric current and the zero-temperature limit is taken. The conductance quantum is g 0 = 4 e 2 / h , taking into account spin and valley degeneracies. U eff denotes the effective voltage difference between the reservoirs (notice that the actual voltage applied may differ from U eff due to charge-screening effects). Similarly, the Fano factor, relating the current variance, I − 〈 I 〉 2 , to the value I − 〈 I 〉 2 Poisson one would measure in the absence of correlations between scattering events (occurring, e.g., in the tunneling limit of T j ≪ 1 for all j-s), is given by
F ( U eff ) = I − 〈 I 〉 2 I − 〈 I 〉 2 Poisson = g 0 G U eff ∫ μ μ + e U eff d ϵ ∑ j T j ( ϵ ) 1 − T j ( ϵ ) ,
where I − 〈 I 〉 2 Poisson = e 〈 I 〉 / Δ t = e G U eff / Δ t with Δ t denoting the time of measurement.
In the linear-response regime ( U eff → 0 ), Eqs. (16) and (17) reduces to
G ( U eff → 0 ) = g 0 ∑ j T j ,
and
F ( U eff → 0 ) = ∑ j T j ( 1 − T j ) ∑ j T j ,
where T j = T j ( μ ) . For the disk geometry, summation range is limited by the number of propagating modes in the inner lead, | j | ⩽ j max = ⌊ K R i ⌋ − 1 2 with ⌊ x ⌋ denoting the floor function of x. (For heavily-doped leads, j max → ∞ .)
As a notable example, we consider the zero-doping limit ( μ → 0 ). In such a case, Eq. (14) simplifies to [45,60]
T j ( μ → 0 ) = 1 cosh 2 [ ( j + Φ / Φ 0 ) ln ( R o / R i ) ] ,
where Φ = π ( R o 2 − R i 2 ) B is the flux piercing the disk area and we have defined Φ 0 = 2 ( h / e ) ln ( R o / R i ) . Assuming the narrow-disk range, R o ≈ R i , we can approximate the sums occurring in Eqs. (18) and (19) by integrals, obtaining
G ≈ G diff = 2 π σ 0 ln ( R o / R i ) and F ≈ F diff = 1 3 .
The above reproduces pseudodiffusive conductance and the shot-noise power for a disk geometry [44]. For larger R o / R i , both characteristics are predicted to show approximately sinusoidal conductance oscillations with the field B [45,60,61].
The case of doped disk, for which one may expect to observe some features of the sub-Sharvin charge transport [35,36], is discussed next.

3. Approximate Conductance and Fano Factor at the Magnetic Field

Before calculating the conductance G and Fano factor F within the mode-matching method described in Sec. Section 2, we first present the approximating formulas for incoherent transport, obtained by adapting the derivation of Ref. [36] for the B > 0 case.

3.1. Corbino Disk in Graphene as a Double Barrier

A key step in the derivation is to observe that, in the multimode regime ( k R i ≫ 1 ) for which one can consider well-defined trajectories, the disk symmetry cause that incident angles θ 1 and θ 2 , corresponding to the interfaces at r = R i and r = R o (see Figure 2) remain constant (up to a sign) after multiple scatterings. Therefore, one can apply the double-contact formula for incoherent transmission [62,63], namely
T incoh = 1 2 π ∫ − π π d ϕ T 1 T 2 2 − T 1 − T 2 + T 1 T 2 − 2 ( 1 − T 1 ) ( 1 − T 2 ) cos ϕ = T 1 T 2 T 1 + T 2 − T 1 T 2 ,
where the transmission probabilities T 1 , T 2 , corresponding to a potential step of infinite height, are given by
T l = 2 cos θ l 1 + cos θ l , l = 1 , 2 ,
and ϕ is assumed to be a random phase acquired during the propagation between r = R i and r = R o (or vice versa). Similarly, we calculate the incoherent squared transmission, useful when evaluating the Fano factor,
T 2 incoh = 1 2 π ∫ − π π d ϕ T 1 T 2 2 − T 1 − T 2 + T 1 T 2 − 2 ( 1 − T 1 ) ( 1 − T 2 ) cos ϕ 2 = ( T 1 T 2 ) 2 ( 2 − T 1 − T 2 + T 1 T 2 ) ( 1 + T 1 T 2 − T 1 T 2 ) 3 .
Next, the incoherent conductance in the linear-response regime is evaluated by inserting T incoh (22) into Eq. (18),
G incoh = G Sharvin T incoh u = sin θ 1 ,
with
G Sharvin = 2 g 0 k R i .
For the Fano factor, one can analogously derive from Eq. (19)
F incoh = 1 − T 2 incoh u = sin θ 1 T incoh u = sin θ 1 .
The summation over 2 k R i modes is approximated in Eqs. (25), (27) by averaging over the variable u = sin θ 1 , within the range of − 1 ⩽ u ⩽ 1 . Explicitly,
T n incoh u = sin θ 1 = 1 2 ∫ u c 1 d u T n incoh , n = 1 , 2 ,
where the lower integration limit ( u c ) is defined via the value of sin θ 1 , below which the trajectory cannot reach the outer interface ( r = R o ). (In other words, for u = sin θ 1 < u c , the geometric derivation to be presented below leads to | sin θ 2 | > 1 .)
The missing elements, necessary to calculate T n incoh u = sin θ 1 in Eq. (28) is the dependence of θ 2 on θ 1 and B [see Eqs. (22), (23), and (24)], as well as the dependence of u c on B. Since we have assumed constant electrostatic potential energy in the disk area, the trajectory between subsequent scatterings (see Figure 2) forms an arc, with the constant radii
r c = ℏ k / ( e B ) = | E | / ( v F e B ) ,
(i.e., the cyclotron radius for massless Dirac particle at B > 0 ), centered at the distance r x from the origin. Now, solving the two triangles with a common edge r x (dashed line) and the opposite vertices in two scattering points, we find
r x 2 = R i 2 + r c 2 + 2 R i r c sin θ 1
(for the triangle containing a scattering point at r = R i ), and
r x 2 = R o 2 + r c 2 − 2 R o r c sin θ 2
(for the triangle containing a scattering point at r = R o ). Together, Eqs. (30) and (31) lead to
sin θ 2 = R o 2 − R i 2 − 2 R i r c u 2 R o r c .
Subsequently, the value of u c in Eq. (28) is given by
u c = − 1 , if B ⩽ B c , 1 R o 2 − R i 2 2 R i r c − R o R i , if B c , 1 < B ⩽ B c , 2 1 , if B > B c , 2 ,
where we have additionally defined
B c , m = 2 ℏ k e R o − ( − 1 ) m R i , m = 1 , 2 .

3.2. The Zero-Field Limit

Typically, averages occurring in Eqs. (25) and (27) need to be evaluated numerically. Analytic expressions are available, e.g., for zero magnetic field [36]
G incoh ( B → 0 ) = G Sharvin ( 2 a + 1 a ) arcsin a + 3 1 − a 2 − π 2 ( a 2 + 2 ) 1 − a 2 , F incoh ( B → 0 ) = 2 a 1 − a 2 ( 53 + 279 a 2 + 88 a 4 ) − 3 π a ( 12 + 82 a 2 + 45 a 4 + a 6 ) + 6 ( 1 + 45 a 2 + 82 a 4 + 12 a 6 ) arcsin a
/ 6 ( 1 − a 2 ) 2 π a ( a 2 + 2 ) − 6 a 1 − a 2 − 2 ( 2 a 2 + 1 ) arcsin a ,
where we have defined the inverse radii ratio a = R i / R o .

3.3. The Zero-Conductance Limit

In the this paper, we focus on the limit of B → B c , 2 − (i.e., B approaching B c , 2 from below), for which G incoh → 0 . Introducing the dimensionless 0 < ε ≪ 1 , one can express the cyclotron diameter, see Eq. (29), as
2 r c = R o − R i + ε R o − R i .
In turn, the value u c , see Eq. (33), can be approximated (up to the leading order in ε ) as
u c ≈ 1 − ε 1 + R o R i .
It is now convenient to define the variable
α = 1 − u ε 1 + R o / R i ,
such that the integration over u c ⩽ u ⩽ 1 , occurring when evaluating T n incoh u = sin θ 1 from Eq. (28), can be replaced by integration over 1 ⩾ α ⩾ 0 . Transmission probabilities ( T 1 , T 2 ) for the interfaces at r = R i and r = R o , see Eqs. (22), (23), (24) and (32), can now be approximated as
T 1 ≈ 2 2 α 1 + 1 a ε 1 / 2 ,
T 2 ≈ 2 2 ( 1 − α ) 1 + a ε 1 / 2 ,
where we have used a = R i / R o again.
Using the above expressions, we can now rewrite the averages occurring in Eq. (28), up to the leading order in ε again, as follows
T incoh u = G incoh G Sharvin ≈ 2 1 + 1 a 2 ε 3 / 2 ∫ 0 1 d α α ( 1 + 1 a ) ( 1 − α ) ( 1 + a ) α ( 1 + 1 a ) + ( 1 − α ) ( 1 + a ) ,
T 2 incoh u ≈ 2 1 + 1 a 2 ε 3 / 2 ∫ 0 1 d α α ( 1 − α ) 2 + 1 a + a α ( 1 + 1 a ) + ( 1 − α ) ( 1 + a ) 3 .
Remarkably, both quantities decays as ∝ ε 3 / 2 , but their ratio, occurring in Eq. (27) for the Fano factor, remains constant (for a given a). The integrals in Eqs. (42) and (43) can be calculated analytically, leading to
F incoh ( B → B c , 2 − ) = 1 − 4 a + 39 a − 58 a 3 / 2 − 23 a 2 − 23 a 5 / 2 − 58 a 3 + 39 a 7 / 2 + 4 a 4 + 69 a 2 − 18 a 3 − 18 a 1 + a a r t a n h a 1 + a + 69 a 2 − 18 a − 18 a 3 1 + a a r t a n h 1 1 + a / ( 1 + a ) 3 1 + a 1 − 3 a + a + 3 a ( 1 + a ) 5 / 2 a r t a n h ( 1 + a ) 1 + a 1 + a + a .
Numerical values of F incoh ( B → 0 ) and F incoh ( B → B c , 2 − ) for selected a = R i / R o are given in Table 1. For F incoh ( B → B c , 2 − ) , we see that the poissonian value of F incoh = 1 , which one could expect due to the vanishing conductance, is reconstructed only for a → 0 (i.e., for R o ≫ R i ). For finite radii ratios, nontrivial values of 0 < F incoh < 1 occurs. Remarkable, for moderate disk proportions ( a ⩾ 0 . 5 ), F incoh ( B → B c , 2 − ) shows very weak dependence on a, decaying by less then 2 % (from F incoh ≈ 0 . 56 at a = 0 . 5 to F incoh ≈ 0 . 55 for a → 1 , with a → 1 representing the narrow-disk limit of R o ≈ R i ).
For this reason, in the following numerical analysis, we fixed the disk radii ratio at a = 0 . 5 (i.e., R o = 2 R i ). We also stress that the derivation presented above, holds true for the parameter ε → 0 + , quantifying the ratio of cyclotron diameter 2 r c to radii difference R o − R i , see Eq. (37). Therefore, it is irrelevant whether one increases the magnetic field at a fixed chemical potential, or reduces the chemical potential at a fixed B > 0 (as long as the system stays in a multimode range, k R i ≫ 1 ).

4. Results and Discussion

Main doubt arising when we consider the applicability of Eq. (44) for real quantum systems concerns the possible role of evanescent waves, totally neglected in our derivation. Obviously, they should not play an important role when the system is highly conducting (such as in the zero-field case [36]); however, since the Fano factor is determined by the ratio of two cumulants, both vanishing for sufficiently high field, it is not fully clear which contribution (from propagating or from evanescent modes) would govern the value of F for B → B c , 2 ? On the other hand, resonances with Landau levels are not expect to play a significant role, as they form very narrow transmission peaks, contributions of which get immediately smeared out beyond the linear-response regime.
In the remaining parts of the paper, compare the results of computer simulation of quantum transport through the disk in graphene, with the predictions for incoherent scattering presented in Sec. Section 3, in attempt to propose an experimental procedure allowing one to extract the nontrivial value of F ≈ 0 . 55 from the data plagued with other contributions.

4.1. The Rectangular Barrier of an Infinite Height

As a first numerical example, we took the limit of V 0 → ∞ and m → ∞ in Eq. (2), for which close-form expressions for transmission probabilities were presented in Sec. Section 2.
In Figure 3, we compare the linear-response conductance G ( U eff → 0 ) , see Eq. (18), with G ( U eff ) calculated from Eq. (16) for a small but nonzero value of U eff = 0 . 01 V, both displayed as functions of the chemical potential. Also in Figure 3, same comparison is presented for the Fano factor F ( U eff ) [see Eqs. (19) and (17)]. It easy to see that prominent, aperiodic oscillations visible for both charge-transfer cumulants in the U eff → 0 limit are significantly reduced even for small U eff > 0 . In fact, for U eff = 0 . 01 V and B > 0 , the values of F incoh calculated from Eq. (27) [black lines] are closely followed by F ( U eff ) obtained from the numerical mode-matching, as long as the former can be defined, i.e., for B < B c , 2 at a given μ . We further notice that the value of μ for which B ≈ B c , 2 and F ( U eff ) ≈ 0 . 56 corresponds to G ( U eff ) ∼ g 0 (up to the order of magnitude). For smaller μ , such that B > B c , 2 and F incoh is undefined, F ( U eff ) saturates near the value ≈ 0 . 75 , apparently below the poissonian limit of F = 1 .
To better understand the nature of the results we now, in Figure 4, go further beyond the linear-response regime, calculating G ( U eff ) and F ( U eff ) for μ = − e U eff / 2 (notice that for an infinite rectangular barrier we have the particle-hole symmetry, and both cumulants are even upon μ ↔ − μ + e U eff ) and displaying them as functions of U eff .
We further introduce the activation voltage U on = U on ( B ) , meaning of which can be understood as follows. The cyclotron diameter, see Eq. (29), naturally defines the range of energies for which 2 r c ( E ) < R o − R i and the system shows G ≈ 0 (up to the evanescent modes). On the other hand, as we have set μ = − e U eff / 2 , the effective voltage defines the energy range of | E | ⩽ e U eff / 2 , being the integration interval in Eqs. (16) and (17). In turn, G ( U eff ) > 0 is expected for U eff ⩾ U on , a value of which can be approximated as
U on , incoh = v F B R o − R i ,
where we have simply rewrite equality 2 r c ( e U on ) = R o − R i neglecting the evanescent modes.
When looking on the conductance spectra illustrated in Figure 4(a) we see, for B > 0 , a wide range of lower U eff for which G ≈ 0 , attached (via a cusp region) to the range of (approximately linearly) increasing G. In order to determine the value of U on ( B ) directly from the conductance spectra G ( U eff ) , we find numerically the value of U on ( 1 ) such that G ( U on ( 1 ) ) = g 0 , and U on ( 2 ) such that G ( U on ( 2 ) ) = 2 g 0 , see the datapoints in Figure 4(a). Then, the linear extrapolation is performed to obtain
U on ( 0 ) = U on ( 1 ) − U on ( 2 ) − U on ( 1 ) G ( U on ( 1 ) ) G ( U on ( 2 ) ) − G ( U on ( 1 ) ) = 2 U on ( 1 ) − U on ( 2 ) ,
such that G ( U on ( 0 ) ) ≈ 0 . The resulting values of U on ( i ) , depicted in Figure 4(c) [datapoints], stay close to U on , incoh obtained from Eq. (45) [dashed line].
Remarkably, the values of the Fano factor corresponding to U eff = U on ( i ) , i = 1 , 2 , see Figure 4(b), are close to F incoh ( B → B c , 2 − ) ≈ 0 . 56 . Similar observation applies for all studied values of B ⩽ 0 . 5 T, see Figure 4(d); a typical deviation does not exceed 5 % .

4.2. Smooth Potential Barriers

In this subsection, we extend our numerical analysis onto smooth potential barriers, defined by choosing 2 ⩽ m < ∞ in Eq. (2). Moreover, the barrier height is now finite, i.e., V 0 = t 0 / 2 = 1 . 35 eV, being not far from the results of some first-principles calculations for graphene-metal structures [64,65]. According to our best knowledge such a model, first proposed in Ref. [35], seems to be the simplest providing qualitatively correct description of the conductance-spectrum asymmetry observed in existing experiments [46,50,53], in which the conductance for μ < 0 is noticeably suppressed, comparing to the μ > 0 range, due to the presence of two circular p-n junctions in the former case. (Such a feature is also correctly reproduced by a simpler model assuming the trapezoidal potential barrier [66], allowing a fully analytic treatment, but this approach produces an artificial conductance maximum near μ = 0 .)
The conductance spectra for five selected values of m are displayed in Figure 5, both for the linear-response regime [see Figure 5(a,c)] and beyond [Figure 5(b,d)]. This time, we have limited out presentation a single value of magnetic field, i.e., B = 0 . 2 T. It must be notice that finite value of V 0 results in small, but visible spectrum asymmetry also for m = ∞ .
The finite-voltage results, G ( U eff ) at μ = − e U eff / 2 , allows as to determine the activation voltage, U on ( B ) , in a similar manner as for an infinite-barrier case (see previous subsection). When attempting to apply the incoherent-scattering approximation to smooth potentials, some modification is required for Eq. (45), which now can be rewritten as
U on , incoh = v F B L diff ( m ) .
In the above, we have introduced the m-dependent effective sample length given by [35,36]
L diff ( m ) = | R o − R i | ℏ v F | R o − R i | V 0 1 / m ,
which reduces to L diff ( ∞ ) = R o − R i for a rectangular barrier, and gives L diff ( m = 2 ) ≪ R o − R i for the parabolic case. In brief, Eq. (48) can be derived by imposing V ( ± L diff / 2 ) = − E diff , where E diff denotes the value of Fermi energy, above which Sharvin conductance overrules the pseudodiffusive conductance, namely,
E diff = ℏ v F R o − R i ≈ 1 meV for R o − R i = 500 nm .
In Figure 6 we show the Fano factor, for same five values of m as previously used for the conductance (see Figure 5), and B = 0 . 2 T, as a function μ in the linear-response limit ( U eff → 0 ), as well as a function U eff for μ = − U eff / 2 (see left or right side of Figure 6, respectively). Again, the aperiodic oscillations almost vanish when entering the nonlinear response regime; in fact, the shape of F max ( U eff ) appears to be much less sensitive to the value of m than the linear-response F ( μ ) . Datapoints on the right side of Figure 6, identifying the values of F ( U on ( i ) ) , i = 1 , 2 , such that G ( U on ( i ) ) = i g 0 (see Figure 5), are available starting from m = 8 (although the deviation from F incoh ( B → B c , 2 − ) ≈ 0 . 56 is significant in such a case), whereas strong asymmetry of F ( μ ) is visible up to m = 32 .
The values of U on ( i ) and the corresponding F ( U on ( i ) ) , for the magnetic fields up to B ⩽ 0 . 5 T, are displayed in Figure 7. It can be noticed that the voltages U on ( i ) , see datapoints in Figure 7(a–d), show relatively good agreement with the approximation given by Eq. (47) [purple solid lines]; in fact, significant deviation from Eq. (45) relevant for the rectangular barrier [black dashed lines] can be noticed for m = 8 only. On the contrary, corresponding Fano factors F ( U on ( i ) ) , see datapoints in Figure 7(e–h), stay close to the value of F incoh ( B → B c , 2 − ) ≈ 0 . 56 only for m = ∞ and m = 128 , showing that the incoherent treatment of the shot-noise power, which we put forward in Sec. Section 3, is applicable only if the potential profiles is close to (but not necessarily perfectly matching) the rectangular shape.

5. Conclusions

We have put forward an analytic description of the shot-noise power in graphene-based disks in high magnetic field and doping. Assuming the incoherent scattering of Dirac fermions between two potential steps of an infinite height, both characterized by a priori nonzero transmission probability due to the Klein tunneling, we find that vanishing conductance should be accompanied by the Fano factor F ≈ 0 . 56 , weakly-dependent on the disk proportions.
Next, the results of analytic considerations are confronted with the outcome of computer simulations, including both rectangular and smooth shapes of the electrostatic potential barrier in the disk area. Calculating both linear-response and finite-voltage transport cumulants, within the zero-temperature Landauer-Büttiker formalism, we point out that the role of evanescent waves (earlier ignored in the analytic approach) is significant in the linear-response regime, however, one should able to detect the quasi-universal F ≈ 0 . 56 noise in a properly designed experiment going beyond the linear response regime. To achieve this goal, the following procedure is proposed: First, the activation voltage (for a fixed magnetic field) needs to be determined, by finding a cusp position on the conductance-versus-voltage plot, above which the conductance grows fast with the voltage (the average chemical potential is controlled by the gate such that the conductance is minimal for a given voltage). Having the activation voltage determined, one measures the noise for such a voltage, expecting the Fano factor to be close to F ≈ 0 . 56 .
We expect that the effect we describe should be observable in ultraclean samples and sub-kelvin temperatures (such as in Ref. [48]); for higher temperatures, hydrodynamic effects may noticeably alter the measurable quantities [53]. Since the noise-related characteristics seem to be generally more sensitive to the potential shape then the conductance (or the thermoelectric properties earlier discussed in Ref. [54]), the experimental study following the scenario presented here may be a suitable way to check whether the flat-potential area of a mesoscopic size is present or not in a given graphene-based structure.

Author Contributions

A.R. designed the algorithm, A.R. and P.W. developed the code and performed preliminary computations, K.R. organized the computations on the PL-Grid supercomputing infrastructure; all authors were involved in data analysis and manuscript preparation.

Acknowledgments

The main part of the work was supported by the National Science Centre of Poland (NCN) via Grant No. 2014/14/E/ST3/00256 (SONATA BIS). Computations were performed using the PL-Grid infrastructure.

Appendix A. Numerical Mode Matching for Smooth Potentials

Here we summarize the numerical approach earlier presented in Ref. [54].
In a typical situation, system of ordinary differential equations for spinor components ( χ a , χ b ) , see Eqs. (4) and (5), needs to be integrated numerically for all j-s. In order to reduce round-off errors that may occur in finite-precision arithmetics due to exponentially growing (or decaying) solutions, one can divide the full interval, R i < r < R o , into M parts, bounded by
R c ( l ) = R i + l R o − R i M < r < R c ( l + 1 ) , with l = 0 , 1 , ⋯ , M − 1 .
(In particular, R c ( 0 ) = R i and R c ( M ) = R o .)
The wave function in the disk area χ j ( disk ) is now given by a series of functions χ j ( l ) for M consecutive intervals given by Eq. (A50). For the l-th interval,
χ j ( l ) = A j ( l ) χ j ( l ) , I + B j ( l ) χ j ( l ) , II ,
where χ j ( l ) , I , χ j ( l ) , II are two linearly independent solutions obtained by integrating Eqs. (4), (5) with two different initial conditions, χ j ( l ) , I r = R i ( l ) = ( 1 , 0 ) T and χ j ( l ) , II r = R i ( l ) = ( 0 , 1 ) T . A j ( l ) and B j ( l ) are complex coefficients (to be determined later).
In particular, for R o = 2 R i = 1000 nm and B < 0 . 5 T considered in this paper, it is sufficient to set M = 20 and employ a standard fourth-order Runge-Kutta (RK4) algorithm with a spatial step of 0 . 5 pm. (For such a choice, the output numerical uncertainties of transmission probabilities T j are smaller then 10 − 7 .)
The matching conditions for the M + 1 interfaces at r = R i , r = R c ( 1 ) , ..., r = R c ( M − 1 ) , and r = R o , can now be written as
χ j ( inner ) ( R i ) = χ j ( 0 ) ( R i ) ,
χ j ( l ) ( R c ( l + 1 ) ) = χ j ( l + 1 ) ( R c ( l + 1 ) ) , l = 0 , ⋯ , M − 2 ,
χ j ( M − 1 ) ( R o ) = χ j ( outer ) ( R o ) ,
and are equivalent to the Cramer’s system of 2 ( M + 1 ) linear equations for the unknowns A j ( 0 ) , B j ( 0 ) , ..., A j ( M − 1 ) , B j ( M − 1 ) , r j , and t j .
Writing down the spinor components appearing in Eqs. (A52), (A53), and (A54) explicitly, we arrive to
− χ j , a out ( R i ) χ j , a ( 0 ) , I ( R c ( 0 ) ) χ j , a ( 0 ) , II ( R c ( 0 ) ) − χ j , b out ( R i ) χ j , b ( 0 ) , I ( R c ( 0 ) ) χ j , b ( 0 ) , II ( R c ( 0 ) ) 0 χ a ( 0 ) , I ( R c ( 1 ) ) χ j , a ( 0 ) , II ( R c ( 1 ) ) 0 χ b ( 0 ) , I ( R c ( 1 ) ) χ j , b ( 0 ) , II ( R c ( 1 ) ) ⋱ ⋱ χ j , a ( M ¯ ) , I ( R c ( M ¯ ) ) χ j , a ( M ¯ ) , II ( R c ( M ¯ ) ) 0 χ j , b ( M ¯ ) , I ( R c ( M ¯ ) ) χ j , b ( M ¯ ) , II ( R c ( M ¯ ) ) 0 χ a ( M ¯ ) , I ( R c ( M ) ) χ j , a ( M ¯ ) , II ( R c ( M ) ) − χ j , a in ( R o ) χ b ( M ¯ ) , I ( R c ( M ) ) χ j , b ( M ¯ ) , II ( R c ( M ) ) − χ j , b in ( R o )
× r j A j ( 0 ) B j ( 0 ) ⋮ A j ( M − 1 ) B j ( M − 1 ) t j = χ j , a in ( R i ) χ j , b in ( R i ) 0 ⋮ 0 ,
where we have defined M ¯ = M − 1 ,
χ j in = H j − 1 / 2 ( 2 ) ( K r ) i H j + 1 / 2 ( 2 ) ( K r ) , χ j out = H j − 1 / 2 ( 1 ) ( K r ) i H j + 1 / 2 ( 1 ) ( K r ) .
For heavily-doped leads ( V 0 → ∞ ) the wave functions given by Eq. (A56) simplify to
χ j ( in ) = e i K r r 1 1 , χ j ( out ) = e − i K r r 1 − 1 ,
with K = | E + V 0 | / ( ℏ v F ) → ∞ .
As the linear systems for different values of j-s are decoupled, standard software packages can be used to find their solutions. We have chosen the double precision LAPACK routine zgesv, see Ref. [67].

References

  1. Datta, S. Quantum Transport: Atom to Transistor, Cambridge University Press: Cambridge, UK, 2005.
  2. Kittel, C. Introduction to Solid State Physics, 8th ed. John Willey and Sons: New York, NY, USA, 2005.
  3. Deaver, Jr., B.S.; Fairbank, W.M. Experimental Evidence for Quantized Flux in Superconducting Cylinders. Phys. Rev. Lett. 1961, 7, 43. [CrossRef]
  4. Doll, R.; Näbauer, M. Experimental Proof of Magnetic Flux Quantization in a Superconducting Ring. Phys. Rev. Lett. 1961, 7, 51. [CrossRef]
  5. Josephson, B.D. The discovery of tunneling supercurrents. Rev. Mod. Phys. 1974, 46, 251. [CrossRef]
  6. Imry, Y. Introduction to Mesoscopic Physics, 2nd ed. Oxford University Press: Oxford, UK, 2002; Chapter 5.
  7. Ando, T.; Matsumoto, Y.; Uemura, Y. Theory of Hall Effect in a Two-Dimensional Electron System. J. Phys. Soc. Jpn. 1975, 39, 279. [CrossRef]
  8. v. Klitzing, K,; Dorda, G.; Pepper, M. New Method for High-Accuracy Determination of the Fine-Structure Constant Based on Quantized Hall Resistance. Phys. Rev. Lett. 1980, 45, 494. [CrossRef]
  9. Laughlin, R.B. Quantized Hall conductivity in two dimensions. Phys. Rev. B 1981, 23, 5632(R). [CrossRef]
  10. Tsui, D.C.; Stormer, H.L.; Gossard, A.C. Two-Dimensional Magnetotransport in the Extreme Quantum Limit. Phys. Rev. Lett. 1982, 48, 1559. [CrossRef]
  11. Novoselov, K.S.; Geim, A.K.; Morozov, S.V.; Jiang, D.; Katsnelson, M.I.; Grigorieva, I.V.; Dubonos, S.V.; Firsov, A.A. Two-dimensional gas of massless Dirac fermions in graphene. Nature 2005, 438, 197. [CrossRef]
  12. Zhang, Y.; Tan, Y.-W.; Stormer, H.L.; Kim, P. Experimental observation of the quantum Hall effect and Berry’s phase in graphene. Nature 2005, 438, 201.
  13. van Wees, B.J.; van Houten, H.; Beenakker, C.W.J.; Williamson, J.G.; Kouwenhoven, L.P.; van der Marel, D.; Foxon, C.T. Quantized conductance of point contacts in a two-dimensional electron gas. Phys. Rev. Lett. 1988, 60, 848. [CrossRef]
  14. Webb, R.A.; Washburn, S.; Umbach, C.P.; Laibowitz, R.B. Observation of h/e Aharonov-Bohm Oscillations in Normal-Metal Rings. Phys. Rev. Lett. 1985, 54, 2696. [CrossRef]
  15. Lee, P.A.; Stone, A.D. Universal Conductance Fluctuations in Metals. Phys. Rev. Lett. 1985, 55, 1622. [CrossRef]
  16. Altshuler, B.L.; Shklovskii, B.I. Repulsion of energy levels and conductivity of small metal samples. Zh. Eksp. Teor. Fiz. 91, 220 1986 [Sov. Phys. – JETP 1986, 64, 127].
  17. Pal, A.N.; Kochat, V.; Ghosh, A. Direct Observation of Valley Hybridization and Universal Symmetry of Graphene with Mesoscopic Conductance Fluctuations. Phys. Rev. Lett. 2012, 109, 196601. [CrossRef]
  18. Hu, Y.; Liu, H.; Jiang, H.; Xie, X.C. Numerical study of universal conductance fluctuations in three-dimensional topological semimetals. Phys. Rev. B 2017, 96, 134201. [CrossRef]
  19. Sharapov, S.G.; Gusynin, V.P.; Beck, H. Transport properties in the d-density-wave state in an external magnetic field: The Wiedemann-Franz law. Phys. Rev. B 2003, 67, 144509. [CrossRef]
  20. Mahajan, R.; Barkeshli, M.; Hartnoll, S.A. Non-Fermi liquids and the Wiedemann-Franz law. Phys. Rev. B 2013, 88, 125107. [CrossRef]
  21. Yoshino, H.; Murata, K. Significant Enhancement of Electronic Thermal Conductivity of Two-Dimensional Zero-Gap Systems by Bipolar-Diffusion Effect. J. Phys. Soc. Jpn. 2015, 84, 024601. [CrossRef]
  22. Crossno, J.; Shi, J.K.; Wang, K.; Liu, X.; Harzheim, A.; Lucas, A.; Sachdev, S.; Kim, P.; Taniguchi, T.; Watanabe, K.; Ohki, T.A.; Fong, K.C. Observation of the Dirac fluid and the breakdown of the Wiedemann-Franz law in graphene. Science 2016, 351, 1058. [CrossRef]
  23. Rycerz, A. Wiedemann–Franz Law for Massless Dirac Fermions with Implications for Graphene. Materials 2021, 14, 2704. [CrossRef]
  24. Tu, Y.-T.; Das Sarma, S. Wiedemann-Franz law in graphene. Phys. Rev. B 2023, 107, 085401. [CrossRef]
  25. Katsnelson, M.I. The Physics of Graphene, 2nd ed.; Cambridge University Press: Cambridge, UK, 2020; Chapter 3. [CrossRef]
  26. Katsnelson, M.I.; Novoselov, K.S.; Geim, A.K. Chiral tunnelling and the Klein paradox in graphene. Nature Phys. 2006, 2, 620. [CrossRef]
  27. Tworzydło, J.; Trauzettel, B.; Titov, M.; Rycerz, A.; Beenakker, C.W.J. Sub-Poissonian shot noise in graphene. Phys. Rev. Lett. 2006, 96, 246802. [CrossRef]
  28. Miao, F.; Wijeratne, S.; Zhang, Y.; Coskun, U.; Bao, W.; Lau, C.N. Phase-Coherent Transport in Graphene Quantum Billiards. Science 2007, 317, 5844. [CrossRef]
  29. Sonin, E.B. Charge transport and shot noise in ballistic graphene sheet. Phys. Rev. B 2008, 77, 233408. [CrossRef]
  30. Danneau, R.; Wu, F.; Craciun, M.F.; Russo, S.; Tomi, M.Y.; Salmilehto, J.; Morpurgo, A.F.; Hakonen, P.J. Shot Noise in Ballistic Graphene. Phys. Rev. Lett. 2008, 100, 196802. [CrossRef]
  31. Laitinen, A.; Paraoanu, G.S.; Oksanen, M.; Craciun, M.F.; Russo, S.; Sonin, E.; Hakonen, P. Contact doping, Klein tunneling, and asymmetry of shot noise in suspended graphene. Phys. Rev. B 2016, 93, 115413. [CrossRef]
  32. Nair, R.R.; Blake, P.; Grigorenko, A.N.; Novoselov, K.S.; Booth, T.J.; Stauber, T.; Peres, N.M.R.; Geim, A.K. Fine Structure Constant Defines Visual Transparency of Graphene. Science 2008, 320, 1308. [CrossRef]
  33. Skulason, H.S.; Gaskell, P.E.; Szkopek, T. Optical reflection and transmission properties of exfoliated graphite from a graphene monolayer to several hundred graphene layers. Nanotechnology 2010, 21, 295709. [CrossRef]
  34. Merthe, D.J.; Kresin, V.V. Transparency of graphene and other direct-gap two-dimensional materials. Phys. Rev. B 2016, 94, 205439. [CrossRef]
  35. Rycerz, A; Witkowski, P. Sub-Sharvin conductance and enhanced shot noise in doped graphene. Phys. Rev. B 2021, 104, 165413.
  36. Rycerz, A; Witkowski, P. Theory of sub-Sharvin charge transport in graphene disks. Phys. Rev. B 2022, 106, 155428.
  37. Sharvin, Yu.V. A possible method for studying Fermi surfaces. Zh. Eksp. Teor. Fiz. 1965, 48, 984 [Sov. Phys. JETP 1965, 21, 655].
  38. Beenakker, C.W.J.; van Houten, H. Quantum transport in semiconductor nanostructures. Solid State Phys. 1991, 44, 1.
  39. Wang, C.R.; Lu, W.-S.; Hao, L.; Lee, W.-L.; Lee, T.-K.; Lin, F.; Cheng, I-C.; Chen, J.-Z. Enhanced thermoelectric power in dual-gated bilayer graphene. Phys. Rev. Lett. 2011, 107, 186602. [PubMed]
  40. Suszalski, D.; Rut, G.; Rycerz, A. Lifshitz transition and thermoelectric properties of bilayer graphene. Phys. Rev. B 2018, 97, 125403. [CrossRef]
  41. Suszalski, D.; Rut, G.; Rycerz, A. Thermoelectric properties of gapped bilayer graphene. J. Phys.: Condens. Matter 2019, 31, 415501. [CrossRef]
  42. Zong, P.; Liang, J.; Zhang, P.; Wan, C.; Wang, Y.; Koumoto, K.; Graphene-Based Thermoelectrics. ACS Appl. Energy Mater. 2020, 3, 2224. [CrossRef]
  43. Jayaraman, A.; Hsieh, K.; Ghawri, B.; Mahapatra, P.S.; Watanabe, K.; Taniguchi, T.; Ghosh, A. Evidence of Lifshitz Transition in the Thermoelectric Power of Ultrahigh-Mobility Bilayer Graphene. Nano Lett. 2021, 21, 1221. [CrossRef]
  44. Rycerz, A.; Recher, P.; Wimmer, M. Conformal mapping and shot noise in graphene. Phys. Rev. B 2009, 80, 125417. [CrossRef]
  45. Rycerz, A. Magnetoconductance of the Corbino disk in graphene. Phys. Rev. B 2010, 81, 121404(R). [CrossRef]
  46. Peters, E.C.; Giesbers, A.J.M.; Burghard, M.; Kern, K. Scaling in the quantum Hall regime of graphene Corbino devices, Appl. Phys. Lett. 2014, 104, 203109.
  47. Abdollahipour, B.; Moomivand, E. Magnetopumping current in graphene Corbino pump. Physica E 2017, 86, 204. [CrossRef]
  48. Zeng, Y.; Li, J.I.A.; Dietrich, S.A.; Ghosh, O.M.; Watanabe, K.; Taniguchi, T.; Hone, J.; Dean, C.R. High-Quality Magnetotransport in Graphene Using the Edge-Free Corbino Geometry, Phys. Rev. Lett. 2019, 122, 137701. [CrossRef]
  49. Suszalski, D.; Rut, G.; Rycerz, A. Mesoscopic valley filter in graphene Corbino disk containing a p-n junction. J. Phys. Mater. 2020, 3, 015006. [CrossRef]
  50. Kamada, M.; Gall, V.; Sarkar, J.; Kumar, M.; Laitinen, A.; Gornyi, I.; Hakonen, P. Strong magnetoresistance in a graphene Corbino disk at low magnetic fields. Phys. Rev. B 2021, 104, 115432. [CrossRef]
  51. Rycerz, A.; Suszalski, D. Graphene disk in a solenoid magnetic potential: Aharonov-Bohm effect without a two-slit-like setup. Phys. Rev. B 2020, 101, 245429. [CrossRef]
  52. Bouhlal, A.; Jellal, A.; Mansouri, M. Quantum tunneling in graphene Corbino disk in a solenoid magnetic potential with wedge disclination. Physica B 2022, 639, 413904. [CrossRef]
  53. Kumar, C.; Birkbeck, J.; Sulpizio, J.A.; Perello, D.; Taniguchi, T.; Watanabe, K.; Reuven, O.; Scaffidi, T.; Stern, A.; Geim, A.K.; Ilani, S. Imaging Hydrodynamic Electrons Flowing Without Landauer-Sharvin Resistance. Nature 2022, 609, 276. [CrossRef]
  54. Rycerz, A.; Rycerz, K.; Witkowski, P. Thermoelectric Properties of the Corbino Disk in Graphene. Materials 2023, 16, 4250. [CrossRef]
  55. For the forthcoming numerical calculations, we set (in the physical units) ℏ vF=0.575214eV·nm, and e/(πℏ)=2067.83T-1·nm-2.
  56. Recher, P.; Nilsson, J.; Burkard, G.; Trauzettel, B. Bound states and magnetic field induced valley splitting in gate-tunable graphene quantum dots. Phys. Rev. B 2009, 79, 085407. [CrossRef]
  57. Abramowitz, M.; Stegun, I.A. (eds.), Handbook of Mathematical Functions; Dover Publications, Inc.: New York, 1965; Chapter 13.
  58. An important step in the derivation of Eq. (uid15) is the recognition of the Wronskian of Hankel functions, cf. Ref. [57]; Chapter 9, Eq. 9.1.17.
  59. Nazarov, Yu.V.; Blanter, Ya.M. Quantum Transport: Introduction to Nanoscience; Cambridge University Press: Cambridge, UK, 2009; Chapter 1.
  60. Katsnelson, M.I. Aharonov-Bohm effect in undoped graphene: Magnetotransport via evanescent waves. Europhys. Lett. 2010, 89, 17001. [CrossRef]
  61. Rycerz, A. Aharonov-Bohm and Relativistic Corbino Effects in Graphene: A Comparative Study of Two Quantum Interference Phenomena. Acta Phys. Polon. A 2012, 121, 1242. [CrossRef]
  62. Datta, S. Electronic Transport in Mesoscopic Systems; Cambridge University Press: Cambridge, UK, 1997; Chap. 3, p. 129. [CrossRef]
  63. We use the identity I ( a , b ) = 1 2 π ∫ - π π du a + b cos u = 1 a 2 − b 2 , for a > | b | , and first derivative of the above over a; see also Gradshteyn, I.S.; Ryzhik, I.M. Table of Integrals, Series, and Products, Seventh Edition; Academic Press: New York, 2007; Eq. 2.553.3.
  64. Giovannetti, G.; Khomyakov, P.A.; Brocks, G.; Karpan, V.M.; van den Brink, J.; Kelly, P.J. Doping Graphene with Metal Contacts. Phys. Rev. Lett. 2008, 101, 026803. [CrossRef]
  65. Cusati, T.; Fiori, G.; Gahoi, A.; Passi, V.; Lemme, M.C.; Fortunelli, A.; Iannaccone, G. Electrical properties of graphene-metal contacts. Sci. Rep. 2017, 7, 5109. [CrossRef]
  66. Paraoanu, G.S. Klein tunneling through the trapezoidal potential barrier in graphene: conductance and shot noise. New J. Phys. 2021, 23, 043027. [CrossRef]
  67. Anderson, E.; Bai, Z.; Bischof, C.; Blackford, S.; Demmel, J.; Dongarra, J.; et al. LAPACK Users’ Guide, 3rd ed.; Society for Industrial and Applied Mathematics: Philadelphia, USA, 1999.
Preprints 104583 g001
Figure 1. (a) Schematic of Corbino disk in graphene, with the inner radius R i and the outer radius R o , contacted by two circular electrodes (dark areas). A voltage source drives a current through the disk. A separate gate electrode (not shown) allows us to tune the carrier concentration around the neutrality point. The coordinate system ( x , y , z ) is also shown. (b) Cross section of the electrostatic potential profile given by Eq. (2) with m → ∞ (i.e., the rectangular barrier) at y = z = 0 . (c) Zoom-in of a single barrier, for x > 0 , showing also the profiles for m = 2 and 8, with symbolic representations of the incident and reflected waves in inner electrode ( x < R i ) and the transmitted wave in outer electrode ( x > R o ) with the amplitudes r and t corresponding to the Fermi energy E > 0 . (d)–(f) Characteristic values of the magnetic field B = ( 0 , 0 , B ) separating different transport regimes. At B = B c , 1 , the cyclotron diameter 2 r c = R o + R i , and the particle leaving the inner lead approaches the outer lead regardless the initial direction (d). At B = B c , 2 , we have 2 r c = R o − R i , and only the trajectory tangent to the inner lead reaches the outer lead (e). For higher fields, classical trajectories do not contribute to the charge transport, which is possibly only if the resonance with Landau level occurs for E ≈ E n LL , with n = 0 , ± 1 , ± 2 , ⋯ (f).
Figure 1. (a) Schematic of Corbino disk in graphene, with the inner radius R i and the outer radius R o , contacted by two circular electrodes (dark areas). A voltage source drives a current through the disk. A separate gate electrode (not shown) allows us to tune the carrier concentration around the neutrality point. The coordinate system ( x , y , z ) is also shown. (b) Cross section of the electrostatic potential profile given by Eq. (2) with m → ∞ (i.e., the rectangular barrier) at y = z = 0 . (c) Zoom-in of a single barrier, for x > 0 , showing also the profiles for m = 2 and 8, with symbolic representations of the incident and reflected waves in inner electrode ( x < R i ) and the transmitted wave in outer electrode ( x > R o ) with the amplitudes r and t corresponding to the Fermi energy E > 0 . (d)–(f) Characteristic values of the magnetic field B = ( 0 , 0 , B ) separating different transport regimes. At B = B c , 1 , the cyclotron diameter 2 r c = R o + R i , and the particle leaving the inner lead approaches the outer lead regardless the initial direction (d). At B = B c , 2 , we have 2 r c = R o − R i , and only the trajectory tangent to the inner lead reaches the outer lead (e). For higher fields, classical trajectories do not contribute to the charge transport, which is possibly only if the resonance with Landau level occurs for E ≈ E n LL , with n = 0 , ± 1 , ± 2 , ⋯ (f).
Preprints 104583 g001
Preprints 104583 g002
Figure 2. Propagation between consecutive scatterings on interfaces at r = R i and r = R o in a uniform magnetic field B c , 1 < B < B c , 2 . A zoom-in shown an arc of single cyclotron orbit centered at r = r x , with its radii r c , and incident angles θ 1 (for r = R i ) and θ 2 (for r = R o ).
Figure 2. Propagation between consecutive scatterings on interfaces at r = R i and r = R o in a uniform magnetic field B c , 1 < B < B c , 2 . A zoom-in shown an arc of single cyclotron orbit centered at r = r x , with its radii r c , and incident angles θ 1 (for r = R i ) and θ 2 (for r = R o ).
Preprints 104583 g002
Preprints 104583 g003
Figure 3. (a,b) Conductance and (c,d) the Fano factor for the Corbino disk in graphene with the radii R o = 2 R i = 1000 nm and the rectangular potential barrier (i.e., V 0 → ∞ and m → ∞ in Eq. (2)) displayed as functions of the chemical potential. The values of magnetic field are B = 0 (red solid lines in all plots), B = 0 . 2 T (green solid lines), and B = 0 . 4 T (blue solid lines). Inset in (a) is a zoom-in, with black dashed lines depicting the incoherent conductance, see Eq. (25). (a) and (c) show the linear-response results, see Eqs. (18) and (19); the datasets in (b) and (d) are obtained from Eqs. (16) and (17) with U eff = 0 . 01 V. Remaining lines in (c,d) [black solid, black dotted, and black dashed] mark the incoherent Fano factor, see Eq. (27); the values of magnetic field are specified for lines in (c), and are the same in (d). (For B = 0 , horizontal lines mark F incoh ( B → 0 ) = 0 . 111074 corresponding to R o = 2 R i , see Table 1.)
Figure 3. (a,b) Conductance and (c,d) the Fano factor for the Corbino disk in graphene with the radii R o = 2 R i = 1000 nm and the rectangular potential barrier (i.e., V 0 → ∞ and m → ∞ in Eq. (2)) displayed as functions of the chemical potential. The values of magnetic field are B = 0 (red solid lines in all plots), B = 0 . 2 T (green solid lines), and B = 0 . 4 T (blue solid lines). Inset in (a) is a zoom-in, with black dashed lines depicting the incoherent conductance, see Eq. (25). (a) and (c) show the linear-response results, see Eqs. (18) and (19); the datasets in (b) and (d) are obtained from Eqs. (16) and (17) with U eff = 0 . 01 V. Remaining lines in (c,d) [black solid, black dotted, and black dashed] mark the incoherent Fano factor, see Eq. (27); the values of magnetic field are specified for lines in (c), and are the same in (d). (For B = 0 , horizontal lines mark F incoh ( B → 0 ) = 0 . 111074 corresponding to R o = 2 R i , see Table 1.)
Preprints 104583 g003
Preprints 104583 g004
Figure 4. (a,b) Minimal conductivity and maximal Fano factor, see Eqs. (16) and (17), corresponding to the chemical potential fixed at μ = − e U eff / 2 , versus the effective voltage. The magnetic field is specified for each line. (c) The activation voltage, defined via G min ( U on ( 1 ) ) = g 0 [blue open circles], G min ( U on ( 2 ) ) = 2 g 0 [red solid squares], or obtained from scaling according to Eq. (46) [green crosses], displayed versus the magnetic field. (d) The Fano factor corresponding to U eff = U on shown in (c). Horizontal dashed lines in (b,d) mark the value of F incoh ( B → B c , 2 − ) = 0 . 557898 for R o = 2 R i (see Table 1). Dashed line in (c) depicts the approximation given in Eq. (45). The remaining system parameters are same as in Figure 3.
Figure 4. (a,b) Minimal conductivity and maximal Fano factor, see Eqs. (16) and (17), corresponding to the chemical potential fixed at μ = − e U eff / 2 , versus the effective voltage. The magnetic field is specified for each line. (c) The activation voltage, defined via G min ( U on ( 1 ) ) = g 0 [blue open circles], G min ( U on ( 2 ) ) = 2 g 0 [red solid squares], or obtained from scaling according to Eq. (46) [green crosses], displayed versus the magnetic field. (d) The Fano factor corresponding to U eff = U on shown in (c). Horizontal dashed lines in (b,d) mark the value of F incoh ( B → B c , 2 − ) = 0 . 557898 for R o = 2 R i (see Table 1). Dashed line in (c) depicts the approximation given in Eq. (45). The remaining system parameters are same as in Figure 3.
Preprints 104583 g004
Preprints 104583 g005
Figure 5. (a) Linear-response conductance as a function of the chemical potential and (b) finite-voltage conductance, for μ = − e U eff / 2 , as a function of the voltage. The magnetic field is B = 0 . 2 T for all plots. The disk radii are the same as in Figure 3, but the barrier height, see Eq. (2)), is now fixed at V 0 = t 0 / 2 = 1 . 35 eV; the parameter m is specified for each line. (c,d) Zoom-in, for low energies, with same datasets as in (a,b). Datapoints in (b,d) mark the values of G ( U on ( i ) ) = i g 0 , i = 1 , 2 , defining the activation voltages U eff = U on ( i ) .
Figure 5. (a) Linear-response conductance as a function of the chemical potential and (b) finite-voltage conductance, for μ = − e U eff / 2 , as a function of the voltage. The magnetic field is B = 0 . 2 T for all plots. The disk radii are the same as in Figure 3, but the barrier height, see Eq. (2)), is now fixed at V 0 = t 0 / 2 = 1 . 35 eV; the parameter m is specified for each line. (c,d) Zoom-in, for low energies, with same datasets as in (a,b). Datapoints in (b,d) mark the values of G ( U on ( i ) ) = i g 0 , i = 1 , 2 , defining the activation voltages U eff = U on ( i ) .
Preprints 104583 g005
Preprints 104583 g006
Figure 6. Left: Linear-response Fano factor as a function of the chemical potential. Right: Finite-voltage Fano factor, for μ = − e U eff / 2 , as a function of the voltage. The magnetic field is B = 0 . 2 T for all plots, the value of exponent m is specified at each plot, and remaining parameters are same as in Figure 5. Horizontal line at each plot marks the value of F incoh ( B → B c , 2 − ) = 0 . 557898 , see Table 1. Datapoints (right) mark the values of F ( U on ( i ) ) , i = 1 , 2 , corresponding to activation voltages U eff = U on ( i ) , for which G ( U on ( i ) ) = i g 0 (see also Figure 5).
Figure 6. Left: Linear-response Fano factor as a function of the chemical potential. Right: Finite-voltage Fano factor, for μ = − e U eff / 2 , as a function of the voltage. The magnetic field is B = 0 . 2 T for all plots, the value of exponent m is specified at each plot, and remaining parameters are same as in Figure 5. Horizontal line at each plot marks the value of F incoh ( B → B c , 2 − ) = 0 . 557898 , see Table 1. Datapoints (right) mark the values of F ( U on ( i ) ) , i = 1 , 2 , corresponding to activation voltages U eff = U on ( i ) , for which G ( U on ( i ) ) = i g 0 (see also Figure 5).
Preprints 104583 g006
Preprints 104583 g007
Figure 7. (a–d) The activation voltage (for the definition, see Figure 4) and (e–h) the corresponding Fano factor for μ = − e U eff / 2 , displayed as functions of the magnetic field [datapoints]. The value of exponent m is specified at each plot; remaining parameters are same as in Figure 5. Lines in (a–d) depict the approximation given by Eq. (47) [purple solid] and Eq. (45) [black dashed] coinciding in the m → ∞ limit. Horizontal lines in (e–h) mark the value of F incoh ( B → B c , 2 − ) = 0 . 557898 , see Table 1.
Figure 7. (a–d) The activation voltage (for the definition, see Figure 4) and (e–h) the corresponding Fano factor for μ = − e U eff / 2 , displayed as functions of the magnetic field [datapoints]. The value of exponent m is specified at each plot; remaining parameters are same as in Figure 5. Lines in (a–d) depict the approximation given by Eq. (47) [purple solid] and Eq. (45) [black dashed] coinciding in the m → ∞ limit. Horizontal lines in (e–h) mark the value of F incoh ( B → B c , 2 − ) = 0 . 557898 , see Table 1.
Preprints 104583 g007
Table 1. Selected numerical values of F incoh ( B → 0 ) , see Eq. (36), and F incoh ( B → B c , 2 − ) , see Eq. (44). Box marks the values for a = 0 . 5 (i.e., R o = 2 R i ) to be compared with the results following from numerical simulations of quantum transport presented Sec. Section 4.
Table 1. Selected numerical values of F incoh ( B → 0 ) , see Eq. (36), and F incoh ( B → B c , 2 − ) , see Eq. (44). Box marks the values for a = 0 . 5 (i.e., R o = 2 R i ) to be compared with the results following from numerical simulations of quantum transport presented Sec. Section 4.
Preprints 104583 i001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

facebook logotwitter logolinkedin logo
MDPI Initiatives
Important Links
Subscribe

Disclaimer

Terms of Use

Privacy Policy

© 2024 MDPI (Basel, Switzerland) unless otherwise stated