Next Article in Journal
Zinc-Ion Storage Mechanism of Polyaniline for Rechargeable Aqueous Zinc-Ion Batteries
Next Article in Special Issue
Quantum Diffusion in the Lowest Landau Level of Disordered Graphene
Previous Article in Journal
Perspectives on Iron Oxide-Based Materials with Carbon as Anodes for Li- and K-Ion Batteries
Previous Article in Special Issue
Simultaneous Extraction of the Grain Size, Single-Crystalline Grain Sheet Resistance, and Grain Boundary Resistivity of Polycrystalline Monolayer Graphene
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Superfluidity of Dipolar Excitons in a Double Layer of αT3 with a Mass Term

1
Physics Department, New York City College of Technology, City University of New York, New York, NY 11201, USA
2
The Graduate School and University Center, City University of New York, New York, NY 10016, USA
3
Department of Physics and Astronomy, Hunter College, City University of New York, New York, NY 10065, USA
4
Donastia International Physics Center (DIPC), P de Manuel Lardizabal, 4, 20018 San Sebastian, Spain
5
US Military Academy at West Point, 606 Thayer Road, West Point, NY 10996, USA
*
Author to whom correspondence should be addressed.
Nanomaterials 2022, 12(9), 1437; https://doi.org/10.3390/nano12091437
Submission received: 13 January 2022 / Revised: 30 March 2022 / Accepted: 12 April 2022 / Published: 22 April 2022
(This article belongs to the Special Issue Graphene for Electronics)

Abstract

