A strongly interacting massive particle (SIMP) is an interesting candidate for dark matter (DM) because its self-interaction cross section can be naturally strong enough to address the astrophysical problem of small-scale structure formation. A simple model was proposed by assuming a monopole condensation, where composite SIMP comes from a "strongly interacting" U(1)$_{\rm d}$ gauge theory. In the original model, the DM relic abundance is determined by the $3\to2$ annihilation process via the Wess-Zumino-Witten term. In this letter, we discuss that the DM relic abundance is naturally determined also by a semi-annihilation process via a kinetic mixing between the hypercharge gauge boson and the dark U(1)$_{\rm d}$ gauge boson (dark photon). The dark photon can be discovered by LDMX-style missing momentum experiments in the near future.

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.

(Mini) split supersymmetry explains the observed Higgs mass and evades stringent constraints, while keeping good features of TeV-scale supersymmetry other than the little hierarchy problem. Such scenarios naturally predict thermal wino dark matter whose mass is around $3 \, {\rm TeV}$. Its non-perturbatively enhanced annihilation is a promising target of indirect detection experiments. It is known that identifying the smallest halos is essential for reducing an uncertainty in interpreting indirect detection experiments. Despite its importance, the smallest halos of thermal wino dark matter have not been well understood and thus are investigated in this work. In particular, we remark on two aspects: 1) the neutral wino is in kinetic equilibrium with primordial plasma predominantly through inelastic processes involving the slightly heavier charged wino; and 2) the resultant density contrast shows larger powers at dark acoustic oscillation peaks than in cold dark matter, which is known as an overshooting phenomenon. By taking them into account, we provide a rigorous estimate of the boost factor. Our result facilitates accurately pinning down thermal wino dark matter through vigorous efforts in indirect detection experiments.

Exothermic scatterings of dark matter (DM) produce DM particles with significant kick velocities inside DM halos. In collaboration with DM self-interaction, the excess kinetic energy of the produced DM particles is distributed to the others, which self-heats the DM particles as a whole. The DM self-heating is efficient towards the halo center, and the heat injection is used to enhance the formation of a uniform density core inside halos. The effect of DM self-heating is expected to be more significant in smaller halos for two reasons: 1) the exothermic cross section times the relative velocity, $\left\langle\sigma_{\rm exo}v_{\rm rel}\right\rangle$, is constant; 2) and the ratio of the injected heat to the velocity dispersion squared gets larger towards smaller halos. For the first time, we quantitatively investigate the core formation from DM self-heating for halos in a wide mass range ($10^{9}$-$10^{15}\,{\rm M}_\odot$) using the gravothermal fluid formalism. Notably, we demonstrate that the core formation is sharply escalating towards smaller halos by taking the self-heating DM (i.e., DM that semi-annihilates and self-interacts) as an example. We show that the sharp escalation of core formation may cause a tension in simultaneously explaining the observed central mass deficit of Milky Way satellites, and field dwarf/low surface brightness spiral galaxies. While the details of the self-heating effect may differ among models, we expect that the sharp halo-mass dependence of the core formation is a general feature of exothermic DM.

We provide a simple UV theory for a Dirac dark matter with a massless Abelian gauge boson. We introduce a single fermion transforming as the $\bf{16}$ representation in the SO(10)$'$ gauge group, which is assumed to be spontaneously broken to SU(5)$'\times$U(1)$'$. The SU(5)$'$ gauge interaction becomes strong at an intermediate scale and then we obtain a light composite Dirac fermion with U(1)$'$ gauge interaction at the low-energy scale. Its thermal relic can explain the observed amount of dark matter consistently with other cosmological and astrophysical constraints. We discuss that a nonzero kinetic mixing between the U(1)$'$ gauge boson and the Hypercharge gauge boson is allowed and the temperature of the visible sector and the dark matter sector can be equal to each other.

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.

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.

