The superfluid $^3$He formed by spin-triplet $p$-wave Cooper pairs is a typical topological superfluid. In the superfluid $^3$He B-phase, several kinds of vortices classified by spatial symmetries $P_1$, $P_2$, and $P_3$ are produced, where $P_1$ is
inversion symmetry, $P_2$ is magnetic reflection symmetry, and $P_3$ is magnetic $pi$-rotation symmetry. We have calculated the vortex bound states by the Bogoliubov-de Gennes theory and the quasiclassical Eilenberger theory, and also clarified symmetry protection of the low energy excitations by the spatial symmetries. On the symmetry protection, $P_3$ symmetry plays a key role which gives two-fold degenerate Majorana zero modes. Then, the bound states in the most symmetric $o$ vortex with $P_1$, $P_2$, and $P_3$ symmetries and in $w$ vortex with $P_3$ symmetry have the symmetry protected degenerate Majorana zero modes. On the other hand, zero energy modes in $v$ vortex, which is believed to be realized in the actual B-phase, are not protected, and in consequence become gapped by breaking axial symmetry. The excitation gap may have been observed as the variation of critical velocity. We have also suggested an experimental setup to create $o$ vortex with Majorana zero modes by a confinement and a magnetic field.
Motivated by recent experiments on heavy fermion materials CeCu$_2$Si$_2$ and UBe$_{13}$, we develop a framework to capture generic properties of multiband superconductors with strong Pauli paramagnetic effect (PPE). In contrast to the single band ca
se, the upper critical field $H_{rm c2}$ can remain second order transition even for strong PPE cases. The expected first order transition is hidden inside $H_{rm c2}$ and becomes a crossover due to the interplay of multibandness. The present theory based on full self-consistent solutions of the microscopic Eilenberger theory explains several mysterious anomalies associated with the crossover and the empty vortex core state which is observed by recent STM experiment on CeCu$_2$Si$_2$.
We investigate the topological aspect of the spin-triplet $f$-wave superconductor UPt$_3$ through microscopic calculations of edge- and vortex-bound states based on the quasiclassical Eilenberger and Bogoliubov-de Gennes theories. It is shown that a
gapless and linear dispersion exists at the edge of the $ab$-plane. This forms a Majorana valley, protected by the mirror chiral symmetry. We also demonstrate that, with increasing magnetic field, vortex-bound quasiparticles undergo a topological phase transition from topologically trivial states in the double-core vortex to zero-energy states in the normal-core vortex. As long as the $d$-vector is locked into the $ab$-plane, the mirror symmetry holds the Majorana property of the zero-energy states, and thus UPt$_3$ preserves topological crystalline superconductivity that is robust against the crystal field and spin-orbit interaction.
We propose a method utilizing edge current to observe Majorana fermions in the surface Andreev bound state for the superfluid $^3$He A- and B-phases. The proposal is based on self-consistent analytic solutions of quasi-classical Greens function with
an edge. The local density of states and edge mass current in the A-phase or edge spin current in the B-phase can be obtained from these solutions. The edge current carried by the Majorana fermions is partially cancelled by quasiparticles (QPs) in the continuum state outside the superfluid gap. QPs contributing to the edge current in the continuum state are distributed in energy even away from the superfluid gap. The effect of Majorana fermions emerges in the depletion of the edge current by temperature within a low-temperature range. The observations that the reduction in the mass current is changed by $T^2$-power in the A-phase and the reduction in the spin current is changed by $T^3$-power in the B-phase establish the existence of Majorana fermions. We also point out another possibility for observing Majorana fermions by controlling surface roughness.
Motivated by a recent angle-resolved thermal conductivity experiment that shows a twofold gap symmetry in the high-field and low-temperature C phase in the heavy-fermion superconductor UPt$_3$, we group-theoretically identify the pairing functions as
$E_{1u}$ with the $f$-wave character for all the three phases. The pairing functions are consistent with the observation as well as with a variety of existing measurements. By using a microscopic quasi-classical Eilenberger equation with the identified triplet pairing function under applied fields, we performed detailed studies of the vortex structures for three phases, including the vortex lattice symmetry, the local density of states, and the internal field distribution. These quantities are directly measurable experimentally by SANS, STM/STS, and NMR, respectively. It is found that, in the B phase of low $H$ and low $T$, the double-core vortex is stabilized over a singular vortex. In the C phase, thermal conductivity data are analyzed to confirm the gap structure proposed. We also give detailed comparisons of various proposed pair functions, concluding that the present scenario of $E_{1u}$ with the $f$-wave, which is an analogue to the triplet planar state, is better than the $E_{2u}$ or $E_{1g}$ scenario. Finally, we discuss the surface topological aspects of Majorana modes associated with the $E_{1u}^f$ state of planar like features.
The total angular momentum associated with the edge mass current flowing at the boundary in the superfluid $^3$He A-phase confined in a disk is proved to be $L=Nhbar/2$, consisting of $L^{rm MJ}=Nhbar$ from the Majorana quasi-particles (QPs) and $L^{
rm cont}=-Nhbar/2$ from the continuum state. We show it based on an analytic solution of the chiral order parameter for quasi-classical Eilenberger equation. Important analytic expressions are obtained for mass current, angular momentum, and density of states (DOS). Notably the DOS of the Majorana QPs is exactly $N_0/2$ ($N_0$: normal state DOS) responsible for the factor 2 difference between $L^{rm MJ}$ and $L^{rm cont}$. The current decreases as $E^{-3}$ against the energy $E$, and $L(T) propto -T^2$. This analytic solution is fully backed up by numerically solving the Eilenberger equation. We touch on the so-called intrinsic angular momentum problem.
