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.

A full gap above the Fermi level: the charge density wave of monolayer VS2


In the standard model of charge density wave (CDW) transitions, the displacement along a single phonon mode lowers the total electronic energy by creating a gap at the Fermi level, making the CDW a metal–insulator transition. Here, using scanning tunneling microscopy and spectroscopy and ab initio calculations, we show that VS2 realizes a CDW which stands out of this standard model. There is a full CDW gap residing in the unoccupied states of monolayer VS2. At the Fermi level, the CDW induces a topological metal-metal (Lifshitz) transition. Non-linear coupling of transverse and longitudinal phonons is essential for the formation of the CDW and the full gap above the Fermi level. Additionally, x-ray magnetic circular dichroism reveals the absence of net magnetization in this phase, pointing to coexisting charge and spin density waves in the ground state.


The many-body ground states of two-dimensional (2D) materials, wherein the reduced dimensionality leads to the enhancement of correlation effects, have been extensively researched in recent years. Of particular interest are the coexistence or competition between charge density waves (CDWs), as found in many 2D transition metal dichalcogenides (TMDCs), with superconducting and magnetic phases1,2. Since these phases can be strongly dependent on the substrate3,4 or the defect density5,6, the intrinsic properties of 2D materials are difficult to determine experimentally. In addition, CDWs themselves are the subject of an ongoing controversy regarding the driving force behind the CDW transition and the exact structure of the electronic system in the CDW phase of 2D materials7,8.

Peierls’ explanation for the CDW in a one-dimensional chain of atoms states that periodic lattice distortions open an electronic gap at the nesting wavevector. This gap at the Fermi level lowers the energy of the occupied states and thus the total energy, while increasing the energy of the unoccupied states that do not contribute to the total energy. Thus, this gapping mechanism requires the gap to be at the Fermi level. However, in many (quasi-)2D cases, CDWs form in the complete or partial absence of Fermi-surface nesting, suggesting that the driving mechanism behind their formation lies beyond a simple electronic disturbance9, and it has been questioned whether the concept of nesting is essential for understanding CDW formation10,11,12. Instead, a strong and wavevector-dependent electron–phonon coupling is often predicted to be the driving force behind the transition7. For these CDWs, spectral reconstructions are not limited to a small energy window around the Fermi energy, but can occur throughout the entire electronic structure, opening the door to novel spectral fingerprints of the CDW. A full gap could occur away from the Fermi energy. However, even for the well-studied strong-coupling TMDCs 2H-NbSe29,13,14 and 1T-TaS27,15,16,17, no experimental verification of a clear CDW gap located away from the Fermi energy has been provided to date. Furthermore, at the Fermi energy, the undistorted phase and the CDW can have different Fermi-surface topologies, with the implication that the transition is a metal–metal Lifshitz transition18.

Metallic 1T-VS2 is not only a promising electrode material in lithium-ion batteries19,20, but also a prototypical d1 system, expected to host strongly correlated physics21. It is stated to be a CDW material22,23 and a candidate for 2D magnetism24,25 with layer-dependent properties26, making it a model system for investigating complex ground states. Although difficult to synthesize, bulk 1T-VS2 has been well studied, with many authors finding a CDW transition at around 305 K when it was prepared via the de-intercalation of Li22,23,27,28,29. However, recent powder samples prepared under high pressure show no CDW transition30. Based on their finding of a phonon instability at \(2/3\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\) corresponding to the experimental CDW wavevector of Li de-intercalated bulk samples23, Gauzzi et al. point out that bulk “VS2 is at the verge of CDW transition”30 but not a CDW material. Due to a similar difficulty in synthesis, the properties of monolayer 1T-VS2 have proven equally elusive31. Theoretical calculations had predicted ferromagnetism and a CDW with a wavevector of \(2/3\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\)21,25. When it was first synthesized however, scanning tunneling microscope (STM) measurements did not reveal a CDW31, presumably due to strong hybridization with the Au(111) substrate, similar to the case of 2H-TaS2 on Au(111)4,32,33.

Here we report the growth of VS2 monolayers on the inert substrate graphene (Gr) on Ir(111) via a two-step molecular beam epitaxy (MBE) synthesis developed for sulfur-based TMDCs34. Using a combination of STM, scanning tunneling spectroscopy (STS), and ab initio density functional theory (DFT) calculations, we determine the spatial and electronic structure of monolayer VS2. We observe a \({{{{{{{\bf{q}}}}}}}}\approx 2/3\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\) CDW as the electronic ground state at 7 K, which remains stable up to room temperature. A full gap in the density of states (DOS), residing completely in the unoccupied states, is measured via STS. From DFT and density functional perturbation theory (DFPT), we find that, although a transverse phonon mode initially becomes unstable in the harmonic approximation, the final CDW has a substantial admixture of longitudinal modes. The calculations are in excellent agreement with experiment, regarding both the electronic structure of the CDW phase and the spatial charge distribution observed on the VS2 islands.

X-ray magnetic circular dichroism (XMCD) measurements at 7 K and 9 T robustly show vanishing total net magnetization. The coupling of the CDW to a spin density wave (SDW), energetically favored in DFT calculations, could explain this observation, providing interesting prospects for future research on the interplay of CDWs and magnetism.


CDW in monolayer VS2

The typical morphology of the MBE-grown monolayer VS2 islands on Gr/Ir(111) is shown in the large-scale STM image in Fig. 1a. The islands were grown by room-temperature deposition of vanadium in a sulfur background pressure of \({P}_{{{{{{{{\rm{S}}}}}}}}}^{{{{{{{{\rm{g}}}}}}}}}\,=\,1\,{\times }\,10^{-8}\ {{{{{{{\rm{mbar}}}}}}}}\) and subsequently annealed at 600 K in the same sulfur pressure. Annealing to temperatures of 800 K and above leads to the formation of a variety of sulfur-depleted phases, which are not under concern here. Similar observations were made by Arnold et al.31, who established monolayer stoichiometric 1T-VS2 on Au(111) by annealing in a sulfiding gas at 670–700 K, while sulfur-depleted monolayer phases form when annealed to the same or higher temperature in the absence of sulfiding species. We also note that depending on growth temperature and sulfur pressure bilayer samples without any monolayer islands evolve.

Fig. 1: Structure of VS2 on Gr/Ir(111) at 7 and 300 K.
figure 1

