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.

Implications of supersymmetrizing the clockwork axions are stuided. Supersymmetry ensures that the saxions and axinos have the same pattern of the coupling hierarchy as the clockwork axions. If we assume supersymmetry breaking is universal over the the clockwork sites, the coupling structure is preserved, while the mass orderings of the saxions and axinos can differ depending on the supersymmetry breaking scale. While the massive saxions and axions quickly decay, the lightest axino can be stable and thus a dark matter candidate. The relic abundance of the axino dark matter from thermal production is mostly determined by decays of the heavier axinos in the normal mass ordering. This exponentially enhances the thermal yield compared to the conventional axino scenarios. Some cosmological issues are discussed.

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 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 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.

Under the expectation that nature is natural, we extend the Standard Model to include SUSY to stabilize the electroweak sector and PQ symmetry to stabilize the QCD sector. Then natural SUSY arises from a Kim-Nilles solution to the SUSY mu problem which allows for a little hierarchy where mu ~ f_a^2/M_P ~ 100-300 GeV while the SUSY particle mass scale m(SUSY)~ 1-10 TeV >> mu. Dark matter then consists of two particles: a higgsino-like WIMP and a SUSY DFSZ axion. The range of allowed axion mass values m(axion) depends on the mixed axion-higgsino relic density. The range of m(axion) is actually restricted in this case by limits on WIMPs from direct and indirect detection experiments. We plot the expected axion detection rate at microwave cavity experiments. The axion-photon-photon coupling is severely diminished by charged higgsino contributions to the anomalous coupling. In this case, the axion may retreat, at least temporarily, back into the regime of near invisibility. From our results, we urge new ideas for techniques which probe both deeper and more broadly into axion coupling versus axion mass parameter space.

Supersymmetric models with radiatively-driven electroweak naturalness require light higgsinos of mass ~ 100-300 GeV. Naturalness in the QCD sector is invoked via the Peccei-Quinn (PQ) axion leading to mixed axion-higgsino dark matter. The SUSY DFSZ axion model provides a solution to the SUSY mu problem and the Little Hierarchy \mu << m_{3/2} may emerge as a consequence of a mismatch between PQ and hidden sector mass scales. The traditional gravitino problem is now augmented by the axino and saxion problems, since these latter particles can also contribute to overproduction of WIMPs or dark radiation, or violation of BBN constraints. We compute regions of the T_R vs. m_{3/2} plane allowed by BBN, dark matter and dark radiation constraints for various PQ scale choices f_a. These regions are compared to the values needed for thermal leptogenesis, non-thermal leptogenesis, oscillating sneutrino leptogenesis and Affleck-Dine leptogenesis. The latter three are allowed in wide regions of parameter space for PQ scale f_a~ 10^{10}-10^{12} GeV which is also favored by naturalness: f_a~ \sqrt{\mu M_P/\lambda_\mu }\sim 10^{10}-10^{12} GeV. These f_a values correspond to axion masses somewhat above the projected ADMX search regions.

By insisting on naturalness in both the electroweak and QCD sectors of the MSSM, the portrait for dark matter production is seriously modified from the usual WIMP miracle picture. In SUSY models with radiatively-driven naturalness (radiative natural SUSY or RNS) which include a DFSZ-like solution to the strong CP and SUSY mu problems, dark matter is expected to be an admixture of both axions and higgsino-like WIMPs. The WIMP/axion abundance calculation requires simultaneous solution of a set of coupled Boltzmann equations which describe quasi-stable axinos and saxions. In most of parameter space, axions make up the dominant contribution of dark matter although regions of WIMP dominance also occur. We show the allowed range of PQ scale f_a and compare to the values expected to be probed by the ADMX axion detector in the near future. We also show WIMP detection rates which are suppressed from usual expectations because now WIMPs comprise only a fraction of the total dark matter. Nonetheless, ton-scale noble liquid detectors should be able to probe the entirety of RNS parameter space. Indirect WIMP detection rates are less propitious since they are reduced by the square of the depleted WIMP abundance.

We examine mixed axion/neutralino cold dark matter production in the SUSY DFSZ axion model where an axion superfield couples to Higgs superfields. We calculate a wide array of axino and saxion decay modes along with their decay temperatures, and thermal and non-thermal production rates. For a SUSY benchmark model with a standard underabundance (SUA) of Higgsino-like dark matter (DM), we find for the PQ scale f_a< 10^{12} GeV that the DM abundance is mainly comprised of axions as the saxion/axino decay occurs before the standard neutralino freeze-out and thus its abundance remains suppressed. For 10^{12}<~ f_a<~ 10^{14} GeV, the saxion/axino decays occur after neutralino freeze-out so that the neutralino abundance is enhanced by the production via decay and subsequent re-annihilation. For f_a>~ 10^{14} GeV, both neutralino dark matter and dark radiation are typically overproduced. For judicious parameter choices, these can be suppressed and the combined neutralino/axion abundance brought into accord with measured values. A SUSY benchmark model with a standard overabundance (SOA) of bino DM is also examined and typically remains excluded due at least to too great a neutralino DM abundance for f_a<~ 10^{15} GeV. For f_a>~ 10^{15} GeV and lower saxion masses, large entropy production from saxion decay can dilute all relics and the SOA model can be allowed by all constraints.

