Interplay of charge density waves, disorder, and superconductivity in 2-TaSe2 elucidated by NMR
Abstract
Single crystals of pristine and 6% Pd-intercalated 2-TaSe2 have been studied by means of 77Se nuclear magnetic resonance (NMR). The temperature dependence of the 77Se spectrum, with an unexpected line narrowing upon Pd intercalation, unravels the presence of correlated local lattice distortions far above the transition temperature of the charge density wave (CDW) order, thereby supporting a strong-coupling CDW mechanism in 2-TaSe2. While, the Knight shift data suggest that the incommensurate CDW transition involves a partial Fermi surface gap opening. As for spin dynamics, the 77Se spin-lattice relaxation rate as a function of temperature shows that a pseudogap behavior dominates the low-energy spin excitations even within the CDW phase, and gets stronger along with superconductivity in the Pd-6% sample. We discuss that CDW fluctuations may be responsible for the pseudogap as well as superconductivity, although the two phenomena are unlikely to be directly linked each other.
1 Introduction
Charge density wave (CDW) order in two dimensions, together with its relationship to superconductivity, has been a central issue in the layered transition metal dichalcogenides (TMDs) [1, 2, 3], and even more so owing to remarkable similarities with high- copper-oxide superconductors (cuprates) [4]. Electronic phase diagrams in many metallic TMDs suggest that the emergence or enhancement of superconductivity is closely related to CDW order, although the relationship between the two phenomena remains unclear. Another particularly interesting feature in TMDs is the presence of a pseudogap regime [5, 6, 7, 8] which involves strange-metal behavior in the normal state. The origin of the pseudogap is often ascribed to a CDW instability [6, 9], an argument also put forward for cuprates [10, 11, 12]. Hence the deep understanding of the nature and origin of the CDWs and the pseudogap in TMDs may provide vital clues to understanding the mechanism for high-temperature superconductivity.
2-TaSe2 is one of the intensely investigated TMDs, as it develops a series of fascinating phases: an unusual metallic state with a pseudogap at high temperatures is followed by an incommensurate CDW (iCDW) transition at K, a lock-in transition into the commensurate CDW (cCDW) at K, and a superconducting (SC) transition at K [13]. The iCDW state may be stabilized by CDW defects termed discommensurations which separate commensurate regions [14, 15, 16, 17]. In this picture, the lock-in transition could be accounted for by the disappearance of discommensurations. However, no consensus has been reached on whether the CDW order is of weak-coupling nature, driven by Fermi-surface nesting [18, 19, 20], or of more local strong-coupling character [21, 22, 23].
Since the two distinct CDW transitions successively occur at much higher temperatures than , it may be possible to establish a relation between CDW and superconductivity by tuning control parameters such as pressure or doping. Recently, it has been demonstrated that Pd intercalation leads to the dramatic enhancement of up to K near an optimal Pd content of % at which the cCDW completely vanishes while the iCDW transition remains robust [24], as shown in Fig. 1(a). Being motivated by the strong effect of Pd intercalation on superconductivity and the commensurate CDW in 2-PdxTaSe2, we carried out 77Se nuclear magnetic resonance (NMR) in pristine and 6% Pd-intercalated 2-TaSe2. Our NMR data suggest that the major driving force for the CDW formation is strong electron-phonon coupling (EPC). They further suggest that Pd intercalation introduces both changes to the electronic structure and random pinning centers, together being responsible for the strong smearing and suppression of the cCDW transition as well as for the strengthening of the pseudogap behavior on top of the substantial increase of .
2 Experimental details
Single crystals of 2-PdxTaSe2 [ (pristine) and 0.06 (Pd-6%)] were grown by the chemical vapor transport method as described in detail in Ref. [24]. It has been confirmed that the Pd intercalation does not alter the structure of the pristine compound, as drawn in the inset of Fig. 1(a). In-plane resistivity was measured by the conventional four probe technique using a conductive silver epoxy in PPMSTM (Quantum Design). Uniform magnetic susceptibility was measured at 5 T applied along the plane in MPMSTM (Quantum Design). Because of the small size of the crystals, the magnetization signal was barely detected in one single crystal so that we stacked 15 single crystals in parallel (a total mass of mg) to obtain . The temperature dependences of and for and are presented in Figs. 1 (b) and (c).

77Se (nuclear spin ) NMR was carried out in pristine and 6% Pd-intercalated 2-TaSe2 single crystals () at an external magnetic field of 15 T, and in the range of temperature 4.2–300 K. The samples were oriented using a goniometer for the accurate alignment along the external field. The 77Se NMR spectra were acquired by a standard spin-echo technique with a typical pulse length 2–3 s. The nuclear spin-lattice relaxation rate was obtained by fitting the recovery of the nuclear magnetization after a saturating pulse to the following function: where is a fitting parameter that is ideally unity. In this study, due to the very weak 77Se NMR signal intensity in tiny single crystals, we have focused on the data for for which the 77Se linewidth is much smaller than that for , as the anisotropy of the physical properties is not our main concern of this work.

3 Results
3.1 NMR linewdith
Figure 2(a) shows 77Se NMR spectra in the pristine and Pd-6% samples as a function of temperature measured at 15 T applied parallel to the plane. As shown in Fig. 2(b), we find that the full width at half maximum (FWHM) of the spectrum, which is a measure of spatial inhomogeneity, shows nearly identical temperature dependence for both samples. A remarkable and unexpected finding is that the FWHM for the Pd-6% sample is considerably reduced compared to the pristine one. Such a narrowing of an NMR line with doping is extremely unusual, because any dopants inevitably introduce chemical disorder typically leading to a much larger NMR line broadening. In fact, the in-plane resistivity significantly increases in the Pd-6% sample, Fig. 1(b), reflecting the increased lattice disorder induced by Pd intercalation. Also, the uniform magnetic susceptibility , Fig. 1(c), reveals a much stronger Curie-like tail for the Pd-6% sample than that for the pristine one, evidencing the increase of localized paramagnetic impurities by Pd intercalation.
We therefore arrive at the conclusion that the spatial inhomogeneity probed by NMR in pristine 2-TaSe2 is largely intrinsic, at least far beyond that originating from impurities. Because of the strong CDW instability in these materials, we think that the intrinsic inhomogeneity arises from local lattice distortions related to CDW phenomena whose amplitude is distributed in space. If so, the fact that the FWHM is reduced in the whole temperature range investigated, Fig. 2(b), suggests that local lattice distortions exist even at temperatures far above , which in turn suggests that the lattice distortions at high temperatures may evolve into the long-ranged CDW state when cooling to below . This is clearly inconsistent with a conventional, i.e., weak-coupling CDW mechanism, in which lattice distortions develop only below the CDW transition temperature. That is, our data indicate that the CDW transition in 2-TaSe2 is primarily driven by local EPC [21, 22, 23] leading to a strong-coupling CDW mechanism with a large fluctuation regime [3, 25, 26], rather than by electronic Fermi-surface nesting of Peierls type [18, 19].
The question may arise here as to how local lattice distortions become static (at least on the NMR time scale) at , as they should rapidly fluctuate in time, and how their amplitude is suppressed by external disorder. One could argue that the presence of any lattice defects will lead to pinning, such that distortions become inhomogeneous but static, realizing local patches of CDW order. Due to the low density of pinning centers, the intrinsic periodicity of the CDW in pristine 2-TaSe2 could be well preserved resulting in the sharp CDW transitions. However, the density of random-field pinning centers rapidly increase upon intercalating Pd, and may strongly compete with the CDW periodicity. This could result in an overall suppression of the distortion amplitude as well as strongly inhomogeneous distortion patterns, which naturally accounts for the nearly indefinable cCDW transition near 6% Pd content [24].
3.2 Knight shift
The Knight shift , which is equivalent to the local spin susceptibility, of the 77Se spectra is presented in Fig. 2(c) as a function of temperature. The weak temperature dependence of reflects the metallic nature of 2-PdxTaSe2, as is independent of temperature in a Pauli metal. Nevertheless, is slightly enhanced with lowering temperature roughly down to , indicating that there is a finite non-Pauli contribution to the total spin susceptibility. For direct comparison, the bulk susceptibilities are drawn as solid curves in Fig. 2(c). For the pristine sample, except a Curie-like upturn at low temperatures, is quite consistent with . The linear relationship between and above permits us to estimate the hyperfine coupling constant, kOe/, and the -independent shift %, using the relation . For the Pd-6% sample, turns out to be governed by a strong Curie-like divergence even up to room temperature overwhelming the Pauli contribution, in contrast to which shows the almost same temperature dependence as the pristine one. This proves that the Curie-like tails indeed arise from localized impurity moments.
Interestingly, looks to be suppressed below for both samples, although there are seldom sufficient data. In great support of this, we note that drops at for the pristine sample. Also, for the Pd-6% sample, a kink was observed at [see Fig. 1(c)], which may be identified as a maximum of if the strong Curie background is relatively weak as for the pristine one. The decrease of the spin susceptibility below suggests that the iCDW transition involves a partial gap opening at the Fermi level. Specifically, the Fermi surface nesting may partially contribute to the CDW transition [27], although it is unlikely the main driving force [21] as discussed above.
3.3 Spin dynamics
We now turn to low-energy spin dynamics probed via the spin-lattice relaxation rate which is very sensitive to quasiparticle excitations near the Fermi level in metals. Figures 3(a) and (b) present and , respectively, as a function of temperature in 2-PdxTaSe2 ( and 0.06) measured at T parallel and perpendicular to . In an ordinary metal, one expects that the relation , or , holds [28]. However, a linear dependence of was not observed, which is better shown in vs. plot in Fig. 3(b). The rapid decrease of with decreasing implies a progressive reduction of low-energy degrees of freedom, or the existence of a pseudogap. In the Pd-6% sample, decreases faster upon cooling than in the pristine one, i.e., the pseudogap behavior becomes stronger. We also measured the data at low temperatures for to check the anisotropy of the pseudogap. The results, Fig. 3(b), show that, while is weakly anisotropic, the pseudogap behavior is nearly isotropic being similarly strengthened by Pd intercalation for both field directions.
Notably, the pseudogap behavior governs the low-energy spin dynamics both above and below , with no significant change across . This suggests that the Fermi surface gapping caused by the iCDW transition observed in the Knight shift measurement, Fig. 2(c), is small compared to the pseudogap. Furthermore, the strengthening of the pseudogap behavior with Pd intercalation contrasts with the weakening of lattice distortions and the CDW. These observations suggest that the pseudogap can be directly linked neither to Fermi surface nesting [8, 20], nor to static lattice distortions above the transition as expected in a strong-coupling CDW scenario [25]. Possibly, the pseudogap is instead related to dynamically fluctuating CDW; such fluctuations of the CDW may also be responsible for the strange metallic behavior characterized by the linear variation of resistivity with temperature [7, 24].

