We point out a unique mechanism to produce the relic abundance for glueball dark matter from a gauged $SU(N)_d$ hidden sector which is bridged to the standard model sector through heavy vectorlike quarks colored under gauge interactions from both sides. A necessary ingredient of our assumption is that the vectorlike quarks, produced either thermally or non-thermally, are abundant enough to dominate the universe for some time in the early universe. They later undergo dark color confinement and form unstable vectorlike-quarkonium states which annihilate decay and reheat the visible and dark sectors. The ratio of entropy dumped into two sectors and the final energy budget in the dark glueballs is only determined by low energy parameters, including the intrinsic scale of the dark $SU(N)_d$, $\Lambda_d$, and number of dark colors, $N_d$, but depend weakly on parameters in the ultraviolet such as the vectorlike quark mass or the initial condition. We call this a cosmic selection rule for the glueball dark matter relic density.

A hidden sector with pure non-abelian gauge symmetry is an elegant and just about the simplest model of dark matter. In this model the dark matter candidate is the lightest bound state made of the confined gauge fields, the dark glueball. In spite of its simplicity, the model has been shown to have several interesting non-standard implications in cosmology. In this work, we explore the gravitational waves from binary boson stars made of self-gravitating dark glueball fields as a natural and important consequence. We derive the dark SU($N$) star mass and radius as functions of the only two fundamental parameters in the model, the glueball mass $m$ and the number of colors $N$, and identify the regions that could be probed by the LIGO and future gravitational wave observatories.

We discuss a minimal solution to the long-standing $(g-2)_\mu$ anomaly in a simple extension of the Standard Model with an extra $Z'$ vector boson that has only flavor off-diagonal couplings to the second and third generation of leptons, i.e. $\mu, \tau, \nu_\mu, \nu_\tau$ and their antiparticles. A simplified model realization, as well as various collider and low-energy constraints on this model, are discussed. We find that the $(g-2)_\mu$-favored region for a $Z'$ lighter than the tau lepton is totally excluded, while a heavier $Z'$ solution is still allowed. Some testable implications of this scenario in future experiments, such as lepton-flavor universality-violating tau decays at Belle 2, and a new four-lepton signature involving same-sign di-muons and di-taus at HL-LHC and FCC-ee, are pointed out. A characteristic resonant absorption feature in the high-energy neutrino spectrum might also be observed by neutrino telescopes like IceCube and KM3NeT.

We investigate the possibility that the dark matter candidate is from a pure non-abelian gauge theory of the hidden sector, motivated in large part by its elegance and simplicity. The dark matter is the lightest bound state made of the confined gauge fields, the hidden glueball. We point out this simple setup is capable of providing rich and novel phenomena in the dark sector, especially in the parameter space of large N. They include self-interacting and warm dark matter scenarios, Bose-Einstein condensation leading to massive dark stars possibly millions of times heavier than our sun giving rise to gravitational lensing effects, and indirect detections through higher dimensional operators as well as interesting collider signatures.

Understanding the spectral and flavor composition of the astrophysical neutrino flux responsible for the recently observed ultra-high energy events at IceCube is of great importance for both astrophysics and particle physics. We perform a statistical likelihood analysis to the 3-year IceCube data and derive the allowed range of the spectral index and flux normalization for various well-motivated physical flavor compositions at source. While most of the existing analyses so far assume the flavor composition of the neutrinos at an astrophysical source to be (1:2:0), it seems rather unnatural to assume only one type of source, once we recognize the possibility of at least two physical sources. Bearing this in mind, we entertain the possibility of a two-component source for the analysis of IceCube data. It appears that our two component hypothesis explains some key features of the data better than a single-component scenario, i.e it addresses the apparent energy gap between 400 TeV to about 1 PeV and easily accommodates the observed track to shower ratio. Given the extreme importance of the flavor composition for the correct interpretation of the underlying astrophysical processes as well as for the ramification for particle physics, this two-component flux should be tested as more data is accumulated.

The recent observation of two PeV events at IceCube, followed by an additional 26 events between 30 - 300 TeV, has generated considerable speculations on its origin, and many exotic New Physics explanations have been invoked. For a reliable interpretation, it is however important to first scrutinize the Standard Model (SM) expectations carefully, including the theoretical uncertainties, mainly due to the parton distribution functions. Assuming a new isotropic cosmic neutrino flux with a simple unbroken power-law spectrum, $\Phi\propto E^{-s}$ for the entire energy range of interest, we find that with $s=1.5$ - 2, the SM neutrino-nucleon interactions are sufficient to explain all the observed events so far, without the need for any beyond the SM explanation. With more statistics, this powerful detector could provide a unique test of the SM up to the PeV scale, and lead to important clues of New Physics.

The possibility of direct detection of light fermionic dark matter in neutrino detectors is explored from a model-independent standpoint. We consider all operators of dimension six or lower which can contribute to the interaction $\bar{f} p \to e^+ n$, where $f$ is a dark Majorana or Dirac fermion. Constraints on these operators are then obtained from the $f$ lifetime and its decays which produce visible $\gamma$ rays or electrons. We find one operator which would allow $\bar{f} p \to e^+ n$ at interesting rates in neutrino detectors, as long as $m_f \lesssim m_{\pi}$. The existing constraints on light dark matter from relic density arguments, supernova cooling rates, and big-bang nucleosynthesis are then reviewed. We calculate the cross-section for $\bar{f} p \to e^+ n$ in neutrino detectors implied by this operator, and find that Super-K can probe the new physics scale $\Lambda$ for this interaction up to ${\cal O}(100 {TeV})$