Physics:Quantum Hall effect

From HandWiki
Short description: Electromagnetic effect in physics

The quantum Hall effect (or integer quantum Hall effect) is a quantized version of the Hall effect which is observed in two-dimensional electron systems subjected to low temperatures and strong magnetic fields, in which the Hall resistance Rxy exhibits steps that take on the quantized values

[math]\displaystyle{ R_{xy} = \frac{V_\text{Hall}}{I_\text{channel}} = \frac{h}{e^2\nu} , }[/math]

where VHall is the Hall voltage, Ichannel is the channel current, e is the elementary charge and h is Planck's constant. The divisor ν can take on either integer (ν = 1, 2, 3,...) or fractional (ν = 1/3, 2/5, 3/7, 2/3, 3/5, 1/5, 2/9, 3/13, 5/2, 12/5,...) values. Here, ν is roughly but not exactly equal to the filling factor of Landau levels. The quantum Hall effect is referred to as the integer or fractional quantum Hall effect depending on whether ν is an integer or fraction, respectively.

The striking feature of the integer quantum Hall effect is the persistence of the quantization (i.e. the Hall plateau) as the electron density is varied. Since the electron density remains constant when the Fermi level is in a clean spectral gap, this situation corresponds to one where the Fermi level is an energy with a finite density of states, though these states are localized (see Anderson localization).[1]

The fractional quantum Hall effect is more complicated; its existence relies fundamentally on electron–electron interactions. The fractional quantum Hall effect is also understood as an integer quantum Hall effect, although not of electrons but of charge-flux composites known as composite fermions. In 1988, it was proposed that there was quantum Hall effect without Landau levels.[2] This quantum Hall effect is referred to as the quantum anomalous Hall (QAH) effect. There is also a new concept of the quantum spin Hall effect which is an analogue of the quantum Hall effect, where spin currents flow instead of charge currents.[3]


The quantization of the Hall conductance ([math]\displaystyle{ G_{xy}= 1/R_{xy} }[/math]) has the important property of being exceedingly precise. Actual measurements of the Hall conductance have been found to be integer or fractional multiples of e2/h to nearly one part in a billion. This phenomenon, referred to as exact quantization, is not really understood but it has sometimes been explained as a very subtle manifestation of the principle of gauge invariance.[4] It has allowed for the definition of a new practical standard for electrical resistance, based on the resistance quantum given by the von Klitzing constant RK. This is named after Klaus von Klitzing, the discoverer of exact quantization. The quantum Hall effect also provides an extremely precise independent determination of the fine-structure constant, a quantity of fundamental importance in quantum electrodynamics.

In 1990, a fixed conventional value RK-90 = 25812.807 Ω was defined for use in resistance calibrations worldwide.[5] On 16 November 2018, the 26th meeting of the General Conference on Weights and Measures decided to fix exact values of h (the Planck constant) and e (the elementary charge),[6] superseding the 1990 value with an exact permanent value RK = h/e2 = 25812.80745... Ω.[7]


The MOSFET (metal-oxide-semiconductor field-effect transistor), invented by Mohamed Atalla and Dawon Kahng at Bell Labs in 1959,[8] enabled physicists to study electron behavior in a nearly ideal two-dimensional gas.[9] In a MOSFET, conduction electrons travel in a thin surface layer, and a "gate" voltage controls the number of charge carriers in this layer. This allows researchers to explore quantum effects by operating high-purity MOSFETs at liquid helium temperatures.[9]

The integer quantization of the Hall conductance was originally predicted by University of Tokyo researchers Tsuneya Ando, Yukio Matsumoto and Yasutada Uemura in 1975, on the basis of an approximate calculation which they themselves did not believe to be true.[10] In 1978, the Gakushuin University researchers Jun-ichi Wakabayashi and Shinji Kawaji subsequently observed the effect in experiments carried out on the inversion layer of MOSFETs.[11]

