Abstract
Conventional two-dimensional superconductivity is destroyed when the critical in-plane magnetic field exceeds the so-called Pauli limit. Some monolayer transition-metal dichalcogenides lack inversion symmetry and the strong spin-orbit coupling leads to a valley-dependent Zeeman-like spin splitting. The resulting spin-valley locking lifts the valley degeneracy and results in a strong enhancement of the in-plane critical magnetic field. In these systems, it was predicted that the density of states in an in-plane field exhibits distinct mirage gaps at finite energies of about the spin-orbit coupling strength, which arise from a coupling of the electron and hole bands at energy larger than the superconducting gap. In this study, we investigate the impact of a triplet pairing channel on the spectral properties, primarily the mirage gap and the superconducting gap, in the clean limit. Notably, in the presence of the triplet pairing channel, the mirage-gap width is reduced for the low magnetic fields. Furthermore, when the temperature is lower than the triplet critical temperature, the mirage gaps survive even in the strong-field limit due to the finite singlet and triplet order parameters. Our work provides insights into controlling and understanding the properties of spin-triplet Cooper pairs.
1 Introduction
Superconductivity in the presence of magnetism has been a topic of great interest within the scientific community for several decades, and is responsible for many exotic properties. Usually, an external magnetic field destroys superconductivity by aligning the electron spins in the same direction. This is typical in conventional superconductors where Cooper pairs are formed by two electrons with opposite spins. The upper critical field is limited by the Pauli paramagnetic effect ; . However, in some superconductors, the Pauli limit can be surpassed by forming Fulde-Ferrell-Larkin-Ovchinnikov (FFLO) states or by creating spin-triplet Cooper pairs so that the parallel-aligned spin configuration in Cooper pairs is not affected by Pauli paramagnetism ; ; . The spin-triplet Cooper pairings are important for spintronics applications as they can generate and control spin currents. They are also proposed as a possible route to realize topological superconductivity which could be used to build robust quantum computers . Therefore, understanding the properties of spin-triplet Cooper pairs is of great interest.
The transition-metal dichalcogenides serve as a platform for both exploring spin-triplet pairing physics and surpassing the Pauli limit in high-field superconductivity. It has been predicted that equal-spin triplet pairs can form in a superconducting few-layer 2H-NbSe2 when an in-plane magnetic field is applied . Unlike conventional superconductors, monolayer transition metal dichalcogenides, such as NbSe2, lack in-plane crystal inversion symmetry ; ; ; . This results in a spin-orbit interaction that generates an effective out-of-plane magnetic field, causing the electron spins to point out of the plane ; . Therefore, it is termed as Ising spin-orbit coupling (ISOC) ; ; , and the corresponding Ising superconductivity was experimentally found in numerous transition-metal dichalcogenides ; ; ; ; ; ; ; ; ; ; ; ; ; ; ; . ISOC is dependent on momentum and has opposite signs at the K and K′ points of the hexagonal Brillouin zone. It prevents the spin directions from being realigned by an externally applied in-plane magnetic field, thus overcoming the Pauli limit to demonstrate high in-plane critical fields ; .
Ising superconductors subjected to an in-plane magnetic field display unique features in their density of states, notably the emergence of additional half-gaps at finite energies of about the ISOC strength . These newly discovered gaps, called mirage gaps, represent a mirroring of the main superconducting gap and are signatures of the equal-spin triplet finite-energy pairing correlations. Their width is determined by the interplay of the in-plane magnetic field and ISOC. In the previous work, the mirage gaps have only been studied in the context of the singlet pairing channel . However, a recent experiment found that the superconducting gap in a few-layer NbSe2 under an in-plane magnetic field was larger than predictions based solely on the singlet-pairing channel . This discrepancy was attributed to the existence of a triplet-pairing channel, in which an equal-spin triplet order parameter couples with a singlet one to enhance the critical magnetic field .
In this work, we investigate the mirage gaps of an Ising superconductor that consists of both singlet and triplet pairing channels ; . For a fixed temperature, the maximal mirage-gap width by varying the magnetic field decreases with increasing the critical temperature of the triplet pairing channel Tct. When the magnetic field is extremely high and the temperature is lower than Tct, the mirage-gap width remains finite due to the nonvanishing spin-singlet and spin-triplet order parameters. This contrasts with the case without considering the triplet-pairing channel where the mirage gaps always disappear at the critical field for the singlet order parameter .
2 Model and formalism
An Ising superconductor with both a singlet order parameter Δs and an equal-spin triplet order parameter Δt can be described by a Bogoliubov–de Gennes Hamiltonian near the K (K′) valley by neglecting the contribution from Γ point . By applying an in-plane magnetic field B, the effective Hamiltonian can be written in the Nambu basis aswhere H0k) is given byHere, s = +1 corresponds to the valley K and s = −1 to K′. The deviation of the momentum from K or K′ is denoted by p. The Pauli matrices σx, σy, and σz act on the spin space and σ0 is the corresponding unit matrix. The dispersion measured from the chemical potential is ξp = p2/2m − μ. The ISOC strength is denoted by βso. The Zeeman term arising from the in-plane magnetic field in the x-direction is Bxσx which absorbs the factor of gμB/2 containing the Landé g factor and the Bohr magneton μB. Note that the ISOC forces the spins to align out of the plane, whereas the in-plane magnetic field aims to align the spins within the plane. The superconducting order parameter can be written as ; .The eigenvalues of the Bogoliubov-de Gennes Hamiltonian are given bywith and P = BxΔs − βsoΔt. The position ϵ0 and the width δ of a mirage gap are obtained from the eigenvalues of the Hamiltonian, being given by ϵ0 = ±(ϵ1 + ϵ2)/2 and δ = ϵ1 − ϵ2, respectively. At ξp = 0 and in the limit of βso ≫Δs, Δt, we have ϵ0 ≈ ρ andThe approximation indicates that the mirage-gap width is influenced by the interplay of the singlet and triplet order parameters.
Since the superconducting gap and the ISOC are much smaller compared to the Fermi energy, it allows us to use the formalism of quasiclassical Green’s function ; ; ; ; . The structure of the Green’s functions can be written asThe bar operation in the above expression is defined as with q ∈ {g0, f0, g, f}. For a homogeneous system in the clean limit, the Eilenberger equation can be written asThe order parameter term consisting of both the singlet and triplet components can be written aswith Δt = (0, isΔt, 0). The Zeeman and ISOC fields are included in the term asThe Pauli matrices τ1, τ2, and τ3 act on the Nambu space and τ0 is the corresponding unit matrix. Using the notation , we can write the Eilenberger equation as the following set of equationsThe normalization condition then gives us the following equationsBy combining the Eilenberger equation with and the normalization condition , we can get the components of the Green’s function. Defining and , we can express g0, g+,x, g−,z, f0, fx, and fy, respectively, asWe also find that , , , and . The remaining Green’s functions fz, g+,z, g−,x and g−,y are zero. We can use Green’s functions f0 in Eq. 18 and fy in Eq. 22 to derive the self-consistent equations for the singlet and the triplet order parameters. The cut-off frequencies Ωs(t) are related to the critical temperatures Tcs(t) and the coupling constants νs(t) by Ωs(t) = 1.764 Tcs(t) sinh(1/νs(t)). Assuming the same cut-off frequency for the singlet and the triplet pairing channels, we obtain νt = νs/[1 + νs ln(Tcs/Tct)]. The singlet and triplet order parameters are coupled with the relationsand
At very high magnetic field and in the limit of Bx ≫ βso, ϵ, Δs, Δt, one can show that , and . This leads to Δs = νsβsoΔt/(νtBx) so that P = (νs/νt − 1)βsoΔt. Therefore, at high magnetic fields, the mirage gaps remain finite as long as the order parameters are finite.
3 Numerical results
In the numerical calculation, the ISOC strength is set as βso = 7Tcs. In Figure 1A, we plot the singlet (Δs) and triplet (Δt) order parameters versus the in-plane magnetic field Bx by varying the triplet critical temperature Tct at a fixed temperature T. The behavior of the order parameters in the cases of T > Tct and T < Tct will be discussed separately in the following. For the case where T > Tct, with increasing magnetic field, the singlet order parameter Δs decreases and finally vanishes at the critical field. This is attributed to the pair-breaking effect of the magnetic field which tries to align the spins of the spin-singlet Cooper pairs in its direction. On the other hand, the in-plane magnetic field induces the equal-spin triplet pairings which are coupled to the singlet pairings. The triplet order parameter Δt first increases with increasing the in-plane magnetic field since the field attempts to align the spins in the plane which is favorable for the formation of triplet pairs. It then decreases due to its coupling with Δs as seen from the self-consistent gap equations. Finally, both Δs and Δt vanish at the same critical magnetic field for T > Tct. For the case where T < Tct, at small magnetic fields, Δs decreases, and Δt begins to increase in a similar manner to the case where T > Tct. The behavior of both order parameters differ from that of T > Tct case at high magnetic fields. Both Δs and Δt reach saturated values instead of vanishing at a critical field. This is because the triplet order parameter is preserved by the in-plane field as the field favors the formation of triplet pairs. Consequently, Δs does not vanish due to its coupling with Δt.
FIGURE 1
Figure 1B shows the mirage-gap width δ for different triplet critical temperatures Tct. The maximal mirage-gap width by varying the magnetic field decreases with increasing Tct. The mirage gaps vanish at the critical in-plane field when T > Tct. This is due to the dependence of mirage-gap width on the order parameters which vanish at the critical field. On the other hand, the mirage gaps are finite even at extremely high in-plane fields in the scenarios where T < Tct. This is because the order parameters remain finite in the high field limit at T < Tct.
4 Conclusion
To conclude, we have studied the effect of a triplet pairing channel on the spectral properties of an Ising superconductor. The presence of a triplet order parameter reduces the maximal mirage-gap width by varying the in-plane magnetic field. Notably, at a temperature lower than Tct, both the order parameters and the mirage gap are finite even at very high fields. In contrast, for the case where the temperature is higher than Tct, both the order parameters and the mirage gaps vanish at the critical field. From an experimental standpoint, it is important to note that the intervalley scattering influences both the order parameters, effectively eliminating equal-spin triplet pairing under moderate intervalley disorder . Our work provides a better understanding about the spin-triplet pairings in Ising superconductors.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Author contributions
SP: Formal Analysis, Investigation, Methodology, Validation, Writing–original draft, Writing–review and editing. GT: Conceptualization, Formal Analysis, Investigation, Project administration, Writing–original draft, Writing–review and editing. WB: Conceptualization, Formal Analysis, Funding acquisition, Methodology, Supervision, Validation, Writing–review and editing.
Funding
We acknowledge funding by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation)–Project-ID 443404566—SPP 2244.
Acknowledgments
We acknowledge useful discussions with M. Aprili, C. Bruder, A. Di Bernardo, and D. Nikolić.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AiL.ZhangE.HuangC.XieX.YangY.JiaZ.et al (2021). Van der Waals ferromagnetic Josephson junctions. Nat. Commun.12, 6580. 10.1038/s41467-021-26946-w
2
AokiD.FlouquetJ. (2012). Ferromagnetism and superconductivity in uranium compounds. J. Phys. Soc. Jpn.81, 011003. 10.1143/JPSJ.81.011003
3
AokiD.HuxleyA.RessoucheE.BraithwaiteD.FlouquetJ.BrisonJ. P.et al (2001). Coexistence of superconductivity and ferromagnetism in URhGe. Nature413, 613–616. 10.1038/35098048
4
BelzigW.WilhelmF. K.BruderC.SchönG.ZaikinA. D. (1999). Quasiclassical Green’s function approach to mesoscopic superconductivity. Superlattices Microstruct.25, 1251–1288. 10.1006/spmi.1999.0710
5
ChandrasekharB. S. (1962). A note on the maximum critical field of high-field superconductors. Appl. Phys. Lett.1, 7–8. 10.1063/1.1777362
6
ChoC.-w.LyuJ.AnL.HanT.LoK. T.NgC. Y.et al (2022). Nodal and nematic superconducting phases in NbSe2 monolayers from competing superconducting channels. Phys. Rev. Lett.129, 087002. 10.1103/physrevlett.129.087002
7
ClogstonA. M. (1962). Upper limit for the critical field in hard superconductors. Phys. Rev. Lett.9, 266–267. 10.1103/PhysRevLett.9.266
8
CostanzoD.ZhangH.ReddyB. A.BergerH.MorpurgoA. F. (2018). Tunnelling spectroscopy of gate-induced superconductivity in MoS2. Nat. Nanotechnol.13, 483–488. 10.1038/s41565-018-0122-2
9
de la BarreraS. C.SinkoM. R.GopalanD. P.SivadasN.SeylerK. L.WatanabeK.et al (2018). Tuning Ising superconductivity with layer and spin-orbit coupling in two-dimensional transition-metal dichalcogenides. Nat. Commun.9, 1427. 10.1038/s41467-018-03888-4
10
DvirT.MasseeF.AttiasL.KhodasM.ApriliM.QuayC. H. L.et al (2018). Spectroscopy of bulk and few-layer superconducting NbSe2 with van der Waals tunnel junctions. Nat. Commun.9, 598. 10.1038/s41467-018-03000-w
11
EilenbergerG. (1968). Transformation of Gorkov’s equation for type II superconductors into transport-like equations. Zeitschrift für Physik A Hadrons Nucl.214, 195–213. 10.1007/BF01379803
12
EschrigM. (2015). Spin-polarized supercurrents for spintronics: a review of current progress. Rep. Prog. Phys.78, 104501. 10.1088/0034-4885/78/10/104501
13
FrigeriP. A.AgterbergD. F.KogaA.SigristM. (2004). Superconductivity without inversion symmetry: MnSi versus CePt3Si. Phys. Rev. Lett.92, 097001. 10.1103/PhysRevLett.92.097001
14
FrolovS. M.ManfraM. J.SauJ. D. (2020). Topological superconductivity in hybrid devices. Nat. Phys.16, 718–724. 10.1038/s41567-020-0925-6
15
HamillA.HeischmidtB.SohnE.ShafferD.TsaiK.-T.ZhangX.et al (2021). Two-fold symmetric superconductivity in few-layer NbSe2. Nat. Phys.17, 949–954. 10.1038/s41567-021-01219-x
16
HuyN. T.GaspariniA.de NijsD. E.HuangY.KlaasseJ. C. P.GortenmulderT.et al (2007). Superconductivity on the border of weak itinerant ferromagnetism in UCoGe. Phys. Rev. Lett.99, 067006. 10.1103/PhysRevLett.99.067006
17
IdzuchiH.PientkaF.HuangK.-F.HaradaK.GülO.ShinY. J.et al (2021). Unconventional supercurrent phase in Ising superconductor Josephson junction with atomically thin magnetic insulator. Nat. Commun.12, 5332–5338. 10.1038/s41467-021-25608-1
18
IlićS.MeyerJ. S.HouzetM. (2017). Enhancement of the upper critical field in disordered transition metal dichalcogenide monolayers. Phys. Rev. Lett.119, 117001. 10.1103/PhysRevLett.119.117001
19
IlicS.MeyerJ. S.HouzetM. (2023). Spectral properties of disordered ising superconductors with singlet and triplet pairing in in-plane magnetic fields. 10.48550/arXiv.2308.02646
20
KangK.JiangS.BergerH.WatanabeK.TaniguchiT.ForróL.et al (2021). Giant anisotropic magnetoresistance in Ising superconductor-magnetic insulator tunnel junctions. 10.48550/arXiv.2101.01327
21
KopninN. (2001). Theory of nonequilibrium superconductivity. Oxford University Press.
22
KuzmanovićM.DvirT.LeBoeufD.IlićS.HaimM.MöckliD.et al (2022). Tunneling spectroscopy of few-monolayer NbSe2 in high magnetic fields: triplet superconductivity and Ising protection. Phys. Rev. B106, 184514. 10.1103/PhysRevB.106.184514
23
LarkinA. I.OvchinnikovY. N. (1969). Quasiclassical method in the theory of superconductivity. JETP 28, 1200–1205. zh. Eksp. Teor. Fiz.55, 2262–2272.
24
LiJ.SongP.ZhaoJ.VaklinovaK.ZhaoX.LiZ.et al (2021). Printable two-dimensional superconducting monolayers. Nat. Mat.20, 181–187. 10.1038/s41563-020-00831-1
25
LuJ. M.ZheliukO.LeermakersI.YuanN. F. Q.ZeitlerU.LawK. T.et al (2015). Evidence for two-dimensional Ising superconductivity in gated MoS2. Science350, 1353–1357. 10.1126/science.aab2277
26
LuJ.ZheliukO.ChenQ.LeermakersI.HusseyN. E.ZeitlerU.et al (2018). Full superconducting dome of strong Ising protection in gated monolayer WS2. Proc. Natl. Acad. Sci. U.S.A.115, 3551–3556. 10.1073/pnas.1716781115
27
MatsudaY.ShimaharaH. (2007). Fulde-Ferrell-Larkin-Ovchinnikov state in heavy fermion superconductors. J. Phys. Soc. Jpn.76, 051005. 10.1143/JPSJ.76.051005
28
MöckliD.KhodasM. (2020). Ising superconductors: interplay of magnetic field, triplet channels, and disorder. Phys. Rev. B101, 014510. 10.1103/PhysRevB.101.014510
29
OhnishiK.KomoriS.YangG.JeonK.-R.Olde OlthofL. A. B.MontielX.et al (2020). Spin-transport in superconductors. Appl. Phys. Lett.116. 10.1063/1.5138905
30
RamiresA. (2022). Nonunitary superconductivity in complex quantum materials. J. Phys. Condens. Matter34, 304001. 10.1088/1361-648X/ac6d3a
31
SaitoY.NakamuraY.BahramyM. S.KohamaY.YeJ.KasaharaY.et al (2016). Superconductivity protected by spin-valley locking in ion-gated MoS2. Nat. Phys.12, 144–149. 10.1038/nphys3580
32
SmidmanM.SalamonM. B.YuanH. Q.AgterbergD. F. (2017). Superconductivity and spin–orbit coupling in non-centrosymmetric materials: a review. Rep. Prog. Phys.80, 036501. 10.1088/1361-6633/80/3/036501
33
SohnE.XiX.HeW.-Y.JiangS.WangZ.KangK.et al (2018). An unusual continuous paramagnetic-limited superconducting phase transition in 2D NbSe2. Nat. Mat.17, 504–508. 10.1038/s41563-018-0061-1
34
TangG.BruderC.BelzigW. (2021). Magnetic field-induced “mirage” gap in an Ising superconductor. Phys. Rev. Lett.126, 237001. 10.1103/PhysRevLett.126.237001
35
WickramaratneD.KhmelevskyiS.AgterbergD. F.MazinI. I. (2020). Ising superconductivity and magnetism in NbSe2. Phys. Rev. X10, 041003. 10.1103/PhysRevX.10.041003
36
XiX.WangZ.ZhaoW.ParkJ.-H.LawK. T.BergerH.et al (2016). Ising pairing in superconducting NbSe2 atomic layers. Nat. Phys.12, 139–143. 10.1038/nphys3538
37
XiaoD.LiuG.-B.FengW.XuX.YaoW. (2012). Coupled spin and valley physics in monolayers of MoS2 and other group-VI dichalcogenides. Phys. Rev. Lett.108, 196802. 10.1103/PhysRevLett.108.196802
38
XingY.ZhaoK.ShanP.ZhengF.ZhangY.FuH.et al (2017). Ising superconductivity and quantum phase transition in macro-size monolayer NbSe2. Nano Lett.17, 6802–6807. 10.1021/acs.nanolett.7b03026
39
ZhouB. T.YuanN. F. Q.JiangH.-L.LawK. T. (2016). Ising superconductivity and Majorana fermions in transition-metal dichalcogenides. Phys. Rev. B93, 180501. 10.1103/PhysRevB.93.180501
40
ZhuZ. Y.ChengY. C.SchwingenschlöglU. (2011). Giant spin-orbit-induced spin splitting in two-dimensional transition-metal dichalcogenide semiconductors. Phys. Rev. B84, 153402. 10.1103/PhysRevB.84.153402
Summary
Keywords
Ising superconductors, spin-orbit coupling, mirage gaps, spin triplet pairing, unconventional superconductivity
Citation
Patil S, Tang G and Belzig W (2023) Spectral properties of a mixed singlet-triplet Ising superconductor. Front. Electron. Mater. 3:1254302. doi: 10.3389/femat.2023.1254302
Received
06 July 2023
Accepted
21 August 2023
Published
01 September 2023
Volume
3 - 2023
Edited by
Sachio Komori, Nagoya University, Japan
Reviewed by
Jianlin Luo, Chinese Academy of Sciences (CAS), China
Kyung-Hwan Jin, Institute for Basic Science (IBS), Republic of Korea
Updates
Copyright
© 2023 Patil, Tang and Belzig.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wolfgang Belzig, wolfgang.belzig@uni-konstanz.de
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.