
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

51608
10.1038/s41467-024-51608-y
Article
Proximity induced charge density wave in a graphene/1T-TaS2 heterostructure
http://orcid.org/0000-0002-3804-555X
Tilak Nikhil 1
http://orcid.org/0000-0002-2122-5888
Altvater Michael 1
Hung Sheng-Hsiung 2
Won Choong-Jae 3
Li Guohong 1
Kaleem Taha 1
Cheong Sang-Wook 1
http://orcid.org/0000-0002-1502-6022
Chung Chung-Hou chung0523@nycu.edu.tw

456
http://orcid.org/0000-0002-2881-3826
Jeng Horng-Tay jeng@phys.nthu.edu.tw

257
http://orcid.org/0000-0002-2516-2749
Andrei Eva Y. eandrei@physics.rutgers.edu

1
1 https://ror.org/05vt9qd57 grid.430387.b 0000 0004 1936 8796 Department of Physics and Astronomy, Rutgers, the State University of New Jersey, Piscataway, New Jersey USA
2 https://ror.org/00zdnkx70 grid.38348.34 0000 0004 0532 0580 Department of Physics, National Tsing Hua University, Hsinchu, Taiwan
3 https://ror.org/04xysgw12 grid.49100.3c 0000 0001 0742 4007 Laboratory for Pohang Emergent Materials and Max Plank POSTECH Center for Complex Phase Materials, Department of Physics, Pohang University of Science and Technology, Pohang, Korea
4 https://ror.org/00se2k293 grid.260539.b 0000 0001 2059 7017 Department of Electrophysics, National Yang Ming Chiao Tung University, Hsinchu, Taiwan
5 https://ror.org/02mfp0b67 grid.468468.0 0000 0000 9060 5564 Physics Division, National Center for Theoretical Sciences, Taipei, Taiwan
6 https://ror.org/00se2k293 grid.260539.b 0000 0001 2059 7017 Center for Theoretical and Computational Physics, National Yang Ming Chiao Tung University, Hsinchu, Taiwan
7 grid.28665.3f 0000 0001 2287 1366 Institute of Physics, Academia Sinica, Taipei, Taiwan
14 9 2024
14 9 2024
2024
15 805629 5 2024
7 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
The proximity-effect, whereby materials in contact appropriate each other’s electronic-properties, is widely used to induce correlated states, such as superconductivity or magnetism, at heterostructure interfaces. Thus far however, demonstrating the existence of proximity-induced charge-density-waves (PI-CDW) proved challenging. This is due to competing effects, such as screening or co-tunneling into the parent material, that obscured its presence. Here we report the observation of a PI-CDW in a graphene layer contacted by a 1T-TaS2 substrate. Using scanning tunneling microscopy (STM) and spectroscopy (STS) together with theoretical-modeling, we show that the coexistence of a CDW with a Mott–gap in 1T-TaS2 coupled with the Dirac-dispersion of electrons in graphene, makes it possible to unambiguously demonstrate the PI-CDW by ruling out alternative interpretations. Furthermore, we find that the PI-CDW is accompanied by a reduction of the Mott gap in 1T-TaS2 and show that the mechanism underlying the PI-CDW is well-described by short-range exchange-interactions that are distinctly different from previously observed proximity effects.

Electronic proximity effects between neighbouring nanoscale materials are an important theme in condensed matter physics. Here, the authors experimentally demonstrate a proximity induced charge density wave in graphene due to the nearby presence of 1T-TaS2.

Subject terms

Nanoscience and technology
Graphene
https://doi.org/10.13039/100000936 Gordon and Betty Moore Foundation (Gordon E. and Betty I. Moore Foundation) FG02-99ER45742 GBMF9453 Tilak Nikhil Andrei Eva Y. https://doi.org/10.13039/100006151 DOE | SC | Basic Energy Sciences (BES) GBMF9453 FG02-99ER45742 Tilak Nikhil Andrei Eva Y. issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

The isolation and manipulation of atomically thin materials provide a ready-made two-dimensional electron system1 whose properties can be tuned by external knobs such as stress, substrate morphology, or doping2–8, to facilitate the emergence of correlated electron phases. Distinct from these external knobs, a powerful approach to manipulate electron correlations is by contact proximity effects. It is well known that proximitizing materials that host correlated electron phases with a normal metal, induces correlations in the metal9. This is a consequence of the quantum mechanical properties of electrons in solids; specifically, the non-local nature of electrons. As quantum particles do not have a well-defined position, electronic states cannot abruptly change from one type of ordering to another at the interface between two materials. Thus, correlated states persist into the normal metal where scattering events begin to destroy the coherence (and vice versa). In the case of 2D materials where scattering is reduced due to their atomically sharp interfaces, proximity-effects are particularly robust allowing proximitized states to persist over long distances. The discovery of graphene and other 2D materials, together with technology enabling the fabrication of 2D heterostructures, has led to the observation of strong proximity effects at the atomic limit, including proximity-induced superconductivity, magnetism10, and spin-orbit effects11–18.

CDW states19,20 that are typically induced by Fermi surface instabilities driven by electron-phonon coupling21 or by correlation-enhanced exchange effects22,23, are quite robust. By contrast, PI-CDW in contacted 3D metals are fragile and easily destroyed by interface defect scattering. The use of 2D heterostructures fabricated in an inert atmosphere to avoid interface damage and contamination has made it possible to address this difficulty by substantially reducing interface scattering24,25. However, unambiguously identifying a PI-CDW in a 2D heterostructure is challenging due to competing mechanisms such as screening, doping, or strain, which can also produce charge modulations. In addition, as is the case of local probes such as STS, co-tunneling into the parent material can obscure the proximity induced CDW. Therefore, in order to demonstrate a PI-CDW one must: (i) demonstrate the presence of a CDW in the contacted material with the same period as that of the parent material (ii) rule out misidentification of electrostatic screening effects as a PI-CDW (iii) rule out misidentification of the parent CDW as the PI-CDW. Thus far studies of 2D heterostructures comprised of graphene and transition metal dichalcogenides (TMD) were unable to satisfy the three necessary criteria to unambiguously demonstrate the existence of a PI-CDW24–27.

