The energy spectrum of electrons on a square lattice under a uniform magnetic field B 0 and a twodimensional magnetic modulation B 1 is calculated. When the flux of B 0 is a rational multiple p/q of the fluxon, we find that the Bloch band is broken into q subbands for even q and into 2q subbands for odd q. Symmetry of the modified Hofstadter spectrum is discussed. ͓S0163-1829͑97͒04239-2͔The electronic structure of tight-binding electrons in two dimensions under a uniform magnetic field has been studied for many years.1-5 The spectrum shows a rich behavior such as the Hofstadter's butterfly, but so far it remains an elusive theoretical result for the experimentalists as an extremely high magnetic field in the order of 10 9 G is required to verify the predictions in a crystal sample. However, recent advances in submicronmeter techniques have made it possible to fabricate a lateral surface superlattice with period of the order of 100 nm using a two-dimensional electron gas ͑2DEG͒ in GaAs/Al x Ga 1Ϫx As heterostructures. In these 2DEG's, one observes interesting magnetotransport phenomena even at moderate magnetic fields, leading to a renewal of interest in the subject.6-12 Some recent highlights are the Weiss oscillations, 6 which is an oscillatory magnetoresistance observed in weak magnetic fields, the transport experiments that reveal some signs of the butterfly 13-15 in a twodimensional electric potential modulation system, and the symmetry breakings found in the theoretical studies of the energy spectrum of Bloch electrons under one-dimensional magnetic modulation by Gumbs et al. 16 To pursue the signitures of the Hofstadter spectrum, we study the energy spectrum of Bloch electrons under two-dimensional magnetic flux modulation with a magnetic field B ជ ϭ͓B 0 ϩ(Ϫ) mϪn B 1 ͔ẑ going through the (m,n) plaquette where r ជ ϭmax ϩnaŷ also labels the lower left-hand corner of the plaquette. Here B 0 is the uniform field and B 1 is the two-dimensional modulating field. This 2D modulating field modifies the Harper equation by the addition of a fluxdependent phase factor. In contrast to one-dimensional modulation, there is no symmetry breaking in the energy spectrum. Furthermore, we find interesting difference in the spectrum for specific values of the uniform field. Consider a charge particle hopping on a square lattice, with hopping amplitude t in the presence of an external twodimensional magnetic field B ជ . Let ͉m,n͘ be the Wannier state localized at site (m,n). The tight-binding Hamiltonian is͉͑m,n͘e i͑ea/បc͒A x ͑m,n͒ ͗mϩ1,n͉ϩ͉m,n͘ ϫe i͑ea/បc͒A y ͑m,n͒ ͗m,nϩ1͉ϩH.c.͒.
͑1͒The vector potential A j (m,n)( jϭx,y) resides on the links. The total flux going through any individual plaquette is ⌽ϭ ͚A j aϭBa 2 ϭa 2 B 0 ϩa 2 B 1 (Ϫ1) (mϪn) . For convenience, we choose the gaugeand define the parameters ␣ϵB 0 a 2 / 0 and 2ϵB 1 a 2 / 0 . Here we consider only rational ␣ϭp/q, so that flux for B 0 through a plaquette is a rational fraction of the flux quantum 0 ϭhc/e. To analyze the symmetry we introduce two translations operat...