In 1980, Klaus von Klitzing, working at the high magnetic field laboratory in Grenoble with silicon-based MOSFET samples developed by Michael Pepper and Gerhard Dorda, made the unexpected discovery that the Hall resistance was exactly quantized.[12][9] For this finding, von Klitzing was awarded the 1985 Nobel Prize in Physics. A link between exact quantization and gauge invariance was subsequently proposed by Robert Laughlin, who connected the quantized conductivity to the quantized charge transport in a Thouless charge pump.[4][13] Most integer quantum Hall experiments are now performed on gallium arsenide heterostructures, although many other semiconductor materials can be used. In 2007, the integer quantum Hall effect was reported in graphene at temperatures as high as room temperature,[14] and in the magnesium zinc oxide ZnO–MgxZn1−xO.[15]

Integer quantum Hall effect

Animated graph showing filling of Landau levels as B changes and the corresponding position on a graph of hall coefficient and magnetic field|Illustrative only. The levels spread out with increasing field. Between the levels the quantum hall effect is seen.

Landau levels

Main page: Physics:Landau quantization

In two dimensions, when classical electrons are subjected to a magnetic field they follow circular cyclotron orbits. When the system is treated quantum mechanically, these orbits are quantized. To determine the values of the energy levels the Schrödinger equation must be solved.

Since the system is subjected to a magnetic field, it has to be introduced as an electromagnetic vector potential in the Schrödinger equation. The system considered is an electron gas that is free to move in the x and y directions, but is tightly confined in the z direction. Then, a magnetic field is applied in the z direction and according to the Landau gauge the electromagnetic vector potential is [math]\displaystyle{ \mathbf{A} = (0, Bx, 0) }[/math] and the scalar potential is [math]\displaystyle{ \phi=0 }[/math]. Thus the Schrödinger equation for a particle of charge [math]\displaystyle{ q }[/math] and effective mass [math]\displaystyle{ m^* }[/math] in this system is:

[math]\displaystyle{ \left \{ \frac{1}{2m^*} \left [ \mathbf{p}-q\mathbf{A}\right ]^2 + V(z) \right \}\psi(x,y,z)=\varepsilon\psi(x,y,z) }[/math]

where [math]\displaystyle{ \mathbf{p} }[/math] is the canonical momentum, which is replaced by the operator [math]\displaystyle{ -i\hbar\nabla }[/math] and [math]\displaystyle{ \varepsilon }[/math] is the total energy.

To solve this equation it is possible to separate it into two equations since the magnetic field just affects the movement along x and y axes. The total energy becomes then, the sum of two contributions [math]\displaystyle{ \varepsilon=\varepsilon_z+\varepsilon_{xy} }[/math]. The corresponding equations in z axis is:

[math]\displaystyle{ \left [- \frac{\hbar^2}{2m^*} {\partial^2\over\partial z^2} + V(z) \right ]u(z)=\varepsilon_zu(z) }[/math]

To simplify things, the solution [math]\displaystyle{ V(z) }[/math] is considered as an infinite well. Thus the solutions for the z direction are the energies [math]\displaystyle{ \varepsilon_z=\frac{n_z^2\pi^2\hbar^2}{2m^*L^2} }[/math], [math]\displaystyle{ n_z=1,2,3... }[/math] and the wavefunctions are sinusoidal. For the [math]\displaystyle{ x }[/math] and [math]\displaystyle{ y }[/math] directions, the solution of the Schrödinger equation can be chosen to be the product of a plane wave in [math]\displaystyle{ y }[/math]-direction with some unknown function of [math]\displaystyle{ x }[/math], i.e., [math]\displaystyle{ \psi_{xy}=u(x)e^{ik_yy} }[/math]. This is because the vector potential does not depend on [math]\displaystyle{ y }[/math] and the momentum operator [math]\displaystyle{ \hat p_y }[/math] therefore commutes with the Hamiltonian. By substituting this Ansatz into the Schrödinger equation one gets the one-dimensional harmonic oscillator equation centered at [math]\displaystyle{ x_{k_y}=\frac{\hbar k_y}{eB} }[/math].

