Eigenstates of the planar magnetic Laplacian with a homogeneous magnetic field form degenerate energy bands, the Landau levels. We discuss the unitary correspondence between states in higher Landau levels and those in the lowest Landau level, where wave functions are holomorphic. We apply this correspondence to many-body systems; in particular, we represent effective Hamiltonians and particle densities in higher Landau levels by using corresponding quantities in the lowest Landau level.
I. INTRODUCTION
The state space of a charged particle moving in a homogeneous magnetic field in a plane orthogonal to the field decomposes into Landau levels, differing in energy by integral multiples of the magnetic field strength. When the position coordinates are expressed as complex numbers in the symmetric gauge, the states in the lowest Landau level form a Bargmann space of holomorphic functions, while wave functions in higher Landau levels involve also powers of the complex conjugate position variables in the standard representation.
As noted by many people for a long time, and emphasized, in particular, in Refs. 1–3, a holomorphic representation of states is not limited to the lowest Landau level, where it has proved to be important for deriving some basic properties, e.g., Refs. 4–8. In fact, there is a natural unitary correspondence between states in different Landau levels, in particular, between higher levels and the lowest one.
In this expository paper, we discuss several ways to arrive at the holomorphic representations and derive explicit formulas for particle densities and effective Hamiltonians in higher Landau levels, expressed in terms of corresponding quantities in the lowest Landau level. The methods have appeared in various disguises in the literature before, but our aim is to present them in a coherent fashion that, we hope, will be found to be useful for students and researchers in quantum Hall physics.
A physically appealing starting point is the decomposition of the position variables into guiding center variables and variables associated with the cyclotron motion of the particle around the guiding centers. While the components of the position operator commute, the other two sets of variables are non-commutative and satisfy canonical commutation relations. They can be represented in terms of creation and annihilation operators for two distinct and mutually commuting harmonic oscillators. One way of arriving at a holomorphic representation of states is an expansion in terms of coherent states for the harmonic oscillator of the guiding center variables.1,3 (These are the same as the “vortex eigenstates” in Refs. 9–11.) The transformation between the position coordinates and the coherent state variables can also be expressed in terms of an integral operator with a kernel that is a modification of the reproducing kernel of a Bargmann space.9
A formally simpler and more direct approach is to use the creation and annihilation operators of the cyclotron oscillator to define unitary mappings between different Landau levels.12 This gives explicit formulas for particle densities of many-body states in one Landau level in terms of polynomials in the Laplacian applied to the corresponding densities in the lowest Landau level. The same formulas can alternatively be obtained by a Fourier transformation, exploiting the factorization of the exponential factor in the guiding center and cyclotron variables, respectively.
The main application of these considerations is in quantum Hall physics.13–17 In this context, an electron gas is confined to two spatial dimensions and submitted to a magnetic field large enough to set the main energy scale. The quantization of the kinetic energy levels then becomes the salient feature. In the full plane, each level is infinitely degenerate, but for a finite area, the degeneracy is proportional to the area times the field strength. For extremely large values of the latter, the lowest Landau level is degenerate enough to accommodate all electrons without violating the Pauli principle. For smaller values of the field, several Landau levels can be completely filled with electrons and become inert in the first approximation. The physics then boils down to the motion of the electrons in the last, partially filled, Landau level.
In both cases, only one Landau level has to be taken into account, and an effective model of widespread use in the literature is given in terms of a Hamiltonian acting on holomorphic functions. We review this first, before describing in more detail the unitary mappings between Landau levels. The remarkable fact is that the dependence of the effective Hamiltonian on the Landau level it corresponds to is quite simple and transparent. An intuitive explanation (albeit not the most direct one from a computational point of view) is that the good variables to use are not the position variables, but rather those of the guiding centers. The Landau level index, which fixes the energy of the cyclotron motion, is encoded in a form factor in Fourier space that modifies external and interaction potentials via a differential operator. In particular, the unitary mappings between Landau levels map multiplication by potentials to operators of the same kind.
One salient feature of the effective operators acting on holomorphic functions is that they naturally suggest variational Ansätze for their ground states, which become exact for certain truncated models. The Laughlin state18,19 is the most emblematic of those, and much of our understanding of the fractional quantum Hall effect rests on its remarkable properties. In Sec. VI, we apply our formulas to Laughlin states in an arbitrary Landau level, computing their density profiles and extending rigidity results from Refs. 4–6 and 8.
II. PROJECTED HAMILTONIANS AND DENSITIES IN QUANTUM HALL PHYSICS
Let us start from the many-body Hamiltonian (in symmetric gauge) for interacting 2D electrons in a constant perpendicular external magnetic field B and a one-body potential V,
Here, w is the radial repulsive pair interaction potential, modeling 3D Coulomb interactions20 in quantum Hall (QH) physics, but more general choices are also of interest. The one-body potential V incorporates trapping in a finite size sample, plus the electrostatic potential generated by impurities. Mathematical conditions on the potentials V and w will be stated below. For convenience, we have subtracted NB from the kinetic part of the energy so that its lowest value is 0. In the sequel vectors, will very often be identified with complex numbers z = x + iy ∈ C.
As appropriate for electrons, we consider the action of H on the fermionic antisymmetric space,
For bosons, one considers the symmetric tensor product instead; this is relevant for rotating cold atomic gases, where the rotation frequency takes over the role of the magnetic field.
In fractional quantum Hall (FQH) physics, the energy scales are set by order of importance: first by the magnetic field, second by the repulsive interactions, and third by the one-body potential. Our discussion in the sequel will reflect this.
A. Landau levels
For large B, it is relevant to restrict particles to live in an eigenspace of . Denote by
the nth Landau level. The lowest level (n = 0) will be denoted by LLL; it is made of analytic × Gaussian functions,
The corresponding fermionic spaces for N particles will be denoted by nLLN and LLLN,
B. Hamiltonians in the LLL
Consider projecting (1) to the LLL. The first term is just a constant, and the others can be expressed using the canonical basis,
Projecting (1) to the LLL leads formally to
where projects21 the relative coordinate ri − rj on the state φm. Similarly, is the operator mapping φℓ to φm, acting on the j variable only.
We assume that the potentials are measurable functions and that the “moments”
are finite for all m. Then, (7) is well defined as a quadratic form on a dense subspace of LLLN. Finiteness for m = 0 means, in particular, that the potentials are in , so derivatives of the potentials are well defined in the sense of distributions. If the moments are uniformly bounded in m and the potentials rotationally symmetric [which implies the absence of terms m ≠ ℓ in (7)], then the corresponding operators are bounded and defined on the whole space.
Usually, in FQH physics, one focuses attention on the interaction term in (7) (i.e., one sets V ≡ 0). There are no off-diagonal terms in it because w is assumed to be radially symmetric. The coefficients are often called “Haldane pseudo-potentials,” (cf. Ref. 22). If w decreases rapidly at infinity, then they also decrease rapidly with increasing m and a basic observation in the theory of the fractional quantum Hall effect (FQHE) is that, if one truncates the sum (7) at m = ℓ − 1, then the Laughlin state
is an exact ground state (L2-normalized by the constant in front). One can then argue and prove to some extent4–6,8 that such functions and natural variants are extremely robust, in particular, to the addition of the external potential V.
For very strong interaction potentials of range much smaller than the magnetic length ∼B−1/2, in particular, if there is a hard core, an expansion in terms of moments as in (7) is not adequate. This situation is analyzed in Ref. 23, which generalizes the paper.24 It is shown that in an appropriate scaling limit, the pseudo-potential operators |φm⟩⟨φm| also emerge, but with renormalized pre-factors involving the scattering lengths of the interaction potentials in the different angular momentum channels, rather than expectation values as in (7).
C. Hamiltonians in higher Landau levels
Consider now a situation where n − 1 Landau levels are filled so that additional electrons must sit in the higher ones because of the Pauli principle. It is a common procedure in the FQH physics community14,17,25 to model this situation using lowest Landau level (LLL) functions again. The basis for this reduction is the following statement, contained in one form or another in a number of sources, in particular, Refs. 1, 3, 14, 17, 26, and 27.
Since we have not assumed any regularity of V and w, except being measurable functions with finite moments, the differentiations in (12) and (13) have, in general, to be understood in the sense of distributions. This poses no problems, however, because the potentials are integrated against densities of wave functions in LLLN, which are smooth functions. Moreover, the densities have the form of polynomial times a Gaussian, so the finiteness of the moments for all m guarantees that the integrals are well defined. In Sec. V B, it will be convenient to assume that the potentials have integrable Fourier transforms, but this is not really an extra restriction because the general case follows by a density argument.
We shall give two proofs of the theorem in Sec. V. Note that the constant we subtract from HnLL is just the magnetic kinetic energy of N particles in the nLL.
What the theorem says is that one can profit from the good properties of the LLL to study phenomena in other Landau levels. This is particularly relevant because the main features are supposed not to depend very much on the potentials Vn, wn entering (7). In particular, the Laughlin states have equivalents in any Landau level (cf. Sec. VI).
Since potential energies are integrals of potentials against particle densities, Theorem II.1 can be seen as a corollary of a general result about particle densities of many body states in different Landau levels. We recall that the k-particle density of an N-particle state with wave function Ψ(r1, …, rN) is, by definition,
If Ψ ∈ nLLN for some n, then is a C∞ function and decreases rapidly at infinity. This is discussed in Sec. V.
Theorem II.1 follows as a corollary if one integrates V(r) against the right-hand side of (16) with k = 1, respectively, w(r1 − r2) with k = 2, and shifts the differentiations to the potentials by partial integration.
Conversely, Theorem II.2 (for k = 1, 2) follows from Theorem II.1 if one regards the potentials as trial functions for the densities.
In the following, we shall define (in several related but distinct ways) the unitary mappings between Landau levels [see (61)] and discuss the proofs of Theorems II.1 and II.2. The physically most appealing way to interpret these unitaries is to see them as replacing the physical coordinates of electrons by the coordinates of the guiding centers of their cyclotron orbits, mathematically implemented through the use of coherent states. Indeed, in the LLL, the position coordinates and the guiding center coordinates are really two different names for the same thing, as will be evident in Sec. IV B. Moreover, the quantum mechanical spread of both coordinates is of the order of the magnetic length ∼B−1/2. The cyclotron radius in the Landau level n has an extra factor . Thus, it is plausible that for large B and small n, the difference between position and guiding center coordinates, and the non-commutativity of the latter, is not of much significance in thermodynamically large systems, i.e., for large N, provided the magnetic length stays much smaller than the interparticle distance.
Although the coherent state approach offers a satisfactory physical picture, it is not always the most convenient one from a computational point of view. This motivates our review of alternate routes to the mappings between levels.
We also take the example of Laughlin states to explain how to deduce the properties of the actual wave functions in nLLN, minimizing effective energies from their representation in the LLLN using the above unitary map. This amounts to saying that the density in guiding center coordinates can to a large extent indeed be identified with the true, physical density in electron coordinates. We believe that this is crucial for the understanding of the efficiency of the correspondence between Landau levels in FQH physics.
III. THE LANDAU HAMILTONIAN AND THE TWO OSCILLATORS
A. The cyclotron oscillator
The magnetic Hamiltonian of a particle of charge q and effective mass m*, moving in a plane with position variables r = (x, y), is
where
is the gauge invariant kinetic momentum with A being the magnetic vector potential and
being the canonical momentum. We assume a homogeneous magnetic field of strength B perpendicular to the plane and choose the symmetric gauge
Moreover, we choose units and signs so that |q| = 1, qB ≡ B > 0, ℏ = 1, and m* = 1. Then,
and the kinetic momentum components satisfy the canonical commutation relations (CCR),
with
being the magnetic length.
In terms of the creation and annihilation operators,
with [a, a†] = 1, the Hamiltonian is
Powers of a† generate normalized eigenstates
with aφ0 = 0 and the energy eigenvalues
In position variables, the corresponding wave functions are
B. Complex notation
With
we can write
Choosing units so that B = 2, or equivalently, defining , this becomes
In addition, the Gaussian factor becomes .
For computations, it is often convenient to use the corresponding operators â†, â, acting on the pre-factors to the Gaussian and defined by
These are
In the sequel, we shall generally use the hat ̂ on operators and functions to indicate that the Gaussian normalization factors are excluded.
Besides the standard definition z = x + iy, other complexifications of are possible and can be useful, as stressed in Ref. 3.
C. The guiding center oscillator
The classical 2D motion of a charged particle in a homogeneous magnetic field consists of a cyclotron rotation around “guiding centers.” The quantization of the cyclotron motion is the physical basis for the energy spectrum (27), and the creation operators a† generate the corresponding harmonic oscillator eigenstates. Every energy eigenvalue is infinitely degenerate due to the different possible positions of the guiding centers.
Quantum mechanically, the dynamics of the guiding centers is described by another harmonic oscillator commuting with the first one. One arrives at this picture by splitting the (gauge invariant) position operator r into the guiding center part R and the cyclotron part
with n being the unit normal vector to the plane. Both R and are gauge invariant, and they commute with each other. On the other hand, the two components of (Rx, Ry) of R do not commute, and likewise for the components of . More precisely, we have
with
and the commutation relations
The creation and annihilation operators for are the same as (24),
Those for the guiding center, on the other hand, are
Note the different signs compared to (39) due to the different signs in (38). We have [b, b†] = 1, and in complex notation,
For B = 2,
and
The splitting (35) corresponds to
While the operators a†, a increase or decrease the Landau level index, the operators b†, b leave each Landau level invariant. Pictorially speaking, we can say that the operators associated with the cyclotron oscillator move states “vertically,” while those associated with the guiding center oscillator move them “horizontally,” i.e., act as ladder operators.
With φ0,0 = φ0 the common, normalized ground state for both oscillators,
the states
form a basis of common eigenstates of the oscillators with φn,0 being the previously defined φn. For fixed n, the states φn,m, m = 0, 1, …, generate the Hilbert space of the nth Landau level, which we shall denote by nLL. The lowest Landau level will be denoted LLL.
Using complex coordinates, the wave functions with n = 0, respectively, m = 0, are
More generally, the wave functions
can be written in terms of associated Laguerre polynomials. They are eigenfunctions of the angular momentum operator in the symmetric gauge (acting on the pre-factor to the Gaussian),
with eigenvalues M = −n + m = −n, −n + 1, … in the nLL. The operators b†, b shift the angular momentum within each Landau level.
IV. EXPRESSIONS OF THE INTER-LEVEL UNITARY MAPS
A. With coherent states
A coherent state associated with the guiding center oscillator in the nLL with parameter is defined in a standard way28,29 as
The overlap of two coherent states is
Moreover,
is the projector on nLL, where is the Lebesgue measure on the plane. Indeed,
and
The coherent states allow an interpretation of nLL as a Bargmann space of analytic functions of the coherent state variable Z: If ψ ∈ nLL, then
is analytic in Z and
has the same L2 norm as ψ because of (52). Thus, the map
is isometric from the nLL to the LLL. From the definition, it is clear that
and
so Un is, in fact, a unitary with
is that used in Theorem II.1.
The function coincides with the LLL wave function of Unψ if Z is identified with the complex position variable z = x + iy. Note, however, that Z is associated with the (non-commutative) components of the guiding center operator R rather than the (commutative) position operator r. By the definition (55), Ψ depends linearly on ψ; the alternative definition Ψ = ⟨ψ|Z, n⟩, which is sometimes used, leads to an anti-unitary correspondence.
B. With integral kernels
Consider the coherent state (50) without the Gaussian normalization factors as as function of ,
The coherent state is an eigenstate of the annihilation operator with eigenvalue Z, so
Furthermore, is an eigenstate of to eigenvalue n, which leads to
with a normalization constant . The full coherent state (50) as a function of Z, z, and , including normalization factors, is thus given by
Inserting this into (55) gives
with
This formula was derived in a different way in Ref. 9 and appears there (in a slightly different notation) as Eq. (34). The inverse map is given by
with
Note that G can be written as
where g is essentially concentrated in a disk of radius and the factor is a phase factor. Recall also that the length unit is .
A further remark is that for n = 0, G is the reproducing kernel in Bargmann space, confirming again that in the LLL, Ψ and ψ are the same function on just with different names for the variables. The phase factor in G is essential for this to hold.
C. With ladder operators
A direct approach to the correspondence nLL ↔ LLL, by-passing the coherent states, starts from (58), noting that
and hence,
Using the representations (33) for the creation and annihilation operators, we conclude that the following holds:
Note that Eq. (75) implies, in particular, that the factor fn(z) to the highest power n of determines uniquely the factors to the lower powers ,
The state is, thus, completely fixed by the holomorphic function fn and the Landau index n.
Incidentally, these considerations also lead to a method for projecting functions to the lowest Landau level.
This is well-known as the recipe “move all factors to the left and replace them by derivatives in z” (see e.g., Ref. 14). For completeness, we give the simple proof:
D. Recap of the different expressions for the unitary maps
Summarizing the contents of this section, we have displayed the following three equivalent ways to represent a state Ψ ∈ nLL by analytic functions in Bargmann space:
Take the scalar product with a coherent state [cf. (55)].
Apply the differential operator -times to the pre-factor of the Gaussian. Equivalently, expand the pre-factor in powers of and keep only the highest power. The inverse mapping, LLL → nLL, is achieved by applying the differential operator to the analytic function representing the state in the LLL.
The last method is formally the simplest, and in Sec. V, we shall use it to discuss particle densities in higher Landau levels in term of their counterparts in the lowest Landau level.
V. PARTICLE DENSITIES AND THE nLL HAMILTONIAN, PROOFS OF THE THEOREMS
We now have all the necessary ingredients to prove Theorems II.1 and II.2. We provide two slightly different approaches.
A. Many body states and particle densities
All considerations in Secs. III and IV carry straightforwardly over to many-body states in symmetric or anti-symmetric tensor powers of single particle states by applying the single particle formulas to each tensor factor.
Let Ψn be a state in nLLN with wave function
Expanding in powers of , we can write
The sum is here over N-tuples (ν1, …νN) such that νk < n for at least one k. The functions fn and are holomorphic, and the latter are, in fact, derivatives of fn [cf. (76)].
The state Ψ0 = UnΨn in LLLN has the wave function
with
The wave function can now be written as
Next, we consider the k-particle density of Ψn, defined by
The density of Ψ0 = UnΨn is given by the same formula with n = 0.
The functions in LLL are holomorphic and decrease at ∞ as ; the latter follows from the fact that the Bargmann kernel (67) with n = 0, which has this decrease, is a reproducing kernel for the Hilbert space LLL. Equation (67) [equivalently, Eq. (73)] also implies that functions in nLL are C∞ in the real position variables and decrease in the same way. This clearly carries over to wave functions in nLLN and the corresponding densities.
To prove Theorem II.2 (which then implies Theorem II.1), we have to compare and . It is, in fact, sufficient to consider the problem for a single variable, i.e., to prove the following lemma:
B. Projected Hamiltonian and guiding center coordinates
We now discuss an alternative road to (12) and (13), providing additional insights. The starting point is the splitting (35) of the position variables in guiding centers and cyclotron motion and the ensuing factorization of matrix elements of exp(iq · r), which enter the Fourier transformed version of (87).
(Plane waves projected in Landau levels). For any , identify eiq·r with the corresponding multiplication operator on , where r is the spatial variable. Let R be the guiding center operator defined in Sec. 3CIII C, Πn be the orthogonal projector on nLL, and Un : nLL → LLL be the inter-LL unitary map.
The equality of the left-hand and right-hand sides of (91) can be seen as a Fourier transformed version of (16) (with k = 1). The identity (91) implies that the norm of Π0eiq·rΠ0 decays faster than any polynomial in |q|. Indeed, on the left-hand side, we have a product of unitaries and projections whose norm is bounded by one. In addition, when q ≠ 0 and n grows, the norm of Πneiq·rΠn decays like n−1/4.
We now provide another proof using guiding center coordinates rather than ladder operators. This also connects with the middle expression in (91).
Some of the following computations can be found in a variety of sources[ e.g., Ref. 14 (Proof of Theorem 3.2) or Ref. 25].
To deduce Eqs. (13) and (12) and, hence, Theorem 2.1, it only remains to write the Fourier decompositions as
and
and use Lemma V.2. The expressions involving Laplacians in Theorem II.1 follow from the Fourier representation −Δ = |q|2. This argument demands absolute integrability of the Fourier transforms, but the general case follows by a density argument.
VI. LAUGHLIN STATES IN HIGHER LANDAU LEVELS
As already mentioned, a crucial approximation in FQH physics is to truncate the Haldane pseudo-potential series in the LLL Hamiltonian (7) to obtain the Laughlin state (9) as an exact ground state of the translation invariant problem V ≡ 0.
In view of Theorem II.1, it is desirable to do the same in a higher Landau level at the level of the effective Hamiltonian and (9) then becomes an exact ground state after the unitary mapping to the LLL. In this section, we explain how the previous considerations allow us to study the properties of the corresponding physical wave function (that is, as expressed in the position coordinates, rather than in the guiding center coordinates).
A. Density estimates on mesoscopic scales
Consider a Laughlin state in the LLLN,
with wave function
Here, ℓ = 1, 3, … for fermions and ℓ = 2, 4, … for bosons. We denote by
the corresponding 1-particle density. According to Laughlin’s plasma analogy,18 the density profile is for large N well approximated by a droplet of radius (ℓN)1/2 and fixed density (πℓ)−1,
Indeed, by a rigorous mean-field analysis, it was proved in Ref. 30 that this approximation holds in the sense of averages over disks of radius Nα with 1/2 > α > 1/4. More generally, the k-particle densities are well approximated in this sense by the k-fold tensor power of the flat density if N → ∞. The more refined analysis of classical Coulomb systems in Refs. 31–33 leads to an extension of this result down to mesoscopic scales Nα for all α > 0. We shall now use results from Ref. 32 to estimate the density of Laughlin states in higher Landau levels.
The Laughlin state corresponding to (108) in the nth Landau level nLLN is
with wave function (cf. Lemma IV.1)
This is, in electronic position variables, the exact ground state of a Hamiltonian obtained by
projecting the physical starting point (1) in the nLLN,
unitarily mapping the result down to an effective Hamiltonian on LLLN using Theorem II.1, and
neglecting the one-body potential Vn and truncating the Haldane pseudo-potential series of the interaction potential wn.
The Hamiltonian obtained this way acts on LLLN, and its exact ground state is a LLL Laughlin state in guiding center variables. Lifting it back up to the nLLN results in (113),
We now vindicate a natural expectation: the density of is very close, for large N, to that of on length scales much larger than the magnetic length (1 in our units). This is because electron coordinates and guiding center coordinates differ only on the scale of a cyclotron orbit, which is much smaller than the thermodynamically large extent of the states themselves.
We shall test the densities with the regularized characteristic functions of disks. Let χ1 be the characteristic function of the unit disc around the origin and for ɛ > 0, and let ηɛ be a function with support in the annulus 1 ≤ |r| ≤ 1 + ɛ such that χ1,ɛ ≔ χ1 + ηɛ is C∞. For R > 0 and , define
In these variables, the extension of the Laughlin state is O(1), and mesoscopic scales are O(N−γ) with 0 < γ < 1/2. The scaled densities are [ is defined in analogy with (110)]
We scale the test functions accordingly and consider . The result on the density and its fluctuations we want to sketch the proof of is as follows:
(Density of Laughlin states on mesoscopic scales).
- For every Landau index n, every fixed ɛ > 0, and all mesoscopic scales r ∼ N−γ with ,(118)
- If r ∼ N−γ, the fluctuation of the linear statistics associated with in the nth Landau level is(119)
B. Rigidity estimates
In Refs. 4–6, 8, 34, and 35, we have investigated the rigidity/stability properties of the LLL Laughlin state. The question is now the response of the Laughlin function to a slight relaxation of the assumptions made in its derivation, namely, that one could in the first approximation neglect the one-body potential and truncate the Haldane pseudo-potential series to a finite order. If one assumes the validity of a certain “spectral gap conjecture” (see Ref. 7, Appendix, and references therein), investigating this question basically means minimizing the one-body energy and the residual part of the interaction within the full ground eigenspace of the truncated interaction energy (cf. degenerate perturbation theory). Our main conclusion was that this problem could be solved to leading order in the large N limit by generating quasi-holes on top of Laughlin’s wave function. We now want to quickly explain how this can be generalized to Laughlin states in higher levels. We discuss only the adaptation of Refs. 5 and 8 for the response to one-body potential. One could consider as well the response to smooth long-range weak interactions as in Ref. 6, but for brevity, we do not write this explicitly.
We take to be a smooth one-body potential, growing polynomially at infinity. We scale it so that it lives on the scale of the Laughlin wave function,
As discussed in the aforementioned references, these assumptions can be relaxed to some extent. The main observation is that after the reduction of the nLL interacting Hamiltonian discussed in Subsection VI A, any multiplication of the LLL Laughlin state by a symmetric analytic function F still yields an exact zero-energy eigenstate in guiding center variables. It is thus relevant to consider the action of the one-body potential VN on the ground-state space of the truncated interaction Hamiltonian. In electron variables, the latter is
where the LLL Laughlin state is the same as in (9). For any many-body wave function , we define its one-particle density as
The variational problem for the response of the Laughlin state to an external potential within the class (126) is now
It is of importance in Laughlin’s theory of the FQHE that one needs to only consider so-called quasi-holes states to solve the above approximately. If one makes this approximation, the minimum energy becomes
The latter energy is obtained by reducing the variational set, so obviously,
What is much less obvious is that this upper bound is optimal in the large N limit.
The n = 0 version of the above was proved in Refs. 5 and 8. The adaptation to higher n follows from the tools therein, together with the representation of discussed at length in Sec. V. We do not give details for brevity. We, however, point out that consequences for minimizing densities also follow so that the density of a (quasi)-minimizer for (128) is approximately flat with value (πℓ)−1 on an open set to be optimized over and quickly drops to 0 outside. This is in accordance with the physical picture of the system responding to external potentials by generating quasi-holes to accommodate their crests. Indeed, the interpretation of the states in (129) is that the zeroes of the analytic function f correspond to the location of quasi-holes in guiding center coordinates.
ACKNOWLEDGMENTS
We had helpful conversations regarding the material of this paper with Thierry Champel, Søren Fournais, and Alessandro Olgiati. We received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation Programme (Grant Agreement CORFRONMAT No. 758620).
REFERENCES
This approach appears already in Ref. 36 where it is attributed to Laughlin.
Although electrons are confined to a 2D interface, they retain their interactions via the 3D Coulomb kernel.
Note that fermionic wave functions do not see the even m terms of (7).