A simple model of dark matter contains a light Dirac field charged under a hidden U(1) gauge symmetry. When a chiral matter content in a strong dynamics satisfies the t'Hooft anomaly matching condition, a massless baryon is a natural candidate of the light Dirac field. One realization is the same matter content as the standard SU(5)$\times$U(1)$_{(B-L)}$ grand unified theory. We propose a chiral [SU(5)$\times$U(1)]$^4$ gauge theory as a unified model of the SM and DM sectors. The low-energy dynamics, which was recently studied, is governed by the hidden U(1)$_4$ gauge interaction and the third-family U(1)$_{(B-L)_3}$ gauge interaction. This model can realize self-interacting dark matter and alleviate the small-scale crisis of collisionless cold dark matter in the cosmological structure formation. The model can also address the semi-leptonic $B$-decay anomaly reported by the LHCb experiment.

A spontaneously broken hidden U(1)$_h$ gauge symmetry can explain both the dark matter stability and the observed relic abundance. In this framework, the light gauge boson can mediate the strong dark matter self-interaction, which addresses astrophysical observations that are hard to explain in collisionless cold dark matter. Motivated by flavoured grand unified theories, we introduce right-handed neutrinos and a flavoured $B - L$ gauge symmetry for the third family U(1)$_{(B-L)_3}$. The unwanted relic of the U(1)$_h$ gauge boson decays into neutrinos via the kinetic mixing with the U(1)$_{(B - L)_3}$ gauge boson. Indirect detection bounds on dark matter are systematically weakened, since dark matter annihilation results in neutrinos. However, the kinetic mixing between U(1)$_{(B - L)_3}$ and U(1)$_Y$ gauge bosons are induced by quantum corrections and leads to an observable signal in direct and indirect detection experiments of dark matter. This model can also explain the baryon asymmetry of the Universe via the thermal leptogenesis. In addition, we discuss the possibility of explaining the lepton flavour universality violation in semi-leptonic $B$ meson decays that is recently found in the LHCb experiment.

The latest Lyman-$\alpha$ forest data severely constrain the conventional warm dark matter solution to small-scale issues in the cold dark matter paradigm. It has been also reported that unconstrained astrophysical processes may address the issues. In response to this situation, we revisit the decaying dark matter solution to the issues, discussing possible signatures to discriminate decaying dark matter from astrophysical processes as a solution to small-scale issues. We consider an axinolike particle (ALPino) decaying into an axionlike particle (ALP) and gravitino with the lifetime around the age of the Universe. The ALPino mass is sub-PeV and slightly ($\Delta m/m\sim 10^{-4}$) larger than the gravitino mass, and thus the dark matter abundance does not alter virtually after the ALPino decays. On the other hand, the gravitino produced from the ALPino decay obtains a kick velocity of $\sim 30 \,{\rm km / s}$, which is sufficiently larger than a circular velocity of dwarf galaxies to impact their dark matter distributions. The Lyman-$\alpha$ forest constraints are relieved since only a small fraction ($\sim10$%) of dark matter experiences the decay at that time. Decaying dark matter is thus promoted to a viable solution to small-scale issues. The ALPino relic abundance is determined predominantly by the decay of the lightest ordinary supersymmetric particle. The monochromatic ALP emission from the ALPino decay is converted to $\sim 50 \,{\rm GeV}$ photon under the Galactic magnetic field. The morphology of the gamma-ray flux shows a distinctive feature of the model when compared to decaying dark matter that directly decays into photons. Once detected, such distinctive signals discriminate the decaying dark matter solution to small-scale issues from unconstrained astrophysical processes.

We present results of a search for unidentified line emission and absorption signals in the 2-12keV energy band of spectra extracted from Perseus Cluster core region observations obtained with the 5eV energy resolution Hitomi Soft X-ray Spectrometer. No significant unidentified line emission or absorption is found. Line flux upper limits (1 sigma per resolution element) vary with photon energy and assumed intrinsic width, decreasing from 100 photons cm$^{-2}$ s$^{-1}$ sr $^{-1}$ at 2keV to $<10$ cm$^{-2}$ s$^{-1}$ sr $^{-1}$ over most of the 5-10 keV energy range for a Gaussian line with Doppler broadening of 640 km/s. Limits for narrower and broader lines have a similar energy dependence and are systematically smaller and larger, respectively. These line flux limits are used to constrain the decay rate of hypothetical dark matter candidates. For the sterile neutrino decay rate, new constraints over the the mass range of 4-24 keV with mass resolution better than any previous X-ray analysis are obtained. Additionally, the accuracy of relevant thermal spectral models and atomic data are evaluated. The Perseus cluster spectra may be described by a composite of multi-temperature thermal and AGN power-law continua. Superposed on these, a few line emission signals possibly originating from unmodeled atomic processes (including Si XIV and Fe XXV) are marginally detected and tabulated. Comparisons with previous X-ray upper limits and future prospects for dark matter searches using high-energy resolution spectroscopy are discussed.

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.

