Assuming that Dark Matter (DM) is made of fermions in the sub-GeV mass range with interactions dominated by electromagnetic moments of higher order, such as the electric and magnetic dipoles or the anapole moment, we show that direct detection experiments searching for atomic ionisation events in xenon targets can shed light on whether DM is a Dirac or Majorana particle. Specifically, we find that between about 45 (120) and 610 (1700) signal events are required to reject Majorana DM in favour of Dirac DM with a statistical significance corresponding to 3 (5) standard deviations. The exact number of DM signal events corresponding to a given significance depends on the relative size of the anapole, magnetic dipole and electric dipole contributions to the expected rate of DM-induced atomic ionisations under the Dirac hypothesis. Our conclusions are based on Monte Carlo simulations and the likelihood ratio test. While the use of asymptotic formulae for the latter is standard in many applications, here it requires a non-trivial extension to the case where one of the hypotheses lies on the boundary of the parameter space. Our results constitute a solid proof of concept about the possibility of using direct detection experiments to reject the Majorana DM hypothesis when the DM interactions are dominated by higher-order electromagnetic moments.

Despite strong evidence for the existence of large amounts of dark matter (DM) in our Universe, there is no direct indication of its presence in our own solar system. All estimates of the local DM density, crucial for all direct DM searches, rely on extrapolating estimates on much larger scales. We demonstrate for the first time the possibility of measuring the local DM density with a direct detection experiment. It relies on the assumption that incoming DM particles frequently scatter on terrestrial nuclei prior to detection, inducing an additional time-dependence of the signal. We show that for sub-GeV DM, with a spin-independent DM-proton cross section $\sim10^{-32}\,\text{cm}^2$, future direct detection experiments should be able to reconstruct the local DM density with 20% uncertainty.

We compute the projected sensitivity to dark matter (DM) particles in the sub-GeV mass range of future direct detection experiments using germanium and silicon semiconductor targets. We perform this calculation within the dark photon model for DM-electron interactions using the likelihood ratio as a test statistic, Monte Carlo simulations, and background models that we extract from recent experimental data. We present our results in terms of DM-electron scattering cross section values required to reject the background only hypothesis in favour of the background plus DM signal hypothesis with a statistical significance, $\mathcal{Z}$, corresponding to 3 or 5 standard deviations. We also test the stability of our conclusions under changes in the astrophysical parameters governing the local space and velocity distribution of DM in the Milky Way. In the best-case scenario, when a high-voltage germanium detector with an exposure of $50$ kg-year and a CCD silicon detector with an exposure of $1$ kg-year and a dark current rate of $1\times10^{-7}$ counts/pixel/day have simultaneously reported a DM signal, we find that the smallest cross section value compatible with $\mathcal{Z}=3$ ($\mathcal{Z}=5$) is about $8\times10^{-42}$ cm$^2$ ($1\times10^{-41}$ cm$^2$) for contact interactions, and $4\times10^{-41}$ cm$^2$ ($7\times10^{-41}$ cm$^2$) for long-range interactions. Our sensitivity study extends and refine previous works in terms of background models, statistical methods, and treatment of the underlying astrophysical uncertainties.

In the leading paradigm of modern cosmology, about 80% of our Universe's matter content is in the form of hypothetical, as yet undetected particles. These do not emit or absorb radiation at any observable wavelengths, and therefore constitute the so-called Dark Matter (DM) component of the Universe. Detecting the particles forming the Milky Way DM component is one of the main challenges for astroparticle physics and basic science in general. One promising way to achieve this goal is to search for rare DM-electron interactions in low-background deep underground detectors. Key to the interpretation of this search is the response of detectors' materials to elementary DM-electron interactions defined in terms of electron wave functions' overlap integrals. In this work, we compute the response of atomic argon and xenon targets used in operating DM search experiments to general, so far unexplored DM-electron interactions. We find that the rate at which atoms can be ionized via DM-electron scattering can in general be expressed in terms of four independent atomic responses, three of which we identify here for the first time. We find our new atomic responses to be numerically important in a variety of cases, which we identify and investigate thoroughly using effective theory methods. We then use our atomic responses to set 90% confidence level (C.L.) exclusion limits on the strength of a wide range of DM-electron interactions from the null result of DM search experiments using argon and xenon targets.

We consider direct-detection searches for sub-GeV dark matter via electron scatterings in the presence of large interactions between dark and ordinary matter. Scatterings both on electrons and nuclei in the Earth's crust, atmosphere, and shielding material attenuate the expected local dark matter flux at a terrestrial detector, so that such experiments lose sensitivity to dark matter above some critical cross section. We study various models, including dark matter interacting with a heavy and ultralight dark photon, through an electric dipole moment, and exclusively with electrons. For a dark-photon mediator and an electric dipole interaction, the dark matter-electron scattering cross-section is directly linked to the dark matter-nucleus cross section, and nuclear interactions typically dominate the attenuation process. We determine the exclusion bands for the different dark-matter models from several experiments - SENSEI, CDMS-HVeV, XENON10, XENON100, and DarkSide-50 - using a combination of Monte Carlo simulations and analytic estimates. We also derive projected sensitivities for a detector located at different depths and for a range of exposures, and calculate the projected sensitivity for SENSEI at SNOLAB and DAMIC-M at Modane. Finally, we discuss the reach to high cross sections and the modulation signature of a small balloon- and satellite-borne detector sensitive to electron recoils, such as a Skipper-CCD. Such a detector could potentially probe unconstrained parameter space at high cross sections for a sub-dominant component of dark matter interacting with a massive, but ultralight, dark photon.

If DM-matter scatterings are assumed to occur in a detector's target material, collisions will naturally take place inside the bulk of planets and stars as well. For large cross sections, these scatterings might occur in the Earth or Sun even prior to the detection. In this thesis, we study the impact of these pre-detection scatterings on direct searches of light DM with the use of Monte Carlo (MC) simulations. By simulating the trajectories and scatterings of many individual DM particles in the Earth or Sun, we determine the local distortions of the statistical properties of DM at any detector caused by elastic DM-nucleus collisions. Scatterings inside the Earth distort the underground DM density and velocity distribution. Any detector moves periodically through these inhomogeneities due to the Earth's rotation, and the expected event rate will vary throughout a sidereal day. Using MC simulations, we can determine the exact amplitude and phase of this diurnal modulation for any experiment. For even higher scattering probabilities, collisions in the overburden above the typically underground detectors start to attenuate the incoming DM flux. The critical cross section above which an experiment loses sensitivity to DM itself is determined for a variety of DM-nucleus and DM-electron scattering experiments and different interaction types. Furthermore, we develop the idea that sub-GeV DM particles can enter the Sun, gain kinetic energy by colliding on hot nuclei and get reflected with great speeds. By deriving an analytic expressions for the particle flux from solar reflection via a single scattering, we demonstrate the prospects of future experiments to probe reflected DM and extend their sensitivity to lower masses than accessible by halo DM alone. Finally, we present first results for MC simulations of solar reflections taking into account the effect of multiple scatterings.

Above a critical dark matter-nucleus scattering cross section any terrestrial direct detection experiment loses sensitivity to dark matter, since the Earth crust, atmosphere, and potential shielding layers start to block off the dark matter particles. This critical cross section is commonly determined by describing the average energy loss of the dark matter particles analytically. However, this treatment overestimates the stopping power of the Earth crust. Therefore the obtained bounds should be considered as conservative. We perform Monte Carlo simulations to determine the precise value of the critical cross section for various direct detection experiments and compare them to other dark matter constraints in the low mass regime. In this region we find parameter space where typical underground and surface detectors are completely blind to dark matter. This "hole" in the parameter space can hardly be closed with an increase in the detector exposure. Dedicated surface or high-altitude experiments may be the only way to directly probe this part of the parameter space.

Sub-GeV halo dark matter that enters the Sun can potentially scatter off hot solar nuclei and be ejected much faster than its incoming velocity. We derive an expression for the rate and velocity distribution of these reflected particles taking into account the Sun's temperature and opacity. We further demonstrate that future direct detection experiments could use these energetic reflected particles to probe light dark matter in parameter space that cannot be accessed via ordinary halo dark matter.

Conventional dark matter direct detection experiments set stringent constraints on dark matter by looking for elastic scattering events between dark matter particles and nuclei in underground detectors. However these constraints weaken significantly in the sub-GeV mass region, simply because light dark matter does not have enough energy to trigger detectors regardless of the dark matter-nucleon scattering cross section. Even if future experiments lower their energy thresholds, they will still be blind to parameter space where dark matter particles interact with nuclei strongly enough that they lose enough energy and become unable to cause a signal above the experimental threshold by the time they reach the underground detector. Therefore in case dark matter is in the sub-GeV region and strongly interacting, possible underground scatterings of dark matter with terrestrial nuclei must be taken into account because they affect significantly the recoil spectra and event rates, regardless of whether the experiment probes DM via DM-nucleus or DM-electron interaction. To quantify this effect we present the publicly available Dark Matter Simulation Code for Underground Scatterings (DaMaSCUS), a Monte Carlo simulator of DM trajectories through the Earth taking underground scatterings into account. Our simulation allows the precise calculation of the density and velocity distribution of dark matter at any detector of given depth and location on Earth. The simulation can also provide the accurate recoil spectrum in underground detectors as well as the phase and amplitude of the diurnal modulation caused by this shadowing effect of the Earth, ultimately relating the modulations expected in different detectors, which is important to decisively conclude if a diurnal modulation is due to dark matter or an irrelevant background.

A well-studied possibility is that dark matter may reside in a sector secluded from the Standard Model, except for the so-called photon portal: kinetic mixing between the ordinary and dark photons. Such interactions can be probed at dark matter direct detection experiments, and new experimental techniques involving detection of dark matter-electron scattering offer new sensitivity to sub-GeV dark matter. Typically however it is implicitly assumed that the dark matter is not altered as it traverses the Earth to arrive at the detector. In this paper we study in detail the effects of terrestrial stopping on dark photon models of dark matter, and find that they significantly reduce the sensitivity of XENON10 and DAMIC. In particular we find that XENON10 only excludes masses in the range (5-3000) MeV while DAMIC only probes (20-50) MeV. Their corresponding cross section sensitivity is reduced to a window of cross sections between $(5\times 10^{-38}-10^{-30})~{\rm cm}^{2}$ for XENON10 and a small window around $\sim 10^{-31}~{\rm cm}^{2}$ for DAMIC. We also examine implications for a future DAMIC run.

We revisit astrophysical constraints on models with an ultra-light gravitino, in particular extending the analysis to more general models and the case of R-parity breaking. These constraints allow to restrict the value of the gravitino mass depending on the masses of scalar moduli. We perform both a frequentist and bayesian analysis and we find comparable results, even if the bayesian analysis is in part affected by volume effects. Only a small window of gravitino masses can be excluded by the SN constraints in a model-independent way, while limits obtained with a definite assumption on the scalar masses result more stringent.