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.

Multiband charge density wave exposed in a transition metal dichalcogenide

This article has been updated


In the presence of multiple bands, well-known electronic instabilities may acquire new complexity. While multiband superconductivity is the subject of extensive studies, the possibility of multiband charge density waves (CDWs) has been largely ignored so far. Here, combining energy dependent scanning tunnelling microscopy (STM) topography with a simple model of the charge modulations and a self-consistent calculation of the CDW gap, we find evidence for a multiband CDW in 2H-NbSe2. This CDW not only involves the opening of a gap on the inner band around the K-point, but also on the outer band. This leads to spatially out-of-phase charge modulations from electrons on these two bands, which we detect through a characteristic energy dependence of the CDW contrast in STM images.


Imposing a new periodicity on a crystal leads to the reorganization of the electronic bands of the parent compound through their back-folding on the new Brillouin zone. New periodicities may be engineered in designer materials, for instance in artificial heterostructures with Moiré minigaps, or emerge due to a structural or electronic phase transition. The charge density wave state is an electronic order where the charge density develops a spatial modulation concomitantly to a periodic distortion of the atomic lattice and the opening of a gap in the quasi-particle spectrum. By reducing the electronic band energy, this gap compensates for the elastic and Coulomb energy costs associated with the formation of the CDW. It also lowers the degeneracy of the electronic states at the crossings of the folded bands. These are the points in the band structure of the parent compound that are connected by the wavevector of the new periodicity. Although a gap should open at all the crossings of the folded bands, previous studies only focused on the primary CDW gap around the Fermi-level, which leads to the highest energy gain of the reconstructed system. The existence of secondary gaps and associated charge modulations (CMs) remains largely unexplored.

In many cases, only a tiny fraction of the electrons are involved in the CDW formation. Therefore, the CDW gap manifests only as a slight reduction of the density of states (DOS)—which can depend on momentum—rather than a full depletion of the DOS. This makes it challenging to measure the CDW gap using spectroscopic probes such as angle-resolved photoemission spectroscopy (ARPES) and scanning tunnelling spectroscopy. This is particularly true for 2H-NbSe21,2,3,4,5 (hereinafter simply NbSe2). However, the effect of the redistributed electrons can be readily detected in topographic STM images, even for minute changes brought upon by the opening of the CDW gap as demonstrated in the following.

NbSe2 is an iconic material of correlated electron physics. It hosts a nearly commensurate charge density wave below TCDW = 33.5 K and a superconducting (SC) order below TSC = 7.2 K6,7,8,9,10. NbSe2 is a layered material with a three-fold symmetric crystal structure around the direction perpendicular to the layers (Fig. 1a). Each unit cell is composed of two slabs of Se–Nb–Se trilayers, where the Se lattices are 60° rotated, while the Nb atoms are aligned on top of each other in a trigonal prismatic coordination with the Se atoms.

Fig. 1: Crystal structure and Fermi surface of NbSe2.

a The three-fold symmetric crystal structure from top and side views33. b The Fermi surface has been calculated using a two-band tight-binding fit to ARPES data3. Inner pockets (red) around Γ and K derive from one band, while the outer pockets (blue) derive from the second band; a small pancake-shaped pocket around Γ originating mainly from Se orbitals has been omitted.

The Fermi surface (FS) of NbSe2 is mainly determined by the bonding and antibonding combinations of the Nb-4d orbitals11,12 leading to double-walled barrel-shaped pockets around the K and Γ points of the hexagonal Brillouin-zone2,3,13,14 (Fig. 1b). The charge-ordered state consists of three CDWs, which form along the three equivalent ΓM directions with wavevectors \((1-\delta )\frac{2}{3}| {{\Gamma }}M|\), where δ ≈ 0.02 and depends on temperature10. In real space, this yields a locally commensurate 3a0 × 3a0 superstructure delimited by discommensurations15,16, where a0 is the atomic periodicity.

