Spin-valley locking in the normal state of a transition-metal dichalocogenide superconductor

Spin-valley locking in the normal state of a
transition-metal dichalocogenide superconductor

L. Bawden SUPA, School of Physics and Astronomy, University of St. Andrews, St. Andrews, Fife KY16 9SS, United Kingdom    S. P. Cooil    F. Mazzola Department of Physics, Norwegian University of Science and Technology (NTNU), N-7491 Trondheim, Norway    J. M. Riley SUPA, School of Physics and Astronomy, University of St. Andrews, St. Andrews, Fife KY16 9SS, United Kingdom Diamond Light Source, Harwell Campus, Didcot, OX11 0DE, United Kingdom    L. J. Collins-McIntyre SUPA, School of Physics and Astronomy, University of St. Andrews, St. Andrews, Fife KY16 9SS, United Kingdom    V. Sunko SUPA, School of Physics and Astronomy, University of St. Andrews, St. Andrews, Fife KY16 9SS, United Kingdom Max Planck Institute for Chemical Physics of Solids, Nöthnitzer Straße 40, 01217 Dresden, Germany    K. Hunvik Department of Physics, Norwegian University of Science and Technology (NTNU), N-7491 Trondheim, Norway    M. Leandersson    C. M. Polley    T. Balasubramanian MAX IV Laboratory, Lund University, P. O. Box 118, 221 00 Lund, Sweden    T. K. Kim    M. Hoesch Diamond Light Source, Harwell Campus, Didcot, OX11 0DE, United Kingdom    J. W. Wells Department of Physics, Norwegian University of Science and Technology (NTNU), N-7491 Trondheim, Norway    G. Balakrishnan Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom    M. S. Bahramy Quantum-Phase Electronics Center and Department of Applied Physics, The University of Tokyo, Tokyo 113-8656, Japan RIKEN center for Emergent Matter Science (CEMS), Wako 351-0198, Japan    P. D. C. King To whom correspondence should be addressed: philip.king@st-andrews.ac.uk SUPA, School of Physics and Astronomy, University of St. Andrews, St. Andrews, Fife KY16 9SS, United Kingdom
July 7, 2019

The metallic transition-metal dichalcogenides (TMDCs) are benchmark systems for studying and controlling intertwined electronic orders in solids, with superconductivity developing upon cooling from a charge density wave state. The interplay between such phases is thought to play a critical role in the unconventional superconductivity of cuprates, Fe-based, and heavy-fermion systems, yet even for the more moderately-correlated TMDCs, their nature and origins have proved highly controversial. Here, we study a prototypical example, -NbSe, by spin- and angle-resolved photoemission and first-principles theory. We find that the normal state, from which its hallmark collective phases emerge, is characterised by quasiparticles whose spin is locked to their valley pseudospin. This results from a combination of strong spin-orbit interactions and local inversion symmetry breaking. Non-negligible interlayer coupling further drives a rich three-dimensional momentum-dependence of the underlying Fermi surface spin texture. Together, these findings necessitate a fundamental re-investigation of the nature of charge order and superconducting pairing in NbSe and related TMDCs.

In combination with broken structural inversion symmetry, spin-orbit coupling (SOC) provides a powerful route to stabilise spin-polarised electronic states without magnetism. This can give rise to electrically-tuneable spin splittings via the Rashba effect, promising new technological developments in spintronics, 1 and underpins the formation of spin-helical Dirac cones at the surfaces of topological insulators. 2 It is strongly desired to realise similar effects in systems where more pronounced electronic interactions drive the emergence of collective phases. In non-centrosymmetric superconductors, for example, spin-splitting driven by strong SOC is expected to induce a mixing of spin-triplet and singlet superconducting order parameters, 3 and offers potential for stabilising topological superconductors. 4 Yet, identifying suitable candidate materials has proved a major challenge to date. Partly this is driven by a relative dearth of non-centrosymmetric metals, while for non-magnetic systems in which the centre of inversion is maintained, robust spin-degeneracies of their electronic states would typically be expected due to the dual constraints of time-reversal and inversion symmetry.

