Skip to main content

Thank you for visiting You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Evidence for equilibrium exciton condensation in monolayer WTe2


We present evidence that the two-dimensional bulk of monolayer WTe2 contains electrons and holes bound by Coulomb attraction—excitons—that spontaneously form in thermal equilibrium. On cooling from room temperature to 100 K, the conductivity develops a V-shaped dependence on electrostatic doping, while the chemical potential develops a step at the neutral point. These features are much sharper than is possible in an independent-electron picture, but they can be accounted for if electrons and holes interact strongly and are paired in equilibrium. Our calculations from first principles show that the exciton binding energy is larger than 100 meV and the radius as small as 4 nm, explaining their formation at high temperature and doping levels. Below 100 K, more strongly insulating behaviour is seen, suggesting that a charge-ordered state forms. The observed absence of charge density waves in this state is surprising within an excitonic insulator picture, but we show that it can be explained by the symmetries of the exciton wavefunction. Therefore, in addition to being a topological insulator, monolayer WTe2 exhibits strong correlations over a wide temperature range.


In a monolayer semimetal, low carrier densities combined with reduced dimensionality provide the conditions for strong correlation effects. One possible form of such correlations is pairing of electrons and holes in the equilibrium state to form excitons. At low temperatures, such excitons could condense to form an excitonic insulator1,2,3,4. However, exciton formation is expected to be easily disrupted by free carriers (which screen the binding interaction) and thus occur only at low temperatures and near charge neutrality. A number of materials have been mooted as excitonic insulator candidates5,6,7,8,9,10, but there is no consensus as to whether any of them truly contains excitons in equilibrium, either as an incoherent gas or as a coherent condensate, except in the case of bilayer heterostructures at high magnetic fields11,12,13,14.

WTe2 is a layered semimetal, but an exfoliated WTe2 monolayer behaves15,16,17 as a two-dimensional (2D) topological insulator, exhibiting helical conducting edge modes, and becomes superconducting when electrostatically doped18,19. The monolayer has the 1T′ structure shown in Fig. 1a. Its bands are spin degenerate due to inversion symmetry and, near the Fermi energy EF, there is a valence (v) band maximum at Γ flanked by two conduction (c) band minima located at kx = ±kΛ, as sketched in Fig. 1b. Some tunnelling spectroscopy measurements20, angle-resolved photoemission20,21 and density functional theory (DFT) calculations20,21,22,23 point to a positive bandgap, Eg, of the order of 50 meV, while others suggest overlapping bands15,24. Noting that, in the photoemission spectra, the v and c band photoemission features are broad enough that they overlap at least somewhat, we use thick lines in the sketch to signify this uncertainty in Eg.

Fig. 1: Bulk conduction measurements on monolayer WTe2.
figure 1

a, 1T′ structure of monolayer WTe2. The x axis is taken to be along the zigzag W chains. b, Schematic Brillouin zone (above) and bands near EF (below). c, The technique used to exclude edge conduction: when the side contacts are grounded, the measured current between top and bottom contacts must flow through the interior and the edge conduction plateau disappears. (Gp ≡ I/V). These measurements were taken on device MW2 at T = 10 K. d, Gp versus gate-induced density ng at a series of temperatures T on the same device. Inset: temperature dependence at positive values of ng. Bottom: regimes of insulating, metallic and superconducting behaviour identified in previous work.

Source data

This band structure immediately invokes the possibility that excitons could occur in equilibrium in monolayer WTe2, and therefore that the insulating state might not be a simple band insulator25,26,27,28. In this Article we argue that the behaviour of the conductivity and the electron chemical potential, even well above 100 K, is impossible to reconcile with an independent particle picture and strongly indicates the presence of excitons in the equilibrium state. Our first-principles calculations of exciton dispersion and Bohr radius support this conclusion. The insulating behaviour below 100 K suggests a charge-ordered state, but in an excitonic insulator one would normally expect charge density waves, no signs of which are seen in scanning tunnelling microscopy or Raman spectroscopy. To explain this, we show that the entanglement of spin, orbital and valley degrees of freedom hides the charge order, as the contributions to the density wave paired through time reversal cancel out.

The measurements were made on exfoliated monolayer WTe2 flakes with platinum contacts, encapsulated by hexagonal boron nitride (hBN), with graphite gates either below or above (Supplementary Section 1). To study the sheet conductivity while excluding edge conduction, we used the approach illustrated in Fig. 1c. As indicated in the insets, a bias V is applied to one contact and the current I flowing to ground through an opposite contact is measured. When the intervening side contacts are grounded, this current must flow through the bulk. The edge current, which produces the plateau at low Vg, is thereby eliminated and the ‘partial conductance’ Gp ≡ I/V reflects the sheet conductivity σ via \({{G_{\rm{p}}^{ - 1}} \approx {\beta /\sigma} + {R_{\rm{c}}}}\), where β is a geometrical factor considerably larger than one and Rc is the contact resistance. The gate-induced areal number density ng is deduced from the voltages applied to the graphite gate(s) and the geometric capacitances.

Figure 1d shows measurements of Gp versus ng and temperature T. These characteristics are not measurably affected by either a normal displacement field or a normal magnetic field of 14 T, confirming the rejection of edge conduction, which is highly sensitive to magnetic field at low temperatures29. (A recent paper30 reports surprising quantum oscillations at low ng, but we have not seen these in any device.) On cooling from room temperature to 100 K, Gp versus np develops a sharp ‘V’ shape centred close to ng = 0. We have seen consistent behaviour across a dozen monolayer devices, although the sharpness of the V varies, probably as a result of variable sample homogeneity. We also note that a similar sharp V occurs in bilayer WTe2 but at temperatures about five times lower31 (Supplementary Section 7). As shown in the inset, for positive ng smaller than a value nce ≈ +5 × 1012 cm−2, Gp decreases monotonically on cooling, whereas for ng > nce it initially increases. For negative ng (hole doping), a similar but less clearcut transition occurs around ncp ≈ −10 × 1012 cm−2. As indicated in the band at the bottom of Fig. 1d, these values of nce and ncp are consistent with the thresholds for metallic behaviour reported in previous work, where it was found that the metallic state at ng > nce becomes superconducting below ~0.8 K.