The 3a0 × 3a0 reconstruction is readily accessible to topographic STM imaging. Its bias-dependent contrast has been the focus of previous studies, with particular emphasis on the contrast inversion expected in a classic Peierls scenario between images acquired above and below the CDW gap17,18, and on the role of defects in stabilizing the CDW5. Sacks et al.18 calculate the bias dependence of the CDW phase in a perturbative approach, considering a single band normal state description of NbSe2. They find that the phase-shift of the CDW component of the local DOS can be very different from the 180° expected in a one-dimensional (1D) case (Supplementary Note VII) when changing the imaging bias across the Fermi level (EF). However, their model does not reproduce the full bias dependence of the CDW amplitude and phase that we find.

Results and discussion

Bias dependent STM topography, CDW phase and amplitude

In Fig. 2a–c, we present a selection from numerous topographic STM images of the same region on a cleaved NbSe2 surface at different sample biases (Vb) between −0.5 V and 0.5 V. They show a triangular atomic lattice with a superimposed 3a0 × 3a0 CDW modulation (see also Supplementary Fig. 1a), consistent with previous STM studies of unstrained bulk NbSe24,5,16,17,19,20,21,22. A defect-free region with a well resolved CDW outlined in red is magnified in Fig. 2d, f, h for each Vb. In order to identify the origin of the bias dependence of the topographic contrast in these images, we separate the atomic lattice and CDW contributions using Fourier filtering (Supplementary Note I). This analysis shows that the bias-dependent STM contrast is due to the changing CDW signal (Fig. 2e, g, i) since the corresponding atomic lattice contrast remains unchanged (see Supplementary Fig. 1).

Fig. 2: Bias-dependent STM imaging of NbSe2 at 1.2 K.

Constant current topography showing the atomic lattice and CDW at aVb = 50 mV, b Vb = − 50 mV and cVb = − 150 mV (Solely for visualization, the STM image in c has been corrected for diagonal running sharp lines (with wavelength much longer than those of the CDW) which arose due to a tiny horizontal tilt of the sample leading to a least significant bit issue in the digital–analogue conversion of the piezo voltage.) with It = 100 pA. d, f and h are magnified images of the areas marked by the red squares in a, b and c, respectively. The overall imaging contrast is very different in these cases, although the atomic lattice appears identical in all images (Supplementary Fig. 1). e, g and i shows the magnified image of the large Fourier-filtered image of the CDWs at the same location as shown in d, f and h, respectively. It demonstrates that the variation of contrast observed in d, f and h is stemming from a variation in the appearance of the CDW at different biases. These appearances can be quantified by a single parameter: the dephasing parameter which describes the relative position of the three CDWs. jl show the spatial variation of the dephasing parameter determined by fitting the CDW modulations of the STM images shown in a-c, respectively. The red squares correspond to the same area that is highlighted in ac and magnified in di. Scalebars: 10 nm in ac and jl; 1 nm in di.

The observed CDW pattern can be modelled as the sum of three plane waves as described in ref 16. While each plane wave has its own phase φi(r), which depends on a selected reference point, the dephasing parameter Θ(r) = φ1(r) + φ2(r) + φ3(r) mod 360 is uniquely defined for each particular CDW pattern, independent of any reference point. Θ(r) represents the internal CDW structure, quantifying the local relative position of the wavefronts of the three CDWs. In Fig. 2j–l, we show Θ(r) corresponding to the STM images in Fig. 2a–c, respectively. They were obtained by fitting the CDW contrast following the method described in ref. 16.

Each bias voltage is characterized by a dominant dephasing parameter (Fig. 2j–l), except in the vicinity of defects discussed later. This visual assessment is confirmed by the peaked histograms of Θ(r) (Supplementary Fig. 2). Fitting a Gaussian to these histograms allows to extract a well-defined dephasing parameter Θ0(Vb) for each imaging bias (Supplementary Fig. 3). For a quantitative analysis of the bias dependence of Θ0, we note that a given local CDW structure is represented by any arbitrary combination of φi(r) summing up to the same dephasing parameter, in particular the one where all three phases are equal. Moreover, the threefold symmetry of the system implies there is no privileged plane wave among the three used to describe the CDW. These observations allow us to map the problem to a one-dimensional (1D) description with an apparent CDW phase φ0(Vb) = Θ0(Vb)/3 (Supplementary Note III), and model φ0(Vb) to understand the bias dependent CDW pattern.