a Large-scale 7 K STM topograph of monolayer VS2 islands with small bilayers present. A height profile along the horizontal black line is shown below the image. b, c Atomically resolved STM images of monolayer VS2 at 7 K (b) and 300 K (c). The Fourier transform of each image is shown as an inset, with the 1 × 1 VS2 structure in red and the superstructure spots indicated in blue. In b, an atomic model for the \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) superstructure is included as an overlay. The model depicts the top sulfur atoms, with their apparent height in STM coded in orange (low) and yellow (high). Measurement parameters: a 80 × 50 nm2, It = 0.8 nA, Vt = −800 meV, b 6 × 6 nm2, It = 0.6 nA, Vt = 400 meV, c 6 × 6 nm2, It = 1.0 nA, Vt = −1000 meV.

The monolayer islands are fully covered by a striped superstructure which is present regardless of island size or defect density and occurs in domains, typically separated by grain boundaries. In the topograph of Fig. 1b, taken at 7 K, the VS2 lattice is resolved, exhibiting the hexagonal arrangement of top layer sulfur atoms as protrusions. We find that monolayer VS2 has a lattice constant of \({a}_{{{{{{{{\rm{VS}}}}}}}}_{2}}={(3.21\,\pm \,0.02)}\,\ \,\,\)Å, in good agreement with the bulk lattice constant of 3.22 Å of 1T-VS227,35. The similarity of the lattice constants indicates also the absence of epitaxial strain, consistent with the random orientation of the VS2 with respect to the Gr.

The stripes of the superstructure have an average periodicity of \((2.28\pm 0.02){a}_{{{{{{{{\rm{VS}}}}}}}}_{2}}\). Close analogues to this structure have previously been observed in stoichiometric monolayer VSe2. There, a superstructure of identical symmetry is attributed to a CDW3,36,37,38 [compare Supplementary Fig. 1]. The superstructure is found to persist up to room temperature, as can be concluded from the STM topograph in Fig. 1c, taken at 300 K. At this temperature, the superstructure appears spontaneously only on larger islands, suggesting that the transition temperature between the superstructure and the undistorted phase is not far above room temperature. Indeed, on smaller islands the STM tip can be used to reversibly switch between the undistorted (1 × 1) structure and the superstructure, shown in Supplementary Fig. 2. This directly excludes the possibility that the superstructure is due to a sulfur-depleted phase. We conclude that the superstructure is most likely a CDW in a stoichiometric monolayer of 1T-VS2.

For the DFT calculations below, the experimental wave pattern must be approximated by a commensurate structure. A close approximation with periodicity \(2.25{a}_{{{{{{\rm{VS}}}}}}_{2}}\) is overlaid on the atomic resolution image in Fig. 1b. It locally matches the incommensurate CDW quite well. The blue box indicates the corresponding \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) unit cell. The Fourier transform of the topograph is shown as inset in Fig. 1b, with the wavevector of the CDW indicated (blue arrow). The same is done for the 300 K topograph in Fig. 1c. Within the margin of error, the wavevector is found to be temperature independent, with \({{{{{{{{\bf{q}}}}}}}}}_{{{{{\rm{CDW}}}}}({{{{\rm{7}}}}}{{{{\rm{\space}}}}}{{{{\rm{K}}}}})}=(0.656\pm 0.006)\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\) and \({{{{{{{{\bf{q}}}}}}}}}_{{{{{\rm{CDW}}}}}({{{{\rm{300}}}}}{{{{\rm{\space}}}}}{{{{\rm{K}}}}})}=(0.65\pm 0.03)\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\). Since the wavevector of the \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) unit cell, \({{{{{{{{\bf{q}}}}}}}}}_{{{\mbox{}}}9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }}= 2/3\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\approx 0.667\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\), is slightly larger than the experimental value, we will in the following also consider another unit cell of size \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\), with a slightly smaller wavevector \({{{{{{{{\bf{q}}}}}}}}}_{{{\mbox{}}}7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }{{\mbox{}}}}=9/14\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\approx 0.643\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\). With the experimental wavevector lying in between \({{{{{{{{\bf{q}}}}}}}}}_{{{\mbox{}}}7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }{{\mbox{}}}}\) and \({{{{{{{{\bf{q}}}}}}}}}_{{{\mbox{}}}9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }{{\mbox{}}}}\), calculations with these two unit cells should capture the essential features of the incommensurate structure and provide a check on any artefacts or errors arising from using them for computational purposes (cf. Supplementary Fig. 3).

Energetics of lattice instabilities

Ab initio DFPT calculations of the acoustic phonon dispersion of undistorted monolayer 1T-VS2 confirm that a structural instability and corresponding tendencies toward CDW formation exist for the experimental wavevector. Figure 2a shows that the longitudinal–acoustic and transverse–acoustic modes feature imaginary frequencies in several parts of the Brillouin zone. In other words, the Born–Oppenheimer energy surface is a downwards-opening parabola for small atomic displacements in the direction of these modes, as seen in Fig. 2b (triangle marks). At the experimental wavevector between \({{{{{{{\bf{q}}}}}}}}=2/3\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\) and \({{{{{{{\bf{q}}}}}}}}=9/14\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\), we find an instability of the transverse–acoustic branch. However, the dominant instability within the harmonic approximation (i.e., DFPT), is located at \({{{{{{{\bf{q}}}}}}}}=1/2\ \overline{{{\Gamma }}{{{{{{{\rm{M}}}}}}}}}\) in the longitudinal–acoustic branch.

Fig. 2: Lattice instabilities in monolayer 1T-VS2 from first principles.
figure 2

a Acoustic phonon dispersion from DFPT. LA, TA, and ZA stand for dominant longitudinal, transverse, and out-of-plane atomic displacements. The insets show selected displacement patterns corresponding to indicated modes. b Total energy from DFT as a function of the displacement amplitude for atomic displacements toward relaxed crystal structures and their projections onto soft phonon modes. ce Relaxed crystal structures on \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\), \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\), and 4 × 4 unit cells from DFT. Vanadium and sulfur atoms are represented by black and yellow dots, their undistorted positions by gray shadows. Arrows represent the projections of the atomic displacements onto soft phonon modes. (Only arrows longer than 2% of the lattice constant are shown.) The contributions of different phonon modes are quantified in the figure titles. The displacements in c, d are drawn to scale, those in e have been magnified by a factor of three for better visibility. The primitive cell of the structure in c, which is in agreement with the results of ref. 21, is outlined in yellow. Dashed lines and crosses mark reflection planes and inversion centers.

