Abstract
A proximity effect at a van der Waals (vdW) interface enables creation of an emergent quantum electronic ground state. Here we demonstrate that an originally superconducting two-dimensional (2D) NbSe2 forms a ferromagnetic ground state with spontaneous spin polarization at a vdW interface with a 2D ferromagnet V5Se8. We investigated the anomalous Hall effect (AHE) of the NbSe2/V5Se8 magnetic vdW heterostructures, and found that the sign of the AHE was reversed as the number of the V5Se8 layer was thinned down to the monolayer limit. Interestingly, the AHE signal of those samples was enhanced with the in-plane magnetic fields, suggesting an additional contribution to the AHE signal other than magnetization. This unusual behavior is well reproduced by band structure calculations, where the emergence of the Berry curvature along the spin-degenerate nodal lines in 2D NbSe2 by the in-plane magnetization plays a key role, unveiling a unique interplay between magnetism and Zeeman-type spin-orbit interaction in a non-centrosymmetric 2D quantum material.
Similar content being viewed by others
Introduction
Van der Waals (vdW) superstructures constructed with a variety of two-dimensional (2D) materials play a central role in modern condensed matter physics and materials science1,2, providing a perspective that has never been discussed so far. A remarkable example is an emergent strong correlation effect in a moiré flat band system made of weakly correlated materials such as graphene and semiconducting transition-metal dichalcogenides (TMDCs)3,4,5, where the motion of electrons is quenched and the low-energy electronic properties are governed by the Coulombic interaction. Another intriguing example is a strong proximity effect at a vdW interface composed of atomically thin quantum materials, where the hybridization of electron wavefunctions leads to generation of a novel quantum electronic ground state that is missing in individual materials, as exemplified by the emergent topological states observed in the graphene/TMDC heterostructures as well as in the NbSe2/CrBr3 heterostructures6,7,8. Those studies also verify a strong electronic coupling at a vdW interface despite that the constituent materials are weakly bonded via the vdW force, suggesting that a vdW interface should provide an ideal material platform for designing and creating novel physical properties and functionalities.
Among various possibilities, fabrication of a magnetic vdW heterostructure is an interesting and important research direction both for fundamental and applied researches2,9,10. Recent studies on the WSe2/CrI3 heterostructures have revealed that a magnetic exchange interaction is present at a vdW interface even though each layer is well separated by a vdW gap9,10, providing the so-called âvalley-Zeeman effectâ11,12,13 at zero field in WSe2 induced by the exchange field from neighboring CrI3 across the vdW interface. In this system, however, the ground state of WSe2 should be still non-magnetic without spontaneous spin polarization, because WSe2 is a semiconductor characterized with the fully occupied band, where the numbers of the up-spin and down-spin electrons below the Fermi level (EF) should be always equal even when proximitized by a ferromagnet. On the other hand, a magnetic proximity effect in a metallic system characterized with the partially occupied band should lead to generation of a novel magnetic ground state with spontaneous spin polarization by producing an imbalance between the numbers of the up-spin and down-spin electrons below EF. Moreover, such a metallic system should enable us to examine the low-energy electronic properties of the system by the anomalous Hall effect (AHE) measurements, through which we should be able to obtain fundamental information on the electronic structure of the system near EF. Taken together, studying a magnetic proximity effect in a metallic vdW system should be of significant importance, but such an attempt has not been reported so far.
In this study, we demonstrate a proximity-induced ferromagnetic ground state in atomically thin NbSe2 at a vdW interface with a 2D ferromagnet, and uncover an intriguing feature of ferromagnetic NbSe2 associated with its unique band structure by the AHE measurements. NbSe2 is one of representative metallic H-type TMDCs showing the charge-density wave (CDW) and the superconducting (SC) transition at low temperature14, while it is magnetically inactive due to highly delocalized nature of 4d electrons in NbSe2. An unprecedented feature of NbSe2 in the context of 2D materials research arises when thinned down to the monolayer limit, where broken in-plane inversion symmetry relevant to the trigonal prismatic structure (see Fig. 1a) and large spin-orbit interaction (SOI) lift the spin degeneracy near EF. This results in generation of the out-of-plane spin-polarized electrons near EF as schematically illustrated in Fig. 1b15,16,17, while the system is non-magnetic due to time-reversal symmetry. Such SOI is termed as âZeeman-type SOIâ or âIsing-type SOIâ, which plays an essential role in the spin-valley locking effect in monolayer H-type TMDCs, providing intriguing physical phenomena both in semiconducting and metallic H-type TMDCs18,19. On the other hand, an interplay between Zeeman-type SOI and magnetism has not been unveiled yet.
As for a 2D ferromagnet, we employed V5Se8, a layered magnet characterized by the 2D VSe2 sheets separated by the layers of the periodically aligned V atoms (see Fig. 1c). This compound is known to be an itinerant antiferromagnet in the bulk form20, whereas it exhibits weak itinerant ferromagnetism with negative AHE in the thin film form, which survives down to the 2D limit21. We recently succeeded in constructing a magnetic vdW heterostructure with an atomically abrupt vdW interface based on NbSe2 and V5Se8 by molecular-beam epitaxy (MBE), and demonstrated that there is a strong proximity effect working at this NbSe2/V5Se8 (Nb/V) interface characterized by the robust out-of-plane magnetic anisotropy and the large enhancement of the transition temperature22. Those results strongly suggest that there should be an interface exchange interaction present in the Nb/V heterostructures and that the electronic state of NbSe2 should be also modulated, although no clear evidence is provided. Here, we systematically investigated the magneto-transport properties of the Nb/V heterostructures, where the number of the V5Se8 layer (N, varied) was set to be thinner than that of the NbSe2 layer (4âL, fixed) so that the transport properties are dominated by NbSe2. The details of the sample fabrication and characterization are described in our previous study22, and brief summary is shown in Methods and Supplementary Note 1.
Results
The AHE of the NbSe2/V5Se8 van der Waals heterostructures
Figure 2a shows the AHE data taken at Tâ=â2âK for the series of samples, where N was systematically reduced from 6âL to 1.2âL. The Nâ=â6âL sample exhibited negative AHE with clear magnetic hysteresis loop, indicating the robust out-of-plane magnetic anisotropy induced in V5Se8 by a strong proximity effect from neighboring NbSe2 as we discussed in the previous study22. On the other hand, as N was decreased, the AHE signal was at first suppressed at around Nâ=â3.0âL and then developed again with the opposite sign below Nâ<â3âL. Such a sign change has never been observed in the V5Se8 individual films down to the 2D limit21, implying essential contribution from NbSe2.
The evolution of the AHE signals depending on N was more systematically investigated by the temperature-dependence measurements. Figure 2b shows the temperature dependence of the AH resistance at the saturated regime (RAH, sat) for the different N samples. For the thick-enough regime (Nâ=â6âL, for example), negative AHE developed monotonously below the transition temperature (~30âK). This behavior is basically the same as those of the V5Se8 individual films (although the transition temperatures are largely enhanced for the heterostructures as mentioned above)21,22, suggesting that negative AHE observed with thick-enough V5Se8 should be attributed to originally ferromagnetic V5Se8. On the other hand, as N was decreased, a positive component started to develop below 25âK, and the sign of the AHE at the lowest temperature was reversed for Nâ<â3âL. Considering that the 3 L-thick V5Se8 individual film shows rather insulating behavior at low temperature (in particular below 10âK) while the 4 L-thick NbSe2 individual film exhibits metallic behavior even below 2 K21,23, it is natural to consider that the electrical conductions of the Nâ<â3âL samples at the low temperature regime are predominantly governed by NbSe2 (see Supplementary Note 2). Given that the AHE is absent in a non-magnetic material, the obtained results suggest that NbSe2 is in a ferromagnetic state, and positive AHE observed with thin-enough V5Se8 should be attributed to ferromagnetically proximitized NbSe2. The sign reversal of the AHE signal shown in Fig. 2a should be therefore originating from the fact that the dominant layer providing larger contribution to the electrical conduction is varied from V5Se8 to NbSe2 as the number of the V5Se8 layer is reduced.
The angle dependence of the AHE of the NbSe2/V5Se8 heterostructure
The AHE arising from NbSe2 is quite surprising in view of the fact that NbSe2 is a well-known superconductor associated with CDW, where no experimental signatures of magnetism have been discussed so far. To get insights into the origin of this novel AHE in the NbSe2/V5Se8 heterostructure system, we performed further experiments on the angle dependence of the AHE for those samples exhibiting positive AHE. Figure 3a displays the AHE of the Nâ=â2.0âL sample taken at Tâ=â2âK with different field angles (θ). We here focus on the high-enough field regime, where the magnetization direction is fixed to the applied field direction. Given that the AHE signal is proportional to the out-of-plane component of the total magnetization, the signal should be reduced when the field is tilted to the in-plane direction by following cos (θ) as schematically illustrated in Fig. 3b. However, the AHE in the present system did not follow this expected behavior. Figure 3c shows the normalized AHE signal as a function of θ taken at μ0Hâ=â9âT, demonstrating a substantial deviation from cos (θ) as highlighted by the yellow-colored area. This suggests that there is an additional contribution to the AHE signal other than magnetization arising when the in-plane magnetization becomes finite, which has never been observed in other 2D quantum material systems. Remarkably, as will be discussed in the following sections, this intriguing angle dependence as well as the positive sign of the AHE signal could be well reproduced by theoretical calculations based on the band structure of ferromagnetically proximitized NbSe2, suggesting their intrinsic origins rather than the extrinsic origins associated with the scattering mechanisms. Moreover, such a good agreement between experiments and theory strongly suggests that NbSe2 forms a ferromagnetic ground state with spontaneous spin polarization at the NbSe2/V5Se8 vdW interface.
Calculations of the band structure of ferromagnetic NbSe2
Now we discuss the origin of the AHE arising from ferromagnetic NbSe2. We here consider the band structure of monolayer NbSe2 assuming that a magnetic proximity effect is limited to the single layer in contact with V5Se8, but a more realistic case with multilayer NbSe2 is discussed in detail in Supplementary Note 7, which provides qualitatively the same results. Figure 4a depicts the band dispersion of monolayer NbSe2 near EF along the Î-K-M direction in the momentum space (and its time-reversal pair), showing spin splitting due to Zeeman-type SOI. Figure 4b maps the magnitude of this spin splitting energy (Îup-down) in the first Brillouin Zone (BZ), which has three-fold symmetry as expected for the trigonal crystal structure. The magnitude of Îup-down at the K and Kâ valleys is exactly the same but opposite in sign because of time-reversal symmetry. Figure 4c illustrates the Fermi surface (FS) of monolayer NbSe2 and the Berry curvature of the lower band (i.e., the band having lower energy at finite k) in the first BZ, ΩLB (k), both of which are symmetric against time-reversal operation. In such a case, the AH conductivity (Ïxy) calculated by the integration of the Berry curvature below EF over the entire BZ becomes exactly zero, and no AHE expected. Another important aspect is that the spin degeneracy is protected along the Î-M direction due to mirror symmetry (see Fig. 1a) as shown by the dashed lines in Fig. 4b, c, which appears only at the Î and M points in the line cut along the Î-K-M direction (Fig. 4a).
Let us now discuss the band structure of ferromagnetic NbSe2 by considering the situation when NbSe2 is subjected to the out-of-plane exchange field. Figure 4d shows the band structure of monolayer NbSe2 with the exchange field (|M|â=â40âmeV) applied parallel to the c-axis (M//c). Now, the up-spin band and the down-spin band are shifted to the opposite directions due to Zeeman effect, and the magnitude of Îup-down at the K and Kâ valleys becomes different (see Fig. 4d). Accordingly, the FS becomes largely distorted (see Fig. 4f), and the numbers of the up-spin and down-spin electrons below EF become unequal, providing a ferromagnetic ground state with spontaneous spin polarization. In such a ferromagnetic state, we expect non-zero Ïxy due to imperfect cancellation of the Berry curvature below EF between two bands.
Figure 4g shows the resultant Ïxy as a function of energy calculated from the band structure of ferromagnetic NbSe2, showing non-zero Ïxy at Eâ=âEF. Moreover, the sign of Ïxy at Eâ=âEF is positive, which is consistent with the experimentally observed positive AHE for the Nâ<â3âL samples, suggesting that positive AHE is indeed originating from ferromagnetic NbSe2. Considering that the K and Kâ valleys are now energetically inequivalent, this ferromagnetic state would possibly be interpreted as a âferrovalleyâ state with spontaneous valley polarization as well24, which is a natural consequence of the spin-valley locking effect in H-type TMDCs associated with Zeeman-type SOI15,16,17. We note that the spin-degenerate nodal lines shown as the dashed lines in Fig. 4b, c are now moved away from the Î-M direction, forming the closed triangular loops surrounding the K valleys (see Fig. 4e, f). The corners of those loops are contacted with the FS near the Î valley (see Fig. 4f), which is essentially important to explain the results of the angle-dependence measurements as will be discussed in the next section.
Calculations of the angle dependence of the AHE of ferromagnetic NbSe2
Figure 5a presents the angle dependence of Ïxy at Eâ=âEF calculated from the band structure of ferromagnetic NbSe2. Surprisingly, a deviation from cos (θ) that we experimentally observed for the Nâ<â3âL samples (see Fig. 3c) was successfully reproduced. This characteristic deviation from cos (θ) could be considered as the enhancement of the AHE signal with the in-plane fields, which could be understood as a consequence of the emergence of the additional Berry curvature with the in-plane magnetization25, where the intersection of the up-spin band and the down-spin band (i.e., the spin-degenerate nodal lines) associated with Zeeman-type SOI plays an essential role. As we explained above, those lines are originally located along the Î-M direction without the exchange field (see Fig. 4b, c), while they are moved to surround the K valleys under the presence of the out-of-plane exchange field (see Fig. 4e, f). Importantly, the location of those nodal lines in the momentum space is determined by the magnitude of the out-of-plane magnetization (Mz), and Mzâ=â40âmeV satisfies the condition that the nodal lines come close to the corner of the FS near the Î valley (see Fig. 4d, f), providing the largest deviation from cos (θ) as will be discussed later.
The top panel of Fig. 5b shows the magnified view of the band structure near one of such nodal lines, corresponding to the dotted rectangular region in Fig. 4d. There is only one crossing point visible in this plot, but this is in reality distributed in the momentum space to surround the K valleys as mentioned above (also schematically shown in Fig. 5c). When the magnetization direction is tilted to the in-plane direction (for example, θâ=â20° is considered in Fig. 5d), the in-plane component of the magnetization (Mxy) becomes finite, which hybridizes the up-spin band and the down-spin band and opens an energy gap along the nodal lines as shown in Fig. 5d, e. Importantly, this gap opening accompanies the generation of the Berry curvature along the nodal lines as shown in the bottom panel of Fig. 5d by the yellow-colored area. Figure 5f illustrates the distribution of this emergent Berry curvature at θâ=â20° defined as δΩLB (k)â=âΩLB (k, 20°)âââΩLB (k, 0°) in the momentum space, which exactly traces the nodal lines (loops) surrounding the K valleys. Also shown is the FS of monolayer NbSe2 in a ferromagnetic state with the exchange field Mzâ=â40âmeV, which is contacted with the peak of the emergent Berry curvature δΩLB (k) near the Î valley. In this situation, the imbalance between Ïxy arising from the lower band and that from the upper band becomes maximum, providing the largest deviation from cos (θ). We however note that such a deviation from cos (θ) is observable in a rather broad range of the exchange field as long as E is fixed at EF, whereas it disappears when E is increased near to the valence band top (see Supplementary Note 6). This suggests that the observed phenomena associated with the emergence of the Berry curvature with the in-plane magnetization are unique to group-V metallic H-type TMDCs, and that the crossing of the FS and the nodal lines near the Î valley plays a key role for the enhancement of the AHE signal with the in-plane fields.
Discussion
The present study demonstrates a proximity-induced ferromagnetic ground state with spontaneous spin polarization in a metallic 2D quantum material NbSe2 originally showing the CDW/SC states. We emphasize that the unique angle dependence of the AHE characterized with the enhancement of the AHE signal with the in-plane fields observed by experiments could be well reproduced by theoretical calculations based on the band structure of ferromagnetic NbSe2, providing firm evidence that NbSe2 forms a ferromagnetic ground state at the interface with V5Se8. Based on the comparison between the experimental and theoretical results, we estimated that the magnetic exchange interaction at the Nb/V interface should be as large as a few tens of millielectronvolt, which is larger than those reported for the WSe2-based magnetic heterostructures9,10,26 and comparable to the value recently reported for the WS2-based heterostructures27. As for a possible origin of the anomalous enhancement of the AHE signal with the in-plane fields, we interpret from theoretical consideration that the generation of the Berry curvature along the spin-degenerate nodal lines in 2D NbSe2 by the in-plane magnetization should play an essential role, resulting from a unique interplay between magnetism and Zeeman-type SOI in 2D NbSe2. Furthermore, recent demonstration of topological superconductivity in the NbSe2/CrBr3 heterostructures implies the possible existence of a topological superconducting state in our all-epitaxial scalable Nb/V heterostructures as well8, which should provide an ideal platform for future topological quantum computing device applications.
Methods
Sample fabrication and characterization
All the Nb/V heterostructures were grown on commercially available insulating Al2O3 (sapphire) (001) substrates (SHINKOSHA Co., Ltd.) by MBE by following our previously established process21,22,23. A substrate was cleaned by ultrasonication in acetone and ethanol, respectively, annealed in air at 1000âoC for three hours to make a surface atomically flat, and then transferred into the ultrahigh vacuum chamber with a base pressure below ~1âÃâ10â7âPa (EIKO Engineering, Ltd.). Prior to the film growth, a substrate surface was treated with the Se flux at 900âoC for an hour. Then, the substrate temperature was decreased down to 450 oC, and the V5Se8 layer was grown at 450âoC. During the growth, V was supplied by an electron beam evaporator with the evaporation rate of ~0.01âà /s, while Se was supplied by a standard Knudsen cell throughout the growth process with the rate of ~2.0âà /s. After the growth of the V5Se8 layer, the NbSe2 layer was formed at the same growth temperature. Nb was supplied by an electron beam evaporator with the rate of ~0.01âà /s as well. The whole growth process was monitored by reflection high energy electron diffraction (RHEED). Clear RHEED intensity oscillations were observed during the growth, confirming the 2D layer-by-layer growth mode for the formation of the V5Se8 layer and the NbSe2 layer. The exact thickness was designed by counting the number of the RHEED intensity oscillation, and confirmed by the x-ray diffraction (XRD) measurement (PANalytical, XâPert MRD) after the growth. The local structure of the obtained heterostructure was characterized by scanning transmission electron microscope (STEM) measurement (JEOL, JEM-ARM200F). All the samples were covered by insulating Se capping layers to protect their surfaces from oxidization. See Supplementary Note 1 for the details.
Transport measurements
All the samples were cut into Hall-bar shape before transport measurements by mechanical scratching through Se capping layers to define the channel regions, which were typically a few hundred micrometers. The electrical transport properties were characterized by Physical Property Measurement System (Quantum Design, PPMS).
Theoretical calculations
The numerical calculations were performed by using density functional theory (DFT) and a multi-orbital tight-binding model. The electronic structure of monolayer NbSe2 was obtained by using quantum-ESPRESSO28, a package of numerical codes for DFT calculations, with the projector augmented wave method. The cut-off energy was set to 50âRy for the plane wave basis and 500âRy for the charge density. The convergence criterion of 10â8âRy was used in the calculations. The lattice constant was estimated to be 3.475âà by using a lattice relaxation code in quantum-ESPRESSO.
A magnetic proximity effect on monolayer NbSe2 was simulated by an exchange potential in a multi-orbital tight binding model describing electronic states in the pristine monolayer crystal. The tight-binding model was constructed on the basis of eleven Wannier orbitals, five d-orbitals in Nb atom and six p-orbitals in top and bottom Se atoms (see Fig. 1a) with spin degrees of freedom. Here, the hopping integrals were computed by using Wannier90 (ref. 29), a code to provide maximally localized Wannier orbitals and hopping integrals between them from a first-principles band structure. In this model, the exchange potential was introduced as an on-site potential represented by the Zeeman coupling \({H}_{z}=-{{{{{\bf{M}}}}}} \cdot {{{{{\boldsymbol{\sigma }}}}}}/2\). The direction of the magnetization was fixed to the yz plane with the tilting angle \(\theta={{{{{\rm{arctan }}}}}}\left({M}_{y}/{M}_{z}\right)\). Here, the coupling constants (\({M}_{y},\;{M}_{z}\)) are the elements of the exchange potential along the y and z axes, and assumed to be independent of orbital characters for simplicity. The Fermi energy \({E}_{{{{{{\rm{F}}}}}}}\) was numerically estimated under the condition of charge neutrality between valence electrons and nuclei.
The Berry curvature was calculated by using a theoretical expression30,
where |Ïnkã and \({E}_{n}=\hslash {\omega }_{n}\) are the wave function and the eigen-energy, respectively, of the electronic state with the band index \(n\) and the wave vector \({{{{{\bf{k}}}}}}\). The anomalous Hall conductivity (Ïxy) was obtained by the integration of the Berry curvatures below EF in the BZ,
with the FermiâDirac distribution function \({f}_{n}.\) The numerical calculation was performed at zero temperature.
Data availability
The data within this paper are available from the corresponding author upon reasonable request.
References
Geim, A. K. & Grigorieva, I. V. Van der Waals heterostructures. Nature 499, 419â425 (2013).
Gong, C. & Zhang, X. Two-dimensional magnetic crystals and emergent heterostructure devices. Sci. (80-.). 363, eaav4450 (2019).
Cao, Y. et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43â50 (2018).
Cao, Y. et al. Correlated insulator behaviour at half-filling in magic-angle graphene superlattices. Nature 556, 80â84 (2018).
Wang, L. et al. Correlated electronic phases in twisted bilayer transition metal dichalcogenides. Nat. Mater. 19, 861â866 (2020).
Avsar, A. et al. Spin-orbit proximity effect in graphene. Nat. Commun. 5, 4875 (2014).
Island, J. O. et al. Spin-orbit-driven band inversion in bilayer graphene by the van der Waals proximity effect. Nature 571, 85â89 (2019).
Kezilebieke, S. et al. Topological superconductivity in a van der Waals heterostructure. Nature 588, 424â428 (2020).
Zhong, D. et al. Van der Waals engineering of ferromagnetic semiconductor heterostructures for spin and valleytronics. Sci. Adv. 3, e1603113 (2017).
Zhong, D. et al. Layer-resolved magnetic proximity effect in van der Waals heterostructures. Nat. Nanotechnol. 15, 187â191 (2020).
Aivazian, G. et al. Magnetic control of valley pseudospin in monolayer WSe2. Nat. Phys. 11, 148â152 (2015).
Srivastava, A. et al. Valley Zeeman effect in elementary optical excitations of monolayer WSe2. Nat. Phys. 11, 141â147 (2015).
Macneill, D. et al. Breaking of valley degeneracy by magnetic field in monolayer MoSe2. Phys. Rev. Lett. 114, 037401 (2015).
Wilson, J. A., Di Salvo, F. J. & Mahajan, S. Charge-density waves and superlattices in the metallic layered transition metal dichalcogenides. Adv. Phys. 24, 117â201 (1975).
Zhu, Z. Y., Cheng, Y. C. & Schwingenschlögl, U. Giant spin-orbit-induced spin splitting in two-dimensional transition-metal dichalcogenide semiconductors. Phys. Rev. B 84, 153402 (2011).
Xiao, D. et al. Coupled spin and valley physics in monolayers of MoS2 and other group-VI dichalcogenides. Phys. Rev. Lett. 108, 196802 (2012).
Xu, X., Yao, W., Xiao, D. & Heinz, T. F. Spin and pseudospins in layered transition metal dichalcogenides. Nat. Phys. 10, 343â350 (2014).
Xi, X. et al. Ising pairing in superconducting NbSe2 atomic layers. Nat. Phys. 12, 139â144 (2016).
Ye, Y. et al. Electrical generation and control of the valley carriers in a monolayer transition metal dichalcogenide. Nat. Nanotechnol. 11, 598â602 (2016).
Kitaoka, Y. & Yasuoka, H. NMR investigations on the spin fluctuations in itinerant antiferromagnets I. V3Se4 and V5Se8. J. Phys. Soc. Jpn. 48, 1460â1469 (1980).
Nakano, M. et al. Intrinsic 2D ferromagnetism in V5Se8 epitaxial thin films. Nano Lett. 19, 8806â8810 (2019).
Matsuoka, H. et al. Spin-orbit-induced Ising ferromagnetism at a van der Waals interface. Nano Lett. 21, 1807â1814 (2021).
Matsuoka, H. et al. Angle dependence of Hc2 with a crossover between the orbital and paramagnetic limits. Phys. Rev. Res. 2, 012064(R) (2020).
Tong, W., Gong, S., Wan, X. & Duan, C. Concepts of ferrovalley material and anomalous valley Hall effect. Nat. Commun. 7, 13612 (2016).
Habe, T. & Koshino, M. Anomalous Hall effect in 2H-phase MX2 transition-metal dichalcogenide monolayers on ferromagnetic substrates (M =Mo, W, and X = S, Se, Te). Phys. Rev. B 96, 085411 (2017).
Zhao, C. et al. Enhanced valley splitting in monolayer WSe2 due to magnetic exchange field. Nat. Nanotechnol. 12, 757â762 (2017).
Norden, T. et al. Giant valley splitting in monolayer WS2 by magnetic proximity effect. Nat. Commun. 10, 4163 (2019).
Giannozzi, P. et al. QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials. J. Phys. Condens. Matter 21, 395502 (2009).
Pizzi, G. et al. Wannier90 as a community code: new features and applications. J. Phys. Condens. Matter 32, 165902 (2020).
Yao, Y. et al. First principles calculation of anomalous Hall conductivity in ferromagnetic bcc Fe. Phys. Rev. Lett. 92, 037204 (2004).
Acknowledgements
We are grateful to M. S. Bahramy and K. Ishizaka for valuable discussions, and also to Y. Wang, Y. Kashiwabara, Y. Majima, Y. Tanaka, S. Yoshida, B. K. Saika, and M. Kawasaki for experimental help. This work was supported by Grants-in-Aid for Scientific Research (Grant Nos. 19H05602, 19H02593, 19H00653, 20H01840, 20H00127, and 21K13888) and A3 Foresight Program from the Japan Society for the Promotion of Science (JSPS), and by PRESTO (Grant No. JPMJPR20AC) and CREST (Grant No. JPMJCR20T3) from Japan Science and Technology Agency.
Author information
Authors and Affiliations
Contributions
H.M. grew and characterized the samples, performed transport measurements, and analyzed the data. T.H. and M.K. contributed to theoretical interpretation. M.N. and Y.I. supervised the study. H.M., T.H., Y.I., and M.N. wrote the manuscript. All the authors discussed the results and commented on the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Peer review
Peer review information
Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.
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/.
About this article
Cite this article
Matsuoka, H., Habe, T., Iwasa, Y. et al. Spontaneous spin-valley polarization in NbSe2 at a van der Waals interface. Nat Commun 13, 5129 (2022). https://doi.org/10.1038/s41467-022-32810-2
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41467-022-32810-2