Plotting φ0(Vb) in Fig. 3a reveals a striking non-monotonic bias dependence, with an inflexion point around −0.15 V and a minimum slightly above the Fermi-level (EF = 0 V). This dependence is robust as long as φ0(Vb) is extracted from topographic STM images away from defects (Supplementary Fig. 4). Close to defects, the dephasing parameter Θ0(Vb) is different and tends to depend much less on imaging bias (Supplementary Fig. 5). This is consistent with earlier findings that defects (and impurities) can act as strong pinning centres23,24 locking the local phase of the CDW or driving the formation of CDW domains25,26.

Fig. 3: Bias dependence of the phase and amplitude in experiments and in modelling.

a and b show the bias dependence of the phase and amplitude of one of the unidirectional CDW extracted from the dephasing parameter and amplitude maps. c and d The phase and amplitude of the best matching simulations in the two-gap model. The data in b and d are normalized to their Vb = 0.5 V value.

The CDW amplitude can be extracted in a similar way to the phase, by fitting the histogram of the local amplitudes ai(r) of each plane wave measured over the entire field of view with a Gaussian, and extracting the peak value ai(Vb). The bias dependence and magnitude of ai(Vb) is nearly the same for all three CDWs (Supplementary Fig. 6b). For the analysis, we consider the average of these three amplitudes at each bias A0(Vb) = (a1(Vb) + a2(Vb) + a3(Vb))/3 plotted in Fig. 3b.

Modelling and calculations

To understand the bias dependence of the CDW amplitude and phase in Fig. 3, we simulate topographic STM traces using a 1D model system (Supplementary Note VIII). In the simplest configuration corresponding to the Peierls reconstruction, we consider the contribution to the tunneling current of a single CM and its associated gap centred on EF (Fig. 4a–c). In this case, traces at the same polarity are always in phase, whereas traces at opposite polarities always show contrast inversion (or a 180° phase shift in the present harmonic model). The latter, often considered as an identifying hallmark of the CDW state27, clearly does not reproduce the data in Fig. 3a.

Fig. 4: One-dimensional model description of the bias-dependent phase of the CDW.

ac There is a single CM and the corresponding gap is centred at the Fermi level. df There are two CMs with two gaps that are centred at different energies (one at the Fermi-level and one below). There is no real-space phase difference between the CMs. gi The same as df except that there is a 120° (2π/3) phase difference -one atom shift- between the two CMs. In all three cases, the first column shows the spatial and energy-dependent CDW local DOS maps. Second column: corresponding simulated topographic traces at selected biases (the actual bias value is shown in the legend). The red and blue topographic traces correspond to positive and negative sample biases, respectively. To clearly see the evolution of the phase and the amplitude the curves are offset vertically in b, e and h such that they all oscillate around the same value. Third column: bias dependence of the phase (left red axis) and DOS (spatially integrated local DOS from the first column, right blue axis).

A single CM can only produce two sets of STM traces differing by contrast inversion in the vicinity of the gap. To generate a more complex bias dependence of the phase, we consider the possibility of a second CM whose associated gap opens in another band and away from EF (Fig. 4d–f). If these two harmonic CMs are in phase (Fig. 4d), the resulting STM traces are either in-phase or 180° out of phase (Fig. 4e), unable to reproduce the data in Fig. 3a. To generate more structures in the bias dependence of the phase, we need to introduce a phase shift between the two CMs (Figs. 4g–i). This leads to a phase that is no longer bi-modal, limited to two values differing by 180° as in Fig. 4c, f. It takes many different values (Fig. 4i), where the precise bias dependence is defined by the magnitude of the two gaps, their position relative to EF and by the relative phase shift between the two CMs. The simulated STM topographic traces in Figs. 4b, e and h also reveal a pronounced bias-dependent imaging amplitude with distinct line-shapes in the three model cases discussed above (Supplementary Fig. 7).

