Measuring high energy cosmic ray electrons/positrons (CRE) provides important means for probing the nearby galactic sources and for the dark matter (DM) indirect detection. In this work, we perform a systematical analysis of the flavor structure of DM annihilations into charged leptons based on the cosmic ray CRE spectra measured by DAMPE, Fermi-LAT, AMS-02, and CALET experiments. We study the annihilations of possible TeV scale DM particles in a nearby subhalo, which is proposed to explain the possible peak-like structure of the DAMPE CRE data. We pay special attention to the possible non-resonant excess (besides the possible peak-like structure) and demonstrate that such non-resonant excess can mainly arise from the decay of muons produced by the DM annihilations in the subhalo. With these we study the flavor composition of the lepton final states of $\chi\chi \to e^+e^-, \mu^+\mu^-, \tau^+\tau^-$ by fitting the CRE data. We further analyze the constraints on the lepton flavor composition using the Fermi-LAT $\gamma$-ray measurements. We find that the flavor composition (favored by the CRE data) is consistent with the Fermi-LAT data at relatively low Galactic latitudes, while the fraction of the final state $\tau^\pm$ is severely bounded.

The XMM-Newton observatory shows evidence with an $11 \sigma$ confidence level for seasonal variation of the X-ray background in the near-Earth environment in the 2-6 keV energy range (Fraser et al. 2014). The interpretation of the seasonal variation given in Fraser et al. 2014 was based on the assumption that solar axions convert to X-rays in the Earth's magnetic field. There are many problems with this interpretation, since the axion-photon conversion must preserve the directionality of the incoming solar axion. At the same time, this direction is avoided by the observations because the XMM-Newton's operations exclude pointing at the Sun and at the Earth. The observed seasonal variation suggests that the signal could have a dark matter origin, since it is very difficult to explain with conventional astrophysical sources. We propose an alternative explanation which involves the so-called Axion Quark Nugget (AQN) dark matter model. In our proposal, dark matter is made of AQNs, which can cross the Earth and emit high energy photons at their exit. We show that the emitted intensity and spectrum is consistent with Fraser et al. 2014, and that our calculation is not sensitive to the specific details of the model. We also find that our proposal predicts a large seasonal variation, on the level of 20-25%, much larger than conventional dark matter models (1-10%). Since the AQN emission spectrum extends up to $\sim$100 keV, well beyond the keV sensitivity of XMM-Newton, we predict the AQN contribution to the hard X-ray and $\gamma$-ray backgrounds in the Earth's environment. The Gamma-Ray Burst Monitor instrument (GBM), aboard the Fermi telescope, is sensitive to the 8 keV-40 MeV energy band. We suggest that the multi-year archival data from the GBM could be used to search for a seasonal variation in the near-Earth environment up to 100 keV as a future test of the AQN framework.

We study the formation of primordial black holes (PBHs) from the collapse of closed domain walls (DWs) which naturally arise in QCD axion models near the QCD scale together with the main string-wall network. The size distribution of the closed DWs is determined by percolation theory, from which we further obtain PBH mass distribution and abundance. Various observational constraints on PBH abundance in turn also constrain axion parameters. Our model prefers axion mass around the meV scale ($f_{a}\sim 10^{9}$ GeV). The corresponding PBHs are in the sublunar-mass window $10^{20}$-$10^{22}$ g (i.e., $10^{-13}$-$10^{-11}M_{\odot}$), one of few mass windows still available for PBHs contributing significantly to dark matter (DM). In our model, PBH abundance could reach $\sim1\%$ or even more of DM, sensitive to the formation efficiency of closed axion DWs.

We consider the formation and evolution of Axion Quark Nugget dark matter particles in the early universe. The goal of this work is to estimate the mass distribution of these objects and assess their ability to form and survive to the present day. We argue that this model allows a broad range of parameter space in which the AQN may account for the observed dark matter mass density, naturally explains a similarity between the "dark" and "visible" components, i.e. $\Omega_{\rm dark}\sim \Omega_{\rm visible}$, and also offer an explanation for a number of other long standing puzzles such as "Primordial Lithium Puzzle" and "the Solar Corona Mystery" among many other cosmological puzzles.

We propose to use atmospheric neutrinos as a powerful probe of new physics beyond the Standard Model via neutrino trident production. The final state with double muon tracks simultaneously produced from the same vertex is a distinctive signal at large Cherenkov detectors. We calculate the expected event numbers of trident production in the Standard Model. To illustrate the potential of this process to probe new physics we obtain the sensitivity on new vector/scalar bosons with coupling to muon and tau neutrinos.

