Using the Fenchel-Eggleston theorem for convex hulls (an extension of the Caratheodory theorem), we prove that any likelihood can be maximized by either a dark matter 1- speed distribution $F(v)$ in Earth's frame or 2- Galactic velocity distribution $f^{\rm gal}(\vec{u})$, consisting of a sum of delta functions. The former case applies only to time-averaged rate measurements and the maximum number of delta functions is $({\mathcal N}-1)$, where ${\mathcal N}$ is the total number of data entries. The second case applies to any harmonic expansion coefficient of the time-dependent rate and the maximum number of terms is ${\mathcal N}$. Using time-averaged rates, the aforementioned form of $F(v)$ results in a piecewise constant unmodulated halo function $\tilde\eta^0_{BF}(v_{\rm min})$ (which is an integral of the speed distribution) with at most $({\mathcal N}-1)$ downward steps. The authors had previously proven this result for likelihoods comprised of at least one extended likelihood, and found the best-fit halo function to be unique. This uniqueness, however, cannot be guaranteed in the more general analysis applied to arbitrary likelihoods. Thus we introduce a method for determining whether there exists a unique best-fit halo function, and provide a procedure for constructing either a pointwise confidence band, if the best-fit halo function is unique, or a degeneracy band, if it is not. Using measurements of modulation amplitudes, the aforementioned form of $f^{\rm gal}(\vec{u})$, which is a sum of Galactic streams, yields a periodic time-dependent halo function $\tilde\eta_{BF}(v_{\rm min}, t)$ which at any fixed time is a piecewise constant function of $v_{\rm min}$ with at most ${\mathcal N}$ downward steps. In this case, we explain how to construct pointwise confidence and degeneracy bands from the time-averaged halo function. Finally, we show that requiring an isotropic ...

We present two different halo-independent methods to assess the compatibility of several direct dark matter detection data sets for a given dark matter model using a global likelihood consisting of at least one extended likelihood and an arbitrary number of Gaussian or Poisson likelihoods. In the first method we find the global best fit halo function (we prove that it is a unique piecewise constant function with a number of down steps smaller than or equal to a maximum number that we compute) and construct a two-sided pointwise confidence band at any desired confidence level, which can then be compared with those derived from the extended likelihood alone to assess the joint compatibility of the data. In the second method we define a "constrained parameter goodness-of-fit" test statistic, whose $p$-value we then use to define a "plausibility region" (e.g. where $p \geq 10\%$). For any halo function not entirely contained within the plausibility region, the level of compatibility of the data is very low (e.g. $p < 10 \%$). We illustrate these methods by applying them to CDMS-II-Si and SuperCDMS data, assuming dark matter particles with elastic spin-independent isospin-conserving interactions or exothermic spin-independent isospin-violating interactions.

We present comparisons of direct dark matter (DM) detection data for light WIMPs with exothermic scattering with nuclei (exoDM), both assuming the Standard Halo Model (SHM) and in a halo model - independent manner. Exothermic interactions favor light targets, thus reducing the importance of upper limits derived from xenon targets, the most restrictive of which is at present the LUX limit. In our SHM analysis the CDMS-II-Si and CoGeNT regions become allowed by these bounds, however the recent SuperCDMS limit rejects both regions for exoDM with isospin-conserving couplings. An isospin-violating coupling of the exoDM, in particular one with a neutron to proton coupling ratio of $-0.8$ (which we call "Ge-phobic"), maximally reduces the DM coupling to germanium and allows the CDMS-II-Si region to become compatible with all bounds. This is also clearly shown in our halo-independent analysis.

We present comparisons of direct detection data for "light WIMPs" with an anapole moment interaction (ADM) and a magnetic dipole moment interaction (MDM), both assuming the Standard Halo Model (SHM) for the dark halo of our galaxy and in a halo-independent manner. In the SHM analysis we find that a combination of the 90% CL LUX and CDMSlite limits or the new 90% CL SuperCDMS limit by itself exclude the parameter space regions allowed by DAMA, CoGeNT and CDMS-II-Si data for both ADM and MDM. In our halo-independent analysis the new LUX bound excludes the same potential signal regions as the previous XENON100 bound. Much of the remaining signal regions is now excluded by SuperCDMS, while the CDMSlite limit is much above them. The situation is of strong tension between the positive and negative search results both for ADM and MDM. We also clarify the confusion in the literature about the ADM scattering cross section.