To go beyond the harmonic approximation, we have performed structural relaxations on appropriate unit cells. The resulting atomic positions are shown in Fig. 2c–e. On the aforementioned \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) and \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) unit cells, which can approximately host an integer multiple of the observed wavelength, the vanadium atoms are displaced from their symmetric positions by up to 8% of the lattice constant, while the positions of the sulfur atoms remain almost unchanged, see Fig. 2c, d. The associated energy gains amount to about 23 meV per VS2 formula unit (cf. ref. 21). The magnitude of these distortions and energy gains is similar to other octahedral TMDCs but exceeds by far what is found in trigonal–prismatic TMDCs7,39. For instance, on the DFT level, the maximum displacement in the \(\sqrt{13}\times \sqrt{13}\) CDW of 1T-NbSe2 is 8.8% of the lattice constant with an energy gain of 57 meV per formula unit40, while in the 3 × 3 CDW of 2H-NbSe2 distortions and energy gain amount to only 2.3% of the lattice constant and 3.7 meV per formula unit41.

The vanadium displacement has components in both the transverse and longitudinal direction (vertical and horizontal in Fig. 2c, d), even though the instability of the phonons at \({{{{{{{\bf{q}}}}}}}}=2/3\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\) and \({{{{{{{\bf{q}}}}}}}}=9/14\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\) is of transverse character (white arrows in Fig. 2c, d). As a consequence, all longitudinal displacement components must stem from non-linear mode–mode coupling beyond the harmonic approximation. The admixture of longitudinal displacement components stems mainly from wavevectors \({{{{{{{\bf{q}}}}}}}}=4/3\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\) and \({{{{{{{\bf{q}}}}}}}}=9/7\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\), which are also commensurate with the \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) and the \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) unit cells, respectively. The admixed longitudinal modes at \({{{{{{{\bf{q}}}}}}}}=4/3\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\) and \({{{{{{{\bf{q}}}}}}}}=9/7\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\) are stable in the harmonic approximation and the non-linear admixture is not related to any nesting or Peierls physics (cf. Supplementary Fig. 4e, f).

We also find a distorted ground state on a 4 × 4 unit cell, see Fig. 2e. This structure is commensurate to the six wavevectors \({{{{{{{\bf{q}}}}}}}}=1/2\ \overline{{{\Gamma }}{{{{{{{\rm{M}}}}}}}}}\), where we have instabilities in the longitudinal–acoustic branch arising from near perfect Fermi-surface nesting, see Supplementary Fig. 4a. However, here the displacements amount to only 4% of the lattice constant with a corresponding energy gain below 4 meV per 1T-VS2 formula unit—much less than what is found for the \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) or \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) CDW structures. Thus, the DFT total energies of the fully relaxed structures are in line with the experimentally observed CDW patterns.

To illustrate the significance of the non-linear mode–mode coupling, in Fig. 2b, we also show the Born–Oppenheimer energy surfaces for displacements toward the relaxed structures (circle marks). The energy curve of the 4 × 4 structure is steeper in the vicinity of the origin. In other words, the 4 × 4 structure wins for small displacements. However, for larger displacements, the structures corresponding to the experimental wavevector reach by far the lowest values. These large energy gains at large displacements are inaccessible without non-linear mode–mode coupling, i.e., without the contribution of stable phonon modes (triangle marks). In the next section, we will address the non-linear regime of the distortions in a quantitative manner.

Non-linear mode–mode coupling

We decompose the entirety of atomic displacements of the relaxed \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) structure as u + v, where u points in the direction of the unstable transverse–acoustic phonon modes at \({{{{{{{\bf{q}}}}}}}}=\pm 9/14\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\) and the orthogonal complement vu combines contributions from all other phonon modes. The unstable modes account for u2/u + v2 ≈ 69% of the total displacement only. In Fig. 2b, we have already seen one-dimensional cross sections of the Born–Oppenheimer energy surface, E(αu) and E(β(u + v)), where α and β are dimensionless scaling factors. Now, we will consider the full 2D Born–Oppenheimer surface spanned by u and v. Figure 3 shows E(xu + yv), where the minimum at x = y = 1 corresponds to the \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) structure and x = y = 0 is the undistorted structure. A fourth-order polynomial fit,

$${\frac{E(x{{{{{{{\bf{u}}}}}}}}+y{{{{{{{\bf{v}}}}}}}})}{{{{{{{\rm{meV/VS}}}}}}}}}_{\!\!\! 2}\approx \, {-25}{x}^{2}+29{y}^{2}+34{x}^{3}-99{x}^{2}y-20x{y}^{2}-12{y}^{3}\\ +0.1{x}^{4}+44{x}^{3}y+13{x}^{2}{y}^{2}+4.8x{y}^{3}+7.2{y}^{4},$$

accurately describes the DFT Born Oppenheimer surface. Here, the first and second line give rise to linear and non-linear forces F = −E, respectively. It turns out that the non-linear part of the forces is dominated by mode–mode coupled terms42,43,44 (dependent on both x and y). The energy reduction stems largely from the x2y and xy2 terms above, which correspond to a shift of the minimum of the potential-energy surface toward finite positive y upon finite displacement in x direction and a softening of the effective spring constant in y direction for finite positive x, respectively. Note that within the harmonic approximation the x2 (y2) term lowers (raises) the energy.

Fig. 3: Born–Oppenheimer energy surface for the \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) structure of 1T-VS2.
figure 3

The axes represent the projection of the full CDW displacement onto the transverse–acoustic (TA) phonon modes at \({{{{{{{\bf{q}}}}}}}}=\pm 9/14\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\) and the orthogonal complement, which combines all other contributing modes. The full CDW displacement is located at the point (1, 1). The forces resulting from this energy surface are predominantly non-linear and coupled in both directions.

The decisive role of mode–mode coupling terms x2y and xy2 distinguishes 1T-VS2 from systems like 2H-NbSe2 or 2H-TaS2, where a single mode can be employed to describe anharmonicities, and distortions along a single effective coordinate suffice to explain the relaxation pattern of the full CDW and associated energy gains (cf. Supplementary Fig. 5).

The non-linear mode–mode coupling also manifests in monolayer 1T-VTe2, which is isoelectronic to monolayer 1T-VS2. Monolayer 1T-VTe2 in experiment realizes a 4 × 4 CDW45 in contrast to monolayer 1T-VS2. In line with experiment, the comparison of DFT total energies in the fully relaxed supercells (Supplementary Table 1) reveals a clear preference of the 4 × 4 structure in 1T-VTe2. At the harmonic level, this is likely related to a shift of the lattice instabilities, especially in the transverse-acoustic branch, toward smaller wavevectors in 1T-VTe2 as compared to 1T-VS2 (cf. Fig. 2a and Supplementary Fig. 6a), which can be traced back to differences in the Fermi surface (cf. Supplementary Figs. 4 and 6b–g). At the harmonic level, a CDW with \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) structure of monolayer 1T-VS2 is not expected, as Supplementary Fig. 6a shows. The small energy gain and still appreciable distortions obtained from the relaxation of a \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) structure of monolayer 1T-VS2 (Supplementary Table 1) despite the stability on the harmonic level suggest that non-linear mode–mode coupling is also effective, here.