:
We predict Bose-Einstein condensation and superfluidity of dipolar excitons, formed by electron-hole pairs in spatially separated gapped hexagonal α T 3 (GHAT3) layers. In the α T 3 model, the AB-honeycomb lattice structure is supplemented with C atoms located at the centers of the hexagons in the lattice. We considered the α T 3 model in the presence of a mass term which opens a gap in the energy-dispersive spectrum. The gap opening mass term, caused by a weak magnetic field, plays the role of Zeeman splitting at low magnetic fields for this pseudospin-1 system. The band structure of GHAT3 monolayers leads to the formation of two distinct types of excitons in the GHAT3 double layer. We consider two types of dipolar excitons in double-layer GHAT3: (a) “A excitons”, which are bound states of electrons in the conduction band (CB) and holes in the intermediate band (IB), and (b) “B excitons”, which are bound states of electrons in the CB and holes in the valence band (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 T c for superfluidity.

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 α T 3 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 α T 3 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 ( α = 1 ) [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 0 α = tan φ 1 . Optically induced dressed states [38], and their tunneling, transport [33,39], and collective properties [40], as well as α T 3 based nanoribbons [41] have been analyzed. The BEC and superfluidity of dipolar magnetoexcitons in α T 3 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 T c 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 α T 3 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 α T 3 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]
H ^ λ = Δ f ( k ) cos ϕ 0 f * ( k ) cos ϕ 0 f ( k ) sin ϕ 0 f * ( k ) sin ϕ Δ ,
where the origin in k -space is defined to be around the K point, k = ( k x , k y ) and tan θ k = k y / k x , ϕ = tan 1 α , f ( k ) = ħ v F λ k x i k y = λ ħ v F k e i λ θ k , with λ = ± 1 being the valley index at the K and K’ points, 2 Δ is the gap in the energy spectrum of a GHAT3 layer due to the mass term in the Hamiltonian. In an α T 3 layer honeycomb lattice, there is an added fermionic hub atom C at the center of each hexagon. Let the hopping integral be t 1 between the hub atom and either an A or B atom on the rim and t 2 between nearest neighbors on the rim of the hexagon. The ratio of these two nearest neighbor hopping terms is denoted as t 2 / t 1 = α , where the parameter α satisfies 0 α 1 . 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 ϵ C B ( k ) is given by the relation [43]
ϵ C B ( k ) Δ + ħ 2 k 2 2 m C B ,
where k = p / ħ and p are the wave vector and momentum of a quasiparticle, m C B is the effective mass of the charge carriers in the conduction band, given by
m C B = 1 + α 2 Δ 2 v F 2 ,
where v F is the Fermi velocity in a GHAT3 layer, and φ = tan 1 α [43]. At small momenta near K and K’ points, the dispersion for the charge carriers in the valence band ϵ V B ( k ) is given by the relation [43]
ϵ V B ( k ) Δ ħ 2 k 2 2 m V B ,
with m V B the effective mass of the charge carriers in the valence band, given by
m V B = 1 + α 2 Δ 2 v F 2 α 2 .
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 α T 3 layer without a mass term, ϵ I B ( k ) is given by the relation [43]
ϵ I B ( k ) ħ 2 k 2 2 m I B ,
where m I B is the effective mass of the charge carriers in the intermediate band, given by
m I B = 1 + α 2 Δ 2 v F 2 1 α 2 .
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 ϵ d 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 V ( r e h ) , where r e h 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 r e h by its projection r on this plane, the potential V ( r e h ) may be expressed as V ( r e h ) = V ( r 2 + D 2 ) , 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 r 1 = ( x 1 , y 1 ) and r 2 = ( x 2 , y 2 ) for the electron and the projection vector of the hole, respectively (where r = r 1 r 2 ), the dipolar exciton can be described by employing a two-body 2D Schrödinger equation with potential V ( r 2 + D 2 ) . 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 α T 3 . 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
H ^ e x = ħ 2 2 m e Δ r 1 ħ 2 2 m h Δ r 2 + V ( r ) ,
where Δ r 1 and Δ r 2 are the Laplacian operators with respect to the components of the vectors r 1 and r 2 , respectively, and m e and m h are the effective masses of the electron and hole, respectively. For CV excitons m e = m C B and m h = m V B ; and for CI excitons m e = m C B and m h = m I B , where m C B , m V B , and m I B 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 μ = m e m h / m e + m h in a V ( r ) potential and motion of the center-of-mass of the exciton with the mass M = m e + m h . We introduce the coordinates of the center-of-mass R of an exciton and the coordinate of the relative motion r of an electron and hole as R = m e r 1 + m h r 2 / m e + m h and r = r 1 r 2 , correspondingly. The Hamiltonian H ^ e x can be represented in the form: H ^ e x = H ^ R + H ^ r , where the Hamiltonian of the motion of the center-of-mass is H ^ R and that of the relative motion of electron and a hole is H ^ r . The solution of the Schrödinger equation for the center-of-mass of an exciton H ^ R ψ ( R ) = E ψ ( R ) is the plane wave ψ ( R ) = e i P · R / ħ with the quadratic energy spectrum E = P 2 / ( 2 M ) , where P 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 ϵ d . The potential energy of electron-hole Coulomb attraction is
V ( r ) = κ e 2 ϵ d r 2 + D 2 ,
where κ = 9 × 10 9 N × m 2 / C 2 , ϵ d is the dielectric constant of the insulator (SiO 2 or h-BN), surrounding the electron and hole GHAT3 monolayers, forming the double layer. For the h-BN barrier we substitute the dielectric constant ϵ d = 4.89 , while for the SiO 2 barrier we substitute the dielectric constant ϵ d = 4.50 . For h-BN insulating layers, ϵ d = 4.89 is the effective dielectric constant, defined as ϵ d = ε ε [14], where ε = 6.71 and ε = 3.56 are the components of the dielectric tensor for h-BN [48]. Assuming r D , we approximate V ( r ) by the first two terms of the Taylor series and obtain
V ( r ) = V 0 + γ r 2 , where V 0 = κ e 2 ϵ d D , γ = κ e 2 2 ϵ d D 3 .
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 H ^ r Ψ ( r ) = E Ψ ( r ) 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
Ψ N L ( r ) = N ! a | L | + 1 n ˜ ! n ˜ ! 2 | L | / 2 sgn ( L ) L r | L | e r 2 / ( 4 a 2 ) × L N | L | ( r 2 / ( 2 a 2 ) ) e i L φ ( 2 π ) 1 / 2 ,
where N = min ( n ˜ , n ˜ ) , L = n ˜ n ˜ , n ˜ , n ˜ = 0 , 1 , 2 , 3 , are the quantum numbers, φ is the polar angle, and a = ħ / 2 2 μ γ 1 / 2 is a Bohr radius of a dipolar exciton. The corresponding energy spectrum is given by
E N L E e ( h ) = V 0 + ( 2 N + 1 + | L | ) ħ 2 γ μ 1 / 2 .
At the lowest quantum state N = L = 0 as it follows from Equation (12) the ground state energy for the exciton is given by
E 00 = V 0 + ħ 2 γ μ 1 / 2 .
The important characteristic of the exciton is the square of the in-plane gyration radius r X 2 . 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 r X 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
r X 2 r 2 = Ψ 00 * ( r ) r 2 Ψ 00 ( r ) d 2 r = 2 π 2 π a 2 0 + r 2 e r 2 2 a 2 r d r = 2 a 2 .
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
μ C V = m C B m V B m C B + m V B = Δ 2 v F 2 ; M C V = m C B + m V B = 1 + α 2 2 Δ 2 v F 2 α 2 .
For CI excitons one has
μ C I = m C B m I B m C B + m I B = 1 + α 2 Δ 2 v F 2 2 α 2 ; M C I = m C B + m I B = 1 + α 2 2 α 2 Δ 2 v F 2 1 α 2 .

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 n A a B A 2 1 and n B a B B 2 1 , where n A ( B ) and a B A ( B ) 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 s = 16 below. The Hamiltonian H ^ of the 2D A and B weakly interacting dipolar excitons can be written as
H ^ = H ^ A + H ^ B + H ^ I ,
where H ^ A ( B ) are the Hamiltonians of A(B) excitons defined as
H ^ A ( B ) = k E A ( B ) ( k ) a k A ( B ) a k A ( B ) + g A A ( B B ) 2 S k l m a k A ( B ) a l A ( B ) a A ( B ) m a A ( B ) k + l m ,
and H ^ I is the Hamiltonian of the interaction between A and B excitons presented as
H ^ I = g A B S k l m a k A a l B a B m a A k + l m ,
where a k A ( B ) and a k A ( B ) are Bose creation and annihilation operators for A(B) dipolar excitons with the wave vector k , correspondingly, S is the area of the system, E A ( B ) ( k ) ϵ A ( B ) = ε ( 0 ) A ( B ) ( k ) + A A ( B ) is the energy spectrum of non-interacting A(B) dipolar excitons, respectively, ε ( 0 ) A ( B ) ( k ) = ħ 2 k 2 / ( 2 M A ( B ) ) , M A ( B ) is an effective mass of non-interacting dipolar excitons, A A ( B ) is the constant, which depends on A(B) dipolar exciton binding energy and the corresponding gap, g A A ( B B ) and g A B 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 U ( R ) = κ e 2 D 2 / ( ϵ d R 3 ) . 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,
G A A = g A A n A = g n A , G B B = g B B n B = g n B , G A B = g A B n A n B = g n A n B , ω A ( k ) = ε ( 0 ) A 2 ( k ) + 2 G A A ε ( 0 ) A ( k ) , ω B ( k ) = ε ( 0 ) B 2 ( k ) + 2 G B B ε ( 0 ) B ( k ) ,
one obtains two modes of the spectrum of Bose collective excitations ε j ( k )
ε j ( k ) = ω A 2 ( k ) + ω B 2 ( k ) + ( 1 ) j 1 ω A 2 ( k ) ω B 2 ( k ) 2 + 4 G A B 2 ε ( 0 ) A ( k ) ε ( 0 ) B ( k ) 2 ,
where j = 1 , 2. In our approach, the condition G A B 2 = G A A G B B holds.
At small momenta p = ħ k , when ε ( 0 ) A ( k ) G A A and ε ( 0 ) B ( k ) G B B , expanding the spectrum of collective excitations ε j ( k ) up to the first order with respect to the momentum p, one obtains two sound modes in the spectrum of the collective excitations ε j ( p ) = c j p , where c j is the sound velocity written as
c j = G A A 2 M A + G B B 2 M B + ( 1 ) j 1 G A A 2 M A G B B 2 M B 2 + G A B 2 M A M B ,
At j = 1 , the spectrum of collective excitations is determined by the non-zero sound velocity c 1 , while at j = 2 the sound velocity vanishes with c 2 = 0 . At large momenta, for the conditions when ε ( 0 ) A ( k ) G A A and ε ( 0 ) B ( k ) G B B , one obtains two parabolic modes of collective excitations with the spectra ε 1 ( k ) = ε ( 0 ) A ( k ) and ε 2 ( k ) = ε ( 0 ) B ( k ) , if M A < M B and if M A > M B with the spectra ε 1 ( k ) = ε ( 0 ) B ( k ) and ε 2 ( k ) = ε ( 0 ) A ( k ) .

4. Superfluidity of the Weakly-Interacting Bose Gas of Dipolar Excitons

Since when j = 2 the sound velocity vanishes, below we take into account only the branch of the spectrum of collective excitations at j = 1 , neglecting the branch at j = 2 . 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 j = 1 is sound-like with the finite sound velocity, c 1 . 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 v c = c 1 . 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 j = 1 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 ρ s ( T ) as ρ s ( T ) = ρ ρ n ( T ) , where ρ = M A n A + M B n B is the total 2D density of the dipolar excitons and ρ n ( T ) denotes the density of the normal component. The density ρ n ( T ) of the normal component can be defined using standard procedure [56]. The assumption that the dipolar exciton system moves with a velocity u implies that the superfluid component moves with the velocity u . 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
J = s d 2 p ( 2 π ħ ) 2 p f ε 1 ( p ) p · u ,
where s = 16 is the spin and valley degeneracy factor, f ε 1 ( p ) ) = exp ε 1 ( p ) / ( k B T ) 1 1 is the Bose-Einstein distribution function for the quasiparticles with the dispersion ε 1 ( p ) , and k B is the Boltzmann constant. Expanding the expression under the integral in Equation (23) up to the first order with respect to p · u / ( k B T ) , one has:
J = s u 2 d 2 p ( 2 π ħ ) 2 p 2 f ε 1 ( p ) ε 1 ( p ) .
The density ρ n of the normal component in the moving weakly-interacting Bose gas of dipolar excitons is defined as [56]
J = ρ n u .
Employing Equations (24) and (25), one derives the normal component density as
ρ n ( T ) = s 2 d 2 p ( 2 π ħ ) 2 p 2 f ε 1 ( p ) ε 1 ( p ) .
At low temperatures k B T M A ( B ) c j 2 , the small momenta ( ε ( 0 ) A ( k ) G A A and ε ( 0 ) B ( k ) G B B ) 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 ε 1 ( k ) = c 1 k with the sound velocity defined by Equation (22). By substituting ε 1 ( k ) = c 1 k into Equation (26), we obtain
ρ n ( T ) = 3 s ζ ( 3 ) 2 π ħ 2 c 1 4 k B 3 T 3 ,
where ζ ( z ) is the Riemann zeta function ( ζ ( 3 ) 1.202 ).
The mean field critical temperature T c of the phase transition at which the superfluidity occurs, implying neglecting the interaction between the quasiparticles, is obtained from the condition ρ s ( T c ) = 0 [56]:
ρ n ( T c ) = ρ = M A n A + M B n B .
At low temperatures k B T M A ( B ) c 1 2 by substituting Equation (27) into Equation (28), one derives
T c = 2 π ħ 2 ρ c 1 4 3 ζ ( 3 ) s k B 3 1 / 3 .
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 T c 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 E b ( α , Δ , D ) for CV and CI excitons as functions of the gap Δ for chosen parameter α = 0.6 and interlayer separations D = 25 nm . According to Figure 2, E b ( α , Δ , D ) 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 E b ( α , Δ , D ) for CV and CI excitons as functions of the parameter α for chosen gap Δ = 0.5 ħ v F / a and interlayer separations D = 25 nm . According to Figure 3, E b ( α , Δ , D ) does not depend on α for CV excitons, whereas it is an increasing function of α for CI excitons. At α 0.7 E b ( α , Δ , D ) for CV excitons is larger than for CI excitons, while at α 0.7 E b ( α , Δ , D ) for CI excitons is larger than for CV excitons.
In Figure 4, we present the results of our calculations for the exciton binding energy E b ( α , Δ , D ) for CV and CI excitons as functions of the interlayer separation D for chosen parameter α = 0.6 and gap Δ = 0.5 ħ v F / a . According to Figure 4, E b ( α , Δ , D ) 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 α = 0.6 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 Δ = 0.5 ħ v F / a 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 α 0.7 , both M and μ for CV excitons are larger than for CI excitons, but when α 0.7 E b ( α , Δ , D ) both M and μ for CV excitons are smaller than for CI excitons.
Figure 7 demonstrates the dependence of the sound velocity c c 1 on the hopping parameter α for chosen Δ = ħ v F / a , interlayer separations D = 25 nm at fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively. According to Figure 7, c does not depend much on α when α 0.5 , while for α 0.5 , the sound velocity c is a decreasing function of α .
In Figure 8, we plot the sound velocity c c 1 versus the gap Δ for chosen parameter α = 0.6 , interlayer separations D = 25 nm for chosen concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 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 c c 1 as a function of the interlayer separation D for hopping parameter α = 0.6 and gap Δ = 0.5 ħ v F / a , for fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 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 c c 1 on the concentrations n A and n B of A and B excitons, respectively for chosen hopping parameter α = 0.6 and gap Δ = 0.5 ħ v F / a , at fixed interlayer separation D = 25 nm . According to Figure 10, the sound velocity c is an increasing function of both concentrations n A and n B .
In Figure 11, we present the mean-field phase transition critical temperature T c ( n A , n B , α , Δ , D ) as a function of the parameter α for chosen gap Δ = 0.5 ħ v F / a , interlayer separations D = 25 nm at the fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively. According to Figure 11, T c is a decreasing function of α at α 0.9 , while at α 0.9 the critical temperature T c is an increasing function of α .
In Figure 12, we present the mean-field phase transition critical temperature T c ( n A , n B , α , Δ , D ) as a function of the gap Δ for chosen parameter α = 0.6 , interlayer separations D = 25 nm at the fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively. According to Figure 12, the criticaltemperature T c is a decreasing function of Δ .
In Figure 13, we demonstrate the mean-field phase transition critical temperature T c ( n A , n B , α , Δ , D ) as a function of the interlayer separation D for chosen parameter α = 0.6 and gap Δ = 0.5 ħ v F / a , at the fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively. According to Figure 13, the critical temperature T c is an increasing function of D.
In Figure 14, we present density plots for the mean-field phase transition critical temperature T c ( n A , n B , α , Δ , D ) as a function of the concentrations n A and n B of A and B excitons, respectively for chosen parameter α = 0.6 and gap Δ = 0.5 ħ v F / a , at the fixed interlayer separation D = 25 nm . According to Figure 14, the the critical temperature T c is an increasing function of both the concentrations n A and n B .
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 T c is a decreasing function of the gap Δ . Therefore, there should be an optimal value for Δ , which would correspond to relatively high T c 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 F = 2 spinor Bose-Einstein condensates of cold atoms in double optical traps [62].