We show that the relic abundance of the minority component of asymmetric dark matter can be very sensitive to the expansion rate of the Universe and the temperature of transition between a non-standard pre-Big Bang Nucleosynthesis cosmological phase and the standard radiation dominated phase, if chemical decoupling happens before this transition. In particular, because the annihilation cross section of asymmetric dark matter is typically larger than that of symmetric dark matter in the standard cosmology, the decrease in relic density of the minority component in non-standard cosmologies with respect to the majority component may be compensated by the increase in annihilation cross section, so that the annihilation rate at present of asymmetric dark matter, contrary to general belief, could be larger than that of symmetric dark matter in the standard cosmology.

We present a halo-independent analysis of direct detection data on "light WIMPs," i.e. weakly interacting massive particles with mass close to or below 10 GeV/c^2. We include new results from silicon CDMS detectors (bounds and excess events), the latest CoGeNT acceptances, and recent measurements of low sodium quenching factors in NaI crystals. We focus on light WIMPs with spin-independent isospin-conserving and isospin-violating interactions with nucleons. For these dark matter candidates we find that a low quenching factor would make the DAMA modulation incompatible with a reasonable escape velocity for the dark matter halo, and that the tension among experimental data tightens in both the isospin-conserving and isospin-violating scenarios. We also find that a new although milder tension appears between the CoGeNT and DAMA annual modulations on one side and the silicon excess events on the other, in that it seems difficult to interpret them as the modulated and unmodulated aspects of the same WIMP dark matter signal.

We derive constraints on parameters of generic dark matter candidates by comparing theoretical predictions with the gamma-ray emission observed by the Fermi-LAT from the region around the Galactic Center. Our analysis is conservative since it simply requires that the expected dark matter signal does not exceed the observed emission. The constraints obtained in the likely case that the collapse of baryons to the Galactic Center is accompanied by the contraction of the dark matter are strong. In particular, we find that for bb and \tau+\tau- or W+W- dark matter annihilation channels, the upper limits on the annihilation cross section imply that the thermal cross section is excluded for a Weakly Interacting Massive Particle (WIMP) mass smaller than about 700 and 500 GeV, respectively. For the \mu+ \mu- channel, where the effect of the inverse Compton scattering is important, depending on models of the Galactic magnetic field the exclusion of the thermal cross-section is for a WIMP mass smaller than about 150 to 400 GeV. The upper limits on the annihilation cross section of dark matter particles obtained are two orders of magnitude stronger than without contraction. In the latter case our results are compatible with the upper limits from the Galactic halo analysis reported by the Fermi-LAT collaboration for the case in which the same conservative approach without modeling of the astrophysical background is employed.

We study how the combined observation of dark matter in various direct detection experiments can be used to determine the phenomenological properties of WIMP dark matter: mass, spin-dependent (SD) and spin-independent (SI) scattering cross section off nucleons. A convenient choice of target materials, including nuclei that couple to dark matter particles through a significantly different ratio of SD vs SI interactions, could break the degeneracies in the determination of those parameters that a single experiment cannot discriminate. In this work we investigate different targets that can be used as scintillating bolometers and could provide complementary information to germanium and xenon detectors. We observe that Al2O3 and LiF bolometers could allow a good reconstruction of the DM properties over regions of the parameter space with a SD scattering cross section as small as 10^(-5) pb and a SI cross section as small as 5x10^(-10) pb for a 50 GeV WIMP. In the case of a CaWO4 bolometer the area in which full complementarity is obtained is smaller but we show that it can be used to determine the WIMP mass and its SI cross section. For each target we study the required exposure and background.

We study the effect that uncertainties in the nuclear spin-dependent structure functions have in the determination of the dark matter (DM) parameters in a direct detection experiment. We show that different nuclear models that describe the spin-dependent structure function of specific target nuclei can lead to variations in the reconstructed values of the DM mass and scattering cross-section. We propose a parametrization of the spin structure functions that allows us to treat these uncertainties as variations of three parameters, with a central value and deviation that depend on the specific nucleus. The method is illustrated for germanium and xenon detectors with an exposure of 300 kg yr, assuming a hypothetical detection of DM and studying a series of benchmark points for the DM properties. We find that the effect of these uncertainties can be similar in amplitude to that of astrophysical uncertainties, especially in those cases where the spin-dependent contribution to the elastic scattering cross-section is sizable.

Very light right-handed (RH) sneutrinos in the Next-to-Minimal Supersymmetric Standard Model can be viable candidates for cold dark matter. We investigate the prospects for their direct detection, addressing their compatibility with the recent signal observed by the CoGeNT detector, and study the implications for Higgs phenomenology. We find that in order to reproduce the correct relic abundance very light RH sneutrinos can annihilate into either a fermion-antifermion pair, very light pseudoscalar Higgses or RH neutrinos. If the main annihilation channel is into fermions, we point out that RH sneutrinos could naturally account for the CoGeNT signal. Furthermore, the lightest Higgs has a very large invisible decay width, and in some cases the second-lightest Higgs too. On the other hand, if the RH sneutrino annihilates mostly into pseudoscalars or RH neutrinos the predictions for direct detection are below the current experimental sensitivities and satisfy the constraints set by CDMS and XENON. We also calculate the gamma ray flux from RH sneutrino annihilation in the Galactic centre, including as an interesting new possibility RH neutrinos in the final state. These are produced through a resonance with the Higgs and the resulting flux can exhibit a significant Breit-Wigner enhancement.