Full CDW gap in the unoccupied states

To better understand this CDW phase, we determined the electronic structure of monolayer VS2 by a combination of STS experiments and simulated dI/dV maps based on the ab initio calculations using the \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) and \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) unit cells. STS spectra were used to locally probe the DOS of monolayer VS2 at 7 K (black line) and 78.5 K (purple line), shown in Fig. 4a. Both spectra were taken with a clean Au tip in the middle of VS2 islands. The most prominent feature is the gap located at about 0.175 eV, which is absent in calculations of undistorted monolayer VS221. At 7 K, the dI/dV signal vanishes completely, corresponding to a full gap in the DOS. At 78.5 K, this gap is not fully open, appearing as a wide depression with a finite value at its minimum. In most other characteristic features the spectra agree qualitatively.

Fig. 4: Spatially and electronically resolved CDW phase in monolayer VS2.
figure 4

a Scanning tunneling spectroscopy (STS) dI/dV spectra taken with a Au tip on monolayer VS2 at 78.5 K (purple) and 7 K (black). The spectra are plotted along with the DFT-calculated density of states (DOS) for the \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) (cyan) and \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) (indigo) CDW phases of monolayer 1T-VS2. b Atomically resolved STM topograph of monolayer VS2 taken at Vt = 175 meV. c, d Fourier-filtered dI/dV conductance maps of the same region as in b, taken at Vt = 75 meV (c) and Vt = 275 meV (d). A linear yellow (maximum) to blue (minimum) color scale is used to depict the dI/dV intensity. The blue box indicates the same location in bd and corresponds to a single \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) unit cell of the CDW. In the same color scale, DFT-simulated dI/dV maps below (c) and above (d) the gap of the charge density wave (CDW) are overlaid as insets. The maps show the integrated DOS from 0 to 137 meV (c) and from 137 to 275 meV (d). Additionally, the Fourier transforms of the conductance maps are shown in the upper-left corners with the 1 × 1 (red) and CDW peaks (blue) highlighted by circles. Measurement parameters: f = 777.7 Hz, a T = 78.5 K, It = 0.3 nA, Vr.m.s. = 6 meV and T = 7 K, It = 0.45 nA, Vr.m.s. = 4 meV, bd T = 7 K, 5.5 × 5.5 nm2, It = 0.3 nA, Vr.m.s. = 10 meV.

While the lack of energy resolution at 78.5 K certainly smears out the spectra and the gap, the reason for its absence is not immediately evident. When discussing the band structure below, it will be seen that the width and existence of the full gap depend on the magnitude of the lattice distortions, which may already be diminished at 78.5 K.

In the same figure, ab initio calculations for the DOS of VS2, structurally relaxed in the \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) (cyan) or \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) (indigo) unit cell, are shown. Both unit cells feature quite similar structures, as expected for close-lying q vectors. Most striking, for both cases a full gap in the unoccupied states is predicted. They only differ in size: 0.13 eV and 0.21 eV for the \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) and \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) unit cell, respectively. The location of the gap matches the STS data. That the width of the gap in the spectrum is smaller than in DFT might stem from the ground-state calculation assumed in DFT, overestimating the vanadium atom displacement at realistic temperatures. Note also that while many of the characteristic features of calculated and measured DOS (peaks, minima) seem to agree, the experimental spectra appear to be compressed with respect to the DFT calculated DOS. This quasiparticle renormalization is indicative of strong electron–electron correlations beyond the approximations of DFT (compare Supplementary Fig. 7).

With theory and experiment largely agreeing on the electronic structure, we turn to the relation between the gap and the CDW measured on the VS2 islands. For that purpose, dI/dV conductance maps were taken on either side of the gap (both in the unoccupied states), in the location shown in Fig. 4b. The maps help to distinguish structural from electronic contributions, providing a close approximation of the spatial distribution of the DOS at the selected energies. As shown in Fig. 4c, d, we find two different DOS distributions on either side of the gap (see Supplementary Fig. 8 for the in-gap DOS). Both distributions are locked into the distorted lattice periodicity. They are out-of-phase, as seen in the blue unit cell drawn in the same location in Fig. 4b–d: The DOS maxima below the gap correspond to DOS minima above the gap and vice versa. This behavior is perfectly analogous to that for a CDW with a symmetric gap around the Fermi level4. Simulated dI/dV maps derived from the DFT DOS for a \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) CDW are shown as an overlay in Fig. 4c, d. In Fig. 4c, the simulation reproduces both the alternating rows of single and zigzag atoms and the DOS minima between the rows. Its counterpart in Fig. 4d shows higher DOS contrast than experiment, but presents the same qualitative features. With the simulated maps based on the displacement patterns of Fig. 2c, the close agreement with experiment emphasizes the need to look beyond the harmonic approximation to understand this type of CDW.

Band structure and Fermi-surface topology

To deepen our understanding of the system, we calculated the spin-degenerate band structure, density of states, and Fermi surface of monolayer 1T-VS2 with DFT. The results are shown in Fig. 5 and the Supplementary Movies. In the undistorted case, we find a single electronic band at the Fermi level, which strongly disperses between M and K and features a Van Hove singularity in the unoccupied states, as shown in Fig. 5a. The Fermi surface, depicted in Fig. 5b, consists of cigar-shaped electron pockets around the M points. For small distortions, partial gaps open at the Fermi level (e.g., between M and K). With increasing amplitude of the distortion, the gaps become larger and the bands are heavily reconstructed also for high energies. Only then, a full gap as observed in STS at 7 K emerges (cf. Supplementary Movie 1).

Fig. 5: Electronic structure of monolayer 1T-VS2 from DFT.
figure 5

Data for the undistorted structure and the \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) CDW is shown in black and blue, respectively. The CDW data has been unfolded to the Brillouin zone of the undistorted structure. Here, the linewidth/saturation corresponds to the overlap of CDW and undistorted wave functions for the same k point. Analogous results for 1T-VSe2 and the \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) CDW are shown in Supplementary Fig. 3. a Electronic band structure and density of states (DOS) for 0 %, 1/3, 2/3, and 100 % of the displacements of the relaxed CDW structure. Please note that since the CDW breaks the C3 symmetry, the chosen path, indicated in b, does not represent the full Brillouin zone. The bands along an extended path are shown in Supplementary Fig. 9. b, c Lifshitz transition. The Fermi surface is shown for displacements toward the relaxed CDW structure (b) and its projection onto unstable transverse-acoustic (TA) phonon modes (c). The Supplementary Movies 1 and 2 show animations of the transitions in b and c, respectively (including bands, DOS, and structures).

The presence of the CDW is therefore in the first place correlated with the gap between M and K, which opens already for small displacements and results in a partial gapping of the total DOS. Presumably it is the associated gain in electronic energy that initially drives the CDW transition. Since the full gap only starts to open at 70% of the final displacements, the experimental observation of a full gap above the Fermi level is an indication that the displacements in the experiments do not fall much below the calculated ones.

At the Fermi level, there is no complete gap even at large distortion, since the downwards-dispersing bands along Γ–M are only slightly shifted downwards and remain above the Fermi level near Γ. On the other hand, the originally flat portion of the band structure between Γ and K now disperses downwards and crosses the Fermi level. The preservation of states near Γ that mask the partial gap at the Fermi level can be understood in terms of band characters and degeneracies, as shown in Supplementary Fig. 10. Altogether, the Fermi surface is reconstructed and not completely destroyed by the lattice distortion. The CDW transition is thus a metal–metal Lifshitz transition with a change in Fermi-surface topology, instead of the usual metal–insulator Peierls transition.

As shown in Fig. 5c, we cannot understand this Fermi-surface reconstruction based on a single unstable mode: The displacements expected from the harmonic approximation (u in Eq. (1)) only induce gaps in two-third of the cigar-shaped electron pockets (cf. Supplementary Movie 2). The other component v couples to segments of the Fermi surface that are not affected by u, i.e., the remaining third of the electron pockets (cf. Supplementary Movie 3). Together, they transform the Fermi surface from multiple cigar-shaped electron pockets around the M points to the single elliptical hole pocket around Γ visible in Fig. 5b. The decomposition of the CDW contains modes at more than one wavevector q, so several approximate Fermi-surface nesting conditions and electron–phonon coupling matrix elements play a role (Supplementary Fig. 4b, c, e, f), enabling the CDW to affect distinct parts of the Fermi surface.

Magnetic properties of monolayer VS2

Prompted by the prediction of ferromagnetism for monolayer 1T-VS2 in its \({{{{{{{\bf{q}}}}}}}}=2/3\ \overline{{{\Gamma }}{{{{{{{\rm{K}}}}}}}}}\) CDW phase21, we also examined the magnetic properties of VS2, by means of X-ray magnetic circular dichroism (XMCD). The monolayer VS2 samples were grown in situ and investigated with STM beforehand to make sure that the same phase and decent coverage were obtained. A STM topograph of the sample investigated by XMCD is shown in Supplementary Fig. 11. The blue curve in Fig. 6a represents the X-ray absorption spectrum averaged over both helicities and external field directions. The overall line shape is very similar to previous bulk crystal measurements23 and clearly fits to a 3d1 configuration46. The red signal in Fig. 6a is the XMCD magnified by a factor of 10, where no signal above the noise level is visible. This implies that the total magnetization vanishes. Sum rule analysis would yield an upper bound of 0.02μB per vanadium atom. Since it cannot be strictly applied to the case of the V2,3 edges47, this analysis yields only a zero-order estimate of the upper bound, but we can safely conclude that neither ferromagnetic nor paramagnetic behavior is present in this system.

Fig. 6: Magnetic properties of monolayer 1T-VS2.
figure 6

a Plotted in blue is the X-ray absorption signal averaged over both helicities and directions of the B field. The corresponding XMCD is shown in red. All measurements have been conducted in B fields of ±9 T and at a temperature of 7 K. b, c Possible SDW pattern with (b) and without CDW (c).

We investigated magnetic order in monolayer VS2 using spin-polarized DFT. We were able to stabilize both ferromagnetic and SDW structures within the \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) unit cell. In fact, magnetically ordered CDW phases are preferred over nonmagnetic CDW phases by energies of the order of 1 meV per VS2 unit. Figure 6b shows the most favorable SDW pattern in the \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) CDW phase. The magnetic moments on vanadium reach ±0.18μB, those on sulfur only ±0.01μB and are thus not shown. While the CDW alone reduces the total energy to −22.7 meV per VS2 unit with respect to the symmetric structure, the SDW lowers this value by another 1.5 to −24.2 meV. Interestingly, without the CDW, a similar SDW with larger local moments of up to ±0.51μB (shown in Fig. 6c) leads to an energy reduction of 7.1 meV. As already suggested by previous calculations of ferromagnetism in the \(9\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) structure of 1T-VS221, there is a competition between the lattice distortion and the formation of local moments. Although a full account of magnetism needs to go beyond the DFT level, in view of the good agreement between our ab initio results and the experimental STS and XMCD data, the formation of coupled CDW–SDW state in 1T-VS2 is plausible. This presumption is further supported by comparison of the calculated DOS in the CDW–SDW state to the STS shown in Supplementary Fig. 12: the SDW formation on top of the CDW leads to a reduction of the gap size and the DOS of the coupled CDW–SDW is even in better agreement with the experiment than the non-spin-polarized CDW DOS.


Both the electronic and magnetic results for VS2 shed some light on the properties of the isoelectronic compound VSe2, which displays a CDW of the same periodicity3,38,48. Our calculations strongly suggest that also for this system non-linear effects are relevant and that a full gap opens in the unoccupied states (compare Supplementary Fig. 3). A full gap at the Fermi level, as proposed for 1T-VSe26,38,48,49, would be unlikely to intrinsically occur for a CDW with the observed wavevector. The strong similarity between our calculations and experimental data, especially for those VSe2 systems where only the \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) CDW is observed36,37, lends credence to our analysis (compare Supplementary Fig. 13). It is possible that the simultaneous occurrence of a 4 × 1 CDW3,6, perhaps due to substrate-induced strain50, causes an additional gap opening near the Fermi level as a result of the interplay between the CDWs. In any case, similar to VS2, the presence of a SDW coupled to the \(7\times \sqrt{3}{{{{{{{\rm{R}}}}}}}}3{0}^{\circ }\) CDW could explain the absence of net magnetization in XMCD experiments37,38,48. Spin-polarized STM or XMLD might be able to detect the magnetic ground state for both VS2 and VSe2.

In conclusion, VS2 defies the common phenomenology of CDW formation, as the complete CDW gap occurs above the Fermi level, there is giant non-linear longitudinal–transverse mode–mode coupling, and the CDW formation is accompanied by a change of the Fermi-surface topology. The unconventional CDW appears to host further electronic correlations as signalled by the quasiparticle renormalization and magnetic-moment formation. In this respect, it is reminiscent of correlated phases in superlattice structures such as Star of David phases7, moiré superlattices51, and doped cuprate superconductors52. In the latter class, lattice anharmonicities are central to boosting superconductivity under THz optical driving44. The case of VS2 presents new terrain: A metal–metal Lifshitz transition from non-linear electron–lattice effects in the strong-coupling regime is intertwined with electronic correlations. We note that the full gap in the DOS, situated within 0.2 eV from the Fermi level, opens up the possibility of inducing a metal–insulator transition upon mild gating or doping (e.g., with Li). Finally, we are convinced that the excellent agreement of experiment and theory for the unconventional CDW of monolayer VS2 with the full gap in the unoccupied states provides a paradigmatic case study of strong-coupling CDWs in general.


The Ir(111) crystal is cleaned by grazing incidence 1.5 keV Ar+ ion exposure and flash annealing to 1500 K. A closed monolayer of single-crystalline Gr on Ir(111) is grown by room temperature exposure of Ir(111) to ethylene until saturation, subsequent annealing to 1300 K, followed by exposure to 200 L ethylene at 1300 K53.

The synthesis of vanadium sulfides on Gr/Ir(111) is based on a two-step MBE approach introduced in detail in ref. 34 for MoS2. In the first step, the sample is held at room temperature and vanadium is evaporated at a rate of FV = 2.5 × 1016 atoms/(m2s) into a sulfur background pressure of \({P}_{{{{{{{{\rm{S}}}}}}}}}^{{{{{{{{\rm{g}}}}}}}}}=1\times 1{0}^{-8}\ {{{{{{{\rm{mbar}}}}}}}}\) built up by thermal decomposition of pyrite inside a Knudsen cell. This results in dendritic TMDC islands of poor epitaxy. To make the islands larger and more compact, the sample is flashed in a sulfur background to 600 K.

The VS2 layers were analyzed by STM, STS, and low-energy electron diffraction (LEED) inside a variable temperature (30–700 K) ultrahigh vacuum apparatus and a low-temperature STM operating at 7 and 78.5 K. The software WSXM54 was used for STM data processing. XMCD measurements have been conducted at the beamline ID32 of the European Synchrotron Radiation Facility (ESRF) in Grenoble, France. The VS2 samples were grown in situ inside the preparation chamber and checked with LEED and STM before X-ray absorption spectroscopy measurements. To be surface sensitive, the measurements were conducted in the total-electron-yield mode under normal incidence. The measurement temperature was 7 K and fields of 9 T were used. The spectra were recorded at the L3,2 edges, i.e., using the dipole allowed transition from 2p states into the 3d shell potentially generating magnetism.

All DFT and DFPT calculations were performed using QUANTUM ESPRESSO55,56. We apply the PBE functional57,58 and norm-conserving pseudopotentials from the PSEUDODOJO table59,60. In the undistorted case, uniform meshes (including Γ) of 12 × 12 q and 24 × 24 k points are combined with a Fermi–Dirac smearing of 300 K. For a fixed unit-cell height of 15 Å, minimizing forces and in-plane pressure to below 1 × 10−5 Ry/Bohr and 0.1 kbar yields a lattice constant of 3.18 Å and a layer height (vertical sulfur–sulfur distance) of 2.93 Å. For the superstructure calculations, appropriate k-point meshes of similar density are chosen, except for the precise total energies quoted in the section about magnetism and in Supplementary Table 1, which required four times as dense meshes. The average lattice constant of superstructures is kept fixed at the value of the symmetric structure. Fourier interpolation to higher k resolutions (1000 × 1000 for calculations of the DOS) and the unfolding of electronic states is based on localized representations generated with WANNIER9061. For the visualization of the unfolded Fermi surfaces, a Fermi–Dirac broadening of 10 meV is used.

Data availability

All the data and methods are present in the main text and the supplementary materials. Any other relevant data are available from the authors upon reasonable request.

Code availability

Codes used in this work are available from the authors upon reasonable request.


  1. 1.

    Ugeda, M. M. et al. Characterization of collective ground states in single-layer NbSe2. Nat. Phys. 12, 92 (2016).

    CAS  Google Scholar 

  2. 2.

    Zheng, B.-X. et al. Stripe order in the underdoped region of the two-dimensional Hubbard model. Science 358, 1155 (2017).

    ADS  MathSciNet  CAS  PubMed  MATH  Google Scholar 

  3. 3.

    Duvjir, G. et al. Emergence of a metal–insulator transition and high-temperature charge-density waves in VSe2 at the monolayer limit. Nano Lett. 18, 5432 (2018).

    ADS  CAS  PubMed  Google Scholar 

  4. 4.

    Hall, J. et al. Environmental control of charge density wave order in monolayer 2H-TaS2. ACS Nano 13, 10210 (2019).

    CAS  PubMed  Google Scholar 

  5. 5.

    Yu, W. et al. Chemically exfoliated VSe2 monolayers with room-temperature ferromagnetism. Adv. Mater. 31, 1903779 (2019).

    CAS  Google Scholar 

  6. 6.

    Chua, R. et al. Can reconstructed Se-deficient line defects in monolayer VSe2 induce magnetism? Adv. Mater. 32, 2000693 (2020).

    CAS  Google Scholar 

  7. 7.

    Rossnagel, K. On the origin of charge-density waves in select layered transition-metal dichalcogenides. J. Phys. Condens. Matter 23, 213001 (2011).

    ADS  CAS  PubMed  Google Scholar 

  8. 8.

    Li, Y. W. et al. Folded superstructure and degeneracy-enhanced band gap in the weak-coupling charge density wave system 2H-TaSe2. Phys. Rev. B 97, 115118 (2018).

    ADS  CAS  Google Scholar 

  9. 9.

    Zhu, X., Guo, J., Zhang, J. & Plummer, E. W. Misconceptions associated with the origin of charge density waves. Adv. Phys. X 2, 622 (2017).

    CAS  Google Scholar 

  10. 10.

    Johannes, M. D. & Mazin, I. I. Fermi surface nesting and the origin of charge density waves in metals. Phys. Rev. B 77, 165135 (2008).

    ADS  Google Scholar 

  11. 11.

    Inosov, D. S. et al. Fermi surface nesting in several transition metal dichalcogenides. New J. Phys. 10, 125027 (2008).

    ADS  Google Scholar 

  12. 12.

    Berges, J., van Loon, E. G. C. P., Schobert, A., Rösner, M. & Wehling, T. O. Ab initio phonon self-energies and fluctuation diagnostics of phonon anomalies: Lattice instabilities from Dirac pseudospin physics in transition metal dichalcogenides. Phys. Rev. B 101, 155107 (2020).

    ADS  CAS  Google Scholar 

  13. 13.

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

    ADS  CAS  PubMed  Google Scholar 

  14. 14.

    Weber, F. et al. Extended phonon collapse and the origin of the charge-density wave in 2H-NbSe2. Phys. Rev. Lett. 107, 107403 (2011).

    ADS  CAS  PubMed  Google Scholar 

  15. 15.

    Dardel, B. et al. Temperature-dependent pseudogap and electron localization in 1T-TaS2. Phys. Rev. B 45, 1462 (1992).

    ADS  CAS  Google Scholar 

  16. 16.

    Dardel, B. et al. Spectroscopic signatures of phase transitions in a charge-density-wave system: 1T-TaS2. Phys. Rev. B 46, 7407 (1992).

    ADS  CAS  Google Scholar 

  17. 17.

    Zwick, F. et al. Spectral consequences of broken phase coherence in 1T-TaS2. Phys. Rev. Lett. 81, 1058 (1998).

    ADS  CAS  Google Scholar 

  18. 18.

    Lifshitz, I. M. Anomalies of electron characteristics of a metal in the high pressure region. Zh. Eksp. Teor. Fiz. 38, 1569 (1960).

    Google Scholar 

  19. 19.

    Ji, Q. et al. Metallic vanadium disulfide nanosheets as a platform material for multifunctional electrode applications. Nano Lett. 17, 4908 (2017).

    ADS  CAS  PubMed  Google Scholar 

  20. 20.

    Li, Q., Zhou, Q., Shi, L., Chen, Q. & Wang, J. Recent advances in oxidation and degradation mechanisms of ultrathin 2D materials under ambient conditions and their passivation strategies. J. Mater. Chem. A 7, 4291 (2019).

    CAS  Google Scholar 

  21. 21.

    Isaacs, E. B. & Marianetti, C. A. Electronic correlations in monolayer VS2. Phys. Rev. B 94, 035120 (2016).

    ADS  Google Scholar 

  22. 22.

    Tsuda, T., Yasuoka, H., Kitaoka, Y. & Di Salvo, F. J. 51V NMR study of the phase transition in 1T-VS2. J. Magn. Magn. Mater. 31-34, 1101 (1983).

    ADS  CAS  Google Scholar 

  23. 23.

    Mulazzi, M. et al. Absence of nesting in the charge-density-wave system 1T-VS2 as seen by photoelectron spectroscopy. Phys. Rev. B 82, 075130 (2010).

    ADS  Google Scholar 

  24. 24.

    Ma, Y. et al. Evidence of the existence of magnetism in pristine VX2 monolayers (X = S, Se) and their strain-induced tunable magnetic properties. ACS Nano 6, 1695 (2012).

    CAS  PubMed  Google Scholar 

  25. 25.

    Zhuang, H. L. & Hennig, R. G. Stability and magnetism of strongly correlated single-layer VS2. Phys. Rev. B 93, 054429 (2016).

    ADS  Google Scholar 

  26. 26.

    Zhang, H., Liu, L.-M. & Lau, W.-M. Dimension-dependent phase transition and magnetic properties of VS2. J. Mater. Chem. A 1, 10821 (2013).

    CAS  Google Scholar 

  27. 27.

    Murphy, D. W., Cros, C., Di Salvo, F. J. & Waszczak, J. V. Preparation and properties of LixVS2 (0 ≤ x ≤ 1). Inorg. Chem. 16, 3027 (1977).

    CAS  Google Scholar 

  28. 28.

    Magonov, S., Contow, H.-J., Hillebrecht, H., Drechlser, M. & Rotter, W. Atomic structure and superstructure of vanadium disulfide by scanning tunneling microscopy. Optik 83, 60 (1989).

    Google Scholar 

  29. 29.

    Sun, X. et al. In situ unravelling structural modulation across the charge-density-wave transition in vanadium disulfide. Phys. Chem. Chem. Phys. 17, 13333 (2015).

    CAS  PubMed  Google Scholar 

  30. 30.

    Gauzzi, A. et al. Possible phase separation and weak localization in the absence of a charge-density wave in single-phase 1T-VS2. Phys. Rev. B 89, 235125 (2014).

    ADS  Google Scholar 

  31. 31.

    Arnold, F. et al. Novel single-layer vanadium sulphide phases. 2D Mater. 5, 045009 (2018).

    CAS  Google Scholar 

  32. 32.

    Sanders, C. E. et al. Crystalline and electronic structure of single-layer TaS2. Phys. Rev. B 94, 081404(R) (2016).

    ADS  Google Scholar 

  33. 33.

    Lefcochilos-Fogelquist, H. M., Albertini, O. R. & Liu, A. Y. Substrate-induced suppression of charge density wave phase in monolayer 1H-TaS2 on Au(111). Phys. Rev. B 99, 174113 (2019).

    ADS  CAS  Google Scholar 

  34. 34.

    Hall, J. et al. Molecular beam epitaxy of quasi-freestanding transition metal disulphide monolayers on van der Waals substrates: a growth study. 2D Mater. 5, 025005 (2018).

    Google Scholar 

  35. 35.

    Feng, J. et al. Metallic few-layered VS2 ultrathin nanosheets: high two-dimensional conductivity for in-plane supercapacitors. J. Am. Chem. Soc. 133, 17832 (2011).

    CAS  PubMed  Google Scholar 

  36. 36.

    Chen, P. et al. Unique gap structure and symmetry of the charge density wave in single-layer VSe2. Phys. Rev. Lett. 121, 196402 (2018).

    ADS  CAS  PubMed  Google Scholar 

  37. 37.

    Coelho, P. M. et al. Charge density wave state suppresses ferromagnetic ordering in VSe2 monolayers. J. Phys. Chem. C 123, 14089 (2019).

    CAS  Google Scholar 

  38. 38.

    Wong, P. K. J. et al. Evidence of spin frustration in a vanadium diselenide monolayer magnet. Adv. Mater. 31, 1901185 (2019).

    Google Scholar 

  39. 39.

    Miller, D. C., Mahanti, S. D. & Duxbury, P. M. Charge density wave states in tantalum dichalcogenides. Phys. Rev. B 97, 045133 (2018).

    ADS  CAS  Google Scholar 

  40. 40.

    Calandra, M. Phonon-assisted magnetic Mott-insulating state in the charge density wave phase of single-layer 1T-NbSe2. Phys. Rev. Lett. 121, 026401 (2018).

    ADS  CAS  PubMed  Google Scholar 

  41. 41.

    Lian, C.-S., Si, C. & Duan, W. Unveiling charge-density wave, superconductivity, and their competitive nature in two-dimensional NbSe2. Nano Lett. 18, 2924 (2018).

    ADS  CAS  PubMed  Google Scholar 

  42. 42.

    Först, M. et al. Nonlinear phononics as an ultrafast route to lattice control. Nat. Phys. 7, 854 (2011).

    Google Scholar 

  43. 43.

    Truitt, P. A., Hertzberg, J. B., Altunkaya, E. & Schwab, K. C. Linear and nonlinear coupling between transverse modes of a nanomechanical resonator. J. Appl. Phys. 114, 114307 (2013).

    ADS  Google Scholar 

  44. 44.

    Mankowsky, R. et al. Nonlinear lattice dynamics as a basis for enhanced superconductivity in YBa2Cu3O6.5. Nature 516, 71 (2014).

    ADS  CAS  PubMed  Google Scholar 

  45. 45.

    Wang, Y. et al. Evidence of charge density wave with anisotropic gap in a monolayer VTe2 film. Phys. Rev. B 100, 241404 (2019).

    ADS  CAS  Google Scholar 

  46. 46.

    Brik, M. G., Ogasawara, K., Ikeno, H. & Tanaka, I. Fully relativistic calculations of the L2,3-edge XANES spectra for vanadium oxides. Eur. Phys. J. B 51, 345 (2006).

    ADS  CAS  Google Scholar 

  47. 47.

    Scherz, A. et al. Relation between L2,3 XMCD and the magnetic ground-state properties for the early 3d element V. Phys. Rev. B 66, 184401 (2002).

    ADS  Google Scholar 

  48. 48.

    Feng, J. et al. Electronic structure and enhanced charge-density wave order of monolayer VSe2. Nano Lett. 18, 4493 (2018).

    ADS  CAS  PubMed  Google Scholar 

  49. 49.

    Bonilla, M. et al. Strong room-temperature ferromagnetism in VSe2 monolayers on van der Waals substrates. Nat. Nanotechnol. 13, 289 (2018).

    ADS  CAS  PubMed  Google Scholar 

  50. 50.

    Si, J. G. et al. Origin of the multiple charge density wave order in 1T-VSe2. Phys. Rev. B 101, 235405 (2020).

    ADS  CAS  Google Scholar 

  51. 51.

    Cao, Y. et al. Correlated insulator behaviour at half-filling in magic-angle graphene superlattices. Nature 556, 80 (2018).

    ADS  CAS  PubMed  Google Scholar 

  52. 52.

    Tranquada, J. M., Sternlieb, B. J., Axe, J. D., Nakamura, Y. & Uchida, S. Evidence for stripe correlations of spins and holes in copper oxide superconductors. Nature 375, 561 (1995).

    ADS  Google Scholar 

  53. 53.

    van Gastel, R. et al. Selecting a single orientation for millimeter sized graphene sheets. Appl. Phys. Lett. 95, 121901 (2009).

    ADS  Google Scholar 

  54. 54.

    Horcas, I. et al. WSXM: a software for scanning probe microscopy and a tool for nanotechnology. Rev. Sci. Instrum. 78, 013705 (2007).

    ADS  CAS  PubMed  Google Scholar 

  55. 55.

    Giannozzi, P. et al. QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials. J. Phys. Condens. Matter 21, 395502 (2009).

    PubMed  Google Scholar 

  56. 56.

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

    CAS  PubMed  Google Scholar 

  57. 57.

    Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 77, 3865 (1996).

    ADS  CAS  Google Scholar 

  58. 58.

    Perdew, J. P., Burke, K. & Ernzerhof, M. Generalized gradient approximation made simple. Phys. Rev. Lett. 78, 1396(E) (1997).

    ADS  Google Scholar 

  59. 59.

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

    ADS  Google Scholar 

  60. 60.

    van Setten, M. J. et al. The PSEUDODOJO: Training and grading a 85 element optimized norm-conserving pseudopotential table. Comput. Phys. Commun. 226, 39 (2018).

    ADS  Google Scholar 

  61. 61.

    Mostofi, A. A. et al. An updated version of wannier90: a tool for obtaining maximally-localised Wannier functions. Comput. Phys. Commun. 185, 2309 (2014).

    ADS  CAS  MATH  Google Scholar 

Download references


This work was funded by the Deutsche Forschungsgemeinschaft (DFG, German Research foundation) through CRC 1238 (project no. 277146847, subprojects A01 and B06). J.B., A.S., and T.Weh. acknowledge financial support by the DFG through EXC 2077, GRK 2247, and SPP 2244. T.Weh. acknowledges support via the European Graphene Flagship Core3 Project (grant agreement 881603). J.B., A.S., and T.Weh. acknowledge computational resources of the North-German Supercomputing Alliance (HLRN). E.v.L. is supported by the Central Research Development Fund of the University of Bremen. L.M.A. acknowledges financial support from CAPES (project no. 9469/13-3).


Open Access funding enabled and organized by Projekt DEAL.

Author information




T.Weh. and T.M. conceived this work and designed the research strategy. J.H. and discovered the CDW and developed, with the assistance of T.Wek., the growth method. C.v.E. and E.P. conducted and analyzed the STS experiments. J.B., A.S., and G.S., with support from E.v.L. and T.Weh., performed ab initio calculations. F.H. and S.K., with support from N.R., K.O., L.A., N.B., K.K., and H.W., performed and analyzed the XMCD experiment. The results were discussed by all authors. C.v.E., J.B., J.H., E.v.L., S.K., T.Weh., and T.M. wrote the manuscript with input from all authors. The first three authors have comparable contributions to this work.

Corresponding author

Correspondence to Camiel van Efferen.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature Communications thanks Matteo Calandra and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.

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

van Efferen, C., Berges, J., Hall, J. et al. A full gap above the Fermi level: the charge density wave of monolayer VS2. Nat Commun 12, 6837 (2021).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI:


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

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