6. Conclusions

This paper is devoted to an investigation of the existence of BEC and superfluidity of dipolar excitons in double layers of GHAT3 which was proposed and analyzed. We have derived the solution of a two-body problem for an electron and a hole for the model Hamiltonian representing double-layer GHAT3. We predict the formation of two types of dipolar excitons, characterized by different binding energies and effective masses, in the double layer of GHAT3. We have calculated the binding energy, effective mass, spectrum of collective excitations, superfluid density, and the mean-field critical temperature of the phase transition to the superfluid state for the two-component weakly interacting Bose gas of A and B dipolar excitons in double-layer GHAT3. We have demonstrated that at fixed exciton density, the mean-field critical temperature for superfluidity of dipolar excitons is decreased as a function of the gap Δ . Our results show that T c is increased as a function of the density n and is decreased as a function of the gap Δ and the interlayer separation D.
The occupancy of the superfluid state at T < T c can result in the existence of persistent dissipationless superconducting oppositely directed electric currents in each GHAT3 layer, forming a double layer. According to the presented results of our calculations, while the external weak magnetic field, responsible for the formation of the gap Δ in the double layer of α T 3 increases the exciton binding energy, the mean-field transition temperature to the superfluid phase is increased as the weak magnetic field and Δ are decreased. Therefore, the dipolar exciton system in a double-layer of GHAT3 can be applied to engineer a switch, where transport properties of dipolar excitons can be tuned by an external weak magnetic field, forming the gap Δ . Varying a weak magnetic field may lead to a phase transition between the superfluid and normal phase, which sufficiently changes the transport properties of dipolar excitons.