In this work, we present evidence of PI-CDW in a graphene layer contacted by a 1T-TaS2 substrate hosting a commensurate CDW (CCDW). This discovery was made possible by the unique coexistence of the CDW and Mott gap in 1T-TaS2 together with the Dirac spectrum of graphene. This exceptional combination, which was not present in any of the previously studied systems24–27, enabled us to unambiguously establish the presence of the PI-CDW by ruling out charge density modulations induced by trivial co-tunneling processes or by screening effects. By comparing STM and STS measurements with first-principle calculations as well as with mean-field theory results, we demonstrate that the charge density modulation of the CCDW in 1T-TaS2 persists within the contacted graphene layer. Further, we identify the effect of charge transfer on the band structure of 1T-TaS2 and reveal that the PI-CDW is described by a mechanism of short-range exchange interactions.

Results and discussion

1T-TaS2 is comprised of hexagonal layers of tantalum atoms coordinated octahedrally by sulfur atoms. Below 550 K a high-temperature metallic phase transits into an incommensurate CDW phase, followed by a nearly commensurate CDW (NCCDW) below 350 K, and a CCDW below 180 K24 (Fig. 1c). The CDW unit cell consists of a 13-atom cluster where 12 atoms displace from their high-temperature equilibrium positions toward the central, 13th Ta atom forming a √13 × √13 reconstructed supercell (Fig. 1b) also known as a David star (DS) structure. The CDW in 1T-TaS2 is driven by strong electron-phonon coupling, which critically suppresses an acoustic phonon mode along the Γ-M direction28,29 (Supplementary Fig. 4) leading to the static displacement of the lattice with the wavevector QCDW, corresponding to the √13 × √13 CDW. This soft phonon mode consists primarily of longitudinal vibrations of Ta atoms with a minor contribution from transverse vibrations of S atoms relative to the phonon propagating direction QCDW. The DS atomic arrangement involves the local lattice contraction around the center of the star, in which the bond lengths between Ta ions are shorter than those between Ta ions outside the DS.Fig. 1 Scanning tunneling microscopy of graphene covered 1T-TaS2.

a Schematic of the device and STM measurement setup. b Cartoon showing a single star-of-David (red) in 1T-TaS2. The black arrows denote the displacement of Ta atoms towards a central Ta atom. c The transition temperatures between the incommensurate CDW (ICCDW), nearly commensurate CDW (NCCDW), and commensurate CDW (CCDW) phases of 1T-TaS2. d Cartoon showing a monolayer graphene layer on top of 1T-TaS2. Here, the angle between the graphene layer and 1T-TaS2 is 13.9° to form a commensurate unit cell (cyan) for DFT calculations. e An STM topography image of graphene-covered 1T-TaS2 measured at 77 K using a tunneling set point of I = 150 pA, Vb = 250 mV. Each bright spot corresponds to a star of David with a measured CDW wavelength of 1.2 nm. This can also be seen in the 2D-FFT of the topography data (inset). f An STM topography image and its 2D FFT (inset) of a bulk 1T-TaS2 sample measured at 77 K using a tunneling set point of I = 150 pA, Vb = 500 mV. The crystals in panels b and d were plotted using VESTA49.

We studied heterostructures of thin 1T-TaS2 flakes (7–50 nm) covered by monolayer graphene (graphene/1T-TaS2) using STM and STS (Fig. 1a). The samples were assembled inside a dry argon-filled glovebox (< 0.1 ppm O2, H2O)24,30 to avoid oxidation of the air-sensitive surface of 1T-TaS2 and subsequently cooled to 77 K inside a home-built STM31,32. At this temperature, 1T-TaS2 is in the CCDW regime. The STM topography (Fig. 1e) measured on a graphene/1T-TaS2 sample using a positive bias set-point (Vb = 250 mV) shows a triangular array of tall spots. The lattice spacing for this array is about 1.2 nm as measured directly from the topography image and confirmed with its fast Fourier transform (FFT) (Fig. 1e, inset). This spacing is equal to the expected CDW wavelength (LCDW) in 1T-TaS2 indicating that each tall spot corresponds to one DS. The topography of the graphene-covered sample thus reflects the distinct arrangement of the DS clusters which are the hallmark of the CDW in 1T-TaS233. The topography is identical to that measured on a bulk 1T-TaS2 sample without a graphene cover under similar tunneling conditions (Fig. 1f). This result also holds true for negative bias set points of Vb < −300 mV. In other words, for sufficiently large bias voltages, of either sign, there is no noticeable difference in the STM topography of 1T-TaS2 with or without a graphene cover. This similarity might suggest that there is little or no interaction between the two materials. However, spectroscopy and low bias-voltage topography scans show that this is not the case.

The high-temperature band structure of bulk 1T-TaS2 consists of an isolated, half-filled spin-degenerate flat band at the Fermi level (EF) (Supplementary Fig. 5). In the CCDW phase, the flat band splits into an occupied spin up lower Hubbard band (LHB) and empty spin down upper Hubbard band (UHB) with a Mott gap of ~ 0.5 eV in between (ab initio calculation Fig. 2c). The other bands corresponding to the CDW distortion appear at energies below the LHB (CDW-VB) and above the UHB (CDW-CB). These spectroscopic features are clearly observed in the experimental dI/dV spectrum (Fig. 2a) measured on a bulk 1T-TaS2 sample, in agreement with previous reports20,23,34,35.Fig. 2 Effects of interaction between graphene and 1T-TaS2.

a, b are the measured dI/dV curves on bulk 1T-TaS2 and graphene/1T-TaS2 respectively. The positions of the Hubbard band peaks and the broad features originating from the CDW distortion are labeled. The red, magenta, and blue dashed lines in (b) indicate the bias voltages at which topography scans (g–i) were measured, respectively. c The calculated DOS for bulk 1T-TaS2. d The calculated DOS for the heterostructure decomposed into the partial DOS of graphene (blue curve) and 1T-TaS2 (red curve). The graphene DOS is multiplied by a factor (39) before adding the 1T-TaS2 DOS (black curve) for comparison with the experiment. This accounts for the smaller tip-graphene separation. See SI for how the factor was estimated. The linear dispersion around the Dirac point of graphene at ~ 0.3 eV above EF indicates charge transfer from graphene to 1T-TaS2. e A map of dI/dV (Vb) vs position across 8 DSs in graphene/1T-TaS2. The intensity of both Hubbard bands is highest at the center of the DSs. f Three representative spectra from the map in (e) demonstrate this intensity variation. g–i STM topography scans measured at bias voltages of − 300 mV, − 150 mV, and + 300 mV respectively indicated by dashed lines in (b). The tunneling resistance was kept constant at 1 GΩ for all scans. The small lateral mismatch between scans is a result of the thermal drift of the sample. All scale bars are 2 nm.

