Electromagnetically neutral dark sector particles may directly couple to the photon through higher dimensional effective operators. Considering electric and magnetic dipole moment, anapole moment, and charge radius interactions, we derive constraints from stellar energy loss in the Sun, horizontal branch and red giant stars, as well as from cooling of the proto-neutron star of SN1987A. We provide the exact formula for in-medium photon-mediated pair production to leading order in the dark coupling, and compute the energy loss rates explicitly for the most important processes, including a careful discussion on resonances and potential double counting between the processes. Stringent limits for dark states with masses below $3\,$keV ($40\,$MeV) arise from red giant stars (SN1987A), implying an effective lower mass-scale of approximately $10^9\,$GeV ($10^7\,$GeV) for mass-dimension five, and $100\,$GeV ($2.5\,$TeV) for mass-dimension six operators as long as dark states stream freely; for the proto-neutron star, the trapping of dark states is also evaluated. Together with direct limits previously derived by us in Chu et al. (2018), this provides the first comprehensive overview of the viability of effective electromagnetic dark-state interactions below the GeV mass-scale.

We examine the cosmology of warm dark matter (WDM), both stable and decaying, from the point of view of structure formation. We compare the matter power spectrum associated to WDM masses of 1.5 keV and 0.158 keV, with that expected for the stable cold dark matter $\Lambda$CDM$\equiv$SCDM paradigm, taken as our reference model. We scrutinize the effects associated to the warm nature of dark matter, as well as the fact that it decays. The decaying warm dark matter (DWDM) scenario is well-motivated, emerging in a broad class of particle physics theories where neutrino masses arise from the spontaneous breaking of a continuous global lepton number symmetry. The majoron arises as a Nambu-Goldstone boson, and picks up a mass from gravitational effects, that explicitly violate global symmetries. The majoron necessarily decays to neutrinos, with an amplitude proportional to their tiny mass, which typically gives it cosmologically long lifetimes. Using N-body simulations we show that our DWDM picture leads to a viable alternative to the $\Lambda$CDM scenario, with predictions that can differ substantially on small scales.

The recently announced results on the 21-cm absorption spectrum by the EDGES experiment can place very stringent limits on dark matter annihilation cross-sections. We properly take into account the heating energy released from dark matter annihilation from the radiation epoch to the 21-cm observation redshifts in the radiative transfer to compute the evolution of the gas temperature. Our results show that the global 21-cm absorption profile is a powerful cosmological probe of the dark matter interactions. For dark matter annihilating into electron-positron pairs, the EDGES results give a more stringent upper limit than the PLANCK result on the annihilation cross section at the lower dark matter mass region.

With recent Lyman-alpha forest data from BOSS and XQ-100, some studies suggested that the lower mass limit on the fuzzy dark matter (FDM) particles is lifted up to $10^{-21}\,\mathrm{eV}$. However, such a limit was obtained by $\Lambda$CDM simulations with the FDM initial condition and the quantum pressure (QP) was not taken into account which could have generated non-trivial effects in large scales structures. We investigate the QP effects in cosmological simulations systematically, and find that the QP leads to further suppression of the matter power spectrum at small scales, as well as the halo mass function in the low mass end. We estimate the suppressing effect of QP in the 1D flux power spectrum of Lyman-alpha forest and compare it with data from BOSS and XQ-100. The rough uncertainties of thermal gas properties in the flux power spectrum model calculation were discussed. We conclude that more systematic studies, especially with QP taken into account, are necessary to constrain FDM particle mass using Lyman-alpha forest.

The Ultra-Light Axion (ULA) is a dark matter candidate with mass $\mathcal{O}(10^{-22})$eV and de Broglie wavelength of order kpc. Such an axion, also called the Fuzzy Dark Matter (FDM), thermalizes via the gravitational force and forms a Bose-Einstein condensate. Recent studies suggested that the quantum pressure from the FDM can significantly affect the structure formation in small scales, thus alleviating the so-called "small-scale crisis". In this paper, we develop a new technique to discretize the quantum pressure and illustrate the interactions among FDM particles in the $N$-body simulation, which accurately simulates the formation of the dark-matter halo and its inner structure in the region outside the softening length. In a self-gravitationally-bound virialized halo, we find a constant density, solitonic core, which is consistent with the theoretical prediction. The existence of the solitonic core reveals the non-linear effect of quantum pressure and impacts the structure formation in the FDM model.