4 Discussion
We now try to synthesize our NMR findings and the phase diagram of Fig. 1(a) into a coherent picture. Our data suggest that the intrinsically strong EPC generates correlated local lattice distortions in the normal state, ruling much of the phenomena in 2-TaSe2. Pd dopants reduce the distortion amplitude and introduce strong random-field-type pinning which destroys the cCDW state. (Note that we reinterpret the data in Ref. [24] so that the cCDW transitions above are not simply suppressed, but also ill-defined due to the increased density of pinning centers which may hinder the lock-in transition to the cCDW state.) By contrast, the superconducting is dramatically enhanced by increased Pd dopants, in the light of the weakly doping-dependent iCDW transition.
Based on the very different role played by disorder for the CDW transitions and superconductivity, we propose that it may be the commensurability of the CDW that competes with superconductivity, rather than CDW ordering itself. Conversely, superconductivity and incommensurate CDW may not be mutually exclusive in nature [29]. Indeed, other superconducting TMDs are also well understood in line with this reasoning. The isostructural 2-NbSe2 and 2-NbS2 [1] in which the cCDW transition is absent, exhibit K and K, respectively, much higher than K and K of the Ta counterparts in which the cCDW is present111Although the lock-in transition has not been reported for 2-TaS2, the local charge density at K is found to be commensurate by NMR [30]. In this respect, the seemingly competing relationship between superconductivity and the CDW in this material, as discussed in Refs. [31, 32], may require a further consideration of the commensurability of the CDW. Furthermore, suppressing in 2-NbSe2 [33] has a negligible effect on , in agreement with a weak correlation between iCDW and superconductivity as in Fig. 1(a). In contrast, suppressing either by doping or by external pressure in 1 polytypes mostly induces a substantial increase of — 1-TaSe2 [34], 1-TaS2 [35, 36], 1-TiSe2 [37, 38], and 1-VSe2 [39], as it precisely did in 2-TaSe2.
Lastly, the fact that both the pseudogap behavior and superconductivity are boosted by Pd intercalation may suggest that the two phenomena are linked in that CDW fluctuations contribute to both. In this case, the very different temperature scales of pseudogap and superconductivity in 2-TaSe2 suggest that the pseudogap in TMDs is a phenomenon distinct from superconductivity, rather than a state of preformed pairs which is often discussed for cuprates [40, 41, 42].
5 Summary
We carried out 77Se NMR measurements in pristine and 6% Pd-intercalated 2-TaSe2 single crystals. Our NMR results, combined with resistivity and uniform magnetic susceptibility measurements, strongly suggest that correlated local lattice distortions exist in the normal state as precursor of CDW formation, and thus a strong-coupling CDW mechanism arising from electron-phonon coupling is likely at play in 2-TaSe2. We argue that the high density of random pinning disorder and modifications to the band-structure caused by Pd intercalation reduce the lattice distortion amplitude and destroy the intrinsic periodicity of the CDW, which accounts for the strong smearing and suppression of the incommensurate-commensurate CDW transition. In contrast, the pseudogap behavior is much enhanced by Pd intercalation, as is the superconducting transition . Based on the subtle effects of external disorder, we propose that quenching the commensurability of the CDW may be a crucial factor for the enhancement of .
Acknowledgments
This work was supported by the National Research Foundation of Korea (NRF) grant funded by the Korea government (MSIT) (NRF-2020R1A2C1003817 and NRF-2019R1A2C2090648) and by the Ministry of Education (2021R1A6C101B418). MV and BB acknowledge support from the Deutsche Forschungsgemeinschaft through SFB 1143 (project-id 247310070) and the Würzburg-Dresden Cluster of Excellence ct.qmat (EXC 2147, project-id 390858490).
References
References
- [1] J.A. Wilson, F.J. Di Salvo, and S. Mahajan. Charge-density waves and superlattices in the metallic layered transition metal dichalcogenides. Adv. Phys., 24:117–201, 1975.
- [2] A. H. Castro Neto. Charge Density Wave, Superconductivity, and Anomalous Metallic Behavior in 2D Transition Metal Dichalcogenides. Phys. Rev. Lett., 86:4382–4385, 2001.
- [3] K. Rossnagel. On the origin of charge-density waves in select layered transition-metal dichalcogenides. J. Phys.: Condens. Matter, 23:213001, 2011.
- [4] C.-W. Chen, J. Choe, and E. Morosan. Charge density waves in strongly correlated electron systems. Rep. Prog. Phys., 79:084505, 2016.
- [5] V. Vescoli, L. Degiorgi, H. Berger, and L. Forró. Dynamics of Correlated Two-Dimensional Materials: The -TaSe2 Case. Phys. Rev. Lett., 81:453–456, Jul 1998.
- [6] R.A. Klemm. Striking similarities between the pseudogap phenomena in cuprates and in layered organic and dichalcogenide superconductors. Physica C, 341-348:839–842, 2000.
- [7] B. Ruzicka, L. Degiorgi, H. Berger, R. Gaál, and L. Forró. Charge Dynamics of 2H-TaSe2 along the Less-Conducting -Axis. Phys. Rev. Lett., 86:4136–4139, 2001.
- [8] S. V. Borisenko, A. A. Kordyuk, A. N. Yaresko, V. B. Zabolotnyy, D. S. Inosov, R. Schuster, B. Büchner, R. Weber, R. Follath, L. Patthey, and H. Berger. Pseudogap and charge density waves in two dimensions. Phys. Rev. Lett., 100:196402, May 2008.
- [9] M. Bovet, D. Popović, F. Clerc, C. Koitzsch, U. Probst, E. Bucher, H. Berger, D. Naumović, and P. Aebi. Pseudogapped Fermi surfaces of and A charge density wave effect. Phys. Rev. B, 69:125117, 2004.
- [10] I. Eremin, M. Eremin, S. Varlamov, D. Brinkmann, M. Mali, and J. Roos. Spin susceptibility and pseudogap in YBa2Cu4O8: An approach via a charge-density-wave instability. Phys. Rev. B, 56:11305–11311, 1997.
- [11] J.-X. Li, C.-Q. Wu, and D.-H. Lee. Checkerboard charge density wave and pseudogap of high- cuprate. Phys. Rev. B, 74:184515, 2006.
- [12] W. D. Wise, M. C. Boyer, K. Chatterjee, T. Kondo, T. Takeuchi, H. Ikuta, Y. Wang, and E. W. Hudson. Charge-density-wave origin of cuprate checkerboard visualized by scanning tunnelling microscopy. Nat. Phys., 4:696–699, 2008.
- [13] D. E. Moncton, J. D. Axe, and F. J. DiSalvo. Neutron scattering study of the charge-density wave transitions in -TaSe2 and -NbSe2. Phys. Rev. B, 16:801–819, 1977.
- [14] W. L. McMillan. Landau theory of charge-density waves in transition-metal dichalcogenides. Phys. Rev. B, 12:1187–1196, 1975.
- [15] W. L. McMillan. Theory of discommensurations and the commensurate-incommensurate charge-density-wave phase transition. Phys. Rev. B, 14:1496–1502, Aug 1976.
- [16] B. H. Suits, S. Couturié, and C. P. Slichter. Confirmation of McMillan’s Concept of “Discommensurations”. Phys. Rev. Lett., 45:194–197, 1980.
- [17] B. H. Suits, S. Couturié, and C. P. Slichter. NMR test of McMillan’s concept of discommensurations in -TaSe2. Phys. Rev. B, 23:5142–5151, 1981.
- [18] Th. Straub, Th. Finteis, R. Claessen, P. Steiner, S. Hüfner, P. Blaha, C. S. Oglesby, and E. Bucher. Charge-Density-Wave Mechanism in -NbSe2: Photoemission Results. Phys. Rev. Lett., 82:4504–4507, 1999.
- [19] Y. W. Li, J. Jiang, H. F. Yang, D. Prabhakaran, Z. K. Liu, L. X. Yang, and Y. L. Chen. Folded superstructure and degeneracy-enhanced band gap in the weak-coupling charge density wave system -. Phys. Rev. B, 97:115118, 2018.
- [20] A. Chikina, A. Fedorov, D. Bhoi, V. Voroshnin, E. Haubold, Y. Kushnirenko, K. H. Kim, and S. Borisenko. Turning charge-density waves into cooper pairs. npj Quantum Mater., 5(1):22, 2020.
- [21] M. D. Johannes and I. I. Mazin. Fermi surface nesting and the origin of charge density waves in metals. Phys. Rev. B, 77:165135, 2008.
- [22] J. Dai, E. Calleja, J. Alldredge, X. Zhu, L. Li, W. Lu, Y. Sun, T. Wolf, H. Berger, and K. McElroy. Microscopic evidence for strong periodic lattice distortion in two-dimensional charge-density wave systems. Phys. Rev. B, 89:165140, 2014.
- [23] V. Petkov, K. Chapagain, S. Shastri, and Y. Ren. Genesis of the periodic lattice distortions in the charge density wave state of -TaSe2. Phys. Rev. B, 101:121114, 2020.
- [24] D. Bhoi, S. Khim, W. Nam, B. S. Lee, C. Kim, B.-G. Jeon, B. H. Min, S. Park, and K. H. Kim. Interplay of charge density wave and multiband superconductivity in -PdxTaSe2. Sci. Rep., 6:24068, 2016.
- [25] W. L. McMillan. Microscopic model of charge-density waves in -TaSe2. Phys. Rev. B, 16:643–650, 1977.
- [26] X. Zhu, J. Guo, J. Zhang, and E. W. Plummer. Misconceptions associated with the origin of charge density waves. Adv. Phys. -X, 2:622–640, 2017.
- [27] D. S. Inosov, V. B. Zabolotnyy, D. V. Evtushinsky, A. A. Kordyuk, B. Büchner, R. Follath, H. Berger, and S .V. Borisenko. Fermi surface nesting in several transition metal dichalcogenides. New J. Phys., 10:125027, 2008.
- [28] G. C. Carter, L. H. Bennett, and D. J. Kahan. Metallic shift in NMR. Pergamon, New York, 1977.
- [29] C. Wen, Y. Xie, Y. Wu, S. Shen, P. Kong, H. Lian, J. Li, H. Xing, and S. Yan. Impurity-pinned incommensurate charge density wave and local phonon excitations in -NbS2. Phys. Rev. B, 101:241404, 2020.
- [30] H. Nishihara, G.A. Scholz, M. Naito, R.F. Frindt, and S. Tanaka. NMR of 181Ta in -TaS2 and -TaSe2 - observation of locally commensurate CDW. J. Magn. Magn. Mater., 31-34:717–718, 1983.
- [31] H. Liu, S. Huangfu, X. Zhang, H. Lin, and A. Schilling. Superconductivity and charge density wave formation in lithium-intercalated -TaS2. Phys. Rev. B, 104:064511, 2021.
- [32] S. Xu, J. Gao, Z. Liu, K. Chen, P. Yang, S. Tian, C. Gong, J. Sun, M. Xue, J. Gouchi, X. Luo, Y. Sun, Y. Uwatoko, H. Lei, B. Wang, and J. Cheng. Effects of disorder and hydrostatic pressure on charge density wave and superconductivity in -TaS2. Phys. Rev. B, 103:224509, 2021.
- [33] M. Leroux, I. Errea, M. Le Tacon, S.-M. Souliou, G. Garbarino, L. Cario, A. Bosak, F. Mauri, M. Calandra, and P. Rodière. Strong anharmonicity induces quantum melting of charge density wave in - under pressure. Phys. Rev. B, 92:140303, 2015.
- [34] Y. Liu, D. F. Shao, L. J. Li, W. J. Lu, X. D. Zhu, P. Tong, R. C. Xiao, L. S. Ling, C. Y. Xi, L. Pi, H. F. Tian, H. X. Yang, J. Q. Li, W. H. Song, X. B. Zhu, and Y. P. Sun. Nature of charge density waves and superconductivity in -TaSe2-xTex. Phys. Rev. B, 94:045131, 2016.
- [35] B. Sipos, A. F. Kusmartseva, A. Akrap, H. Berger, L. Forró, and E. Tutiš. From Mott state to superconductivity in -TaS2. Nature Mater., 7:960–965, 2008.
- [36] R. Ang, Y. Miyata, E. Ieki, K. Nakayama, T. Sato, Y. Liu, W. J. Lu, Y. P. Sun, and T. Takahashi. Superconductivity and bandwidth-controlled Mott metal-insulator transition in -TaS2-xSex. Phys. Rev. B, 88:115145, 2013.
- [37] E. Morosan, H. W. Zandbergen, B. S. Dennis, J. W. G. Bos, Y. Onose, T. Klimczuk, A. P. Ramirez, N. P. Ong, and R. J. Cava. Superconductivity in CuxTiSe2. Nature Phys., 2:544–550, 2006.
- [38] A. F. Kusmartseva, B. Sipos, H. Berger, L. Forró, and E. Tutiš. Pressure Induced Superconductivity in Pristine -TiSe2. Phys. Rev. Lett., 103:236401, 2009.
- [39] S. Sahoo, U. Dutta, L. Harnagea, A. K. Sood, and S. Karmakar. Pressure-induced suppression of charge density wave and emergence of superconductivity in -VSe2. Phys. Rev. B, 101(1):014514, jan 2020.
- [40] V. J. Emery and S. A. Kivelson. Importance of phase fluctuations in superconductors with small superfluid density. Nature, 374:434–437, 1995.
- [41] P. W. Anderson. Physics of the Resonating Valence Bond (Pseudogap) State of the Doped Mott Insulator: Spin-Charge Locking. Phys. Rev. Lett., 96:017001, 2006.
- [42] Y. Kohsaka, C. Taylor, P. Wahl, A. Schmidt, J. Lee, K. Fujita, J. W. Alldredge, K. McElroy, J. Lee, H. Eisaki, S. Uchida, D.-H. Lee, and J. C. Davis. How Cooper pairs vanish approaching the Mott insulator in Bi2Sr2CaCu2O8+δ. Nature, 454:1072–1078, 2008.