Author Contributions

Conceptualization, O.L.B. and G.G.; computations, G.P.M. and P.F.; writing reviewing and editing, O.L.B. and G.G.; supervision, O.L.B., G.G. and P.F.; project administration, G.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by U.S. ARO grant No. W911NF1810433. G.G. would like to acknowledge the support from the Air Force Research Laboratory (AFRL) through Grant No. FA9453-21-1-0046.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lozovik, Y.E.; Yudson, V.I. Feasibility of superfluidity of paired spatially separated electrons and holes; a new superconductivity mechanism. Sov. Phys. JETP Lett. 1975, 22, 274. [Google Scholar]
  2. Lozovik, Y.E.; Yudson, V.I. A new mechanism for superconductivity: Pairing between spatially separated electrons and holes. Sov. Phys. JETP Lett. 1976, 44, 389. [Google Scholar]
  3. Snoke, D.W. Spontaneous Bose coherence of excitons and polaritons. Science 2002, 298, 1368. [Google Scholar] [CrossRef] [PubMed]
  4. Butov, L.V. Condensation and pattern formation in cold exciton gases in coupled quantum wellsm. J. Phys. Condens. Matter 2004, 16, R1577. [Google Scholar] [CrossRef]
  5. Eisenstein, J.P.; MacDonald, A.H. Bose–Einstein condensation of excitons in bilayer electron systems. Nature 2004, 432, 691. [Google Scholar] [CrossRef] [Green Version]
  6. Snoke, D.W. Quantum Gases: Finite Temperature and Non-equilibrium Dynamics; Cold Atom Series; Proukakis, N.P., Gardiner, S.A., Davis, M.J., Szymanska, M.H., Eds.; Imperial College Press: London, UK, 2013; Volume 1, p. 419. [Google Scholar]
  7. Saberi-Pouya, S.; Conti, S.; Perali, A.; Croxall, A.F.; Hamilton, A.R.; Peeters, F.M.; Neilson, D. Experimental conditions for the observation of electron-hole superfluidity in GaAs heterostructures. Phys. Rev. B 2020, 101, 140501(R). [Google Scholar] [CrossRef] [Green Version]
  8. Berman, O.L.; Lozovik, Y.E.; Gumbs, G. Dynamical equation for an electron-hole pair condensate in a system of two graphene layers. Phys. Rev. B 2008, 77, 155433. [Google Scholar] [CrossRef] [Green Version]
  9. Lozovik, Y.E.; Sokolik, A.A. Electron-hole pair condensation in a graphene bilayer. JETP Lett. 2008, 87, 55. [Google Scholar] [CrossRef]
  10. Lozovik, Y.E.; Sokolik, A.A. Multi-band pairing of ultrarelativistic electrons and holes in graphene bilayer. Phys. Lett. A 2009, 374, 326. [Google Scholar] [CrossRef]
  11. Bistritzer, R.; MacDonald, A.H. Influence of disorder on electron-hole pair condensation in graphene bilayers. Phys. Rev. Lett. 2008, 101, 256406. [Google Scholar] [CrossRef] [Green Version]
  12. Berman, O.L.; Kezerashvili, R.Y.; Ziegler, K. Superfluidity of dipole excitons in the presence of band gaps in two-layer graphene. Phys. Rev. B 2012, 85, 035418. [Google Scholar] [CrossRef] [Green Version]
  13. Perali, A.; Neilson, D.; Hamilton, A.R. High-temperature superfluidity in double-bilayer graphene. Phys. Rev. Lett. 2013, 110, 146803. [Google Scholar] [CrossRef]
  14. Fogler, M.M.; Butov, L.V.; Novoselov, K.S. High-temperature superfluidity with indirect excitons in van der Waals heterostructures. Nat. Commun. 2014, 5, 4555. [Google Scholar] [CrossRef] [Green Version]
  15. Wu, F.-C.; Xue, F.; MacDonald, A.H. Theory of two-dimensional spatially indirect equilibrium exciton condensates. Phys. Rev. B 2015, 92, 165121. [Google Scholar] [CrossRef] [Green Version]
  16. Berman, O.L.; Kezerashvili, R.Y. High-temperature superfluidity of the two-component Bose gas in a transition metal dichalcogenide bilayer. Phys. Rev. B 2016, 93, 245410. [Google Scholar] [CrossRef] [Green Version]
  17. Berman, O.L.; Kezerashvili, R.Y. Superfluidity of dipolar excitons in a transition metal dichalcogenide double layer. Phys. Rev. B 2017, 96, 094502. [Google Scholar] [CrossRef] [Green Version]
  18. Conti, S.; Van der Donck, M.; Perali, A.; Peeters, F.M.; Neilson, D. Doping-dependent switch from one- to two-component superfluidity in coupled electron-hole van der Waals heterostructures. Phys. Rev. B 2020, 101, 220504(R). [Google Scholar] [CrossRef]
  19. Berman, O.L.; Gumbs, G.; Kezerashvili, R.Y. Bose-Einstein condensation and superfluidity of dipolar excitons in a phosphorene double layer. Phys. Rev. B 2017, 96, 014505. [Google Scholar] [CrossRef] [Green Version]
  20. Saberi-Pouya, S.; Zarenia, M.; Perali, A.; Vazifehshenas, T.; Peeters, F.M. High-temperature electron-hole superfluidity with strong anisotropic gaps in double phosphorene monolayers. Phys. Rev. B 2018, 97, 174503. [Google Scholar] [CrossRef] [Green Version]
  21. Mazuz-Harpaz, Y.; Cohen, K.; Leveson, M.; West, K.; Pfeiffer, L.; Khodas, M.; Rapaport, R. Dynamical formation of a strongly correlated dark condensate of dipolar excitons. Proc. Nat. Acad. Sci. USA 2019, 116, 18328. [Google Scholar] [CrossRef] [Green Version]
  22. Raoux, A.; Morigi, M.; Fuchs, J.-N.; Piéchon, F.; Montambaux, G. From dia- to paramagnetic orbital susceptibility of massless Fermions. Phys. Rev. Lett. 2014, 112, 026402. [Google Scholar] [CrossRef] [Green Version]
  23. Sutherland, B. Localization of electronic wave functions due to local topology. Phys. Rev. B 1986, 34, 5208. [Google Scholar] [CrossRef]
  24. Illes, E.; Carbotte, J.P.; Nicol, E.J. Hall quantization and optical conductivity evolution with variable Berry phase in the αT3 model. Phys. Rev. B 2015, 92, 245410. [Google Scholar] [CrossRef] [Green Version]
  25. Islam, S.K.F.; Dutta, P. Valley-polarized magnetoconductivity and particle-hole symmetry breaking in a periodically modulated αT3 lattice. Phys. Rev. B 2017, 96, 045418. [Google Scholar] [CrossRef] [Green Version]
  26. Illes, E.; Nicol, E.J. Magnetic properties of the αT3 model: Magneto-optical conductivity and the Hofstadter butterfly. Phys. Rev. B 2016, 94, 125435. [Google Scholar] [CrossRef] [Green Version]
  27. Dey, B.; Ghosh, T.K. Photoinduced valley and electron-hole symmetry breaking in αT3 lattice: The role of a variable Berry phase. Phys. Rev. B 2018, 98, 075422. [Google Scholar] [CrossRef] [Green Version]
  28. Dey, B.; Ghosh, T.K. Floquet topological phase transition in the αT3 lattice. Phys. Rev. B 2019, 99, 205429. [Google Scholar] [CrossRef] [Green Version]
  29. Biswas, T.; Ghosh, T.K. Dynamics of a quasiparticle in the αT3 model: Role of pseudospin polarization and transverse magnetic field on zitterbewegung. J. Phys. Condens. Matter 2018, 30, 075301. [Google Scholar] [CrossRef] [Green Version]
  30. Kovacs, A.D.; David, G.; Dora, B.; Cserti, J. Frequency-dependent magneto-optical conductivity in the generalized αT3 model. Phys. Rev. B 2017, 95, 035414. [Google Scholar] [CrossRef] [Green Version]
  31. Biswas, T.; Ghosh, T.K. Magnetotransport properties of the αT3 model. J. Phys. Condens. Matter 2016, 28, 495302. [Google Scholar] [CrossRef] [Green Version]
  32. Oriekhov, D.O.; Gusynin, V.P. RKKY interaction in a doped pseudospin-1 fermion system at finite temperature. arXiv 2020, arXiv:2001.00272. [Google Scholar] [CrossRef]
  33. Huang, D.; Iurov, A.; Xu, H.-Y.; Lai, Y.-C.; Gumbs, G. Interplay of Lorentz-Berry forces in position-momentum spaces for valley-dependent impurity scattering in αT3 lattices. Phys. Rev. B 2019, 99, 245412. [Google Scholar] [CrossRef] [Green Version]
  34. Li, Y.; Kita, S.; Munoz, P.; Reshef, O.; Vulis, D.I.; Yin, M.; Loncar, M.; Mazur, E. On-chip zero-index metamaterials. Nat. Photon 2015, 9, 738. [Google Scholar] [CrossRef]
  35. Xu, H.-Y.; Huang, L.; Huang, D.H.; Lai, Y.-C. Geometric valley Hall effect and valley filtering through a singular Berry flux. Phys. Rev. B 2017, 96, 045412. [Google Scholar] [CrossRef] [Green Version]
  36. Leykam, D.; Andreanov, A.; Flach, S. Artificial flat band systems: From lattice models to experiments. Adv. Phys. X 2018, 3, 677. [Google Scholar] [CrossRef]
  37. Sherafati, M.; Satpathy, S. Analytical expression for the RKKY interaction in doped graphene. Phys. Rev. B 2011, 84, 125416. [Google Scholar] [CrossRef] [Green Version]
  38. Iurov, A.; Gumbs, G.; Huang, D. Peculiar electronic states, symmetries, and Berry phases in irradiated αT3 materials. Phys. Rev. B 2019, 99, 205135. [Google Scholar] [CrossRef] [Green Version]
  39. Iurov, A.; Zhemchuzhna, L.; Dahal, D.; Gumbs, G.; Huang, D. Quantum-statistical theory for laser-tuned transport and optical conductivities of dressed electrons in αT3 materials. Phys. Rev. B 2020, 101, 035129. [Google Scholar] [CrossRef] [Green Version]
  40. Weekes, N.; Iurov, A.; Zhemchuzhna, L.; Gumbs, G.; Huang, D. Generalized WKB theory for electron tunneling in gapped αT3 lattices. Phys. Rev. B 2021, 103, 165429. [Google Scholar] [CrossRef]
  41. Iurov, A.; Zhemchuzhna, L.; Gumbs, G.; Huang, D.; Fekete, P.; Anwar, F.; Dahal, D.; Weekes, N. Tailoring plasmon excitations in αT3 armchair nanoribbons. Sci. Rep. 2021, 11, 20577. [Google Scholar] [CrossRef]
  42. Abranyos, Y.; Berman, O.L.; Gumbs, G. Superfluidity of dipolar excitons in doped double-layered hexagonal lattice in a strong magnetic field. Phys. Rev. B 2020, 102, 155408. [Google Scholar] [CrossRef]
  43. Balassis, A.; Gumbs, G.; Roslyak, O. Temperature-Induced Plasmon Excitations for the αT3 Lattice in Perpendicular Magnetic Field. Nanomaterials 2021, 11, 1720. [Google Scholar] [CrossRef]
  44. Balassis, A.; Dahal, D.; Gumbs, G.; Iurov, A.; Huang, D.; Roslyak, O. Magnetoplasmons for the αT3 model with filled Landau levels. J. Phys. Condens. Matter 2020, 32, 485301. [Google Scholar] [CrossRef]
  45. Malcolm, J.D.; Nicol, E.J. Frequency-dependent polarizability, plasmons, and screening in the two-dimensional pseudospin-1 dice lattice. Phys. Rev. B 2016, 93, 165433. [Google Scholar] [CrossRef] [Green Version]
  46. Illes, E. Properties of the αT3 Model. Ph.D. Thesis, University of Guelph, Guelph, ON, Canada, 2017. [Google Scholar]
  47. Wunsch, B.; Stauber, T.; Sols, F.; Guinea, F. Dynamical polarization of graphene at finite doping. New J. Phys. 2006, 8, 318. [Google Scholar] [CrossRef]
  48. Cai, Y.; Zhang, L.; Zeng, Q.; Cheng, L.; Xu, Y. Infrared reflectance spectrum of BN calculated from first principles. Solid State Commun. 2007, 141, 262. [Google Scholar] [CrossRef] [Green Version]
  49. Maksym, P.A.; Chakraborty, T. Quantum dots in a magnetic field: Role of electron-electron interactions. Phys. Rev. Lett. 1990, 65, 108. [Google Scholar] [CrossRef]
  50. Iyengar, A.; Wang, J.; Fertig, H.A.; Brey, L. Excitations from filled Landau levels in graphene. Phys. Rev. B 2007, 75, 125430. [Google Scholar] [CrossRef] [Green Version]
  51. Berman, O.L.; Kezerashvili, R.Y.; Kolmakov, G.V.; Lozovik, Y.E. Turbulence in a Bose-Einstein condensate of dipolar excitons in coupled quantum wells. Phys. Rev. B 2012, 86, 045108. [Google Scholar] [CrossRef] [Green Version]
  52. Lifshitz, E.M.; Pitaevskii, L.P. Statistical Physics, Part 2; Pergamon Press: Oxford, UK, 1980. [Google Scholar]
  53. Tommasini, P.; de Passos, E.J.V.; de Toledo Piza, A.F.R.; Hussein, M.S.; Timmermans, E. Bogoliubov theory for mutually coherent condensates. Phys. Rev. A 2003, 67, 023606. [Google Scholar] [CrossRef]
  54. Sun, B.; Pindzola, M.S. Bogoliubov modes and the static structure factor for a two-species Bose–Einstein condensate. J. Phys. B 2010, 43, 055301. [Google Scholar] [CrossRef]
  55. Abrikosov, A.A.; Gorkov, L.P.; Dzyaloshinskii, I.E. Methods of Quantum Field Theory in Statistical Physics; Prentice-Hall: Englewood Cliffs, NJ, USA, 1963. [Google Scholar]
  56. Pitaevskii, L.; Stringari, S. Bose-Einstein Condensation; Clarendon Press: Oxford, UK, 2003. [Google Scholar]
  57. López Ríos, P.; Perali, A.; Needs, R.J.; Neilson, D. Evidence from quantum Monte Carlo simulations of large-gap superfluidity and BCS-BEC crossover in double electron-hole layers. Phys. Rev. Lett. 2018, 120, 177701. [Google Scholar] [CrossRef] [Green Version]
  58. Kosterlitz, J.M.; Thouless, D.J. Ordering, metastability and phase transitions in two-dimensional systems. J. Phys. C 1973, 6, 1181. [Google Scholar] [CrossRef]
  59. Nelson, D.R.; Kosterlitz, J.M. Universal jump in the superfluid density of two-dimensional superfluids. Phys. Rev. Lett. 1977, 39, 1201. [Google Scholar] [CrossRef] [Green Version]
  60. Vyas, V.M.; Panigrahia, P.K.; Banerji, J. A scheme to observe universal breathing mode and Berezinskii–Kosterlitz–Thouless phase transition in a two-dimensional photon gas. Phys. Lett. A 2014, 378, 1434. [Google Scholar] [CrossRef]
  61. Ji, A.-C.; Liu, W.M.; Song, J.L.; Zhou, F. Dynamical creation of fractionalized vortices and vortex lattices. Phys. Rev. Lett. 2008, 101, 010402. [Google Scholar] [CrossRef] [Green Version]
  62. Qi, R.; Yu, X.-L.; Li, Z.B.; Liu, W.M. Non-Abelian Josephson effect between two spinor Bose-Einstein condensates in double optical traps. Phys. Rev. Lett. 2009, 102, 185301. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Schematic illustration of a dipolar excitonin a pair of GHAT3 double layers embedded in an insulating material.
Figure 1. Schematic illustration of a dipolar excitonin a pair of GHAT3 double layers embedded in an insulating material.
Nanomaterials 12 01437 g001
Figure 2. The exciton binding energy E b ( α , Δ , D ) for CV and CI excitons as functions of the gap Δ for chosen parameter α = 0.6 and interlayer separations D = 25 nm . The lattice constant of α T 3 is a = 2.46 .
Figure 2. The exciton binding energy E b ( α , Δ , D ) for CV and CI excitons as functions of the gap Δ for chosen parameter α = 0.6 and interlayer separations D = 25 nm . The lattice constant of α T 3 is a = 2.46 .
Nanomaterials 12 01437 g002
Figure 3. The exciton binding energy E b ( α , Δ , D ) for CV and CI excitons as functions of the parameter α for chosen gap Δ = 0.5 ħ v F / a and interlayer separations D = 25 nm .
Figure 3. The exciton binding energy E b ( α , Δ , D ) for CV and CI excitons as functions of the parameter α for chosen gap Δ = 0.5 ħ v F / a and interlayer separations D = 25 nm .
Nanomaterials 12 01437 g003
Figure 4. The exciton binding energy E b ( α , Δ , D ) for CV and CI excitons as functions of the interlayer separation D for chosen parameter α = 0.6 and gap Δ = 0.5 ħ v F / a .
Figure 4. The exciton binding energy E b ( α , Δ , D ) for CV and CI excitons as functions of the interlayer separation D for chosen parameter α = 0.6 and gap Δ = 0.5 ħ v F / a .
Nanomaterials 12 01437 g004
Figure 5. The effective masses of a dipolar exciton for CV and CI excitons as functions of the gap Δ for chosen α = 0.6 for (a) center-of-mass exciton mass M on the left panel and (b) reduced exciton mass μ , on the right.
Figure 5. The effective masses of a dipolar exciton for CV and CI excitons as functions of the gap Δ for chosen α = 0.6 for (a) center-of-mass exciton mass M on the left panel and (b) reduced exciton mass μ , on the right.
Nanomaterials 12 01437 g005
Figure 6. The effective masses of dipolar excitons for CV and CI excitons as functions of the hopping parameter α for chosen gap Δ = 0.5 ħ v F / a for (a) center-of-mass exciton mass M in the left panel and (b) reduced exciton mass μ , on the right.
Figure 6. The effective masses of dipolar excitons for CV and CI excitons as functions of the hopping parameter α for chosen gap Δ = 0.5 ħ v F / a for (a) center-of-mass exciton mass M in the left panel and (b) reduced exciton mass μ , on the right.
Nanomaterials 12 01437 g006
Figure 7. Plot of the sound velocity c c 1 versus α for chosen gap Δ = ħ v F / a , interlayer separations D = 25 nm at fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively.
Figure 7. Plot of the sound velocity c c 1 versus α for chosen gap Δ = ħ v F / a , interlayer separations D = 25 nm at fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively.
Nanomaterials 12 01437 g007
Figure 8. The sound velocity c c 1 versus the gap Δ for chosen parameter α = 0.6, interlayer separations D = 25 nm at the fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively.
Figure 8. The sound velocity c c 1 versus the gap Δ for chosen parameter α = 0.6, interlayer separations D = 25 nm at the fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively.
Nanomaterials 12 01437 g008
Figure 9. The sound velocity c c 1 versus the interlayer separation D for chosen parameter α = 0.6 and gap Δ = 0.5 ħ v F / a , at fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively.
Figure 9. The sound velocity c c 1 versus the interlayer separation D for chosen parameter α = 0.6 and gap Δ = 0.5 ħ v F / a , at fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively.
Nanomaterials 12 01437 g009
Figure 10. The sound velocity c c 1 versus the concentrations n A and n B of A and B excitons, respectively, for chosen parameter α = 0.6 and gap Δ = 0.5 ħ v F / a , at the fixed interlayer separation D = 25 nm .
Figure 10. The sound velocity c c 1 versus the concentrations n A and n B of A and B excitons, respectively, for chosen parameter α = 0.6 and gap Δ = 0.5 ħ v F / a , at the fixed interlayer separation D = 25 nm .
Nanomaterials 12 01437 g010
Figure 11. The mean-field phase transition critical temperature T c ( n A , n B , α , Δ , D ) versus the parameter α for chosen gap Δ = 0.5 ħ v F / a , interlayer separations D = 25 nm at the fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively.
Figure 11. The mean-field phase transition critical temperature T c ( n A , n B , α , Δ , D ) versus the parameter α for chosen gap Δ = 0.5 ħ v F / a , interlayer separations D = 25 nm at the fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively.
Nanomaterials 12 01437 g011
Figure 12. The mean-field phase transition critical temperature T c ( n A , n B , α , Δ , D ) versus the gap Δ for chosen parameter α = 0.6 , interlayer separations D = 25 nm at the fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively.
Figure 12. The mean-field phase transition critical temperature T c ( n A , n B , α , Δ , D ) versus the gap Δ for chosen parameter α = 0.6 , interlayer separations D = 25 nm at the fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively.
Nanomaterials 12 01437 g012
Figure 13. The mean-field phase transition critical temperature T c ( n A , n B , α , Δ , D ) versus the interlayer separation D for chosen parameter α = 0.6 and gap Δ = 0.5 ħ v F / a , at fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively.
Figure 13. The mean-field phase transition critical temperature T c ( n A , n B , α , Δ , D ) versus the interlayer separation D for chosen parameter α = 0.6 and gap Δ = 0.5 ħ v F / a , at fixed concentrations n A = 50 × 10 11 cm 2 and n B = 50 × 10 11 cm 2 of A and B excitons, respectively.
Nanomaterials 12 01437 g013
Figure 14. Density plot for the mean-field phase transition critical temperature T c ( n A , n B , α , Δ , D ) versus the concentrations n A and n B of A and B excitons, respectively, for chosen parameter α = 0.6 and gap Δ = 0.5 ħ v F / a , at the fixed interlayer separation D = 25 nm .
Figure 14. Density plot for the mean-field phase transition critical temperature T c ( n A , n B , α , Δ , D ) versus the concentrations n A and n B of A and B excitons, respectively, for chosen parameter α = 0.6 and gap Δ = 0.5 ħ v F / a , at the fixed interlayer separation D = 25 nm .
Nanomaterials 12 01437 g014
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Berman, O.L.; Gumbs, G.; Martins, G.P.; Fekete, P. Superfluidity of Dipolar Excitons in a Double Layer of αT3 with a Mass Term. Nanomaterials 2022, 12, 1437. https://doi.org/10.3390/nano12091437

AMA Style

Berman OL, Gumbs G, Martins GP, Fekete P. Superfluidity of Dipolar Excitons in a Double Layer of αT3 with a Mass Term. Nanomaterials. 2022; 12(9):1437. https://doi.org/10.3390/nano12091437

Chicago/Turabian Style

Berman, Oleg L., Godfrey Gumbs, Gabriel P. Martins, and Paula Fekete. 2022. "Superfluidity of Dipolar Excitons in a Double Layer of αT3 with a Mass Term" Nanomaterials 12, no. 9: 1437. https://doi.org/10.3390/nano12091437

APA Style

Berman, O. L., Gumbs, G., Martins, G. P., & Fekete, P. (2022). Superfluidity of Dipolar Excitons in a Double Layer of αT3 with a Mass Term. Nanomaterials, 12(9), 1437. https://doi.org/10.3390/nano12091437

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop