Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

Unidirectional Kondo scattering in layered NbS2

Abstract

Crystalline defects can modify quantum interactions in solids, causing unintuitive, even favourable, properties such as quantum Hall effect or superconducting vortex pinning. Here we present another example of this notion—an unexpected unidirectional Kondo scattering in single crystals of 2H-NbS2. This manifests as a pronounced low-temperature enhancement in the out-of-plane resistivity and thermopower below 40 K, hidden for the in-plane charge transport. The anomaly can be suppressed by the c-axis-oriented magnetic field, but is unaffected by field applied along the planes. The magnetic moments originate from layers of 1T-NbS2, which inevitably form during the growth, undergoing a charge-density-wave reconstruction with each superlattice cell (David-star-shaped cluster of Nb atoms) hosting a localised spin. Our results demonstrate the unique and highly anisotropic response of a spontaneously formed Kondo-lattice heterostructure, intercalated in a layered conductor.

Introduction

Layered van der Waals materials, such as transition metal dichalcogenides (TMDs), have attracted major interest, thanks to their rich variety of ground states and the possibility of their exfoliation down to an atomically thin level, which remarkably modifies their electronic properties1,2. Recent observations of intriguing physics in artificially assembled heterostructures highlight the importance of interlayer interactions. Examples include the outstanding stability of interlayer excitons in semiconducting TMDs3, and strongly correlated states in twisted bilayer systems4. Relevant aspects of the inter-plane coupling can be deduced by probing out-of-plane charge transport, even in bulk materials5. However, enforcing the current flow strictly along the c axis can be rather challenging due to the crystals’ common flake-like appearance and their propensity for delamination. Such a pitfall can distort the measurement results by orders of magnitude, as demonstrated in our recent study of microstructured samples of 1T-TaS2 with a well-defined current flow6. This observation motivates a careful re-examination of the out-of-plane charge transport properties in this class of materials by adopting the latest state-of-the-art for quantum matter microfabrication7.

Here we present data on the out-of-plane electrical resistivity of bulk monocrystalline 2H-NbS2. This material is one of the three known structural variants of layered NbS2. The two other polytypes are 3R and 1T, the latter occurring only in atomically thin form8,9. As illustrated in Fig. 1a, 1T-NbS2 consists of corner-sharing octahedral NbS6 cells. Layers of 2H- and 3R-NbS2 (Fig. 1b, c) both contain NbS6 units of trigonal prismatic geometry, but exhibit different stacking configurations. The 1T polytype has been attracting interest recently as a candidate for realising a two-dimensional magnetic system10,11. 2H-NbS2 has been actively featured in the literature due to a superconductivity below 6 K, proposed to have a multiband character. It also does not show any charge-density-wave (CDW) order, which is uncommon for metallic TMDs12,13,14,15,16,17. Another distinguishing feature of 2H-NbS2 is its non-trivial synthesis procedure. This polytype is thermodynamically stable in a relatively narrow range of temperatures and reactant stoichiometries18,19,20. Crystals formed during high-temperature growth must be rapidly quenched in order to capture 2H-NbS2 in a metastable room-temperature state. However, X-ray diffraction studies have shown that the resultant material has up to 18% of pairs of neighbouring layers stacked in a 3R-like manner21,22. Additionally, diffuse X-ray scattering experiments22 revealed weak traces of the \(\sqrt {13} \times \sqrt {13}\) CDW reconstruction, which appears as a triangular superlattice of David-star-shaped clusters defined by 13 Nb atoms23. Such a reconstruction is not expected for pure 2H-NbS2 or 3R-NbS2. Earlier theoretical investigations have predicted 1T-NbS2 to be particularly prone to developing such a CDW order10,11. One can therefore conclude that single crystals of 2H-NbS2 contain rare, atomically thin inclusions of the 1T polytype.

Fig. 1: Polytypes of NbS2.
figure 1

Crystalline lattices of 1T-NbS2 (a), 2H-NbS2 (b) and 3R-NbS2 (c). The corresponding 1-, 2- and 3-layer unit cells are marked with black wireframes. The compass shows the directions of the principal axes for all three structures. 2H-NbS2 and 3R-NbS2 share the same structure in the ab plane, but have different stacking of layers along the c axis. For 3R-NbS2, only one Nb atom with 6 nearest S atoms are shown for the middle two layers, for a clearer illustration of the stacking. Images produced with VESTA57.

Our study of 2H-NbS2 revealed a remarkably strong low-temperature anomaly in the compound’s out-of-plane resistivity (ρc), manifesting as a minimum at around 40 K, followed by a pronounced upturn upon further cooling. The feature is simultaneously invisible in the in-plane resistivity (ρab), and shows a highly anisotropic response to magnetic field. Neither 2H-NbSe2 nor 3R-NbS2 display such an anomaly, implying that the phenomenon is linked to the structural defects specific to 2H-NbS2. 1T-NbS2, layers of which are one of such defects, were predicted to form a lattice of unpaired localised spins located at the centre of each David-star CDW superlattice cluster10,11. We argue that planes of magnetic moments, hosted by the inclusions of 1T-NbS2, cause a Kondo effect observable only when the current flows across these planes.

Results

Resistivity anisotropy

Optimisation of sample geometry with focused ion beam (FIB) micro-milling greatly improves charge transport study precision6,7. Using this approach, we shaped single crystals into samples with well-defined, few micron-thick and -wide current channels, oriented along the two principal directions: normal and parallel to the atomic planes (Fig. 2a shows a sample of 2H-NbS2 produced this way). Such a design allowed simultaneous measurement of both ρab and ρc via the four-point technique. Probing ρc on two segments of different surface-to-volume resulted in mutually consistent values, allowing us to ensure that our results were not distorted by the presence of surface-related effects or macroscopic defects.

Fig. 2: Interlayer charge dynamics in 2H-NbS2.
figure 2

a Scanning electron microscope image of a 2H-NbS2 sample, structured with focused ion beam for accurate resistivity anisotropy measurements. False colouring is used: purple – crystal, beige – gold film. The scale bar in the bottom left is 20 µm long. The current sourcing (I) and voltage probing (V) electrodes are labelled. b Plots of the in-plane (ρab, blue) and out-of-plane (ρc, red) resistivities of 2H-NbS2 against temperature (T). The dashed red line stands for the out-of-plane resistivity of 2H-NbSe2, scaled by a factor of 2. Resistivity anisotropy of 2H-NbS2 is plotted in the inset. c Seebeck coefficients of 2H-NbS2 for the out-of-plane (Sc) and in-plane (Sab) directions as functions of temperature, measured on bulk single crystals (note that Sab is negative). The setup for measuring Sc is illustrated schematically. The crystal was approximately 1 mm long in the c axis direction (indicated in the drawing), and 2–3 mm long laterally. The value of Sc is the ratio of the voltage across the sample (Vs) and the thermal gradient across it, determined from the differential thermocouple voltage (VDTC). The sample sat between two copper plates, which homogenised temperature at its two faces and was electrically decoupled from the heatsink by a thin sapphire plate.

Figure 2b shows the plots of ρab and ρc of 2H-NbS2 against temperature (T), as well as their ratio in the inset. Note that in contrast to the earlier study which reports an anisotropy of the order of 1000 (ref. 24), our measured value was as low as 10 at room temperature, monotonically increasing to 180 on cooling. As it was shown for the case of 1T-TaS2 using finite element simulations6, such an overestimate by the previous study could be a result of a non-optimised measurement geometry and an incorrect prior assumption that the anisotropy is very large. While ρab has a conventional metallic temperature dependence, ρc is also metallic, but shows a few noteworthy features. First, the residual out-of-plane resistivity is very high, presumably due to a significant concentration of static defects. Second, ρc approaches saturation in the high-temperature region. This flattening of resistivity may be attributed to the mean free paths decreasing to the point of becoming comparable to the interlayer separation, a concept known as the Mott-Ioffe-Regel limit25. Third, at low temperatures, ρc displays a minimum at around 40 K, with a major upturn at lower temperatures. No corresponding feature exists in ρab (in agreement with the previous results26). All studied samples of 2H-NbS2 showed qualitatively identical behaviour, with slight differences in the absolute values of resistivity—related to slight impurity content variations—and temperatures of the minimum distributed in the 30–40 K range.

We compared ρc of 2H-NbS2 to that of the isostructural and isovalent compound 2H-NbSe2 (dashed line in Fig. 2b). The latter material did not exhibit a similar low-temperature anomaly. Based on the nominal lattice parameters, density functional theory calculations predict that the two compounds will have nearly identical electronic band structures (see Supplementary Note 1 and Supplementary Figs. 1 and 2). We therefore conclude that the upturn of ρc of 2H-NbS2 is not intrinsic to the nominal structure of the compound, but is caused by crystalline lattice defects.

Seebeck coefficient anisotropy

The Seebeck coefficient (S) is a useful quantity for sensitively probing energy landscape variations near Fermi level. The open-circuit voltage generated by the thermal gradient is unaffected by the presence of static defects such as vacancies or non-magnetic stacking faults. On the other hand, the Seebeck coefficient is a function of the energy dependence of the conduction electron scattering rate. This is then strongly affected by the occurrence of resonance peaks in the density of states close to the chemical potential. Seebeck coefficient of 2H-NbS2 revealed a prominent peak at approximately 15 K, appearing only for the out-of-plane thermal gradient (Fig. 2c). The phonon drag phenomenon produces a similar feature in the temperature dependence. However, it manifests only in conductors with long phonon and charge carrier mean free paths—such as semimetals or extremely pure metals—where momentum-conserving scattering is dominant27. This is highly unlikely for 2H-NbS2, as its high content of static defects should clearly favour momentum-relaxing scattering. Furthermore, absence of the corresponding Sab peak rules out the phonon drag from the possible Sc anomaly origins. An alternative interpretation of the peak, the Kondo effect, will be discussed further below.

