We develop an innovative technique for studying inhomogeneous phases with a spontaneous broken symmetry. The method relies on the knowledge of the exact form of the free energy in the homogeneous phase and on a specific gradient expansion of the order parameter. We apply this method to quark matter at vanishing temperature and large chemical potential, which is expected to be relevant for astrophysical considerations. The method is remarkably reliable and fast as compared to performing the full numerical diagonalization of the quark Hamiltonian in momentum space, and is designed to improve the standard Ginzburg-Landau expansion close to the phase transition points. For definiteness we focus on inhomogeneous chiral symmetry breaking, accurately reproducing known results for 1D and 2D modulations and examining novel crystalline structures, as well. Consistently with previous results, we find that the energetically favored modulation is the so-called 1D real kink crystal. We propose a qualitative description of the pairing mechanism to motivate this result.

We consider the capture of dark matter (DM) in neutron stars via scattering on hadronic targets, including neutrons, protons and hyperons. We extend previous analyses by including momentum dependent form factors, which account for hadronic structure, and incorporating the effect of baryon strong interactions in the dense neutron star interior, rather than modelling the baryons as a free Fermi gas. The combination of these effects suppresses the DM capture rate over a wide mass range, thus increasing the cross section for which the capture rate saturates the geometric limit. In addition, variation in the capture rate associated with the choice of neutron star equation of state is reduced. For proton targets, the use of the interacting baryon approach to obtain the correct Fermi energy is essential for an accurate evaluation of the capture rate in the Pauli-blocked regime. For heavy neutron stars, which are expected to contain exotic matter, we identify cases where DM scattering on hyperons contributes significantly to the total capture rate. Despite smaller neutron star capture rates, compared to existing analyses, we find that the projected DM-nucleon scattering sensitivity greatly exceeds that of nuclear recoil experiments for a wide DM mass range.

Neutron star models with maximum mass close to 2 M⊙ reach high central densities, which may activate nucleonic and hyperon direct Urca neutrino emission. To alleviate the tension between fast theoretical cooling rates and thermal luminosity observations of moderately magnetized, isolated thermally-emitting stars (with Lγ ≳ 1031 erg s−1 at t ≳ 105.3 yr), some internal heating source is required. The power supplied by the internal heater is estimated for both a phenomenological source in the inner crust and Joule heating due to magnetic field decay, assuming different superfluidity models and compositions of the outer stellar envelope. It is found that a thermal power of W(t) ≈ 1034 erg s−1 allows neutron star models to match observations of moderately magnetized, isolated stars with ages t ≳ 105.3 yr. The requisite W(t) can be supplied by Joule heating due to crust-confined initial magnetic configurations with (i) mixed poloidal-toroidal fields, with surface strength Bdip = 1013 G at the pole of the dipolar poloidal component and ∼90 per cent of the magnetic energy stored in the toroidal component; and (ii) poloidal-only configurations with Bdip = 1014 G.

The dissipation of intense crustal electric currents produces high Joule heating rates in cooling neutron stars. Here it is shown that Joule heating can counterbalance fast cooling, making it difficult to infer the presence of hyperons (which accelerate cooling) from measurements of the observed thermal luminosity Lγ. Models with and without hyperon cores match Lγ of young magnetars (with poloidal-dipolar field Bdip ≳ 1014 G at the polar surface and Lγ ≳ 1034 erg s−1 at t ≲ 105 yr) as well as mature, moderately magnetized stars (with Bdip ≲ 1014 G and 1031 erg s−1 ≲ Lγ ≲ 1032 erg s−1 at t ≳ 105 yr). In magnetars, the crustal temperature is almost independent of hyperon direct Urca cooling in the core, regardless of whether the latter is suppressed or not by hyperon superfluidity. The thermal luminosities of light magnetars without hyperons and heavy magnetars with hyperons have Lγ in the same range and are almost indistinguishable. Likewise, Lγ data of neutron stars with Bdip ≲ 1014 G but with strong internal fields are not suitable to extract information about the equation of state as long as hyperons are superfluid, with maximum amplitude of the energy gaps of the order ≈1 MeV.