It was recently pointed out that semi-annihilating dark matter (DM) may experience a novel temperature evolution dubbed as self-heating. Exothermic semi-annihilation converts the DM mass to the kinetic energy. This yields a unique DM temperature evolution, $T_{\chi} \propto 1 / a$, in contrast to $ T_{\chi} \propto 1 / a^{2}$ for free-streaming non-relativistic particles. Self-heating continues as long as self-scattering sufficiently redistributes the energy of DM particles. In this paper, we study the evolution of cosmological perturbations in self-heating DM. We find that sub-GeV self-heating DM leaves a cutoff on the subgalactic scale of the matter power spectrum when the self-scattering cross section is $\sigma_{\rm self} / m_{\chi} \sim {\cal O} (1) \,{\rm cm}^{2} /{\rm g}$. Then we present a particle physics realization of the self-heating DM scenario. The model is based on recently proposed strongly interacting massive particles with pion-like particles in a QCD-like sector. Pion-like particles semi-annihilate into an axion-like particle, which is thermalized with dark radiation. The dark radiation temperature is smaller than the standard model temperature, evading the constraint from the effective number of neutrino degrees of freedom. It is easily realized when the dark sector is populated from the standard model sector through a small coupling.

The freeze-out of dark matter (DM) depends on the evolution of the DM temperature. The DM temperature does not have to follow the standard model one, when the elastic scattering is not sufficient to maintain the kinetic equilibrium. We study the temperature evolution of the semi-annihilating DM, where a pair of the DM particles annihilate into one DM particle and another particle coupled to the standard model sector. We find that the kinetic equilibrium is maintained solely via semi-annihilation until the last stage of the freeze-out. After the freeze-out, semi-annihilation converts the mass deficit to the kinetic energy of DM, which leads to non-trivial evolution of the DM temperature. We argue that the DM temperature redshifts like radiation as long as the DM self-interaction is efficient. We dub this novel temperature evolution as self-heating. Notably, the structure formation is suppressed at subgalactic scales like keV-scale warm DM but with GeV-scale self-heating DM if the self-heating lasts roughly until the matter-radiation equality. The long duration of the self-heating requires the large self-scattering cross section, which in turn flattens the DM density profile in inner halos. Consequently, self-heating DM can be a unified solution to apparent failures of cold DM to reproduce the observed subgalactic scale structure of the Universe.

We explore the phenomenology of having a second epoch of dark matter annihilation into dark radiation long after the standard thermal freeze-out. Such a hidden reannihilation process could affect visible sectors only gravitationally. As a concrete realization we consider self-interacting dark matter (SIDM) with a light force mediator coupled to dark radiation. We demonstrate how resonantly Sommerfeld enhanced cross sections emerge to induce the reannihilation epoch. The effect is a temporally local modification of the Hubble expansion rate and we show that the Cosmic Microwave Background (CMB) measurements -- as well as other observations -- have a high sensitivity to observe this phenomenon. Special attention is given to the model region where late kinetic decoupling and strong self-interactions can alleviate several small-scale problems in the cold dark matter paradigm at the same time. Interestingly, we find that reannihilation might here also simultaneously lower the tension between CMB and low-redshift astronomical observations of $H_{0}$ and $\sigma_{8}$. Moreover, we identify reannihilation as a clear signature to discriminate between the phenomenologically otherwise almost identical vector and scalar mediator realizations of SIDM.

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.

We propose a new realization of strongly interacting massive particles (SIMP) as self-interacting dark matter, where SIMPs couple to the Standard Model sector through an axion-like particle. Our model gets over major obstacles accompanying the original SIMP model, such as a missing mechanism of kinetically equilibrating SIMPs with the SM plasma as well as marginal perturbativity of the chiral Lagrangian density. Remarkably, the parameter region realizing $\sigma_{\rm self}/m_{\rm DM} \simeq 0.1 \textrm{--} 1 \, {\rm cm}^{2}/{\rm g}$ is within the reach of future beam dump experiments such as the Search for Hidden Particles (SHiP) experiment.