Below 100 K, as T decreases Gp collapses over an increasingly wide range of ng. This insulating behaviour, which allows the edge conduction to dominate in normal geometries, is consistent with the findings in ref. 26. We used microwave impedance microscopy32 (MIM; Supplementary Section 2) on devices with no top gate to confirm that this is not a contact effect, as well as to detect any cracks in the monolayer WTe2 that could invalidate the measurements. Figure 2a is a MIM image of device MW10 at 11 K. Red dashed lines mark the edges of the monolayer WTe2 flake, and the wiggly bright lines are cracks. Figure 2b shows the MIM-derived conductivity σMIM, measured in the centre of the white dashed square. Like Gp, it collapses over a range of ng that grows as T falls, indicated by the dotted white contour which is drawn at σMIM ≈ 0.1 μS. The drop-off of Gp at low temperatures, even at large ng, as is evident in Fig. 1d, can be explained by the fact that the monolayer adjacent to the metal contacts is partially screened from the gates and so is less doped and remains insulating. In addition, the contact pattern in MW10 was aligned with the crystal axes, as determined by Raman spectroscopy (Fig. 2a inset and Supplementary Section 3), allowing us to compare conductivities along the x and y axes (Fig. 2c,d and Supplementary Section 4). We see that there is substantial gate-dependent anisotropy, with the lowest conductivity occurring for p-doping parallel to the x axis. This is consistent with the direction in which the valence band edge has a large effective mass.

Fig. 2: Local conductivity imaging and anisotropy.
figure 2

a, MIM image of device MW10. Here, the uncalibrated imaginary MIM signal is plotted for the purpose of identifying conducting features. Multiple conducting cracks are seen running from upper left to lower right. Inset: polar plot of the Raman peak P11 (210 cm−1) intensity, used to determine the crystal axes. b, Conductivity versus ng and T at 2 GHz deduced from MIM measurements in the centre of the white dotted square in a. c, Partial conductance parallel to the x axis (left) and y axis (right), measured in the white dotted square using the configurations shown in the insets.

To measure the chemical potential, we study a device including a separately contacted graphene sheet in parallel with the WTe2 monolayer, as shown in Fig. 3a. Briefly, the WTe2 is almost an equipotential because it has finite conductivity and carries no current. With both the graphene and bottom gate grounded, a voltage Vg is applied to the top gate relative to the WTe2 and the voltage VW on the WTe2 is adjusted to bring the graphene conductance to a minimum. This keeps the graphene neutral and maintains zero electric field beneath the WTe2. The electrostatic potential in the WTe2 is thus in effect fixed to that of the graphene, so the change in VW is due to the change in chemical potential, Δμ = −eΔVW, associated with the gate-induced charge density −eng = εrε0Vg/d, where d is the distance to the gate εr is the hBN dielectric constant and ε0 the vacuum permitivity. From VW versus Vg we thereby obtain μ(ng), choosing the zero of μ at each temperature for convenience. Figure 3b shows measurements of both μ (black) and Gp (red) versus ng made on device MW12. As usual, Gp forms a sharp V as a function of ng at 100 K. Meanwhile, μ exhibits a step at the centre of the V that is still discernible at room temperature and which grows on cooling, saturating at ~40 meV in height below ~50 K. The same behaviour was seen in two devices (Supplementary Section 5).

Fig. 3: Chemical potential measurements and comparison with the single-particle model.
figure 3

a, Schematic of the device structure used to measure the chemical potential versus doping. In parallel with (below) the monolayer WTe2 is a graphene sheet, which is maintained at its Dirac point so that the electric field beneath the WTe2 is zero. The doping ng and chemical potential μ are then obtained from the voltages as shown (see text). The hBN dielectric between the layers is not shown. b, Measurements of μ (black) and conductance Gp (red) versus ng on monolayer WTe2 device MW12. The length of the green bars indicates the thermal energy, kT. c, Single-particle density of states D(E) that reproduces the low-temperature μng behaviour. d, Calculated μ (black line), electron density n (red dotted line) and hole density p (blue dashed line) using this D(E) at the same temperatures as the measurements.

The variations of Gp and μ with ng are impossible to reconcile with a single-particle picture, as can be shown by an elementary calculation: in such a picture, μ and ng are related by \({n_{\rm{g}}} = {{{\int}_{ - \infty }^{ + \infty }} {D\left( E \right)f\left( E \right){\rm{d}}E}}\), where D(E) is the total electron density of states and f(E) = [1 + exp{(E − μ)/kT}]−1. To match the variation of μ with ng at low temperatures, for energies in the relevant range, D(E) must have roughly the form shown in Fig. 3c: a constant value, D0, in the conduction band to give a uniform slope for ng > 0; a gap, Eg; and a large spike at the valence band edge (E = 0). The latter is needed to pin μ to the valence band edge and thus make it independent of ng when ng < 0, to be consistent with the data at 25 K and 10 K. In Fig. 3d we plot the chemical potential calculated at the same temperatures as the measurements in Fig. 3b, using this D(E) with best-fit parameters D0 = 3.7 × 1011 cm−2 meV−1 and Eg = 43 meV. At 150 K the step is washed out; in fact, no choice of D(E) can yield a distinct step in μ whose height is less than kT as is needed to match the measurements at T ≥ 150 K. In particular, states in the gap will only smear the step more. We also plot (dotted lines) the calculated populations of electrons in the conduction band, n, and holes in the valence band, p = ng − n, to contrast their thermally smeared dependence on ng with the sharp V shape seen in the conductance when T ≥ 100 K.

The contradictions between the single-particle picture and the observed dependence of μ and G on ng and T can be largely resolved simply by positing that some electrons and holes are bound as neutral excitons with density nx so that the conductivity is σ=μee(n − nx) +μhe(p − nx), the sum of the contributions of n − nx free electrons and p − nx free holes with respective mobilities μe and μh. In Fig. 4 we compare predictions based on this equation with data from Fig. 3 (Supplementary Section 6 provides more details). When nx = 0 we just have σ=μeen+μhep, which is thermally smeared for any choice of D(E). To illustrate this, in Fig. 4a we plot p and n calculated using the same D(E) shown in Fig. 3c, at 100 K, and in Fig. 4b we plot the calculated σ for nx = 0 (blue dotted), using a mobility ratio chosen to obtain the best match to the measured conductance at 100 K (black line). However, when nx takes its maximal value, determined by the number of minority carriers nx = min(n, p), then for ng > 0 we have σ=μee(n − p) =μeeng while for ng < 0 we have σ=μhe(p − n) =μheng. The result is a sharp, asymmetric V shape (red dashed line) that matches the measurements much better. Note that, in this limit, where only the unbalanced gate-induced charge is free to move, the detailed form of D(E) becomes immaterial.

Fig. 4: Electron–hole pairing and conductivity.
figure 4

