In neutron star matter, there exist $^{1}S_{0}$ superfluids in lower density in the crust while $^{3}P_{2}$ superfluids are believed to exist at higher density deep inside the core. In the latter, depending on the temperature and magnetic field, either the uniaxial nematic phase, the D$_{2}$-biaxial nematic phase, or the D$_{4}$-biaxial nematic phase appears. In this paper, we discuss a mixture of the $^{1}S_{0}$ and $^{3}P_{2}$ superfluids and find their coexistence. Adopting the loop expansion and the weak-coupling approximation for the interaction between two neutrons, we obtain the Ginzburg-Landau (GL) free energy in which both of the $^{1}S_{0}$ and $^{3}P_{2}$ condensates are taken into account by including the coupling terms between them. We analyze the GL free energy and obtain the phase diagram for the temperature and magnetic field. We find that the $^{1}S_{0}$ superfluid excludes the $^{3}P_{2}$ superfluid completely in the absence of magnetic field, they can coexist for weak magnetic fields, and the $^{1}S_{0}$ superfluid is expelled by the $^{3}P_{2}$ superfluid at strong magnetic fields, thereby proving the robustness of $^{3}P_{2}$ superfluid against the magnetic field. We further show that the D$_{4}$-BN phase covers the whole region of the $^{3}P_{2}$ superfluidity as a result of the coupling term, in contrast to the case of a pure $^{3}P_{2}$ superfluid studied before in which the D$_{4}$-BN phase is realized only under strong magnetic fields. Thus, the D$_{4}$-BN phase is topologically the most interesting phase, e.g., admitting half-quantized non-Abelian vortices relevant not only in magnetars but also in ordinary neutron stars.

We study surface effects of neutron $^{3}P_{2}$ superfluids in neutron stars. $^{3}P_{2}$ superfluids are in uniaxial nematic (UN), D$_{2}$ biaxial nematic (BN), or D$_{4}$ BN phase, depending on the strength of magnetic fields from small to large. We suppose a neutron $^{3}P_{2}$ superfluid in a ball with a spherical boundary. Adopting a suitable boundary condition for $^{3}P_{2}$ condensates, we solve the Ginzburg-Landau equation to find several surface properties for the neutron $^{3}P_{2}$ superfluid. First, the phase on the surface can be different from that of the bulk, and symmetry restoration or breaking occurs in general on the surface. Second, the distribution of the surface energy density has an anisotropy depending on the polar angle in the sphere, which may lead to the deformation of the geometrical shape of the surface. Third, the order parameter manifold induced on the surface, which is described by two-dimensional vector fields induced on the surface from the condensates, allows topological defects (vortices) on the surface, and there must exist such defects even in the ground state thanks to the Poincar\'{e}-Hopf theorem: although the numbers of the vortices and antivortices depend on the bulk phases, the difference between them is topologically invariant (the Euler number $\chi=2$) irrespective of the bulk phases. These vortices, which are not extended to the bulk, are called boojums in the context of liquid crystals and helium-3 superfluids. The surface properties of the neutron $^{3}P_{2}$ superfluid found in this paper may provide us useful information to study neutron stars.

We work out domain walls in neutron $^{3}P_{2}$ superfluids realized in the core of neutron stars. Adopting the Ginzburg-Landau (GL) theory as a bosonic low-energy effective theory, we consider configurations of domain walls interpolating ground states, i.e., the uniaxial nematic (UN), D$_{2}$-biaxial nematic (D$_{2}$-BN), and D$_{4}$-biaxial nematic (D$_{4}$-BN) phases in the presence of zero, small and large magnetic fields, respectively. We solve the Euler-Lagrange equation from the GL free energy density, and calculate surface energy densities of the domain walls. We find that one extra Nambu-Goldstone mode is localized in the vicinity of a domain wall in the UN phase while a U(1) symmetry restores in the vicinity of one type of domain wall in the D$_{2}$-BN phase and all domain walls in the D$_{4}$-BN phase. Considering a pile of domain walls in the neutron stars, we find that the most stable configurations are domain walls perpendicular to the magnetic fields piled up in the direction along the magnetic fields in the D$_{2}$-BN and D$_{4}$-BN phases. We estimate the energy released from the deconstruction of the domain walls in the edge of a neutron star, and show that it can reach an astrophysical scale such as glitches in neutron stars.

The Ginzburg-Landau (GL) effective theory is a useful tool to study a superconductivity or superfluidity near the critical temperature, and usually the expansion up to the 4th order in terms of order parameters is sufficient for the description of the second-order phase transition. In this paper, we discuss the GL equation for the neutron $^{3}P_{2}$ superfluidity relevant for interior of neutron stars. We derive the GL expansion up to the 8th order in the condensates and find that this order is necessary for the system to have the unique ground state, unlike the ordinary cases. Starting from the $LS$ potential, which provides the dominant attraction between two neutrons at the high density, we derive the GL equation in the path-integral formalism, where the auxiliary field method and the Nambu-Gor'kov representation are used. We present the detailed description for the trace calculation necessary in the derivation of the GL equation. As numerical results, we show the phase diagram of the neutron $^{3}P_{2}$ superfluidity on the plane spanned by the temperature and magnetic field, and find that the 8th order terms lead to a first-order phase transition, whose existence was predicted in the Bogoliubov-de Gennes equation but has not been found thus far within the framework of the GL expansion up to the 6th order.The first-order phase transition will affect the interior structures inside the neutron stars.

We discuss the effect of a strong magnetic field on neutron $^{3}P_{2}$ superfluidity. Based on the attraction in the $^{3}P_{2}$ pair of two neutrons, we derive the Ginzburg-Landau equation in the path-integral formalism by adopting the bosonization technique and leave the next-to-leading order in the expansion of the magnetic field $B$. We determine the $(T,B)$ phase diagram with temperature $T$, comprising three phases: the uniaxial nematic (UN) phase for $B=0$, D$_{2}$-biaxial nematic (BN) and D$_{4}$-BN phases in finite $B$ and strong $B$ such as magnetars, respectively, where D$_{2}$ and D$_{4}$ are dihedral groups. We find that, compared with the leading order in the magnetic field known before, the region of the D$_{2}$-BN phase in the $(T,B)$ plane is extended by the effect of the next-to-leading-order terms of the magnetic field. We also present the thermodynamic properties, such as heat capacities and spin susceptibility, and find that the spin susceptibility exhibits anisotropies in the UN, D$_{2}$-BN, and D$_{4}$-BN phases. This information will be useful to understand the internal structures of magnetars.

We discuss neutron $^{3}P_{2}$ phases in the core of neutron stars in strong magnetic field (magnetars). The neutron $^{3}P_{2}$ pairing provides a wide variety of condensates, such as the uniaxial nematic and (D$_{2}$ and D$_{4}$) biaxial nematic, with different symmetries stemming from the combinations of spin and momentum. Based on the spin-orbital angular momentum coupling and the spin-magnetic field coupling of the neutrons, we derive the Ginzburg-Landau equation containing higher order terms of the magnetic field. We investigate the phase diagram of the neutron $^{3}P_{2}$ superfluidity, and find that the D$_{2}$ biaxial nematic phase is extended by the higher order terms of the magnetic field. We also discuss the thermodynamic properties, the heat capacity and the spin susceptibility.

Quark-hadron continuity was proposed as crossover between hadronic matter and quark matter without a phase transition, based on the matching of the symmetry and excitations in both phases. In the limit of a light strange-quark mass, it connects hyperon matter and the color-flavor-locked (CFL) phase exhibiting color superconductivity. Recently, it was proposed that this conjecture could be generalized in the presence of superfluid vortices penetrating both phases (arXiv:1803.05115 [hep-ph]), and it was suggested that one hadronic superfluid vortex in hyperon matter could be connected to one non-Abelian vortex (color magnetic flux tube) in the CFL phase. Here, we argue that their proposal is consistent only at large distances; instead, we show that three hadronic superfluid vortices must combine with three non-Abelian vortices with different colors with the total color magnetic fluxes canceled out, where the junction is called a colorful boojum. We rigorously prove this in both a macroscopic theory based on the Ginzburg-Landau description in which symmetry and excitations match (including vortex cores), and a microscopic theory in which the Aharonov-Bohm phases of quarks around vortices match.

The quark-hadron crossover conjecture was proposed as a continuity between hadronic matter and quark matter with no phase transition. It is based on matching of symmetry and excitations in both the phases. It connects hyperon matter and color-flavor locked (CFL) phase of color superconductivity in the limit of light strange quark mass. We study generalization of this conjecture in the presence of topological vortices. We propose a picture where hadronic superfluid vortices in hyperon matter could be connected to non-Abelian vortices (color magnetic flux tubes) in the CFL phase during this crossover. We propose that three hadronic superfluid vortices must join together to three non-Abelian vortices with different color fluxes with the total color magnetic fluxes canceled out, where the junction is called a colorful boojum.

In this review, we discuss various properties of topological solitons in dense QCD matter, with a particular emphasis on the CFL phase exhibiting superfluidity and superconductivity, and their phenomenological implications in terms of the effective field theories such as the Ginzburg-Landau theory, the chiral Lagrangian, or the Bogoliubov--de Gennes equation. The most fundamental topological excitations are non-Abelian vortices, which are 1/3 quantized superfluid vortices and color magnetic flux tubes. They are created at a phase transition or a rotation such compact stars. The intervortex-interaction is repulsive and consequently a vortex lattice is formed. Bosonic and fermionic zero-energy modes are trapped in the vortex core and propagate along it as gapless excitations. The former consists of translational zero modes (a Kelvin mode) with a quadratic dispersion and CP(2) Nambu-Goldstone gapless modes with a linear dispersion, while the latter is the triplet Majorana fermion zero modes. The low-energy effective theory of the bosonic zero modes is a non-relativistic free complex scalar field and a CP(2) model in 1+1 dimensions. The effects of strange quark mass, electromagnetic interactions and non-perturbative quantum corrections are taken into account. Colorful boojums at the CFL interface, quantum color magnetic monopole confined by vortices, which supports the notion of quark-hadron duality, and Yang-Mills instantons inside a vortex as lumps are discussed. The interactions between a vortex and quasi-particles such as phonons, gluons, mesons, and photons are studied. A vortex lattice is shown to behave as a cosmic polarizer. Non-Abelian vortices are shown to behave as a novel kind of non-Abelian anyons. For the chiral symmetry breaking, we discuss fractional and integer axial domain walls, Abelian and non-Abelian axial vortices, axial wall-vortex composites, and Skyrmions.

We study the internal structure of a non-Abelian vortex in color superconductivity with respect to quark degrees of freedom. Stable non-Abelian vortices appear in the Color-Flavor-Locked phase whose symmetry SU(3)_{c+L+R} is further broken to SU(2)_{c+L+R} x U(1)_{c+L+R} at the vortex cores. Microscopic structure of vortices at scales shorter than the coherence length can be analyzed by the Bogoliubov-de Gennes (B-dG) equation (rather than the Ginzburg-Landau equation). We obtain quark spectra from the B-dG equation by treating the diquark gap having the vortex configuration as a background field. We find that there are massless modes (zero modes) well-localized around a vortex, in the triplet and singlet states of the unbroken symmetry SU(2)_{c+L+R} x U(1)_{c+L+R}. The velocities v_i of the massless modes (i=t,s for triplet and singlet) change at finite chemical potential \mu, and decrease as \mu becomes large. Therefore, low energy excitations in the vicinity of the vortices are effectively described by 1+1 dimensional massless fermions whose velocities are reduced v_i<1.