High-field magnetotransport

The out-of-plane resistivity anomaly of 2H-NbS2 demonstrated a particularly curious response to magnetic fields. As can be seen in Fig. 3a, the transverse and longitudinal out-of-plane magnetoresistances of the material are strikingly different. Applying the field along the c axis suppresses the resistivity upturn, shifting the minimum to lower temperature. Yet even at 63 T the anomaly is still present. Consequently, at 50 K and below, ρc decreases when magnetic field is increased, with signs of saturation appearing around 50 T (Fig. 3b). In contrast, in transverse magnetic field, ρc behaves as a more typical orbital magnetoresistance, common to metals. It is positive, and about three times weaker in magnitude than the longitudinal one (Fig. 3c) without significantly affecting the shape of the upturn in ρc(T).

Fig. 3: Out-of-plane magnetotransport in 2H-NbS2.
figure 3

a Out-of-plane resistivity (ρc) as a function of temperature (T) for various longitudinal (main plot) and transverse (inset) magnetic fields (B). Solid lines and markers represent the measured data. Dashed lines in the main plot are fits according to the numerical renormalisation group theory of Kondo effect. The fitted model includes a field-independent contribution due to the residual temperature-independent resistance as well as the electron–phonon scattering, described by the Bloch-Grüneisen formula (dash-dot line). Dashed lines in the inset are guides for the eye. The two datasets were collected using different samples, which explains the slight resistivity mismatch for zero field. b, c Relative changes in ρc, under out-of-plane (b) and in-plane (c) magnetic fields.

In order to emphasise the observed phenomenon’s highly anisotropic nature, we also report the in-plane magnetotransport of 2H-NbS2 up to 14 T. The magnetoresistance is weak for all field directions, as can be seen in Fig. 4a, but similarly to ρc, ρab is also reduced by the field along the c axis (clearly depicted in Fig. 4b), which could be a trace of the same anomaly. The in-plane field-dependence of ρab (Fig. 4c) is likely governed by orbital effects, like in the case of ρc.

Fig. 4: In-plane magnetotransport in 2H-NbS2.
figure 4

a In-plane resistivity (ρab) as a function of temperature (T) at zero and 14 T magnetic field (B). The change of ρab between zero field and 14 T (out-of-plane field) is plotted in the inset on a logarithmic T scale. The black line in the inset corresponds to the logarithmic temperature dependence. b, c In-plane magnetoresistance for the out-of-plane (b) and in-plane (c) field.

Discussion

The question of the origin of the anomalies of ρc in 2H-NbS2 will now be addressed. A number of phenomena could result in a finite-temperature resistivity minimum in a metal. Resistivity upturns can be caused by electron–electron interactions in the presence of static disorder28,29,30. However, the corresponding quantum mechanical correction is either weakly enhanced by a magnetic field, or is effectively field-independent. A closely linked phenomenon of weak localisation (WL) is also known to produce an additional contribution to resistivity at low temperatures31. In this scenario, when a series of scattering events cause an electron to follow a closed path, quantum interference favours the net backward scattering over the forward one. Magnetic flux threading these scattering loops shifts the phases of the wave functions, diminishing the effect. In our case, when electrons are scattered between different planes, the closed paths should have comparable projections along the in-plane and out-of-plane directions. However, the upturn is only influenced by the c-axis-oriented field, contradicting the WL-based interpretation. A metal–insulator transition or conduction based on a thermally activated hopping between defects32 would cause a divergence of ρc at the lowest temperatures, which was not the case. The possibility of quantum tunnelling playing a significant role is ruled out based on a linear relation between current and voltage (Supplementary Note 2 and Supplementary Fig. 3). Resistivity upturns have also been observed in strongly doped cuprate superconductors33,34. In those materials, the effect is believed to be caused by scattering from magnetic droplets forming around non-magnetic impurities. This interpretation, however, relies on the existence of strong electronic correlations, and therefore does not apply to our system.

Finite-temperature resistivity minimum in a metal is also a well-known signature of the Kondo effect, a scattering of conduction electrons of dilute localised magnetic moments35. Besides the upturn, the characteristic features of the phenomenon, observable in charge transport, include a negative curvature of ρc(T) at the lowest temperatures and a suppression of the upturn by magnetic field, which causes spin-flip scattering to become inelastic36,37. The observed peak in the Seebeck coefficient is also characteristic to dilute and concentrated spin systems, including Kondo lattices38,39. It originates from the resonant scattering in the Kondo channel at the Fermi level. Above the Kondo temperature, TK, the resonance is smeared out, and depending on the specifics of a system, the peak in S appears at a temperature between 0.3 and 0.9 TK (refs. 39,40). When temperature is low enough, the localised spins are screened and the excitations obey simple power laws, like those of a Fermi liquid. For example, S varies as T/TK for T/TK < 0.1–0.15 for several typical Kondo alloys in the dilute, single impurity limit41.