The broad parameter space of our 1D model makes it challenging to run a self-converging fit to the data. Visually optimizing the size and position of the two gaps in Fig. 4i, we find a range of parameters (Supplementary Note X) simultaneously reproducing the experimental bias dependent CDW phase and amplitude data remarkably well (Fig. 3). As for the relative phase between the two CMs, it is chosen to minimize the Coulomb interaction of the CMs and to conform with the strong commensuration energy that locally locks them to the lattice. Reducing the Coulomb energy is obtained by introducing a phase shift between the two CMs, which can only be ±120° (Fig. 4g) to satisfy the lock-in criterion with the lattice given the 3a0 periodicity of the CDW.

In the following, we turn to theoretically modelling multiple CDW gaps on different bands in NbSe2. We deploy self-consistent calculations to include the CDW gap within the random phase approximation on the two-dimensional two-band tight-binding fit to the NbSe2 band structure constrained by ARPES3. The corresponding FS shown in Fig. 1b consists of inner (red) and outer (blue) bands originating from symmetric and antisymmetric combinations of the Nb d\({}_{3{z}^{2}-{r}^{2}}\) orbitals. The model (see “Methods”) was previously shown to accurately reproduce the full range of experimental measurements on the charge-ordered state12,28. The resulting DOS for the gapped and ungapped cases in each band are shown in Fig. 5a. To emphasize the DOS reduction accompanying the CDW phase transition, we plot the difference between the gapped and ungapped DOS for each band in Fig. 5b.

Fig. 5: Density of state of the two bands in the self-consistent calculations.

a The red (blue) indicates the band making up the inner (outer) pockets; solid/dashed indicates the ungapped/gapped band structure, where gaps are included self-consistently at the mean-field level. b Difference between the gapped and ungapped cases in a, showing small DOS suppressions at energies away from EF on both bands in addition to the gap at EF on the inner band.

Our theoretical modelling shows a clear gap on the inner band at EF, consistent with the gap measured by ARPES around the K-point3. Interestingly, Fig. 5b reveals further DOS reductions, for example near 100 meV on the inner band and −50 meV on the outer band. These features are indicative of CDW gaps opening away from EF in addition to the (primary) gap at EF, supporting the simple model we propose to understand the bias dependence of the CDW appearance in STM images of NbSe2. According to Fig. 5b, there could even be more than two gaps. Consequently, we have included up to n = 8 gaps to our 1D model. However, the agreement with the data is similar for n = 2 and n = 3 (Supplementary Note XI), and we see no improvements adding more gaps.

In summary, the remarkable match between the bias dependence of the CDW contrast in STM topography and the simple 1D model proposed here provides compelling evidence that the CDW in NbSe2 is composed of at least two out-of-phase CMs on the inner and outer bands. While a 180° phase shift between these two CMs would minimize the Coulomb energy, the complex bias dependence of the CDW amplitude and phase observed by STM can only be reproduced when considering also the commensuration energy. This highlights the importance of the coupling of charge order to the lattice, which manifests in the formation of discommensurations15,16 and ultimately enables the observation of the multiband CMs uncovered here. The present study further highlights the power of topographic imaging to gain unique insight into detailed features of the CDW too faint to be detected accurately by tunneling spectroscopy. The formation of multiple modulations in response to new periodicities of a primary transition directly observed here, is extremely general and should in principle be present in all charge (and spin) density wave materials, and suggests new directions to explore in the physics of spatially modulated electronic orders.


Crystal growth and STM measurements

Single crystals of NbSe2 were grown via iodine-assisted chemical vapour transport and cleaved in-situ at room temperature. STM experiments were done in UHV (base pressure below 2  10−10 mbar) using tips mechanically cut from a PtIr wire and conditioned in-situ on a clean Ag(111) single crystal. The bias voltage was applied to the sample. STM images were recorded in constant current mode with at least 64 pixel/nm resolution. Details of the CDW amplitude and phase fitting procedure can be found in ref 16.

Self-consistent calculations

