Indirect detection of dark matter via its annihilation products is a key technique in the search for dark matter in the form of weakly interacting massive particles (WIMPs). Strong constraints exist on the annihilation of WIMPs to highly visible Standard Model final states such as photons or charged particles. In the case of s-wave annihilation, this typically eliminates thermal relic cross sections for dark matter of mass below $\cal{O}$(10) GeV. However, such limits typically neglect the possibility that dark matter may annihilate to assumed invisible or hard-to-detect final states, such as neutrinos. This is a difficult paradigm to probe due to the weak neutrino interaction cross section. Considering dark matter annihilation in the Galactic halo, we study the prospects for indirect detection using the Hyper-Kamiokande (HyperK) neutrino experiment, for dark matter of mass below 1 GeV. We undertake a dedicated simulation of the HyperK detector, which we benchmark against results from the similar Super-Kamiokande experiment and HyperK physics projections. We provide projections for the annihilation cross-sections that can be probed by HyperK for annihilation to muon or neutrino final states, and discuss uncertainties associated with the dark matter halo profile. For neutrino final states, we find that HyperK is sensitive to thermal annihilation cross-sections for dark matter with mass around 20 MeV, assuming an NFW halo profile. We also discuss the effects of Gadolinium enhancement, and prospects for improving the reach at low mass.

Neutron stars provide a cosmic laboratory to study the nature of dark matter particles and their interactions. Dark matter can be captured by neutron stars via scattering interactions, in which kinetic energy is transferred to the star. This can have a number of observational consequences, such as the heating of old neutron stars to infra-red temperatures. Previous treatments of the capture process have employed various approximation or simplifications. We present here an improved treatment of dark matter capture, valid for a wide dark matter mass range, that correctly incorporates all relevant physical effects. These include gravitational focusing, a fully relativistic scattering treatment, Pauli blocking, neutron star (NS) opacity and multi-scattering effects. We provide general expressions that enable the exact capture rate to be calculated numerically, and derive simplified expressions that are valid for particular interaction types or mass regimes and that greatly increase the computational efficiency. Our formalism is applicable to the scattering of dark matter from any NS constituents, or to the capture of dark matter in other compact objects.

In this paper we study the presence of generalized blind spots, i.e. regions of the parameter space where the spin-independent cross section for dark matter direct detection is suppressed, in the context of a generic 2HDM and a minimal fermionic Higgs-portal dark sector. To this end, we derive analytical expressions for the couplings of the dark matter to the light and heavy Higgses, and thus for the blind spot solutions. Unlike the case of a standard Higgs sector, blind spots can occur even without a cancellation between different contributions, while keeping unsuppressed and efficient the annihilation processes in the early Universe involving Higgs states. As a consequence, the allowed parameter space is dramatically enhanced.

Dark matter particles will be captured in neutron stars if they undergo scattering interactions with nucleons or leptons. These collisions transfer the dark matter kinetic energy to the star, resulting in appreciable heating that is potentially observable by forthcoming infrared telescopes. While previous work considered scattering only on nucleons, neutron stars contain small abundances of other particle species, including electrons and muons. We perform a detailed analysis of the neutron star kinetic heating constraints on leptophilic dark matter. We also estimate the size of loop induced couplings to quarks, arising from the exchange of photons and Z bosons. Despite having relatively small lepton abundances, we find that an observation of an old, cold, neutron star would provide very strong limits on dark matter interactions with leptons, with the greatest reach arising from scattering off muons. The projected sensitivity is orders of magnitude more powerful than current dark matter-electron scattering bounds from terrestrial direct detection experiments.

Semi-analytic models are best suited to compare galaxy formation and evolution theories with observations. These models rely heavily on halo merger trees, and their realistic features (i.e., no drastic changes on halo mass or jumps on physical locations). Our aim is to provide a new framework for halo merger tree generation that takes advantage of the results of large volume simulations, with a modest computational cost. We treat halo merger tree construction as a matrix generation problem, and propose a Generative Adversarial Network that learns to generate realistic halo merger trees. We evaluate our proposal on merger trees from the EAGLE simulation suite, and show the quality of the generated trees.