In contrast, we show here from spin- and angle-resolved photoemission (ARPES) measurements that the normal state of the centrosymmetric TMDC superconductor -NbSe (hereafter denoted NbSe) hosts a strong layer-resolved and momentum-dependent spin-polarisation of its electronic states at and in the vicinity of the Fermi level. We attribute this as a consequence of a recently-realised form of spin polarisation that can emerge in globally centrosymmetric materials in which constituent structural units nonetheless break inversion symmetry. 5; 6; 7 Together with first-principles calculations, we show how this drives a critical and complex interplay of interlayer interactions and spin-orbit coupling in NbSe. This yields a rich underlying spin-polarised electronic landscape from which charge order and superconductivity emerge upon cooling.

Figure 1: Superconductivity and charge density wave order in -NbSe. (a) Centrosymmetric bulk crystal structure (side view) of -NbSe. This is formed by stacking non-centrosymmetric layers of symmetry with 180 relative rotations, restoring the bulk inversion centre. (b) Resistivity measurements show clear signatures of charge-density wave formation at  K and superconductivity at  K (magnified in the insets). (c) Normal-state Fermi surface measured by ARPES with  eV (left-hand-side) and  eV (right-hand-side);  meV. This consists of two Nb-derived barrels centred around the zone-corner () points, two Nb-derived barrels at the Brillouin zone centre, and an additional central diffuse pocket (most visible at  eV) predominantly derived from Se  orbitals.


Bulk electronic properties of NbSe. We consider here exclusively the layered polymorph (Fig. 1(a)). Each layer forms a graphene-like honeycomb structure with Nb occupying the “A” sublattice and two Se atoms situated on the “B” sublattice. These lie out of the basel plane, equidistant above and below the transition metal. The unit cell contains two such layers, stacked along the -axis with 180 in-plane rotation. Our resistivity measurements (Fig. 1(b)) from single-crystal NbSe samples show a metallic temperature-dependence. They additionally exhibit a pronounced hump at a temperature of  K indicative of charge-density wave (CDW) formation,8 as well as a sharp superconducting transition at  K. The corresponding normal-state Fermi surface is shown in Fig. 1(c), as measured by ARPES. There are two barrels centred around each zone-corner () point which are strongly trigonally-warped. Two further barrels are centred at the point, the inner of which is hexagonal while the outer exhibits additional warping. From our first-principles density-functional theory calculations, and consistent with previous studies, 9; 10; 11 we assign all four of these Fermi surface sheets as being predominantly derived from Nb orbitals. Additional spectral weight at the zone centre is evident in our ARPES measurements for selected photon energies, which we attribute as a fifth, highly three-dimensional, Fermi surface sheet of predominantly Se orbital character.

Figure 2: Spin-polarised bulk electronic structure. (a) Dispersion measured by ARPES ( eV) along the direction. (b)-(f) Spin-resolved EDCs at the momenta marked in (a), revealing a strong spin polarisation of these electronic states along .
Figure 3: Spin-valley locked Fermi surfaces. (a) Dispersion measured ( eV) along the direction, together with the corresponding spin-polarisation of an MDC at the Fermi level and (b-d) EDCs at the momenta marked in (a). These reveal how the sign of the spin polarisation for each zone-corner Fermi surface sheet becomes locked to the valley degree of freedom, as shown schematically in (e).

This sheet also contributes diffuse filled-in intensity, due to the finite out-of-plane momentum () resolution of ARPES, close to the zone centre in measurements of the electronic dispersions along in-plane high-symmetry directions (Fig. 2(a)). The Nb-dominated states, on the other hand, yield clear spectral features close to their Fermi crossings. Kinks in their measured dispersion and a decrease in linewidth near the Fermi level point to relatively-strong electron-phonon coupling in this system. 12; 13; 11 It is these Nb-derived states that are known to host the largest energy gaps at the Fermi level arising from the CDW and superconducting instabilities in this system, 15; 14; 16; 13 although the origins of these have proved highly controversial. 8; 18; 15; 10; 14; 19; 16; 20; 13; 21; 22; 11; 17 To the best of our knowledge, all prior theoretical treatments assume these orders emerge from an electronic liquid of trivially spin-degenerate character. Indeed, standard expectations of group theory would say this must be the case for the centrosymmetric space group () of -NbSe.

