Neutral excitation and bulk gap of fractional quantum Hall liquids in disk geometry∗

Wu-Qing Yang(杨武庆),Qi Li(李骐),Lin-Peng Yang(杨林鹏),and Zi-Xiang Hu(胡自翔),†

1Department of Physics,Chongqing University,Chongqing 401331,China

2Department of Physics,South University of Science and Technology of China,Shenzhen 518055,China

3School of Physics and Technology,Wuhan University,Wuhan 430072,China

Keywords:fractional quantum Hall,bulk states,edge states,magneto-roton,dipolar-dipolar interaction

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,suchasDMRG[2-5]orMPS.[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.[7-9]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,[11-17]quasiparticle interference,[18-28]and edge reconstruction.[29-33]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 inthe bulk.[35-37]There aremorepara meters 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 magnetoroton 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 V1interaction 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,Bz= εjk∂jAk,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 r=R+.The cyclotron coordinates are de fined by the canonical momentum:

It is easy to know

Fig.1.The spectrum of V1Hamiltonian 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 Mtotal＞135 are the edge excitation states.

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 H=(a†a+1/2)¯hωcand 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 Neelectrons filling Norborbitals.Because of the rotational symmetry along z axis,the total angular momentum Mtotal=∑imiis conserved and mostly used in the exact diagonalization.Figure1shows the energy spectrum with V1Hamiltonian 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 statehas total angular momentum Mtotal=ML=3Ne(Ne-1)/2=135 for 10 electrons.The low-lying excited states that only exist in the subspace Mtotal＞MLare the chiral edge states.Here we should notice that the edge excitations have zero energies in an in finite plane for the model Hamiltonian.The finite number of orbitals results in finite energies for partially edge states while Mtotal＞ML+Norb-3Ne+2.

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 con figuration “100100···1001001”which has edge counting “1,1,2,3,5,...”related to the chiral Luttinger liquid theory of the edge.However,some compressible charge density wave state,such as the largest antisqueezed pattern “111···000···111”,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 de fined 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 de fined 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.

Fig.2.The spectrum of the COM operator and energy spectrum for 3 electrons in 20 orbitals withV1interaction.

In the energy spectrum in Mtotalsubspace 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 in finite potential barrier.The Laughlin state is the densest zero energy state for Neelectrons in 3Ne-2 orbitals for V1interaction,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 de fined 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=B†B:

For an in finite plane with translational invariance,the[,H]=0 and thus MCOMis a good quantum number.We expect that the energy spectrum in different MCOMsubspaces 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 theperator in a system with 20 guiding center orbitals with the quantum number Mtotal.The eigenvalues of theare well quantized to integers while Mtotal＜Norb.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 truncatedsubspace.The results for 3 electrons in 20 orbitals are shown in Fig.2(b).It shows that the spectrum is highly degenerated for differentsubspaces.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 forV1interaction as the way we would expect.Second,the first excited state de fines the bulk gap.As shown in Fig.3,we compare the bulk gap for theV1interaction 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 V1interaction.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 Mtotalfor 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 M0-Mtotalwhere M0is the total angular momentum for ground state,i.e.,M0=3Ne(Ne-1)/2,the magneto-roton comes out.Thus the magneto-roton excitation actually hidden in the Mtotal＜M0subspace in Fig.1.As a comparison,we plot the energy spectrum for the same V1Hamiltonian 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.3.The comparison of the band gap for the V1interaction(a)and ground state energy for Coulomb interaction(b)between disk and sphere geometries.

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 of52Cr atom[42]or the degenerate quantum gas of40K87Rb,[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 x-z plane.The polarized dipole interaction(in unit of d2/(4πε0l3),where d is the dipole moment of the neutral atom and ε0is the vacuum permittivity)is

In this section,we want to provide complementary numericalresults about this phase transition in the disk geometry.In this case,we have one additional parameter,the con fining potential strength γ which characterizes the relative strength of con fining 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 V1+λV1,2.Therefore,as the V1model 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.

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 con fining potential γm in the m′-th orbital.

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 α ≃ 53°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 pro files 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 pro files 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 andtheir pro files 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 V1model as shownin 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 if lling 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 in finite 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.