The Dark Matter Particle Explorer (DAMPE) satellite detector announced its first result for measuring the cosmic-ray electron/positron (CRE) energy spectrum up to 4.6TeV, including a tentative peak-like event excess at (1.3-1.5)TeV. In this work, we uncover a significant hidden excess in the DAMPE CRE spectrum over the energy range (0.6-1.1)TeV, which has a non-peak-like structure. We propose a new mechanism to explain this excess by a set of 1.5TeV $\mu^\pm$ events with subsequent decays into $e^\pm$ plus neutrinos. For explaining this new excess together with the peak excess around 1.4TeV, we demonstrate that the {\it flavor structure} of the original lepton final-state produced by dark matter (DM) annihilations (or other mechanism) should have a composition ratio $N_e : (N_\mu +\frac{1}{6}N_\tau) = 1 : y$, with $y \simeq 2.6-10.8$. For lepton portal DM models, this puts important constraint on the lepton-DM-mediator couplings $\lambda_e : (\lambda_\mu^4 + \frac{1}{6}\lambda_\tau^4)^{\frac{1}{4}} = 1 : y^{\frac{1}{4}}$ with a narrow range $y^{\frac{1}{4}} \simeq 1.3-1.8$.

We study a dark matter (DM) model offering a very natural explanation of two (naively unrelated) problems in cosmology: the observed relation $\Omega_{\rm DM}\sim\Omega_{\rm visible}$ and the observed asymmetry between matter and antimatter in the Universe, known as the "baryogenesis" problem. In this framework, both types of matter (dark and visible) have the same QCD origin, form at the same QCD epoch, and both proportional to one and the same dimensional parameter of the system, $\Lambda_{\rm QCD}$, which explains how these two, naively distinct, problems could be intimately related, and could be solved simultaneously within the same framework. More specifically, the DM in this model is composed by two different ingredients: the (well- studied) DM axions and (less-studied) the quark nuggets made of matter or antimatter. The focus of the present work is the quantitative analysis of the relation between these two distinct components contributing to the dark sector of the theory determined by $\Omega_{\rm DM}\equiv [\Omega_{\rm DM}(\rm nuggets)+ \Omega_{\rm DM}(\rm axion)]$. We argue that the nugget's DM component always traces the visible matter density, i.e. $\Omega_{\rm DM}(\rm nuggets)\sim\Omega_{\rm visible}$ and this feature is not sensitive to the parameters of the system such as the axion mass $m_a$ or the misalignment angle $\theta_0$. It should be contrasted with conventional axion production mechanism due to the misalignment when $\Omega_{\rm DM}(\rm axion)$ is highly sensitive to the axion mass $m_a$ and the initial misalignment angle $\theta_0$. We also discuss the constraints on this model related to the inflationary scale $H_I$, non-observation of the isocurvature perturbations, $r_T < 0.12$, and also, varies axions search experiments.

We consider a dark matter (DM) model offering a very natural explanation of the observed relation, $\Omega_{\rm dark} \sim \Omega_{\rm visible}$. This generic consequence of the model is a result of the common origin of both types of matter (DM and visible) which are formed during the same QCD transition. The masses of both types of matter in this framework are proportional to one and the same dimensional parameter of the system, $\Lambda_{\rm QCD}$. The focus of the present work is the detail study of the dynamics of the $\cal{CP}$-odd coherent axion field $a(x)$ just before the QCD transition. We argue that the baryon charge separation effect on the largest possible scales inevitably occurs as a result of merely existence of the coherent axion field in early Universe. It leads to preferential formation of one species of nuggets on the scales of the visible Universe where the axion field $a(x)$ is coherent. A natural outcome of this preferential evolution is that only one type of the visible baryons remains in the system after the nuggets complete their formation. This represents a specific mechanism on how the baryon charge separation mechanism (when the Universe is neutral, but visible part of matter consists the baryons only) replaces the conventional "baryogenesis" scenarios.

Neutrino oscillation data strongly support mu-tau symmetry as a good approximate flavor symmetry of the neutrino sector, which has to appear in any viable theory for neutrino mass-generation. The mu-tau breaking is not only small, but also the source of Dirac CP-violation. We conjecture that both discrete mu-tau and CP symmetries are fundamental symmetries of the seesaw Lagrangian (respected by interaction terms), and they are only softly broken, arising from a common origin via a unique dimension-3 Majorana mass-term of the heavy right-handed neutrinos. From this conceptually attractive and simple construction, we can predict the soft mu-tau breaking at low energies, leading to quantitative correlations between the apparently two small deviations \theta_{23} - 45^o and \theta_{13} - 0^o. This nontrivially connects the on-going measurements of mixing angle \theta_{23} with the upcoming experimental probes of \theta_{13}. We find that any deviation of \theta_{23} - 45^o must put a lower limit on \theta_{13}. Furthermore, we deduce the low energy Dirac and Majorana CP violations from a common soft-breaking phase associated with mu-tau breaking in the neutrino seesaw. Finally, from the soft CP breaking in neutrino seesaw we derive the cosmological CP violation for the baryon asymmetry via leptogenesis. We fully reconstruct the leptogenesis CP-asymmetry from the low energy Dirac CP phase and establish a direct link between the cosmological CP-violation and the low energy Jarlskog invariant. We predict new lower and upper bounds on the \theta_{13} mixing angle, 1^o < \theta_{13} < 6^o. In addition, we reveal a new hidden symmetry that dictates the solar mixing angle \theta_12 by its group-parameter, and includes the conventional tri-bimaximal mixing as a special case, allowing deviations from it.