The search for dark matter weakly interacting massive particles with noble liquids has probed masses down and below a GeV/c^2. The ultimate limit is represented by the experimental threshold on the energy transfer to the nuclear recoil. Currently, the experimental sensitivity has reached a threshold equivalent to a few ionization electrons. In these conditions, the contribution of a Bremsstrahlung photon or a so-called Migdal electron due to the sudden acceleration of a nucleus after a collision might be sizable. In the present work, we use a Bayesian approach to study how these effects can be exploited in experiments based on liquid argon detectors. In particular, taking inspiration from the DarkSide-50 public spectra, we develop a simulated experiment to show how the Migdal electron and the Bremsstrahlung photon allow to push the experimental sensitivity down to masses of 0.1 GeV/c^2, extending the search region for dark matter particles of previous results. For these masses we estimate the effect of the Earth shielding that, for strongly interacting dark matter, makes any detector blind. Finally, we show how the sensitivity scales for higher exposure.

Proportional electroluminescence (EL) in noble gases is used in two-phase detectors for dark matter search to record (in the gas phase) the ionization signal induced by particle scattering in the liquid phase. The "standard" EL mechanism is considered to be due to noble gas excimer emission in the vacuum ultraviolet (VUV). In addition to it, there are two alternative mechanisms, producing light in the visible and near infrared (NIR) range. The first is due to bremsstrahlung of electrons scattered on neutral atoms (so-called neutral bremsstrahlung, NBrS). The second, responsible for electron avalanche scintillations in the NIR at higher electric fields, is due to transitions between excited atomic states. In this work, we have for the first time demonstrated two alternative techniques of the optical readout of two-phase argon detectors, in the visible and NIR range, using a SiPM matrix and electroluminescence due to either neutral bremsstrahlung or avalanche scintillations effects. The amplitude yield and position resolution were measured for these readout techniques, which allowed to assess the detection threshold for electron and nuclear recoils in two-phase argon detectors for dark matter search. It should be emphasized that this is the first practical application of the NBrS effect in detection science.

In a recent paper [arXiv:2002.00459], Buttazzo et al. show how the annually modulated rate of the DAMA experiments can be possibly interpreted as an artefact due to the interplay between a time-dependent background and the method to account for it. In this work, we compare this hypothesis against the sinusoidal dark matter signal as proposed by the DAMA collaboration. We produce in a Bayesian approach a quantitative comparison of how much the experimental observations are in support of each hypothesis. Our conclusions are that the odds against the hypothesis of a time varying background being responsible for the annual modulation are decreased by a Bayes factor larger than 10^8 after considering the public available data of the DAMA/NaI and DAMA/LIBRA experiments. In this work we also elaborate on general aspects of the analysis procedure. Indeed, in order to optimise the background subtraction procedure, the DAMA collaboration only considers data-taking cycles with a duration of roughly one year. We argue that any data-taking cycle is informative, and we propose a procedure to include this effect, as well as the possibility to include a slowly varying component for the background.

Large liquid argon detectors offer one of the best avenues for the detection of galactic weakly interacting massive particles (WIMPs) via their scattering on atomic nuclei. The liquid argon target allows exquisite discrimination between nuclear and electron recoil signals via pulse-shape discrimination of the scintillation signals. Atmospheric argon (AAr), however, has a naturally occurring radioactive isotope, $^{39}$Ar, a $\beta$ emitter of cosmogenic origin. For large detectors, the atmospheric $^{39}$Ar activity poses pile-up concerns. The use of argon extracted from underground wells, deprived of $^{39}$Ar, is key to the physics potential of these experiments. The DarkSide-20k dark matter search experiment will operate a dual-phase time projection chamber with 50 tonnes of radio-pure underground argon (UAr), that was shown to be depleted of $^{39}$Ar with respect to AAr by a factor larger than 1400. Assessing the $^{39}$Ar content of the UAr during extraction is crucial for the success of DarkSide-20k, as well as for future experiments of the Global Argon Dark Matter Collaboration (GADMC). This will be carried out by the DArT in ArDM experiment, a small chamber made with extremely radio-pure materials that will be placed at the centre of the ArDM detector, in the Canfranc Underground Laboratory (LSC) in Spain. The ArDM LAr volume acts as an active veto for background radioactivity, mostly $\gamma$-rays from the ArDM detector materials and the surrounding rock. This article describes the DArT in ArDM project, including the chamber design and construction, and reviews the background required to achieve the expected performance of the detector.

