We investigate the cosmology of the minimal model of neutral naturalness, the mirror Twin Higgs. The softly-broken mirror symmetry relating the Standard Model to its twin counterpart leads to significant dark radiation in tension with BBN and CMB observations. We quantify this tension and illustrate how it can be mitigated in several simple scenarios that alter the relative energy densities of the two sectors while respecting the softly-broken mirror symmetry. In particular, we consider both the out-of-equilibrium decay of a new scalar as well as reheating in a toy model of twinned inflation, Twinflation. In both cases the dilution of energy density in the twin sector does not merely reconcile the existence of a mirror Twin Higgs with cosmological constraints, but predicts contributions to cosmological observables that may be probed in current and future CMB experiments. This raises the prospect of discovering evidence of neutral naturalness through cosmology rather than colliders.

An excess in gamma-rays emanating from the galactic centre has recently been observed in the Fermi-LAT data. This signal can be interpreted as resulting from WIMP annihilation, with the spectrum well-fit by dark matter annihilating dominantly into either bottom-quark or Higgs pairs. Supersymmetric models provide a well-motivated framework to study the implications of this signal in these channels. With a neutralino dark matter candidate, the gamma-ray excess cannot be easily accommodated in the minimal supersymmetric model, which in any case requires tuning below the percent level to explain the observed Higgs mass. Instead we are naturally led to consider the next-to-minimal model with a singlet superfield. This not only allows for the annihilation channel into bottom-quark pairs to be implemented, but also provides new possibilities for annihilation into Higgs-pseudoscalar pairs. We show that the fit to the gamma-ray excess for the Higgs-pseudoscalar channel can be just as good as for annihilation into bottom-quark pairs. Moreover, in the parameter range of interest, the next-to-minimal supersymmetric model solves the mu-problem and can explain the 125 GeV Higgs mass with improved naturalness. We also consider an extension by adding a right-handed neutrino superfield with the right-handed sneutrino acting as a dark matter candidate. Interestingly, this allows for the annihilation into pseudoscalar pairs which also provide a good fit to the gamma-ray excess. Furthermore, in the case of a neutralino LSP, the late decay of a sneutrino NLSP can non-thermally produce the observed relic abundance. Finally, the WIMP annihilation into scalar pairs allows for the possibility of detecting the Higgs or pseudoscalar decay into two photons, providing a smoking-gun signal of the model.