We consider freeze-in production of 7 keV axino dark matter (DM) in the supersymmetric Dine-Fischler-Srednicki-Zhitnitsky (DFSZ) model in light of the 3.5 keV line excess. The warmness of such 7 keV DM produced from the thermal bath, in general, appears in tension with Ly-$\alpha$ forest data, although a direct comparison is not straightforward. This is because the Ly-$\alpha$ forest constraints are usually reported on the mass of the conventional warm dark matter (WDM), where large entropy production is implicitly assumed to occur in the thermal bath after WDM particles decouple. The phase space distribution of freeze-in axino DM varies depending on production processes and axino DM may alleviate the tension with the tight Ly-$\alpha$ forest constraints. By solving the Boltzmann equation, we first obtain the resultant phase space distribution of axinos produced by 2-body decay, 3-body decay, and 2-to-2 scattering, respectively. The reduced collision term and resultant phase space distribution are useful for studying other freeze-in scenarios as well. We then calculate the resultant linear matter power spectra for such axino DM and directly compare them with the linear matter power spectra for the conventional WDM. In order to demonstrate realistic axino DM production, we consider benchmark points with the Higgsino next-to-lightest supersymmetric particle (NLSP) and wino NLSP. In the case of the Higgsino NLSP, the phase space distribution of axinos is colder than that in the conventional WDM case, so the most stringent Ly-$\alpha$ forest constraint can be evaded with mild entropy production from saxion decay inherent in the supersymmetric DFSZ axion model.

Thermal evolution of neutron stars (NSs) is studied as a probe of physics beyond the standard model. We first review the standard cooling theory of NSs in detail, with an emphasis on the roles of nucleon superfluidity. Then we discuss non-standard evolution with axion and dark matter (DM); axion production enhances the cooling while DM accretion leads to heating of NSs. To evaluate the effect of DM heating, we need to compare it to the rotochemical heating, which is caused by the out-of-equilibrium beta processes in a NS. In the dissertation, we also investigate the rotochemical heating in the presence of both neutron and proton superfluidity.

Recent observations have found several candidates for old warm neutron stars whose surface temperatures are above the prediction of the standard neutron star cooling scenario, and thus require some heating mechanism. Motivated by these observations, we study the non-equilibrium beta process in the minimal cooling scenario of neutron stars, which inevitably occurs in pulsars. This out-of-equilibrium process yields the late time heating in the core of a neutron star, called the rotochemical heating, and significantly changes the time evolution of the neutron star surface temperature. To perform a realistic analysis of this heating effect, we include the singlet proton and triplet neutron pairing gaps simultaneously in the calculation of the rate and emissivity of this process, where the dependence of these pairing gaps on the nucleon density is also taken into account. We then compare the predicted surface temperature of neutron stars with the latest observational data. We show the simultaneous inclusion of both proton and neutron gaps is advantageous for the explanation of the old warm neutron stars since it enhances the heating effect. It is then found that the observed surface temperatures of the old warm millisecond pulsars, J2124-3358 and J0437-4715, are explained for various choices of nucleon gap models. The same setup is compatible with the observed temperatures of ordinary pulsars including old warm ones, J0108-1431 and B0950+08, by choosing the initial rotational period of each neutron star accordingly. In particular, the upper limit on the surface temperature of J2144-3933 can be satisfied if its initial period is $\gtrsim 10\,\mathrm{ms}$.

A feebly interacting massive particle (FIMP), contrasting with a weakly interacting massive particle (WIMP), is an intriguing dark matter (DM) candidate. Light (keV-scale) FIMP DM is of particular interest: its radiative decay leaves a line signal in x-ray spectra; and it is warm dark matter (WDM) and alters the galactic-scale structure formation of the Universe from that with WIMP DM. Once a possible x-ray line is reported (e.g., $3.5 \,\mathrm{keV}$ line and $7 \,\mathrm{keV}$ FIMP DM is inferred), one has to check whether or not this FIMP DM is compatible with the structure formation. Here is an issue: the structure formation constraint on WDM is often reported in terms of the so-called thermal WDM mass $m_{\mathrm{WDM}}$, which cannot be directly applied to FIMP parameters. In this paper, we introduce a benchmark FIMP model that represents well a broad class of FIMP models. A big advantage of this benchmark is that we can derive the analytic formula of the non-thermal phase space distribution of FIMPs produced from freeze-in processes. By further deriving a certain "warmness" quantity, we can analytically map $m_{\mathrm{WDM}}$ to FIMP parameters. Our analytic map indicates that $7 \,\mathrm{keV}$ FIMP DM, without entropy production or a degenerate spectrum, is in tension with the latest Lyman-$\alpha$ forest data. Our analytic map will be very useful for future updates of observational constraints and reports of x-ray lines. It is also very easy to incorporate our analytic formula into a Boltzmann solver so that a linear matter power spectrum is readily accessible. Our benchmark model will facilitate FIMP searches and particle physics model-building.

Dark matter (DM) particles in the Universe accumulate in neutron stars (NSs) through their interactions with ordinary matter. It has been known that their annihilation inside the NS core causes late-time heating, with which the surface temperature becomes a constant value of $T_s \simeq (2-3) \times 10^3$ K for the NS age $t \gtrsim 10^{6-7}$ years. This conclusion is, however, drawn based on the assumption that the beta equilibrium is maintained in NSs throughout their life, which turns out to be invalid for rotating pulsars. The slowdown in the pulsar rotation drives the NS matter out of beta equilibrium, and the resultant imbalance in chemical potentials induces late-time heating, dubbed as rotochemical heating. This effect can heat a NS up to $T_s \simeq 10^6$ K for $t \simeq 10^{6-7}$ years. In fact, recent observations found several old NSs whose surface temperature is much higher than the prediction of the standard cooling scenario and is consistent with the rotochemical heating. Motivated by these observations, in this letter, we reevaluate the significance of the DM heating in NSs, including the effect of the rotochemical heating. We then show that the signature of DM heating can still be detected in old ordinary pulsars, while it is concealed by the rotochemical heating for old millisecond pulsars. To confirm the evidence for the DM heating, however, it is necessary to improve our knowledge on nucleon pairing gaps as well as to evaluate the initial period of the pulsars accurately. In any cases, a discovery of a very cold NS can give a robust constraint on the DM heating, and thus on DM models. To demonstrate this, as an example, we also discuss the case that the DM is the neutral component of an electroweak multiplet, and show that an observation of a NS with $T_s \lesssim 10^3$ K imposes a stringent constraint on such a DM candidate.

Galactic-scale structure is of particular interest since it provides important clues to dark matter properties and its observation is improving. Weakly interacting massive particles (WIMPs) behave as cold dark matter on galactic scales, while beyond-WIMP candidates suppress galactic-scale structure formation. Suppression in the linear matter power spectrum has been conventionally characterized by a single parameter, the thermal warm dark matter mass. On the other hand, the shape of suppression depends on the underlying mechanism. It is necessary to introduce multiple parameters to cover a wide range of beyond-WIMP models. Once multiple parameters are introduced, it becomes harder to share results from one side to the other. In this work, we propose adopting neural network technique to facilitate the communication between the two sides. To demonstrate how to work out in a concrete manner, we consider a simplified model of light feebly interacting massive particles.

The observed rapid cooling of the neutron star (NS) located at the center of the supernova remnant Cassiopeia A (Cas A) can be explained in the minimal NS cooling scenario. This consequence may be changed if there exists an extra cooling source, such as axion emission. In this work, we study the Cas A NS cooling in the presence of axion emission, taking account of the temperature evolution in the whole life of the Cas A NS. We obtain a lower limit on the axion decay constant, $f_a\gtrsim (5-7)\times 10^8$ GeV, if the star has an envelope with a thin carbon layer. This is as strong as existing limits imposed by other astrophysical observations such as SN1987A.

We construct a self-interacting dark matter model that could simultaneously explain the observed muon anomalous magnetic moment. It is based on a gauged U$(1)_{L_{\mu} - L_{\tau}}$ extension of the standard model, where we introduce a vector-like pair of fermions as the dark matter candidate and a new Higgs boson to break the symmetry. The new gauge boson has a sizable contribution to muon $(g-2)$, while being consistent with other experimental constraints. The U$(1)_{L_{\mu} - L_{\tau}}$ Higgs boson acts as a light force carrier, mediating dark matter self-interactions with a velocity-dependent cross section. It is large enough in galaxies to thermalize the inner halo and explain the diverse rotation curves and diminishes towards galaxy clusters. Since the light mediator dominantly decays to the U$(1)_{L_{\mu} - L_{\tau}}$ gauge boson and neutrinos, the astrophysical and cosmological constraints are weak. We study the thermal evolution of the model in the early Universe and derive a lower bound on the gauge boson mass.

We point out that 7 keV axino dark matter (DM) in the R-parity violating (RPV) supersymmetric (SUSY) Dine-Fischler-Srednicki-Zhitnitsky model can simultaneously reproduce the 3.5keV X-ray excess, and evade stringent constraints from the Ly-alpha forest data. Peccei-Quinn symmetry breaking naturally generates both axino interactions with minimal SUSY standard model particles and RPV interactions. The RPV interaction introduces an axino-neutrino mixing and provides axino DM as a variant of sterile neutrino DM, whose decay into a monochromatic photon can be detected by X-ray observations. Axinos, on the other hand, are produced by freeze-in processes of thermal particles in addition to the Dodelson-Widrow mechanism of sterile neutrinos. The resultant phase space distribution tends to be colder than the Fermi-Dirac distribution. The inherent entropy production from late-time saxion decay makes axinos even colder. The linear matter power spectrum satisfies even the latest and strongest constraints from the Ly-alpha forest data.