In condensed matter physics, Hofstadter's butterfly describes the spectral properties of non-interacting two-dimensional electrons in a magnetic field in a lattice. The fractal, self-similar nature of the spectrum was discovered in the 1976 Ph.D. work of Douglas Hofstadter and is one of the early examples of computer graphics. The name reflects the visual resemblance of the figure on the right to a swarm of butterflies flying to infinity.
Hofstadter described the structure in 1976 in an article on the energy levels of Bloch electrons in magnetic fields. It gives a graphical representation of the spectrum of Harper's equation at different frequencies. The intricate mathematical structure of this spectrum was independently discovered by Soviet physicist Mark Azbel in 1964 (the Azbel-Hofstadter model), but Azbel did not plot the structure as a geometrical object.
Written while Hofstadter was at the University of Oregon, his paper was influential in directing further research. It predicted on theoretical grounds that the allowed energy level values of an electron in a two-dimensional square lattice, as a function of a magnetic field applied to the system, formed what is now known as a fractal set. That is, the distribution of energy levels for small scale changes in the applied magnetic field recursively repeat patterns seen in the large-scale structure. "Gplot", as Hofstadter called the figure, was described as a recursive structure in his 1976 article in Physical Review B, written before Benoit Mandelbrot's newly coined word "fractal" was introduced in an English text. Hofstadter also discusses the figure in his 1979 book Gödel, Escher, Bach. The structure became generally known as "Hofstadter's butterfly".
In 1997 the Hofstadter butterfly was reproduced in experiments with microwave guide equipped by an array of scatterers. Similarity between the mathematical description of the microwave guide with scatterers and Bloch's waves in magnetic field allowed the reproduction of the Hofstadter butterfly for periodic sequences of the scatterers.
In 2001, Christian Albrecht, Klaus von Klitzing and coworkers realized an experimental setup to test Thouless et al.'s predictions about Hofstadter's butterfly with a two-dimensional electron gas in a supperlattice potential.
In 2013, three separate groups of researchers independently reported evidence of the Hofstadter butterfly spectrum in graphene devices fabricated on hexagonal boron nitride substrates. In this instance the butterfly spectrum results from interplay between the applied magnetic field and the large scale moiré pattern that develops when the graphene lattice is oriented with near zero-angle mismatch to the boron nitride.
In September 2017, John Martinis’s group at Google, in collaboration with the Angelakis group at CQT Singapore, published results from a simulation of 2D electrons in a magnetic field using interacting photons in 9 superconducting qubits. The simulation recovered Hofstadter's butterfly, as expected.
In his original paper, Hofstadter considers the following derivation: a charged quantum particle in a two-dimensional square lattice, with a lattice spacing , is described by a periodic Schrödinger equation, under a static homogeneous magnetic field restricted to a single Bloch band. For a 2D square lattice, the tight binding energy dispersion relation is
where is the energy function, is the crystal momentum, and is an empirical parameter. The magnetic field , where the magnetic vector potential, can be taken into account by using Peierls substitution, replacing the crystal momentum with the canonical momentum , where is the particle momentum operator and is the charge of the particle ( for the electron, is the elementary charge). For convenience we choose the gauge .
Using that is the translation operator, so that , where and is the particle's two-dimensional wave function. One can use as an effective Hamiltonian to obtain the following time-independent Schrödinger equation:
Considering that the particle can only hop between points in the lattice, we write , where are integers. Hofstadter makes the following ansatz: , where depends on the energy, in order to obtain Harper's equation (also known as almost Mathieu operator for ):
Hofstadter's butterfly is the resulting plot of as a function of the flux ratio , where is the set of all possible that are a solution to Harper's equation.
Solutions to Harper's equation and Wannier treatment
Due to the cosine function's properties, the pattern is periodic on with period 1 (it repeats for each quantum flux per unit cell). The graph in the region of between 0 and 1 has reflection symmetry in the lines and . Note that is necessarily bounded between -4 and 4.
Harper's equation has the particular property that the solutions depend on the rationality of . By imposing periodicity over , one can show that if (a rational number), where and are distinct prime numbers, there are exactly energy bands. For large , the energy bands converge to thin energy bands corresponding to the Landau levels.
where and are integers, and is the density of states at a given . Here counts the number of states up to the Fermi energy, and corresponds to the levels of the completely filled band (from to ). This equation characterizes all the solutions of Harper's equation. Most importantly, one can derive that when is an irrational number, there are infinitely many solution for .
The union of all forms a self-similar fractal that is discontinuous between rational and irrational values of . This discontinuity is nonphysical, and continuity is recovered for a finite uncertainty in  or for lattices of finite size. The scale at which the butterfly can be resolved in a real experiment depends on the system's specific conditions.
Phase diagram, conductance and topology
This section needs expansion. You can help by adding to it. (September 2020)
The phase diagram of electrons in a two-dimensional square lattice, as a function of magnetic field, chemical potential and temperature, has infinitely many phases. Thouless and coworkers showed that each phase is characterized by an integral Hall conductance, where all integer values are allowed. These integers are known as Chern numbers.
- Hofstadter, Douglas R. (1976). "Energy levels and wavefunctions of Bloch electrons in rational and irrational magnetic fields". Physical Review B. 14 (6): 2239–2249. Bibcode:1976PhRvB..14.2239H. doi:10.1103/PhysRevB.14.2239.
- Avron J, Osadchy D., and Seiler R. (2003). "A topological look at the quantum Hall effect". Physics Today. 53: 38. doi:10.1063/1.1611351.
- Harper, P.G. (1955). "Scaling analysis of quasiperiodic systems: Generalized harper model". Proceedings of the Physical Society. 68: 874.
- Azbel', Mark Ya. (1964). "Energy Spectrum of a Conduction Electron in a Magnetic Field". Journal of Experimental and Theoretical Physics. 19 (3): 634–645.
- Thouless D. , Kohmoto M, Nightngale and M. den-Nijs (1982). "Quantized Hall conductance in a two dimensional periodic potential". Physical Review Letters. 49 (6): 405–408. Bibcode:1982PhRvL..49..405T. doi:10.1103/PhysRevLett.49.405.
- Kuhl, U.; Stöckmann, H.-J. (13 April 1998). "Microwave realization of the Hofstadter butterfly". Physical Review Letters. 80 (15): 3232–3235. Bibcode:1998PhRvL..80.3232K. doi:10.1103/PhysRevLett.80.3232.
- Albrecht, C.; Smet, J. H.; von Klitzing, K.; Weiss, D.; Umansky, V.; Schweizer, H. (2001-01-01). "Evidence of Hofstadter's Fractal Energy Spectrum in the Quantized Hall Conductance". Physical Review Letters. 86 (1): 147–150. doi:10.1103/PhysRevLett.86.147. ISSN 0031-9007.
- Dean, C. R.; Wang, L.; Maher, P.; Forsythe, C.; Ghahari, F.; Gao, Y.; Katoch, J.; Ishigami, M.; Moon, P.; Koshino, M.; Taniguchi, T.; Watanabe, K.; Shepard, K. L.; Hone, J.; Kim, P. (30 May 2013). "Hofstadter's butterfly and the fractal quantum Hall effect in moiré superlattices". Nature. 497 (7451): 598–602. arXiv:1212.4783. Bibcode:2013Natur.497..598D. doi:10.1038/nature12186. PMID 23676673.
- Ponomarenko, L. A.; Gorbachev, R. V.; Yu, G. L.; Elias, D. C.; Jalil, R.; Patel, A. A.; Mishchenko, A.; Mayorov, A. S.; Woods, C. R.; Wallbank, J. R.; Mucha-Kruczynski, M.; Piot, B. A.; Potemski, M.; Grigorieva, I. V.; Novoselov, K. S.; Guinea, F.; Fal’ko, V. I.; Geim, A. K. (30 May 2013). "Cloning of Dirac fermions in graphene superlattices". Nature. 497 (7451): 594–597. arXiv:1212.5012. Bibcode:2013Natur.497..594P. doi:10.1038/nature12187. hdl:10261/93894. PMID 23676678.
- Hunt, B.; Sanchez-Yamagishi, J. D.; Young, A. F.; Yankowitz, M.; LeRoy, B. J.; Watanabe, K.; Taniguchi, T.; Moon, P.; Koshino, M.; Jarillo-Herrero, P.; Ashoori, R. C. (2013). "Massive Dirac fermions and Hofstadter butterfly in a van der Waals heterostructure". Science. 340 (6139): 1427–1430. arXiv:1303.6942. Bibcode:2013Sci...340.1427H. doi:10.1126/science.1237240. PMID 23686343.
- Roushan, P.; Neill, C.; Tangpanitanon, J.; Bastidas, V. M.; Megrant, A.; Barends, R.; Chen, Y.; Chen, Z.; Chiaro, B.; Dunsworth, A.; Fowler, A.; Foxen, B.; Giustina, M.; Jeffrey, E.; Kelly, J.; Lucero, E.; Mutus, J.; Neeley, M.; Quintana, C.; Sank, D.; Vainsencher, A.; Wenner, J.; White, T.; Neven, H.; Angelakis, D. G.; Martinis, J. (2017-12-01) [2017-09-20]. "Spectroscopic signatures of localization with interacting photons in superconducting qubits" [Spectral signatures of many-body localization with interacting photons]. Science. 358 (6367): 1175–1179. arXiv:1709.07108. doi:10.1126/science.aao1401. ISSN 0036-8075. PMID 29191906.
- Wannier, G. H. (1978-08-01). "A Result Not Dependent on Rationality for Bloch Electrons in a Magnetic Field". Physica Status Solidi (b). 88 (2): 757–765. doi:10.1002/pssb.2220880243.
- Analytis, James G.; Blundell, Stephen J.; Ardavan, Arzhang (May 2004). "Landau levels, molecular orbitals, and the Hofstadter butterfly in finite systems". American Journal of Physics. 72 (5): 613–618. doi:10.1119/1.1615568. ISSN 0002-9505.