In contrast, we directly observe pronounced spin polarisations of the underlying electronic states in spin-resolved energy distribution curves (spin-EDCs, see Methods) measured along the direction, shown in Fig. 2(c-f). This is particularly apparent close to the saddle point of these bands (Fig. 2(d)), where two clearly-separated peaks can be observed in spin-EDCs. The measured spin polarisation is almost entirely out-of-plane (Fig. 2(d) and Supplemental Fig. S1), with a sign that reverses between the two bands. We attribute this as arising from local inversion symmetry breaking within the individual layers that make up the bulk crystal structure. 5; 6; 23 For sufficiently weak interlayer interactions (a point we return to below), pronounced spin-orbit coupling characteristic of the transition metal can lift the spin-degeneracy of the states localised within each layer, that thus strongly “feel” the local inversion asymmetry. A layer-dependent sign change of this spin polarisation, mediated by the rotated layer stacking of the bulk crystal structure, would act to restore overall spin degeneracy as required by global inversion and time-reversal symmetries. Nonetheless, the electronic states retain strong layer-resolved spin polarisations. Depth-averaging probes would largely be insensitive to these, but photoemission is a surface sensitive technique. Incoherent superposition of photoelectrons emitted from neighbouring layers of the unit cell would lead to some suppression of the measured spin polarisation, while interference effects can further complicate this picture 5. Nonetheless, the extreme surface sensitivity at the photon energies used here, with an inelastic mean free path on the order of the inter-layer separation, renders us predominantly sensitive to the spin texture of the top-most layer of the unit cell.

Spin-valley locking. Our measurements show that such non-trivial spin textures in NbSe persist up to the Fermi level. This is evident in spin-resolved EDCs along the direction (Fig. 2(c-f)) as well as the spin polarisation of EDCs and a Fermi-level momentum distribution curve (MDC) along the direction (Fig. 3). The latter clearly reveals how the spin polarisation reverses sign for each pair of Fermi surface sheets centred on neighbouring and points. This is a natural consequence of time-reversal symmetry. Here, this results in a coupling of the spin to the so-called valley index, the quantum number which distinguishes - and -centred Fermi surfaces in NbSe (Fig. 3(e)). Such spin-valley coupling has recently been extensively investigated for the band extrema in monolayers of semiconducting TMDCs, 24; 25; 26 where it has not only been shown to lead to new physics, such as a valley Hall effect, 27 but also to offer potential for devices exploiting the valley pseudospin. 29; 28 Our observations here point to a pronounced role of spin-valley coupling also for the low-energy quasiparticle excitations of the metallic -structured TMDCs.

Figure 4: Interplay of interlayer interactions and intralayer inversion symmetry breaking. (a,b) Density-functional theory calculations of the (a) out-of-plane and (b) in-plane spin polarisation of the three-dimensional Fermi surface of NbSe projected onto the first layer of the unit cell. (c,d) Calculated electronic structure along (c) --- () and (d) --- () with and without SOC. (e) Corresponding influence of SOC on the Fermi surface contours in the --- and --- planes (shown throughout the full three-dimensional Brillouin zone in Supplemental Fig. S2). (f) Our experimental ARPES measurements of such -dependent Fermi surfaces ( to eV) are in good general agreement with the theoretical calculations including SOC, supporting a SOC-mediated suppression of interlayer hopping in the --- plane.

Critically, it is from this spin-polarised Fermi sea that electron-hole and electron-electron pairing interactions drive the formation of CDW order and superconductivity. The largest CDW gaps in NbSe are located on the zone-corner spin-valley locked Fermi surfaces. 13 The CDW wave vector is entirely in-plane, 30 making such hidden layer-dependent spin polarisations relevant: we thus conclude that these cannot be viewed as purely charge density wave states with no role of the spin degree of freedom. The underlying spin textures can also be expected to have important implications for superconductivity. Recent measurements on electrically-gated MoS, which is known to host a strong spin-valley locking in its band structure, 25; 26 have found evidence for unconventional so-called Ising superconductivity. 32; 31 In this, the pronounced spin-orbit field that induces the underlying spin texture of the electronic states pins the spin of Cooper pair electrons in the out-of-plane direction. This leads to upper critical fields dramatically exceeding the Pauli paramagnetic limit for a magnetic field applied within the ab-plane. Similar phenomenology has recently been reported for mono- and few-layer NbSe33 entirely consistent with our direct observation of spin-valley locking in this compound.

