Supersolid of indirect excitons in electron-hole quantum Hall systems
Abstract
We investigate the ground state of a balanced electron-hole system in the quantum Hall regime using mean-field theory and obtain a rich phase diagram as a function of interlayer distance and the filling factor within a layer. We identify an excitonic condensate phase, an excitonic supersolid phase, as well as uncorrelated Wigner crystal states. We find that balanced electron-hole system exhibits a supersolid phase a wide range of filling factors, with different crystal structure ground states. We obtain the ground state stiffness in the excitonic phases and show that the phase transitions from a uniform condensate to a supersolid is accompanied by a marked change in the stiffness. Our results provide the first semi-quantitative determination excitonic supersolid phase diagram and properties.
Introduction: Over the past decade, Bose-Einstein condensation of indirect excitons Gri in optically pumped and doped electron-hole double quantum wells has been extensively explored but ; sno ; lai ; Butov (2004); Sea The concurrent theoretical analysis has largely focused on either a uniform condensate, where a macroscopic number of excitons occupy the zero-momentum state, or the condensation in a trap where the translational symmetry is broken explicitly she ; Littlewood and Zhu (1996); ga ; pal ; kee . On the other hand, recent observations of the supersolid phase in He Kim ; kim have revived the interest in and questions about excitonic Bose-Einstein condensates with spontaneously broken symmetries and their properties sas (a); Phillips and Balatsky (2007). (“Supersolid” phase of cold atoms in optical lattices has been extensively discussed, although, in that case, the translational symmetry is explicitly broken by the optical lattice.)
Lozovik and Berman Lozovik and Berman (1998) first discussed the coherent charge-density-wave (CCDW) ground state of indirect excitons in electron-hole double quantum wells. Recently, based on general principles, it was shown that electron-hole systems must support a supersolid phase of excitons due to their dipolar repulsion Joglekar et al. (2006). Although a simple qualitative analysis implies the existence of a supersolid phase in electron-hole bilayer systems, the quantitative determination of the phase-boundary and supersolid properties is made difficult by the dispersion of electron (hole) bands.
In this paper, we study an electron-hole system in a strong magnetic field, where the kinetic energy of carriers is quenched ja . We consider a system with equal electron and hole filling factors . CCDW ground states of such a system have been investigated Chen and Quinn (1991, 1992); however, such states are destabilized by fluctuations Côté and Fertig (2000); Brey and Fertig (2000). We obtain the ground state phase diagram in the plane (Fig. 1). Particle-hole symmetry in the lowest Landau level maps the electron-hole system at onto an electron-hole system at and implies that the phase diagram is symmetric around . We verify that this exact relation is satisfied by our mean-field results. The phase diagram of a closely related system - bilayer quantum Hall system where carriers in both layers have the same polarity, near total filling factor - has been extensively explored ja . These systems have provided signatures of excitonic condensation in interlayer tunneling and counterflow experiments at small ja and in the presence of a bias voltage tut .
For an electron-hole quantum Hall system, we find that (Fig. 1): i) for small (region I), the ground state is a uniform excitonic condensate irrespective of the total filling factor ii) for large (region IV), the ground state consists of uncorrelated triangular Wigner crystals in the electron layer and the hole layer iii) for , the ground state is a supersolid of excitons, i.e. a state with spontaneous interlayer phase coherence and spontaneously broken translational symmetry. This supersolid has either a triangular lattice structure (region II) or an anisotropic lattice structure (region III). To study the robustness of phase coherence, we calculate the ground state stiffness by considering the mean-field energy dependence on an in-plane magnetic field. We find that the phase transitions with increasing - from a uniform excitonic condensate to a supersolid to Wigner crystals - are accompanied by marked changes in the stiffness.
In the following section, we briefly sketch the details of mean-field calculations. Then we discuss the phase diagram (Fig. 1) focusing on the excitonic supersolid phases, and the dependence of stiffness on the interlayer distance. We conclude the section with a comment on electron-hole quantum Hall systems at total filling factors and . In the last section, we summarize our results.
Microscopic Model: Let us consider a bilayer system with electrons as carriers in the top layer and holes in the bottom layer, in a uniform magnetic field . The magnetic field normal to the layers quenches the kinetic energy of the carriers, whereas the in-plane magnetic field allows us to study mean-field states with a winding interlayer phase. We assume a vanishing interlayer tunneling amplitude, keeping in mind that the in-plane magnetic field cannot be gauged away for any nonzero interlayer tunneling. The Hamiltonian for the system in the lowest Landau level approximation is given by Chen and Quinn (1991)
(1) |
where denotes electron (hole), is the area of the sample, is the repulsive intralayer Coulomb interaction ( 10 is the dielectric constant of the semiconductor) and is the interlayer attractive Coulomb interaction. The momentum-space density operator in second quantization is
(2) |
where () is the annihilation (creation) operator for a particle in the lowest Landau level state , is the degeneracy of a single Landau level, and is the magnetic length. We define the excitonic condensate operator as
(3) |
and recall that single-particle wave-functions for holes are obtained by complex-conjugation from the single-particle wavefunctions for electrons, . Following standard procedure, we obtain the Hartree-Fock Hamiltonian
(4) |
The first term, contains the intralayer Hartree and exchange, and interlayer Hartree contributions respectively (), and second term denotes the excitonic condensate contribution. These dimensionless contributions are given by , , , and
(5) |
To obtain the self-consistent density matrices, we introduce a two-component operator and define the matrix Green’s function Fetter and Walecka (2003)
(6) |
The electron and hole density matrices as well as the (complex) excitonic order parameter are then determined from the equal-time limit () of this Green’s function matrix. In the Hartree-Fock approximation, the equation of motion for matrix in the frequency space is given by Chen and Quinn (1991); Cot ; Zheng and Fertig (1995)
(7) |
where is the identity matrix, , represents the relative displacement of single-particle wavefunctions in the two layers due to the Lorenz drift caused by the in-plane field Zheng and Fertig (1995), and we have defined a modified Green’s function
(8) |
The self-energies in Eq. (Supersolid of indirect excitons in electron-hole quantum Hall systems) are given by with . The modified Green’s function can be expressed as
(9) |
where are the eigenvectors of the self-energy matrix with eigenvalue ,
(10) |
We notice that, in contrast to a similar analysis for bilayer quantum Hall systems Cot ; Zheng and Fertig (1995) - the chemical potential appears with opposite signs in the electron and hole self-energies. The Hartree-Fock ground state energy is calculated using the self-consistent density matrix using Eq.(4).
Equations (8)-(9) provide the set of equations that are iteratively solved. We check that the resulting density matrix satisfies the sum-rule Cot
(11) |
Results: In order to determine the phase diagram in
the plane, we compare the energies of uniform excitonic
condensate (, ), Wigner
crystal in each layer (,), and supersolid (,
) in the absence of an
in-plane magnetic field . We use simplified anisotropic
lattice with two primitive lattice vectors
citebrey2000
and , and define the lattice anisotropy
. The lattice constant is determined by the constraint
that the unit cell contains one electron-hole pair, and the optimal
value of is obtained by minimizing the mean-field energy.
Note that the triangular lattice ( 1.155) and
a stripe lattice () are its special cases. (We
have found that square lattices have higher energy in the supersolid
and uncorrelated Wigner crystal phases.)
Figure 1 shows the ground state phase diagram. At small values of , we find a uniform excitonic condensate as the ground state over the entire range of total filling factor (region I). For such a state, the mean-field equations can be analytically solved and we obtain the layer densities and the excitonic condensate order parameter . These analytical results satisfy Eq. (11) and imply that the uniform phase coherence order parameter is the same for systems with filling factors and . The monotonic increase in with decreasing Chen and Quinn (1991) is consistent with observed strengthening of excitonic condensate phase in bilayer quantum Hall systems with layer imbalance tut . At large values of (which is the same as for ), the ground state is uncorrelated triangular Wigner crystals in the electron layer and the hole layer (region IV).
For , we find that the ground state is a supersolid at intermediate values of distance . In region II, the supersolid has a triangular lattice structure () whereas in region III, the lattice is anisotropic. (The anisotropic supersolid with small is similar to the CCDW solutions discussed in Ref.Chen and Quinn (1991)). The boundary between two regions corresponds to a line of first-order phase transitions from a triangular lattice to an anisotropic lattice. The inset in Fig. 1 shows the ground state lattice anisotropy along the line . For , is constant, as expected for a triangular lattice; near , drops discontinuously and then reduces monotonically with . This discontinuity in indicates a first order transition from a triangular to an anisotropic supersolid. The reverse transition, from region III to region II, can be induced by increasing the interlayer distance at a given value of . Figure 2 shows the mean-field energy as function of lattice anisotropy for increasing . We see that the optimal value of jumps discontinuously from (anisotropic lattice) to (triangular lattice). The triangular lattice structure in regions II and IV suggests that the phase transition between them is a continuous phase transition. On the other hand, different lattice structures in region III and IV imply that the transition between them will be a first order transition.
An important property of a Bose-Einstein condensate, uniform or otherwise, is the energy cost associated with spatial variations of its phase. A uniform in-plane magnetic field induces such a variation in a uniform excitonic condensate sas (b). We calculate the stiffness, defined as the variation of ground state energy with an in-plane magnetic field not
(12) |
from the ground state energy. Figure 3 shows dependence of the stiffness on interlayer distance at We see a change in the stiffness slope at , where the system changes from a uniform excitonic condensate to a triangular supersolid and , where it changes from the supersolid to uncorrelated Wigner crystals, along with the expected monotonic decrease with . Since the transport properties (in superfluid or supersolid phase) depend on the stiffness, generically, we expect that the phase transitions will be visible in transport measurements.
An electron-hole system at is mapped on to a system at after a particle-hole transformation on both layers. Figure 4 shows that the mean-field density matrices, obtained separately for and , indeed satisfy and . Therefore, it is enough to restrict ourselves to .
Discussion: The supersolid phase, with two spontaneously broken continuous symmetries, has been a source of extensive investigations in He Kim ; kim ; sas (a), but has not been experimentally realized in any other system. Our results predict that electron-hole quantum Hall systems exhibit a rich phase diagram including a supersolid phase robust over filling factors . We show that the supersolid phase exhibits a nonzero stiffness and obtain the dependence of stiffness on interlayer distance. It is well known that mean-field approach overestimates the stability ordered states. However, it is not clear whether it overestimates the stability of one type of order (excitonic condensation) over another (broken translational symmetry); therefore, we believe that our conclusions - existence of the supersolid phase in the intermediate distance regime and the accompanying change in the phase stiffness - are generically valid.
Our results present an example of a supersolid whose properties can be investigated starting from a microscopic Hamiltonian; in particular, the study of low-energy excitations and transport properties in the supersolid phase is possible. We point out that this excitonic system, although analytically tractable, is significantly different from the zero-field system in which, at low exciton densities, the excitons behave as non-interacting bosons. Nonetheless, an experimental verification (or falsification) of our predictions will provide a better understanding of the supersolid phase and the properties of excitons in the quantum Hall regime.
Acknowledgments: It is a pleasure to thank Allan MacDonald for useful discussion.
References
- (1) L. V. Keldysh, in Bose-Einstein Condensation, edited by A. Griffin, D. W. Snoke, and S. Stringari (Cambridge University Press, Cambridge, UK, 1995) and references therein.
- (2) L.V. Butov, A.C. Gossard, and D.S. Chemla, Nature (London) 417, 47 (2002); ibid. 418, 751 (2002).
- (3) D. Snoke, S. Denew, Y. Liu, L.N. Pfeiffer, and K.W. West, Nature (London) 418, 754 (2002).
- (4) C.W. Lai, J. Zoch, A.C. Gossard, and D.S. Chemla, Science 303, 503 (2004).
- Butov (2004) L. V. Butov, J. Phys.: Condens. Matter 16, R1577 (2004).
- (6) J. A. Seamons et al., App. Phys. Lett. 90, 052103 (2007).
- (7) S.I. Shevchenko, Phys. Rev. Lett. 72, 3242 (1992).
- Littlewood and Zhu (1996) P. B. Littlewood and X. Zhu, Phys. Scr. T68, 56 (1996).
- (9) S. Vignale and A.H. MacDonald, Phys. Rev. Lett. 76, 2786 (1996).
- (10) S.De Palo, F. Rapisarda, and G. Senatore, Phys. Rev. Lett. 88, 206401 (2002).
- (11) J. Keeling, L.S. Levitov, and P.B. Littlewood, Phys. Rev. Lett. 92, 176402 (2004).
- (12) E. Kim and M. W. H. Chan, Nature 427, 225 (2004); Science 305 1941(2004).
- (13) X. Lin, A.C. Clark, and M.W.H. Chan, Nature 449, 1025 (2007).
- sas (a) S. Sasaki et al., Science 313, 1098 (2006).
- Phillips and Balatsky (2007) P. Phillips and A. V. Balatsky, Science 316, 1435 (2007).
- Lozovik and Berman (1998) Y. E. Lozovik and O. L. Berman, Physica Scripta 58, 86 (1998).
- Joglekar et al. (2006) Y. N. Joglekar, A. V. Balatsky, and S. Das Sarma, Phys. Rev. B 74, 233302 (2006).
- (18) J.P. Eisenstein and A.H. MacDonald, Nature (London) 432, 691 (2004); J.P. Eisenstein, Science 305, 950 (2004) and references therein.
- Chen and Quinn (1991) X. M. Chen and J. J. Quinn, Phys. Rev. Lett. 67, 895 (1991).
- Chen and Quinn (1992) X. M. Chen and J. J. Quinn, Phys. Rev. B 45, 11 054 (1992).
- Côté and Fertig (2000) R. Côté and H. A. Fertig, Phys. Rev. B 62, 1993 (2000).
- Brey and Fertig (2000) L. Brey and H. A. Fertig, Phys. Rev. B 62, 10268 (2000).
- (23) E. Tutuc et al., Phys. Rev. Lett. 91, 076802 (2003).
- Fetter and Walecka (2003) A. L. Fetter and J. D. Walecka, Quantum Theory of Many-Particle System (Dover Publications, Inc., NewYork, 2003).
- (25) R. Côté and A. H. Macdonald, Phys. Rev. Lett. 65, 2662 (1990); Phys. Rev. B 44, 8759 (1991).
- Zheng and Fertig (1995) L. Zheng and H. A. Fertig, Phys. Rev. B 52, 12282 (1995).
- sas (b) A.V. Balatsky, Y.N. Joglekar, and P.B. Littlewood, Phys. Rev. Lett. 93, 266801 (2006).
- (28) This stiffness, for a spatially homogeneous system, corresponds to the cost of phase fluctuations. Our calculations, with a supersolid ground state, include contributions from the phase variations and relative density displacements.