In this paper, as a preparation of developing data analysis procedures for using 3-dimensional information offered by directional Dark Matter (DM) detection experiments, we study the patterns of the angular distribution of the Monte Carlo-generated 3-D velocity of halo Weakly Interacting Massive Particles (WIMPs) as well as apply the Bayesian fitting technique to reconstruct the radial distribution of the 3-D WIMP velocity. Besides the diurnal modulation of the angular WIMP velocity distribution, the so-called "directionality" of DM signals proposed in literature, we will also demonstrate possible "annual" modulations of both of the angular and the radial distributions of the 3-D WIMP velocity. Our Bayesian reconstruction results of (the annual modulation of) the radial WIMP velocity distribution will also be discussed in detail. For readers' reference, the angular distribution patterns of the 3-D WIMP velocity in the "laboratory (location)-dependent" reference frames of several underground laboratories are given in the Appendix.

Cosmic ray electron (CRE) acceleration and cooling are important physical processes in astrophysics. We develop an approximative framework to treat CRE physics in the parallel smoothed particle hydrodynamics code Gadget-3. In our methodology, the CRE spectrum of each fluid element is approximated by a single power-law distribution with spatially varying amplitude, upper cut-off, lower cut-off, and spectral index. We consider diffusive shock acceleration to be the source of injection, and oppositely the sinking processes is attributed to synchrotron radiation, inverse Compton scatters, and Coulomb scatters. The adiabatic gains and losses are also included. We show that our formalism produces the energy and pressure with an accuracy of $ > 90\%$ for a free cooling CRE spectrum. Both slope and intensity of the radio emission computed from the CRE population given by our method in cosmological hydro-simulation coincide well with observations, and our results also show that relaxed clusters have lower fluxes. Finally, we investigate several impacts of the CRE processes on the cosmological hydro-simulation, we find that: (1) the pressure of the CRE spectrum is very small and can be ignored in hydro-simulation, (2) the impacts of the CRE processes on the gas phase-space state of hydro-simulation is up to $3\%$, (3) the CRE processes induce a $5\%$ influence on the mass function in the mass range $10^{12} -10^{13} h^{-1} M_{\odot}$, (4) The gas temperature of massive galaxy cluster is influenced by the CRE processes up to $\sim 10\%$.

In this article, I present our (Monte Carlo) simulation results of the angular distribution of the 3-dimensional WIMP velocity, in particular, a possible "annual" modulation and the diurnal modulation proposed in literature, for different underground laboratories.

The overwhelming foreground contamination is one of the primary impediments to probing the Epoch of Reionization (EoR) through measuring the redshifted 21 cm signal. Among various foreground components, radio halos are less studied and their impacts on the EoR observations are still poorly understood. In this work, we employ the Press-Schechter formalism, merger-induced turbulent re-acceleration model, and the latest SKA1-Low layout configuration to simulate the SKA "observed" images of radio halos. We calculate the one-dimensional power spectra from simulated images and find that radio halos can be about $10^4$, $10^3$ and $10^{2.5}$ times more luminous than the EoR signal on scales of $0.1\,\text{Mpc}^{-1} < k < 2\,\text{Mpc}^{-1}$ in the 120-128, 154-162, and 192-200 MHz bands, respectively. By examining the two-dimensional power spectra inside properly defined EoR windows, we find that the power leaked by radio halos can still be significant, as the power ratios of radio halos to the EoR signal on scales of $0.5\,\text{Mpc}^{-1} \lesssim k \lesssim 1\,\text{Mpc}^{-1}$ can be up to about 230-800%, 18-95%, and 7-40% in the three bands, when the 68% uncertainties caused by the variation of the number density of bright radio halos are considered. Furthermore, we find that radio halos located inside the far side-lobes of the station beam can also impose strong contamination within the EoR window. In conclusion, we argue that radio halos are severe foreground sources and need serious treatments in future EoR experiments.