Genetic Algorithms (GAs) are explored as a tool for probing new physics with high dimensionality. We study the 19-dimensional pMSSM, including experimental constraints from all sources and assessing the consistency of potential signals of new physics. We show that GAs excel at making a fast and accurate diagnosis of the cross-compatibility of a set of experimental constraints in such high dimensional models. In the case of the pMSSM, it is found that only ${\cal O}(10^4)$ model evaluations are required to obtain a best fit point in agreement with much more costly MCMC scans. This efficiency allows higher dimensional models to be falsified, and patterns in the spectrum identified, orders of magnitude more quickly. As examples of falsification, we consider the muon anomalous magnetic moment, and the Galactic Centre gamma-ray excess observed by Fermi-LAT, which could in principle be explained in terms of neutralino dark matter. We show that both observables cannot be explained within the pMSSM, and that they provide the leading contribution to the total goodness of the fit, with $\chi^2_{\delta a_\mu^{\mathrm{SUSY}}}\approx12$ and $\chi^2_{\rm GCE}\approx 155$, respectively.

The latest experimental results from the LHC and dark matter (DM) searches suggest that the parameter space allowed in supersymmetric theories is subject to strong reductions. These bounds are especially constraining for scenarios entailing light DM particles. Previous studies have shown that light neutralino DM in the Minimal Supersymmetric Standard Model (MSSM), with parameters defined at the electroweak scale, is still viable when the low energy spectrum of the model features light sleptons, in which case, the relic density constraint can be fulfilled. In view of this, we have investigated the viability of light neutralinos as DM candidates in the MSSM, with parameters defined at the grand unification scale. We have analysed the optimal choices of non-universalities in the soft supersymmetry-breaking parameters for both, gauginos and scalars, in order to avoid the stringent experimental constraints. We show that light neutralinos, with a mass as low as 25 GeV, are viable in supergravity scenarios if the gaugino mass parameters at high energy are very non universal, while the scalar masses can remain of the same order. These scenarios typically predict a very small cross section of neutralinos off protons and neutrons, thereby being very challenging for direct detection experiments. However, a potential detection of smuons and selectrons at the LHC, together with a hypothetical discovery of a gamma-ray signal from neutralino annihilations in dwarf spheroidal galaxies could shed light on this kind of solutions. Finally, we have investigated the naturalness of these scenarios, taking into account all the potential sources of tuning. Besides the electroweak fine-tuning, we have found that the tuning to reproduce the correct DM relic abundance and that to match the measured Higgs mass can also be important when estimating the total degree of naturalness.

We study spectral features in the gamma-ray emission from dark matter (DM) annihilation in the Next-to-Minimal Supersymmetric Standard Model (NMSSM), with either neutralino or right-handed (RH) sneutrino DM. We perform a series of scans over the NMSSM parameter space, compute the DM annihilation cross section into two photons and the contribution of box-shaped features, and compare them with the limits derived from the Fermi-LAT search for gamma-ray lines using the latest Pass 8 data. We implement the LHC bounds on the Higgs sector and on the masses of supersymmetric particles as well as the constraints on low-energy observables. We also consider the recent upper limits from the Fermi-LAT satellite on the continuum gamma-ray emission from dwarf spheroidal galaxies (dSphs). We show that in the case of the RH sneutrino the constraint on gamma-ray spectral features can be more stringent than the dSph bounds. This is due to the Breit-Wigner enhancement near the ubiquitous resonances with a CP even Higgs and the contribution of scalar and pseudoscalar Higgs final states to box-shaped features. By contrast, for neutralino DM, the di-photon final state is only enhanced in the resonance with a $Z$ boson and box-shaped features are even more suppressed. Therefore, the observation of spectral features could constitute a discriminating factor between both models. In addition, we compare our results with direct DM searches, including the SuperCDMS and LUX limits on the elastic DM-nucleus scattering cross section and show that some of these scenarios would be accessible to next generation experiments. Thus, our findings strengthen the idea of complementarity among distinct DM search strategies.

We show that the right-handed (RH) sneutrino in the NMSSM can account for the observed excess in the Fermi-LAT spectrum of gamma rays from the Galactic Centre, while fulfilling all the current experimental constraints from the LHC as well as from direct and indirect dark matter searches. We have explored the parameter space of this scenario, computed the gamma ray spectrum for each phenomenologically viable solution and then performed a chi^2 fit to the excess. Unlike previous studies based on model independent interpretations, we have taken into account the full annihilation spectrum, without assuming pure annihilation channels. Furthermore, we have incorporated limits from direct detection experiments, LHC bounds and also the constraints from Fermi-LAT on dwarf spheroidal galaxies (dSphs) and gamma ray spectral lines. In addition, we have estimated the effect of the most recent Fermi-LAT reprocessed data (Pass 8). In general, we obtain good fits to the GCE when the RH sneutrino annihilates mainly into pairs of light singlet-like scalar or pseudoscalar Higgs bosons that subsequently decay in flight, producing four-body final states and spectral features that improve the goodness of the fit at large energies. The best fit (chi^2=20.8) corresponds to a RH sneutrino with a mass of 64 GeV which annihilates preferentially into a pair of light singlet-like pseudoscalar Higgs bosons (with masses of order 60 GeV). Besides, we have analysed other channels that also provide good fits to the excess. Finally, we discuss the implications for direct and indirect detection searches paying special attention to the possible appearance of gamma ray spectral features in near future Fermi-LAT analyses, as well as deviations from the SM-like Higgs properties at the LHC. Remarkably, many of the scenarios that fit the GCE can also be probed by these other complementary techniques.

We present a detailed analysis of the local evolution of 206 Lagrangian Volumes (LVs) selected at high redshift around galaxy seeds, identified in a large-volume $\Lambda$ cold dark matter ($\Lambda$CDM) hydrodynamical simulation. The LVs have a mass range of $1 - 1500 \times 10^{10} M_\odot$. We follow the dynamical evolution of the density field inside these initially spherical LVs from $z=10$ up to $z_{\rm low} = 0.05$, witnessing highly non-linear, anisotropic mass rearrangements within them, leading to the emergence of the local cosmic web (CW). These mass arrangements have been analysed in terms of the reduced inertia tensor $I_{ij}^r$, focusing on the evolution of the principal axes of inertia and their corresponding eigendirections, and paying particular attention to the times when the evolution of these two structural elements declines. In addition, mass and component effects along this process have also been investigated. We have found that deformations are led by dark matter dynamics and they transform most of the initially spherical LVs into prolate shapes, i.e. filamentary structures. An analysis of the individual freezing-out time distributions for shapes and eigendirections shows that first most of the LVs fix their three axes of symmetry (like a skeleton) early on, while accretion flows towards them still continue. Very remarkably, we have found that more massive LVs fix their skeleton earlier on than less massive ones. We briefly discuss the astrophysical implications our findings could have, including the galaxy mass-morphology relation and the effects on the galaxy-galaxy merger parameter space, among others.

Recent results from direct and indirect searches for dark matter (DM) have motivated the study of particle physics models that can provide weakly interacting massive particles (WIMPs) in the mass range 1-50 GeV. Viable candidates for light WIMP DM must fulfil stringent constraints. On the one hand, the observation at the LHC of a Higgs boson with Standard Model properties set an upper bound on the coupling of light DM particles to the Higgs, thereby making it difficult to reproduce the correct relic abundance. On the other hand, the recent results from direct searches in the CDMSlite, SuperCDMS and LUX experiments have set upper constraints on the DM scattering cross section. In this paper, we investigate the viability of light right-handed sneutrino DM in the Next-to-Minimal Supersymmetric Model (NMSSM) in the light of these constraints. To this aim, we have carried out a scan in the NMSSM parameter space, imposing experimental bounds on the Higgs sector and low-energy observables, such as the muon anomalous magnetic moment and branching ratios of rare decays. We demonstrate that the enlarged Higgs sector of the NMSSM, together with the flexibility provided by the RH sneutrino parameters, make it possible to obtain viable RH sneutrino DM with a mass as light as 2 GeV. We have also considered the upper bounds on the annihilation cross section from Fermi LAT data on dwarf spheroidal galaxies, and extracted specific examples with mass in the range 8-50 GeV that could account for the apparent low-energy excess in the gamma-ray emission at the Galactic Centre. Then, we have computed the theoretical predictions for the elastic scattering cross-section of RH sneutrinos. Finally, after imposing the recent bounds from SuperCDMS and LUX, we have found a wide area of the parameter space that could be probed by future low-threshold direct detection experiments.