a, Total electron (n) and hole (p) densities calculated for T = 100 K using the single-particle density of states D(E) shown in Fig. 3c. The difference is n − p=ng and the shaded area min(n,p) is the maximum possible density nx of non-interacting equilibrium excitons. b, Measured conductance of device MW12 (black solid line) at 100 K compared with the calculated conductance for nx = min(n,p) (red dashed line) and for no excitons, nx = 0 (blue dotted line). c, As in b but at 150 K. The additional green dash-dotted line was obtained by allowing nx to vary in kinetic equilibrium (see text).

As T is increased from 100 K, the conductance at ng = 0 rises and the sides of the V become shallower. This is illustrated in Fig. 4c, where we replot the conductance at 150 K (black line). The behaviour remains highly incongruous with the single-particle model (nx = 0, blue dotted line), and is again more similar to the calculation for the case of maximal nx with slightly decreased mobilities (red dashed line). However, there is now a discrepancy in the form of a vertical shift, equivalent to an extra gate-independent contribution to the conductance. The vertical shift cannot be accounted for by varying nx; to illustrate this we also plot (green dash-dotted line) the result of assuming that nx varies with gate voltage according to a chemical equilibrium condition, nx = K(n − nx)(p − nx), where K is an equilibrium constant. It also cannot be reproduced by using a single-electron spectrum, for example with overlapping bands. A more sophisticated treatment of the correlated-electron system may therefore be needed to understand this aspect of the behaviour.

Excitons that persist in equilibrium at 100 K and at doping levels above 1 × 1012 cm−2 must have binding energy much larger than the thermal energy of ~10 meV and small size to survive screening by free charges. To see whether this is plausible, we solved the exciton (Bethe–Salpeter) equation of motion from first principles, building on the DFT band structure (Fig. 5a) and including spin–orbit effects in a non-perturbative way (Methods). The resulting excitation energy versus momentum q is shown in Fig. 5b. As the binding energy only weakly depends on Eg because the gap is indirect2,10, and the value of Eg is uncertain, we tuned the DFT hybrid functional to make Eg vanish. The dielectric function was evaluated in the random phase approximation, and the uncertainty induced by the numerical discretization of k space is shown by the error bars. The excitation energy is negative for all q, ranging from −100 meV for direct excitons at q = 0 to a minimum of −330 meV for indirect excitons made of a hole at Γ and an electron at Λ. In Fig. 5c,d we plot the spatial profile of an exciton in the centre-of-mass frame. The exciton radius is as small as 4 nm. This is comparable with the typical electron separation at the critical doping \({{{n_{\rm{ce}}}^{{ - 1/2}}} = {4.5\,{\rm{nm}}}}\), suggesting that excitons could play a role in the insulator–metal transition at nce.

Fig. 5: Calculated exciton properties.
figure 5

a, Band structure along Γ–Λ–X in Fig. 1b, obtained from first principles (DFT-PBE0 level, Methods). The lines are guides to the eye. b, Excitation energy of the lowest exciton versus centre-of-mass momentum q = (q, 0) along the same cut, calculated by solving the Bethe–Salpeter equation. The thick curve is a guide to the eye. The excitation energy is negative over the whole range of q, implying instability of the ground state due to the spontaneous generation of excitons. Error bars are estimated by extrapolation of Brillouin zone sampling (Methods and Supplementary Fig. 7). c,d, Wavefunction of the lowest exciton with q = 0 (c) and q = kΛ (d), respectively. The plot shows the conditional probability distribution of the electron for the hole located at the red dot. The plot in d is an average over three degenerate states. e, Simulated behaviour of μ versus ng (black solid line) for an excitonic insulator within the spinful two-band model (Methods). The density of free electrons (red dashed lines) and holes (blue dashed lines) is computed self-consistently as the population of the c and v bands, respectively, which are renormalized by the presence of the condensate.

Source data

In interpreting the situation below 100 K in terms of an excitonic insulator, formed by condensation of the excitons, we face two problems. First, at low temperatures, insulating behaviour sets in over a wide doping range, approaching ncp < ng < nce. Conventional neutral excitonic insulator theory provides no mechanism for localizing the unbalanced charge. Nevertheless, it seems possible that the Coulomb interaction, which has a long range for these doping values, could stabilize both the excitonic phase10 and the Wigner crystallization of unbound carriers, whose effective mass is enhanced by the opening of a many-body gap. Second, the exciton condensate is naively expected to exhibit a charge density wave (CDW) with wavevector kΛ, but CDWs are not seen in tunnelling microscopy20,24,33 and our detailed temperature-dependent Raman spectra show no evidence of any CDW transition (Supplementary Section 3). One possible explanation is that the condensate is made of direct excitons having q = 0, as predicted25 for the T′ phase of monolayer MoS2. We have checked that this leads to no substantial symmetry breaking, due to the anisotropic character of the WTe2 band structure (Methods and Supplementary Fig. 8). However, this state would have higher energy than a condensate made of indirect excitons with finite q. A more likely possibility is that the peculiar symmetries of excitons with q = ±kΛ prevent the condensate from exhibiting charge order.

A mean-field theory of this condensate (Methods), which builds on the knowledge of indirect excitons obtained from first principles and takes into account both spin and valley degrees of freedom, shows that there is no CDW with momentum q = kΛ. We find that the charge order is hidden by the combined effect of time-reversal symmetry and spin–orbit interaction, which is ultimately related to the topological properties of WTe2 (ref. 15). The CDW may be unveiled by breaking time-reversal symmetry; practically, one may split the hole states that sustain the density wave through Zeeman coupling with the magnetic field, as the hole spin is polarized along x close to Γ. Finally, the ground state exhibits a spin density wave of momentum q = kΛ, plus a weak charge modulation of period q = 2kΛ. These features were also found in ref. 28 by adding an intervalley scattering term to a two-band model, in the absence of which different ordered states would have been degenerate. In contrast, the hidden order we predict here is inherent to all possible spin symmetries of the condensate, and hence robust.

We also consider the behaviour of μ versus ng. We start by putting forward the following heuristic argument. At low temperatures (T 50 K), where every minority carrier is paired, for ng < 0 when an electron is added to the system it pairs with a hole, reducing the addition energy by the exciton binding energy and leading to diverging compressibility and self-consistent pinning of μ (in the charge-ordered state). For ng > 0, the added electron does not pair and so μ is not affected by the binding energy; the result is a step in μ at ng = 0 whose height is related to the binding energy. The decrease in the height of the step for T 100 K could be because the exciton binding weakens due to screening by the free carriers.