[math]\displaystyle{ \left [- \frac{\hbar^2}{2m^*} {\partial^2\over\partial x^2}+\frac{1}{2}m^*\omega_{\rm c}^2(x-l_B^2k_y)^2\right ]u(x)=\varepsilon_{xy}u(x) }[/math]

where [math]\displaystyle{ \omega_{\rm c}=\frac{eB}{m^*} }[/math] is defined as the cyclotron frequency and [math]\displaystyle{ l_B^2=\frac{\hbar}{eB} }[/math] the magnetic length. The energies are:

[math]\displaystyle{ \varepsilon_{xy}\equiv\varepsilon_{n_x}=\hbar \omega_{\rm c} \left ( n_x+\frac{1}{2} \right ) }[/math], [math]\displaystyle{ n_x=1,2,3... }[/math]

And the wavefunctions for the motion in the [math]\displaystyle{ xy }[/math] plane are given by the product of a plane wave in [math]\displaystyle{ y }[/math] and Hermite polynomials attenuated by the gaussian function in [math]\displaystyle{ x }[/math], which are the wavefunctions of a harmonic oscillator.

From the expression for the Landau levels one notices that the energy depends only on [math]\displaystyle{ n_x }[/math], not on [math]\displaystyle{ k_y }[/math]. States with the same [math]\displaystyle{ n_x }[/math] but different [math]\displaystyle{ k_y }[/math] are degenerate.

Density of states

At zero field, the density of states per unit surface for the two-dimensional electron gas taking into account degeneration due to spin is independent of the energy

[math]\displaystyle{ n_{\rm 2D}=\frac{m^*}{\pi \hbar^2} }[/math].

As the field is turned on, the density of states collapses from the constant to a Dirac comb, a series of Dirac [math]\displaystyle{ \delta }[/math] functions, corresponding to the Landau levels separated [math]\displaystyle{ \Delta\varepsilon_{xy}=\hbar \omega_{\rm c} }[/math]. At finite temperature, however, the Landau levels acquire a width [math]\displaystyle{ \Gamma=\frac{\hbar}{\tau_i} }[/math] being [math]\displaystyle{ \tau_i }[/math] the time between scattering events. Commonly it is assumed that the precise shape of Landau levels is a Gaussian or Lorentzian profile.

Another feature is that the wave functions form parallel strips in the [math]\displaystyle{ y }[/math] -direction spaced equally along the [math]\displaystyle{ x }[/math]-axis, along the lines of [math]\displaystyle{ \mathbf{A} }[/math]. Since there is nothing special about any direction in the [math]\displaystyle{ xy }[/math]-plane if the vector potential was differently chosen one should find circular symmetry.

Given a sample of dimensions [math]\displaystyle{ L_x \times L_y }[/math] and applying the periodic boundary conditions in the [math]\displaystyle{ y }[/math]-direction [math]\displaystyle{ k=\frac{2\pi}{L_y}j }[/math] being [math]\displaystyle{ j }[/math] an integer, one gets that each parabolic potential is placed at a value [math]\displaystyle{ x_k=l_B^2k }[/math].

Parabolic potentials along the [math]\displaystyle{ x }[/math]-axis centered at [math]\displaystyle{ x_k }[/math] with the 1st wave functions corresponding to an infinite well confinement in the [math]\displaystyle{ z }[/math] direction. In the [math]\displaystyle{ y }[/math]-direction there are travelling plane waves.

The number of states for each Landau Level and [math]\displaystyle{ k }[/math] can be calculated from the ratio between the total magnetic flux that passes through the sample and the magnetic flux corresponding to a state.

[math]\displaystyle{ N_B=\frac{\phi}{\phi_0}=\frac{BA}{BL_y\Delta x_k}=\frac{A}{2\pi l_B^2}\begin{array}{lcr}& l_B&\\&=&\\ &&\end{array} \frac{AeB}{2\pi \hbar}\begin{array}{lcr}& \omega_{\rm c}&\\&=&\\ &&\end{array} \frac{m^*\omega_{\rm c}A}{2\pi \hbar} }[/math]

