Two-fold symmetric superconductivity in few-layer NbSe2

The strong Ising spin-orbit coupling in certain two-dimensional transition metal dichalcogenides can profoundly affect the superconducting state in few-layer samples. For example, in NbSe2, this effect combines with the reduced dimensionality to stabilize the superconducting state against magnetic fields up to similar to 35 T, and could lead to topological superconductivity. Here we report a two-fold rotational symmetry of the superconducting state in few-layer NbSe2 under in-plane external magnetic fields, in contrast to the three-fold symmetry of the lattice. Both the magnetoresistance and critical field exhibit this two-fold symmetry, and it also manifests deep inside the superconducting state in NbSe2/CrBr3 superconductor-magnet tunnel junctions. In both cases, the anisotropy vanishes in the normal state, demonstrating that it is an intrinsic property of the superconducting phase. We attribute the behaviour to the mixing between two closely competing pairing instabilities, namely the conventional s-wave instability typical of bulk NbSe2 and an unconventional d- or p-wave channel that emerges in few-layer NbSe2. Our results demonstrate the unconventional character of the pairing interaction in few-layer transition metal dichalcogenides and highlight the exotic superconductivity in this family of two-dimensional materials.