The first stars to form in the Universe may be powered by the annihilation of weakly interacting dark matter particles. These so-called dark stars, if observed, may give us a clue about the nature of dark matter. Here we examine which models for particle dark matter satisfy the conditions for the formation of dark stars. We find that in general models with thermal dark matter lead to the formation of dark stars, with few notable exceptions: heavy neutralinos in the presence of coannihilations, annihilations that are resonant at dark matter freeze-out but not in dark stars, some models of neutrinophilic dark matter annihilating into neutrinos only and lighter than about 50 GeV. In particular, we find that a thermal DM candidate in standard Cosmology always forms a dark star as long as its mass is heavier than about 50 GeV and the thermal average of its annihilation cross section is the same at the decoupling temperature and during the dark star formation, as for instance in the case of an annihilation cross section with a non-vanishing s-wave contribution.

A TeV scale decaying dark matter chiral multiplet N is introduced in addition to the minimal supersymmetric standard model(MSSM). For a calculable abundance of N, we introduce heavy axino decaying to N and MSSM particles including the lightest supersymmetric particle(LSP). In the scenario that heavy axino, once dominating the energy density of the universe, decays after the LSP decouples, it is possible to estimate the relative cosmic abundances of N and the LSP. Dimension 6 interactions allow the lifetime of the fermionic or the bosonic superpartner of N in the 10^{27} s range to be compatible with the recent astrophysical bounds. A diagrammatic strategy obtaining a suppression factor 1/M^2 is also given.

Recently, there are two hints arising from physics beyond the standard model. One is a possible energy loss mechanism due to emission of very weakly interacting light particles from white dwarf stars, with a coupling strength ~ 0.7x10^{-13}, and another is the high energy positrons observed by the PAMELA satellite experiment. We construct a supersymmetric flipped-SU(5) model, SU(5)xU(1)_X with appropriate additional symmetries, [U(1)_H]_{gauge}x[U(1)_RxU(1)_\Gamma]_{global}xZ_2, such that these are explained by a very light electrophilic axion of mass 0.5 meV from the spontaneously broken U(1)_\Gamma and two component cold dark matters from Z_2 parity. We show that in the flipped-SU(5) there exists a basic mechanism for allowing excess positrons through the charged SU(2) singlet leptons, but not allowing anti-proton excess due to the absence of the SU(2) singlet quarks. We show the discovery potential of the charged SU(2) singlet E at the LHC experiments by observing the electron and positron spectrum. With these symmetries, we also comment on the mass hierarchy between the top and bottom quarks.

We present the dark matter (DM) extension (by N) of the minimal supersymmetric standard model to give the recent trend of the high energy positron spectrum of the PAMELA/HEAT experiments. If the trend survives by future experiments, the MSSM needs to be extended. Here, we minimally extend the MSSM with one more DM component N together with a heavy lepton E, and introduce the coupling e_R E_R^c N_R. This coupling naturally appears in the flipped SU(5) GUT models. This N_{DM}MSSM contains the discrete symmetry Z_6, and for some parameter ranges there result two DM components. For the MSSM fields, the conventional R-parity, which is a subgroup of Z_6, is preserved. We also present the needed parameter ranges of these additional particles.

We propose a possible explanation for the recently observed anomalous 511 keV line with a new "millicharged" fermion. This new fermion is light [${\cal O}({\rm MeV})$]. Nevertheless, it has never been observed by any collider experiments by virtue of its tiny electromagnetic charge $\epsilon e$. In particular, we constrain parameters of this millicharged particle if the 511 keV cosmic $\gamma$-ray emission from the galactic bulge is due to positron production from this new particle.

We improve the estimate of the axion CDM energy density by considering the new values of current quark masses, the QCD phase transition effect and a possible anharmonic effect.

We present a possible explanation of the recently observed 511 keV $\gamma$-ray anomaly with a new ``millicharged'' fermion. The new fermion is light (${\cal O}({\rm MeV})$) but has never been observed by any collider experiments mainly because of its tiny electromagnetic charge $\epsilon e$. We show that constraints from its relic density in the Universe and collider experiments allow a parameter range such that the 511 keV cosmic $\gamma$-ray emission from the galactic bulge may be due to positron production from this millicharged fermion.