Measurements of ρc under high pressure, presented in Supplementary Note 3 and Supplementary Fig. 4, show that the upturn remains extremely robust up to the highest achieved pressure of 1.9 GPa. Applying the pressure weakly shifts the minimum of ρc up in temperature. This is consistent with the behaviour expected from Kondo systems42,43.

We therefore argue that scattering of magnetic impurities is the most fitting explanation of our observations. The temperature dependence of ρc in 2H-NbS2 is consistent with the one expected from the numerical renormalisation group (NRG) theory calculations for Kondo effect36,44, as illustrated by the fit in Fig. 3a. We modelled ρc with a sum of three contributions: a temperature-independent residual resistivity ρ0, an electron–phonon scattering term ρe–p (captured by the Bloch-Grüneisen formula) and the Kondo term ρK (the only magnetic-field-dependent term), for which we used the common empirical expression closely following the results of the NRG theory36,44:

$$\rho _{\mathrm{K}}\left( T \right) = \rho _{{{{\mathrm{K0}}}}}\left( {1 + \left( {2^{1/\alpha } - 1} \right)\left( {T/T_{\mathrm{K}}} \right)^2} \right)^{ - \alpha }.$$
(1)

The fitting procedure is described in more detail in Supplementary Note 4, with the help of Supplementary Figs. 5 and 6 and Supplementary Table 1.

This explanation immediately raises a question regarding the nature of our system’s magnetic impurities. The standard scenario where magnetic atoms are uniformly distributed clearly does not fit our picture. Doping 2H-NbS2 with Fe results in the upturn observable in ρab as well as the disappearance of superconductivity45. Additionally, the undoped material does not display a corresponding signature in the heat capacity46. Lack of pronounced anomaly effects on ρab implies that the responsible defects take the form of sparse planes, extending along the layers. When the current then flows along the layers, only a small fraction of the conduction electrons move in close proximity to these planes. But for the out-of-plane current flow, effectively all charge carriers have to pass through them, resulting in a particularly strong influence. Although we observed planar irregularities in the crystalline lattice via transmission electron microscopy, their atomic structure could not be determined due to a limited resolution (see Supplementary Note 5 and Supplementary Figs. 7 and 8). Taking a closer look at the in-plane magnetoresistance for the c-axis-oriented field reveals that the difference between ρab at 0 and 14 T increases linearly with respect to ln(T) between 30 and 15 K (Fig. 4a inset), which further supports our hypothesis. This means that a minute contribution of Kondo scattering is present in ρab, but it is not strong enough to change the sign of the gradient of ρab(T).

The presence of 1T-NbS2 layers, evidenced by the characteristic CDW signatures22, offers a fascinating interpretation of our findings, illustrated in Fig. 5. As we mentioned in the introduction, the \(\sqrt {13} \times \sqrt {13}\) CDW order, associated with the 1T polytype, forms a triangular superlattice of David-star-shaped clusters defined by 13 Nb atoms23. The electronic structure of monolayer 1T-NbS2 as well as 1T-NbSe2 in such a configuration has been predicted to contain one very flat band around the Fermi level. This makes the materials susceptible to electronic instabilities like Mott localisation, with a concomitant magnetic order10,47,48. The referenced works found the ferromagnetic insulating state as the most stable, although others have proposed that such triangular lattices can host antiferromagnetic spin-liquid phases49,50. These magnetic planes play the role of scatterers in the Kondo effect, however, their concentration appears to be too low for a detection via magnetometry measurements. The described scenario is conceptually similar to the Kondo effect occurring in artificially fabricated magnetic tunnel junctions51,52, yet in our case the phenomenon is observed in a spontaneously formed system. The same kind of Kondo interaction has been very recently observed in scanning tunnelling spectroscopy studies of the 2H/1T or 1H/1T heterostructures of NbSe2 (ref. 53), TaS2 (ref. 54) and TaSe2 (ref. 55), grown by molecular beam epitaxy. One outstanding question is the anomaly’s markedly different response to the two orientations of magnetic field. This difference is probably coming from the localised electron’s highly anisotropic g-factor, causing a very small spin splitting (less than TK), but could also be related to the magnetic ordering. Sizeable anisotropy of the g-factor is expected for systems with strong spin–orbit coupling, such as TMDs50.

Fig. 5: 1T-NbS2 inclusions in 2H-NbS2.
figure 5

Schematic visualisation of the proposed interpretation of the observed out-of-plane charge transport anomaly. The 2H-NbS2 crystal (structure in the centre) contains inclusion layers of 1T-NbS2 (highlighted in red). 1T-NbS2 undergoes a \(\sqrt {13} \times \sqrt {13}\) charge-density-wave (CDW) reconstruction. The Nb atoms in the CDW state are arranged into David-star-shaped clusters, superlattice of which is depicted on the right (only Nb atoms are shown, the dashed line marks the unit cell after the reconstruction). Each cluster contains an unpaired localised spin at the centre (red arrows on the right). The orientations of spins in the illustration are arbitrary and are not meant to suggest any particular ordering. This array of localised magnetic moments causes the itinerant electrons in the 2H-NbS2 bulk to experience Kondo scattering during the out-of-plane current flow (blue arrows).

Since 2H-NbS2 is known to contain frequent 3R-like stacking faults, it is natural to ask whether the anomaly is somehow caused by the inclusions of 3R-NbS2. We measured the latter compound’s out-of-plane resistivity, and while the corresponding temperature dependence was surprisingly found to be non-metallic, the extremely weak reaction of the interlayer conduction to the longitudinal magnetic field (Δρc/ρc ≈ 0.1% at 14 T) was incompatible with the behaviour observed in 2H-NbS2 (see Supplementary Note 6 as well as Supplementary Fig. 9 for the relevant data on 3R-NbS2). The abundance of these stacking faults could explain the high residual component of ρc. The current understanding is that the poor conductivity of 3R-NbS2 is not intrinsic, but rather originates from the disorder due to self-intercalated Nb atoms56. However, in 2H-NbS2, the abundance of stacking faults results in high residual component of ρc and good in-plane metallicity.

In summary, we have demonstrated that a delicate alternation of the interlayer crystalline structure of 2H-NbS2, by introducing different polymorphs of the same atomic composition, dramatically affects the material’s physical properties. In particular, the crystal’s 2H stacking is occasionally disrupted by the 1T layers, which undergo a CDW instability. This then results in a triangular superlattice of David-star-shaped clusters, each hosting a lone spin at the centre. Such a texture of localised magnetic moments can be seen as a two-dimensional Kondo lattice, immersed into the metallic bulk of the 2H polytype. When the electric field or thermal gradient are then applied along the c axis, electronic transport shows a pronounced Kondo effect manifesting as anomalies in the out-of-plane resistivity and Seebeck coefficient. But when they are applied within the plane, there is no sign of spin-dependent scattering. The observation of this highly anisotropic phenomenon occurring in a naturally formed heterostructure has been made possible, thanks to the careful tailoring of the crystal by FIB. Our work therefore shows the importance of adopting new experimental techniques in studying novel electronic materials, especially highly anisotropic Van der Waals structures.

Methods

Focused ion beam microfabrication

Microstructured samples were extracted from monocrystalline flakes of TMDs. The starting crystals had the lateral size of the order of 1 mm and were at least 100 µm thick. After identifying a clean region on a crystal’s surface, free of cracks or buckling, a rectangular lamella was defined by milling away the surrounding material using an FEI Helios G4 Xe plasma FIB microscope. The typical dimensions of a lamella were around 120 µm × 60 µm × 5 µm (with up to 20% variations in lengths between different samples), with the intermediate dimension corresponding to the extent along the c axis of a crystal. The milling current for this stage was 60 nA, with 30 kV column voltage. An FEI Helios G3 Ga FIB microscope was then used for polishing the surface of the lamella with a 1 nA beam in order to ensure the parallelism of the two largest faces. After extraction, the lamella was glued to a sapphire substrate with a tiny amount of Araldite Rapid epoxy, keeping the external face exposed. Besides anchoring the lamella, the epoxy also formed a meniscus around it that smoothly connected the substrate’s surface to the lamella’s exposed face. The setup was then sputter coated with a 100 nm layer of gold. Next, the Ga FIB milling at 10 nA was used for defining the probing electrodes by selectively removing the sputtered gold layer, and for patterning the lamella in order to form the current channel and voltage probing points. The procedure was concluded with polishing the exposed side faces of the sample with a 1 nA ion beam in order to clean the surface of the re-deposited material and define the final dimensions of the device. Since the entire bottom face of the sample was rigidly attached to the substrate, differential thermal contraction and compressibility were expected to produce inhomogeneous stresses throughout the lamella. In our study, these stresses did not have a significant influence on the measured data. More detailed information about the FIB-assisted sample preparation can be found in the relevant review paper7 and references therein.

Resistivity measurements

Resistivity was measured via the four-point technique with direct or alternating excitation currents in the 20–40 µA range. Temperature sweeps’ rate was limited to 1 K/min for the ambient pressure and of 0.5 K/min for the high-pressure measurements in order to reduce the thermal lag and gradients.

Resistivity at high pressure was measured using a piston cylinder cell produced by C&T Factory. Daphne oil 7474 was used as a pressure-transmitting medium. Pressure was determined from the changes in resistance and superconducting transition temperature of a sample of Pb located next to the 2H-NbS2 sample.

Measurements in high magnetic fields were conducted at the high magnetic field facilities in Grenoble (up to 34 T DC field) and Toulouse (up to 63 T pulsed field). Quantum Design PPMS was used for measurements in fields up to 14 T.

Seebeck coefficient measurements

Seebeck coefficient was measured using an in-house setup. For the in-plane Seebeck coefficient measurement, a thin and long sample was mounted on a ceramic bar. One end of the bar was connected to the thermal bath, while the other one had a resistive heater attached. A differential thermocouple was used to measure the temperature difference across the sample. The out-of-plane Seebeck coefficient measurement was performed using a setup displayed in Fig. 2c and described in the corresponding caption.

