Spin-orbital Yu-Shiba-Rusinov states in single Kondo molecular magnet

Studies of single-spin objects are essential for designing emergent quantum states. We investigate a molecular magnet Tb$_2$Pc$_3$ interacting with a superconducting Pb(111) substrate, which hosts unprecedented Yu-Shiba-Rusinov (YSR) subgap states, dubbed spin-orbital YSR states. Upon adsorption of the molecule on Pb, the degeneracy of its lowest unoccupied molecular orbitals (LUMO) is lifted, and the lower LUMO forms a radical spin via charge transfer. This leads to Kondo screening and subgap states. Intriguingly, the YSR states display two pairs of resonances with clearly distinct behavior. The energy of the inner pair exhibits prominent inter and intra molecular variation, and it strongly depends on the tip height. The outer pair, however, shifts only slightly. As is unveiled through theoretical calculations, the two pairs of YSR states originate from the ligand spin and charge-fluctuating higher LUMO, coexisting in a single molecule, but only weakly coupled presumably due to different spatial distribution. Our work paves the way for understanding complex many-body excitations and constructing molecule-based topological superconductivity.