Thus the density of states per unit surface is

[math]\displaystyle{ n_B=\frac{m^*\omega_{\rm c}}{2\pi \hbar} }[/math].

Note the dependency of the density of states with the magnetic field. The larger the magnetic field is, the more states are in each Landau level. As a consequence, there is more confinement in the system since less energy levels are occupied.

Rewriting the last expression as [math]\displaystyle{ n_B=\frac{\hbar \omega_{\rm c}}{2} \frac{m^*}{\pi \hbar^2} }[/math] it is clear that each Landau level contains as many states as in a 2DEG in a [math]\displaystyle{ \Delta\varepsilon=\hbar \omega_{\rm c} }[/math].

Given the fact that electrons are fermions, for each state available in the Landau levels it corresponds two electrons, one electron with each value for the spin [math]\displaystyle{ s=\pm\frac{1}{2} }[/math]. However, if a large magnetic field is applied, the energies split into two levels due to the magnetic moment associated with the alignment of the spin with the magnetic field. The difference in the energies is [math]\displaystyle{ \Delta E = \pm \frac{1}{2}g\mu_{\rm B}B }[/math] being [math]\displaystyle{ g }[/math] a factor which depends on the material ([math]\displaystyle{ g=2 }[/math] for free electrons) and [math]\displaystyle{ \mu_{\rm B} }[/math] the Bohr magneton. The sign [math]\displaystyle{ + }[/math] is taken when the spin is parallel to the field and [math]\displaystyle{ - }[/math] when it is antiparallel. This fact called spin splitting implies that the density of states for each level is reduced by a half. Note that [math]\displaystyle{ \Delta E }[/math] is proportional to the magnetic field so, the larger the magnetic field is, the more relevant is the split.

Density of states in a magnetic field, neglecting spin splitting. (a)The states in each range [math]\displaystyle{ \hbar \omega_{\rm c} }[/math] are squeezed into a [math]\displaystyle{ \delta }[/math]-function Landau level.(b) Landau levels have a non-zero width [math]\displaystyle{ \Gamma }[/math] in a more realistic picture and overlap if [math]\displaystyle{ \hbar \omega_{\rm c}\lt \Gamma }[/math]. (c) The levels become distinct when [math]\displaystyle{ \hbar \omega_{\rm c}\gt \Gamma }[/math].

In order to get the number of occupied Landau levels, one defines the so-called filling factor [math]\displaystyle{ \nu }[/math] as the ratio between the density of states in a 2DEG and the density of states in the Landau levels.

[math]\displaystyle{ \nu=\frac{n_{\rm 2D}}{n_B}=\frac{hn_{\rm 2D}}{eB} }[/math]

In general the filling factor [math]\displaystyle{ \nu }[/math] is not an integer. It happens to be an integer when there is an exact number of filled Landau levels. Instead, it becomes a non-integer when the top level is not fully occupied. Since [math]\displaystyle{ n_B\propto B }[/math], by increasing the magnetic field, the Landau levels move up in energy and the number of states in each level grow, so fewer electrons occupy the top level until it becomes empty. If the magnetic field keeps increasing, eventually, all electrons will be in the lowest Landau level ([math]\displaystyle{ \nu\lt 1 }[/math]) and this is called the magnetic quantum limit.

Occupation of Landau levels in a magnetic field neglecting the spin splitting, showing how the Fermi level moves to maintain a constant density of electrons. The fields are in the ratio [math]\displaystyle{ 2:3:4 }[/math] and give [math]\displaystyle{ \nu=4,\frac{8}{3} }[/math] and [math]\displaystyle{ 2 }[/math].

Longitudinal resistivity

