1 Introduction
Universal conductance fluctuation (UCF) is one of the striking features in mesoscopic electronic transport [
1-
3]. It is a quantum interference effect that occurs in systems with sample size being less than the phase coherence length. And the UCF amplitude only depends on the dimensionality and system symmetries [
4-
6]. According to random matrix theory [
7], there are usually three ensembles due to different symmetries: (1) Circular orthogonal ensemble (COE), where both time-reversal and spin-rotation symmetries are preserved; (2) Circular unitary ensemble (CUE), where the time-reversal symmetry is broken but the spin-rotation symmetry ia preserved; (3) Circular symplectic ensemble (CSE), where the time-reversal symmetry is preserved but the spin-rotational symmetry is broken. These three ensembles can be indexed by
, corresponding respectively to COE, CUE and CSE. For quantum dot, quasi-one dimensional, and two dimensional systems, the UCF at the diffusive regimes are respectively
, where
[
4-
7].
So far, there have been tremendous focus on the UCF for the integer dimensions [
8-
22]. However, there is very limited study on the UCF in fractal dimension [
23], which is also named as Hausdorff dimension measuring the dimensionality of self-similar structures [
24]. As the simplest models of regular fractal structure, Sierpinski gasket and carpet have attracted numerous attentions. The related research mainly focused on the electronic transport properties [
25-
34], and the topological properties in Sierpinski lattices [
35]. For example, Schwalm
et al. investigated the conductance distribution in the presence of random bond disorder via length scaling theory and described the curves by two parameters that are independent of the model system [
29]; spin-related transport study of fractal conductors proposed a spin filter based on Sierpinski carpet [
32,
33]; and the nonzero Chern number in Sierpinski carpet under external magnetic field [
35] demonstrates the extraordinary electronic properties in fractal dimension. In addition to the previous theoretical proposals, the realistic experimental realization of Sierpinski triangle fractals either by molecular self-assembling or by atomic manipulation [
36,
37] arises some interesting questions: Does there exist a UCF in fractal dimension? How does the symmetry affect the conductance fluctuation?
Fractals are self-similarity patterns and can be iterated infinitely. According to the characteristics of fractal, it can be classified into different types, e.g., exact/quasi/statistical/qualitative self-similar patterns. As a representative exact self-similar fractal, we choose Sierpinski carpet to address the above questions. We provide a systematic investigation on the electronic transport properties of Sierpinski carpets in the presence of disorders. Unlike integral dimensions, we find that the Sierpinski carpets do not have a UCF at the COE situation; When the time-reversal symmetry or the spin-rotational symmetry is broken, we show that UCF can be reached at the CUE or CSE system, i.e., ; And the corresponding conductance distributions in the diffusive region for both CUE and CSE systems obey nearly the same distribution, i.e., normal distribution at the small conductance region but log-normal distribution at the large conductance region.
2 Model System
In our numerical calculations, the tight-binding model Hamiltonian of square lattice system can be expressed as [
35]:
where is the electron creation (annihilation) operator with spin α at site i, and t is the hopping energy between nearest neighbor sites. The effect of external magnetic field is considered by a Peierls phase factor in the first term, where is vector potential with Landau gauge and is the quantum of flux. The second term describes the removed region of Sierpinski carpet with a large on-site potential , which has been elaborately tested for the convergence of conductance. The third term represents Rashba spin-orbit coupling with coupling strength , σ are the Pauli matrices and dij stands for the unit vector pointing from site j to site i. The last term corresponds to the Anderson type disorder with being uniformly distributed within the range of , where W characterizes the disorder strength.
In our consideration, the two-terminal electronic transport was evaluated by employing the Landauer-Büttiker formula [
36] to calculate the conductance
G, i.e.,
. Here
are the retarded and advanced Green’s functions of the central disordered region, and
are the line-width functions coupling left and right terminals to the central region, respectively. In the presence of finite disorder strength, the conductance fluctuation can be obtained by
, where
represents the ensemble average over different disordered samples at the same strength
W. In our calculation,
t is set to be the unit of energy
E, disorder strength
W and spin-orbit coupling strength
; and the magnetic field is measured by magnetic flux
.
3 Fractal structures and density of states
Figure 1(a) displays the system configurations of Sierpinski carpet at different iteration steps. For example, “R2" indicates the system configuration after the second iteration step [see the middle panel of Fig. 1(a)]. The white parts are set to be empty. Due to the nature of fractal, the structure can be infinitely iterated and the empty regions would become more and more intensive. Figure 1(b) displays density of states as a function of Fermi energy E in the presence of a series of magnetic flux ϕ. For pristine square lattice (i.e., in black line, R0), Landau level splitting gradually appears along with the increase of magnetic flux ϕ. Meanwhile, the interval of Landau level becomes increasing with the increase of ϕ, agreeing well with the theoretical expectations. Compared with that in pristine square lattice system, the density of states of Sierpinski carpet decreases for the whole range of energy. And the more iteration steps of the Sierpinski carpet, the smaller of density of states. At the first iteration step (i.e., the R1 case where the structure is similar to square lattice except for an empty region in the center), the density of states exhibits the similar characteristic as that of the pristine square lattice. However, one can see that the density of states for structures with larger iteration steps are much more robust against weak magnetic field, implying that stronger magnetic field is required to split their energy levels. This means that the special structure of Sierpinski carpet results in distinct electronic properties. Due to the symmetric feature of density of states about the Fermi energy, we choose a series of within the range of 0 to 4 in our numerical calculation.
4 Conductance fluctuations in systems with different symmetries
4.1 Absence of universal conductance fluctuation in COE symmetry
First, We investigate the average conductance and the fluctuation in the Sierpinski carpet with COE symmetry (). Figures 2(a)-2(c) display the average conductance as a function of the disorder strength W at different energies by evaluating over 3000 samples at each point. In the absence of disorder, i.e.,, one can find that decays rapidly as the iteration step increases, e.g., at , the average conductance was decreased to be about 10% from ~63 e2/h to ~6.3 e2/h as the iteration step varies from R1 to R3. This can be attributed to the strong electron back-scattering from the boundaries of the empty regions. In the presence of disorders (i.e., ), the average conductance decreases along with the increase of the disorder strength. Figures. 2(d)-2(f) describe the conductance fluctuation as a function of disorder strength W. As disorder strength W increases, the conductance fluctuation first increases to its maximal plateau, and then gradually decreases. If disorder strength W is strong enough, all electrons in the system would be completely localized to reach the Anderson insulating state, and thus the corresponding conductance fluctuation becomes zero. One can find that the peaks of conductance fluctuation that reached its maximum are energy dependent at fixed iteration step, e.g., at R1 [see Fig. 2(d)], the peak of conductance fluctuation drops from 1.07 e2/h to 0.93 e2/h as the Fermi energy increases from 0.2t to 3.2t. Moreover, the dashed orange line in the figures represents the peak of conductance fluctuation at (R1, ). From this dashed orange line, we can find that the peak tends to decrease as the iteration step increases, e.g., at , the peaks corresponding to R1, R2, and R3 are 1.00 e2/h, 0.95 e2/h, and 0.88 e2/h, respectively. From above discussion, we can find that the peaks of conductance fluctuation are dependent on the iteration step and Fermi energy. Therefore, the conductance fluctuation in the COE () symmetry of Sierpinski carpet does not exhibit a universal feature.
4.2 Universal conductance fluctuation in CUE symmetry
Then, we turn to study the conductance fluctuation of CUE system that breaks time-reversal symmetry (i.e.,
). Figures 3(a)-3(c) display the average conductance as a function of disorder strength
W. In the absence of disorder (i.e.,
), the conductance G at different Fermi energies decreases as iteration step increases. And at fixed magntic flux and Fermi energy, the average conductance also decreases as the disorder strength
W increases, similar to that in the COE system. But at the same system parameters, the average conductance in CUE system is smaller than that in COE system due to the enhancement of localized states induced by external magnetic field. At fixed iteration step, e.g., R
1 [see the lines with (
) of (0.012, 1.4), (0.024, 1.4), and (0.048, 1.4) in Fig. 3(a)], the average conductance quickly decreases as the magnetic flux increases for weak disorder strength
W. When the disorder strength
W exceeds certain critical value, such a dependence weakens and results in an saturation of average conductance. Figures. 3(d)-3(f) display the conductance fluctuation
as a function of disorder strength
W. As illustrated by the dashed orange line, there exists a maximum of
around 0.78 ± 0.01 (
e2/
h) for different system configurations, magnetic fluxes, or Fermi energies, suggesting a universal feature of conductance fluctuation [
37]. Moreover, the critical disorder strength
that reaches the maximum of conductance fluctuation for different Fermi energies exists in a tiny range, and decreases as the iteration step increases, e.g.,
for R
1, R
2, and R
3, respectively. To confirm it is indeed a universal value, we further study the conductance fluctuation
for even larger system size. Figure 5(a) displays the maximum of conductance fluctuation in a wide range of system size from
to
. Due to the finite-size effect, the maximum of conductance fluctuation gradually decreases as the system size increases, and reaches a convergence at larger system size, i.e., the plateau is around 0.74 ± 0.01 (
e2/
h) when the system size is larger than
. Therefore, different from COE system where the fluctuation does not manifests a universal feature, there exists a universal conductance fluctuation 0.74 ± 0.01 (
e2/
h) in CUE system with time-reversal symmetry breaking.
4.3 Universal conductance fluctuations in CSE symmetry
Next, we move to explore the conductance fluctuation of systems with CSE symmetry (i.e., ). Figures 4(a)-(c) display the average conductance as a function of disorder strength W. One can find some universal features of the variation of average conductance in COE, CUE, and CSE systems, i.e., the average conductance decreases with the increase of iteration steps, Fermi energies, or disorder strength. Moreover, the average conductance increases as the spin-orbit coupling strength increases [see the curves with () of (0.2, 1.7), (0.4, 1.7), and (0.6, 1.7) in Figs. 4(b) or (c)]. Figures 4(d)-4(f) display the conductance fluctuation as a function of disorder strength W. The dashed orange line in these figures denotes the maximum of around 0.75 ± 0.01 (e2/h) at different iteration steps, spin-orbit coupling strengths, and Fermi energies. All these together indicate the presence of universal conductance fluctuation in the CSE system. In addition, the critical disorder strength at fixed decreases as the iteration step increases (e.g., in the case of and , , 4.74, and 3.68 for R1, R2, and R3 respectively). By further investigating the maxima of conductance fluctuation in a wide range of system size from to [see Fig. 5(b)], one can find that the system size shows little influence on the conductance fluctuation plateau, which is around 0.74 ± 0.01 (e2/h). Therefore, we conclude that there also exists a universal conductance fluctuation in systems with CSE symmetry. In addition, the difference of the size effect of UCF in CUE and CSE indicates that time-reversal symmetry has significant influence on electronic transport in a disordered system with different size. A phenomenological picture of band structure can be used to explain the difference. In a magnetic field, electronic band will form Landau level and this is a global influence on the band structure. The number of Landau levels depend on the strength of magnetic field, as well as system size. On the contrary, bands structure will open band gaps at the crossing points when Rahsba spin-orbit coupling is considered and the system size have little influence on the electronic band structure. This is a local influence on the band structure. Therefore, UCF in CUE decrease gradually to a saturation as the increase of system size; while it seems unchanged in CSE.
5 Conductance distributions in both CUE and CSE
From above discussions, one can reach that the Sierpinski carpet does not exhibit a universal feature of conductance fluctuations in COE system, but reaches a universal conductance fluctuation of 0.74 ± 0.01 (e2/h) for both CUE (time-reversal symmetry breaking) and CSE (spin-rotation symmetry breaking) symmetries. To further reveal the physical mechanism of the universal conductance fluctuations in Sierpinski carpet, we study the conductance distribution around the critical disorder strength . Figure 6 displays the conductance distribution near the critical disorder strength in CUE system (). One can find that the general function of the conductance distribution is independent of Fermi energies or iteration steps, i.e., the left regime obeys the Gaussian distribution function, while the right regime obeys the log-normal distribution function. In the CSE systems, we also find that the conductance near the critical disorder strength follows the same distribution function as that in the case of CUE symmetry [see Fig. 7]. Therefore, one can conclude that in the Sierpinski carpet, both CUE and CSE systems share the same universal conductance fluctuations and conductance distribution near the critical disorder strength.
As displayed in Fig. 8, we further study the conductance distribution beyond the critical disorder strength
. In CUE system [see Fig. 8(a)-(c)], the conductance distribution belongs to Gaussian function when disorder strength
W less than
[see Fig. 8(a)], indicating the system stays in the ballistic regime. While it follows the log-normal distribution when disorder strength
W is larger than
[see Fig. 8(c)], implying that the system enters the localized regime [
5,
12]. In particular, the conductance distribution experiences a transition from Gaussian to log-normal function around
as mentioned above. In CSE system [see Fig. 8(d)-(f)], the variation of conductance distribution around the critical disorder strength
behaves similar as that in CUE system. Therefore, from the analysis of conductance distribution, we can conclude that the universal conductance fluctuation in Sierpinski carpet occurs in the crossover region between ballistic and Anderson localization regime in both CUE and CSE systems.
6 Summary
In this work, we theoretically investigate the electronic transport properties (e.g., conductance and its fluctuation) of disordered two-terminal device in Sierpinski carpet for circular orthogonal, circular unitary, and circular symplectic ensembles. Our results show that there exist universal conductance fluctuations in both CUE and CSE systems, while not in the COE system [
41]. By further studying the conductance distribution around the critical disorder strength
in CUE and CSE systems, we find that both of systems share the same distribution function, which obeys Gaussian distribution when
or follows log-normal distribution when
, indicating that the universal conductance fluctuations in both CUE and CSE systems occur in the crossover region between ballistic and Anderson localized regimes. In integer dimensions, the universal conductance fluctuations in metal-insulator crossover regime have been studied by Qiao
et al. [
14]. Compared to aforementioned study, universal conductance fluctuations in Sierpinski carpet have different values, indicating the difference between the universal conductance fluctuations in fractal and integer dimensions. Our finding in Sierpinski carpet can be extended to study other fractal structures that have different Hausdorff dimensions.
Higher Education Press and Springer-Verlag GmbH Germany, part of Springer Nature 2019