We carried out a 22-orbital Slater–Koster tight-binding fit29 to the NbSe2 band structure, constrained by ARPES measurements2,3 and local density approximation numerical calculations11. This provided not only the band structure, but the orbital composition of the bands. We found that the two bands crossing the Fermi level are predominantly composed of the niobium \({d}_{3{z}^{2}-{r}^{2}}\) orbitals (>60% across the Brillouin zone). A third, small, pancake-shaped pocket centred on Γ derives primarily from the selenium p-orbitals and is therefore not expected to mix significantly with the other bands, in-keeping with experimental observations that it plays no role in the charge ordering. We then re-fit the two bands of interest using only the two relevant orbitals; the fit was indistinguishable from the phenomenological fit to ARPES data provided in ref. 3.

The Coulomb interaction can be neglected in NbSe2, as the large DOSs at the Fermi level leads to strong screening. This remains true down to the lowest temperatures (above the SC transition at 7.2 K) since the CDW gap only opens on small regions of the FS. The relevant interaction is the electron-phonon coupling, for which we constructed an analytic expression following Ref. 30. It has long been suggested that the CDW in NbSe2 originates not from FS nesting, but from the dependence of the electron–phonon coupling on the momentum transfer in the phonon-mediated electron–electron scattering11. Our calculation includes the dependence of the coupling on the ingoing and outgoing electron momenta, as well as the orbital composition of the electronic bands scattered between. Only by taking all of these factors into account were we able to find a consistent explanation of the full range of experimental observations, including ARPES2,3, scanning tunneling spectroscopy/microscopy4, and inelastic X-ray scattering31. Our model has only one free parameter, the overall magnitude of the electron–phonon coupling, which we fixed using TCDW = 33.5 K.

We modelled the effect of the CDW on the electronic band structure using the random phase approximation. We employed the Nambu–Gor’kov method to work within the gapped phase; this method consists of promoting the electronic Green’s function to a 9 × 9 matrix, representing the tripling of the real-space unit cell in both lattice directions induced by the CDW formation. The CDW gap appears in off-diagonal elements, and diagonalisation then results in a gapped electronic band structure. We solved for the CDW gap self-consistently at high-symmetry points in the Brillouin zone, and used the results to constrain a tight-binding fit for the gap structure. We assumed the gap to be independent of energy. The CDW gap serves as an order parameter, and so our model naturally accounts for long-range CDW correlations. Further details of the method are given in refs. 12,28.

We found that a CDW gap opens at the Fermi level on the outer pocket centred on the K-point, along the MK line, in agreement with ARPES2,3. However, since the order parameter is non-zero at all points in the Brillouin zone, and at all energies, gaps also open wherever band crossings are introduced. This is the origin of the multiband CDW, evidenced by the suppression of DOS at energies below EF seen in Fig. 5. We calculated the DOS at different energies, with and without the CDW gap, by summing over the Brillouin zone the spectral function found from the electronic Green’s function.

Data availability

The data that support the findings of this study has been deposited in the Yareta repository (

Code availability

Computer codes are available upon reasonable request and preferably within a collaboration.

Change history

  • 11 November 2021

    In this article the journal title in the ‘Peer review information’ was incorrectly given as Nature Aging but should have been Nature Communications. The original article has been corrected.


  1. 1.

    Wang, C., Giambattista, B., Slough, C. G., Coleman, R. V. & Subramanian, M. A. Energy gaps measured by scanning tunneling microscopy. Phys. Rev. B 42, 8890 (1990).

    ADS  CAS  Article  Google Scholar 

  2. 2.

    Borisenko, S. V. et al. Two energy gaps and fermi-surface “arcs” in NbSe2. Phys. Rev. Lett. 102, 166402 (2009).

    ADS  CAS  Article  Google Scholar 

  3. 3.

    Rahn, D. J. et al. Gaps and kinks in the electronic structure of the superconductor 2H-NbSe2 from angle-resolved photoemission at 1 K. Phys. Rev. B 85, 224532 (2012).

    ADS  MathSciNet  Article  Google Scholar 

  4. 4.

    Soumyanarayanan, A. et al. Quantum phase transition from triangular to stripe charge order in NbSe2. Proc. Natl Acad. Sci. USA 110, 1623 (2013).

    ADS  CAS  Article  Google Scholar 

  5. 5.

    Arguello, C. J. et al. Visualizing the charge density wave transition in 2H-NbSe2 in real space. Phys. Rev. B 89, 235115 (2014).

    ADS  Article  Google Scholar 

  6. 6.

    Revolinsky, E., Spiering, G. A. & Beerntsen, D. J. Superconductivity in the niobium-selenium system. J. Phys. Chem. Solids 26, 1029 (1965).

    ADS  CAS  Article  Google Scholar 

  7. 7.

    Wilson, J. A., DiSalvo, F. J. & Mahajan, S. Charge-density waves and superlattices in the metallic layered transition metal dichalcogenides. Adv. Phys. 24, 117 (1975).

    ADS  CAS  Article  Google Scholar 

  8. 8.

    Harper, J. M. E., Geballe, T. H. & Di Salvo, F. J. Heat capacity of 2H-NbSe2 at the charge density wave transition. Phys. Lett. A 54, 27 (1975).

    ADS  Article  Google Scholar 

  9. 9.

    Long, J. R., Bowen, S. P. & Lewis, N. E. Anomalous resistivity of iodine-free 2H-NbSe2. Solid State Commun. 22, 363 (1977).

    ADS  CAS  Article  Google Scholar 

  10. 10.

    Moncton, D. E., Axe, J. D. & DiSalvo, F. J. Study of superlattice formation in 2H-NbSe2 and 2H-TaSe2 by neutron scattering. Phys. Rev. Lett. 34, 734 (1975).

    ADS  CAS  Article  Google Scholar 

  11. 11.

    Johannes, M. D., Mazin, I. I. & Howells, C. A. Fermi-surface nesting and the origin of the charge-density wave in NbSe2. Phys. Rev. B 73, 205102 (2006).

    ADS  Article  Google Scholar 

  12. 12.

    Flicker, F. & van Wezel, J. Charge order from orbital-dependent coupling evidenced by NbSe2. Nat. Commun. 6, 7034 (2015).

    ADS  CAS  Article  Google Scholar 

  13. 13.

    Corcoran, R. et al. Quantum oscillations in the mixed state of the type II superconductor 2H-NbSe2. J. Phys.: Condens. Matter 6, 4479 (1994).

    ADS  CAS  Google Scholar 

  14. 14.

    Rossnagel, K. et al. Fermi surface of 2H-NbSe2 and its implications on the charge-density-wave mechanism. Phys. Rev. B 64, 235119 (2001).

    ADS  Article  Google Scholar 

  15. 15.

    McMillan, W. L. Theory of discommensurations and the commensurate-incommensurate charge-density-wave phase transition. Phys. Rev. B 14, 1496 (1976).

    ADS  CAS  Article  Google Scholar 

  16. 16.

    Pásztor, Á. et al. Holographic imaging of the complex charge density wave order parameter. Phys. Rev. Res. 1, 033114 (2019).

    Article  Google Scholar 

  17. 17.

    Mallet, P., Sacks, W., Roditchev, D., Défourneau, D. & Klein, J. Spatial and energy variation of the local density of states in the charge density wave phase of 2H-NbSe2. J. Vac. Sci., Technol. B: Microelectron. Nanometer Struct. 14, 1070 (1996).

    ADS  CAS  Article  Google Scholar 

  18. 18.

    Sacks, W., Roditchev, D. & Klein, J. Voltage-dependent STM image of a charge density wave. Phys. Rev. B 57, 13118 (1998).

    ADS  CAS  Article  Google Scholar 

  19. 19.

    Gao, S. et al. Atomic-scale strain manipulation of a charge density wave. Proc. Natl Acad. Sci. USA 115, 6986 (2018).

    ADS  CAS  Article  Google Scholar 

  20. 20.

    Martinez-Castro, J. et al. Scanning tunneling microscopy of an air sensitive dichalcogenide through an encapsulating layer. Nano Lett. 18, 6696 (2018).

    ADS  CAS  Article  Google Scholar 

  21. 21.

    Gye, G., Oh, E. & Yeom, H. W. Topological landscape of competing charge density waves in 2H-NbSe2. Phys. Rev. Lett. 122, 016403 (2019).

    ADS  CAS  Article  Google Scholar 

  22. 22.

    Guster, B. et al. Coexistence of elastic modulations in the charge density wave state of 2H-NbSe2. Nano Lett. 19, 3027 (2019).

    ADS  CAS  Article  Google Scholar 

  23. 23.

    Fukuyama, H. Pinning in Peierls-Fröhlich state and conductivity. J. Phys. Soc. Jpn. 41, 513 (1976).

    ADS  CAS  Article  Google Scholar 

  24. 24.

    Fukuyama, H. & Lee, P. A. Dynamics of the charge-density wave. I. Impurity pinning in a single chain. Phys. Rev. B 17, 535 (1978).

    ADS  CAS  Article  Google Scholar 

  25. 25.

    Hildebrand, B. et al. Short-range phase coherence and origin of the 1T-TiSe2 charge density wave. Phys. Rev. B 93, 125140 (2016).

    ADS  Article  Google Scholar 

  26. 26.

    Novello, A. M. et al. Stripe and short range order in the charge density wave of 1T-Cux TiSe2. Phys. Rev. Lett. 118, 017002 (2017).

    ADS  CAS  Article  Google Scholar 

  27. 27.

    Spera, M. et al. Insight into the charge density wave gap from contrast inversion in topographic STM images. Phys. Rev. Lett. 125, 267603 (2020).

    ADS  CAS  Article  Google Scholar 

  28. 28.

    Flicker, F. & van Wezel, J. Charge order in NbSe2. Phys. Rev. B 94, 235135 (2016).

    ADS  Article  Google Scholar 

  29. 29.

    Slater, J. C. & Koster, G. F. Simplified LCAO method for the periodic potential problem. Phys. Rev. 94, 1498 (1954).

    ADS  CAS  Article  Google Scholar 

  30. 30.

    Varma, C. M., Blount, E. I., Vashishta, P. & Weber, W. Electron-phonon interactions in transition metals. Phys. Rev. B 19, 6130 (1979).

    ADS  CAS  Article  Google Scholar 

  31. 31.

    Weber, F. et al. Optical phonons and the soft mode in 2H-NbSe2. Phys. Rev. B 87, 245111 (2013).

    ADS  Article  Google Scholar 

  32. 32.

    Pásztor, Á. et al. Dataset Multiband charge density wave exposed in a transition metal dichalcogenide. (2021).

  33. 33.

    Momma, K. & Izumi, F. Vesta 3 for three-dimensional visualization of crystal, volumetric and morphology data. J. Appl. Crystallogr. 44, 1272 (2011).

    CAS  Article  Google Scholar 

Download references


We acknowledge A. F. Morpurgo and I. Maggio-Aprile for inspiring discussions. We thank A. Guipet and G. Manfrini for their technical support in the STM laboratories. This project was supported by Div.2 (Grant No. 182652) and Sinergia (Grant No. 147607) of the Swiss National Science Foundation. F. F. acknowledges support from the Astor Junior Research Fellowship of New College, Oxford.

Author information




C.R. designed the experiment. M.S. and A.S. took care of the STM experiments. Á.P. performed the data analysis. F.F. and J.v.W. performed the self-consistent calculations. Á.P. did the one-dimensional model simulations. C.B. and E.G. synthesized the bulk crystals. Á.P., A.S., F.F., J.v.W. and C.R. wrote the paper. All authors contributed to the scientific discussions and manuscript revisions.

Corresponding authors

Correspondence to Árpád Pásztor or Christoph Renner.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review informationNature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. This article has been peer reviewed as part of Springer Nature’s Guided Open Access initiative.

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

Supplementary information

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

Pásztor, Á., Scarfato, A., Spera, M. et al. Multiband charge density wave exposed in a transition metal dichalcogenide. Nat Commun 12, 6037 (2021).

Download citation


By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.


Quick links