When applying the foreground removal methods to uncover the faint cosmological signal from the epoch of reionization (EoR), the foreground spectra are assumed to be smooth. However, this assumption can be seriously violated in practice since the unresolved or mis-subtracted foreground sources, which are further complicated by the frequency-dependent beam effects of interferometers, will generate significant fluctuations along the frequency dimension. To address this issue, we propose a novel deep-learning-based method that uses a 9-layer convolutional denoising autoencoder (CDAE) to separate the EoR signal. After being trained on the SKA images simulated with realistic beam effects, the CDAE achieves excellent performance as the mean correlation coefficient ($\bar{\rho}$) between the reconstructed and input EoR signals reaches $0.929 \pm 0.045$. In comparison, the two representative traditional methods, namely the polynomial fitting method and the continuous wavelet transform method, both have difficulties in modelling and removing the foreground emission complicated with the beam effects, yielding only $\bar{\rho}_{\text{poly}} = 0.296 \pm 0.121$ and $\bar{\rho}_{\text{cwt}} = 0.198 \pm 0.160$, respectively. We conclude that, by hierarchically learning sophisticated features through multiple convolutional layers, the CDAE is a powerful tool that can be used to overcome the complicated beam effects and accurately separate the EoR signal. Our results also exhibit the great potential of deep-learning-based methods in future EoR experiments.

We present a morphological classification of 14,245 radio active galactic nuclei (AGNs) into six types, i.e., typical Fanaroff--Riley Class I / II (FRI/II), FRI/II-like bent-tailed, X-shaped radio galaxy, and ringlike radio galaxy, by designing a convolutional neural network (CNN) based autoencoder, namely MCRGNet, and applying it to a labeled radio galaxy (LRG) sample containing 1442 AGNs and an unlabeled radio galaxy (unLRG) sample containing 14,245 unlabeled AGNs selected from the Best--Heckman sample. We train MCRGNet and implement the classification task by a three-step strategy, i.e., pre-training, fine-tuning, and classification, which combines both unsupervised and supervised learnings. A four-layer dichotomous tree is designed to classify the radio AGNs, which leads to a significantly better performance than the direct six-type classification. On the LRG sample, our MCRGNet achieves a total precision of $\sim 93\%$ and an averaged sensitivity of $\sim 87\%$, which are better than those obtained in previous works. On the unLRG sample, whose labels have been human-inspected, the neural network achieves a total precision of $\sim 80\%$. Also, using the Sloan Digital Sky Survey (SDSS) Data Release 7 (DR7) to calculate the $r$-band absolute magnitude ($M_\mathrm{opt}$) and using the flux densities to calculate the radio luminosity ($L_\mathrm{radio}$), we find that the distributions of the unLRG sources on the $L_\mathrm{radio}$--$M_\mathrm{opt}$ plane do not show an apparent redshift evolution and could confirm with a sufficiently large sample that there could not exist an abrupt separation between FRIs and FRIIs as reported in some previous works.