It is possible to relate the filling factor to the resistivity and hence, to the conductivity of the system. When [math]\displaystyle{ \nu }[/math] is an integer, the Fermi energy lies in between Landau levels where there are no states available for carriers, so the conductivity becomes zero (it is considered that the magnetic field is big enough so that there is no overlap between Landau levels, otherwise there would be few electrons and the conductivity would be approximately [math]\displaystyle{ 0 }[/math]). Consequently, the resistivity becomes zero too (At very high magnetic fields it is proven that longitudinal conductivity and resistivity are proportional).[16]

Instead, when [math]\displaystyle{ \nu }[/math] is a half-integer, the Fermi energy is located at the peak of the density distribution of some Landau Level. This means that the conductivity will have a maximum .

This distribution of minimums and maximums corresponds to ¨quantum oscillations¨ called Shubnikov–de Haas oscillations which become more relevant as the magnetic field increases. Obviously, the height of the peaks are larger as the magnetic field increases since the density of states increases with the field, so there are more carriers which contribute to the resistivity. It is interesting to notice that if the magnetic field is very small, the longitudinal resistivity is a constant which means that the classical result is reached.

Longitudinal and transverse (Hall) resistivity, [math]\displaystyle{ \rho_{xx} }[/math] and [math]\displaystyle{ \rho_{xy} }[/math], of a two-dimensional electron gas as a function of magnetic field. Both vertical axes were divided by the quantum unit of conductance [math]\displaystyle{ e^2/h }[/math] (units are misleading). The filling factor [math]\displaystyle{ \nu }[/math] is displayed for the last 4 plateau.

Transverse resistivity

From the classical relation of the transverse resistivity [math]\displaystyle{ \rho_{xy}=\frac{B}{en_{\rm 2D}} }[/math] and substituting [math]\displaystyle{ n_{\rm 2D}=\nu \frac{eB}{h} }[/math] one finds out the quantization of the transverse resistivity and conductivity:

[math]\displaystyle{ \rho_{xy}=\frac{h}{\nu e^2}\Rightarrow \sigma=\nu \frac{e^2}{h} }[/math]

One concludes then, that the transverse resistivity is a multiple of the inverse of the so-called conductance quantum [math]\displaystyle{ e^2/h }[/math]. Nevertheless, in experiments a plateau is observed between Landau levels, which indicates that there are in fact charge carriers present. These carriers are localized in, for example, impurities of the material where they are trapped in orbits so they can not contribute to the conductivity. That is why the resistivity remains constant in between Landau levels. Again if the magnetic field decreases, one gets the classical result in which the resistivity is proportional to the magnetic field.

Photonic quantum Hall effect

The quantum Hall effect, in addition to being observed in two-dimensional electron systems, can be observed in photons. Photons do not possess inherent electric charge, but through the manipulation of discrete optical resonators and coupling phases or on-site phases, an artificial magnetic field can be created.[17][18][19][20][21] This process can be expressed through a metaphor of photons bouncing between multiple mirrors. By shooting the light across multiple mirrors, the photons are routed and gain additional phase proportional to their angular momentum. This creates an effect like they are in a magnetic field.


The integers that appear in the Hall effect are examples of topological quantum numbers. They are known in mathematics as the first Chern numbers and are closely related to Berry's phase. A striking model of much interest in this context is the Azbel–Harper–Hofstadter model whose quantum phase diagram is the Hofstadter butterfly shown in the figure. The vertical axis is the strength of the magnetic field and the horizontal axis is the chemical potential, which fixes the electron density. The colors represent the integer Hall conductances. Warm colors represent positive integers and cold colors negative integers. Note, however, that the density of states in these regions of quantized Hall conductance is zero; hence, they cannot produce the plateaus observed in the experiments. The phase diagram is fractal and has structure on all scales. In the figure there is an obvious self-similarity. In the presence of disorder, which is the source of the plateaus seen in the experiments, this diagram is very different and the fractal structure is mostly washed away.