We report measurements of the charged daughter fraction of $^{218}$Po as a result of the $^{222}$Rn alpha decay, and the mobility of $^{218}$Po$^+$ ions, using radon-polonium coincidences from the $^{238}$U chain identified in 532 live-days of DarkSide-50 WIMP-search data. The fraction of $^{218}$Po that is charged is found to be 0.37$\pm$0.03 and the mobility of $^{218}$Po$^+$ is (8.6$\pm$0.1)$\times$10$^{-4}$$\frac{\text{cm}^2}{\text{Vs}}$.

The DarkSide-50 direct-detection dark matter experiment is a dual-phase argon time projection chamber operating at Laboratori Nazionali del Gran Sasso. This paper reports on the blind analysis of a (16,660+-270) kg d exposure using a target of low-radioactivity argon extracted from underground sources. We find no events in the dark matter selection box and set a 90% C.L. upper limit on the dark matter-nucleon spin-independent cross section of 1.14E-44 cm^2 (3.78E-44 cm^2, 3.43E-43 cm^2) for a WIMP mass of 100 GeV/c^2 (1 TeV/c^2, 10 TeV/c^2).

We present new constraints on sub-GeV dark matter particles scattering off electrons in argon based on an analysis of ionization signal data from the DarkSide-50 detector.

We present the results of a search for dark matter WIMPs in the mass range below 20 GeV/c^2 using a target of low-radioactivity argon. The data were obtained using the DarkSide-50 apparatus at Laboratori Nazionali del Gran Sasso (LNGS). The analysis is based on the ionization signal, for which the DarkSide-50 time projection chamber is fully efficient at 0.1 keVee. The observed rate in the detector at 0.5 keVee is about 1.5 events/keVee/kg/day and is almost entirely accounted for by known background sources. We obtain a 90% C.L. exclusion limit above 1.8 GeV/c^2 for the spin-independent cross section of dark matter WIMPs on nucleons, extending the exclusion region for dark matter below previous limits in the range 1.8-6 GeV/c^2.

We measure the power spectrum of the galaxy distribution in the ESO Slice Project (ESP) galaxy redshift survey. We develope a technique to describe the survey window function analytically, and then deconvolve it from the measured power spectrum using a variant of the Lucy method. We test the whole deconvolution procedure on ESP mock catalogues drawn from large N-body simulations, and find that it is reliable for recovering the correct amplitude and shape of $P(k)$ at $k> 0.065 h$ Mpc$^{-1}$. In general, the technique is applicable to any survey composed by a collection of circular fields with arbitrary pattern on the sky, as typical of surveys based on fibre spectrographs. The estimated power spectrum has a well-defined power-law shape $k^n$ with $n\simeq -2.2$ for $k\ge 0.2 h$ Mpc$^{-1}$, and a smooth bend to a flatter shape ($n\simeq -1.6$) for smaller $k$'s. The smallest wavenumber, where a meaningful reconstruction can be performed ($k\sim 0.06 h$ Mpc$^{-1}$), does not allow us to explore the range of scales where other power spectra seem to show a flattening and hints for a turnover. We also find, by direct comparison of the Fourier transforms, that the estimate of the two-point correlation function $\xi(s)$ is much less sensitive to the effect of a problematic window function as that of the ESP, than the power spectrum. Comparison to other surveys shows an excellent agreement with estimates from blue-selected surveys. In particular, the ESP power spectrum is virtually indistinguishable from that of the Durham-UKST survey over the common range of $k$'s, an indirect confirmation of the quality of the deconvolution technique applied.