By eschewing finetuning from the electroweak and QCD sectors of supersymmetry (natural supersymmetry or SUSY), and by invoking the Kim-Nilles solution to the SUSY mu problem, one is lead to models wherein the dark matter is comprised of a mixture of axions and higgsino-like WIMPs. Over a large range of Peccei-Quinn breaking scale f_a~ 10^9-10^{12} GeV, one then expects about 90-95% axion dark matter. In such a scenario, both axion and WIMP direct detection may be expected.

The supersymmetrized DFSZ axion model is highly motivated not only because it offers solutions to both the gauge hierarchy and strong CP problems, but also because it provides a solution to the SUSY mu problem which naturally allows for a Little Hierarchy. We compute the expected mixed axion-neutralino dark matter abundance for the SUSY DFSZ axion model in two benchmark cases-- a natural SUSY model with a standard neutralino underabundance (SUA) and an mSUGRA/CMSSM model with a standard overabundance (SOA). Our computation implements coupled Boltzmann equations which track the radiation density along with neutralino, axion (produced thermally (TH) and via coherent oscillations (CO)), saxion (TH- and CO-produced), axino and gravitino densities. In the SUSY DFSZ model, axions, axinos and saxions go through the process of freeze-in-- in contrast to freeze-out or out-of-equilibrium production as in the SUSY KSVZ model-- resulting in thermal yields which are largely independent of the re-heat temperature. We find the SUA case with suppressed saxion-axion couplings (\xi=0) only admits solutions for PQ breaking scale f_a~< 5\times 10^{12} GeV where the bulk of parameter space tends to be axion-dominated. For SUA with allowed saxion-axion couplings (\xi =1), then f_a values up to ~ 2\times 10^{14} GeV are allowed. For the SOA case, almost all of SUSY DFSZ parameter space is disallowed by a combination of overproduction of dark matter, overproduction of dark radiation or violation of BBN constraints. An exception occurs at very large f_a~ 10^{15}-10^{16} GeV where large entropy dilution from CO-produced saxions leads to allowed models.

A variety of supersymmetric models give rise to a split mass spectrum characterized by very heavy scalars but sub-TeV gauginos, usually with a wino-like LSP. Such models predict a thermally-produced underabundance of wino-like WIMP dark matter so that non-thermal DM production mechanisms are necessary. We examine the case where theories with a wino-like LSP are augmented by a Peccei-Quinn sector including an axion-axino-saxion supermultiplet in either the SUSY KSVZ or SUSY DFSZ models and with/without saxion decays to axions/axinos. We show allowed ranges of PQ breaking scale f_a for various cases which are generated by solving the necessary coupled Boltzmann equations. We also present results for a model with radiatively-driven naturalness but with a wino-like LSP.

If one simultaneously invokes the SUSY solution to the gauge hierarchy problem and the PQ solution to the strong CP problem, then one might expect mixed axion/neutralino dark matter (DM), i.e. two dark matter particles. In this case, dark matter production involves oscillation production of axions, thermal and non-thermal neutralino production and thermal and non-thermal production of axinos and saxions. The latter particles may feed additional neutralinos into the cosmic soup, or dilute all relics via additional entropy production, or feed dark radiation into the cosmos. Depending on the SUSY and PQ model, the PQ scale can easily exist from 10^9-10^{16} GeV. Axions may lie well beyond the range currently being probed by experiment. In the case of relic WIMPs, theoretical predictions must be rescaled for the possibility that WIMPs only comprise a fraction of the local DM abundance.

Recent analyses of WMAP9 data show that dark radiation-- parametrized by the apparent number of additional neutrinos \Delta N_{eff} contributing to the cosmic expansion-- is consistent with the Standard Model and bounded from above by about \Delta N_{eff} ~< 0.5 at 95% CL. We consider the mixed axion/neutralino cold dark matter scenario which arises in R-parity conserving supersymmetric (SUSY) models wherein the strong CP problem is solved by hadronic axions with a concommitant axion(a)/saxion(s)/axino(\ta) supermultiplet. Our new results include improved calculations of thermal axion and saxion production and include effects of saxion decay to axinos and axions. We show that the above bound on \Delta N_{eff} is easily satisfied if saxions are mainly thermally produced and m_{LSP} < m_{\ta} \lesssim m_s. However, if the dominant mechanism of saxion production is through coherent oscillations, the WMAP9 data provides a strong bound on saxion production followed by saxion decays to axions. Furthermore we show that scenarios with mixed neutralino/axion dark matter are highly constrained by combined WMAP9, BBN and Xe-100 constraints. In particular, supersymmetric models with a standard overabundance of neutralino dark matter are excluded for all values of the Peccei-Quinn breaking scale. Next generation WIMP direct detection experiments may be able to discover or exclude mixed axion-neutralino CDM scenarios where $s\to aa$ is the dominant saxion decay mode.

We study the cosmological impact of the supersymmetric DFSZ axion model. Extending recent works, we first provide a comprehensive analysis of thermal production of the DFSZ axino considering all the possible scattering, decay and inverse decay processes depending on various mass parameters and the reheat temperature. Although it is hard for the DFSZ axino to be in thermal equilibrium, its coupling is still large enough to generate huge axino population which can turn into overabundant neutralino density. We examine the neutralino parameter space to identify the dark matter property depending on the Peccei-Quinn scale. As the Peccei-Quinn scale becomes higher resulting in longer axino lifetime, the neutralino dark matter appears in a lighter Higgsino-like LSP region or a more restricted Bino-like LSP region allowing a resonant annihilation through a CP-odd Higgs boson to meet stronger reannihilation.

We discuss the effective interactions of axion supermultiplet, which might be important for analyzing the cosmological aspect of supersymmetric axion model. Related to axino cosmology, it is stressed that three seemingly similar but basically different quantities, the Wilsonian axino-gluino-gluon coupling, the 1PI axino-gluino-gluon amplitude, and the PQ anomaly coefficient, should be carefully distinguished from each other for correct analysis of the thermal production of axinos in the early Universe. It is then noticed that the 1PI axino-gluino-gluon amplitude at energy scale p in the range M_\Phi < p < v_{PQ} is suppressed by (M_\Phi ^2)/(p^2) in addition to the well-known suppression by p/(16\pi^2 v_{PQ}), where M_\Phi is the mass of the heaviest PQ-charged and gauge-charged matter supermultiplet in the model, which can be well below the PQ scale v_{PQ}. As a result, axino production at temperature T>M_\Phi is dominated by the production by matter supermultiplet, not by the production by gauge supermultiplet. Still the axino production rate is greatly reduced if M_\Phi << v_{PQ}, which would make the subsequent cosmology significantly altered. This would be most notable in the supersymmetric DFSZ model in which M_\Phi corresponds to the Higgsino mass which is around the weak scale, however a similar reduction is possible in the KSVZ model also. We evaluate the relic axino density for both the DFSZ and KSVZ models while including the axino production in the processes involving the heaviest PQ-charged and gauge-charged matter supermultiplet.

We discuss dark matter decay mediated by exotically charged particles ("exotics") in a supersymmetric model with two dark matter (DM) components: One is the (bino-like) lightest supersymmetric particle (LSP) \chi, and the other is a newly introduced meta-stable neutral singlet $N$. $N$ decays to \chi e^+e^- via a dimension 6 operator induced by a penguin-type one loop diagram with the life time of 10^{26} sec., explaining energetic cosmic e^\pm excess observed recently by PAMELA and ATIC/PPB-BETS. The superheavy masses of exotics (\sim 10^{15-16} GeV) are responsible for the longevity of $N$. The superpartner of $N$ develops the vacuum expectation value (VEV) of order TeV so that the DM $N$ achieves the desired mass of 2 TeV. By the VEV, the U(1)_R symmetry is broken to the discrete Z_2 symmetry, which is identified with the matter parity in the minimal supersymmetric standard model (MSSM). Since we have the two DM components, even extremely small amount of $N$ [O(10^{-10}) < (n_N/n_\chi)] could account for the observed positron flux with relatively light exotics' masses [10^{12} GeV < M_{exo.} < 10^{16} GeV].

Recently, there are two hints arising from physics beyond the standard model. One is a possible energy loss mechanism due to emission of very weakly interacting light particles from white dwarf stars, with a coupling strength ~ 0.7x10^{-13}, and another is the high energy positrons observed by the PAMELA satellite experiment. We construct a supersymmetric flipped-SU(5) model, SU(5)xU(1)_X with appropriate additional symmetries, [U(1)_H]_{gauge}x[U(1)_RxU(1)_\Gamma]_{global}xZ_2, such that these are explained by a very light electrophilic axion of mass 0.5 meV from the spontaneously broken U(1)_\Gamma and two component cold dark matters from Z_2 parity. We show that in the flipped-SU(5) there exists a basic mechanism for allowing excess positrons through the charged SU(2) singlet leptons, but not allowing anti-proton excess due to the absence of the SU(2) singlet quarks. We show the discovery potential of the charged SU(2) singlet E at the LHC experiments by observing the electron and positron spectrum. With these symmetries, we also comment on the mass hierarchy between the top and bottom quarks.

We improve the estimate of the axion CDM energy density by considering the new values of current quark masses, the QCD phase transition effect and a possible anharmonic effect.