$SU(2)_L \times SU(2)_R$ gauge symmetry requires three right-handed neutrinos ($ N _i $), one of which, $N_1$, can be sufficiently stable to be dark matter. In the early universe, $ W _R $ exchange with the Standard Model thermal bath keeps the right-handed neutrinos in thermal equilibrium at high temperatures. $N_1$ can make up all of dark matter if they freeze-out while relativistic and are mildly diluted by subsequent decays of a long-lived and heavier right-handed neutrino, $N_2$. We systematically study this parameter space, constraining the symmetry breaking scale of $SU(2)_R$ and the mass of $N_1$ to a triangle in the $(v_R,M_1)$ plane, with $v_R = (10^6 - 3 \times 10^{12})$ GeV and $M_1 = (2\, {\rm keV} - 1 \, {\rm MeV)}$. Much of this triangle can be probed by signals of warm dark matter, especially if leptogenesis from $N_2$ decay yields the observed baryon asymmetry. The minimal value of $v_R$ is increased to $10^8 \, {\rm GeV}$ for doublet breaking of $SU(2)_R$, and further to $10^9 \, {\rm GeV}$ if leptogenesis occurs via $N_2$ decay, while the upper bound on $M_1$ is reduced to 100 keV. In addition, there is a component of hot $N_1$ dark matter resulting from the late decay of $N_2 \rightarrow N_1 \ell^+ \ell^-$ that can be probed by future cosmic microwave background observations. Interestingly, the range of $v_R$ allows both precision gauge coupling unification and the Higgs Parity understanding of the vanishing of the Standard Model Higgs quartic at scale $v_R$. Finally, we study freeze-in production of $N_1$ dark matter via the $W_R$ interaction, which allows a much wider range of $(v_R,M_1)$.

An exact parity replicates the Standard Model giving a Mirror Standard Model, SM $\leftrightarrow$ SM$'$. This "Higgs Parity" and the mirror electroweak symmetry are spontaneously broken by the mirror Higgs, $\left\langle H'\right\rangle = v' \gg \left\langle H\right\rangle$, yielding the Standard Model Higgs as a Pseudo-Nambu-Goldstone Boson of an approximate $SU(4)$ symmetry, with a quartic coupling $\lambda_{SM}(v') \sim 10^{-3}$. Mirror electromagnetism is unbroken and dark matter is composed of $e'$ and $\bar{e}'$. Direct detection may be possible via the kinetic mixing portal, and in unified theories this rate is correlated with the proton decay rate. With a high reheat temperature after inflation, the $e'$ dark matter abundance is determined by freeze-out followed by dilution from decays of mirror neutrinos, $\nu' \rightarrow \ell H$. Remarkably, this requires $v' \sim (10^8 - 10^{10})$ GeV, consistent with the Higgs mass, and a Standard Model neutrino mass of $(10^{-2} - 10^{-1})$ eV, consistent with observed neutrino masses. The mirror QCD sector exhibits a first order phase transition producing gravitational waves that may be detected by future observations. Mirror glueballs decay to mirror photons giving dark radiation with $\Delta N_{\rm eff} \sim 0.03 - 0.4$. With a low reheat temperature after inflation, the $e'$ dark matter abundance is determined by freeze-in from the SM sector by either the Higgs or kinetic mixing portal.

We study interactions of cosmological relics, $X$, of mass $m$ and electric charge $qe$ in the galaxy, including thermalization with the interstellar medium, diffusion through inhomogeneous magnetic fields and Fermi acceleration by supernova shock waves. We find that for $m \mathop{}_{\textstyle \sim}^{\textstyle <} 10^{10} q \; ~{\rm GeV}$, there is a large flux of accelerated $X$ in the disk today, with a momentum distribution $\propto 1/p^{2.5}$ extending to $(\beta p)_{max} \sim 5 \times10^4 q \; ~{\rm GeV}$. Even though acceleration in supernova shocks is efficient, ejecting $X$ from the galaxy, $X$ are continually replenished by diffusion into the disk from the halo or confinement region. For $m \mathop{}_{\textstyle \sim}^{\textstyle >} 10^{10} q \; ~{\rm GeV}$, $X$ cannot be accelerated above the escape velocity within the lifetime of the shock. The accelerated $X$ form a component of cosmic rays that can easily reach underground detectors, as well as deposit energies above thresholds, enhancing signals in various experiments. We find that nuclear/electron recoil experiments place very stringent bounds on $X$ at low $q$; for example, $X$ as dark matter is excluded for $q$ above $10^{-9}$ for any $m$. For larger $q$, stringent bounds on the fraction of dark matter that can be $X$ are set by Cherenkov and ionization detectors. Nevertheless, very small $q$ is highly motivated by the kinetic mixing portal, and we identify regions of $(m,q)$ that can be probed by future experiments.