This scenario is supported by simulating the behaviour of an excitonic insulator in the presence of free charge carriers (Fig. 5e), with both μ and the many-body gap being computed self-consistently (Methods). Here the gap has a purely excitonic origin, as the starting non-interacting phase is a semimetal. The step in μ remains clearly visible up to 100 K, in contrast with the smeared profile of the independent-electron model (Fig. 3d). This is a peculiar consequence of exciton condensation, as electrons injected for ng > 0 fill in the lowest c states blockading the formation of e–h pairs—an effect due to the Pauli exclusion principle and suppressed with temperature34. At even higher T the step is smeared anyway, because the condensate is depopulated by thermal excitations and the excitonic gap melts. The simulation of Fig. 5e was performed for direct excitons for the sake of illustration25, but indirect excitons will exhibit the same qualitative features.

In conclusion, examination of the conductivity and thermodynamics of monolayer WTe2 provides evidence that neutral excitons are present in thermal equilibrium, not only at low temperatures where they probably form a collective excitonic insulator state, but also even at temperatures above 100 K when they coexist with free carriers. The size of the step in chemical potential and theoretical calculations suggest that the exciton binding energy is at least 40 meV.


Ground-state calculations from first principles

We obtained the ground-state electronic structure within DFT with a plane wave basis set, as implemented in the Quantum ESPRESSO package35. We fixed a kinetic energy cutoff of 80 Ry for the wavefunctions and used fully relativistic norm-conserving pseudopotentials36 to include the spin–orbit interaction. We optimized the lattice parameters and atomic positions using the PBE exchange-correlation functional, the final cell parameters being a = 3.52 Å and b = 6.29 Å. We set the cell side along z to 15 Å. We obtained the band structure using a PBE0 pseudopotential, for which we considered a small fraction of exact exchange, 2%.

Excitation energies and exciton wavefunctions from first principles

We calculated the excitation energies of excitons as well as the dispersion of the lowest-energy exciton within the framework of many-body perturbation theory37,38,39, by solving the Bethe–Salpeter equation through the YAMBO code40,41 and including spin–orbit interaction in a non-perturbative way42. We considered the PBE0 electronic structure as a starting point and calculated the static screening in the direct term within the random phase approximation, with inclusion of local field effects, and we employed the Tamm–Dancoff approximation for the Bethe–Salpeter Hamiltonian. To avoid spurious interactions among layers, we employed a truncated Coulomb cutoff technique43. We obtained converged excitation energies considering, respectively, two valence and two conduction bands in the Bethe–Salpeter matrix, the irreducible Brillouin zone being sampled up to a 48 × 24 × 1 k-point grid. We extrapolated the excitation energy of the lowest exciton with momentum q = 0 to the limit of a dense k-point grid, as shown in Supplementary Fig. 7.

Excitonic insulator ground state from indirect excitons and synopsis of theory

From the first-principles solution of Bethe–Salpeter equation we find that the lowest-energy exciton with q = kΛ (or −kΛ) is three-fold degenerate and separated by 20 meV from a non-degenerate first excited state. This energy splitting is due to the residual exchange interaction, present despite the strong spin–orbit coupling. These excitons are made of electrons and holes that populate, respectively, the lowest c and highest v band, with momentum k lying near the ΓΛ line. Along this high-symmetry line, both c and v Bloch spinors, each doubly degenerate, may be chosen as the eigenstates of the two-fold screw-axis rotation, which has complex conjugated values as irreducible representations. We label these Bloch spinors as ζ= ±i and use them to write explicitly the exciton wavefunctions within a minimal two-band model, which includes both spin and orbital degrees of freedom (Spinful two-band model). In addition to the rotational symmetry, we may classify the excitons according to their triplet- or singlet-like character, that is, whether they respectively maximize or minimize the electron–hole spatial overlap. In fact, Bloch spinors labelled by ζ transform like spins polarized along the x axis under the two-fold rotation, even if their actual spin polarization away from Γ is zero. Specifically, the exciton wavefunctions that are even with respect to the screw-axis rotation are

$${{\left| \pm \right\rangle} = {\frac{1}{{\sqrt 2 }}{\mathop {\sum }\limits_{k_x,k_y}}{ \psi ^ \pm \left( {{{\bf{k}}}} \right)}{\left[ {{c_i^{\dagger}} \left( {k_x + k_{{{\mathrm{{\Lambda}}}}},\,k_y} \right){v_i}\left( {{{\bf{k}}}} \right) \pm {c_{ - i}^{\dagger}} \left( {k_x + k_{{{\mathrm{{\Lambda}}}}},\,k_y} \right){v_{ - i}}({{{\bf{k}}}})} \right]}{\left| 0 \right\rangle},}}$$

where +/− stands for singlet/triplet-like symmetry, \({{c_\zeta ^{\dagger}} \left( {{{\bf{k}}}} \right)}\) creates and \({v_\zeta \left( {{{\bf{k}}}} \right)}\) destroys an electron of momentum k and screw-axis symmetry ζ in the c and v band, respectively, |0〉 is the non-interacting ground state with all v states filled and c states empty, and ψ is the exciton wavefunction in reciprocal space in the centre-of-mass frame, which is nodeless and even in ky (and, approximately, in kx): ψ(kx, ky) =ψ(kx, −ky).

In the excitonic insulator phase, the condensate is macroscopically occupied by excitons, hence the expectation value of the operator that creates an electron–hole pair, \({\left\langle {{c_\zeta ^ +} {v_{\zeta \prime }}}\right\rangle}\), has finite magnitude ϕζζ and arbitrary phase θζζ; that is, the complex wavefunction of the condensate is ϕζζexp(ζζ) (here \(\left\langle \ldots \right\rangle\) is the average over the many-body ground state). As excitons with q = (kΛ, 0) and (−kΛ, 0) are degenerate, they may separately condense, so the wavefunction of the condensate has two valley components of equal magnitude, whose respective phases are related by time-reversal symmetry10. Furthermore, both condensate and exciton wavefunctions share the same screw-axis and spin-like symmetries3, as well as the same parity in k space. These fundamental constraints relate all components of the condensate to a unique wavefunction magnitude, ϕ(k), and phase, θ. If the excitonic ground state has triplet-like symmetry, one has

