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

51631
10.1038/s41467-024-51631-z
Article
Nematic Ising superconductivity with hidden magnetism in few-layer 6R-TaS2
Liu Shao-Bo 1
Tian Congkuan 12
Fang Yuqiang 3
Rong Hongtao 4
http://orcid.org/0000-0001-9549-8136
Cao Lu 5
Wei Xinjian 2
Cui Hang 1
Chen Mantang 1
http://orcid.org/0000-0002-0758-1763
Chen Di 2
Song Yuanjun 2
Cui Jian 2
Li Jiankun 2
Guan Shuyue 1
http://orcid.org/0000-0001-8705-6055
Jia Shuang 1
http://orcid.org/0000-0002-3930-8294
Chen Chaoyu 4
http://orcid.org/0000-0002-1595-8807
He Wenyu 6
http://orcid.org/0000-0003-0526-5473
Huang Fuqiang huangfq@sjtu.edu.cn

3
http://orcid.org/0000-0003-4842-4921
Jiang Yuhang yuhangjiang@ucas.ac.cn

5
http://orcid.org/0000-0002-9034-3642
Mao Jinhai 7
Xie X. C. 189
http://orcid.org/0000-0003-0501-6290
Law Kam Tuen 10
http://orcid.org/0000-0002-9485-1759
Chen Jian-Hao chenjianhao@pku.edu.cn

12911
1 https://ror.org/02v51f717 grid.11135.37 0000 0001 2256 9319 International Center for Quantum Materials, School of Physics, Peking University, Beijing, China
2 https://ror.org/04nqf9k60 grid.510904.9 0000 0004 9362 2406 Beijing Academy of Quantum Information Sciences, Beijing, China
3 https://ror.org/0220qvk04 grid.16821.3c 0000 0004 0368 8293 School of Materials Science and Engineering, Shanghai Jiao Tong University, Shanghai, China
4 https://ror.org/049tv2d57 grid.263817.9 0000 0004 1773 1790 Shenzhen Institute for Quantum Science and Engineering and Department of Physics, Southern University of Science and Technology, Shenzhen, China
5 https://ror.org/05qbk4x57 grid.410726.6 0000 0004 1797 8419 College of Materials Science and Optoelectronic Technology, Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing, China
6 https://ror.org/030bhh786 grid.440637.2 0000 0004 4657 8879 School of Physical Science and Technology, ShanghaiTech University, Shanghai, China
7 https://ror.org/05qbk4x57 grid.410726.6 0000 0004 1797 8419 CAS Center for Excellence in Topological Quantum Computation, University of Chinese Academy of Sciences, Beijing, China
8 https://ror.org/013q1eq08 grid.8547.e 0000 0001 0125 2443 Institute for Nanoelectronic Devices and Quantum Computing, Fudan University, Shanghai, China
9 grid.59053.3a 0000000121679639 Hefei National Laboratory, Hefei, China
10 grid.24515.37 0000 0004 1937 1450 Department of Physics, Hong Kong University of Science and Technology, Hong Kong, China
11 https://ror.org/02v51f717 grid.11135.37 0000 0001 2256 9319 Key Laboratory for the Physics and Chemistry of Nanodevices, Peking University, Beijing, China
31 8 2024
31 8 2024
2024
15 75699 2 2024
14 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/.
In van der Waals heterostructures (vdWHs), the manipulation of interlayer stacking/coupling allows for the construction of customizable quantum systems exhibiting exotic physics. An illustrative example is the diverse range of states of matter achieved through varying the proximity coupling between two-dimensional (2D) quantum spin liquid (QSL) and superconductors within the TaS2 family. This study presents a demonstration of the intertwined physics of spontaneous rotational symmetry breaking, hidden magnetism, and Ising superconductivity (SC) in the three-fold rotationally symmetric, non-magnetic natural vdWHs 6R-TaS2. A distinctive phase emerges in 6R-TaS2 below a characteristic temperature (T*) of approximately 30 K, which is characterized by a remarkable set of features, including a giant extrinsic anomalous Hall effect (AHE), Kondo screening, magnetic field-tunable thermal hysteresis, and nematic magneto-resistance. At lower temperatures, a coexistence of nematicity and Kondo screening with Ising superconductivity is observed, providing compelling evidence of hidden magnetism within a superconductor. This research not only sheds light on unexpected emergent physics resulting from the coupling of itinerant electrons and localized/correlated electrons in natural vdWHs but also emphasizes the potential for tailoring exotic quantum states through the manipulation of interlayer interactions.

6R-TaS2 is a natural van der Waals heterostructure formed by 1H- and 1T-phase TaS2 layers, which can individually exhibit Ising superconductivity, correlated states and charge density waves. Here, the authors show experimental evidence of emergent nematic Ising superconductivity with simultaneous hidden magnetism (extrinsic anomalous Hall effect and Kondo screening) in 6R-TaS2 under 30 K.

Subject terms

Superconducting properties and materials
Magnetic properties and materials
Electronic properties and materials
Phase transitions and critical phenomena
https://doi.org/10.13039/501100001809 National Natural Science Foundation of China (National Science Foundation of China) 11934001 92265106 11774010 11921005 Chen Jian-Hao Shanghai Rising-Star Program (23QA1410700)issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

2D transition metal dichalcogenides (TMDs) have attracted considerable interest as fascinating quantum materials, exhibiting intricate interplay between charge, spin, and lattice degrees of freedom with strong spin-orbit coupling (SOC)1,2. This unique combination leads to the emergence of multiple intertwined orders with similar energies, giving rise to exotic physics that holds great potential for next-generation quantum devices1–3. In particular, the assembly of TMDs into vdWHs offers exciting prospects for realizing unexpected electronic ground states that are not present in the individual constituent layers2,4–8. Some of these states include nematicity, unconventional Ising superconductivity and chiral spin liquids4,6,9–11.

Ising superconductivity in TMDs, such as 2H-MoS2, is generally in-plane isotropic12. A few exceptions, such as few-layer 1Td-MoTe2, 2H-NbSe2 and 2M-WS2, possess in-plane anisotropic Ising-like superconductivity13–15. The mechanisms underlying this anisotropy have been under debate, with proposals such as anisotropic SOC13, competing superconducting order parameters14 and spin-orbit-parity coupling15. Hence, data from additional in-plane anisotropic Ising superconductors would help to unravel this puzzle. Another important physical effect, the AHE, typically found in magnetic systems with spin-dependent scattering16, in topological materials with nontrivial Berry phase17, or in systems with non-coplanar spin structures18,19, could be an indication of time-reversal symmetry breaking (TRSB). Recent reports of giant extrinsic AHE and TRSB in non-magnetic kagome metals AV3Sb520–23 and ScV6Sn624,25 suggest the involvement of chiral charge order induced loop currents. Thus, unveiling AHE in non-magnetic systems is of great importance to advancing scientific understanding in this field.

