PHYSICAL REVIEW B VOLUME 58, NUMBER 16 Numerical study of localization in antidot lattices 15 OCTOBER 1998-II Seiji Uryu and Tsuneya Ando Institute for Solid State Physics, University of Tokyo, 7-22-1 Roppongi, Minato-ku, Tokyo 106-8666, Japan Received 12 June 1998 Localization effects in antidot lattices in weak magnetic fields are numerically studied with the use of a Thouless-number method. In hexagonal antidot lattices, both conductance and inverse localization length oscillate as a function of a magnetic flux with the same period as an Al tshuler-aronov-spivak oscillation, in qualitative agreement with recent experiments. S0163-1829 98 01440-4 I. INTRODUCTION An antidot lattice is a two-dimensional electron system modulated by a strong repulsive periodic potential. Recently, a metal-insulator transition was observed in hexagonal antidot lattices in magnetic fields. 1 The temperature dependence of the resistivity exhibits Mott s variable-range hopping and its analysis gave an interesting conclusion that the localization length itself oscillates with the period 0 /2, a half of the flux quantum 0 ch/e, as a function of the total magnetic flux passing through a unit cell. The purpose of this paper is to study the Anderson localization in antidot lattices in a weak magnetic field. Studies of localization effects in antidot lattices began with observation of Al tshuler-aronov-spivak AAS oscillation 2 with the period 0 /2 as a function of in a weak localization regime. 3 5 Numerical calculations revealed that fluctuations in the antidot potential play an essential role in suppressing irregular Aharonov-Bohm oscillations and giving a regular AAS oscillation in realistic antidot lattices where the mean-free path is much larger than a lattice constant. 6,7 In this paper, we shall numerically study the localization effect in antidot lattices based on a Thouless-number method. Difficulties in theoretical treatment of antidot lattices with a large size and realistic potential can be overcome by a scattering S matrix method proposed previously. 8 In Sec. II, the models and the method of numerical calculations are described. Results of explicit calculations are presented in Sec. III and discussed in Sec. IV. A summary is given in Sec. V. II. MODEL AND METHOD In this paper, we shall use a conventional Thoulessnumber method in order to study localization effects. 9,10 This method is useful also in magnetic fields. 11 According to a scaling argument, 12 the Thouless number g(l) is a dimensionless conductance defined by the ratio of the strength of effective coupling V(L) and energy difference W(L) of two systems with size L when the systems are combined, i.e., g L V L W L. 2.1 There can be various ways to estimate the effective coupling strength V(L) and the average energy difference W(L). In the present work, the effective coupling strength is estimated from a geometric average of the curvature of energy bands 2 E/ k 2 for periodic systems having a unit cell with size L and the energy difference is determined from the average energy spacing given by D(E) 1, where D(E) is the density of states at energy E. We then have g L 2 E k 2 D E. 2.2 The energy bands are calculated in an S-matrix method, 8 in which an antidot lattice is replaced by an array of twodimensional quantum-wire junctions characterized by S matrices. It was demonstrated previously 8 that the inclusion of all traveling modes and only a few evanescent modes in the wire region between nearest-neighbor antidots gives accurate energy bands. This reduces the effective matrix size required in the calculation of energy levels considerably and makes calculations in a lattice having a very large unit cell possible. We shall choose quantum-wire junctions illustrated in Fig. 1 a for a square lattice and Fig. 1 b for a hexagonal lattice. The area of a junction for a hexagonal lattice is twice as large as a unit cell. The S matrix for a junction is calculated in a nearest-neighbor tight-binding model based on a recursive Green s-function technique. 13 In the S-matrix method, at a given energy, several k x values are obtained for each k y or several k y values for each k x on an equienergy line in the k space. The second derivative of the energy bands is calculated numerically by varying the energy slightly. The curvature of the bands at each energy for a sample system can be estimated from an average over several such points in the k space. In actual numerical calculations, g(l) is calculated for n n systems with n 4, 8, 12, and 16 for a square lattice and for 2n n systems with n 2, 4, 6, and 8 for a hexagonal lattice. The inverse localization length is determined by fitting the results to g L g 0 exp L, 2.3 with L na and 2na for the square and hexagonal lattice, respectively, where a is a lattice constant. We use the following model potential V(r) of an antidot with its center at the origin: 0163-1829/98/58 16 /10583 6 /$15.00 PRB 58 10 583 1998 The American Physical Society
10 584 SEIJI URYU AND TSUNEYA ANDO PRB 58 is expected that other fluctuations like in the position and shape of an antidot are effectively included by the dot-size fluctuations. We consider antidots with several different diameters d i (i 1,...,n d ), whose distribution is characterized by average d and fluctuation d f, given by d 1 n d i d i, d f 2 1 n d i d i d 2. 2.7 FIG. 1. Schematic illustration of quantum-wire junctions for a a square antidot lattice and b a hexagonal antidot lattice. Dashed lines indicated boundaries of Wigner-Seitz cells. V r U 0 F r, within a Wigner-Seitz cell, where for a square lattice F r cos x a cos y and for a hexagonal lattice 6,7 F r cos a 1 r a 2 cos a 2 r a 2 2 a 2.4, 2.5 cos a 4 /3 1 a 2 r, a 2 2.6 with U 0 being a maximum of the potential, a parameter describing steepness of the potential, a 1 ( 3a/2,a/2), and a 2 (0,a). Note that these potentials vanish identically at the boundary of the Wigner-Seitz cell. The diameter d of the antidot is defined as V(r) E F with r (d/2,0) or (0,d/2) for the square lattice and r (d/2,0)a 1 or (0,d/2)a 2 for the hexagonal lattice, where E F is the Fermi energy. According to self-consistent calculations of energy levels in quantum wires in two-dimensional systems fabricated at a GaAs/Al x Ga 1 x As heterostructure, the electrostatic potential is well approximated by a parabolic potential and the width of the region where the potential increases from the bottom of the Fermi energy is of the same order as the Fermi wavelength for typical electron concentrations. 14,15 This leads to roughly 1 for the aspect ratio d/a 0.5 between an antidot diameter and a lattice constant in the case d 1000 Å, corresponding to usual experiments. Two kinds of disorders are taken into account, those associated with random impurities and fluctuations in the antidot potential. The former is replaced by short-range scatterers and their strength is characterized by the mean free path l e of the two-dimensional system in the absence of the antidot potential. As an example of fluctuations in the antidot potential we shall consider those of the antidot diameter. 6,7 It We denote the average diameter of an antidot d as d for simplicity. In the case of the square antidot lattice, we have to prepare n 4 d junctions because each junction contains four antidots. In the hexagonal case, the number of different kinds of junctions is n 5 d because of the presence of the central antidot. In order to construct a system with size n n, we should combine n 2 junctions in such a way that the antidots of adjacent junctions have the same sizes and the potential is continuous across all boundaries of junctions. Further, we have to prepare several samples with different short-range disorders for each junction in order to include the effects of impurity scattering. For simplicity, we shall use 1 independent of the antidot diameter. This is likely to give channels between neighboring antidots wider than those in actual systems when the aspect ratio is larger, i.e., d/a 0.5. In actual antidot lattices, the potential in the channel region is likely to be raised higher because of the overlapping of the potential of neighboring antidots with the increase of d/a. It is not easy to include effects of such potential overlapping in the hexagonal lattice in particular, because of the necessity to increase the number of different junctions considerably the sample number becomes n 7 d instead of n 5 d ). We consider the case a/ F 3.77, where F is the Fermi wavelength. This corresponds to the electron concentration n s 2.2 10 11 cm 2 for a 2000 Å, which is comparable to n s 1.8 10 11 cm 2 in the experiments. 1 We shall consider two cases d f /a 0.035 and 0.07. This amount of fluctuation is necessary for the explanation 6,7 of the AAS oscillation observed experimentally 4,5 and is likely to be present in actual systems in which the experiments were carried out. 1 We use l e /a 4, which is consistent with a real system with high mobility where the mean free path is much larger than the lattice constant. Further, we have n d 5 for the square lattice and n d 3 for the hexagonal lattice. The number of energy levels used for the determination of the Thouless number is 2000 3000 typically. III. NUMERICAL RESULTS A. Hexagonal antidot lattices Figure 2 shows some examples of the calculated Thouless number g as a function of a magnetic field in hexagonal antidot lattices with size 4 2 (L 4a). They show an AAS oscillation with the period 0 /2 as a function of the total magnetic flux passing through a unit cell given by 3Ba 2 /2, where B is the strength of the magnetic field. The amplitude of the oscillation g normalized by the Thouless number at zero magnetic field g 0 is about g/g 0 0.19 almost independent of d f. This is consistent with the results g/g 0 0.15 calculated for similar systems in a dif-
PRB 58 NUMERICAL STUDY OF LOCALIZATION IN ANTIDOT... 10 585 FIG. 2. Calculated Thouless numbers in hexagonal antidot lattices with d/a 0.7 and the size 4 2. The horizontal axis is the flux 3Ba 2 /2 passing through a unit cell normalized by the flux quantum 0. FIG. 4. Magnetic-field dependence of a the Thouless number and b the inverse localization length in a hexagonal lattice with d/a 0.7 and d f /a 0.035. FIG. 3. Calculated dependence of a the Thouless number and b the inverse localization length on the aspect ratio in the absence of a magnetic field. In a the histograms in the region 0.475 d/a 0.525, for example, indicate the Thouless number for a system with size L 4a, 8a,12a, and 16a from the left to right in antidot lattices with d/a 0.5. Similarly, the histograms in the regions 0.525 d/a 0.575, 0.575 d/a 0.625, 0.625 d/a 0.675, 0.675 d/a 0.725, 0.725 d/a 0.775, and 0.775 d/a 0.825 indicate the system-size dependence for d/a 0.55, 0.6, 0.65, 0.7, 0.75, and 0.8, respectively. The inverse localization length for each aspect ratio is calculated from the linear fitting shown in the figure. In b the horizontal solid line indicates the Fermi energy and the dotted lines show the energy of onedimensional subbands in the quantum wire between neighboring antidots. ferent method, i.e., through a direct calculation of the conductance from the transmission probability. 6,7 Figure 3 shows calculated results of a the Thouless number and b the inverse localization length in antidot lattices with d f /a 0.035 in the absence of a magnetic field. It contains the energy of one-dimensional subband in the quantum-wire region between nearest-neighbor antidots in the absence of the disorder. The figure shows that the localization is quite sensitive to the aspect ratio d/a. In fact, the localization length is reduced from 1 50a for d/a 0.6 to 1 8a for d/a 0.8. Figure 4 shows the magnetic-field dependence of the Thouless number and the inverse localization length for d/a 0.7 and d f /a 0.035. The Thouless number for each system size oscillates with the period of about 0 /2 like an AAS oscillation. The localization length also exhibits a clear oscillation with the same period in good qualitative agreement with the experiments. 1 The localization length oscillates in the range between 20a 70a. Figure 5 shows calculated results for d/a 0.7 and d f /a 0.07. The Thouless number for each system size oscillates with the period of about 0 /2 and at the same time the localization length oscillates with the same period. Due to the presence of larger fluctuations, the localization has become strong in comparison with the results shown in Fig. 4. The amount of the oscillation of the localization length is slightly smaller, i.e., 1 lies in the range between 20a 50a. Figure 6 shows the calculated results for d/a 0.8 and d f /a 0.035. In this case, the channel width between neighboring antidots is small and only a single traveling channel is present see Fig. 3 b. The localization effect is enhanced considerably and the AAS oscillation of the localization length is reduced significantly. In fact, the localization length varies only in the range 8a 10a and the tendency corresponding to a negative magnetoresistance that the localization effect becomes weaker with the increase of the magnetic field can be more clearly identified.
10 586 SEIJI URYU AND TSUNEYA ANDO PRB 58 FIG. 5. Magnetic-field dependence of a the Thouless number and b the inverse localization length in a hexagonal lattice with d/a 0.7 and d f /a 0.07. B. Square antidot lattices FIG. 7. Magnetic-field dependence of a the Thouless number and b the inverse localization length in a square lattice with d/a 0.8 and d f /a 0.035. Figure 7 shows some examples of the results for a square lattice. The parameters are d/a 0.8 and d f /a 0.035, which are the same as those of Fig. 6. In the square lattice the Thouless number is much larger and the inverse localization length is significantly smaller than those of the hexagonal lattice shown in Fig. 6. The Thouless number for each system size oscillates very weakly with the period of about 0 /2 as a function of Ba 2 like an AAS oscillation in agreement with a previous numerical simulation 6,7 and with the experiments. 3 An AAS oscillation cannot clearly be identified in the inverse localization length, while a monotonic decrease of the localization effect with the increase of the magnetic field is much more apparent.. IV. DISCUSSIONS FIG. 6. Magnetic-field dependence of a the Thouless number and b the inverse localization length in a hexagonal lattice with d/a 0.8 and d f /a 0.035. The origin of the AAS oscillation is known to be quantum interferences of a path encircling an antidot with its timereversal path in the weak localization regime. In fact, the quantum correction of the conductivity obtained perturbationally from so-called maximally crossed diagrams takes a maximum when the flux passing through the area is surrounded by the path is an integer multiple of a half of the flux quantum. The AAS oscillation of the conductance does not necessarily mean the presence of the oscillation of the localization length itself, however. The localization depends strongly on the symmetry of the system. In usual systems in the absence of a magnetic field, the Hamiltonian is chosen as a real symmetric matrix and the corresponding wave function is real except for an unimportant phase factor. In this orthogonal case the effective coupling V(L) is given by a real number and an effective dimension of V(L) is given by 1. In the presence of a magnetic field, the Hamiltonian is given by a complex unitary matrix and the corresponding wave function becomes complex. In this unitary case, V(L) is given by a complex number and consequently 2. In the presence of a strong
PRB 58 NUMERICAL STUDY OF LOCALIZATION IN ANTIDOT... 10 587 spin-orbit interaction in the absence of a magnetic field, states are always doubly degenerate Kramers degeneracy and V(L) is given by a quarternion number. In this symplectic case we have 4. This is the exponent characterizing level repulsion effects in the random-matrix theory. 16 It is clear that the localization effect becomes weaker with the increase of. In fact, the localization length in quasi one-dimensional systems is known to be approximately given by N c l e, 17 20 where N c is the number of channels and l e is the mean free path. The AAS oscillation of the conductivity is closely related to the fact that the system belongs to an orthogonal ensemble for n 0 /2 with an integer n and to a unitary ensemble in other cases. This symmetry change can lead to an oscillation of the localization length with a period of 0 /2 in a tight-binding lattice also. However, a numerical study showed a dominant 0 oscillation rather than 0 /2 oscillation. 19 Tight-binding lattices exhibit self-similar spectra called Hofstadter s butterfly 21 and the effective band width itself oscillates as a function of the magnetic flux with period 0. In antidot lattices, on the other hand, the oscillation in the total bandwidth as a function of the flux is not important because the butterfly spectrum is destroyed due to strong mixing among many different bands. This is presumably the reason that an AAS oscillation in the localization length is not appreciable in tight-binding lattices but is clearly observed in antidot lattices. The present results show that the amplitude of localization-length oscillation for strong localization shown in Fig. 6 is much smaller than that for weaker localization shown in Fig. 4. This is consistent with the fact that the AAS oscillation requires interferences between the paths encircling an antidot in different directions. In fact, with the increase in the localization effect the probability amplitude for an electron to circle around an antidot is reduced. In the strong localization limit where the localization length is smaller than the antidot period, in particular, the oscillation amplitude should be reduced considerably because the interference between paths rounding an antidot is almost completely suppressed. The large difference between the square and hexagonal lattices is considered to arise from the nature of the electron motion in these lattices. In fact, electrons are more frequently scattered by antidots and perform more complicated trajectories in hexagonal lattices than in square lattices. As a results, electrons circle around a dot more frequently in hexagonal lattices than in square lattices, leading to a possible enhancement in the modulation of interference effects by a magnetic flux as well as a reduction of the conductance. According to our results on hexagonal antidot lattices with d f /2 0.035, the oscillation amplitude of the inverse localization length normalized by that at zero magnetic field 0 is / 0 0.37 for d/a 0.7. This is in reasonable agreement with experimental results / 0 0.20. 1 However, the absolute magnitude of the localization length obtained in the present calculation is 0 1 20a, which is much larger than 0 1 a estimated experimentally. The localization length can be much smaller for larger aspect ratio, but the amplitude of the AAS oscillation of the localization length becomes smaller with the increase of the localization effect ( / 0 0.11 for d/a 0.8 for which 0 1 8a). When the confinement is strong, the system is considered as a quantum-dot array rather than an antidot lattice. The hopping of electrons to neighboring dots can be considerably suppressed by a Coulomb blockade if effects of electronelectron interactions are considered. In a weak localization regime, effects of electron-electron interactions give another quantum correction to the conductivity. 22 Further, when such a charging effect is important in a strong localization regime, the system exhibits a Coulomb gap. 23 In the case, however, the temperature dependence of hopping resistivity obeys (T) 0 exp (T 0 /T) 1/2 different from Mott s variable-range hopping observed in the experiments. 1 V. SUMMARY The Anderson localization in antidot lattices has been numerically studied with the use of a Thouless-number method. The localization is very sensitive to the aspect ratio between the antidot diameter and the lattice constant. In a hexagonal antidot lattice, the Thouless number oscillates with the period equal to that of the AAS oscillation and the localization length oscillates also with the same period. The AAS oscillation of the localization length has a large amplitude when the localization length is much larger than the lattice constants but becomes smaller in the stronger localization regime. In a square lattice, the AAS oscillation of the Thouless number itself is quite small and that of the localization length has not been identified clearly. ACKNOWLEDGMENTS We would like to thank S. Ishizaka, T. Nakanishi, F. Nihey, K. Asano, and H. Matsumura for useful discussions. This work was supported in part by the Japan Society for the Promotion of Science Research for the Future Program No. JSPS-RFTF96P00103. One of the authors S.U. has been supported by the Japan Society of the Promotion of Science for Young Scientists. Some numerical calculations were performed on FACOM VPP500 in the Supercomputer Center, Institute for Solid State Physics, University of Tokyo. 1 F. Nihey, M. A. Kastner, and K. Nakamura, Phys. Rev. B 55, 4085 1997. 2 B. L. Al tshuler, A. G. Aronov, and B. Z. Spivak, Pis ma Zh. Eksp. Teor. Fiz. 33, 101 1981 JETP Lett. 33, 94 1981. 3 G. M. Gusev, Z. D. Kvon, L. V. Litvin, Yu. V. Nastaushev, A. K. Kalagin, and A. I. Toropov, Pis ma Zh. Eksp. Teor. Fiz. 55, 129 1991 JETP Lett. 55, 123 1992. 4 K. Nakamura, S. Ishizaka, and F. Nihey, Physica B 197, 144 1994. 5 F. Nihey, S. W. Hwang, and K. Nakamura, Phys. Rev. B 51, 4649 1995. 6 T. Nakanishi and T. Ando, Phys. Rev. B 54, 8021 1996.
10 588 SEIJI URYU AND TSUNEYA ANDO PRB 58 7 T. Nakanishi and T. Ando, Physica B 227, 127 1996. 8 S. Uryu and T. Ando, Phys. Rev. B 53, 13 613 1996. 9 J. T. Edwards and D. J. Thouless, J. Phys. C 5, 807 1972. 10 D. C. Licciardello and D. J. Thouless, J. Phys. C 8, 4157 1975. 11 T. Ando, J. Phys. Soc. Jpn. 52, 1740 1983 ; 53, 3101 1984 ; 53, 3126 1984. 12 E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. V. Ramakrishnan, Phys. Rev. Lett. 42, 673 1979. 13 T. Ando, Phys. Rev. B 44, 8017 1991. 14 S. E. Laux, D. J. Frank, and F. Stern, Surf. Sci. 196, 101 1988. 15 T. Suzuki and T. Ando, J. Phys. Soc. Jpn. 62, 2986 1993 ; Physica B 201, 345 1994 ; 227, 46 1996. 16 F. J. Dyson, J. Math. Phys. 3, 140 1962. 17 O. N. Dorokhov, Zh. Eksp. Teor. Fiz. 85, 1040 1983 Sov. Phys. JETP 58, 606 1983. 18 K. B. Efetov and A. I. Larkin, Zh. Eksp. Teor. Fiz. 85, 764 1983 Sov. Phys. JETP 58, 444 1983. 19 J.-L. Pichard, M. Sanquer, K. Slevin, and P. Debray, Phys. Rev. Lett. 65, 1812 1990. 20 H. Tamura, Ph.D. thesis, University of Tokyo, 1992. 21 D. R. Hofstadter, Phys. Rev. B 14, 2239 1976. 22 See, for example, H. Fukuyama, Proc. Theor. Phys. Suppl. 84, 47 1985, and references cited therein. 23 A. L. Efros and B. I. Shklovskii, J. Phys. C 8, L49 1975.