1. Introduction
The many-particle systems of dipolar (indirect) excitons, formed by spatially separated electrons and holes, in semiconductor coupled quantum wells (CQWs) and novel two-dimensional (2D) materials have been the subject of numerous experimental and theoretical studies. These systems are attractive in large part due to the possibility of Bose-Einstein condensation (BEC) and superfluidity of dipolar excitons, which can be observed as persistent electrical currents in each quantum well, and also through coherent optical properties [
1,
2,
3,
4,
5]. Recent progress in theoretical and experimental studies of BEC and superfluidity of dipolar excitons in CQWs have been reviewed in [
6]. Electron-hole superfluidity in double layers can occur not only in the BEC regime, but also in the Bardeen-Cooper-Schrieffer (BCS)-BEC crossover regime [
7].
A number of experimental and theoretical investigations have been devoted to the BEC of electron-hole pairs, formed by spatially separated electrons and holes in a double layer formed by parallel graphene layers. These investigations were reported in [
8,
9,
10,
11,
12,
13]. Both BEC and superfluidity of dipolar excitons in double layers of transition-metal dichalcogenides (TMDCs) [
14,
15,
16,
17,
18] and phosphorene [
19,
20] have been discussed, because the exciton binding energies in novel 2D semiconductors are quite large. Possible BEC in a long-lived dark spin state of 2D dipolar excitons has been experimentally observed for GaAs/AlGaAs semiconductor CQWs [
21].
Recently, the electronic properties of the
lattice have been the subject of the intensive theoretical and experimental investigations due to its surprising fundamental physical properties as well as its promising applications in solid-state devices [
22,
23,
24,
25,
26,
27,
28,
29,
30,
31,
32,
33,
34,
35]. For a review of artificial flat band systems, see [
36]. Raoux, et al. [
22] proposed that an
lattice could be assembled from cold fermionic atoms confined to an optical lattice by means of three pairs of laser beams for the optical dice lattice (
) [
37]. This structure consists of an AB-honeycomb lattice (the rim) like that in graphene which is combined with C atoms at the center/hub of each hexagon. A parameter
represents the ratio of the hopping integral between the rim and the hub to that around the rim of the hexagonal lattice. By dephasing one of the three pairs of laser beams, one could vary the parameter
. Optically induced dressed states [
38], and their tunneling, transport [
33,
39], and collective properties [
40], as well as
based nanoribbons [
41] have been analyzed. The BEC and superfluidity of dipolar magnetoexcitons in
double layers in a strong uniform perpendicular magnetic field were proposed in [
42].
We present the conditions for BEC and superfluidity of a two-component weakly interacting Bose gas of dipolar excitons, formed by electron-hole pairs in spatially separated GHAT3 layers. An applied weak magnetic field to this pseudospin-1 monolayer system results in a Zeeman-type splitting of the energy subbands [
43]. This dispersion relation consists of three bands: CB, IB, and VB. We consider two types of dipolar excitons in a double-layer of GHAT3: (a) “A excitons”, formed as bound states of electrons in CB and holes in IB, and (b) “B excitons”, formed as bound states of electrons in CB and holes in VB. The binding energy of A and B dipolar excitons is calculated. For a two-component weakly interacting Bose gas of dipolar excitons in a GHAT3 double layer, we obtain the energy dispersion of collective excitations, the sound velocity, the superfluid density, and the mean-field critical temperature
for superfluidity.
Our paper is organized in the following way. In
Section 2, the two-body problem for an electron and a hole, spatially separated in two parallel GHAT3 monolayers, is formulated, and the effective masses and binding energies are obtained for two types of dipolar excitons. The spectrum of collective excitations and the sound velocity for the two-component weakly interacting Bose gas of dipolar excitons in the double layer of GHAT3 are derived in
Section 3. In
Section 4 the superfluidity of the weakly interacting Bose gas of dipolar excitons in the double layer of GHAT3 is predicted, and the mean-field critical temperature of the phase transition is obtained. The results of our calculations are discussed in
Section 5. In
Section 6 our conclusions are reported.
2. Dipolar Excitons in a Double Layer of with a Mass Term
We will consider charge carriers in the conduction band, valence band, and the intermediate band, which corresponds to the flat band in an
layer without a mass term. In the presence of a weak magnetic field, the low-energy Hamiltonian of the charge carriers in a GHAT3 monolayer at the K and K’ points are given by [
43]
where the origin in
-space is defined to be around the K point,
and
,
,
, with
being the valley index at the K and K’ points,
is the gap in the energy spectrum of a GHAT3 layer due to the mass term in the Hamiltonian. In an
layer honeycomb lattice, there is an added fermionic hub atom C at the center of each hexagon. Let the hopping integral be
between the hub atom and either an A or B atom on the rim and
between nearest neighbors on the rim of the hexagon. The ratio of these two nearest neighbor hopping terms is denoted as
, where the parameter
satisfies
. The largest value when
is 1 is for the dice lattice, whereas its value of 0 corresponds to graphene for decoupled hub from rim atoms [
43].
At small momenta near K and K’ points, the dispersion for the charge carriers in the conduction band
is given by the relation [
43]
where
and
are the wave vector and momentum of a quasiparticle,
is the effective mass of the charge carriers in the conduction band, given by
where
is the Fermi velocity in a GHAT3 layer, and
[
43]. At small momenta near K and K’ points, the dispersion for the charge carriers in the valence band
is given by the relation [
43]
with
the effective mass of the charge carriers in the valence band, given by
At small momenta near K and K’ points, the dispersion for the charge carriers in the intermediate band, corresponding to the flat band in an
layer without a mass term,
is given by the relation [
43]
where
is the effective mass of the charge carriers in the intermediate band, given by
It is worth noting that there are spin degeneracy and valley degeneracy for the energy of the charge carriers in a GHAT3 layer.
In the system under consideration in this paper, electrons are confined in a 2D GHAT3 monolayer, while an equal number of positive holes are located in a parallel GHAT3 monolayer at a distance
D away as demonstrated in
Figure 1. This electron-hole system in two parallel GHAT3 layers is treated as a 2D system without interlayer hopping. Due to the absence of tunneling of electrons and holes between different GHAT3 monolayers, electron-hole recombination is suppressed by a dielectric barrier with dielectric constant
that separates the GHAT3 monolayers. Therefore, the dipolar excitons, formed by electrons and holes, located in two different GHAT3 monolayers, have a longer lifetime than direct excitons. The electron and hole are attracted via electromagnetic interaction
where
is the distance between the electron and hole, and they could form a bound state, i.e., an exciton, in three-dimensional (3D) space. Therefore, to determine the binding energy of the exciton a two-body problem in restricted 3D space has to be solved. However, if one projects the electron position vector onto the GHAT3 plane with holes and replaces the relative coordinate vector
by its projection
on this plane, the potential
may be expressed as
where
r is the relative distance between the hole and the projection of the electron position vector onto the GHAT3 plane with holes. A schematic illustration of the dipolar exciton in a GHAT3 double layer is presented in
Figure 1. By introducing in-plane coordinates
and
for the electron and the projection vector of the hole, respectively (where
), the dipolar exciton can be described by employing a two-body 2D Schrödinger equation with potential
. So that the restricted 3D two-body problem can be reduced to a 2D two-body problem on a GHAT3 layer with the holes.
The dipolar excitons with spatially separated electrons and holes in two parallel GHAT3 monolayers can be created by laser pumping with an applied external voltage. While an electron in the conduction band and a hole in the valence or intermediate band are excited due to absorption of a photon, voltages are applied with opposite signs to confine electrons on one layer and holes on another so that dipoles point in one direction only.
In our case, “both” the energy bands and the exciton modes referred to the K-point, not one to the
point and the other to the K point. We note that in the dispersion equations appearing in [
43,
44,
45,
46] the origin of the k-space was specified to be around the K point, (and not the
point) as did several authors investigating
. So, our choice of origin not being the center of the Brillouin zone has precedence. For graphene, the plasmon dispersion relation and low-energy bands, presented by [
47] were both consistently measured from the K point taken as the origin and not the center of the Brillouin zone.
We consider excitons, formed by an electron and a hole from the same valley, because an electron and a hole from different valleys cannot be excited by absorption of photon due to conservation of momentum. The reason is that photons carry momenta much smaller than the difference between K and K’ in reciprocal space.
The effective Hamiltonian of an electron and a hole, spatially separated in two parallel GHAT3 monolayers with the interlayer distance
D has the following form
where
and
are the Laplacian operators with respect to the components of the vectors
and
, respectively, and
and
are the effective masses of the electron and hole, respectively. For CV excitons
and
; and for CI excitons
and
, where
,
, and
are given by Equations (
3), (
5) and (
7), correspondingly. The problem of the in-plane motion of an interacting electron and hole forming the exciton in a GHAT3 double layer can be reduced to that of one particle with the reduced mass
in a
potential and motion of the center-of-mass of the exciton with the mass
. We introduce the coordinates of the center-of-mass
of an exciton and the coordinate of the relative motion
of an electron and hole as
and
, correspondingly. The Hamiltonian
can be represented in the form:
, where the Hamiltonian of the motion of the center-of-mass is
and that of the relative motion of electron and a hole is
. The solution of the Schrödinger equation for the center-of-mass of an exciton
is the plane wave
with the quadratic energy spectrum
, where
is the momentum of the center-of-mass of an exciton.
We consider electrons and holes to be located in GHAT3 parallel layers, embedded in a dielectric with the dielectric constant
. The potential energy of electron-hole Coulomb attraction is
where
,
is the dielectric constant of the insulator (SiO
or
h-BN), surrounding the electron and hole GHAT3 monolayers, forming the double layer. For the
h-BN barrier we substitute the dielectric constant
, while for the SiO
barrier we substitute the dielectric constant
. For
h-BN insulating layers,
is the effective dielectric constant, defined as
[
14], where
and
are the components of the dielectric tensor for
h-BN [
48]. Assuming
, we approximate
by the first two terms of the Taylor series and obtain
The similar approach has been applied for excitons in TMDC double layers [
16,
17]. The solution of the Schrödinger equation for the relative motion of an electron and a hole
with the potential (
10) is reduced to the problem of a 2D harmonic oscillator with the exciton reduced mass
. Following [
49,
50] one obtains the radial Schrödinger equation and the solution for the eigenfunctions for the relative motion of an electron and a hole in a GHAT3 double layer in terms of associated Laguerre polynomials, which can be written as
where
,
,
are the quantum numbers,
is the polar angle, and
is a Bohr radius of a dipolar exciton. The corresponding energy spectrum is given by
At the lowest quantum state
as it follows from Equation (
12) the ground state energy for the exciton is given by
The important characteristic of the exciton is the square of the in-plane gyration radius
. It allows one to estimate the condition when the excitonic gas is dilute enough. One can obtain the square of the in-plane gyration radius
of a dipolar exciton [
14], which is expressed as the average squared projection of an electron-hole separation onto the plane of a GHAT3 monolayer
We consider dipolar excitons, formed by an electron in the conduction band and a hole in the valence band (CV excitons) and formed by and electron in the conduction band and a hole in the intermediate valence band (CI excitons). For CV excitons one has
3. The Collective Excitations Spectrum and Superfluidity for the Two-Component System of Dipolar Excitons
We consider the dilute limit for dipolar exciton gas in a GHAT3 double layer, when
and
, where
and
are the concentration and effective exciton Bohr radius for A(B) dipolar excitons, correspondingly. In the dilute limit, dipolar A and B excitons are formed by electron-hole pairs with the electrons and holes spatially separated in two different GHAT3 layers. We will treat the two-component weakly interacting Bose gas of dipolar excitons in a GHAT3 double layer by applying the approach analogous to the one used for dipolar excitons in a transition metal dichalcogenide (TMDC) double layer [
16,
17].
Since the dipolar excitons, formed by the charge carriers in different valleys, are characterized by the same energy, the exciton states are degenerate with respect to the valley degree of freedom. Therefore, we consider the Hamiltonian of the weakly interacting Bose gas of dipolar excitons, formed in a single valley. We will take into account the degeneracy of the exciton states with respect to spin and valley degrees of freedom by the introducing the spin and valley degeneracy factor
below. The Hamiltonian
of the 2D A and B weakly interacting dipolar excitons can be written as
where
are the Hamiltonians of A(B) excitons defined as
and
is the Hamiltonian of the interaction between A and B excitons presented as
where
and
are Bose creation and annihilation operators for A(B) dipolar excitons with the wave vector
, correspondingly,
S is the area of the system,
is the energy spectrum of non-interacting A(B) dipolar excitons, respectively,
,
is an effective mass of non-interacting dipolar excitons,
is the constant, which depends on A(B) dipolar exciton binding energy and the corresponding gap,
and
are the interaction constants for the repulsion between two A dipolar excitons, two B dipolar excitons and for the interaction between A and B dipolar excitons, respectively.
In dilute system with large interlayer separation
D, two dipolar excitons, located at distance
R, repel each other via the dipole-dipole interaction potential
. Following the procedure described in [
51], the interaction parameters for the exciton-exciton repulsion in very dilute systems can be obtained implying the exciton-exciton dipole-dipole repulsion exists only at the distances between excitons greater than the distance from the exciton to the classical turning point.
The many-particle Hamiltonian for a weakly interacting Bose gas can be diagonalized within the Bogoliubov approximation [
52], replacing the product of four operators in the interaction term with the product of two operators. The Bogoliubov approximation is valid if one assumes that most of the particles belong to BEC. In this case, in the Hamiltonian one can keep only the terms responsible for the interactions between the condensate and non-condensate particles, while the terms describing the interactions between non-condensate particles are neglected.
Following the procedure, described in [
16,
17], applying the Bogoliubov approximation [
52], generalized for a two-component weakly interacting Bose gas [
53,
54] and introducing the following notation,
one obtains two modes of the spectrum of Bose collective excitations
where
, 2. In our approach, the condition
holds.
At small momenta
, when
and
, expanding the spectrum of collective excitations
up to the first order with respect to the momentum
p, one obtains two sound modes in the spectrum of the collective excitations
, where
is the sound velocity written as
At , the spectrum of collective excitations is determined by the non-zero sound velocity , while at the sound velocity vanishes with . At large momenta, for the conditions when and , one obtains two parabolic modes of collective excitations with the spectra and , if and if with the spectra and .
4. Superfluidity of the Weakly-Interacting Bose Gas of Dipolar Excitons
Since when
the sound velocity vanishes, below we take into account only the branch of the spectrum of collective excitations at
, neglecting the branch at
. According to [
52,
55], it is clear that we need a finite sound velocity for superfluidity. Since the branch of the collective excitations at zero sound velocity for the collective excitations corresponds to the zero energy of the quasiparticles (which means that no quasiparticles are created with zero sound velocity), this branch does not lead to the dissipation of energy resulting in finite viscosity and, therefore, does not influence the Landau critical velocity. This is the reason for eliminating the zero sound velocity case in our considerations here. The weakly-interacting gas of dipolar excitons in the double layer of GHAT3 satisfies the Landau criterion for superfluidity [
52,
55], because at small momenta, the energy spectrum of the quasiparticles in the weakly-interacting gas of dipolar excitons at
is sound-like with the finite sound velocity,
. In the moving weakly-interacting gas of dipolar excitons the quasiparticles are created at velocities above the velocity of sound, and the critical velocity for superfluidity reads as
. The difference between the ideal Bose gas and two-component weakly interacting Bose gas of dipolar excitons is that while the spectrum of ideal Bose gas has no branch with finite sound velocity, the dipolar exciton system under consideration has one branch in the spectrum of collective excitations with finite sound velocity at
due to exciton-exciton interaction. Therefore, at low temperatures, the two-component system of dipolar excitons exhibits superfluidity due to exciton-exciton interactions, while the ideal Bose gas does not demonstrate superfluidity.
We defined the density of the superfluid component
as
, where
is the total 2D density of the dipolar excitons and
denotes the density of the normal component. The density
of the normal component can be defined using standard procedure [
56]. The assumption that the dipolar exciton system moves with a velocity
implies that the superfluid component moves with the velocity
. The energy dissipation at nonzero temperatures
T is characterized by the occupancy of quasiparticles in this system. Since the density of quasiparticles is small at low temperatures, the gas of quasiparticles can be treated as an ideal Bose gas. In order to obtain the density of the superfluid component, one can define the total mass flow for a Bose gas of quasiparticles in the frame, in which the superfluid component is assumed to be at rest, as
where
is the spin and valley degeneracy factor,
is the Bose-Einstein distribution function for the quasiparticles with the dispersion
, and
is the Boltzmann constant. Expanding the expression under the integral in Equation (
23) up to the first order with respect to
, one has:
The density
of the normal component in the moving weakly-interacting Bose gas of dipolar excitons is defined as [
56]
Employing Equations (
24) and (
25), one derives the normal component density as
At low temperatures
, the small momenta (
and
) make the dominant contribution to the integral on the right-hand side of Equation (
26). The quasiparticles with such small momenta are characterized by the sound spectrum
with the sound velocity defined by Equation (
22). By substituting
into Equation (
26), we obtain
where
is the Riemann zeta function (
).
The mean field critical temperature
of the phase transition at which the superfluidity occurs, implying neglecting the interaction between the quasiparticles, is obtained from the condition
[
56]:
At low temperatures
by substituting Equation (
27) into Equation (
28), one derives
While Bose-Einstein condensation occurs at absolute zero even in a two-dimensional (2D) system, it is well known that in a 2D bosonic system, Bose-Einstein condensation does not occur at finite temperature, and only the quasi-long-range order appears. In this paper, we have obtained the mean-field critical temperature
of the phase transition at which superfluidity appears without claiming BEC in a 2D system at finite temperature. In this work, we have considered BEC only at absolute zero temperature. The similar approach has been applied for excitons in TMDC double layers [
16,
17].
5. Discussion
In this section we now discuss the results of our calculations. In
Figure 2, we present the results for the exciton binding energy
for CV and CI excitons as functions of the gap
for chosen parameter
and interlayer separations
. According to
Figure 2,
is an increasing function of
, whereas for CV excitons the exciton binding energy is slightly larger than that for CI excitons.
In
Figure 3, we present our results for the exciton binding energy
for CV and CI excitons as functions of the parameter
for chosen gap
and interlayer separations
. According to
Figure 3,
does not depend on
for CV excitons, whereas it is an increasing function of
for CI excitons. At
for CV excitons is larger than for CI excitons, while at
for CI excitons is larger than for CV excitons.
In
Figure 4, we present the results of our calculations for the exciton binding energy
for CV and CI excitons as functions of the interlayer separation
D for chosen parameter
and gap
. According to
Figure 4,
is a decreasing function of
D, whereas for CV excitons the exciton binding energy is slightly larger than for CI excitons.
In
Figure 5, we present plots of the effective masses for CV and CI dipolar excitons as functions of the gap
for chosen
for (a) center-of-mass exciton mass
M on the left-hand side and (b) reduced exciton mass
, on the right. According to
Figure 5, both
M and
for the CV and CI excitons are increasing functions of
, while for CV excitons both
M and
are slightly larger than for CI excitons.
Figure 6 shows the effective masses of a dipolar exciton for CV and CI excitons as functions of
for chosen
for (a) center-of-mass exciton mass
M in the left panel and (b) reduced exciton mass
, on the right. According to
Figure 6, for CV excitons
M is a decreasing function of
, whereas
does not depend on
. For CI excitons, both
M and
increase as
is increased. For
, both
M and
for CV excitons are larger than for CI excitons, but when
both
M and
for CV excitons are smaller than for CI excitons.
Figure 7 demonstrates the dependence of the sound velocity
on the hopping parameter
for chosen
, interlayer separations
at fixed concentrations
and
of A and B excitons, respectively. According to
Figure 7,
c does not depend much on
when
, while for
, the sound velocity
c is a decreasing function of
.
In
Figure 8, we plot the sound velocity
versus the gap
for chosen parameter
, interlayer separations
for chosen concentrations
and
of A and B excitons, respectively. According to
Figure 8, the sound velocity
c is a decreasing function of
.
In
Figure 9, we show the sound velocity
as a function of the interlayer separation
D for hopping parameter
and gap
, for fixed concentrations
and
of A and B excitons, respectively. According to
Figure 9, the sound velocity
c is an increasing function of
D.
In
Figure 10, we illustrate the dependence of the sound velocity
on the concentrations
and
of A and B excitons, respectively for chosen hopping parameter
and gap
, at fixed interlayer separation
. According to
Figure 10, the sound velocity
c is an increasing function of both concentrations
and
.
In
Figure 11, we present the mean-field phase transition critical temperature
as a function of the parameter
for chosen gap
, interlayer separations
at the fixed concentrations
and
of A and B excitons, respectively. According to
Figure 11,
is a decreasing function of
at
, while at
the critical temperature
is an increasing function of
.
In
Figure 12, we present the mean-field phase transition critical temperature
as a function of the gap
for chosen parameter
, interlayer separations
at the fixed concentrations
and
of A and B excitons, respectively. According to
Figure 12, the criticaltemperature
is a decreasing function of
.
In
Figure 13, we demonstrate the mean-field phase transition critical temperature
as a function of the interlayer separation
D for chosen parameter
and gap
, at the fixed concentrations
and
of A and B excitons, respectively. According to
Figure 13, the critical temperature
is an increasing function of
D.
In
Figure 14, we present density plots for the mean-field phase transition critical temperature
as a function of the concentrations
and
of A and B excitons, respectively for chosen parameter
and gap
, at the fixed interlayer separation
. According to
Figure 14, the the critical temperature
is an increasing function of both the concentrations
and
.
At a formal level, the weakly interacting Bose gas of A and B dipolar excitons in a GHAT3 double layer are similar to the two-component weakly interacting Bose gas of trapped cold atoms in a planar harmonic trap. The spectrum of collective excitations in the Bogoliubov approximation for dipolar excitons in a GHAT3 double layer is similar to one for a two-component BEC of trapped cold atoms, studied in [
53,
54].
The gap parameter
, has a dual role, since it appears as chemical potential in the Hamiltonian, as also in the mass of the excitons through the band curvature. According to
Figure 2 and
Figure 12, the dipolar exciton binding energy is an increasing function of the gap
, while is the mean-field phase transition temperature
is a decreasing function of the gap
. Therefore, there should be an optimal value for
, which would correspond to relatively high
at the relatively high dipolar exciton binding energy. The latter condition provides the formation of the superfluid phase by the relatively stable dipolar excitons.
Note that electron-hole superfluids can be formed not only in the BEC regime but also in the BCS-BEC crossover regime [
7]. Quantum Monte Carlo simulations analyzing the BCS-BEC crossover regime for electron-hole systems have been performed [
57]. In this paper, we concentrate on the dilute electron-hole system, which corresponds to the BEC, which matches experimentally achievable densities in the electron-hole systems in 2D materials. BCS regime requires higher concentrations beyond the model of weakly interacting Bose gas. The studies of the BCS regime, and BEC-BCS crossover for an electron-hole superfluid in a GHAT3 double layer, seem to be a promising direction for future studies.
The considered system of dipolar excitons in a GHAT3 double layer has also a strong similarity, with photon condensation in a cavity. The collective modes and possibility of the Kosterlitz-Thouless phase transition to the superfluid phase [
58,
59] has been studied for a photon condensation in a cavity in [
60]. If we consider only one type of excitons in a GHAT3 double layer, assuming the concentration of the excitons of another type to be zero, the expressions for the spectrum of collective excitations reported in this paper can be reduced to the expressions similar to [
60].
The Kosterlitz-Thouless phase transition to the superfluid phase [
58,
59] can be inferred from the variation of the superfluid density, which has been computed in this paper.
Note that, in this paper, we did not consider vortices, as within the mean-field approximation it was assumed that the number of quasiparticles are relatively low. However, beyond the mean-field approximation, it is possible to consider the properties of vortices in the system of dipolar excitons. Thus, the dynamical creation of fractionalized vortices and vortex lattices can be considered by applying the approach, developed for the BEC of cold atoms in [
61].
The Josephson phenomena for two trapped condensates of dipolar excitons can be studied by applying an approach similar to the one developed for non-Abelian Josephson effect between two
spinor Bose-Einstein condensates of cold atoms in double optical traps [
62].