Concerning physical mechanisms, impurities and/or particular states (e.g., edge currents) are important for both the 'integer' and 'fractional' effects. In addition, Coulomb interaction is also essential in the fractional quantum Hall effect. The observed strong similarity between integer and fractional quantum Hall effects is explained by the tendency of electrons to form bound states with an even number of magnetic flux quanta, called composite fermions.

The Bohr atom interpretation of the von Klitzing constant

The value of the von Klitzing constant may be obtained already on the level of a single atom within the Bohr model while looking at it as a single-electron Hall effect. While during the cyclotron motion on a circular orbit the centrifugal force is balanced by the Lorentz force responsible for the transverse induced voltage and the Hall effect, one may look at the Coulomb potential difference in the Bohr atom as the induced single atom Hall voltage and the periodic electron motion on a circle as a Hall current. Defining the single atom Hall current as a rate a single electron charge [math]\displaystyle{ e }[/math] is making Kepler revolutions with angular frequency [math]\displaystyle{ \omega }[/math]

[math]\displaystyle{ I = \frac{\omega e}{2\pi}, }[/math]

and the induced Hall voltage as a difference between the hydrogen nucleus Coulomb potential at the electron orbital point and at infinity:

[math]\displaystyle{ U=V_\text{C}(\infty) - V_\text{C}(r) = 0 - V_\text{C}(r) = \frac{e}{4\pi\epsilon_0 r} }[/math]

One obtains the quantization of the defined Bohr orbit Hall resistance in steps of the von Klitzing constant as

[math]\displaystyle{ R_\text{Bohr}(n) = \frac{U}{I} = n\frac{h}{e^2} }[/math]

which for the Bohr atom is linear but not inverse in the integer n.

Relativistic analogs

Relativistic examples of the integer quantum Hall effect and quantum spin Hall effect arise in the context of lattice gauge theory.[22][23]