The spectral features of graphene/1T-TaS2 (Fig. 2b), like those of bare 1T-TaS2, show two peaks at ~ − 270 meV and + 70 meV, identified as the LHB and UHB respectively. This identification is supported by the intensity variation of these features as a function of distance from the DS centers (Fig. 2e, f). The main differences between the 1T-TaS2 and graphene/1T-TaS2 spectra are the shift of the Fermi level and the linearly dispersing states within the Mott gap. The former can be understood in terms of the interfacial charge transfer caused by the work function difference between graphene (ΦG ~ 4.6 eV36) and 1T-TaS2 (Φ1T-TaS2 ~ 5.2 eV37). Consequently, the Fermi level moves from the LHB top to the UHB bottom, while the graphene Dirac point shifts to a higher energy of about 0.3 eV above EF, indicating hole doping. The linearly dispersing Dirac cone states give rise to a linear band within the Mott gap of the 1T-TaS2. While the DFT model reproduces the presence of graphene-attributed density of states at the Fermi level, we find that the experimental dI/dV spectrum at higher energies is qualitatively different from the calculated DOS. The discrepancy arises from the exponential sensitivity of the signal on the distance between the STM tip and electronic states within the sample structure. Taking into account the decay of the wavefunction from the graphene surface to the 1T-TaS2 substrate, we find that the contribution of the graphene layer to the calculated DOS should be multiplied by ~ 39 compared to that of the 1T-TaS2 layer (Supporting Information). The weighted sum (black) of the projected DOS of 1T-TaS2 and that of graphene multiplied by 39 (blue) reproduces the asymmetric V-shaped DOS measured experimentally (Fig. 2d). While STM measurements of van der Waals heterostructures are common, de-convolving the different layer contributions to spectroscopy is not. Our procedure provides a computational method to accurately simulate STS measurements of 2D material heterostructures based on DFT models.

Another notable difference between the 1T-TaS2 and the graphene/1T-TaS2 spectra is that the separation between the Hubbard peaks in graphene/1T-TaS2 (~ 338 mV) is smaller by about 31% compared to the bare 1T-TaS2 (~ 491 mV). We attribute this to the reduction of the on-site Coulomb repulsion (U) in the graphene/1T-TaS2 sample due to screening by the highly mobile carriers in graphene. These carriers suppress the localized picture of the narrow Hubbard d-bands in 1T-TaS2 and thus reduce the on-site U value of Ta. Comparing the ab initio calculated band structures (supplementary Fig. 6) for different U values with the STS results suggests that the Hubbard U is lowered from 2.27 eV to 1.70 eV. Using U = 1.70 eV in our ab initio calculations (Fig. 2d) we find that the calculated DOS reproduces the main features of the observed STS. Our observation reveals important microscopic details about the effect of screening on correlated insulators that were not initially captured by the DFT model. Similar discrepancies between the measured and calculated DOS in heterostructures of graphene and sister TMDs (1T-NbSe2 and 1T-TaSe238–43) may be resolved by taking into account the screening effects of the conducting substrates.

In addition to the local height variation in the sample, STM topography also depends on the total number of states in the energy window [EF,EF+eVb]. For CDW materials, the local DOS variation dominates the topography signal. Since the integrated DOS of 1T-TaS2 is orders of magnitude larger than that of graphene for large Vb (Fig. 3b), we estimate that there is substantial tunneling current into both the graphene and 1T-TaS2, despite the larger z separation between the tip and the 1T-TaS2 layer (Fig. 3a). Consequently, at high bias voltages the topography scans on graphene/1T-TaS2 resemble those on 1T-TaS2, as shown for Vb = −300 mV (Fig. 2g) and Vb = + 300 mV (Fig. 2i). Note that these two energy windows include the LHB and UHB of the 1T-TaS2, respectively.Fig. 3 Modeling co-tunneling and electrostatic screening.

a A schematic cross-section of the tunneling junction. The mean graphene-1T-TaS2 separation is ~ 0.55 nm. Co-tunneling into the graphene and 1T-TaS2 layers was estimated by measuring the tunneling current decay constant z0. b The integrated DOS of the graphene (blue curve) and the 1T-TaS2 (red curve) in the energy window [EF, EF + eVb] were calculated using the DOS from DFT calculations. The integrated DOS for 1T-TaS2 is about an order of magnitude larger than that of graphene for large Vb. For Vb inside the Mott gap, the graphene DOS dominates. c A schematic that considers the effect of screening of the CDW in 1T-TaS2 by the electrons in graphene. The bottom subplot shows the electron density (ρ) in the 1T-TaS2 as a function of position. The resulting electric potential (U(r)) denoted by the orange line, tunes the local charge density of the graphene layer. This is denoted by the shift of the charge neutrality point of graphene’s Dirac cones (blue) at two representative positions. The Fermi level is denoted by the dashed black line. The red squiggly lines denote the Hubbard bands of the 1T-TaS2. The cyan region shows the states within the energy window [EF, EF + eVb] for the two positions. There are more states available for tunneling in the graphene above the center of the DSs compared to the sides. This should result in a triangular CDW like pattern in topography measured at these tunneling conditions. This is at odds with the observed hexagonal honeycomb pattern. d A dI/dV spectrum measured on a bilayer graphene/1T-TaS2 area which shows that the CNP is clearly around 0.3 eV. The Hubbard bands are greatly suppressed compared to the monolayer graphene/1T-TaS2 owing to the increased separation of the topmost graphene layer from the 1T-TaS2, additional screening, and the lack of direct hybridization between the two materials.