The non-magnetic TMDs 6R-TaS2 represents a promising natural vdWH26,27, with alternating 2D layers of trigonal-prismatic 1H-TaS2 and octahedral 1T-TaS2. Single layer 1H-TaS2 is shown to exhibit SOC-driven Ising superconductivity28 with a conventional s-wave order parameter29–31, while single layer 1T-TaS2 is proven to be a correlated Mott insulator hosting a gapless QSL ground state32,33. Below 200 K, single layer 1T-TaS2 undergoes a13×13 charge-density wave (CDW) transition, leading to the formation of a “start-of-David” deformation of the Ta atoms with localized electrons on an emergent frustrated triangular superlattice32–34. The electron spins of this superlattice are devoid of magnetic ordering due to strong frustration and quantum fluctuations, creating a QSL5,34 (Fig. 1b). Thus, differently stacked 1H and 1 T layers in the TaS2 family become a fruitful source of exotic physics. Recent spectroscopic studies have revealed Kondo coupling between the itinerant (1H) and the localized (1 T) electrons in MBE grown single layer 1H/1 T heterostructures of TaS2 and TaSe25,8. However, in another natural vdWHs 4Hb-TaS2, Kondo coupling and nematic Ising superconductivity is absent deep in the superconducting state6,35–38. Furthermore, transport signature of intertwined nematicity and local time-reversal symmetry breaking in 1H/1 T TaS2 vdWHs has been elusive so far.Fig. 1 Natural van der Waals heterostructures (vdWHs), quantum spin liquid (QSL) candidate and nematic Ising superconductivity (SC) in few-layer 6R-TaS2.

a Crystal structure of 6R-TaS2 in three different angle views, exhibiting three-fold rotational (C3) symmetry along the c axis with broken inversion symmetry. The blue dashed lines represent a triangular lattice. The black dashed boxes represent a unit cell. b Schematics of a 1 T/1H bilayer comprising one-third of the 6R-TaS2 unit cell. The localized moments in the 1 T layer interacts with the Ising superconducting 1H layer through Kondo coupling (JK) or proximity effect. c Resistivity vs. T at zero magnetic field for 6R-TaS2 device #1 for T from 0.3 K to 360 K. Inset shows the zoom-in view of the superconducting transition. Here, TCDW and Tconset represents the charge-density-waves transition and superconducting onset transition temperatures, respectively. d Polar plot of angular-dependent normalized in-plane resistance Rab/RN for device #1 under various T near the Tconset at B// = 5 T, where the magnetic field angle α with respect to x axis ranges from 0° to 360°. Here, B// represents the in-plane magnetic field, and RN represents the normal state resistance at temperatures just above the superconducting transition. e Normalized in-plane upper critical field Hc2///Hp vs. reduced temperature T/Tc along the in-plane directions of α = 90° (solid dots) and 0° (empty dots). The brown dashed line denotes the normalized Pauli limit Hp.

In this work, we provide solid experimental evidence of the emergence of a hidden order in 6R-TaS2 thin flakes below a characteristic temperature T* ~ 30 K. Such hidden order simultaneously breaks local time-reversal symmetry and in-plane rotational symmetry, and persists below the superconducting transition temperature. The coexistence of the hidden order with Ising superconductivity produces strongly nematic Ising superconductivity with simultaneous Kondo screening in the 1T-layer. The origin of the nematic Ising superconductivity is narrowed down to anisotropic Ising SOC13, unconventional mixed pairing channels14, or nematicity from the normal state; while the hidden magnetism is attributed to proximity coupling between 1H and 1 T TaS2 layers4,39. Our study provides a platform in the exploration of intertwined physics of nematicity, hidden magnetism and unconventional Ising superconductivity in non-magnetic TMDs vdWHs. A summary of the variety of interesting properties of the TaS2 family with different interlayer coupling can be found in Supplementary Fig. 1.

Results

Nematic Ising superconductivity in thin flake 6R-TaS2

The crystal structure of 6R-TaS2 consists of alternating 1H-TaS2 and 1T-TaS2 layers, which exhibits three-fold rotational (C3) symmetry along the c axis with broken inversion symmetry26 (Fig. 1a). The magneto-transport properties of superconducting and normal state 6R-TaS2 were studied in various sample thicknesses. The device fabrication process and characterization can be found in Methods and Supplementary Figs. 2–6. Three 6R-TaS2 thin flake devices, device #1 (with a thickness of 43.5 nm), device #2 (20.7 nm), and device #3 (7.9 nm), as well as bulk crystals, were investigated.

We found that all the thin flake devices exhibit 2D nematic Ising superconductivity. Data from device #1 is shown in Fig. 1c–e as well as in Supplementary Figs. 7 and 8. Similar data from devices #2 and #3 can be found in Supplementary Figs. 8–11. Figure 1c plots the temperature-dependent resistivity of device #1 at zero magnetic field. Several kinks at around 300 K represent the various CDW transitions of the sample26,27, and the superconductivity transition at around 2 K is shown in the inset. Figure 1d is the polar plot of the normalized in-plane resistance Rab/RN of device #1 from below Tc to above Tc with rotating in-plane magnetic field B// = 5 T, where α ranges from 0° to 360° is the magnetic field angle with respect to the x axis. Two-fold (C2) nematicity in the superconducting state is evident, with minima of Rab/RN appearing at α = 90° and 270° and maxima at α = 0° and 180°. Angular-dependent in-plane upper critical field Hc2//(α) (Supplementary Fig. 7e) and angular-dependent critical current Ic(α) (Supplementary Fig. 7f) of device #1 also shows the same nematic behavior of Rab/RN. Furthermore, this anisotropy is pronounced in the superconducting state (T<Tconset~2K), while it becomes undetectable at B// ≤ 5 T immediately above Tconset (Fig. 1c). The above evidence unambiguously indicates that the two-fold anisotropy arises from nematic superconducting state rather than from resistance anisotropy9.

Figure 1e shows the Hc2///Hp versus reduced temperature T/Tc for α = 0° and 90° in 6R-TaS2 devices #1-#3, where Hp is the Pauli limit. In all the devices, Hc2// exceeds the Pauli limit for both B// directions, demonstrating nematic Ising superconductivity in the materials. The in-plane anisotropy obtained from Hc2,max///Hc2,min// ranges from ~2.6 to 6.5, which is much larger than previous results of Ising-like superconductors13–15. Additionally, nematic behavior is observed throughout the superconducting state of 6R-TaS2, in contrast to the isotropic behavior observed in the ultra-low T regime of another vdWHs 4Hb-TaS26. We have ruled out trivial origins of the C2 behaviors, including current induced vortex motion40, strain stress and misalignment of magnetic field (details in Supplementary Figs. 13, 14), leaving the only possibility of nematic Ising superconductivity. Together with evidence shown later in the text, we argue that this state emerges with the hidden order formed at the normal state of the material.

Nematicity and hidden magnetism below T* ~ 30 K