Data availability

The data that support the findings of this study are available from the authors (K.S., E.M., and L.F.) upon reasonable request. Data used for generating the plots can be found at https://doi.org/10.6084/m9.figshare.16720621.

References

  1. 1.

    Novoselov, K. S. et al. Two-dimensional gas of massless Dirac fermions in graphene. Nature 438, 197–200 (2005).

    CAS  Google Scholar 

  2. 2.

    Ellis, J. K., Lucero, M. J. & Scuseria, G. E. The indirect to direct band gap transition in multilayered MoS2 as predicted by screened hybrid density functional theory. Appl. Phys. Lett. 99, 261908 (2011).

    Google Scholar 

  3. 3.

    Wang, Z. et al. Evidence of high-temperature exciton condensation in two-dimensional atomic double layers. Nature 574, 76–80 (2019).

    CAS  Google Scholar 

  4. 4.

    Cao, Y. et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43–50 (2018).

    CAS  Google Scholar 

  5. 5.

    Basov, D. N. & Timusk, T. Electrodynamics of high-Tc superconductors. Rev. Mod. Phys. 77, 721–779 (2005).

    CAS  Google Scholar 

  6. 6.

    Martino, E. et al. Preferential out-of-plane conduction and quasi-one-dimensional electronic states in layered 1T-TaS2. npj 2D Mater. Appl. 4, 7 (2020).

  7. 7.

    Moll, P. J. W. Focused ion beam microstructuring of quantum matter. Annu. Rev. Condens. Matter Phys. 9, 147–162 (2018).

    Google Scholar 

  8. 8.

    Carmalt, C. J., Manning, T. D., Parkin, I. P., Peters, E. S. & Hector, A. L. Formation of a new (1T) trigonal NbS2 polytype via atmospheric pressure chemical vapour deposition. J. Mater. Chem. 14, 290 (2004).

    CAS  Google Scholar 

  9. 9.

    Shimakawa, M., Kawachi, K., Nishikawa, S. & Hayashi, K. Structural stability of the 1T structure on transition-metal dichalcogenides. J. Solid State Chem. 129, 242–249 (1997).

    CAS  Google Scholar 

  10. 10.

    Tresca, C. & Calandra, M. Charge density wave and spin 1/2 insulating state in single layer 1T-NbS2. 2D Mater. 6, 035041 (2019).

    CAS  Google Scholar 

  11. 11.

    Wang, W. et al. Charge density wave instability and pressure-induced superconductivity in bulk 1T-NbS2. Phys. Rev. B 102, 155115 (2020).

    CAS  Google Scholar 

  12. 12.

    Guillamón, I. et al. Superconducting density of states and vortex cores of 2H-NbS2. Phys. Rev. Lett. 101, 166407 (2008).

    Google Scholar 

  13. 13.

    Yan, R. et al. Thickness dependence of superconductivity in ultrathin NbS2. Appl. Phys. Express 12, 023008 (2019).

    Google Scholar 

  14. 14.

    Leroux, M., Rodière, P., Cario, L. & Klein, T. Anisotropy and temperature dependence of the first critical field in 2H-NbS2. Phys. B Condens. Matter 407, 1813–1815 (2012).

    CAS  Google Scholar 

  15. 15.

    Leroux, M. et al. Anharmonic suppression of charge density waves in 2H-NbS2. Phys. Rev. B 86, 155125 (2012).

    Google Scholar 

  16. 16.

    Lin, D. et al. Patterns and driving forces of dimensionality-dependent charge density waves in 2H-type transition metal dichalcogenides. Nat. Commun. 11, 2406 (2020).

    CAS  Google Scholar 

  17. 17.

    Heil, C. et al. Origin of superconductivity and latent charge density wave in NbS2. Phys. Rev. Lett. 119, 087003 (2017).

    Google Scholar 

  18. 18.

    Jellinek, F., Brauer, G. & Müller, H. Molybdenum and niobium sulphides. Nature 185, 376–377 (1960).

    CAS  Google Scholar 

  19. 19.

    Fisher, W. G. & Sienko, M. J. Stoichiometry, structure, and physical properties of niobium disulfide. Inorg. Chem. 19, 39–43 (1980).

    CAS  Google Scholar 

  20. 20.

    Witteveen, C. et al. Polytypism and superconductivity in the NbS2 system. Dalton Trans. 50, 3216–3223 (2021).

    CAS  Google Scholar 

  21. 21.

    Katzke, H. Stacking disorder in 2H-NbS2 and its intercalation compounds Kx(H2O)yNbS2. I. Description and model calculations of stacking faults in the host lattice NbS2. Z. Kristallogr. Cryst. Mater. 217, 127–130 (2002).

  22. 22.

    Leroux, M., Cario, L., Bosak, A. & Rodière, P. Traces of charge density waves in NbS2. Phys. Rev. B 97, 195140 (2018).

    CAS  Google Scholar 

  23. 23.

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

    CAS  Google Scholar 

  24. 24.

    Pfalzgraf, B. W. & Spreckels, H. The anisotropy of the upper critical field Hc2 and electrical resistivity in 2H-NbS2. J. Phys. C. 20, 4359–4367 (1987).

    CAS  Google Scholar 

  25. 25.

    Gunnarsson, O., Calandra, M. & Han, J. E. Colloquium: saturation of electrical resistivity. Rev. Mod. Phys. 75, 1085–1099 (2003).

    Google Scholar 

  26. 26.

    Naito, M. & Tanaka, S. Electrical transport properties in 2H-NbS2, -NbSe2, -TaS2 and -TaSe2. J. Phys. Soc. Jpn. 51, 219–227 (1982).

    CAS  Google Scholar 

  27. 27.

    Goldsmid, H. J. Introduction to Thermoelectricity (Springer Science & Business Media, 2009).

  28. 28.

    Aronov, B. L. & Al’tshuler, A. G. Contribution to the theory of disordered metals in strongly doped semiconductors. Zh. Eksp. Teor. Fiz. 77, 2028–2044 (1979).

    Google Scholar 

  29. 29.

    Lee, P. A. & Ramakrishnan, T. V. Disordered electronic systems. Rev. Mod. Phys. 57, 287–337 (1985).

    CAS  Google Scholar 

  30. 30.

    Xu, Y., Zhang, J., Cao, G., Jing, C. & Cao, S. Low-temperature resistivity minimum and weak spin disorder of polycrystalline La2∕3Ca1∕3MnO3 in a magnetic field. Phys. Rev. B 73, 224410 (2006).

  31. 31.

    Altshuler, B. L., Khmel’nitzkii, D., Larkin, A. I. & Lee, P. A. Magnetoresistance and Hall effect in a disordered two-dimensional electron gas. Phys. Rev. B 22, 5142–5153 (1980).

    CAS  Google Scholar 

  32. 32.

    Forro, L., Ilakovac, V., Cooper, J. R., Ayache, C. & Henry, J.-Y. Out-of-plane conductivity of YBa2Cu3O7-δ. Phys. Rev. B 46, 6626–6629 (1992).

    CAS  Google Scholar 

  33. 33.

    Chen, W., Andersen, B. M. & Hirschfeld, P. J. Theory of resistivity upturns in metallic cuprates. Phys. Rev. B 80, 134518 (2009).

    Google Scholar 

  34. 34.

    Alloul, H., Bobroff, J., Gabay, M. & Hirschfeld, P. J. Defects in correlated metals and superconductors. Rev. Mod. Phys. 81, 45–108 (2009).

    CAS  Google Scholar 

  35. 35.

    Kondo, J. Resistance minimum in dilute magnetic alloys. Prog. Theor. Phys. 32, 37–49 (1964).

    CAS  Google Scholar 

  36. 36.

    Costi, T. A. et al. Kondo decoherence: finding the right spin model for iron impurities in gold and silver. Phys. Rev. Lett. 102, 056802 (2009).

    CAS  Google Scholar 

  37. 37.

    Hanl, M. et al. Iron impurities in gold and silver: comparison of transport measurements to numerical renormalization group calculations exploiting non-Abelian symmetries. Phys. Rev. B 88, 075146 (2013).

    Google Scholar 

  38. 38.

    Hewson, A. C. The Kondo Problem to Heavy Fermions (Cambridge University Press, 2003).

  39. 39.

    Zlatić, V., Costi, T. A., Hewson, A. C. & Coles, B. R. Thermoelectric power of concentrated Kondo systems. Phys. Rev. B 48, 16152–16155 (1993).

    Google Scholar 

  40. 40.

    Bickers, N. E., Cox, D. L. & Wilkins, J. W. Self-consistent large-N expansion for normal-state properties of dilute magnetic alloys. Phys. Rev. B 36, 2036–2079 (1987).

    CAS  Google Scholar 

  41. 41.

    Cooper, J. R., Vucic, Z. & Babic, E. The thermoelectric power of AlMn alloys. J. Phys. F. Met. Phys. 4, 1489–1500 (1974).

    CAS  Google Scholar 

  42. 42.

    Schilling, J. S. & Holzapfel, W. B. Effect of pressure on the Kondo temperature of Cu:Fe – existence of a universal resistivity curve. Phys. Rev. B 8, 1216–1227 (1973).

    CAS  Google Scholar 

  43. 43.

    Crone, J. & Schilling, J. Effect of pressure on the Kondo temperatures of Au(Fe) and Au(Mn). Solid State Commun. 17, 791–794 (1975).

    CAS  Google Scholar 

  44. 44.

    Goldhaber-Gordon, D. et al. From the Kondo regime to the mixed-valence regime in a single-electron transistor. Phys. Rev. Lett. 81, 5225–5228 (1998).

    CAS  Google Scholar 

  45. 45.

    Nobukane, H., Tabata, Y., Kurosawa, T., Sakabe, D. & Tanda, S. Coexistence of the Kondo effect and spin glass physics in Fe-doped NbS2. J. Phys. Condens. Matter 32, 165803 (2020).

    CAS  Google Scholar 

  46. 46.

    Kačmarčı́k, J. et al. Specific heat measurements of a superconducting NbS2 single crystal in an external magnetic field: energy gap structure. Phys. Rev. B 82, 014518 (2010).

    Google Scholar 

  47. 47.

    Pasquier, D. & Yazyev, O. V. Charge density wave phase, Mottness, and ferromagnetism in monolayer 1T-NbSe2. Phys. Rev. B 98, 045114 (2018).

    CAS  Google Scholar 

  48. 48.

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

    CAS  Google Scholar 

  49. 49.

    Iqbal, Y., Hu, W.-J., Thomale, R., Poilblanc, D. & Becca, F. Spin liquid nature in the Heisenberg J1-J2 triangular antiferromagnet. Phys. Rev. B 93, 144411 (2016).

    Google Scholar 

  50. 50.

    Law, K. T. & Lee, P. A. 1T-TaS2 as a quantum spin liquid. Proc. Natl Acad. Sci. USA 114, 6996–7000 (2017).

    CAS  Google Scholar 

  51. 51.

    Lee, K. I. et al. Kondo effect in magnetic tunnel junctions. Phys. Rev. Lett. 98, 107202 (2007).

    CAS  Google Scholar 

  52. 52.

    Zheng, C., Shull, R. D., Chen, P. J. & Pong, P. W. T. Kondo effect in magnetic tunnel junctions with an AlOx tunnel barrier. Phys. Lett. A 380, 2237–2241 (2016).

    CAS  Google Scholar 

  53. 53.

    Liu, M. et al. Monolayer 1T-NbSe2 as a correlated magnetic insulator. Preprint at https://arxiv.org/abs/2103.12948 (2021).

  54. 54.

    Vaňo, V. et al. Artificial heavy fermions in a van der Waals heterostructure. Preprint at https://arxiv.org/abs/2103.11989 (2021).

  55. 55.

    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).

    CAS  Google Scholar 

  56. 56.

    Youbi, Z. E. et al. Fermiology and electron–phonon coupling in the 2H and 3R polytypes of NbS2. Phys. Rev. B 103, 155105 (2021).

    Google Scholar 

  57. 57.

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

    CAS  Google Scholar 