$${{\left\langle {c_\zeta ^{\dagger} \left( {k_x \pm k_{{{\mathrm{{\Lambda}}}}},\,k_y} \right)v_{\zeta \prime }\left( {{{\bf{k}}}} \right)} \right\rangle} = {\pm \left[ {\sigma _z} \right]_{\zeta \zeta \prime }\exp \left( { \pm i\theta } \right)\phi ^ - ({{{\bf{k}}}}).}}$$

For the singlet-like ground state, one replaces the 2 × 2 Pauli matrix σz with the identity matrix 1, the right-hand side of the equation then reading 1exp(±)ϕ+(k). The charge/spin density wave of the excitonic insulator is dictated by the interband contribution to the expectation value of the corresponding density operator, which is proportional to the left-hand side of the above equation. Therefore, we may assess the occurrence of charge order in the excitonic phase without actually computing ϕ±(k), which is given by a gap equation10 that depends on Eg. Indeed, ϕ only provides the intensity of the density modulation, whereas θ rigidly shifts the density wave with respect to the frame origin.

Spinful two-band model

The spinful two-band model provides the non-interacting, doubly degenerate c and v Bloch states of crystal momentum k that comply with the symmetry group of monolayer WTe2 (the T′ structure is centrosymmetric and non-symmorphic). Within the 4D spin/orbital space, the Hamiltonian is a 4 × 4 Hamiltonian matrix, HQSH(k), whose off-diagonal elements are the spin–orbit interaction terms. Here the orbital degree of freedom identifies the c and v Bloch states at Γ, whose energies are ‘inverted’ with respect to the usual order of bulk semiconductors. These two states have been variously identified in the literature as a pair of orbitals having either opposite15,25 or like29,44 parities under spatial inversion, leading to two different forms of the spin–orbit interaction: we label the corresponding model Hamiltonians, respectively, as HQSH,1(k) and HQSH,2(k). We find that the charge order is hidden in the condensate of indirect excitons, regardless of the model. The reason is that the screw-axis rotation (around the W atom chain direction) maintains the same form within the spin/orbital space, entangling the degrees of freedom of its eigenstates, labelled by ζ= ±i. This is pivotal to the discussion of the main text.

In the following we detail the Hamiltonians HQSH,1(k) and HQSH,2(k). The first model, HQSH,1(k), is taken from refs. 15,25 with minor adjustments. The Hamiltonian, written in the present reference frame with the W atom chain parallel to the x axis, reads

$$\begin{array}{rcl}{H_{{\rm{QSH}},{1}}}\left( {{{\bf{k}}}} \right) &=& {{\frac{1}{2}}\left[ {{{\it{\epsilon }}_u\left( {{{\bf{k}}}} \right)} + {{\it{\epsilon }}_g\left( {{{\bf{k}}}} \right)}} \right]{\bf{1}}_\tau} \otimes {{\bf{1}}_\sigma} + {\frac{1}{2}\left[ {{{\it{\epsilon }}_u\left( {{{\bf{k}}}} \right)} - {{\it{\epsilon }}_g\left( {{{\bf{k}}}} \right)}} \right]\tau _z} \otimes {{\bf{1}}_\sigma}\\ &&+ {\hbar v_2k_x\tau _x} \otimes {\sigma _y},\end{array}$$

where v2 is the spin–orbit coupling parameter, τx, τy, τz and σx, σy, σz are 2 × 2 Pauli matrices in orbital and spin space, respectively, and the 2 × 2 unit matrices are 1τ and 1σ. The diagonal matrix elements, εl(k) with l = u, g, are the band energies in the absence of spin–orbit interaction, which are inverted at Γ and cross at the points ±kΛ of the ΓX line in the Brillouin zone. The energies εl(kx, ky) are even with respect to both kx and ky axes. The corresponding Bloch states transform like px and dxz orbitals at Γ. The functional dependence of εl on k differs from the effective-mass expression given in ref. 25, but its precise form is irrelevant to the discussion of the main text, solely based on symmetry arguments. We discard the off-diagonal spin–orbit term linear in σx considered in ref. 25, the Hamiltonian being now even in ky, HQSH,1(kx, ky) = HQSH,1(kx, −ky). This choice agrees with the evidence that the spin–orbit field lies in the xz plane29, as also predicted by ref. 45. Furthermore, we find that the spin–orbit term proportional to σy provides a good matching with the strongly anisotropic DFT bands, as we checked through comparison with our own first-principles calculations. In the spin/orbital space, the inversion operator reads \({{{\bf{I}}}} = {{- \tau _z} \otimes {1_\sigma}}\) and the screw-axis rotation around the x axis is \({{{{\bf{C}}}}_{2x}} = {{i\tau _z} \otimes {\sigma _x\exp (ik_xa/2)}}\), with a being the lattice constant in the direction parallel to the W chains.

The second model, HQSH,2(k), is taken from ref. 29 and builds on the four-band tight-binding Hamiltonian proposed in ref. 44 to improve the matching between model and DFT bands. The c and v orbital states are Wannier functions, respectively an antibonding combination of \({d_{{x^2} - {y^2}}}\)-type orbitals localized on W atoms (energy εW(k)) and a bonding superposition of px-type orbitals localized on Te atoms (εTe(k)). These Wannier functions have the same parities under inversion but opposite parities under the screw-axis two-fold rotation, like the Bloch states at Γ of our own DFT calculations. The Hamiltonian is

$$\begin{array}{rcl}H_{{\rm{QSH}},2}\left( {{{\bf{k}}}} \right) &=& {\frac{1}{2}\left[ {{{\it{\epsilon }}_{\rm{W}}\left( {{{\bf{k}}}} \right)} + {{\it{\epsilon }}_{\rm{Te}}\left( {{{\bf{k}}}} \right)}} \right]1_\tau} \otimes {1_\sigma} + {\frac{1}{2}\left[ {{{\it{\epsilon }}_{\rm{W}}\left( {{{\bf{k}}}} \right)} - {{\it{\epsilon }}_{\rm{Te}}\left( {{{\bf{k}}}} \right)}} \right]{\tau _z}} \otimes {1_\sigma}\\ &&+ {\lambda _{\rm{SO}}\tau _y} \otimes {\sigma _y},\end{array}$$

where λSO > 0 is the spin–orbit coupling parameter, which is independent of k. For the sake of simplicity, here we have neglected the spin–orbit term proportional to σz proposed by ref. 29. Within the envelope function approximation, the band energies εW(k) and εTe(k) are provided by the tight-binding calculation of ref. 44 and are even with respect to both kx and ky axes. The inversion operator now reads \({{{\bf{I}}}} = {{{\bf{1}}_\tau} \otimes {{\bf{1}}_\sigma}}\) whereas the screw-axis rotation around the x axis is again \({{{{\bf{C}}}}_{2x}} = {i\tau _z} \otimes {\sigma _x\exp (ik_xa/2)}\).