Above the superconducting temperature and below the CDW transition temperature, a hidden order is found in 6R-TaS2 below T* ~ 30 K, which is characterized by the concurrent emergence of nematicity and hidden magnetism. Figure 2a shows the derivative resistance (dRab(T)/dT) of device #1 at zero magnetic field, exhibiting an abrupt change at around 30 K. Figure 2b presents the longitudinal MR of device #1 under 14 T of out-of-plane magnetic field B⊥, together with the two-fold anisotropic MR (AMR = ΔR/Rmin(α)) under 14 T of in-plane magnetic field B//. Although MR and AMR are from completely different measurement configurations, we found that they could be easily scaled together, and both exhibit an upturn at around 30 K. Figure 2c plots the Hall resistivity of device #1 at various temperatures from 3 K to 80 K, with nonlinear Hall curves apparent at low temperatures. Since such nonlinear Hall curves mainly exhibits two sections with different slops (Hall resistance RH), we plot RH of the low and high magnetic field regions in Fig. 2d, where a bifurcation of the two RH is found below T* ~ 30 K. Figure 2e shows the anomalous Hall conductivity σAHE versus temperature, which again exhibits an upturn below T* ~ 30 K. Figure 2f depicts σAHE insensitive to the angles θ between the magnetic field and the c-axis (up to 75° for T = 8 K, and > 45° for T = 21 K and 25.5 K), proving that the observed signal is indeed AHE20,21. Similar data from devices #2, #3 and #bulk can be found in Supplementary Figs. 18–20.Fig. 2 Hidden order with intertwined physics of nematicity, anomalous Hall effect (AHE), and magnetic field-tunable thermal hysteresis around T* ~ 30 K for device #1.

a T-dependent derivative in-plane resistance (dRab/dT) at zero magnetic field, showing an abrupt change at around T* ~ 30 K. Inset shows measurement configuration, with the S and D represents source and drain electrode, respectively. b T-dependent longitudinal magnetoresistance (MR) under B⊥ = 14 T (extracted from Supplementary Fig. 15) and anisotropic magnetoresistance (AMR) under B// = 14 T (extracted from Supplementary Fig. 16), exhibiting a concurrent upturn at T* ~ 30 K. Here B⊥ represents the vertical magnetic field. c Hall resistivity ρxy(B⊥) at various T from 3 K to 80 K, with nonlinear Hall curves apparent below T* ~ 30 K. Inset shows measurement configuration. d T-dependent Hall coefficientRH≡dρxy/dB extracted from both low and high magnetic field regions, exhibiting a bifurcation below T* ~ 30 K. e T-dependent anomalous Hall conductivity σAHE (main panel) and magnetic field-dependent σAHE (inset) at various T, exhibiting an upturn in σAHE below T* ~ 30 K. The Hall conductivityσxy=−ρxy/(ρxx2+ρxy2), then σAHE obtained by subtracting the local linear ordinary Hall conductivity background. Here ρxx is the longitudinal resistivity. f σAHE(θ)/σAHE(0∘) vs. θ, where θ is the angle between the magnetic field and the c axis under various T (extracted from Supplementary Fig. 17). Inset shows measurement configuration. g–i Magnetic field-tunable thermal hysteresis. g T-dependent resistance with decreasing and increasing T under B⊥ = 0 T and 8 T. Inset is a zoom-in plot of the thermal hysteresis at zero field. h T-dependent ΔR (difference in resistance between warming and cooling ramps), exhibiting the temperature range (between ~11 K to 30 K) of the thermal hysteresis. i B⊥-dependent integral area of the hysteresis loop and B⊥-dependent longitudinal MR at T = 15 K (extracted from Supplementary Fig. 21). The loop area is clearly independent of the MR.

To our knowledge, no AHE has been reported in other non-magnetic TMDs. And the AHE in 6R-TaS2 is not associated with a CDW transition, in contrast to the CDW-induced AHE in non-magnetic kagome metals AV3Sb520,21 and ScV6Sn625. Furthermore, below T* ~ 30 K, a magnetic field-tunable thermal hysteresis loop (tunable by magnetic field up to ~20% and independent of the MR) is observed in the normal state resistance curve (Fig. 2g–i and Supplementary Figs. 21 and 22), suggesting a first-order phase transition which is not associated with conventional structural transitions, but rather, with an electronic symmetry breaking process arising from orbital41–43 or spin44,45 degrees of freedom. Indeed, there has been no report of structural or conventional magnetic phase transitions in 6R-TaS2 ~30 K26,27,46–48. The above results unveil a hidden order with intertwined nematicity and hidden magnetism below T* ~ 30 K in 6R-TaS2. The coexistence of the hidden order with Ising superconductivity produces strongly nematic Ising superconductivity with simultaneous Kondo screening in the 1T layer.

Giant extrinsic AHE and Kondo resonance in 6R-TaS2

AHE reflects the anomalous transverse velocity of charge carriers via extrinsic (side-jump or skew-scattering) and intrinsic mechanisms16. Figure 3a plots σAHE versus σxx of device #1, #2, and #3 together with data from other materials in the literature, spanning over the side-jump (σAHE∝σxx1.6), intrinsic (σAHE∝ constant), and skew-scattering (σAHE∝σxx) regimes16. The linear dependence of σAHE versus σxx is not obvious with a log-log scale in Fig. 3a, and can be clearly seen in Fig. 3b with a linear scale, similar to previous results in MnGe19, gated CsV3Sb549 and doped Fe16,50, pointing to the extrinsic skew-scattering origin16,19,49,50. Note that AHE is rarely found in non-magnetic materials, with limited examples such as kagome metals AV3Sb5 and ScV6Sn6 (chiral-CDW effects)20,21,25, Weyl semimetal ZrTe5 (Berry curvature effects)17, and 6R-TaS2 in this work.Fig. 3 Extrinsic AHE with linear scaling relation and Kondo resonance in non-magnetic transition metal dichalcogenides (TMDs) 6R-TaS2.

a Map of AHE (σAHE vs. σxx) for various materials in double logarithmic coordinates, spanning over the side-jump (σAHE∝σxx1.6), intrinsic (σAHE∝ constant), and skew-scattering (σAHE∝σxx) regimes. The solid lines are the linear fitting of σAHE vs. σxx for 6R-TaS2 device #1, #2 and #3. The dashed line is the linear fitting of σAHE vs. σxx for MnGe19. b The linear dependence of σAHE vs. σxx for 6R-TaS2 samples with a linear scale. c T-dependent dI/dV spectrum of the 1H layer at zero magnetic field (Vs = −5 mV, I = 400 pA). d Zero-field scanning tunneling spectroscopy (STS) of the superconducting gap at T = 0.3 K and atomically resolved topography of the 1H layer (inset). The red dot shows the location where the STS spectra were measured. e Kondo resonance in the 1 T layer from 0.3 K (superconducting state) to up to 12 K (normal state) (Vs = −50 mV, I = 400 pA). f Kondo resonance of the 1 T layer at T = 0.3 K. Inset shows scanning tunneling microscopy (STM) topography of the 1 T layer with the 13×13 charge-density-wave (CDW) pattern. The red dot shows the location where the STS spectra were measured.

In order to elucidate the origin of AHE, we performed scanning tunneling spectroscopy (STS) on 6R-TaS2 bulk crystals with 1H and 1 T termination at temperatures ranging from 0.3 K to 12 K. The STS spectra on the 1H termination (Fig. 3c, d) reveal a normal metal to superconductor phase transition at ~2.5 K with a V-shaped nodal-like pairing gap Δ ~ 0.53 meV, proving 6R-TaS2 as a strong-coupling superconductor (2Δ/kBTc ~ 4.7), in contrast to 4Hb-TaS2 as a weak-coupling superconductor (2Δ/kBTc ~ 3.2)11,35,37. More interestingly, no clear superconducting gap is found in the 1 T termination (Fig. 3e, f), while Kondo screening emerges on each of the CDW centers below ~12 K (normal state) and persists down to 0.3 K (deep inside the superconducting state). The Kondo resonance in 6R-TaS2 is due to the screening of the localized spins on the CDW sites by the metallic 1H layer underneath5, possibly constituting a Kondo lattice like the artificial 1H/1 T heterojunctions5, which is further supported by the spatial dI/dV spectroscopic mapping of Kondo resonance at 4.2 K and 0.3 K (Supplementary Figs. 26 and 27). The quantitative analysis of Kondo resonance and Kondo coupling JK can be seen in Supplementary Fig. 25.

Based on all the findings presented above, we propose that spin-dependent coupling between the 1 T and 1H layers in 6R-TaS2 induces the unconventional hidden magnetism in 6R-TaS2 below T*; the spin state of the QSL in the 1 T layer may be significantly altered, i.e., the highly fluctuating spins may acquire additional broken symmetries akin to the chiral spin liquid phase4,18,51,52. It serves as magnetic scattering impurities without overall net magnetization, giving rise to the large extrinsic AHE in the normal states and hidden magnetism persisting down to the superconducting state in 6R-TaS2.

Phase diagram of nematicity, hidden magnetism and Ising SC

A global phase diagram of the nematicity, hidden magnetism and unconventional Ising superconductivity in 6R-TaS2 crystals is constructed with our experimental findings (Fig. 4a). The phase diagram can be divided into three regions: the isotropic region (gray color), the nematic normal state region with hidden magnetism (blue color), and the nematic Ising superconductivity region with Kondo screening (red color). The first region starts from high temperature down to 30 K, where the crystal is effectively featureless after the CDW transitions at ~305 K, with isotropic in-plane MR (anisotropy <0.5%) down to 30 K. The second region is the normal state below T* ~ 30 K, where nematicity and hidden magnetism emerge simultaneously. Here nematicity is characterized by in-plane AMR as defined in Fig. 2b; hidden magnetism is characterized by 1) large extrinsic AHE under B⊥, 2) magnetic field-tunable thermal hysteresis and 3) Kondo screening as shown in Figs. 3 and 4. The third region is the nematic Ising superconductivity coexisting with Kondo screening at the lower left corner of the phase diagram, defined as the resistivity of the sample below 90% of normal state resistivity RN (purple dots in Fig. 4a), under temperature-dependent critical in-plane fields along α = 90°. Since the orientation of the nematicity is the same for the normal state and the superconducting state, we can conclude that nematicity of the hidden order coexist with Ising superconductivity and substantially modified the superconducting gap in different in-plane directions, producing up to ~10,000% of anisotropy in resistivity. More intriguingly, when a small out-of-plane magnetic field is applied to suppress the superconductivity, the AHE effect is immediately visible without any intermediate state (Supplementary Fig. 28), together with the Kondo screening detected by STS below the superconductivity temperature, proving the coexistence of hidden magnetism and nematic Ising superconductivity in the material.Fig. 4 Phase diagram of nematicity, hidden magnetism and Ising superconductivity in 6R-TaS2.

a B// - T phase diagram of 6R-TaS2 (device #1), divided into three regions: the isotropic region (gray color), the weakly nematic with hidden magnetism region (blue color), and the nematic Ising superconductivity with hidden magnetism region (red color). The purple dots are the boundary between nematic Ising superconducting state and weakly nematic normal state, representing the upper critical field along α = 90° determined by the 90% RN criterion (extracted from Supplementary Fig. 6g). The yellow ribbon around T* is the boundary between the weakly nematic and isotropic normal state regions. The color scale represents nematicity defined as [Rmax(α)−Rmin(α)]/Rmin(α)(%) and extracted from the AMR (See Supplementary Fig. 16). b–d The representative in-plane AMR curves R(α) at 5 K and 10 T (blue color region), 1.6 K and 2 T (red color region), 40 K and 10 T (gray color region), respectively. e In-plane AMR curves at 5 K measured with B// from 1 T to 14 T.

Some key points are worth noting from the phase diagram (Fig. 4a). First, the hidden phase in 6R-TaS2 (T* ~ 30 K), driven by the particular coupling between the 1 T and 1H layers in the 6 R phase, is different from the spin liquid phase proposed in bulk 1T-TaS2 which persists up to 200 K32–34; Second, 6R-TaS2 is the first material in the TaS2 family to exhibit an AHE, which interestingly echoes with theoretical prediction of AHE in triangular lattice systems with Kondo coupling of itinerant electrons and a chiral spin texture53; Third, both the nematic normal state and nematic Ising superconductivity has never been reported in the TaS2 family, which are absent in both the single layer 1H-TaS228 and in ultra-low temperature limit of superconducting 4Hb-TaS26. Fourth, the coexistence of superconductivity and Kondo resonance in 6R-TaS2 has not been observed in either the artificial 1H/1 T heterojunctions5 (without SC gap on its 1H layer) or the natural vdWHs 4Hb-TaS2 (without Kondo resonance on its 1 T layer)36. We argue that although Kondo coupling may occur in all cases (Supplementary Fig. 1), subtle differences in the coupling can lead to large variation in the spin states of the 1 T layer as well as the electronic nematicity in the 1H layer, leading to marked differences in the QSL phases and other exotic physics. Thus manipulating interlayer stacking/coupling of vdWHs may help to construct customizable quantum systems with promising properties (Supplementary Fig. 1)54. Fifth, we note that the origin of the hidden magnetism in 4Hb-TaS2 is still under debate, including theoretical mechanisms related to CSL52, vison-vortex nucleation with ℤ2 topological order55, and type-II heavy Fermi liquids56. For 6R-TaS2, we cannot rule out these three possible mechanisms for the origin of hidden magnetism, as the complex interactions between the 1H and 1 T layers may lead to various possibilities for the spin state of 1 T. On the other hand, nematicity57 has been observed both in the normal state of cuprates58, iron-based59, kagome superconductors60, as well as in the superconducting state of magic-angle graphene9, doped Bi2Se361 and NbSe214. Although the microscopic origin is still unclear, nematicity may coexist, cooperate, or compete with other orders9,57,62. The discovery of a nematic electronic state and hidden magnetism in the normal phase 6R-TaS2 that cooperate with the Ising superconductivity points to the unconventional nature of this hidden order in the material.

Discussion

In conclusion, a hidden order is found in 6R-TaS2 below T* ~ 30 K that is characterized by simultaneous emergence of nematicity and hidden magnetism. The development of hidden magnetism may be related to significantly altered QSL state which acquire additional broken symmetries. Magneto-transport and scanning tunneling spectroscopy data strongly suggest a coexistence of hidden magnetism and nematic Ising superconductivity in 6R-TaS2. The entangled physics of unconventional nematic Ising SC, strong-coupling SC, hidden magnetism, Kondo screening, and 13×13 CDW in 6R-TaS2 may provide a fertile ground for the exploration of pair density waves63,64, chiral superconductivity11, Yu-Shiba-Rusinov (YSR)–like bound states7, spin triplet superconductivity65 and more. This work unveils the potential of natural van der Waals heterostructures as a promising platform for exploring the intertwined and exotic physics.

Methods

Single crystals growth

Single crystal samples of 6R-TaS2 have been prepared by phase transition of 1T-TaS2 to 6 R at 800 °C in an inert atmosphere26,27. The pure 1T-TaS2 single crystals were firstly synthesized by the chemical vapor transport (CVT) method using the raw materials of tantalum powder (99.9%,) and sulfur powder (99.999%), combined in a molar ratio of 1:2. Meanwhile, the transport agent iodine spheres (99.99%) weighing 0.25 g/cm3 were placed into a crucible. Then, the crucible was sealed in a quartz tube under a high vacuum (pressure below 10−5 Torr). This assembly was then placed in a two-zone furnace with 1274 K for the end containing the mixture and 1224 K for the growth end. After about 200 h, the quartz tube was rapidly quenched into ice water. Then 1T-TaS2 single crystals can be obtained. At last, the 1T-TaS2 single crystals were further heated at 800°C in evacuated quartz tubes to obtain the 6R-TaS2 single crystal samples.

Single crystals characterizations

The sample has been characterized by diffraction (XRD) studies in θ–2θ geometry using lab-based sources. A Quantum Design Magnetic Property Measurement System (MPMS-3) was used to measure the temperature- and field-dependent magnetization of the samples. The thickness of the various samples was measured using Atomic Force Microscopy (AFM).

Device fabrication

Al2O3-assisted exfoliation technique was used to obtain thin flakes of 6R-TaS2 crystals with a thickness of down to 7.9 nm (more details in Supplementary Fig. 3). Standard e-beam lithography was used to pattern electrodes, followed by e-beam evaporation of Ti (5 nm) and Au (100 nm). The device fabrication process was carried out in an inert atmosphere and vacuum to minimize sample oxidation, and samples were briefly exposed to air only under PMMA capping layer protection. More details of device fabrication process can be found in Supplementary Fig. 4.

Transport measurements

Transport measurements were conducted at temperatures between 0.3 K and 360 K with magnetic fields up to 14 T using an Oxford Teslatron cryostat and a Quantum Design PPMS. Lock-in amplifiers were used to measure longitudinal resistance (Rxx) and Hall resistance (Rxy) at a frequency of 17.77 Hz with an AC current of 10 μA for nano-devices and 2 mA for the bulk devices. Changing the magnetic field direction was achieved by rotating the sample holder. To eliminate the influence of the slight Hall signals on the raw data of angular dependence of resistivity, the resistivity taken at every angle has been averaged with positive and negative magnetic fields. The in-plane resistivity was measured by the standard four-electrode method. Additionally, the c-axis resistance for the bulk sample was measured by the four-electrode method with the Corbino-shape-like configuration. Various resistivity/resistance quantities have been used in this experiment, and their definitions are summarized in the following: Rab represents the resistance in the ab plane, Rc is the out-of-plane resistance, Rxx and Rxy represent longitudinal resistance and transverse Hall resistance, respectively; while ρxx and ρxy represent longitudinal resistivity and transverse Hall resistivity, respectively.

STM/STS measurements

Bulk 6R-TaS2 samples were cleaved in an ultrahigh vacuum chamber at room temperature, and then immediately inserted into the STM head for further topography and spectrum measurements in continuous ultrahigh vacuum. STM measurements were performed using PtIr tips that were well-calibrated on Au(111) surface. All the dI/dV spectrum were taken using the standard lock-in technique (frequency f = 973 Hz) with a small AC modulation (Vmod = 0.05 mV) added to the DC bias.

Supplementary information

Supplementary Information

Peer Review File

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-51631-z.

Acknowledgements

This project has been supported by the National Key R&D Program of China (Grant No. 2019YFA0308402 (J.C.), 2019YFA0307800 (J.M.)), the Innovation Program for Quantum Science and Technology (Grant No. 2021ZD0302403 (X.C.X. and J.C.)), the National Natural Science Foundation of China (Grant Nos. 11934001 (J.C.), 92265106 (J.C.), 11921005 (J.C.), 11974347 (J.M.), 12074377 (Y.J.)). Shanghai Rising-Star Program (Grant Nos. 23QA1410700 (Y.F.)). Fundamental Research Funds for the Central Universities (J.M. and Y.J.). J.-H.C. acknowledges technical support from Peking Nanofab.

Author contributions

J.-H.C. and S.L. conceived the idea and designed the experiments; S.L. performed all the transport measurements. S.L. and Y.S. performed the AFM measurements; C.T. performed XRD and SQUID measurements with the help of S.J. and S.G.; Y.F., F.H., H.R. and C.C. provided high-quality crystals; S.L. and H.C. fabricated the devices with the help of X.W. and C.T.; C.T., J.C. and J.L. aided the transport measurements; Y.J., J.H.M. and L.C. performed the STM/STS measurements; X.C.X., W.H. and K.T.L. provided theoretical analysis; S.L. and J.-H.C. analyzed the data and wrote the manuscript; M.C. and D.C. aided the modification of the manuscript; all authors commented and modified the manuscript.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

Data availability

Source data are provided in this paper. Data for figures that support the current study are available at 10.7910/DVN/IUPOMC. All raw data generated during the current study are available from the corresponding authors upon request.

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: Shao-Bo Liu, Congkuan Tian, Yuqiang Fang, Hongtao Rong.
==== Refs
References

1. Manzeli S Ovchinnikov D Pasquier D Yazyev OV Kis A 2D transition metal dichalcogenides Nat. Rev. Mater. 2017 2 17033 10.1038/natrevmats.2017.33
Manzeli, S., Ovchinnikov, D., Pasquier, D., Yazyev, O. V. & Kis, A. 2D transition metal dichalcogenides. Nat. Rev. Mater. 2, 17033 (2017).10.1038/natrevmats.2017.33
2. Novoselov K Mishchenko A Carvalho A Neto AC 2D materials andvan der Waals heterostructures Science 2016 353 aac9439 10.1126/science.aac9439 27471306
Novoselov, K., Mishchenko, A., Carvalho, A. & Neto, A. C. 2D materials andvan der Waals heterostructures. Science 353, aac9439 (2016).27471306 10.1126/science.aac9439
3. Lu JM Evidence for two-dimensional Ising superconductivity in gated MoS2 Science 2015 350 1353 1357 10.1126/science.aab2277 26563134
Lu, J. M. et al. Evidence for two-dimensional Ising superconductivity in gated MoS2. Science 350, 1353–1357 (2015).26563134 10.1126/science.aab2277
4. Persky E Magnetic memory and spontaneous vortices in a van der Waals superconductor Nature 2022 607 692 696 10.1038/s41586-022-04855-2 35896649
Persky, E. et al. Magnetic memory and spontaneous vortices in a van der Waals superconductor. Nature 607, 692–696 (2022).35896649 10.1038/s41586-022-04855-2
5. Vaňo V Artificial heavy fermions in a van der Waals heterostructure Nature 2021 599 582 586 10.1038/s41586-021-04021-0 34819682
Vaňo, V. et al. Artificial heavy fermions in a van der Waals heterostructure. Nature 599, 582–586 (2021).34819682 10.1038/s41586-021-04021-0
6. Silber I Two-component nematic superconductivity in 4Hb-TaS2 Nat. Commun. 2024 15 824 10.1038/s41467-024-45169-3 38280890
Silber, I. Two-component nematic superconductivity in 4Hb-TaS2. Nat. Commun. 15, 824 (2024).38280890 10.1038/s41467-024-45169-3
7. Liu MK Monolayer 1T-NbSe2 as a 2D-correlated magnetic insulator Sci. Adv. 2021 7 eabi6339 10.1126/sciadv.abi6339 34797708
Liu, M. K. Monolayer 1T-NbSe2 as a 2D-correlated magnetic insulator. Sci. Adv. 7, eabi6339 (2021).34797708 10.1126/sciadv.abi6339
8. 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
9. Cao Y Nematicity and competing orders in superconducting magic-angle graphene Science 2021 372 264 271 10.1126/science.abc2836 33859029
Cao, Y. Nematicity and competing orders in superconducting magic-angle graphene. Science 372, 264–271 (2021).33859029 10.1126/science.abc2836
10. Chen MY Chen XY Yang H Du ZY Wen H-H Superconductivity with twofold symmetry in Bi2Te3/FeTe0.55Se0.45 heterostructures Sci. Adv. 2018 4 eaat1084 10.1126/sciadv.aat1084 29888330
Chen, M. Y., Chen, X. Y., Yang, H., Du, Z. Y. & Wen, H.-H. Superconductivity with twofold symmetry in Bi2Te3/FeTe0.55Se0.45 heterostructures. Sci. Adv. 4, eaat1084 (2018).29888330 10.1126/sciadv.aat1084
11. Ribak A Chiral superconductivity in the alternate stacking compound 4Hb-TaS2 Sci. Adv. 2020 6 eaax9480 10.1126/sciadv.aax9480 32258393
Ribak, A. et al. Chiral superconductivity in the alternate stacking compound 4Hb-TaS2. Sci. Adv. 6, eaax9480 (2020).32258393 10.1126/sciadv.aax9480
12. Wan P Orbital Fulde–Ferrell–Larkin–Ovchinnikov state in an Ising superconductor Nature 2023 619 46 51 10.1038/s41586-023-05967-z 37225992
Wan, P. et al. Orbital Fulde–Ferrell–Larkin–Ovchinnikov state in an Ising superconductor. Nature 619, 46–51 (2023).37225992 10.1038/s41586-023-05967-z
13. Cui J Transport evidence of asymmetric spin–orbit coupling in few-layer superconducting 1Td-MoTe2 Nat. Commun. 2019 10 2044 10.1038/s41467-019-09995-0 31053717
Cui, J. et al. Transport evidence of asymmetric spin–orbit coupling in few-layer superconducting 1Td-MoTe2. Nat. Commun. 10, 2044 (2019).31053717 10.1038/s41467-019-09995-0
14. Hamill A Two-fold symmetric superconductivity in few-layer NbSe2 Nat. Phys. 2021 17 949 954 10.1038/s41567-021-01219-x
Hamill, A. et al. Two-fold symmetric superconductivity in few-layer NbSe2. Nat. Phys. 17, 949–954 (2021).10.1038/s41567-021-01219-x
15. Zhang E Spin–orbit–parity coupled superconductivity in atomically thin 2M-WS2 Nat. Phys. 2023 19 106 113 10.1038/s41567-022-01812-8
Zhang, E. et al. Spin–orbit–parity coupled superconductivity in atomically thin 2M-WS2. Nat. Phys. 19, 106–113 (2023).10.1038/s41567-022-01812-8
16. Nagaosa N Sinova J Onoda S MacDonald AH Ong NP Anomalous Hall effect Rev. Mod. Phys. 2010 82 1539 1592 10.1103/RevModPhys.82.1539
Nagaosa, N., Sinova, J., Onoda, S., MacDonald, A. H. & Ong, N. P. Anomalous Hall effect. Rev. Mod. Phys. 82, 1539–1592 (2010).10.1103/RevModPhys.82.1539
17. Liang T Anomalous Hall effect in ZrTe5 Nat. Phys. 2018 14 451 455 10.1038/s41567-018-0078-z
Liang, T. et al. Anomalous Hall effect in ZrTe5. Nat. Phys. 14, 451–455 (2018).10.1038/s41567-018-0078-z
18. Machida Y Nakatsuji S Onoda S Tayama T Sakakibara T Time-reversal symmetry breaking and spontaneous Hall effect without magnetic dipole order Nature 2010 463 210 213 10.1038/nature08680 20010605
Machida, Y., Nakatsuji, S., Onoda, S., Tayama, T. & Sakakibara, T. Time-reversal symmetry breaking and spontaneous Hall effect without magnetic dipole order. Nature 463, 210–213 (2010).20010605 10.1038/nature08680
19. Fujishiro Y Giant anomalous Hall effect from spin-chirality scattering in a chiral magnet Nat. Commun. 2021 12 317 10.1038/s41467-020-20384-w 33436576
Fujishiro, Y. et al. Giant anomalous Hall effect from spin-chirality scattering in a chiral magnet. Nat. Commun. 12, 317 (2021).33436576 10.1038/s41467-020-20384-w
20. Yang S-YGiant unconventional anomalous Hall effect in the metallic frustrated magnet candidate, KV3Sb5 Sci. Adv. 2020 6 eabb6003 10.1126/sciadv.abb6003 32789181
Yang, S.-Y. Giant unconventional anomalous Hall effect in the metallic frustrated magnet candidate, KV3Sb5. Sci. Adv. 6, eabb6003 (2020).32789181 10.1126/sciadv.abb6003
21. Yu FH Concurrence of anomalous Hall effect and charge density wave in a superconducting topological kagome metal Phys. Rev. B 2021 104 L041103 10.1103/PhysRevB.104.L041103
Yu, F. H. et al. Concurrence of anomalous Hall effect and charge density wave in a superconducting topological kagome metal. Phys. Rev. B 104, L041103 (2021).10.1103/PhysRevB.104.L041103
22. Mielke C Time-reversal symmetry-breaking charge order in a kagome superconductor Nature 2022 602 245 250 10.1038/s41586-021-04327-z 35140387
Mielke, C. et al. Time-reversal symmetry-breaking charge order in a kagome superconductor. Nature 602, 245–250 (2022).35140387 10.1038/s41586-021-04327-z
23. Guo C Switchable chiral transport in charge-ordered kagome metal CsV3Sb5 Nature 2022 611 461 466 10.1038/s41586-022-05127-9 36224393
Guo, C. et al. Switchable chiral transport in charge-ordered kagome metal CsV3Sb5. Nature 611, 461–466 (2022).36224393 10.1038/s41586-022-05127-9
24. Guguchia Z Hidden magnetism uncovered in a charge ordered bilayer kagome material ScV6Sn6 Nat. Commun. 2023 14 7796 10.1038/s41467-023-43503-9 38016982
Guguchia, Z. Hidden magnetism uncovered in a charge ordered bilayer kagome material ScV6Sn6. Nat. Commun. 14, 7796 (2023).38016982 10.1038/s41467-023-43503-9
25. Yi C Quantum oscillations revealing topological band in kagome metal ScV6Sn6 Phys. Rev. B 2024 109 035124 10.1103/PhysRevB.109.035124
Yi, C. Quantum oscillations revealing topological band in kagome metal ScV6Sn6. Phys. Rev. B 109, 035124 (2024).10.1103/PhysRevB.109.035124
26. Achari A Alternating superconducting and charge density wave monolayers within bulk 6R-TaS2 Nano Lett 2022 22 6268 6275 10.1021/acs.nanolett.2c01851 35857927
Achari, A. et al. Alternating superconducting and charge density wave monolayers within bulk 6R-TaS2. Nano Lett 22, 6268–6275 (2022).35857927 10.1021/acs.nanolett.2c01851
27. Pal S Charge density wave and superconductivity in 6R-TaS2 Phys. B 2023 669 415266 10.1016/j.physb.2023.415266
Pal, S. Charge density wave and superconductivity in 6R-TaS2. Phys. B 669, 415266 (2023).10.1016/j.physb.2023.415266
28. de la Barrera SC Tuning Ising superconductivity with layer and spin–orbit coupling in two-dimensional transition-metal dichalcogenides Nat. Commun. 2018 9 1427 10.1038/s41467-018-03888-4 29650994
de la Barrera, S. C. et al. Tuning Ising superconductivity with layer and spin–orbit coupling in two-dimensional transition-metal dichalcogenides. Nat. Commun. 9, 1427 (2018).29650994 10.1038/s41467-018-03888-4
29. Kvashnin Y Coexistence of superconductivity and charge density waves in tantalum disulfide: experiment and theory Phys. Rev. Lett. 2020 125 186401 10.1103/PhysRevLett.125.186401 33196259
Kvashnin, Y. et al. Coexistence of superconductivity and charge density waves in tantalum disulfide: experiment and theory. Phys. Rev. Lett. 125, 186401 (2020).33196259 10.1103/PhysRevLett.125.186401
30. Li L Superconducting order from disorder in 2H-TaSe2-xSx npj Quantum Mater. 2017 2 1 7 10.1038/s41535-017-0016-9
Li, L. et al. Superconducting order from disorder in 2H-TaSe2-xSx. npj Quantum Mater. 2, 1–7 (2017).10.1038/s41535-017-0016-9
31. Wan, Z. Unconventional superconductivity in chiral molecule-TaS2 hybrid superlattices. Nature 632, 69–74 (2024).
32. Law KT Lee PA 1T-TaS2 as a quantum spin liquid Proc. Natl. Acad. Sci. USA 2017 114 6996 7000 10.1073/pnas.1706769114 28634296
Law, K. T. & Lee, P. A. 1T-TaS2 as a quantum spin liquid. Proc. Natl. Acad. Sci. USA 114, 6996–7000 (2017).28634296 10.1073/pnas.1706769114
33. He W-Y Xu XY Chen G Law KT Lee PA Spinon Fermi surface in a cluster mott insulator model on a triangular lattice and possible application to 1T−TaS2 Phys. Rev. Lett. 2018 121 046401 10.1103/PhysRevLett.121.046401 30095934
He, W.-Y., Xu, X. Y., Chen, G., Law, K. T. & Lee, P. A. Spinon Fermi surface in a cluster mott insulator model on a triangular lattice and possible application to 1T−TaS2. Phys. Rev. Lett. 121, 046401 (2018).30095934 10.1103/PhysRevLett.121.046401
34. Klanjšek M A high-temperature quantum spin liquid with polaron spins Nat. Phys. 2017 13 1130 1134 10.1038/nphys4212
Klanjšek, M. et al. A high-temperature quantum spin liquid with polaron spins. Nat. Phys. 13, 1130–1134 (2017).10.1038/nphys4212
35. Nayak AK Evidence of topological boundary modes with topological nodal-point superconductivity Nat. Phys. 2021 17 1413 1419 10.1038/s41567-021-01376-z
Nayak, A. K. et al. Evidence of topological boundary modes with topological nodal-point superconductivity. Nat. Phys. 17, 1413–1419 (2021).10.1038/s41567-021-01376-z
36. Wen C Roles of the narrow electronic band near the fermi level in 1T-TaS2-related layered materials Phys. Rev. Lett. 2021 126 256402 10.1103/PhysRevLett.126.256402 34241511
Wen, C. et al. Roles of the narrow electronic band near the fermi level in 1T-TaS2-related layered materials. Phys. Rev. Lett. 126, 256402 (2021).34241511 10.1103/PhysRevLett.126.256402
37. Shen, S. et al. Coexistence of quasi-two-dimensional superconductivity and tunable Kondo Lattice in a van der Waals Superconductor. Chin. Phys. Lett. 39, 077401 (2022).
38. Nayak AK First-order quantum phase transition in the hybrid metal–Mott insulator transition metal dichalcogenide 4Hb-TaS2 Proc. Natl. Acad. Sci. USA 2023 120 e2304274120 10.1073/pnas.2304274120 37856542
Nayak, A. K. First-order quantum phase transition in the hybrid metal–Mott insulator transition metal dichalcogenide 4Hb-TaS2. Proc. Natl. Acad. Sci. USA 120, e2304274120 (2023).37856542 10.1073/pnas.2304274120
39. Song X-Y Vishwanath A Zhang Y-H Doping the chiral spin liquid: topological superconductor or chiral metal Phys. Rev. B 2021 103 165138 10.1103/PhysRevB.103.165138
Song, X.-Y., Vishwanath, A. & Zhang, Y.-H. Doping the chiral spin liquid: topological superconductor or chiral metal. Phys. Rev. B 103, 165138 (2021).10.1103/PhysRevB.103.165138
40. Wang Y-L Parallel magnetic field suppresses dissipation in superconducting nanostrips Proc. Natl Acad. Sci. USA 2017 114 E10274 E10280 10.1073/pnas.1619550114 29133405
Wang, Y.-L. et al. Parallel magnetic field suppresses dissipation in superconducting nanostrips. Proc. Natl Acad. Sci. USA 114, E10274–E10280 (2017).29133405 10.1073/pnas.1619550114
41. Liang YG Tuning the hysteresis of a metal-insulator transition via lattice compatibility Nat. Commun. 2020 11 3539 10.1038/s41467-020-17351-w 32669544
Liang, Y. G. et al. Tuning the hysteresis of a metal-insulator transition via lattice compatibility. Nat. Commun. 11, 3539 (2020).32669544 10.1038/s41467-020-17351-w
42. Haverkort MW Orbital-assisted metal-insulator transition in VO2 Phys. Rev. Lett. 2005 95 196404 10.1103/PhysRevLett.95.196404 16384003
Haverkort, M. W. et al. Orbital-assisted metal-insulator transition in VO2. Phys. Rev. Lett. 95, 196404 (2005).16384003 10.1103/PhysRevLett.95.196404
43. Wei Z Two superconductor-insulator phase transitions in the spinel oxide Li1 ± xTi2O4−δ induced by ionic liquid gating Phys. Rev. B 2021 103 L140501 10.1103/PhysRevB.103.L140501
Wei, Z. et al. Two superconductor-insulator phase transitions in the spinel oxide Li1 ± xTi2O4−δ induced by ionic liquid gating. Phys. Rev. B 103, L140501 (2021).10.1103/PhysRevB.103.L140501
44. Liuwan Zhang CI Amlan Biswas RLG Lozanne AD Direct observation of percolation in a manganite thin film Science 2002 298 805 10.1126/science.1077346 12242450
Liuwan Zhang, C. I., Amlan Biswas, R. L. G. & Lozanne, A. D. Direct observation of percolation in a manganite thin film. Science 298, 805 (2002).12242450 10.1126/science.1077346
45. Catalano S Tailoring the electronic transitions of NdNiO3 films through (111)pc oriented interfaces APL Mater. 2015 3 062506 10.1063/1.4919803
Catalano, S. et al. Tailoring the electronic transitions of NdNiO3 films through (111)pc oriented interfaces. APL Mater. 3, 062506 (2015).10.1063/1.4919803
46. Thompson AH The synthesis and properties of 6R-TaS2 Solid State Commun 1975 17 1115 1117 10.1016/0038-1098(75)90267-7
Thompson, A. H. The synthesis and properties of 6R-TaS2. Solid State Commun 17, 1115–1117 (1975).10.1016/0038-1098(75)90267-7
47. Figueroa E Physical properties of 6R-TaS2 J. Solid State Chem. 1992 114 486 490 10.1006/jssc.1995.1073
Figueroa, E. Physical properties of 6R-TaS2. J. Solid State Chem. 114, 486–490 (1992).10.1006/jssc.1995.1073
48. Liu, Y. et al. Tuning the charge density wave and superconductivity in 6R-TaS2−xSex. J. Appl. Phys. 117, 163912 (2015).
49. Zheng G Electrically controlled superconductor-to-failed insulator transition and giant anomalous Hall effect in kagome metal CsV3Sb5 nanoflakes Nat. Commun. 2023 14 678 10.1038/s41467-023-36208-6 36755031
Zheng, G. et al. Electrically controlled superconductor-to-failed insulator transition and giant anomalous Hall effect in kagome metal CsV3Sb5 nanoflakes. Nat. Commun. 14, 678 (2023).36755031 10.1038/s41467-023-36208-6
50. Shiomi Y Onose Y Tokura Y Extrinsic anomalous Hall effect in charge and heat transport in pure iron, Fe0.997Si0.003, and Fe0.97Co0.03 Phys. Rev. B 2009 79 100404 10.1103/PhysRevB.79.100404
Shiomi, Y., Onose, Y. & Tokura, Y. Extrinsic anomalous Hall effect in charge and heat transport in pure iron, Fe0.997Si0.003, and Fe0.97Co0.03. Phys. Rev. B 79, 100404 (2009).10.1103/PhysRevB.79.100404
51. Wen XG Wilczek F Zee A Chiral spin states and superconductivity Phys. Rev. B 1989 39 11413 11423 10.1103/PhysRevB.39.11413
Wen, X. G., Wilczek, F. & Zee, A. Chiral spin states and superconductivity. Phys. Rev. B 39, 11413–11423 (1989).10.1103/PhysRevB.39.11413
52. Lin, S.-Z. Kondo enabled transmutation between spinons and superconducting vortices: Origin of magnetic memory in 4Hb-TaS2. Phys. Rev. Res. 6, 023224 (2024).
53. Martin, I. & Batista, C. D. Itinerant electron-driven chiral magnetic ordering and spontaneous quantum hall effect in triangular lattice models. Phys. Rev. Lett. 101, 156402 (2008).
54. Andrei EY The marvels of moiré materials Nat. Rev. Mater. 2021 6 201 206 10.1038/s41578-021-00284-1
Andrei, E. Y. et al. The marvels of moiré materials. Nat. Rev. Mater. 6, 201–206 (2021).10.1038/s41578-021-00284-1
55. Luo, R. L. & Chen, G. Is spontaneous vortex generation in superconducting 4Hb-TaS2 from vison-vortex nucleation with Z2 topological order? Quantum Frontiers 3, 10 (2024).
56. K¨onig EJ Type-II heavy Fermi liquids and the magnetic memory of 4Hb-TaS2 Phys. Rev. Res. 2024 6 L012058 10.1103/PhysRevResearch.6.L012058
K¨onig, E. J. Type-II heavy Fermi liquids and the magnetic memory of 4Hb-TaS2. Phys. Rev. Res. 6, L012058 (2024).10.1103/PhysRevResearch.6.L012058
57. Fradkin E Kivelson SA Lawler MJ Eisenstein JP Mackenzie AP Nematic fermi fluids in condensed matter physics Annu. Rev. Condens. Matter Phys. 2010 1 153 178 10.1146/annurev-conmatphys-070909-103925
Fradkin, E., Kivelson, S. A., Lawler, M. J., Eisenstein, J. P. & Mackenzie, A. P. Nematic fermi fluids in condensed matter physics. Annu. Rev. Condens. Matter Phys. 1, 153–178 (2010).10.1146/annurev-conmatphys-070909-103925
58. Tranquada JM Sternlieb BJ Axe JD Nakamura Y Uchida S Evidence for stripe correlations of spins and holes in copper oxide Nature 1995 375 561 563 10.1038/375561a0
Tranquada, J. M., Sternlieb, B. J., Axe, J. D., Nakamura, Y. & Uchida, S. Evidence for stripe correlations of spins and holes in copper oxide. Nature 375, 561–563 (1995).10.1038/375561a0
59. Chu J-H In-plane resistivity anisotropy in an underdoped iron arsenide Science 2010 329 824 826 10.1126/science.1190482 20705856
Chu, J.-H. In-plane resistivity anisotropy in an underdoped iron arsenide. Science 329, 824–826 (2010).20705856 10.1126/science.1190482
60. Xiang Y Twofold symmetry of c-axis resistivity in topological kagome superconductor CsV3Sb5 with in-plane rotating magnetic field Nat. Commun. 2021 12 6727 10.1038/s41467-021-27084-z 34795303
Xiang, Y. et al. Twofold symmetry of c-axis resistivity in topological kagome superconductor CsV3Sb5 with in-plane rotating magnetic field. Nat. Commun. 12, 6727 (2021).34795303 10.1038/s41467-021-27084-z
61. Matano K Kriener M Segawa K Ando Y Zheng G-Q Spin-rotation symmetry breaking in the superconducting state of CuxBi2Se3 Nat. Phys. 2016 12 852 854 10.1038/nphys3781
Matano, K., Kriener, M., Segawa, K., Ando, Y. & Zheng, G.-Q. Spin-rotation symmetry breaking in the superconducting state of CuxBi2Se3. Nat. Phys. 12, 852–854 (2016).10.1038/nphys3781
62. Fernandes RM Chubukov AV Schmalian J What drives nematic order in iron-based superconductors? Nat. Phys. 2014 10 97 104 10.1038/nphys2877
Fernandes, R. M., Chubukov, A. V. & Schmalian, J. What drives nematic order in iron-based superconductors? Nat. Phys. 10, 97–104 (2014).10.1038/nphys2877
63. Chen H Roton pair density wave in a strong-coupling kagome superconductor Nature 2021 599 222 228 10.1038/s41586-021-03983-5 34587621
Chen, H. et al. Roton pair density wave in a strong-coupling kagome superconductor. Nature 599, 222–228 (2021).34587621 10.1038/s41586-021-03983-5
64. Zhao H Smectic pair-density-wave order in EuRbFe4As4 Nature 2023 618 940 945 10.1038/s41586-023-06103-7 37380689
Zhao, H. et al. Smectic pair-density-wave order in EuRbFe4As4. Nature 618, 940–945 (2023).37380689 10.1038/s41586-023-06103-7
65. Gu Q Detection of a pair density wave state in UTe2 Nature 2023 618 921 927 10.1038/s41586-023-05919-7 37380691
Gu, Q. et al. Detection of a pair density wave state in UTe2. Nature 618, 921–927 (2023).37380691 10.1038/s41586-023-05919-7