Surprisingly, the topography pattern measured at bias voltages within the Mott gap of 1T-TaS2 is hexagonal instead of the triangular pattern observed outside the Mott gap, but with the same period (Fig. 2h). We note that for bias voltages within the Mott gap of 1T-TaS2, the STM topography includes only states belonging to graphene (Fig. 3b) because tunneling into the 1T-TaS2 layer is disallowed. Therefore, since a hexagonal lattice is the dual of a triangular lattice, we conclude that the electronic density redistribution in graphene due to its proximity to 1T-TaS2 forms a CDW pattern which is out-of-phase with the CDW in the underlying 1T-TaS2. This result excludes an interpretation of the data in terms of co-tunneling into the 1T-TaS2 substrate. The hexagonal pattern was observed for multiple energy windows inside the 1T-TaS2 Mott gap but never outside it (Supplementary Fig. 10)

Note that this out-of-phase CDW cannot be explained by a simple screening picture (Fig. 3c). To wit: if one models the CDW in 1T-TaS2 as a periodic modulation of electron density (ρ(r)) which peaks at the center of the DSs (Fig. 3c, bottom), then the electrostatic potential (U(r)) near the surface of 1T-TaS2 will also peak near the center of the DSs (orange line in Fig. 3c). A graphene layer placed on top of 1T-TaS2 will experience this spatially varying potential leading to spatial doping variation, i.e., shifting of the local CNP. This CNP variation is in addition to the average shift of the CNP (~ 0.3 eV) due to the work function difference between the two materials as discussed above. If we consider an energy window in the Mott gap of the 1T-TaS2, then there are more graphene states (cyan regions in Fig. 3b) near the center of the DS compared to the region between neighboring DSs. This should result in a triangular CDW pattern in STM topography, which is in phase with the underlying 1T-TaS2 CDW. This is contrary to our observation of a hexagonal pattern, thus ruling out this simple picture. These results provide strong evidence for the existence of a PI-CDW in graphene.

We also studied an area of the sample covered with bilayer graphene. In this area, the topography showed almost no CDW pattern. dI/dV spectra (Fig. 3d) showed a bilayer graphene-like V shape with a pronounced dip near 0.3 V, indicating the CNP. Very weak features corresponding to the position of the Hubbard bands were also seen. These observations indicate that the PI-CDW relies on very short-range coupling between the graphene and 1T-TaS2.

We performed first-principle electronic structure calculations using the projector augmented wave (PAW) approach within the framework of density functional theory (DFT) as implemented in the Vienna ab initio Simulation Package (VASP)44,45. The exchange-correlation is described in the Perdew-Burke-Ernzerhof (PBE) form of generalized gradient approximation (GGA)46,47. To take the strong correlation of Ta d-electrons into consideration, we perform GGA plus on-site U (GGA + U) calculations with U = 2.27 eV for monolayer 1T-TaS2 in accordance with previous DFT calculations33,48. The lattice structure is theoretically optimized with the atomic forces converging within 0.01 eV/Ang.

Next, we calculated the GGA + U band structure of a monolayer 1T-TaS2 covered with graphene using the reduced value U = 1.7 eV (Fig. 4d) to reproduce the observed spacing between the Hubbard bands. To enable this calculation, we considered a twist angle of 13.9 degrees between the graphene and 1T-TaS2 layers. This forms a commensurate superstructure (Fig. 4a) because 5×aG≈LCDW, where aG = 0.246 nm is the lattice constant of graphene (Supplementary Fig. 8). We note that both the 1T-TaS2 Hubbard bands and the graphene Dirac cone are preserved (Fig. 4d) in the band structure. At the intersections between the Dirac cone with the Hubbard bands, four small gaps of ~ 10 meV emerge (Fig. 4d insets). These gaps, which originate from the weak interlayer couplings, indicate finite hybridization between the Dirac and Hubbard states, giving rise to the proximity effect between graphene and 1T-TaS2 as illustrated in Fig. 4b. In addition, we found a gap associated with the CDW distortion (SI Fig. 9), which was too small (~ 1.7 meV) to detect experimentally.Fig. 4 Ab initio band structure calculations.

a Top view of graphene/1T-TaS2. Blue, brown, and yellow spheres indicate C, Ta, and S atoms, respectively. Black, blue, and brown rhombuses show the 5 × 5 graphene/√13 × √13 1T-TaS2 supercell, graphene 1 × 1 unit cell, and 1T-TaS2 1 × 1 unit cell, respectively. The calculations were carried out for a ~ 13.9° twist angle between graphene and 1T-TaS2. b Side view of graphene/1T-TaS2 overlaid with the charge density map (bubble-shaped color contour indicating the electron density) corresponding to the states at the two crossing points of the Dirac cone and Lower Hubbard band. The slight overlap between graphene and 1T-TaS2 electron clouds gives rise to the interlayer coupling and proximity effect. Here, blue, red, and green spheres represent C, Ta, and S atoms, respectively. c GGA + U spin resolved band structure of √13 × √13 CDW reconstructed 1T-TaS2 with U = 2.27 eV. d GGA + U spin resolved band structure of graphene/1T-TaS2 using the phenomenological value U = 1.70 eV. The red and blue curves are the spin-up and down bands, respectively. Owing to the charge transfer from graphene to 1T-TaS2, the Fermi level (zero energy) moves from the lower Hubbard band to the upper Hubbard band, and the graphene-associated Dirac point is shifted to ~ 0.3 eV above the Fermi level indicating hole doping. Insets on the right are zoomed in views of the bands near the upper and lower Hubbard bands. e–g The integrated DOS of graphene/1T-TaS2 in the energy range from EF to − 300 mV, − 200 mV, and + 300 mV, respectively, relative to that of freestanding graphene in corresponding energy ranges. DSs are overlaid for clarity. For the energy above the UHB (e) and below the LHB (g), the integrated DOS forms a triangular lattice. While within the Mott gap (f), the in-gap DOS forms a hexagonal lattice, in good agreement with our experimental observations in Fig. 2(g–i).

Since different carbon atoms in the heterostructure unit-cell have a different registry with respect to the underlying CDW, it is natural to expect that the hybridization between graphene pz-bands and 1T-TaS2 d-bands varies as a function of position within the supercell. This can be seen in the calculated charge density map (Fig. 4b). To show the charge modulation of graphene in the proximity CDW state (Fig. 4e–g) for comparison with the experimental topography maps (Fig. 2h–g), we calculate the integrated DOS in the energy windows [EF, EF + eVb ] for Vb = − 0.3 V, − 0.2 V and + 0.3 V by subtracting the counterparts of freestanding graphene in the corresponding energy ranges, i.e., in the energy windows [− 0.3 eV, − 0.6 eV], [− 0.3 eV, − 0.5 eV] and [− 0.3 eV, 0.0 eV] with 0.0 eV being the Dirac point (EF) of freestanding graphene. Consistent with the experimental results, we find triangular, hexagonal, and triangular CDW patterns below, within, and above the Mott gap region, respectively, providing further support for the out-of-phase PI-CDW in the graphene layer.

The PI-CDW is an effect driven by periodically varying band hybridization. In addition to the DFT calculations, we further modeled the graphene/1T-TaS2 system using a simplified mean-field Hamiltonian (detailed in the SI) based on the existing CDW order in 1T-TaS2 with a weak charge transfer between the graphene and 1T-TaS2 layers. The charge transfer is described by a short-ranged exchange term produced by second-order perturbation of the weak interlayer hopping. Using a mean-field decoupling of this exchange term, a PI-CDW order (much weaker than that in 1T-TaS2), is generated in the graphene layer. We find that the CDW in graphene is out of phase with that in the 1T-TaS2 layer, consistent with both the DFT and the experimental results. In addition, the mean-field theory predicts that the amplitude of the CDW in graphene is much weaker than that in 1T-TaS2, again consistent with the DFT and the STM results. These findings provide strong support for the CDW proximity effect in graphene. We emphasize that the above mechanism based on charge transfer is distinct from previously realized proximity effects, including superconducting, magnetic, and spin-orbit proximity effects.

In summary, using STM/STS together with DFT calculations and mean field Hamiltonian to study graphene/1T-TaS2, we demonstrated the existence of a PI-CDW in graphene and elucidated details of the coupling between the two systems. This was made possible by the unique properties of these two materials, the coexistence of the CDW and Mott gap in 1T-TaS2 together with the Dirac spectrum of electrons in graphene, which enabled us to distinguish the PI-CDW in graphene from the contribution of the host CDW, as well as to rule out a screening-induced charge modulation. We propose a model based on the short-range exchange interaction between carriers in graphene and in 1T-TaS2 which captures the main features of the CDW proximity effect. Concomitant with the PI-CDW in graphene, we observe a reduction in the Mott gap in 1T-TaS2, indicating the presence of proximity-induced mid-gap carriers that screen the Mott–Hubbard interaction. The graphene/1T-TaS2 system provides a sensitive and discriminating probe of contact-induced contributions such as charge transfer, screening, and exchange interactions at the atomic level. Future studies using techniques such as STM/STS and spin-resolved STM to focus on these effects in heterostructures comprised of TMDs and metallic layers could provide new insights into correlated insulator physics and quantum interactions at vdW interfaces, crucial for identifying emergent phenomena, and for understanding quantum transport through 2D materials embedded in 3D structures.

Methods

Sample fabrication and STM measurements

Samples were fabricated by mechanical exfoliation of graphene and separately 1T-TaS2 flakes inside an argon-filled glovebox. 1T-TaS2 flakes were exfoliated from bulk 1T-TaS2 crystals (samples were purchased from 2D Semiconductors or grown by iodine chemical vapor transport) and transferred onto a passivated SiO2-capped degenerately doped Si wafer. The graphene and 1T-TaS2 flakes were aligned vertically, brought into close contact with micromanipulators under an optical microscope, and then heated to promote adhesion. Standard electron beam lithography and electrode deposition (4–5 nm Ti/40–50 nm Au) were used to make electrical contact with the sample. After removing the PMMA mask, the resulting heterostructure was annealed (180 − 220 °C) in hydrogen/argon (10%: 90%) to remove polymer residues followed by AFM tip sweeping. STM and STS were performed using a homebuilt STM at 77 K in high-vacuum < 10−5 Torr. To locate the micron-size samples we employed a technique using the STM tip (mechanically cut Pt/Ir) as a capacitive antenna32. Two such samples were measured by STM using several mechanically cut Pt/Ir wire tips. In addition, we also studied a 1T-TaS2 flake transferred to a pre-patterned gold electrode using PDMS polymer. This sample was transferred to the STM system with < 15 mins of exposure to the atmosphere.

Supplementary information

Supplementary Information

Peer Review File

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-51608-y.

Acknowledgements

N.T. and E.Y.A. acknowledge support from the Department of Energy grant DOE-FG02-99ER45742 and the Gordon and Betty Moore Foundation EPiQS initiative grant GBMF9453; MAA was supported by the National Science Foundation grant EFRI 1433307; G.L. was supported by Rutgers University; C.J.W. was supported by The National Research Foundation of Korea (NRF), Ministry of Science and ICT(No. 2022M3H4A1A04074153); S.W.C. was supported by the Center for Quantum Materials Synthesis (cQMS), funded by the Gordon and Betty Moore Foundation’s EPiQS initiative grant GBMF10104, and by Rutgers University, C.-H.C. was supported by MOST (Grant NO.: 107-2112-M-009-010-MY3, 110-2112-M-A49-018-MY3), NSTC (Grant NO. 112-2124-M-A49-003), and NCTS of Taiwan, R.O.C., JHT acknowledges support from the Ministry of Science and Technology, Taiwan under grant: MOST 109-2112-M-007 −034 -MY3, and from NCHC, CINC-NTU, AS-iMATE-109-13, and CQT-NTHU-MOE, Taiwan.

Author contributions

N.T. and M.A. built the STM, performed the experimental studies, and analysis, and wrote the paper; S.-H.H. and H.-T.J. performed the first principles calculations; C.-J.W. grew one of the TaS2 samples; G.L. built the STM and helped with data collection; T.K. helped with sample fabrication; S.-W.C. helped obtain the TaS2 sample; C.-H.C. carried out the mean-field calculations; E.Y.A. directed the project, analyzed data, and wrote the paper.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewer(s) for their contribution to the peer review of this work. A peer review file is available.

Data availability

The experimental data generated in this study have been deposited in the Zenodo database under the accession code 10.5281/zenodo.12745631.

Competing interests

The authors declare no competing interests.

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

These authors contributed equally: Nikhil Tilak, Michael Altvater, Sheng-Hsiung Hung.
==== Refs
References

1. Geim AK Grigorieva IV Van der Waals heterostructures Nature 2013 499 419 425 10.1038/nature12385 23887427
Geim, A. K. & Grigorieva, I. V. Van der Waals heterostructures. Nature 499, 419–425 (2013).23887427 10.1038/nature12385
2. Li G Observation of Van Hove singularities in twisted graphene layers Nat. Phys. 2010 6 109 109 10.1038/nphys1463
Li, G. et al. Observation of Van Hove singularities in twisted graphene layers. Nat. Phys. 6, 109–109 (2010).10.1038/nphys1463
3. Levy N Strain-induced pseudo-magnetic fields greater than 300 tesla in graphene nanobubbles Science 2010 329 544 547 10.1126/science.1191700 20671183
Levy, N. et al. Strain-induced pseudo-magnetic fields greater than 300 tesla in graphene nanobubbles. Science 329, 544–547 (2010).20671183 10.1126/science.1191700
4. Carrillo-Bastos R Strained fold-assisted transport in graphene systems Phys. Rev. B 2016 94 125422 125422 10.1103/PhysRevB.94.125422
Carrillo-Bastos, R. et al. Strained fold-assisted transport in graphene systems. Phys. Rev. B 94, 125422–125422 (2016).10.1103/PhysRevB.94.125422
5. Jiang Y Visualizing strain-induced pseudomagnetic fields in graphene through an hBN magnifying glass Nano Lett. 2017 17 2839 2843 10.1021/acs.nanolett.6b05228 28409936
Jiang, Y. et al. Visualizing strain-induced pseudomagnetic fields in graphene through an hBN magnifying glass. Nano Lett. 17, 2839–2843 (2017).28409936 10.1021/acs.nanolett.6b05228
6. Cao Y Correlated insulator behaviour at half-filling in magic-angle graphene superlattices Nature 2018 556 80 84 10.1038/nature26154 29512654
Cao, Y. et al. Correlated insulator behaviour at half-filling in magic-angle graphene superlattices. Nature 556, 80–84 (2018).29512654 10.1038/nature26154
7. Mao J Evidence of flat bands and correlated states in buckled graphene superlattices Nature 2020 584 215 220 10.1038/s41586-020-2567-3 32788735
Mao, J. et al. Evidence of flat bands and correlated states in buckled graphene superlattices. Nature 584, 215–220 (2020).32788735 10.1038/s41586-020-2567-3
8. Andrei EY MacDonald AH Graphene bilayers with a twist Nat. Mater. 2020 19 1265 1275 10.1038/s41563-020-00840-0 33208935
Andrei, E. Y. & MacDonald, A. H. Graphene bilayers with a twist. Nat. Mater. 19, 1265–1275 (2020).33208935 10.1038/s41563-020-00840-0
9. Meissner H Superconductivity of contacts with interposed barriers Phys. Rev. 1960 117 672 680 10.1103/PhysRev.117.672
Meissner, H. Superconductivity of contacts with interposed barriers. Phys. Rev. 117, 672–680 (1960).10.1103/PhysRev.117.672
10. Hauser JJ Magnetic proximity effect Phys. Rev. 1969 187 580 583 10.1103/PhysRev.187.580
Hauser, J. J. Magnetic proximity effect. Phys. Rev. 187, 580–583 (1969).10.1103/PhysRev.187.580
11. Zhao R Tuning phase transitions in 1T-TaS2 via the substrate Nano Lett. 2017 17 3471 3477 10.1021/acs.nanolett.7b00418 28463560
Zhao, R. et al. Tuning phase transitions in 1T-TaS2 via the substrate. Nano Lett. 17, 3471–3477 (2017).28463560 10.1021/acs.nanolett.7b00418
12. Frano A Long-range charge-density-wave proximity effect at cuprate/manganate interfaces Nat. Mater. 2016 15 831 834 10.1038/nmat4682 27322824
Frano, A. et al. Long-range charge-density-wave proximity effect at cuprate/manganate interfaces. Nat. Mater. 15, 831–834 (2016).27322824 10.1038/nmat4682
13. Li AJ Zhu X Stewart GR Hebard AF Bi-2212/1T-TaS2 Van der Waals junctions: Interplay of proximity induced high-Tcsuperconductivity and CDW order Sci. Rep. 2017 7 4639 10.1038/s41598-017-04645-1 28680063
Li, A. J., Zhu, X., Stewart, G. R. & Hebard, A. F. Bi-2212/1T-TaS2 Van der Waals junctions: Interplay of proximity induced high-Tcsuperconductivity and CDW order. Sci. Rep. 7, 4639 (2017).28680063 10.1038/s41598-017-04645-1
14. Dreher P Proximity effects on the charge density wave order and superconductivity in single-layer NbSe2 ACS Nano 2021 15 19430 19438 10.1021/acsnano.1c06012 34846856
Dreher, P. et al. Proximity effects on the charge density wave order and superconductivity in single-layer NbSe2. ACS Nano 15, 19430–19438 (2021).34846856 10.1021/acsnano.1c06012
15. Du, X., Skachko, I. & Andrei, E. Y. Josephson current and multiple Andreev reflections in graphene SNS junctions. Phys. Rev. B 77, 184507 (2008).
16. Heersche HB Jarillo-Herrero P Oostinga JB Vandersypen LM Morpurgo AF Bipolar supercurrent in graphene Nature 2007 446 56 59 10.1038/nature05555 17330038
Heersche, H. B., Jarillo-Herrero, P., Oostinga, J. B., Vandersypen, L. M. & Morpurgo, A. F. Bipolar supercurrent in graphene. Nature 446, 56–59 (2007).17330038 10.1038/nature05555
17. Avsar A Spin-orbit proximity effect in graphene Nat. Commun. 2014 5 4875 4875 10.1038/ncomms5875 25255743
Avsar, A. et al. Spin-orbit proximity effect in graphene. Nat. Commun. 5, 4875–4875 (2014).25255743 10.1038/ncomms5875
18. Liang X The magnetic proximity effect and electrical field tunable valley degeneracy in MoS2/EuS van der Waals heterojunctions Nanoscale 2017 9 9502 9509 10.1039/C7NR03317F 28660948
Liang, X. et al. The magnetic proximity effect and electrical field tunable valley degeneracy in MoS2/EuS van der Waals heterojunctions. Nanoscale 9, 9502–9509 (2017).28660948 10.1039/C7NR03317F
19. Keimer B Kivelson SA Norman MR Uchida S Zaanen J From quantum matter to high-temperature superconductivity in copper oxides Nature 2015 518 179 186 10.1038/nature14165 25673411
Keimer, B., Kivelson, S. A., Norman, M. R., Uchida, S. & Zaanen, J. From quantum matter to high-temperature superconductivity in copper oxides. Nature 518, 179–186 (2015).25673411 10.1038/nature14165
20. Lin H Scanning tunneling spectroscopic study of monolayerT-TaS 2 and 1 T-TaSe 2 Nano Res. 2020 13 133–137 133–137 10.1007/s12274-019-2584-4
Lin, H. et al. Scanning tunneling spectroscopic study of monolayerT-TaS 2 and 1 T-TaSe 2. Nano Res. 13, 133–137–133–137 (2020).10.1007/s12274-019-2584-4
21. Johannes MD Mazin II Fermi surface nesting and the origin of charge density waves in metals Phys. Rev. B 2008 77 165135 165135 10.1103/PhysRevB.77.165135
Johannes, M. D. & Mazin, I. I. Fermi surface nesting and the origin of charge density waves in metals. Phys. Rev. B 77, 165135–165135 (2008).10.1103/PhysRevB.77.165135
22. Overhauser AW Exchange and correlation instabilities of simple metals Phys. Rev. 1968 167 691 698 10.1103/PhysRev.167.691
Overhauser, A. W. Exchange and correlation instabilities of simple metals. Phys. Rev. 167, 691–698 (1968).10.1103/PhysRev.167.691
23. Chen C-W Choe J Morosan E Charge density waves in strongly correlated electron systems Rep. Prog. Phys. 2016 79 084505 10.1088/0034-4885/79/8/084505 27376547
Chen, C.-W., Choe, J. & Morosan, E. Charge density waves in strongly correlated electron systems. Rep. Prog. Phys. 79, 084505 (2016).27376547 10.1088/0034-4885/79/8/084505
24. Altvater MA Charge density wave vortex lattice observed in graphene-passivated 1T-TaS2 by ambient scanning tunneling microscopy Nano Lett. 2021 21 6132 6138 10.1021/acs.nanolett.1c01655 34231367
Altvater, M. A. et al. Charge density wave vortex lattice observed in graphene-passivated 1T-TaS2 by ambient scanning tunneling microscopy. Nano Lett. 21, 6132–6138 (2021).34231367 10.1021/acs.nanolett.1c01655
25. Zhang Z Watanabe K Taniguchi T LeRoy BJ Local characterization and engineering of proximitized correlated states in graphene/NbSe2 vertical heterostructures Phys. Rev. B 2020 102 085429 10.1103/PhysRevB.102.085429
Zhang, Z., Watanabe, K., Taniguchi, T. & LeRoy, B. J. Local characterization and engineering of proximitized correlated states in graphene/NbSe2 vertical heterostructures. Phys. Rev. B 102, 085429 (2020).10.1103/PhysRevB.102.085429
26. Rahnejat KC Charge density waves in the graphene sheets of the superconductor CaC6 Nat. Commun. 2011 2 558 558 10.1038/ncomms1574 22127054
Rahnejat, K. C. et al. Charge density waves in the graphene sheets of the superconductor CaC6. Nat. Commun. 2, 558–558 (2011).22127054 10.1038/ncomms1574
27. Chen, Y. et al. Visualizing the Anomalous Charge Density Wave States in Graphene/NbSe2 Heterostructures. Adv. Mater. 32, 2003746 (2020).
28. Chen C Singh B Lin H Pereira VM Reproduction of the charge density wave phase diagram in 1T-TiSe2 exposes its excitonic character Phys. Rev. Lett. 2018 121 226602 226602 10.1103/PhysRevLett.121.226602 30547625
Chen, C., Singh, B., Lin, H. & Pereira, V. M. Reproduction of the charge density wave phase diagram in 1T-TiSe2 exposes its excitonic character. Phys. Rev. Lett. 121, 226602–226602 (2018).30547625 10.1103/PhysRevLett.121.226602
29. Zhang Evidence for a quasi-one-dimensional charge density wave in CuTe by angle-resolved photoemission spectroscopy Phys. Rev. Lett. 2018 121 206402 206402 10.1103/PhysRevLett.121.206402 30500247
Zhang et al. Evidence for a quasi-one-dimensional charge density wave in CuTe by angle-resolved photoemission spectroscopy. Phys. Rev. Lett. 121, 206402–206402 (2018).30500247 10.1103/PhysRevLett.121.206402
30. Altvater MA Observation of a topological defect lattice in the charge density wave of 1T-TaS2 Appl. Phys. Lett. 2021 119 121601 10.1063/5.0059662
Altvater, M. A. et al. Observation of a topological defect lattice in the charge density wave of 1T-TaS2. Appl. Phys. Lett. 119, 121601 (2021).10.1063/5.0059662
31. Li G Luican A Andrei EY Scanning tunneling spectroscopy of graphene on graphite Phys. Rev. Lett. 2009 102 176804 10.1103/PhysRevLett.102.176804 19518809
Li, G., Luican, A. & Andrei, E. Y. Scanning tunneling spectroscopy of graphene on graphite. Phys. Rev. Lett. 102, 176804 (2009).19518809 10.1103/PhysRevLett.102.176804
32. Li G Luican A Andrei EY Self-navigation of a scanning tunneling microscope tip toward a micron-sized graphene sample Rev. Sci. Instrum. 2011 82 073701 10.1063/1.3605664 21806183
Li, G., Luican, A. & Andrei, E. Y. Self-navigation of a scanning tunneling microscope tip toward a micron-sized graphene sample. Rev. Sci. Instrum. 82, 073701 (2011).21806183 10.1063/1.3605664
33. Qiao, S. et al. Mottness collapse in 1T-${{\mathrm{TaS_{2-x}}}Se_x}$ transition-metal dichalcogenide: An interplay between localized and itinerant orbitals. Phys. Rev. X 7, 041054 (2017).
34. Ma, L. et al. A metallic mosaic phase and the origin of Mott-insulating state in 1T-TaS2. Nat. Commun. 7, 10956 (2016).
35. Wang Z Surface-limited superconducting phase transition on 1 T-TaS2 ACS Nano 2018 12 12619 12628 10.1021/acsnano.8b07379 30403840
Wang, Z. et al. Surface-limited superconducting phase transition on 1 T-TaS2. ACS Nano 12, 12619–12628 (2018).30403840 10.1021/acsnano.8b07379
36. Yu Y-J Tuning the graphene work function by electric field effect Nano Lett. 2009 9 3430 3434 10.1021/nl901572a 19719145
Yu, Y.-J. et al. Tuning the graphene work function by electric field effect. Nano Lett. 9, 3430–3434 (2009).19719145 10.1021/nl901572a
37. Shimada T Ohuchi FS Parkinson BA Work function and photothreshold of layered metal dichalcogenides Jpn. j. Appl. Phys. 1994 33 2696 2698 10.1143/JJAP.33.2696
Shimada, T., Ohuchi, F. S. & Parkinson, B. A. Work function and photothreshold of layered metal dichalcogenides. Jpn. j. Appl. Phys. 33, 2696–2698 (1994).10.1143/JJAP.33.2696
38. Chen Y Strong correlations and orbital texture in single-layer 1T-TaSe2 Nat. Phys. 2020 16 218 224 10.1038/s41567-019-0744-9
Chen, Y. et al. Strong correlations and orbital texture in single-layer 1T-TaSe2. Nat. Phys. 16, 218–224 (2020).10.1038/s41567-019-0744-9
39. Liu Z-Y Charge transfer gap tuning via structural distortion in monolayer 1T-NbSe2 Nano Lett. 2021 21 7005 7011 10.1021/acs.nanolett.1c02348 34350759
Liu, Z.-Y. et al. Charge transfer gap tuning via structural distortion in monolayer 1T-NbSe2. Nano Lett. 21, 7005–7011 (2021).34350759 10.1021/acs.nanolett.1c02348
40. Nakata Y Robust charge-density wave strengthened by electron correlations in monolayer 1T-TaSe2 and 1T-NbSe2 Nat. Commun. 2021 12 5873 10.1038/s41467-021-26105-1 34620875
Nakata, Y. et al. Robust charge-density wave strengthened by electron correlations in monolayer 1T-TaSe2 and 1T-NbSe2. Nat. Commun. 12, 5873 (2021).34620875 10.1038/s41467-021-26105-1
41. Ruan W Evidence for quantum spin liquid behaviour in single-layer 1T-TaSe2 from scanning tunnelling microscopy Nat. Phys. 2021 17 1154 1161 10.1038/s41567-021-01321-0
Ruan, W. et al. Evidence for quantum spin liquid behaviour in single-layer 1T-TaSe2 from scanning tunnelling microscopy. Nat. Phys. 17, 1154–1161 (2021).10.1038/s41567-021-01321-0
42. Nakata Y Monolayer 1T-NbSe2 as a mott insulator NPG Asia Mater. 2016 8 e321 e321 10.1038/am.2016.157
Nakata, Y. et al. Monolayer 1T-NbSe2 as a mott insulator. NPG Asia Mater. 8, e321–e321 (2016).10.1038/am.2016.157
43. Liu M Monolayer 1T-NbSe2 as a 2D-correlated magnetic insulator Sci. Adv. 2011 7 eabi6339 10.1126/sciadv.abi6339
Liu, M. et al. Monolayer 1T-NbSe2 as a 2D-correlated magnetic insulator. Sci. Adv. 7, eabi6339 (2011).10.1126/sciadv.abi6339
44. Kresse G Ab initio molecular dynamics for liquid metals J. Non Cryst. Solids 1995 192-193 222 229 10.1016/0022-3093(95)00355-X
Kresse, G. Ab initio molecular dynamics for liquid metals. J. Non Cryst. Solids 192-193, 222–229 (1995).10.1016/0022-3093(95)00355-X
45. Kresse G Furthmüller J Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set Comput. Mater. Sci. 1996 6 15 50 10.1016/0927-0256(96)00008-0
Kresse, G. & Furthmüller, J. Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set. Comput. Mater. Sci. 6, 15–50 (1996).10.1016/0927-0256(96)00008-0
46. Perdew JP Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation Phys. Rev. B 1992 46 6671 6687 10.1103/PhysRevB.46.6671
Perdew, J. P. et al. Atoms, molecules, solids, and surfaces: Applications of the generalized gradient approximation for exchange and correlation. Phys. Rev. B 46, 6671–6687 (1992).10.1103/PhysRevB.46.6671
47. Perdew JP Wang Y Pair-distribution function and its coupling-constant average for the spin-polarized electron gas Phys. Rev. B 1992 46 12947 12954 10.1103/PhysRevB.46.12947
Perdew, J. P. & Wang, Y. Pair-distribution function and its coupling-constant average for the spin-polarized electron gas. Phys. Rev. B 46, 12947–12954 (1992).10.1103/PhysRevB.46.12947
48. Darancet, P., Millis, A. J. & Marianetti, C. A. Three-dimensional metallic and two-dimensional insulating behavior in octahedral tantalum dichalcogenides. Phys. Revi. B 90, 045134 (2014).
49. Momma, K. & Izumi, F. VESTA 3 for three-dimensional visualization of crystal, volumetric and morphology data. J. Appl. Cryst. 44, 1272–1276 (2011).