Eigenstates of the screw-axis rotation

The eigenvectors of HQSH,1(k) that were explicitly given in ref. 25 are spin-polarized along the direction perpendicular to the W atom chains. Here we use the notation \({{\left| {{{{\bf{k}}}},\lambda} \right\rangle} }\) to identify these eigenvectors, which belong to either the c or v band and have spin polarization λ = ↑,↓. In the main text we introduced an alternative, equally legitimate set of eigenvectors, which are simultaneous eigenstates of HQSH,1(k) and C2x (the latter with eigenvalues ζ= ±i), by applying a unitary rotation of the basis for any wavevector k. Explicitly,

$${{\left| {{{{\bf{k}}}},\zeta = - i} \right\rangle} = {\left( {i\left| {{{{\bf{k}}}}, \uparrow } \right\rangle + \left| {{{{\bf{k}}}}, \downarrow } \right\rangle } \right)}/{\surd 2}},$$
$${{\left| {{{{\bf{k}}}},\zeta = + i} \right\rangle} = {\left( { - i\left| {{{{\bf{k}}}}, \uparrow } \right\rangle + \left| {{{{\bf{k}}}}, \downarrow } \right\rangle } \right)/\surd 2}},$$


$${{{{\bf{C}}}}_{2x}\left| {{{{\bf{k}}}},\zeta = - i} \right\rangle} = {- i\exp (ik_xa/2)\left| {{{{\bf{k}}}},\zeta = - i} \right\rangle}$$


$${{{{{\bf{C}}}}_{2x}}{\left| {{{{\bf{k}}}},\zeta = + i} \right\rangle} = {i\exp (ik_xa/2)\left| {{{{\bf{k}}}},\zeta = {+ i}} \right\rangle}},$$

as may be checked by direct substitution (followed by rotation of the original frame axis of ref. 25). Importantly, the states ζ= ±i are not spin-polarized (except at k = 0), because the spin and orbital degrees of freedom are now entangled. Note that \({\left| {{{{\bf{k}}}},\zeta = - i} \right\rangle}\) and \({\left| { - {{{\bf{k}}}},\zeta = + i} \right\rangle}\) are time-reversal mates, with \({{\varTheta}\left| {{{{\bf{k}}}},\zeta = - i} \right\rangle = i\left| { - {{{\bf{k}}}},\zeta = + i} \right\rangle}\), the time-reversal operator being \({{\varTheta} = {i{\bf{1}}_\tau \otimes \sigma _yK}}\) (K is the complex conjugation operator).

The simultaneous eigenvectors of HQSH,2(k) and C2x that we use to build the excitonic insulator ground state (within the envelope function approximation) are as follows. The c band state with k = (0, kΛ) and ζ= +i is \({( {- 1} + {i})/{(2\sqrt 2 )}{[1,1,1,- 1]}}\), and the v band state with k = 0 and ζ= +i is (approximately) \({{\left( {1/\sqrt 2 } \right)}{[0, - 1,0,1]}}\). The states with ζ= −i as well as those with k = (0, −kΛ) are obtained through time-reversal and inversion transformations. Here the first and third (second and fourth) components of the 4D vector, \({[{u_{{\rm{W}}, \uparrow }},{v_{{\rm{Te}}, \uparrow }},{u_{{\rm{W}}, \downarrow }},{v_{{\rm{Te}}, \downarrow }}]}\), correspond to a spinor whose orbital part is a Wannier function localized on W (Te) atoms in the crystal unit cell.

Condensate of indirect excitons and charge/spin order

We assess the charge (spin) order of a permanent condensate of indirect excitons of momentum q = (±kΛ, 0) within the multivalley framework developed in ref. 10. This approach, in turn, relies on the scheme to decouple the equations of motion for Green functions of ref. 46. The theory, which deals with spinless electrons, may be straightforwardly generalized to spinors labelled by the screw-axis symmetry =±i. Indeed, for any given electron and hole species ζ and ζ′, the structure of the equations of motion for Green functions that govern the condensate component \({\left\langle {{c_\zeta ^ +} {v_{\zeta \prime }}} \right\rangle}\) remains the same, because ζ electrons pair with ζ′ holes only. As far as pairing is concerned, ζ spinors behave as if they were spinless fermions. Therefore, the k-dependent spinless Bogoliubov–Valatin-like creation operator that defines the excitonic insulator ground state in equations (20) and (21) of ref. 10 is simply replicated for ζ = ±i, provided one specializes equation (21) to the present case of two valley components and chooses the condensate phases as shown above (Excitonic insulator ground state from indirect excitons and synopsis of theory). Finally, we compute the expectation value of the charge (spin) density operator over the excitonic ground state, after making the spin/orbital structure of ζ Bloch states given in the previous section explicit, the derivation being lengthy but straightforward. As discussed above, to simply assess the occurrence of charge (spin) order without computing the density wave modulation intensity, it is not necessary to evaluate explicitly the coherence coefficients \({u_k^0}\) and \({v_k^0}\) that occur in equation (21) of ref. 10 (these are provided by the self-consistent gap in equation (3)).

Condensate of direct excitons and simulation of Fig. 5e

We obtain the results shown in Fig. 5e within the spinful two-band model HQSH,1(k) described above, the band energies εl(k) being parametrized through comparison with our own DFT calculations (l = u, g). The non-interacting ground state, in the absence of spin–orbit interaction, is taken to be a semimetal with a band overlap of 38 meV. The energy parametrization relies partly on the tight-binding model of ref. 44 (table III therein) and partly on ad hoc parameters. In detail, we correct the tight-binding energies by adding the terms \({{2\left( {{t_l^\prime} - {t_l^{\prime\prime}} } \right){\cos 2k_x}} + {{2t_l^{\prime\prime}} \cos 2\left| {{{\bf{k}}}} \right|}}\), with \({t_g^\prime} = {0.149}\,{\rm{eV}}\), \({t_u^\prime} = {0.075}\,{\rm{eV}}\) (ref. 23), \({t_g^{\prime\prime}} = {0.049}\,{\rm{eV}}\) and \({t_u^{\prime\prime}} = {0.055}\). The spin–orbit parameter is v2 = 1 × 1014 Å s−1, and the strength of the Coulomb interaction is fixed by the 2D polarizability, α2D = 5.5 (the notation of ref. 25). To compute the chemical potential μ versus charge density ng in the many-body excitonic phase, we adapt the theory of ref. 25, which deals with a condensate of direct excitons in an intrinsic semiconductor, to the case of a doped system, by means of a fully self-consistent calculation of both the excitonic order parameter, ΔX(k), and μ. Furthermore, we assume that the Coulomb interaction remains long-ranged for any doping value, which is supported by the evidence that charge carriers localize in a wide doping interval at low T. The free carriers populating the renormalized bands of the excitonic insulator are conventionally taken as non-interacting, which is the origin of the unphysical behaviour of μ for small, positive values of ng at T = 25 K (dμ/dng is negative close to the axis origin in Fig. 5e).