Download references

Acknowledgements

We would like to express gratitude to Dr. Osor S. Barišić (Institute of Physics in Zagreb), Prof. Andrew Mackenzie (MPI CPfS Dresden), Prof. Fakher Assaad (University of Würzburg), Prof. Frederic Mila (EPFL), Dr. Reza Zamani (EPFL) and particularly Prof. John Cooper (University of Cambridge) for valuable discussions and feedback. We acknowledge the support of Dr. Gaetan Giriat (EPFL) concerning the instrumentation and high-pressure cells, Dr. Maja Bachmann (MPI CPfS Dresden) for her assistance with FIB microfabrication, Dr. Wen Hua (David) Bi and Davor Tolj (EPFL) for their aid with the characterisation of crystals, Dr. Diego Pasquier (EPFL) for auxiliary numerical calculations. We acknowledge the support of the European Magnetic Field Laboratory (EMFL) for access to a 34 T static magnet at LNCMI-CNRS in Grenoble (Proposal: GMA04-217), and access to a 70 T pulsed magnet at LNCMI-CNRS in Toulouse (Proposal: TSC05-119). This study has been funded by the Swiss National Science Foundation through its SINERGIA network MPBH and grants No. 200021_175836 and PP00P2_170544. C. Putzke and P.J.W.M. acknowledge the support by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (grant agreement No 715730) and the Max-Planck-Society. S.N.Z., Q.S.W. and O.V.Y. acknowledge support from NCCR Marvel.

Author information

Affiliations

Authors

Contributions

E.M. and K.S. prepared and conducted resistivity and Seebeck coefficient measurements and wrote the manuscript together with L.F. C.P., M.K. and P.M. assisted with the FIB fabrication process. H.B. synthesised the crystals used in the study. D.L. assisted with resistivity measurements at the high DC magnetic field facility in Grenoble. M.L. and C.P. conducted resistivity measurements at the pulsed magnetic field facility in Toulouse. A.A. secured the magnet time for the experiments in Toulouse. H.K. and C.K. conducted the TEM study. S.Z., Q.W. and O.Y. provided theoretical support and the DFT data. L.F. is the project leader.

Corresponding authors

Correspondence to Edoardo Martino, László Forró or Konstantin Semeniuk.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

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

Supplementary information

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Martino, E., Putzke, C., König, M. et al. Unidirectional Kondo scattering in layered NbS2. npj 2D Mater Appl 5, 86 (2021). https://doi.org/10.1038/s41699-021-00265-6

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1038/s41699-021-00265-6

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing