We discuss the effect of a strong magnetic field on neutron 3 P2 superfluidity. Based on the attraction in the 3 P2 pair of two neutrons, we derive the Ginzburg-Landau equation in the pathintegral 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, D2-biaxial nematic (BN) and D4-BN phases in finite B and strong B such as magnetars, respectively, where D2 and D4 are dihedral groups. We find that, compared with the leading order in the magnetic field known before, the region of the D2-BN phase in the (T, B) plane is extended by the effect of the nextto-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, D2-BN, and D4-BN phases. This information will be useful to understand the internal structures of magnetars.
I. INTRODUCTIONNeutron stars are interesting astrophysical objects whose phenomena are induced by combinations of fundamental forces, i.e., strong interaction, weak interaction, electromagnetic interaction and gravitation (see Refs. [1,2] for recent reviews). Studying neutron stars by several different signals can open a new approach to unveil the interiors in neutron stars. It was epoch-making that gravitational waves from neutron star merger were observed directly [3]. One of the most important properties of neutron stars is accompanied by a strong magnetic field. It is known that the magnitude of the magnetic field is around 10 12 G in standard neutron stars, and, in magnetars, it can reach about 10 15 G at the surface and may reach even larger values in the inside (see Refs. [4,5] for reviews on magnetars). As for the origin of the strong magnetic fields, researchers have studied several mechanisms such as spin-dependent interactions between neutrons [6-9], 1 the pion domain wall [11,12], the spin polarization in quark-matter core [13][14][15] and so on. In neutron matter, the influence of the magnetic field on a neutron is provided by a finite magnetic moment, µ n = −(γ n /2)σ with the gyromagnetic ratio of a neutron, γ n = 1.2 × 10 −13 MeV/T (in natural units, = c = 1), and the Pauli matrices for the neutron spin σ. The interaction between the neutron and the magnetic field (B) supplies the energy splitting between the spin-up state and the spin-down state for a neutron, −µ n ·B. For example, the strong magnetic field |B| ∼ 10 15 G in magnetars gives the mass splitting about O(10) keV (notice the unit relation, 1 T = 10 4 G).The dominant ingredient in neutron stars is the neutron matter, and it exists as superfluidity with a small admixture of superconducting protons and normal electrons (see Refs. [16,17] for recent reviews). In relation to the observations of neutron stars, it is considered that the neutron superfluidity affects relaxation time afte...