We choose the bright compact group HCG 62, which was found to exhibit both excess X-ray emission and high Fe abundance to the southwest of its core, as an example to study the impact of mergers on chemical enrichment in the intragroup medium. We first reanalyze the high-quality Chandra and XMM-Newton archive data to search for the evidence for additional SN II yields, which is expected as a direct result of the possible merger-induced starburst. We reveal that, similar to the Fe abundance, the Mg abundance also shows a high value in both the innermost region and the southwest substructure, forming a high-abundance plateau, meanwhile all the SN Ia and SN II yields show rather flat distributions in $>0.1r_{200}$ in favor of an early enrichment. Then we carry out a series of idealized numerical simulations to model the collision of two initially isolated galaxy groups by using the TreePM-SPH GADGET-3 code. We find that the observed X-ray emission and metal distributions, as well as the relative positions of the two bright central galaxies with reference to the X-ray peak, can be well reproduced in a major merger with a mass ratio of 3 when the merger-induced starburst is assumed. The `best-match' snapshot is pinpointed after the third pericentric passage when the southwest substructure is formed due to gas sloshing. By following the evolution of the simulated merging system, we conclude that the effects of such a major merger on chemical enrichment are mostly restricted within the core region when the final relaxed state is reached.

In this paper, we revisit our model-independent methods developed for reconstructing properties of Weakly Interacting Massive Particles (WIMPs) by using measured recoil energies from direct Dark Matter detection experiments directly and take into account more realistically non-negligible threshold energy. All expressions for reconstructing the mass and the (ratios between the) spin-independent and the spin-dependent WIMP-nucleon couplings have been modified. We focus on low-mass (m_chi <~ 15 GeV) WIMPs and present the numerical results obtained by Monte Carlo simulations. Constraints caused by non-negligible threshold energy and technical treatments for improving reconstruction results will also be discussed.

In this paper, we investigate the modification of our expressions developed for the model-independent data analysis procedure of the reconstruction of the (time-averaged) one-dimensional velocity distribution of Galactic Weakly Interacting Massive Particles (WIMPs) with a non-negligible experimental threshold energy. Our numerical simulations show that, for a minimal reconstructable velocity of as high as O(200) km/s, our model-independent modification of the estimator for the normalization constant could provide precise reconstructed velocity distribution points to match the true WIMP velocity distribution with a <~ 10% bias.

We report on the optical spectroscopic analysis of a sample of 99 low-luminosity quasi-stellar objects (LLQSOs) at $z\leq 0.06$ base the Hamburg/ESO QSO survey (HES). The LLQSOs presented here offer the possibility of studying the faint end of the QSO population at smaller cosmological distances and, therefore, in greater detail. A small number of our LLQSO present no broad component. Two sources show double broad components, whereas six comply with the classic NLS1 requirements. As expected in NLR of broad line AGNs, the [S{\sc{ii}}]$-$based electron density values range between 100 and 1000 N$_{e}$/cm$^{3}$. Using the optical characteristics of Populations A and B, we find that 50\% of our sources with H$\beta$ broad emission are consistent with the radio-quiet sources definition. The remaining sources could be interpreted as low-luminosity radio-loud quasar. The BPT-based classification renders an AGN/Seyfert activity between 50 to 60\%. For the remaining sources, the possible star burst contribution might control the LINER and HII classification. Finally, we discuss the aperture effect as responsible for the differences found between data sets, although variability in the BLR could play a significant role as well.

In this paper, we give a detailed user's guide to the AMIDAS (A Model-Independent Data Analysis System) package and website, which is developed for online simulations and data analyses for direct Dark Matter detection experiments and phenomenology. Recently, the whole AMIDAS package and website system has been upgraded to the second phase: AMIDAS-II, for including the new developed Bayesian analysis technique. AMIDAS has the ability to do full Monte Carlo simulations as well as to analyze real/pseudo data sets either generated by another event generating programs or recorded in direct DM detection experiments. Moreover, the AMIDAS-II package can include several "user-defined" functions into the main code: the (fitting) one-dimensional WIMP velocity distribution function, the nuclear form factors for spin-independent and spin-dependent cross sections, artificial/experimental background spectrum for both of simulation and data analysis procedures, as well as different distribution functions needed in Bayesian analyses.

In this paper, we extended our earlier work on the reconstruction of the (time-averaged) one-dimensional velocity distribution of Galactic Weakly Interacting Massive Particles (WIMPs) and introduce the Bayesian fitting procedure to the theoretically predicted velocity distribution functions. In this reconstruction process, the (rough) velocity distribution reconstructed by using raw data from direct Dark Matter detection experiments directly, i.e. measured recoil energies, with one or more different target materials, has been used as "reconstructed-input" information. By assuming a fitting velocity distribution function and scanning the parameter space based on the Bayesian analysis, the astronomical characteristic parameters, e.g. the Solar and Earth's Galactic velocities, will be pinned down as the output results. Our Monte-Carlo simulations show that this Bayesian scanning procedure could reconstruct the true (input) WIMP velocity distribution function pretty precisely with negligible systematic deviations of the reconstructed characteristic Solar and Earth's velocities and 1 sigma statistical uncertainties of <~ 20 km/s. Moreover, for the use of an improper fitting velocity distribution function, our reconstruction process could still offer useful information about the shape of the velocity distribution. In addition, by comparing these estimates to theoretical predictions, one could distinguish different (basic) functional forms of the theoretically predicted one-dimensional WIMP velocity distribution function with 2 sigma to 4 sigma confidence levels.

In this script handbook, we collect the basic (and partially upgraded) PHP scripts used for building the AMIDAS website (http://pisrv0.pit.physik.uni-tuebingen.de/darkmatter/amidas/), an online interactive simulation/data analysis system for direct Dark Matter detection experiments and phenomenology. In this (1.73205) version, we add more materials and improve the scripts for offering a more convenient, comfortable and user-friendly environment on interactive scientific computing websites. Some basic, often used commands of (X)HTML, CSS, JavaScript, HTML DOM, and PHP are also given in Appendix. Online demonstrations and downloadable template scripts are given on http://www.tir.tw/iswb/ss2012/templates.php.

In this paper, we introduce model-independent data analysis procedures for identifying inelastic WIMP-nucleus scattering as well as for reconstructing the mass and the mass splitting of inelastic WIMPs simultaneously and separately. Our simulations show that, with O(50) observed WIMP signals from one experiment, one could already distinguish the inelastic WIMP scattering scenarios from the elastic one. By combining two or more data sets with positive signals, the WIMP mass and the mass splitting could even be reconstructed with statistical uncertainties of less than a factor of two.

In this talk, I presented effects of small, but non-negligible unrejected background events on the determinations of WIMP couplings/cross sections.

In this talk I have presented the data analysis results of extracting properties of halo WIMPs: the mass and the (ratios between the) spin-independent and spin-dependent couplings/cross sections on nucleons by the AMIDAS website by taking into account possible unrejected background events in the analyzed data sets. Although non-standard astronomical setup has been used to generate pseudodata sets for our analyses, it has been found that, without prior information/assumption about the local density and velocity distribution of halo Dark Matter, these WIMP properties have been reconstructed with ~ 2% to <~ 30% deviations from the input values.

In this talk I have presented the data analysis results of extracting properties of halo WIMPs: the mass and the (ratios between the) spin-independent and spin-dependent couplings/cross sections on nucleons by the AMIDAS website. Although non-standard astronomical setup has been used to generate pseudodata sets for our analyses, it has been found that, without prior information/assumption about the local density and velocity distribution of halo Dark Matter, these WIMP properties have been reconstructed with ~5% to < ~40% deviations from the input values.

In our work on the development of model-independent data analysis methods for determining ratios between different couplings/cross sections of Weakly Interacting Massive Particles (WIMPs) by using measured recoil energies from direct Dark Matter detection experiments directly, it was assumed that the analyzed data sets are background-free, i.e., all events are WIMP signals. In this article, as a more realistic study, we take into account a fraction of possible residue background events, which pass all discrimination criteria and then mix with other real WIMP-induced events in our data sets. Our simulations show that, assuming that the spin-dependent (SD) WIMP-nucleus interaction dominates over the spin-independent (SI) one, the maximal acceptable fraction of residue background events in the analyzed data sets for determining the ratio of the SD WIMP coupling on neutrons to that on protons is ~ 20% - 40%; whereas considering a general combination of the SI and SD WIMP interactions, the maximal acceptable background ratio for determining the ratio between two SD WIMP couplings as well as the ratios of the SD cross section on protons (neutrons) to the SI one is ~ 10% - 20%. Moreover, by considering different forms of background spectrum, we find that only background events in the lowest energy ranges could affect the reconstructions (significantly); those in high energy ranges would almost not change the reconstructed ratios or only very slightly.

Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for Dark Matter. For understanding the nature of WIMPs and identifying them among new particles produced at colliders (hopefully in the near future), determinations of their mass and couplings on nucleons from direct Dark Matter detection experiments are essential. Based on our model-independent method for determining the WIMP mass from experimental data, I present a way to also estimate the spin-independent (SI) WIMP-nucleon coupling by using measured recoil energies directly. This method is independent of the velocity distribution of halo WIMPs as well as (practically) of the as yet unknown WIMP mass. In a background-free environment, for a WIMP mass of ~ 100 GeV the SI WIMP-nucleon coupling could in principle be estimated with an uncertainty of ~ 15% by using 2 (or 3) x 50 events from experiments.

In our work on the development of a model-independent data analysis method for estimating the spin-independent (SI) scalar coupling of Weakly Interacting Massive Particles (WIMPs) on nucleons by using measured recoil energies from direct Dark Matter detection experiments directly, it was assumed that the analyzed data sets are background-free, i.e., all events are WIMP signals. In this article, as a more realistic study, we take into account a fraction of possible residue background events, which pass all discrimination criteria and then mix with other real WIMP-induced events in our data sets. Our simulations show that, for the estimation of the SI WIMP-nucleon coupling, the maximal acceptable fraction of residue background events in the analyzed data set of O(50) total events is ~ 10% - 20%. For a WIMP mass of 100 GeV and 20% residue background events, the systematic deviation of the reconstructed SI WIMP coupling (with a reconstructed WIMP mass) would in principle be ~ +13% with a statistical uncertainty of ~ ^{+21%}_{-14%} (~ -3.3%^{+18%}_{-13%} for background-free data sets).

Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for Dark Matter. So far the usual procedure for constraining the WIMP-nucleon cross sections in direct Dark Matter detection experiments have been to fit the predicted event rate based on some model(s) of the Galactic halo and of WIMPs to experimental data. One has to assume whether the spin-independent (SI) or the spin-dependent (SD) WIMP-nucleus interaction dominates, and results of such data analyses are also expressed as functions of the as yet unknown WIMP mass. In this article, I introduce methods for extracting information on the WIMP-nucleon cross sections by considering a general combination of the SI and SD interactions. Neither prior knowledge about the local density and the velocity distribution of halo WIMPs nor about their mass is needed. Assuming that an exponential-like shape of the recoil spectrum is confirmed from experimental data, the required information are only the measured recoil energies (in low energy ranges) and the number of events in the first energy bin from two or more experiments.

We reexamine the model-independent data analysis methods for extracting properties of Weakly Interacting Massive Particles (WIMPs) by using data (measured recoil energies) from direct Dark Matter detection experiments directly and, as a more realistic study, consider a small fraction of residue background events, which pass all discrimination criteria and then mix with other real WIMP-induced signals in the analyzed data sets. In this talk, the effects of residue backgrounds on the determination of the WIMP mass as well as the spin-independent WIMP coupling on nucleons will be discussed.

We reexamine the model-independent data analysis methods for extracting properties of Weakly Interacting Massive Particles (WIMPs) by using data (measured recoil energies) from direct Dark Matter detection experiments directly and, as a more realistic study, consider a small fraction of residue background events, which pass all discrimination criteria and then mix with other real WIMP-induced signals in the analyzed data sets. In this talk, the effects of residue backgrounds on the determination of the mass of halo Dark Matter particle as well as on the reconstruction of its one-dimensional velocity distribution function will be discussed.

In the earlier work on the development of a model-independent data analysis method for determining the mass of Weakly Interacting Massive Particles (WIMPs) by using measured recoil energies from direct Dark Matter detection experiments directly, it was assumed that the analyzed data sets are background-free, i.e., all events are WIMP signals. In this article, as a more realistic study, we take into account a fraction of possible residue background events, which pass all discrimination criteria and then mix with other real WIMP-induced events in our data sets. Our simulations show that, for the determination of the WIMP mass, the maximal acceptable fraction of residue background events in the analyzed data sets of O(50) total events is ~20%, for background windows of the entire experimental possible energy ranges, or in low energy ranges; while, for background windows in relatively higher energy ranges, this maximal acceptable fraction of residue background events can not be larger than ~10%. For a WIMP mass of 100 GeV with 20% background events in the windows of the entire experimental possible energy ranges, the reconstructed WIMP mass and the 1-sigma statistical uncertainty are ~97 GeV^{+61%}_{-35%} (~94 GeV^{+55%}_{-33%} for background-free data sets).

In our earlier work on the development of a model-independent data analysis method for reconstructing the (moments of the) time-averaged one-dimensional velocity distribution function of Weakly Interacting Massive Particles (WIMPs) by using measured recoil energies from direct Dark Matter detection experiments directly, it was assumed that the analyzed data sets are background-free, i.e., all events are WIMP signals. In this article, as a more realistic study, we take into account a fraction of possible residue background events, which pass all discrimination criteria and then mix with other real WIMP-induced events in our data sets. Our simulations show that, for the reconstruction of the one-dimensional WIMP velocity distribution, the maximal acceptable fraction of residue background events in the analyzed data set(s) of O(500) total events is ~10% - 20%. For a WIMP mass of 50 GeV with a negligible uncertainty and 20% residue background events, the deviation of the reconstructed velocity distribution would in principle be ~7.5% with a statistical uncertainty of ~18% (~19% for a background-free data set).

In this article I review model-independent procedures for extracting properties of Weakly Interacting Massive Particles (WIMPs) from direct Dark Matter detection experiments. Neither prior knowledge about the velocity distribution function of halo Dark Matter particles nor about their mass or cross sections on target nucleus is needed. The unique required information is measured recoil energies from experiments with different detector materials.

In this article I review two data analysis methods for determining the mass (and eventually the spin-independent cross section on nucleons) of Weakly Interacting Massive Particles with positive signals from direct Dark Matter detection experiments: a maximum likelihood analysis with only one experiment and a model-independent method requiring at least two experiments. Uncertainties and caveats of these methods will also be discussed.

The AMIDAS website has been established as an online interactive tool for running simulations and analyzing data in direct Dark Matter detection experiments. At the first phase of the website building, only some commonly used WIMP velocity distribution functions and elastic nuclear form factors have been involved in the AMIDAS code. In order to let the options for velocity distribution as well as for nuclear form factors be more flexible, we have extended the AMIDAS code to be able to include user-uploaded files with their own functions. In this article, I describe the preparation of files of user-defined functions onto the AMIDAS website. Some examples will also be given.

Following our long-term work on development of model-independent data analysis methods for reconstructing the one-dimensional velocity distribution function of halo WIMPs as well as for determining their mass and couplings on nucleons by using data from direct Dark Matter detection experiments directly, we combined the simulation programs to a compact system: AMIDAS (A Model-Independent Data Analysis System). For users' convenience an online system has also been established at the same time. AMIDAS has the ability to do full Monte Carlo simulations, faster theoretical estimations, as well as to analyze (real) data sets recorded in direct detection experiments without modifying the source code. In this article, I give an overview of functions of the AMIDAS code based on the use of its website.

In this talk we present data analysis methods for reconstructing the mass and couplings of Weakly Interacting Massive Particles (WIMPs) by using directly future experimental data (i.e., measured recoil energies) from direct Dark Matter detection. These methods are independent of the model of Galactic halo as well as of WIMPs. The basic ideas of these methods and the feasibility and uncertainties of applying them to direct detection experiments with the next generation detectors will be discussed.

Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for Dark Matter. For understanding the properties of WIMPs and identifying them among new particles produced at colliders (hopefully in the near future), determinations of their mass and their couplings on nucleons from direct Dark Matter detection experiments are essential. Based on our method for determining the WIMP mass model-independently from experimental data, we present a way to also estimate the spin-independent (SI) WIMP-nucleon coupling by using measured recoil energies directly. This method isindependent of the as yet unknown velocity distribution of halo WIMPs. In spite of the uncertainty of the local WIMP density (of a factor of ~ 2), at least an upper limit on the SI WIMP-nucleon coupling could be given, once two (or more) experiments with different target nuclei obtain positive signals. In a background-free environment, for a WIMP mass of 100 GeV its SI coupling on nucleons could in principle be estimated with a statistical error of only ~ 15% with just 50 events from each experiment.

Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for Dark Matter. We developed a model-independent method for determining the WIMP mass by using data (i.e., measured recoil energies) of direct detection experiments. Our method is independent of the as yet unknown WIMP density near the Earth, of the form of the WIMP velocity distribution, as well as of the WIMP-nucleus cross section. It requires however positive signals from at least two detectors with different target nuclei. At the first phase of this work we found a systematic deviation of the reconstructed WIMP mass from the real one for heavy WIMPs. Now we improved this method so that this deviation can be strongly reduced for even very high WIMP mass. The statistical error of the reconstructed mass has also been reduced. In a background-free evironment, a WIMP mass of ~ 50 GeV could in principle be determined with an error of ~ 35% with only 2 times 50 events.

Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for Dark Matter. We develop a model-independent method for determining the mass $m_\chi$ of the WIMP by using data (i.e., measured recoil energies) of direct detection experiments. Our method is independent of the as yet unknown WIMP density near the Earth, of the form of the WIMP velocity distribution, as well as of the WIMP-nucleus cross section. However, it requires positive signals from at least two detectors with different target nuclei. In a background-free environment, $m_\chi \sim 50$ GeV could in principle be determined with an error of $\sim 35%$ with only $2 \times 50$ events; in practice upper and lower limits on the recoil energy of signal events, imposed to reduce backgrounds, can increase the error. The method also loses precision if $m_\chi$ significantly exceeds the mass of the heaviest target nucleus used.

Weakly interacting massive particles (WIMPs) are one of the leading candidates for Dark Matter. So far we can use direct Dark Matter detection to estimate the mass of halo WIMPs only by fitting predicted recoil spectra to future experimental data. Here we develop a model-independent method for determining the WIMP mass by using experimental data directly. This method is independent of the as yet unknown WIMP density near the Earth as well as of the WIMP-nuclear cross section and can be used to extract information about WIMP mass with O(50) events.

Weakly Interacting Massive Particles (WIMPs) are one of the leading candidates for Dark Matter. Currently, the most promising method to detect WIMPs is the direct detection of the recoil energy deposited in a low-background laboratory detector due to elastic WIMP-nucleus scattering. So far the usual procedure has been to predict the event rate of direct detection of WIMPs based on some model(s) of the Galactic halo from cosmology and of WIMPs from elementary particle physics. The aim of this work is to invert this process. In this thesis I present methods which allow to reconstruct (the moments of) the WIMP velocity distribution function as well as to determine the WIMP mass from the recoil energy spectrum as well as from experimental data directly. The reconstruction of the velocity distribution function has been further extended to take into account the annual modulation of the event rate. Moreover, the reconstruction of the amplitude of the annual modulation of the velocity distribution and an alternative, better way for confirming the annual modulation of the event rate have been discussed. On the other hand, the determination of the WIMP mass by combining two (or more) experiments with different detector materials has been developed. All formulae and expressions given here are not only independent of the model of Galactic halo but also of that of WIMPs. This means that we need neither the as yet unknown WIMP density near the Earth nor the WIMP-nucleus cross section. The only information which we need is the measured recoil energies and their measuring times.

Weakly interacting massive particles (WIMPs) are one of the leading candidates for dark matter. Currently, the most promising method to detect many different WIMP candidates is the direct detection of the recoil energy deposited in a low-background laboratory detector due to elastic WIMP-nucleus scattering. So far the usual procedure has been to predict the event rate of direct detection of WIMPs based on some model(s) of the galactic halo. The aim of our work is to invert this process. That is, we study what future direct detection experiment can teach us about the WIMP halo. As the first step we consider a time-averaged recoil spectrum, assuming that no directional information exists. We develop a method to construct the (time-averaged) one-dimensional velocity distribution function from this spectrum. Moments of this function, such as the mean velocity and velocity dispersion of WIMPs, can also be obtained directly from the recoil spectrum. The only input needed in addition to a measured recoil spectrum is the mass of the WIMP; no information about the scattering cross section or WIMP density is required.