In bulk NbSe, the spin-layer locking, and also multi-band nature as compared to gated MoS, raises further prospects for stabilising a delicate balance between different pairing states. In the normal state, we find that the -centred Fermi surface barrels, which are known to support a modulated superconducting gap in the bulk, 13; 34 are also strongly spin-polarised along the direction (Fig. 2(c)). Their spin polarisation, however, is completely suppressed along (Fig. 2(b)). Intriguingly, the largest superconducting gaps of this Fermi surface determined in Ref. 13 are located at the regions of strongest spin polarisation evident here. Pronounced superconducting gaps are also known to occur on the zone-corner Fermi surfaces, 13; 15; 34 which, as demonstrated above, host strong spin-valley locking even in the bulk. Given the significant influence of SOC, a mixing of spin-triplet and spin-singlet order parameters could potentially be expected. 3 Forming even pseudospin-singlet Cooper pairs from the strongly spin-polarised Fermi surfaces could necessitate a phase locking between the order parameter of neighbouring Fermi surface sheets. Moreover, the -axis coherence length in NbSe is much greater than the interlayer separation. 35; 36 This raises the tantalising prospect that the inherent coupling between the spin and layer pseudospins reported here could be tuned to drive an instability to an odd-parity pair density wave state, where the sign of the superconducting gap becomes tied to the layer index. 7; 37 The proximity of bulk NbSe and similar compounds to such phases requires further theoretical exploration, and will depend sensitively on the relative importance of inter- and intra-band as well as interlayer pairing interactions.

Three-dimensional spin structure. Our density-functional theory calculations (Fig. 4) already reveal a key role of interlayer interactions in mediating and controlling the underlying spin texture of the normal-state bulk Fermi surface. The calculated spin-polarisation projected onto the first layer of the unit cell is shown throughout the full three-dimensional Brillouin zone in Fig. 4(a,b). The momentum-dependent spin polarisations are determined by the effective spin-orbit field, , where is a momentum-dependent scaling factor, is the net electrostatic potential gradient, and is the crystal momentum. The horizontal mirror symmetry of each NbSe structural unit about the transition-metal plane ( of the point group; Fig. 1(a)) ensures that is entirely within the plane. Due to the 180 relative rotation of neighbouring NbSe monolayers in the bulk crystal structure, , and thus , has opposite sign for successive layers in the unit cell. This causes the sign of the spin polarisation to reverse at each momentum-space point when projected onto layer 2 vs. layer 1 of the unit cell (Supplemental Fig. S3). This confirms that local inversion asymmetry within each NbSe layer drives the formation of the spin-polarised states observed here.

For electronic states whose wavefunctions are completely delocalised over both layers of the unit cell, the spin-orbit field from each layer cancels. Such states are consequently unpolarised, restoring the conventional expectations for a centrosymmetric space group. This can be observed for the highly three-dimensional “pancake” Fermi surface at the Brillouin zone centre here. In contrast, for electronic states predominantly localised within individual layers of the unit cell, strong layer-dependent spin-orbit fields mediate large layer-resolved spin polarisations. This is ideally realised at the three-dimensional Brillouin zone boundary along (). In a tight-binding picture, for this , interlayer hoppings within a unit cell and between neighbouring unit cells are out of phase with each other and thus cancel. 10 This can be directly visualised by comparing the dispersion of the electronic states along the --- and --- directions: neglecting SOC, a single four-fold degenerate band crosses the Fermi level in the plane (Fig. 4(c)), whereas this is split into a pair of two-fold degenerate bands by interlayer interactions for , as evident in Fig. 4(d).

With its forbidden interlayer coupling, the electronic structure for the plane is thus formally equivalent to that of an isolated monolayer. The spin-orbit field is consequently maximised, driving the largest (%) layer-dependent spin-polarisations. These are purely out-of-plane (Fig. 4(a,b)), as must lie entirely along at the Brillouin zone boundary along . Along the - direction, however, a vertical mirror plane in the crystal structure forbids any out-of-plane component of . The spin-orbit field must therefore be strictly zero along -. This enforces a touching, and hence spin-degeneracy, of the zone-centre Fermi surface barrels along this direction, which are otherwise strongly spin-polarised throughout the plane.

For other through the Brillouin zone, the spin-orbit field strength can be partially suppressed by finite interlayer coupling. Moreover, the -component of the momentum induces a non-zero component of in the plane. Together, this causes not only the magnitude of the Fermi surface spin polarisation, but also its vectorial spin texture, to develop a strong dependence on both the in- and out-of-plane momentum (Fig. 4(a,b)). Within the --- plane, the calculated Fermi surface crossings are relatively strongly dispersive in (Fig. 4(e)). This is consistent with our experimental measurements of the -dependent Fermi surface (Fig. 4(f)) and points to significant interlayer coupling. Even away from , the Fermi surface spin polarisations therefore become strongly suppressed in the entire vicinity of the --- plane.