Although poor for small dynamic scales, the Particle-Mesh (PM) model allows in astrophysics good insight for large dynamic scales at low computational cost. Furthermore, it is possible to employ a very high number of particles to get high mass resolution. These properties could be exploited by suitable parallelization of the algorithm. In PM there are two types of data: the particle data, i.e. position and velocity, which are stored in one-dimensional arrays of N elements, and the grid data, i.e. density and force, which are stored in three-dimensional arrays $M\times M\times M$ in size. Since individual particles can change cell under the action of gravitational force, parallelization is a real challenge on parallel machine and must account for the machine architecture. We have implemented a dynamic work distribution through agenda parallelism. By subdividing the work in small tasks, the implementation is well balanced, scalable and efficient also for clustered particle distributions. In this contribution we describe this efficient, load balanced, parallel implementation of PM algorithm on Cray T3E at CINECA and show its performances on cosmological simulation results.

The possible cosmological effects of primordial fluctuation corrections to the evolution equation of matter obtained from the Wheeler--De Witt equation are explored. In particular, both the metric and a scalar matter field are expanded around their homogeneous values and the corrections induced on the scalar field fluctuation spectrum are perturbatively estimated. Finally, results of a preliminary numerical simulation to investigate the effects on large--scale structure formation are presented.

We analyse the large--scale velocity field obtained by N--body simulations of cold dark matter (CDM) models with non--Gaussian primordial density fluctuations, considering models with both positive and negative primordial skewness in the density fluctuation distribution. We study the velocity probability distribution and calculate the dependence of the bulk flow, one--point velocity dispersion and Cosmic Mach Number on the filtering size. We find that the sign of the primordial skewness of the density field provides poor discriminatory power on the evolved velocity field. All non--Gaussian models here considered tend to have lower velocity dispersion and bulk flow than the standard Gaussian CDM model, while the Cosmic Mach Number turns out to be a poor statistic in characterizing the models. Next, we compare the large--scale velocity field of a composite sample of optically selected galaxies as described by the Local Group properties, bulk flow, velocity correlation function and Cosmic Mach Number, with the velocity field of mock catalogues extracted from the N--body simulations. The comparison does not clearly permit to single out a best model: the standard Gaussian model is however marginally preferred by the maximum likelihood analysis.

We investigate the evolution of the skewness of the distribution of density fluctuations in CDM models with both Gaussian and non--Gaussian initial fluctuations. We show that the method proposed by Coles \& Frenk (1991), which uses the skewness of galaxy counts to test the hypothesis of Gaussian primordial density fluctuations, is a potentially powerful probe of initial conditions. As expected, the mass distribution in models with initially non--Gaussian fluctuations shows systematic departures from the Gaussian behaviour on intermediate to large scales. We investigate the effect of peculiar velocity distortions and normalisation upon the relationship between skewness and variance. These effects are generally small for the models we consider. Comparing our results to the QDOT measurements of the skewness, we find that our initially positive--skew models are clearly excluded by this analysis, but the available data do not rule out the negative--skew models.

We generate artificial Lick maps using $N$--body simulations and compare the angular correlation function, $w(\vartheta)$, measured from the simulations with the APM correlation. For the Gaussian CDM model, neither the standard biassed model nor a more evolved model (as suggested by the COBE data), reproduce the correlations on large angular scales. We come to a similar conclusion about CDM models with positively skewed initial fluctuation distributions. In contrast, models with initially negatively skewed fluctuations produce a $w(\vartheta)$ that declines much more gently on large scales. Such models are therefore in principle capable of reconciling the lack of large--scale power of the CDM spectrum with the observed clustering of APM galaxies.