We investigate anomalous strong lens systems, particularly the effects of weak lensing by structures in the line of sight, in models with long-lived electrically charged massive particles (CHAMPs). In such models, matter density perturbations are suppressed through the acoustic damping and the flux ratio of lens systems are impacted, from which we can constrain the nature of CHAMPs. For this purpose, first we perform $N$-body simulations and develop a fitting formula to obtain non-linear matter power spectra in models where cold neutral dark matter and CHAMPs coexist in the early Universe. By using the observed anomalous quadruple lens samples, we obtained the constraints on the lifetime ($\tau_{\rm Ch}$) and the mass density fraction ($r_{\rm Ch}$) of CHAMPs. We show that, for $r_{\rm Ch}=1$, the lifetime is bounded as $\tau_{\rm Ch} < 0.96\,$yr (95% confidence level), while a longer lifetime $\tau_{\rm Ch} = 10\,$yr is allowed when $r_{\rm Ch} < 0.5$ at the 95% confidence level. Implications of our result for particle physics models are also discussed.

We argue that the acoustic damping of the matter power spectrum is not a generic feature of the kinetic decoupling of dark matter, but even the enhancement can be realized depending on the nature of the kinetic decoupling when compared to that in the standard cold dark matter model. We consider a model that exhibits a ${\it sudden}$ kinetic decoupling and investigate cosmological perturbations in the ${\it standard}$ cosmological background numerically in the model. We also give an analytic discussion in a simplified setup. Our results indicate that the nature of the kinetic decoupling could have a great impact on small scale density perturbations.

We examine the possibility that dark matter consists of charged massive particles (CHAMPs) in view of the cosmic microwave background (CMB) anisotropies. The evolution of cosmological perturbations of CHAMP with other components is followed in a self-consistent manner, without assuming that CHAMP and baryons are tightly coupled. We incorporate for the first time the "kinetic re-coupling" of the Coulomb scattering, which is characteristic of heavy CHAMPs. By a direct comparison of the predicted CMB temperature/polarization auto-correlations in CHAMP models and the observed spectra in the Planck mission, we show that CHAMPs leave sizable effects on CMB spectra if they are lighter than $10^{11}\,{\rm GeV}$. Our result can be applicable to any CHAMP as long as its lifetime is much longer than the cosmic time at the recombination ($\sim 4 \times 10^{5}\, {\rm yr}$). An application to millicharged particles is also discussed.

Dark Matter (DM) models providing possible alternative solutions to the small- scale crisis of standard cosmology are nowadays of growing interest. We consider DM interacting with light hidden fermions via well motivated fundamental operators showing the resultant matter power spectrum is suppressed on subgalactic scales within a plausible parameter region. Our basic description of the evolution of cosmological perturbations relies on a fully consistent first principles derivation of a perturbed Fokker-Planck type equation, generalizing existing literature. The cosmological perturbation of the Fokker-Planck equation is presented for the first time in two different gauges, where the results transform into each other according to the rules of gauge transformation. Furthermore, our focus lies on a derivation of a broadly applicable and easily computable collision term showing important phenomenological differences to other existing approximations. As one of the main results and concerning the small-scale crisis, we show the equal importance of vector and scalar boson mediated interactions between DM and light fermions.

The rotation curves of spiral galaxies exhibit a diversity that has been difficult to understand in the cold dark matter (CDM) paradigm. We show that the self-interacting dark matter (SIDM) model provides excellent fits to the rotation curves of a sample of galaxies with asymptotic velocities in the 25 to 300 km/s range that exemplify the full range of diversity. We only assume the halo concentration-mass relation predicted by the CDM model and a fixed value of the self-interaction cross section.In dark matter dominated galaxies, thermalization due to self-interactions creates large cores and reduces dark matter densities. In contrast, thermalization leads to denser and smaller cores in more luminous galaxies, and naturally explains the flat rotation curves of the highly luminous galaxies. Our results demonstrate that the impact of the baryons on the SIDM halo profile and the scatter from the assembly history of halos as encoded in the concentration-mass relation can explain the diverse rotation curves of spiral galaxies.

We provide a variant model of strongly interacting massive particle (SIMP), where composite dark matter comes from a strongly interacting U(1) theory. We first explain a non-Abelian version of the model with an additional singlet field, which is mixed with the Higgs field to maintain the kinetic equilibrium between the hidden and Standard Model (SM) sectors. The mixing leads to signals that would be detected by future collider experiments, direct DM detection experiments, and beam-dump experiments. Then we investigate a U(1) theory with a scalar monopole, where U(1) charged particles are confined by monopole condensation. In this model, the radial component of monopole can mix with the Higgs field, so that we do not need to introduce the additional singlet field.

Light gravitinos of mass $\lesssim \mathcal{O} (10)$ eV are of particular interest in cosmology, offering various baryogenesis scenarios without suffering from the cosmological gravitino problem. The gravitino may contribute considerably to the total matter content of the Universe and affect structure formation from early to present epochs. After the gravitinos decouple from other particles in the early Universe, they free-stream and consequently suppress density fluctuations of (sub-)galactic length scales. Observations of structure at the relevant length-scales can be used to infer or constrain the mass and the abundance of light gravitinos. We derive constraints on the light gravitino mass using the data of cosmic microwave background (CMB) lensing from Planck and of cosmic shear from the Canada France Hawaii Lensing Survey, combined with analyses of the primary CMB anisotropies and the signature of baryon acoustic oscillations in galaxy distributions. The obtained constraint on the gravitino mass is $m_{3/2} < 4.7$ eV (95% C.L.), which is substantially tighter than the previous constraint from clustering analysis of Ly-$\alpha$ forests.

Recently it has been claimed that the warm dark matter (WDM) model cannot at the same time reproduce the observed Lyman-{\alpha} forests in distant quasar spectra and solve the small-scale issues in the cold dark matter (CDM) model. As an alternative candidate, it was shown that the mixed dark matter (MDM) model that consists of WDM and CDM can satisfy the constraint from Lyman-{\alpha} forests and account for the "missing satellite problem" as well as the reported 3.5 keV anomalous X-ray line. We investigate observational constraints on the MDM model using strong gravitational lenses. We first develop a fitting formula for the nonlinear power spectra in the MDM model by performing N-body simulations and estimate the expected perturbations caused by line-of-sight structures in four quadruply lensed quasars that show anomaly in the flux ratios. Our analysis indicates that the MDM model is compatible with the observed anomaly if the mass fraction of the warm component is smaller than 0.47 at the 95% confidence level. The MDM explanation to the anomalous X-ray line and the small-scale issues is still viable even after this constraint is taken into account.

We perform a set of cosmological simulations of structure formation in a mixed dark matter (MDM) model. Our model is motivated by the recently identified $3.5\,{\rm keV}$ X-ray line, which can be explained by the decay of non-resonantly produced sterile neutrinos accounting for $20$-$60\%$ of the dark matter in the Universe. These non-resonantly produced sterile neutrinos have a sizable free-streaming length and hence behave effectively as warm dark matter (WDM). Assuming the rest of dark matter is composed of some cold dark matter (CDM) particles, we follow the coevolution of a mixed WDM plus CDM cosmology. Specifically, we consider the models with the warm component fraction of $r_{\rm warm}=0.25$ and $0.50$. Our MDM models predict that the comoving Jeans length at the matter-radiation equality is close to that of the thermally produced warm dark matter model with particle mass $m_{\rm WDM}=2.4\,{\rm keV}$, but the suppression in the fluctuation power spectrum is weaker. We perform large $N$-body simulations to study the structure of non-linear dark halos in the MDM models. The abundance of substructure is significantly reduced in the MDM models, and hence the so-called small-scale crisis is mitigated. The cumulative maximum circular velocity function (CVF) of at least one halo in the MDM models is in good agreement with the CVFs of the observed satellites in the Milky Way and the Andromeda galaxy. We argue that the MDM models open an interesting possibility to reconcile the reported $3.5\,{\rm keV}$ line and the internal structure of galaxies.