We have checked that the observable effects related to the breaking of inversion symmetry, due to the condensation of excitons with q = 0, are negligible for WTe2, in contrast to the case of T′-MoS2. This is related to the limited extent of the excitonic order parameter, ΔX(k), along the Brillouin zone direction perpendicular to the ΓΛ line, which is in turn caused by the strong anisotropy of the non-interacting c and v bands close to Λ (compare Fig. 5a with Supplementary Fig. 8). In particular, the real part of ΔX(k) is negligible, so the c and v bands, renormalized by electron–hole pairing, each remain doubly degenerate.

Data availability

The datasets generated during and/or analysed during the current study are available from the corresponding authors on reasonable request. Source data are provided with this paper.

Code availability

Many-body perturbation theory calculations were performed using the codes YAMBO ( and Quantum ESPRESSO (, which are both open-source software. Results for the two-band model were obtained through custom Fortran codes that are available from M.R. upon reasonable request.


  1. Keldysh, L. & Kopaev, Y. Possible instability of semimetallic state toward Coulomb interaction. Sov. Phys. Sol. State 6, 2219–2224 (1964).

    Google Scholar 

  2. Cloizeaux, J. D. Exciton instability and crystallographic anomalies in semiconductors. J. Phys. Chem. Solids 26, 259–266 (1965).

    ADS  Google Scholar 

  3. Jerome, D., Rice, T. M. & Kohn, W. Excitonic insulator. Phys. Rev. 158, 462–475 (1967).

    ADS  Google Scholar 

  4. Halperin, B. I. & Rice, T. M. Possible anomalies at a semimetal–semiconductor transition. Rev. Mod. Phys. 40, 755–766 (1968).

    ADS  Google Scholar 

  5. Rohwer, T. et al. Collapse of long-range charge order tracked by time-resolved photoemission at high momenta. Nature 471, 490–493 (2011).

    ADS  Google Scholar 

  6. Kogar, A. et al. Signatures of exciton condensation in a transition metal dichalcogenide. Science 358, 1314–1317 (2017).

    ADS  Google Scholar 

  7. Lu, Y. F. et al. Zero-gap semiconductor to excitonic insulator transition in Ta2NiSe5. Nat. Commun. 8, 14408 (2017).

    ADS  Google Scholar 

  8. Werdehausen, D. et al. Coherent order parameter oscillations in the ground state of the excitonic insulator Ta2NiSe5. Sci. Adv. 4, eaap8652 (2018).

    ADS  Google Scholar 

  9. Varsano, D. et al. Carbon nanotubes as excitonic insulators. Nat. Commun. 8, 1461 (2017).

    ADS  Google Scholar 

  10. Ataei, S. S., Varsano, D., Molinari, E. & Rontani, M. Evidence of ideal excitonic insulator in bulk MoS2 under pressure. Proc. Natl Acad. Sci. 118, e2010110118 (2021).

    Google Scholar 

  11. Eisenstein, J. P. & MacDonald, A. H. Bose–Einstein condensation of excitons in bilayer electron systems. Nature 432, 691–694 (2004).

    ADS  Google Scholar 

  12. Nandi, D., Finck, A. D. K., Eisenstein, J. P., Pfeiffer, L. N. & West, K. W. Exciton condensation and perfect Coulomb drag. Nature 488, 481–484 (2012).

    ADS  Google Scholar 

  13. Li, J. I. A., Taniguchi, T., Watanabe, K., Hone, J. & Dean, C. R. Excitonic superfluid phase in double bilayer graphene. Nat. Phys. 13, 751–755 (2017).

    Google Scholar 

  14. Liu, X., Watanabe, K., Taniguchi, T., Halperin, B. I. & Kim, P. Quantum Hall drag of exciton condensate in graphene. Nat. Phys. 13, 746–750 (2017).

    Google Scholar 

  15. Qian, X., Liu, J., Fu, L. & Li, J. Quantum spin Hall effect in two-dimensional transition metal dichalcogenides. Science 346, 1344–1347 (2014).

    ADS  Google Scholar 

  16. Fei, Z. et al. Edge conduction in monolayer WTe2. Nat. Phys. 13, 677–682 (2017).

    Google Scholar 

  17. Wu, S. et al. Observation of the quantum spin Hall effect up to 100 kelvin in a monolayer crystal. Science 359, 76–79 (2018).

    ADS  MathSciNet  MATH  Google Scholar 

  18. Sajadi, E. et al. Gate-induced superconductivity in a monolayer topological insulator. Science 362, 922–925 (2018).

    ADS  Google Scholar 

  19. Fatemi, V. et al. Electrically tunable low-density superconductivity in a monolayer topological insulator. Science 362, 926–929 (2018).

    ADS  MathSciNet  Google Scholar 

  20. Tang, S. et al. Quantum spin Hall state in monolayer 1T′-WTe2. Nat. Phys. 13, 683–687 (2017).

    Google Scholar 

  21. Cucchi, I. et al. Microfocus laser-angle-resolved photoemission on encapsulated mono-, bi- and few-layer 1T′-WTe2. Nano Lett. 19, 554–560 (2019).

    ADS  Google Scholar 

  22. Zheng, F. et al. On the quantum spin Hall gap of monolayer 1T′-WTe2. Adv. Mater. 28, 4845–4851 (2016).

    Google Scholar 

  23. Ok, S. et al. Custodial glide symmetry of quantum spin Hall edge modes in monolayer WTe2. Phys. Rev. B 99, 121105 (2019).

    ADS  Google Scholar 

  24. Song, Y.-H. et al. Observation of Coulomb gap in the quantum spin Hall candidate single-layer 1T′-WTe2. Nat. Commun. 9, 4071 (2018).

    ADS  Google Scholar 

  25. Varsano, D., Palummo, M., Molinari, E. & Rontani, M. A monolayer transition-metal dichalcogenide as a topological excitonic insulator. Nat. Nanotechnol. 15, 367–372 (2020).

    ADS  Google Scholar 

  26. Jia, Y. et al. Evidence for a monolayer excitonic insulator. Nat. Phys. (2021).

  27. Lee, P. A. Quantum oscillations in the activated conductivity in excitonic insulators: possible application to monolayer WTe2. Phys. Rev. B 103, L041101 (2021).

    ADS  Google Scholar 

  28. Kwan, Y. H., Devakul, T., Sondhi, S. L. & Parameswaran, S. A. Theory of competing excitonic orders in insulating WTe2 monolayers. Phys. Rev. B 104, 125133 (2021).

    ADS  Google Scholar 

  29. Zhao, W. et al. Determination of the spin axis in quantum spin Hall insulator monolayer WTe2. Preprint at (2021).

  30. Wang, P. et al. Landau quantization and highly mobile fermions in an insulator. Nature 589, 225–229 (2021).

    ADS  Google Scholar 

  31. Fei, Z. et al. Ferroelectric switching of a two-dimensional metal. Nature 560, 336–339 (2018).

    ADS  Google Scholar 

  32. Shi, Y. et al. Imaging quantum spin Hall edges in monolayer WTe2. Sci. Adv. 5, eaat8799 (2019).

    ADS  Google Scholar 

  33. Lupke, F. et al. Proximity-induced superconducting gap in the quantum spin Hall edge state of monolayer WTe2. Nat. Phys. 16, 526–530 (2020).

    Google Scholar 

  34. Zittartz, J. Anisotropy effects in the excitonic insulator. Phys. Rev. 162, 752–758 (1967).

    ADS  Google Scholar 

  35. Giannozzi, P. et al. Advanced capabilities for materials modelling with Quantum ESPRESSO. J. Phys. Condens. Matter 29, 465901 (2017).

    Google Scholar 

  36. Hamann, D. R. Optimized norm-conserving Vanderbilt pseudopotentials. Phys. Rev. B 88, 085117 (2013).

    ADS  Google Scholar 

  37. Onida, G., Reining, L. & Rubio, A. Electronic excitations: density-functional versus many-body Green’s-function approaches. Rev. Mod. Phys. 74, 601–659 (2002).

    ADS  Google Scholar 

  38. Hybertsen, M. S. & Louie, S. G. Electron correlation in semiconductors and insulators: band gaps and quasiparticle energies. Phys. Rev. B 34, 5390–5413 (1986).

    ADS  Google Scholar 

  39. Strinati, G. Application of the Green’s functions method to the study of the optical properties of semiconductors. Riv. Nuovo Cim. 11, 1–86 (1988).

    ADS  Google Scholar 

  40. Marini, A., Hogan, C., Grüning, M. & Varsano, D. Yambo: an ab initio tool for excited state calculations. Comput. Phys. Commun. 180, 1392–1403 (2009).

    ADS  Google Scholar 

  41. Sangalli, D. et al. Many-body perturbation theory calculations using the yambo code. J. Phys. Condens. Matter 31, 325902 (2019).

    Google Scholar 

  42. Elliott, J. D. et al. Surface susceptibility and conductivity of MoS2 and WSe2 monolayers: a first-principles and ellipsometry characterization. Phys. Rev. B 101, 045414 (2020).

    ADS  Google Scholar 

  43. Rozzi, C. A., Varsano, D., Marini, A., Gross, E. K. U. & Rubio, A. Exact Coulomb cutoff technique for supercell calculations. Phys. Rev. B 73, 205119 (2006).

    ADS  Google Scholar 

  44. Muechler, L., Alexandradinata, A., Neupert, T. & Car, R. Topological nonsymmorphic metals from band inversion. Phys. Rev. X 6, 041069 (2016).

    Google Scholar 

  45. Garcia, J. H. et al. Canted persistent spin texture and quantum spin Hall effect in WTe2. Phys. Rev. Lett. 125, 256603 (2020).

    ADS  Google Scholar 

  46. Monney, C. et al. Spontaneous exciton condensation in 1T-TiSe2: BCS-like approach. Phys. Rev. B 79, 045116 (2009).

    ADS  Google Scholar 

Download references


The experiments and analysis were supported primarily by NSF DMR awards nos. MRSEC 1719797 and EAGER 1936697. Materials synthesis at the University of Washington was partially supported by the Gordon and Betty Moore Foundation’s EPiQS Initiative, grant no. GBMF6759 to J.-H.C. X.H. and Y.-T.C. acknowledge support from the NSF under award no. DMR-2004701 and a Hellman Fellowship award. The development of the devices was supported as part of Programmable Quantum Materials, an Energy Frontier Research Center funded by the US Department of Energy (DOE), Office of Science, Basic Energy Sciences (BES), under award no. DE-SC0019443. On the theory side, S.S.A., D.V., E.M. and M.R. were supported in part by the MaX European Centre of Excellence (‘MAterials design at the eXascale’, funded by the European Union H2020-INFRAEDI-2018-1 programme, grant no. 824143. D.V., E.M. and M.R. were also supported by the Italian national programme PRIN2017 no. 2017BZPKSZ ‘Excitonic insulator in two-dimensional long-range interacting systems (EXC-INS)’. M.P. was supported by the INFN Time2Quest project and by the European Union H2020-MSCA-RISE project DiSeTCom, grant no. 823728. S.S.A., M.P., D.V., E.M. and M.R. acknowledge access to the Marconi supercomputing system based at CINECA, Italy, through the PRACE and ISCRA programmes.

Author information




B.S., W.Z., T.P., Z.F. and E.R. made the devices and carried out the transport measurements. P.M. and J.-H.C. provided the crystals. J.C. and X.X. performed the Raman measurements. X.H. and Y.-T.C. carried out the MIM measurements. S.S.A., D.V., M.P., E.M. and M.R. developed the theory of excitonic effects and wrote the theoretical sections. D.H.C. led the experiments and their analysis. All authors contributed to all processes and to the manuscript.

Corresponding authors

Correspondence to Massimo Rontani or David H. Cobden.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature Physics thanks Vitor Pereira and the other, anonymous, reviewer(s) for their contribution to the peer review of this work.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Supplementary Information

Supplementary Figs. 1–8 and Discussion.

Source data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Sun, B., Zhao, W., Palomaki, T. et al. Evidence for equilibrium exciton condensation in monolayer WTe2. Nat. Phys. 18, 94–99 (2022).

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI:

Further reading


Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing