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.

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 investigate models involving a vector-like quark X, which forms a 750 GeV bound state and reproduces the observed diphoton signals at the LHC, in connection with other excesses in the LHC data. An exotic hypercharge of -4/3 is required to fit the signal cross section, which indicates that there is additional particle(s) that mediates the decay of X in the full theory. We find that, introducing an SU(2) doublet vector-like quark of mass around 600 GeV in our UV-complete framework can accommodate not only the diphoton but also the on-Z excess (and potentially a slight excess in the monojet events). Our models also include a dark matter candidate. The most useful way to probe the models at the LHC is via monojet searches. The relic dark matter density is largely determined by coannihilation effects, and indirect detection of dark matter annihilation signals is the alternative and complementary probe of our models.

In light of positive signals reported by the CDMS-II Si experiment and the recent results of the LUX and SuperCDMS experiments, we study isospin-violating dark matter scenarios assuming that the interaction of the dark matter is mediated by colored particles. We investigate the phenomenology of the model, including collider searches, flavor and CP phenomenology. A minimal possible scenario includes scalar dark matter and new vector-like colored fermions with masses of O(1) TeV as mediators. Such a scenario may be probed at the 14 TeV LHC, while flavor and CP constraints are stringent and severe tuning in the couplings is unavoidable. We also found that, as an explanation of the CDMS-II Si signal, isospin-violating fermionic dark matter models with colored scalar mediators are disfavored by the LHC constraints.

Axino as the superpartner of axion that solves the strong CP problem can be a good candidate of dark matter. Inspired by the 3.5 keV X-ray line signal found to be originated from galaxy clusters and Andromeda galaxy, we study axino models with R-parity violations, and point out that axino dark matter with trilinear R-parity violations is an attractive scenario that reproduces the X-ray line. The Peccei-Quinn scale is required to be $f_a \sim{\cal O}(10^{9}-10^{11})$ GeV for trilinear R-parity violating couplings $\lambda \sim {\cal O} (10^{-3}-10^{-1})$ in order to explain the line signal. Moreover, the right-handed stau is predicted to be light, i.e. $\sim{\cal O}(100)$ GeV, and thus can be looked for at the LHC. Cosmological aspects of the model are also discussed in this study.

We study radiative decay of gravitino dark matter with trilinear R-parity violations. We show that the branching ratio of the decay of gravitino into monochromatic photon can be large enough to explain the observed gamma-ray line from the Galactic centre in the Fermi-LAT data without producing too much continuum gamma-ray and anti-proton flux. This scenario is realized when the mass of sfermions and the trilinear R-parity violating coupling are $O(1-10)$ TeV and $O(10^{-7}-10^{-6})$ respectively.

We show that decaying axino dark matter with R-parity violation can explain the observed excess of the 130GeV gamma-ray line from the Galactic center in the Fermi data. The branching fraction of the axino decay into monochromatic photons can be O(1), and constraints from continuum gamma-rays and the anti-proton flux are ameliorated. The Peccei-Quinn scale of $O(10^{13}-10^{14})$\,GeV and the R-parity violation parameter of $O(10^{-12}-10^{-11})$ are cosmologically favored.