Supersymmetric theories, including the minimal supersymmetric standard model, usually contain many scalar fields whose potentials are absent in the exact supersymmetric limit and within the renormalizable level. Since their potentials are vulnerable to the finite energy density of the Universe through supergravity effects, these flat directions have nontrivial dynamics in the early Universe. Recently, we have pointed out that a flat direction may have a positive Hubble induced mass term during inflation whereas a negative one after inflation. In this case, the flat direction stays at the origin of the potential during inflation and then obtain a large vacuum expectation value after inflation. After that, when the Hubble parameter decreases down to the mass of the flat direction, it starts to oscillate around the origin of the potential. In this paper, we investigate the dynamics of the flat direction with and without higher dimensional superpotentials and show that topological defects, such as cosmic strings and domain walls, form at the end of inflation and disappear at the beginning of oscillation of the flat direction. We numerically calculate their gravitational signals and find that the observation of gravitational signals would give us information of supersymmetric scale, the reheating temperature of the Universe, and higher dimensional operators.

The energy spectrum of high-energy neutrinos reported by the IceCube collaboration shows a dip between 400 TeV and 1 PeV. One intriguing explanation is that high-energy neutrinos scatter with the cosmic neutrino background through a $\sim$ MeV mediator. Taking the density matrix approach, we develop a formalism to study the propagation of PeV neutrinos in the presence of the new neutrino interaction. If the interaction is flavored such as the gauged $L_\mu-L_\tau$ model we consider, the resonant collision may not suppress the PeV neutrino flux completely. The new force mediator may also contribute to the number of effectively massless degrees of freedom in the early universe and change the diffusion time of neutrinos from the supernova core. Astrophysical observations such as Big Bang Nucleosynthesis and supernova cooling provide an interesting test for the explanation.

The discovery of the primordial B-mode polarisation by the BICEP2 experiment indicates inflation with a relatively high energy scale. Taking this indication into account, we propose consistent scenarios to account for the observed baryon and dark matter densities in gravity and gauge mediated supersymmetry breaking models. The baryon asymmetry is explained by the Afflck-Dine mechanism, which requires relatively low reheating temperature to avoid a sizable baryonic isocurvature perturbation. The low reheating temperature then requires non-thermal production of dark matter to account for the correct relic density of dark matter. Our scenarios can account for the observations of baryon and dark matter density in gravity and gauge mediation and predict some parameters, including the mass of dark matter.

We investigate the sources of the Hubble-induced mass for a flat direction in supersymmetric theories and show that the sign of the Hubble-induced mass generally changes just after the end of inflation. This implies that global cosmic strings generally form after the end of inflation in most supersymmetric models, including the Minimal Supersymmetric Standard Model. The cosmic strings emit gravitational waves whose frequency corresponds to the Hubble scale, until they disappear when the Hubble parameter decreases down to the soft mass of the flat direction. As a result, the peak frequency of gravitational waves is related to the supersymmetric scale. The observation of this gravitational wave signal will give us information of supersymmetric scale and reheating temperature.

We show that the baryon-dark matter coincidence problem is solved in the CMSSM. The baryons and dark matter are generated simultaneously through the late-time decay of non-topological solitons, Q-balls, which are formed after the Affleck-Dine baryogenesis. A certain relation between the universal scalar mass, $m_0$, and the universal gaugino mass, $M_{1/2}$, is required to solve the coincidence problem, marginally depending on the other CMSSM parameters, and the result is consistent with the observation of the 126-GeV Higgs boson.

We explore the discovery potential of light gravitino mass m_{3/2} by combining future cosmology surveys and collider experiments. The former probe the imprint of light gravitinos in the cosmic matter density field, whereas the latter search signatures of a supersymmetry breaking mechanism. Free-streaming of light gravitinos suppresses the density fluctuations at galactic and sub-galactic length scales, where weak gravitational lensing can be used as a powerful probe. We perform numerical simulations of structure formation to quantify the effect. We then run realistic ray-tracing simulations of gravitational lensing to measure the cosmic shear in models with light gravitino. We forecast the possible reach of future wide-field surveys by Fisher analysis; the light gravitino mass can be determined with an accuracy of m_{3/2}=4\pm 1 eV by a combination of the Hyper Suprime Cam survey and cosmic microwave background anisotropy data obtained by Planck satellite. The corresponding accuracy to be obtained by the future Large Synoptic Survey Telescope is \delta m_{3/2}=0.6 eV. Data from experiments at Large Hadron Collider at 14 TeV will provide constraint at m_{3/2} \simeq 5 eV in the minimal framework of gauge-mediated supersymmetry breaking (GMSB) model. We conclude that a large class of the GMSB model can be tested by combining the cosmological observations and the collider experiments.

We study phenomenological aspects of the bino-wino co-annihilation scenario in high-scale supersymmetry breaking models. High-scale SUSY breaking scenarios are considered to be promising possibility after the discovery of the Higgs boson with a mass around 126 GeV. In this paper, we discuss the bino lightest supersymmetric particle (LSP) accompanied by the at most around 30 GeV heavier wino. With the suitable mass splitting between the bino and the wino, the bino LSP has the correct relic abundance of dark matter. For the smaller mass splitting, the late-time decay of the gravitino can provide the correct abundance of the bino dark matter. It is extremely challenging to find signals from the bino dark matter in direct and indirect detections. By utilizing multi-jets plus missing transverse momentum events at the LHC, we can constraint the gluino mass and thus probe the bino-wino co-annihilation scenario indirectly. The collider experiment, however, can not search the bino dark matter directly. In this paper, we suggest the direct probe of the bino dark matter. We show that the bino dark matter leaves imprints on the small-scale matter power spectrum when the bino dark matter is produced by the decay of the gravitino. The non-thermal bino dark matter behaves as mixed (cold+warm) dark matter.

The observed baryon and dark matter densities are equal up to a factor of 5. This observation indicates that the baryon asymmetry and dark matter have the same origin. The Affleck-Dine baryogenesis is one of the most promising mechanisms in this context. Q balls, which are often formed in the early Universe associated with the Affleck-Dine baryogenesis, decay both into supersymmetric particles and into quarks. Recently, it was pointed out that annihilation of squarks into quarks gives a dominant contribution to the Q-ball decay rate and the branching ratio of Q-ball decay into supersymmetric particles changes from the previous estimate. In this paper, the scenario of baryon and dark matter cogenesis from Q ball in gravity mediation is revisited in respect of the improved Q-ball decay rates. It is found that the successful cogenesis takes place when a wino with mass 0.3-1 TeV is dark matter.

We study the formation of non-linear structures in Warm Dark Matter (WDM) models and in a Long-Lived Charged Massive Particle (CHAMP) model. CHAMPs with a decay lifetime of about 1 yr induce characteristic suppression in the matter power spectrum at subgalactic scales through acoustic oscillations in the thermal background. We explore structure formation in such a model. We also study three WDM models, where the dark matter particles are produced through the following mechanisms: i) WDM particles are produced in the thermal background and then kinematically decoupled; ii) WDM particles are fermions produced by the decay of thermal heavy bosons; and iii) WDM particles are produced by the decay of non-relativistic heavy particles. We show that the linear matter power spectra for the three models are all characterised by the comoving Jeans scale at the matter-radiation equality. Furthermore, we can also describe the linear matter power spectrum for the Long-Lived CHAMP model in terms of a suitably defined characteristic cut-off scale k_{Ch}, similarly to the WDM models. We perform large cosmological N-body simulations to study the non-linear growth of structures in these four models. We compare the halo mass functions, the subhalo mass functions, and the radial distributions of subhalos in simulated Milky Way-size halos. We study the models with k_{cut}=51, 410, 820 h/Mpc, and confirm that these statistics are indeed similar between the different WDM models and the Long-Lived CHAMP model. The result suggests that the cut-off scale k_{cut} not only characterises the linear power spectra but also can be used to predict the non-linear clustering properties. The radial distribution of subhalos in Milky Way-size halos is consistent with the observed distribution for k_{cut}~50-800 h/Mpc; such models resolve the so-called "missing satellite problem".

We study how "warm" the wino dark matter is when it is non-thermally produced by the decays of the gravitino in the early Universe. We clarify the energy distribution of the wino at the decay of the gravitino and the energy loss process after their production. By solving the Boltzmann equation, we show that a sizable fraction of the wino dark matter can be "warm" for the wino mass m_{\tilde w} \sim 100-500 GeV. The "warmness" of the wino dark matter leaves imprints on the matter power spectra and may provide further insights on the origin of dark matter via the future 21 cm line survey. Our calculations can be applied to other non-thermal wino production scenarios such as the wino dark matter produced by the decay of the moduli fields.