Trial wavefunctions for the Goldstone mode in \nu=1/2+1/2 quantum Hall bilayers

Based on the known physics of the excitonic superfluid or 111 state of the quantum Hall \nu=1/2+1/2 bilayer, we create a simple trial wavefunction ansatz for constructing a low energy branch of (Goldstone) excitations by taking the overall ground state and boosting one layer with respect to the other. This ansatz works extremely well for any interlayer spacing. For small d this is simply the physics of the Goldstone mode, whereas for large d this is a reflection of composite fermion physics. We find hints that certain aspects of composite fermion physics persist to low d whereas certain aspects of Goldstone mode physics persist to high d. Using these results we show nonmonotonic behavior of the Goldstone mode velocity as a function of d.

Based on the known physics of the excitonic superfluid or 111 state of the quantum Hall ν = 1/2 + 1/2 bilayer, we create a simple trial wavefunction ansatz for constructing a low energy branch of (Goldstone) excitations by taking the overall ground state and boosting one layer with respect to the other. This ansatz works extremely well for any interlayer spacing. For small d this is simply the physics of the Goldstone mode, whereas for large d this is a reflection of composite fermion physics. We find hints that certain aspects of composite fermion physics persist to low d whereas certain aspects of Goldstone mode physics persist to high d. Using these results we show nonmonotonic behavior of the Goldstone mode velocity as a function of d.
There are some results, however, that are extremely well established theoretically. In the limit where d becomes very small, it is known that the Halperin 111 trial wavefunction becomes exact [4]. In a more BCS-like language this wavefunction can be expressed as [3] |111 = X (c † X,↑ + c † X,↓ )|0 (1) where ↑ and ↓ indicate the layer index (we assume the real spin is frozen throughout this paper) and X constitutes the orbital index within the Lowest Landau level (chosen to be the x-directed momentum in Landau gauge for example) [20]. The BCS-like form of Eq. 1 allows one to consider long wavelength Goldstone excitations of the form [3,21] These modes are expected to form a linearly dispersing low energy branch with energy proportional to k for small k. Physically, this Goldstone mode corresponds to superflow -one layer being boosted with respect to the other. Both a linearly dispersing mode [22] and excitonic superflow [7,9,[23][24][25] were observed experimentally in this system. Away from the d → 0 limit the form of the bilayer ground state is not known exactly. However, so long as we remain in the same phase of matter, there will continue to be a linearly dispersing Goldstone mode in the long wavelength limit. An approximate expression for this Goldstone mode can be obtained from the ground state wavefunction at any d simply by boosting one layer with respect to the other. One purpose of the current paper is to test this technique of generating trial wavefunctions for the long wavelength Goldstone modes.
We note that this technique is not expected to be exact away from d = 0, but for small d is expected to be quite accurate. In a conventional picture of superfluidity, one might imagine that it would be better to find a way to boost the superfluid fraction while leaving the "normal" fraction unboosted (only at d = 0 is the system entirely super in some sense [13]). Nonetheless, our technique appears to work quite well even away from d = 0.
One might expect that once the system is no longer in the 111 phase of matter (roughly d > 1.5 magnetic lengths), our technique for generating excited states would fail. However, this turns out not to be the case. First of all, at intermediate d there may exist an interlayer paired state as discussed in Ref. 18. Such a paired state would also have a Goldstone mode that could be generated from the ground state by boosting one layer in exactly the same way.
However, even at very large d when such a pairing phase is either absent or pairing is extremely weak, our scheme for generating excited states still works surprisingly well. To understand why this is so, we realize that at large enough d each layer is essentially independent. To first approximation, each layer forms a composite fermion Fermi liquid, which for finite size system has finite momentum (except when the number of electrons exactly fills a shell). The two Fermi liquids are weakly coupled and can combine their momenta to form an overall zero momentum ground state. But since the coupling between the two layers is weak, it costs very little energy to instead form a state of overall finite momentum -which can be interpreted as boosting one layer with respect to the other in comparison to the ground state.
In the absence of any superfluid order parameter (at large d) it is probably not strictly appropriate to refer to this low energy mode as a Goldstone mode. However, since this mode may evolve continuously into the Goldstone mode at smaller d we will abuse nomenclature and continue to call it a Goldstone mode. (If, as conjectured in Ref. 18, the bilayer is actually paired out to large d, then the usage remains correct). Throughout this paper we will work with a spherical geometry. In this case, boosting one layer with respect to the other corresponds to applying the angular momentum raising operator L + to one layer but not the other (call this operator L +,↑ meaning that it is applied to the ↑ layer only). In the appendix we show that if we start with any L = 0 state of the entire system, applying (L +,↑ ) J generates a bilayer state with overall angular momentum L = L z = J. Our technique is then to use exact diagonalization to generate the L = 0 ground state of the bilayer system, which is used to obtain the trial wavefunction for the excited state (3) In turn, we compare this trial state to the exact excited states with angular momentum L = J.
Our numerical work is based on exact diagonalization of the Coulomb Hamiltonian for a bilayer system on the sphere [26]. We simplify the problem to exclude issues related to a finite tunneling amplitude between the layers, Landau-level mixing, or spin (which we assume is polarized), and model each layer as an ideal 2D plane without considering its width into the third dimension. At fillings smaller than one per layer, the Hamiltonian is thus given by the projection of the Coulomb interaction into the lowest Landau-level where sums run over all particles with the given pseudospin. The interactions are parametrized by the layer separation d that is measured in units of the magnetic length 0 = c/eB. All lengths given in this paper should be understood to be measured in units of 0 , where this is not explicitly indicated.
In our exact diagonalization calculations, we focus specifically on the density-balanced bilayer system with N ↑ = N ↓ = N/2 and consider the state at the shift of the 111 state, namely N φ = N − 1. We obtain the two Overlaps of the trial states for the Goldstone mode at d = 0 (top) and d = 3 (bottom) with the exact state, as in Fig. 1, plotted here as a function of wavevector k ∼ L/ N φ /2. Despite the size of the Hilbert spaces increasing significantly between the smallest and largest system shown, the overlap remains roughly constant or maybe slightly increases with N . The region of high overlap extends roughly up to k ≈ 2 −1 0 . The failure of this approach for d = 3 and N = 12 is discussed in detail in the text.
lowest-lying eigenvalues and eigenvectors in each sector of angular momentum. This is most easily achieved using a projected Lanczos algorithm [27] which uses an additional projection to the lowest energy subspace of minimal angular momentum after each multiplication with the Hamiltonian. In a given sector with fixed L z , this procedure therefore directly yields eigenstates of L 2 with the eigenvalue L = L z .
To assess the accuracy of the trial states (3) for the Goldstone mode, we consider their overlap with the lowest energy state in each sector of angular momentum L 2 and for layer separations d = 0 . . . 3 0 , in steps of 1 2 0 [28]. We note again that the the trial states are generated by applying the operator (L +,↑ ) J to the exact groundstate at L = 0. Only at d = 0 is an exact analytical expression of the groundstate known -the 111-state. At other values of d, the numerical groundstates from exact diagonalization are used, although very accurate trial wavefunctions are also known [19]. The results are summarized in Fig. 1, which also indicates overlaps with the first excited state in addition to the overlaps with lowest energy state in each sector. At d = 0, the ansatz (3) is very successful, describing excited states up to high angular momentum accurately. This is shown in more detail in Fig. 2 a), which displays the magnitude of the overlap as a function of wavevector k. Surprisingly, for d = 0 the overlap is very consistent with system size at given k, even though the Hilbert space dimension increases strongly with N . A very good description with overlaps above 0.8 is given up to a wavevector of k ∼ 2 −1 0 . Turning back to Fig. 1, we now focus on the overlaps at finite values of layer sep-aration. At d = 0.5 0 , our trial states obtain significant overlaps only with the first excited state at L = 1. This is not due to a disappearance of the linearly dispersing mode however. Rather, a level crossing appears with distinct excitations occurring at energies less than that of the Goldstone mode, as can be seen from the significant overlaps with the first excited excited state in the sector of L = 2. At d > 1.5 0 this overlap with the first excited state again disappears, signalling the presence of additional low-lying excitations of a nature different from the linearly dispersing Goldstone mode.
Finally at the largest value of d = 3 0 , the ansatz for the boosted trial wavefunctions becomes more accurate than at intermediate d, signalling the possible emergence of a distinct mode of low-lying excitations. Based on the overlap with the trial states, we can point out that there are strong finite size effects in the physics at large d. As with d = 0, Fig. 2 b) displays the numerical values of the overlaps at d = 3. These data single out the system with N = 12 particles as particularly poorly described by these trial states. Here, the physics at large d is clearly dominated by the shell filling effects of composite fermions. As N ↑ = N ↓ = 6 electrons per layer precisely fill the lowest two shells of composite fermion orbitals in one quantum of effective flux, this system size is aliased with the situation where each layer forms its own incompressible ν = 2/5 state with angular momentum zero in each layer -making it impossible to form higher angular momenta states by boosting one layer with respect to each other.
Before proceeding further, note the rather unusual feature that, excepting N = 12, the trial states give higher overlaps for larger systems at d = 3 0 . This unusual behaviour is related to a different manifestation of the composite fermion shell filling effect. As we have shown in Ref. 18, the groundstate at large layer separation is a state in which each layer individually obeys Hund's rule and maximizes the angular momentum per layer L σ = L ↑ = L ↓ , while both layers are combined into a total L = 0 state. Without modifying the correlations inside each layer, the same states with L σ per layer can be paired into excited states with subsequently larger angular momenta, up to a maximum L max = 2L σ . For the system sizes with partially filled composite fermion shells in each layer, one obtains the values of L σ = 3 2 for N = 10, L σ = 5 2 for N = 14, and L σ = 3 2 + 5 2 = 4 for N = 16 particles. We therefore expect a low-lying mode of excitations with angular momenta up to L max = 3, L max = 5, and L max = 8, respectively. Indeed, upon inspection of the spectra, such a mode can be identified.
As an example, Fig. 3  Once the termination of the low energy branch at L max = 8 is identified at large d, it becomes apparent that this feature of a jump in the spectrum at L max exists at all values of layer separation shown, with the exception of the SU(2) invariant case of vanishing d. Note that this termination is a feature of composite fermion physics, explained by successively filling the lowest shells of these composite particles, while obeying Hund's rule. The observation that composite fermion physics intervenes at very small layer separation had been made previously by the current authors. While the 111-state can be regarded as a condensate of composite bosons, it was shown that an accurate description of the groundstate requires a mixedfluid description of both composite bosons and composite fermions at any finite layer separation [16,19]. The identified jump may constitute evidence for the mixed-fluid picture in the excitation spectrum, but more study of these excitation will certainly be required.
A linearly dispersing mode at small k exists at all values of the layer separation. The states which were shown in Fig. 1 to have large overlap with the trial states (3) are highlighted by blue circles in the spectra shown in Fig. 3. These states very accurately come to lie on a single line, which is true especially for the 111-state at d = 0, but also for the intermediate layer separations such as d = 1 0 , where the first excited state at L = 2 lies in the continuation of the line through the points at L = 0 and L = 1 and is shown to be associated to the Goldstone mode by its overlap. Judging by the spectra, it is also very suggestive that multiple level crossings occur at larger values of L, for example at L = 2 for d = 1.5 0 . Finally, between d = 2 0 and d = 3 0 , a change occurs in the association of the low-lying mode with the Goldstone mode trial states. For d = 2 0 only the states up to L = 2 have a good overlap, and the remaining states of the already well formed band of low-lying states in the exact spectrum are of a different nature. The transition to low overlaps occurs at a point where this band has a visible kink and flattens out. Finally at d = 3 0 , this low lying band has good overlaps up to much higher momenta, which, as we have discussed above, is a reflection of two approximately uncoupled composite fermion Fermi seas which each maximizes its own angular momentum according to Hund's rule.
Given the existence of a linearly dispersing mode over the whole range of layer separations, we now consider how its velocity changes with d. Comparing the results obtained for different system sizes, a relatively strong dependence of the velocity v = ∂E/∂k is evident. Applying charging corrections to take account of the shift in the charge-density of the system [29] does not suffice to absorb these effects. Thus, in addition to measuring energies in units of the rescaled magnetic length The velocity has significant finite-size effects on the sphere, therefore we extract an estimate of its value in the thermodynamic limit by extrapolating the finite size values over the inverse system size (see inset).
the thermodynamic limit, as shown in the inset of Fig. 4. The data for N = 12 and large layer separation are easily identified as outliers, due to the shell filling effects discussed above. Generally, the slope is decreasing with system size. However, at small finite layer separations the opposite scaling takes effect. The resulting estimates for the mode velocity in the thermodynamic limit shown in the main graph of Fig. 4 therefore show a pronounced maximum near d = 0 , which has about twice the magnitude as the value found at d = 0. Beyond this point, the mode velocity decreases monotonically and drops to about 1/10 th the value of the 111-state at layer separation d = 3 0 .
A previous experiment [22] has probed the velocity of the neutral mode at layer separations near the transition into the incompressible phase at large d, obtaining a combined best fit of v = 1.4 × 10 4 m/s for data at the three layer separations d 1 = 1.61 0 , d 2 = 1.71 0 and d 3 = 1.76 0 . Based on linear extrapolation between our numerical data at d = 1.5 0 and d = 2 0 , the corresponding estimates are v(d 1 ) = 1.14 × 10 4 m/s, v(d 2 ) = 1.21 × 10 4 m/s, and v(d 3 ) = 1.33 × 10 4 m/s, all slightly smaller but within about 20% of the proposed fit to the experimentally obtained values. Had the data in Ref. 22 been fitted separately at each layer separation, the velocity at d = 1.71 0 would be estimated to be about 10% smaller than that at d = 1.61 0 , roughly reflecting the ratio of our predicted values. The data at d = 1.76 0 appears to be rather noisier, probably due to the vicinity to the phase transition, and would be difficult to fit on its own. We suggest that a significant enhancement of the linear mode velocity should be seen deeper inside the inter-layer coherent phase at smaller layer separation.
To summarize our results, we use the ansatz Eq. 3 to construct trial wavefunctions for a low energy branch of excitations based on the exact ground state wavefunction. This ansatz is accurate at all interlayer spacings d when k is small, and it is accurate at all k when either d is small or d is large (so long as we do not have a filled shell configuration, whereupon only k = 0 is in this low energy branch). We find hints that certain aspects of the composite fermion physics persist to low d whereas certain aspects of the Goldstone mode physics persist to high d. Applying these results to the analysis of our numerical data we show nonmonotonic behavior of the Goldstone mode velocity as a function of the layer separation d. It would be interesting to look for this nonmonotonicity of the Goldstone mode velocity experimentally.
GM acknowledges support from Trinity Hall Cambridge, as well as from the Institute of Complex Adaptive Matter (ICAM). SHS thanks the Aspen Center for Physics for its hospitality.

Appendix: Angular Momentum of Two Coupled Subsystems
In this appendix, we will use the standard angular momentum notation |l, m to indicate eigenstates of the L 2 and L z operators.
We consider two subsystems ↑ and ↓ with corresponding angular momentum operators L ↑ and L ↓ . These two subsystems combine to form the total system with angular momentum operator Our objective is to show that given an eigenstate of the total system with |l = 0, m = 0 application of (L +↑ ) J to this system will produce an eigenstate of the total system with |l = J, m = J . To achieve this, it is sufficient to show that L +↑ |J, J ∼ |J + 1, J + 1 Obviously applying L +↑ to any state increments its overall L z eigenvalue by one (or kills the state), so all that remains is to show that applying L +↑ to |J, J results in an eigenstate L = J +1 of L 2 , i.e, results in an eigenvalue of L 2 being given by (J + 1)(J + 2). Using Eq. 5 it is just a matter of some algebra to show that [L +↑ , L 2 ] = 2L z↑ L +↓ − 2L +↑ L z↓ ) = 2(L + − L +↑ )(L z + 1) − 2L z↓ L + We then apply both sides of this equation to the state |J, J . Noting that L + kills |J, J , we obtain