Neutral excitation and bulk gap of fractional quantum Hall liquids in disk geometry
Yang Wu-Qing1, Li Qi2, 3, Yang Lin-Peng1, Hu Zi-Xiang1, †
Department of Physics, Chongqing University, Chongqing 401331, China
Department of Physics, South University of Science and Technology of China, Shenzhen 518055, China
School of Physics and Technology, Wuhan University, Wuhan 430072, China

 

† Corresponding author. E-mail: zxhu@cqu.edu.cn

Abstract

For the numerical simulation of the fractional quantum Hall (FQH) effects on a finite disk, the rotational symmetry is the only symmetry that is used in diagonalizing the Hamiltonian. In this work, we propose a method of using the weak translational symmetry for the center of mass of the many-body system. With this approach, the bulk properties, such as the energy gap and the magneto-roton excitation are consistent with those in the closed manifolds like the sphere and torus. As an application, we consider the FQH phase and its phase transition in the fast rotated dipolar fermions. We thus demonstrate the disk geometry having versatility in analyzing the bulk properties beside the usual edge physics.

1. Introduction

The fractional quantum Hall effect (FQHE), a topological quantum state of matter which was experimentally realized in two-dimensional (2D) electron gas placed in a low temperature and strong magnetic field environment,[1] has attracted strong interest of theoretical and experimental physicists due to its embedded nature of the electronic topology and strong correlation. Since the kinetic energy of electrons has been frozen by a strong magnetic field, the FQHE system is typically strongly correlated and cannot be treated by the perturbation approach. The most powerful numerical tool for studying the FQHE is exact diagonalizing a microscopic Hamiltonian for small number of electrons or other advanced numerical methods, such as DMRG[25] or MPS.[6] Theoretically, the numerical calculation can be applied in different geometries for different purposes. For the case of compact geometries without edge, such as putting electrons on the surface of a sphere or a torus, one always considers the bulk topological properties of the FQH states, i.e., the ground state topological degeneracy and magneto-roton excitation.[79] The open boundary systems, such as the cylinder and disk geometries, are always aimed for the exploration of the edge physics,[10] such as the edge tunneling,[1117] quasiparticle interference,[1828] and edge reconstruction.[2933] Thanks for the bulk–edge correspondence, the topological properties of the FQH state can be unveiled from both the bulk and edge perspectives.

For the compact geometries, the proposal of the entanglement spectrum by Li and Haldane[34] supplies a way of detecting the edge physics of the FQH droplets by artificially bi-partite the system. The entanglement spectrum is actually the eigenvalue spectrum of the reduced density matrix for subsystem after truncating the rest part of the system. It reflects the bulk topology via the counting of the so-called conformal edge states. However, the entanglement spectrum cannot give the quantitative properties of the edge excitation of the FQH droplet, such as the edge velocities or edge reconstruction. In a Hall bar sample, the existence of an edge is unavoidable except the Corbino geometry which provides direct access to the longitudinal conductivity in the bulk.[3537] There are more parameters such as the strength of the background confinement and edge potential to tune the system. With these knobs, the characteristics of the FQH edge are frequently explored. The edge excitation is gapless, which overwhelms the bulk excitation in the low-energy sector. The bulk excitation, such as the magneto-roton, is rarely discussed in disk geometry and therefore, the topological phase transition accompanying gap closing in the bulk has been incompletely discussed.

In this paper, we give a way for digging out the magneto-roton excitation of the FQH liquid in disk geometry and exactly determining the bulk energy gap by using the degrees of freedom in the center of mass system. For the model Hamiltonian with V1 interaction which gives the unambiguous Laughlin state and its low-lying excitations, we find the spectrum of magneto-roton excitation is exactly matched with that from the sphere geometry. The energy gap is found to be less sensitive to the finite size effect in our approach. As an example, we consider the dipolar-interacted neutral atoms in a fast rotated trap.[38] A phase transition from the FQH region to molecular phase as tilting the dipolar angle is characterized by the gap closing. In the anisotropic FQH state before gap closing, we observed multiple branches in the magneto-roton excitation. The rest of this paper is organized as following: In Section 2, we review the model of the FQH system and the edge spectrum in orbital space on a disk. The center of mass diagonalization is introduced in Section 3. We make a comparison of the energy gap with the spherical geometry and the magneto-roton excitation with torus geometry. As an example of application, we consider the dipole–dipole interaction between the neutral atoms in a fast rotated trap in Section 4, in which the bulk gap closing manifests the phase transition and the multi-branch of the magneto-roton excitation indicates the anisotropy of the FQH state. Conclusions and discussions are arranged in the last section.

