Orbital-resolved vortex core states in FeSe Superconductors: calculation based on a three-orbital model

We study electronic structure of vortex core states of FeSe superconductors based on a t$_{2g}$ three-orbital model by solving the Bogoliubov-de Gennes(BdG) equation self-consistently. The orbital-resolved vortex core states of different pairing symmetries manifest themselves as distinguishable structures due to different quasi-particle wavefunctions. The obtained vortices are classified in terms of the invariant subgroups of the symmetry group of the mean-field Hamiltonian in the presence of magnetic field. Isotropic $s$ and anisotropic $s$ wave vortices have $G_5$ symmetry for each orbital, whereas $d_{x^2-y^2}$ wave vortices show $G^{*}_{6}$ symmetry for $d_{xz/yz}$ orbitals and $G^{*}_{5}$ symmetry for $d_{xy}$ orbital. In the case of $d_{x^2-y^2}$ wave vortices, hybridized-pairing between $d_{xz}$ and $d_{yz}$ orbitals gives rise to a relative phase difference in terms of gauge transformed pairing order parameters between $d_{xz/yz}$ and $d_{xy}$ orbitals, which is essentially caused by a transformation of co-representation of $G^{*}_{5}$ and $G^{*}_{6}$ subgroup. The calculated local density of states(LDOS) of $d_{x^2-y^2}$ wave vortices show qualitatively similar pattern with experiment results. The phase difference of $\frac{\pi}{4}$ between $d_{xz/yz}$ and $d_{xy}$ orbital-resolved $d_{x^2-y^2}$ wave vortices can be verified by further experiment observation.