For the --- plane, however, we find that the dispersion of the Fermi surface crossings are significantly reduced by the inclusion of spin-orbit coupling in the calculations (Fig. 4(e)), in keeping with our experimental measurements of more two-dimensional Fermi contours for this plane (Fig. 4(f)). This reduction in dispersion is achieved by a lifting, via SOC, of the four-fold degeneracy that would otherwise be present along the - line, and is accompanied by the emergence of strong spin-polarisation of the Fermi surface crossings for . This is in contrast to the - line, where the energetic degeneracy of the Fermi crossings is protected by their symmetry-enforced spin degeneracy. The reduction of -dispersion within the --- plane is equivalent to a spin-orbit-mediated suppression of interlayer hopping here. 28 This allows relatively strong spin polarisations to be maintained in the vicinity of this plane, with only a moderate suppression of the spin-orbit field strength away from .

For momenta close to the - line, the in-plane momentum is always much larger than the out-of-plane component. This causes to remain predominately aligned along . The spin polarisations of the zone-corner Fermi surfaces are thus largely out-of-plane throughout the full Brillouin zone (Fig. 4(a)). For the Nb-derived zone-centre Fermi surface barrels, however, much stronger (up to %) in-plane components develop (Fig. 4(b)). The in-plane spin texture is largely radial to the Fermi surface (Supplemental Fig. S4), and switches sign about the plane. This indicates a non-zero component of within the plane, with a direction that is tied to the sign of . It thus confirms that an in-plane component of the spin texture arises due to finite out-of-plane momentum, where no symmetry constraint exists to cancel the in-plane component of the spin-orbit field for low-symmetry momenta with . The existence of such hidden in-plane spin textures is a unique property of the bulk compound. Indeed, the in-plane spin component goes strictly to zero at the Brillouin zone boundary along , enforced by symmetry, and recovering a monolayer-like purely out-of-plane spin texture for . The bulk system thus hosts a rich intertwining of spin, orbital, and layer degrees of freedom, mediating a three-dimensional nature of its spin texture that can be expected to further modulate its pairing interactions. 31


Taken together, our results show how the combination of interlayer hopping and intra-layer inversion symmetry breaking can lead to particularly rich momentum-dependent spin textures of metallic TMDCs. Suppression of Fermi surface spin polarisation away from the Brillouin zone boundary along will be reduced with decreased interlayer interactions. Indeed, upper critical fields, already known to exceed the Pauli limit for bulk NbSe38 are dramatically enhanced in other -TMDC superconductors as a function of increasing interlayer separation.39 This points to a susceptibility of the bulk systems to Ising superconductivity similar to that recently observed in isolated monolayers. 33 The delicate balance between interlayer hopping and spin-orbit coupling strength in NbSe makes this an ideal material for understanding and ultimately controlling the role of layer-dependent spin polarisations on the collective states and phases of transition-metal dichalcogenides.

ARPES: Spin-resolved ARPES measurements were performed at the I3 beamline of MAX-IV Laboratory, Sweden, and ARPES measurements were performed using the I05 beamline of Diamond Light Source, UK. Measurements were performed at temperatures of  K using -polarised synchrotron light. Scienta R4000 hemispherical electron analysers were utilised for all measurements. This was additionally fitted with a mini-Mott detector scheme for the spin-ARPES measurements, configured to simultaneously probe the out-of-plane and in-plane (along the analyser slit direction) component of the photoelectron spin. 40 The finite spin-detection efficiency was corrected using a Sherman function of 40 and the spin-resolved EDCs determined according to

with , , the measured intensity on the individual detectors in the Mott scattering chamber, corrected by a relative detector efficiency calibration, and the total spin polarisation,

Calculations: Density-functional theory (DFT) calculations including spin-orbit coupling were performed using the modified Becke-Johnson exchange potential and Perdew-Burke-Ernzerhof correlation functional implemented in the WIEN2K program. 41 A 20x20x10 -mesh was employed. The DFT results were downfolded using maximally localised Wannier functions, 42; 43; 44 employing Nb orbitals and Se orbitals as basis functions. The resulting tight-binding Hamiltonian allows a direct extraction of the spin and layer projections of the electronic structure.

Acknowledgements We thank M.R. Lees for assistance with the transport measurements and C. Hooley, P. Wahl, A.P. Mackenzie, S. Lee, and B. Braunecker for useful discussions. We gratefully acknowledge support from the Engineering and Physical Sciences Research Council, UK (work at St Andrews under Grant Nos. EP/I031014/1 and EP/M023427/1 and work at Warwick under Grant No. EP/M028771/1) and the International Max Planck partnership. PDCK acknowledges support from the Royal Society through a University Research Fellowship. MSB was supported by the Grant-in-Aid for Scientific Research (S) (No. 24224009) from the Ministry of Education, Culture, Sports, Science and Technology (MEXT) of Japan. LB, JR, and VS acknowledge studentship funding from EPSRC through grant nos. EP/G03673X/1, EP/L505079/1, and EP/L015110/1, respectively. The experiments at MAX IV Laboratory were made possible through funding from the Swedish Research Council and the Knut and Alice Wallenberg Foundation. We also thank Diamond Light Source for access to beamline I05 (proposal no. SI11383) that contributed to the results presented here.


  • (1) Manchon, A. Koo, H. C., Nitta, J., Frolov, S. M. and Duine, R. A., New perspectives for Rashba spin-orbit coupling, Nature Mater. 14, 871-882 (2015).
  • (2) Hasan, M. Z. and Kane, C. L., Colloquium: Topological insulators, Rev. Mod. Phys. 85, 3045 (2010).
  • (3) Gor’kov, L. P. and Rashba, E. I., Superconducting 2D system with lifted spin degeneracy: Mixed singlet-triplet state, Phys. Rev. Lett. 87, 037004 (2001).
  • (4) Sato, M. and Fujimoto, S., Topological phases of noncentrosymmetric superconductors: Edge states, Majorana fermions, and non-Abelian statistics, Phys. Rev. B 79, 094504 (2009).
  • (5) Riley, J. M. et al., Direct observation of spin-polarized bulk bands in an inversion-symmetric semiconductor, Nature Phys. 10, 835-839 (2014).
  • (6) Zhang, X., Liu, Q., Luo, J.-W., Freeman, A. J., and Zunger, A., Hidden spin polarization in inversion-symmetric bulk crystals, Nature Phys. 10, 387-393 (2014).
  • (7) Sigrist, M. et al., Superconductors with staggered non-centrosymmetricity, J. Phys. Soc. Jpn. 83, 061014 (2014).
  • (8) Wilson, J. A., Di Salvo, F. J., and Mahajan, S., Charge-density waves in metallic, layered, transition-metal dichalcogenides, Phys. Rev. Lett. 32, 882 (1974).
  • (9) Rossnagel, K. et al., Fermi surface of 2-NbSe and its implications on the charge-density-wave mechanism, Phys. Rev. B 64, 235119 (2001).
  • (10) Johannes, M. D, Mazin, I. I, and Howells, C. A., Fermi-surface nesting and the origin of the charge-density wave in 2-NbSe, Phys. Rev. B 73, 205102 (2006).
  • (11) Flicker, F. and van Wezel, J., Charge order from orbital-dependent coupling evidenced by NbSe, Nature Commun. 6, 7034 (2015).
  • (12) Weber, F. et al., Electron-phonon coupling and the soft phonon mode in TiSe, Phys. Rev. Lett 107, 107403 (2011).
  • (13) Rahn, D. J. et al., Gaps and kinks in the electronic structure of the superconductor 2-NbSe from angle-resolved photoemission at  K Phys. Rev. B 85, 224532 (2012).
  • (14) Kiss, T. et al., Charge-order-maximized momentum-dependent superconductivity, Nature Phys. 3, 720-725 (2007).
  • (15) Yokoya, T., Kiss, T., Chainani, A., Shin, S., Nohara, M. and Takagi, H., Fermi surface sheet-dependent superconductivity in 2-NbSe, Science 294, 2518-2520 (2001).
  • (16) Borisenko, S. V. et al., Two energy gaps and Fermi-surface “arcs” in NbSe, Phys. Rev. Lett. 102, 166402 (2009).
  • (17) Xi, X. et al., Strongly enhanced charge-density-wave order in monolayer NbSe, Nature Nano. 10, 765-769 (2015).
  • (18) Rice, T. M. and Scott, G. K., New mechanism for a charge-density-wave instability, Phys. Rev. Lett. 35, 120 (1975).
  • (19) Shen, D. W. et al., Primary role of the barely occupied states in the charge density wave formation of NbSe, Phys. Rev. Lett. 101, 226406 (2008).
  • (20) Rossnagel, K., On the origin of charge-density waves in select layered transition-metal dichalcogenides, J. Phys. Cond. Mat. 23, 21 (2011).
  • (21) Soumyanarayanan, A. et al., Quantum phase transition from triangular to stripe charge order in NbSe P. Nat. Acad. Sci. USA 110, 1623-1627 (2013).
  • (22) Laverock, J. et al., -resolved susceptibility function of 2-TaSe from angle-resolved photoemission, Phys. Rev. B 88, 035108 (2013).
  • (23) Riley, J. M. et al., Negative electronic compressibility and tunable spin splitting in WSe, Nature Nano. 10, 1043-1047 (2015).
  • (24) Xiao, D., Liu, G.-B., Feng, W., Xu, X. and Yao, W., Coupled spin and valley physics in monolayers of MoS and other group-VI dichalcogenides, Phys. Rev. Lett. 108, 196802 (2012).
  • (25) Mak, K. F., He, K., Shan, J. and Heinz, T. F., Control of valley polarization in monolayer MoS by optical helicity Nature Nano. 7, 494-498 (2012).
  • (26) Zeng, H., Dai, J., Yao, W., Xiao, D. and Cui, X., Valley polarization in MoS monolayers by optical pumping, Nature Nano. 7, 490-493 (2012).
  • (27) Mak, K. F., McGill, K. L., Park, J., and McEuen, P. L., The valley Hall effect in MoS transistors, Science 344, 1489-1492 (2014).
  • (28) Gong, Z. et al., Magnetoelectric effects and valley-controlled spin quantum gates in transition metal dichalcogenide bilayers, Nature Comm. 4, 2053 (2013).
  • (29) Xu, X, Yao, W, Xiao, D. and Heinz, T. F., Spin and pseudospins in layered transition metal dichalcogenides, Nature Phys. 10, 343-350 (2014).
  • (30) Moncton, D. E. et al., Neutron scattering study of the charge-density wave transitions in -TaSe and -NbSe, Phys. Rev. B 16, 801 (1977).
  • (31) Saito, Y. et al., Superconductivity protected by spin-valley locking in gate-tuned MoS, Nature Phys. 12, 144-149 (2016).
  • (32) Lu, J. M. et al., Evidence for two-dimensional Ising superconductivity in gated MoS, Science 350, 1353-1357 (2015).
  • (33) Xi, X. et al., Ising pairing in superconducting NbSe atomic layers, Nature Phys. 12, 139-143 (2016).
  • (34) Huang, C. L. et al., Experimental evidence for a two-gap structure of superconducting NbSe: A specific-heat study in external magnetic fields, Phys. Rev. B 76, 212504 (2007).
  • (35) Sanchez, D., Junod, A., Muller, J., Berger, H. and Lévy, F., Specific heat of 2-NbSe in high magnetic fields, Physica B 204, 167-175 (1995).
  • (36) Nader, A. and Monceau, P., Critical field of 2-NbSe down to  mK, SpringerPlus 3, 16 (2014).
  • (37) Yoshida, T., Sigrist, M. and Yanase, Y., Pair-density wave states through spin-orbit coupling in multilayer superconductors, Phys. Rev. B 86, 134514 (2012).
  • (38) Foner, S. and McNiff Jr, E. J., Upper critical fields of layered superconducting NbSe at low temperature, Phys. Lett. A 45, 429-430 (1973).
  • (39) Klemm, R. A., Pristine and intercalated transition metal dichalcogenide superconductors, Physica C 514, 86-94 (2015).
  • (40) Berntsen, M. H. et al., A spin- and angle-resolving photoelectron spectrometer, Rev. Sci. Instrum. 81, 035104 (2010).
  • (41) Blaha, P. et al., WIEN2K package, Version 10.1 (2010); available at http://www.wien2k.at.
  • (42) Souza, I. et al., Maximally localized Wannier functions for entangled energy bands, Phys. Rev. B 65, 035109 (2001).
  • (43) Mostofi, A. A. et al., Wannier90: A Tool for obtaining maximally localised Wannier functions, Comp. Phys. Commun. 178, 685Ð699 (2008).
  • (44) Kuneš, J. et al., WIEN2WANNIER: From linearized augmented plane waves to maximally localized Wannier functions, Comp. Phys. Commun. 181, 1888-1895 (2010).
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description