2. Cyclotron motion in a magnetic field

For an electron in a magnetic field along z direction, , the single particle Hamiltonian is

The coordinates of an electron in a magnetic field can be decomposed into the guiding center coordinates and cyclotron coordinates, namely . The cyclotron coordinates are defined by the canonical momentum:
It is easy to know

With these two types of coordinates, we can construct corresponding ladder operators

they obey the following commutation relations:
The Hamiltonian can be written in a diagonal form of one-dimensional (1D) harmonic oscillator and the Hilbert space can be expanded by the above two ladder operators
The index N labels the Landau level (LL) and m labels the degenerated orbitals in each LL. After projecting into the lowest Landau level, only the orbital angular momentum quantum number m remains and the system can be treated as a quasi-1D system for Ne electrons filling Norb orbitals. Because of the rotational symmetry along z axis, the total angular momentum is conserved and mostly used in the exact diagonalization. Figure 1 shows the energy spectrum with V1 Hamiltonian in the language of the Haldane’s pseudopotential[39] for 10 electrons in 30 orbitals. This is the model Hamiltonian for the celebrated Laughlin state[40] at filling factor 1/3. The Laughlin state has total angular momentum for 10 electrons. The low-lying excited states that only exist in the subspace are the chiral edge states. Here we should notice that the edge excitations have zero energies in an infinite plane for the model Hamiltonian. The finite number of orbitals results in finite energies for partially edge states while .

Fig. 1. The spectrum of V1 Hamiltonian for 10 electrons in 30 orbitals in disk geometry. The unique zero energy state at Mtotal = 135 is the Laughlin state and the low-lying states in are the edge excitation states.
3. Center of mass diagonalization

The information from the edge spectrum of the FQH state, especially for the realistic Coulomb interaction, sometimes does not completely specify the FQH state itself. For example, we know the incompressible Laughlin phase at 1/3 filling is accompanied by the root configuration “ ” which has edge counting “ ” related to the chiral Luttinger liquid theory of the edge. However, some compressible charge density wave state, such as the largest anti-squeezed pattern “ ”, could also be the global ground state and has exactly the same counting numbers because of the integer quantum Hall-like edge near the edge. Therefore, the edge spectrum is just one aspect of the topological order for the FQH liquids. On the other hand, the topological quantum ground state is always defined by a bulk gap which is the origination of the incompressibility for the FQH liquids. The ground state phase transition is displayed by the gap closing. Besides the gapless edge excitations on the boundary, the FQH liquids are embedded some bulk excitations. The most studied ones are the quasihole and quasielectron excitations which can be realized by a local potential[41] that mimics the STM or AFM tip in experiments. Relatively little progress has been made in understanding the neutral excitations since the seminal work of Girvin, MacDonald, and Platzman[7] introduced the single-mode approximation (SMA) to describe the lowest excitation in terms of a neutral density wave, or the “magnetoroton”. For numerical simulation, the neutral excitation can be observed as an in-gap mode in the spectrum of a translational-invariant torus geometry in both directions or rotational-invariant sphere in the total angular momentum L = 0 subspace. The bulk gap of the FQH liquid is defined as the gap between the roton minimum and the ground state. Thus the precise definition of the bulk gap, or finding out the roton minimum, is important to locate the phase transition of the ground state.

In the energy spectrum in Mtotal subspace on a disk, due to the gapless state of the FQH edge excitation, we cannot exactly locate the bulk gap closing except the changing of the global ground state due to edge reconstruction[31] in the lower energy sector. For electrons on a finite disk, the translational symmetry is broken near the edge and the open boundary can be regarded as an infinite potential barrier. The Laughlin state is the densest zero energy state for Ne electrons in orbitals for V1 interaction, thus all the electrons prefer to form a droplet at the center of the disk. If the ground state is a homogeneous liquid state, we would expect that the center of mass (COM) locates at the center of disk, or the angular momentum for COM should be zero which is similar to the total angular momentum L = 0 for the spherical geometry. The COM ladder operator for a many-body state is defined as the summation of all the operators in each guiding center orbitals.

Therefore, the COM for a given state can be calculated as an expect value of its number operator :

For an infinite plane with translational invariance, the and thus MCOM is a good quantum number. We expect that the energy spectrum in different MCOM subspaces should be degenerate. The presence of the edge for a finite disk will lift this degeneracy when the COM operator has contributions from the edge. As shown in Fig. 2(a), we diagonalize the operator in a system with 20 guiding center orbitals with the quantum number Mtotal. The eigenvalues of the are well quantized to integers while . If we make a truncation to keep all the well-quantized states in the spectrum of , the electron Hamiltonian can be approximately diagonalized in each truncated subspace. The results for 3 electrons in 20 orbitals are shown in Fig. 2(b). It shows that the spectrum is highly degenerated for different subspaces. Some high energy excited states gradually vanish as increasing MCOM. Therefore, from the experience of the spherical geometry, we speculate that the spectrum in MCOM=0 subspace contains all the information in the bulk. First of all, the ground state energy is exactly zero for V1 interaction as the way we would expect. Second, the first excited state defines the bulk gap. As shown in Fig. 3, we compare the bulk gap for the V1 interaction and the ground state energy per particle for Coulomb interaction with that on sphere for several system sizes. It is found that the extrapolated values for the thermodynamic limit are consistent for different geometries. However, it is obvious that both the energy gap (V1) and single particle energy (Coulomb) in disk geometry are always smaller than that on sphere. Moreover, the gap has much fewer finite size effects in disk geometry. It reaches the saturation value even at a very small size for 5 electrons with V1 interaction. At last, except the bulk gap, we can also extract the neutral magneto-roton excitation from the energy spectrum in MCOM=0 subspace. Because the rotational invariance is always there, we calculate the expected value of Mtotal for each of states in MCOM=0 subspace and replot the spectrum as a function of Mtotal. The results are shown in Fig. 4(a). While we plot the spectrum by where M0 is the total angular momentum for ground state, i.e., , the magneto-roton comes out. Thus the magneto-roton excitation actually hidden in the subspace in Fig. 1. As a comparison, we plot the energy spectrum for the same V1 Hamiltonian in sphere (with the total angular momentum L) and torus geometry (with the total translational momentum k) as shown in Figs. 4(b) and 4(c) respectively. The two closed manifolds clearly show the magneto-roton excitation branch in the spectrum. As we can see, the curves for the neutral excitations are almost the same although the system size we can do in the disk geometry is smaller. Because the full diagonalization of the COM matrix to construct the MCOM=0 basis limits the system size.

Fig. 2. The spectrum of the COM operator and energy spectrum for 3 electrons in 20 orbitals with V1 interaction.
Fig. 3. The comparison of the band gap for the V1 interaction (a) and ground state energy for Coulomb interaction (b) between disk and sphere geometries.
Fig. 4. The magneto-roton spectrum on disk (a), sphere (b), and torus (c) for V1 Hamiltonian. The energy gap is labeled by the difference between ground state and the roton minimum.

Therefore, as we can see, all the bulk properties of the FQH liquid could be extracted from a finite disk. All we just need is to diagonalize the electron Hamiltonian in the MCOM=0 subspace. The ground state, the gap, and the neutral excitation in the bulk can be obtained. Because of no curvature, the gap in disk geometry has much fewer finite size effects.

4. Phase transition of the dipolar atoms in FQH regime

In the Bose–Einstein condensation of 52Cr atom[42] or the degenerate quantum gas of 40K87Rb,[43] the interaction can be described as dipole–dipole interaction with the s-wave collisional interaction vanishing for spin-polarized fermions. We assume all the dipoles are polarized in the same direction and without loss of generality, say xz plane. The polarized dipole interaction (in unit of , where d is the dipole moment of the neutral atom and is the vacuum permittivity) is

with being the direction of the polarized dipole (parametrized by the angle α) as sketched in Fig. 5. When the dipole fermions are placed in a fast rotated trap and the rotating frequency is approaching the trap frequency along z axis, the Coriolis force experienced by the atoms results in an effective magnetic field. Then the FQH states are expected in the fast rotating limit.[38,4447] By using our recently developed generalized pseudopotential description of the anisotropic FQH state,[48,49] we studied[50] the anisotropic FQH state and its quantum phase transition as tilting the dipolar angle α. In torus geometry, it was shown that the gap of the Laughlin-like state at 1/3 filling was closed at when the thickness was set to q = 0.01. Previous studies in disk geometry[51] could not gain this critical value because of the edge effects.

Fig. 5. The model for the dipolar fermions. α is the angle between the direction of the dipole moment and z axis. The system is placed in a confining potential in the -th orbital.

In this section, we want to provide complementary numerical results about this phase transition in the disk geometry. In this case, we have one additional parameter, the confining potential strength γ which characterizes the relative strength of confining potential with respect to interaction. In our calculation, we always choose the value of γ to maximize the energy gap. The pseudopotential analysis[50] tells us that the dipole–dipole interaction in the lowest Landau level can be described by a model Hamiltonian . Therefore, as the V1 model in Section 3, we expect the finite size effect is small in the dipole–dipole interaction and the results for small system should be good enough to predict the physics in the thermodynamical limit. On the other hand, the rotational invariance is broken while the dipole moment has a component in the plane, we thus need to include the orbitals as many as possible in order to gain convergence of the energy.

In Fig. 6(a), we plot the energy spectrum in the MCOM=0 subspace for 5 particles in 15 orbitals with confinement γ = 0.005. It shows clearly that the bulk gap between the ground state and the first excited state gradually decreases as increasing the angle α. The gap closes at while q = 0.01 which is exactly the same as that from the calculation in torus geometry.[50] To trace the ground state, we calculate the largest wave function overlap between the Laughlin state and the eigenstates in each subspace. Predictably, the ground state before gap closing has the largest overlap with the model wave function. Here we just consider the isotropic Laughlin state as a reference, it is shown in Fig. 6(b) that the overlap is almost 100% even though the tilted angle α reaches to 30° which demonstrates that the ground state is still isotropic. If we continue to trace the highest overlap state in the spectrum, as outlined in blue, it looks like the energy of the Laughlin-like state flows up and becomes the highest energy state finally. When α = 90°, we find the highest state has overlap ∼75% with Laughlin state. In Figs. 6(c) and 6(d), we plot the particle density and their profiles along x/y directions for the ground state at α = 0° and highest state at α = 90°. When α = 0°, the droplet is isotropic as expected. For α = 90°, the shape of the density is elliptical and has different profiles along x and y directions. We should note that the reference state in the overlap calculation is the isotropic Laughlin state which does not contain the geometric metric.[52] This anisotropic Laughlin-like state should have larger overlap with the anisotropic model wave function.[53]

Fig. 6. (a) The energy spectrum for 5 dipolar fermions in 15 ortitals with the confinement potential γ=0.005 as varying α. The states labelled by dark blue points have the maximum overlap with Laughlin state. (b) The largest overlap along the dark blue points. The particle density and their profiles along x/y directions for the ground state at α=0° (c) and the highest energy state at α=90° (d).

In spite of the fact that COM diagonalization can provide the bulk gap and the outline of the neutral magneto-roton excitation, since the angular momentum is discreted, we do not find a way to get the continuous spectrum for the roton excitation as that in torus (The reason that one can get the continuous spectrum on torus is that the angle between two primitive lattice vector can smoothly be tuned). Here, as an attachment, we plot the magneto-roton spectrum for the dipolar fermions at 1/3 filling in the LLL in torus geometry. As shown in Fig. 7, when α = 0°, the spectrum looks similar to the V1 model as shown in Fig. 4(c) except that the long-range interaction makes the points scattered for large k part. When α = 30° at which the ground state still looks isotropic since it has almost 100% overlap with the isotropic model wave function as shown in Fig. 6(b). However, the neutral excitation in the bulk as shown in Fig. 7(b) has some multi-branch structures. This result is similar to the case while considering the effects of anisotropic effective mass on the FQH state.[54] It is reasonable since we know both the anisotropic effective mass and the anisotropic interaction could be equally treated as an intrinsic metric in the geometric description of the FQH states.[52]

Fig. 7. The energy spectrum on torus for 10 dipolar fermions at 1/3 filling in LLL at α=0° (a) and α=30° (b).
5. Summaries and conclusions

In conclusion, we propose a numerical method to extract the bulk properties of the FQH liquids, such as the energy gap and the neutral magneto-roton excitation in the disk geometry which is mostly used to study the edge physics in the guiding center orbitals. The main point in our approach is constructing the electron Hamiltonian in the subspace of the COM angular momentum. Although the COM angular momentum is conserved in an infinite plane, we find the diagonalization of the interaction Hamiltonian in a truncated MCOM=0 subspace provides all the information in the bulk. By comparing the results in closed geometries, such as sphere and torus which naturally provide the bulk properties, we find the extrapolated energy gap in thermodynamic limit is consistent both for the model Hamiltonian and the Coulomb interaction. More importantly, because of no curvature in the plane, the finite size effect is smaller than that on a sphere. The magneto-roton excitation appears while all the eigen-energies in the MCOM=0 subspace are plotted as a function of the Mtotal. It looks very similar to the spectrum on the sphere while considering the total angular momentum L. As an application, we consider the fast rotated dipolar fermions in quantum Hall regime. The critical point of the gap closing for the phase transition induced by tilting the dipole angle is consistent to the previous study on torus. We also find that the Laughlin state in the COM spectrum flows to be the highest energy state when all the dipoles are parallel to plane although it comes anisotropic. The multi-branch structure of the magneto-roton excitation in the anisotropic FQH state is also consistent to the previous study of the anisotropic effective mass. The weakness of this method is the system size has an upper limit due to the full diagonalization of the COM matrix. To avoid confusion, here we should note that the geometries we discussed above are the theoretical models in the numerical calculations which has nothing to do with the sample geometries in experiments.

Reference
[1] Tsui D C Störmer H L Gossard A C 1982 Phys. Rev. Lett. 48 1559
[2] Shibata N Yoshioka D 2001 Phys. Rev. Lett. 86 5755
[3] Feiguin A E Rezayi E Nayak C Das Sarma S 2008 Phys. Rev. Lett. 100 166803
[4] Zhao J Z Sheng D N Haldane F D M 2011 Phys. Rev. 83 195135
[5] Hu Z X Papic Z Johri S Bhatt R N Schmitteckert P 2012 Phys. Lett. 376 2157
[6] Zaletel M P Mong R S K Pollmann F 2013 Phys. Rev. Lett. 110 236801
[7] Girvin S M MacDonald A H Platzman P M 1985 Phys. Rev. Lett. 54 581
[8] Girvin S M MacDonald A H Platzman P M 1986 Phys. Rev. 33 2481
[9] Yang B Hu Z X Papic Z Haldane F D M 2012 Phys. Rev. Lett. 108 256807
[10] Wen X G 1992 Int. J. Mod. Phys. 6 1711
[11] Kan C L Fisher M P A 1992 Phys. Rev. 46 15233
[12] Moon K Yi H Kane C L Girvin S M Fisher M P A 1993 Phys. Rev. Lett. 71 4381
[13] Kane C L Fisher M P A 1994 Phys. Rev. Lett. 72 724
[14] de Chamon C Wen X G 1993 Phys. Rev. Lett. 70 2605
[15] Lin X Dillard C Kastner M A Pfeiffer L N West K W 2012 Phys. Rev. 85 165321
[16] Li Q Jiang N Wan X Hu Z X 2018 J. Phys.: Condens. Matter 30 255601
[17] Hu Z X Lee K H Rezayi E H Wan X Yang K 2011 New J. Phys. 13 035020
[18] Fendley P Fisher M P A Nayak C 2007 Phys. Rev. 75 045317
[19] Fradkin E Nayak C Tsvelik A Wilczek F 1998 Nucl. Phys. 516 704
[20] Das Sarma S Freedman M Nayak C 2005 Phys. Rev. Lett. 94 166802
[21] Stern A Halperin B I 2006 Phys. Rev. Lett. 96 016802
[22] Rosenow B Halperin B I Simon S H Stern A 2008 Phys. Rev. Lett. 100 226803
[23] Bishara B Nayak C 2008 Phys. Rev. 77 165302
[24] Bonderson P Shtengel K Slingerland J K 2006 Phys. Rev. Lett. 97 016401
[25] Willett R L Pfeiffer L N West K W 2009 Proc. Natl. Acad Sci. 106 8853
[26] Willett R L Pfeiffer L N West K W 2010 Phys. Rev. 82 205301
[27] Willett R L, Pfeiffer L N and West K W 2013 arXiv: 1301.2594
[28] Lin X Du R R Xie X C 2014 National Science Review 1 564
[29] MacDonald A H 1990 Phys. Rev. Lett. 64 220
[30] de Chamon C Wen X G 1994 Phys. Rev. 49 8227
[31] Wan X Rezayi E H Yang K 2003 Phys. Rev. 68 125307
[32] Wang J H Meir Y Gefen Y 2013 Phys. Rev. Lett. 111 246803
[33] Sabo R Gurman I et al 2017 Nat. Phys. 13 491
[34] Li H Haldane F D M 2008 Phys. Rev. Lett. 101 010504
[35] Laughlin R B 1981 Phys. Rev. 23 5632
[36] Dolgopolov V T Shashkin A A Zhitenev N B Dorozhkin S I von Klitzing K 1992 Phys. Rev. 46 12560
[37] Zhu Z Y Wang P J Fu H L Pfeiffer L N West K W Du R R Lin X 2018 Physica 95 1
[38] Cooper N R 2008 Adv. Phys. 57 539
[39] Haldane F D M 1983 Phys. Rev. Lett. 51 605
[40] Laughlin R B 1983 Phys. Rev. Lett. 50 1395
[41] Hu Z X Wan X Schmitteckert P 2008 Phys. Rev. 77 075331
[42] Griesmaier A Werner J Hensler S Stuhler J Pfau T 2005 Phys. Rev. Lett. 94 160401
[43] Myatt C J Burt E A Ghrist R W Cornell E A Wieman E A 1997 Phys. Rev. Lett. 78 586
[44] Fetter A L 2009 Rev. Mod. Phys. 81 647
[45] Baranov M A Osterloh K Lewenstein M 2005 Phys. Rev. Lett. 94 070404
[46] Osterloh K Barberán N Lewenstein M 2007 Phys. Rev. Lett. 99 160403
[47] Baranov M A Fehrmann H Lewenstein M 2008 Phys. Rev. Lett. 100 200402
[48] Yang B Hu Z X Lee C H Papic Z 2017 Phys. Rev. Lett. 118 146403
[49] Yang L P Li Q Hu Z X 2018 Chin. Phys. 27 087306
[50] Hu Z X Li Q Yang L P Yang W Q Jiang N Qiu R Z Yang B 2018 Phys. Rev. 97 035140
[51] Qiu R Z Kou S P Hu Z X Wan X Yi S 2011 Phys. Rev. 83 063633
[52] Haldane F D M 2011 Phys. Rev. Lett. 107 116801
[53] Qiu R Z Haldane F D M Wan X Yang K Yi S 2012 Phys. Rev. 85 115308
[54] Yang B Papic Z Rezayi E H Bhatt R N Haldane F D M 2012 Phys. Rev. 85 165318