See also


  1. Editorial (2020-07-29). "The quantum Hall effect continues to reveal its secrets to mathematicians and physicists" (in en). Nature 583 (7818): 659. doi:10.1038/d41586-020-02230-7. PMID 32728252. Bibcode2020Natur.583..659.. 
  2. F. D. M. Haldane (1988). "Model for a Quantum Hall Effect without Landau Levels: Condensed-Matter Realization of the 'Parity Anomaly'". Physical Review Letters 61 (18): 2015–2018. doi:10.1103/PhysRevLett.61.2015. PMID 10038961. Bibcode1988PhRvL..61.2015H. 
  3. Ezawa, Zyun F. (2013). Quantum Hall Effects: Recent Theoretical and Experimental Developments (3rd ed.). World Scientific. ISBN 978-981-4360-75-3. 
  4. 4.0 4.1 R. B. Laughlin (1981). "Quantized Hall conductivity in two dimensions". Phys. Rev. B 23 (10): 5632–5633. doi:10.1103/PhysRevB.23.5632. Bibcode1981PhRvB..23.5632L. 
  5. "2018 CODATA Value: conventional value of von Klitzing constant". The NIST Reference on Constants, Units, and Uncertainty. NIST. 20 May 2019. Retrieved 2019-05-20. 
  6. "26th CGPM Resolutions". 
  7. "2018 CODATA Value: von Klitzing constant". The NIST Reference on Constants, Units, and Uncertainty. NIST. 20 May 2019. Retrieved 2019-05-20. 
  8. "1960 - Metal Oxide Semiconductor (MOS) Transistor Demonstrated". The Silicon Engine (Computer History Museum). 
  9. 9.0 9.1 9.2 Lindley, David (15 May 2015). "Focus: Landmarks—Accidental Discovery Leads to Calibration Standard". Physics 8. doi:10.1103/physics.8.46. 
  10. Tsuneya Ando; Yukio Matsumoto; Yasutada Uemura (1975). "Theory of Hall effect in a two-dimensional electron system". J. Phys. Soc. Jpn. 39 (2): 279–288. doi:10.1143/JPSJ.39.279. Bibcode1975JPSJ...39..279A. 
  11. Jun-ichi Wakabayashi; Shinji Kawaji (1978). "Hall effect in silicon MOS inversion layers under strong magnetic fields". J. Phys. Soc. Jpn. 44 (6): 1839. doi:10.1143/JPSJ.44.1839. Bibcode1978JPSJ...44.1839W. 
  12. K. v. Klitzing; G. Dorda; M. Pepper (1980). "New method for high-accuracy determination of the fine-structure constant based on quantized Hall resistance". Phys. Rev. Lett. 45 (6): 494–497. doi:10.1103/PhysRevLett.45.494. Bibcode1980PhRvL..45..494K. 
  13. D. J. Thouless (1983). "Quantization of particle transport". Phys. Rev. B 27 (10): 6083–6087. doi:10.1103/PhysRevB.27.6083. Bibcode1983PhRvB..27.6083T. 
  14. K. S. Novoselov; Z. Jiang; Y. Zhang; S. V. Morozov; H. L. Stormer; U. Zeitler; J. C. Maan; G. S. Boebinger et al. (2007). "Room-temperature quantum Hall effect in graphene". Science 315 (5817): 1379. doi:10.1126/science.1137201. PMID 17303717. Bibcode2007Sci...315.1379N. 
  15. Tsukazaki, A.; Ohtomo, A.; Kita, T.; Ohno, Y.; Ohno, H.; Kawasaki, M. (2007). "Quantum Hall effect in polar oxide heterostructures". Science 315 (5817): 1388–91. doi:10.1126/science.1137430. PMID 17255474. Bibcode2007Sci...315.1388T. 
  16. Davies J.H.. The physics of low-dimension. 6.4 Uniform magnetic Field; 6.5 Magnetic Field in a Narrow Channel, 6.6 The Quantum Hall Effect. ISBN 9780511819070. 
  17. Raghu, S.; Haldane, F. D. M. (2008-09-23). "Analogs of quantum-Hall-effect edge states in photonic crystals" (in en). Physical Review A 78 (3): 033834. doi:10.1103/PhysRevA.78.033834. ISSN 1050-2947. 
  18. Fang, Kejie; Yu, Zongfu; Fan, Shanhui (November 2012). "Realizing effective magnetic field for photons by controlling the phase of dynamic modulation" (in en). Nature Photonics 6 (11): 782–787. doi:10.1038/nphoton.2012.236. ISSN 1749-4885. 
  19. Schine, Nathan; Ryou, Albert; Gromov, Andrey; Sommer, Ariel; Simon, Jonathan (June 2016). "Synthetic Landau levels for photons" (in en). Nature 534 (7609): 671–675. doi:10.1038/nature17943. ISSN 0028-0836. PMID 27281214. Bibcode2016Natur.534..671S. 
  20. Minkov, Momchil; Savona, Vincenzo (2016-02-20). "Haldane quantum Hall effect for light in a dynamically modulated array of resonators" (in en). Optica 3 (2): 200. doi:10.1364/OPTICA.3.000200. ISSN 2334-2536. 
  21. Dutt, Avik; Lin, Qian; Yuan, Luqi; Minkov, Momchil; Xiao, Meng; Fan, Shanhui (2020-01-03). "A single photonic cavity with two independent physical synthetic dimensions" (in en). Science 367 (6473): 59–64. doi:10.1126/science.aaz3071. ISSN 0036-8075. 
  22. D. B. Kaplan (1992). "A Method for simulating chiral fermions on the lattice". Physics Letters B288 (3–4): 342–347. doi:10.1016/0370-2693(92)91112-M. Bibcode1992PhLB..288..342K. 
  23. M. F. L. Golterman; K. Jansen; D. B. Kaplan (1993). "Chern-Simons currents and chiral fermions on the lattice". Physics Letters B301 (2–3): 219–223. doi:10.1016/0370-2693(93)90692-B. Bibcode1993PhLB..301..219G. 

Further reading