Black holes formed in the early universe, prior to the formation of stars, can exist as dark matter and also contribute to the black hole merger events observed in gravitational waves. We set a new limit on the abundance of primordial black holes (PBHs) by considering interactions of PBHs with the interstellar medium, which result in the heating of gas. We examine generic heating mechanisms, including emission from the accretion disk, dynamical friction, and disk outflows. Using the data from the Leo T dwarf galaxy, we set a new cosmology-independent limit on the abundance of PBHs in the mass range $\mathcal{O}(1) M_{\odot}-10^7 M_{\odot}$.

A simple extension of the Standard Model consists of a scalar field that can potentially constitute the dark matter (DM). Significant attention has been devoted to probing light $\mathcal{O}(\lesssim 10~\rm{eV})$ scalar DM, with a multitude of experimental proposals based on condensed matter systems as well as novel materials. However, the previously overlooked effective in-medium mixing of light scalars with longitudinal plasmons can strongly modify the original sensitivity calculations of the direct detection experiments. We implement the in-medium effects for scalar DM detection, using thermal field theory techniques, and show that the reach of a large class of direct DM detection experiments searching for light scalars is significantly reduced. This development identifies setups based on Dirac materials and tunable plasma haloscopes as particularly promising for scalar DM detection. Further, we also show that scalars with significant boost with respect to halo DM, such as those produced in the Sun, decay of other particles or by cosmic rays, will not suffer from in-medium suppression. Hence, multi-tonne direct DM detection experiments, such as those based on xenon or argon, also constitute a favorable target. We also discuss scalar mediated DM scattering.

In this addendum to the article JCAP 1912 (2019) no.12, 047 (arXiv:1909.13328) on the cosmological dependence of non-resonantly produced sterile neutrinos we discuss, using an analytic treatment, the parameter regions of large active-sterile neutrino mixing angles where sterile neutrinos can approach thermalization. We show that these additional considerations affect only large active-sterile neutrino mixing already rejected by different limits. Hence, the allowed sterile neutrino parameter regions are unaffected.

Dark photons (DPs) produced in the early Universe are well motivated dark matter (DM) candidates. We show that the recently proposed tunable plasma haloscopes are particularly advantageous for DP searches. While in-medium effects suppress the DP signal in conventional searches, plasma haloscopes make use of metamaterials that enable resonant absorption of the DP by matching its mass to a tunable plasma frequency and thus enable efficient plasmon production. Using thermal field theory, we confirm the in-medium DP absorption rate within the detector. This scheme allows to competitively explore a significant part of the DP DM parameter space in the DP mass range of $6-400$ $\mu$eV. If a signal is observed, the observation of a daily or annual modulation of the signal would be crucial to clearly identify the signal as due to DP DM and could shed light on the production mechanism.

We consider the cosmological effects of sterile neutrinos with the masses of $150- 450$ MeV. The decay of sterile neutrinos changes the thermal history of the Universe and affects the energy density of radiation at the recombination and the big bang nucleosynthesis (BBN) results. We derive severe constraints on the parameters of sterile neutrinos from the primordial abundances of helium-4 and deuterium. We also find that in a particular model the constraints can be considerably relaxed by assuming a large lepton asymmetry in the active neutrinos. In this case, the consistent parameters result in $N_{\mathrm{eff}} \simeq 3.2- 3.4$ and can alleviate the Hubble tension.

The detection of a sterile neutrino could constitute the first observation of a particle that could have been produced before Big-Bang Nucleosynthesis (BBN), and could provide information about the yet untested pre-BBN era. The cosmological evolution in this era could be drastically different than typically assumed in what constitutes the standard cosmology, as happens in a variety of motivated particle models. In this work we assess the sensitivity to different pre-BBN cosmologies in which entropy is conserved of 0.01 eV to 1 MeV mass sterile neutrinos produced in the early Universe via resonant active-sterile oscillations, which requires a large lepton asymmetry. We identify mass ranges where it is possible to have two populations of the same sterile neutrino, one with a colder and one with a hotter momentum spectra, which is in principle an observable effect. Furthermore, we show the regions in mass and mixing where fully resonant production (i.e. simultaneously coherent and adiabatic) can occur. We find that in several of the cosmologies we consider, including the standard one, for a lepton asymmetry larger than $\sim10^{-4}$ fully resonantly produced sterile neutrinos in the eV-mass range can evade all cosmological constraints.

We discuss how a laboratory detection of a sterile neutrino not only would constitute a fundamental discovery of a new particle, but could also provide an indication of the evolution of the Universe before Big-Bang Nucleosynthesis (BBN), a fundamental discovery in cosmology. These "visible" sterile neutrinos could be detected in experiments such as KATRIN/TRISTAN and HUNTER in the keV mass range and PTOLEMY, KATRIN and reactor neutrino experiments in the eV mass range. Standard assumptions are usually made to compute the relic abundance and momentum distribution of particles produced before the temperature of the Universe was 5 MeV, an epoch from which there are no observed remnants thus far. However, non-standard pre-BBN cosmologies based on other assumptions that are equally in agreement with all existing data can arise in some theoretical models. We revisit the production of 0.01 eV to 1 MeV sterile neutrinos via non-resonant active-sterile flavor oscillations in several pre-BBN cosmologies. We give general equations for models in which the expansion of the Universe is parametrized by its amplitude and temperature power and where entropy is conserved, which include kination and scalar tensor models as special cases.

A laboratory detection of a sterile neutrino could provide the first indication of the evolution of the Universe before Big-Bang Nucleosynthesis (BBN), an epoch yet untested. Such "visible" sterile neutrinos are observable in upcoming experiments such as KATRIN/TRISTAN and HUNTER in the keV mass range and PTOLEMY and others in the eV mass range. A set of standard assumptions is typically made about cosmology before the temperature of the Universe was 5 MeV. However, non-standard pre-BBN cosmologies based on alternative assumptions could arise in motivated theoretical models and are equally in agreement with all existing data. We revisit the production of sterile neutrinos of mass 0.01 eV to 1 MeV in two examples of such models: scalar-tensor and low reheating temperature pre-BBN cosmologies. In both of them, the putative 3.5 keV X-ray signal line corresponds to a sterile neutrino with a mixing large enough to be tested in upcoming laboratory experiments. Additionally, the cosmological/astrophysical upper limits on active-sterile neutrino mixings are significantly weaker than in the standard pre-BBN cosmology, which shows that these limits are not robust. For example, in the scalar-tensor case the potential signal regions implied by the LSND and MiniBooNE short-baseline as well as the DANSS and NEOS reactor experiments are entirely not bound by cosmological restrictions. Our work highlights that sterile neutrinos may constitute a sensitive probe of the pre-BBN epoch.

This is an idiosyncratic account of the main results presented at the 31st Rencontres de Blois, which took place from June 2nd to June 7th, 2019 in the Castle of Blois, France.

Local Universe observations find a value of the Hubble constant $H_0$ that is larger than the value inferred from the Cosmic Microwave Background and other early Universe measurements, assuming known physics and the $\Lambda$CDM cosmological model. We show that additional radiation in active neutrinos produced just before Big Bang Nucleosynthesis by an unstable sterile neutrino with mass $m_s=$ O(10) MeV can alleviate this discrepancy. The necessary masses and couplings of the sterile neutrino, assuming it mixes primarily with $\nu_{\tau}$ and/or $\nu_{\mu}$ neutrinos, are within reach of Super-Kamiokande as well as upcoming laboratory experiments such as NA62.

As dark matter (DM) direct detection experiments continue to improve their sensitivity they will inevitably encounter an irreducible background arising from coherent neutrino scattering. This so-called "neutrino floor" may significantly reduce the sensitivity of an experiment to DM-nuclei interactions, particularly if the recoil spectrum of the neutrino background is approximately degenerate with the DM signal. This occurs for the conventionally considered spin-independent (SI) or spin-dependent (SD) interactions. In such case, an increase in the experiment's exposure by multiple orders of magnitude may not yield any significant increase in sensitivity. The typically considered SI and SD interactions, however, do not adequately reflect the whole landscape of the well-motivated DM models, which includes other interactions. Since particle DM has not been detected yet in laboratories, it is essential to understand and maximize the detection capabilities for a broad variety of possible models and signatures. In this work we explore the impact of the background arising from various neutrino sources on the discovery potential of a DM signal for a large class of viable DM-nucleus interactions and several potential futuristic experimental settings, with different target elements. For some momentum suppressed cross sections, large DM particle masses and heavier targets, we find that there is no suppression of the discovery limits due to neutrino backgrounds. Further, we explicitly demonstrate that inelastic scattering, which could appear in models with multicomponent dark sectors, would help to lift the signal degeneracy associated with the neutrino floor. This study could assist with mapping out the optimal DM detection strategy for the next generation of experiments.

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 an updated halo-dependent and halo-independent analysis of viable light WIMP dark matter candidates which could account for the excess observed in CDMS-II-Si. We include recent constraints from LUX, PandaX-II, and PICO-60, as well as projected sensitivities for XENON1T, SuperCDMS SNOLAB, LZ, DARWIN, DarkSide-20k, and PICO-250, on candidates with spin-independent isospin conserving and isospin-violating interactions, and either elastic or exothermic scattering. We show that there exist dark matter candidates which can explain the CDMS-II-Si data and remain very marginally consistent with the null results of all current experiments, however such models are highly tuned, making a dark matter interpretation of CDMS-II-Si very unlikely. We find that these models can only be ruled out in the future by an experiment comparable to LZ or PICO-250.

Directional dark matter detection attempts to measure the direction of motion of nuclei recoiling after having interacted with dark matter particles in the halo of our Galaxy. Due to Earth's motion with respect to the Galaxy, the dark matter flux is concentrated around a preferential direction. An anisotropy in the recoil direction rate is expected as an unmistakable signature of dark matter. The average nuclear recoil direction is expected to coincide with the average direction of dark matter particles arriving to Earth. Here we point out that for a particular type of dark matter, inelastic exothermic dark matter, the mean recoil direction as well as a secondary feature, a ring of maximum recoil rate around the mean recoil direction, could instead be opposite to the average dark matter arrival direction. Thus, the detection of an average nuclear recoil direction opposite to the usually expected direction would constitute a spectacular experimental confirmation of this type of dark matter.

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.

Cosmological observations indicate that most of the matter in the Universe is Dark Matter. Dark Matter in the form of Weakly Interacting Massive Particles (WIMPs) can be detected directly, via its elastic scattering off target nuclei. Most current direct detection experiments only measure the energy of the recoiling nuclei. However, directional detection experiments are sensitive to the direction of the nuclear recoil as well. Due to the Sun's motion with respect to the Galactic rest frame, the directional recoil rate has a dipole feature, peaking around the direction of the Solar motion. This provides a powerful tool for demonstrating the Galactic origin of nuclear recoils and hence unambiguously detecting Dark Matter. Furthermore, the directional recoil distribution depends on the WIMP mass, scattering cross section and local velocity distribution. Therefore, with a large number of recoil events it will be possible to study the physics of Dark Matter in terms of particle and astrophysical properties. We review the potential of directional detectors for detecting and characterizing WIMPs.

Earth's rotation about the Sun produces an annual modulation in the expected scattering rate at direct dark matter detection experiments. The annual modulation as a function of the recoil energy $E_\text{R}$ imparted by the dark matter particle to a target nucleus is expected to vary depending on the detector material. However, for most interactions a change of variables from $E_\text{R}$ to $v_\text{min}$, the minimum speed a dark matter particle must have to impart a fixed $E_\text{R}$ to a target nucleus, produces an annual modulation independent of the target element. We recently showed that if the dark matter-nucleus cross section contains a non-factorizable target and dark matter velocity dependence, the annual modulation as a function of $v_\text{min}$ can be target dependent. Here we examine more extensively the necessary conditions for target-dependent modulation, its observability in present-day experiments, and the extent to which putative signals could identify a dark matter-nucleus differential cross section with a non-factorizable dependence on the dark matter velocity.

We study how gravitational focusing (GF) of dark matter by the Sun affects the annual and biannual modulation of the expected signal in non-directional direct dark matter searches, in the presence of dark matter substructure in the local dark halo. We consider the Sagittarius stream and a possible dark disk, and show that GF suppresses some, but not all, of the distinguishing features that would characterize substructure of the dark halo were GF neglected.

Due to Earth's revolution around the Sun, the expected scattering rate in direct dark matter searches is annually modulated. This modulation is expected to differ between experiments when given as a function of recoil energy $E_\text{R}$, e.g. due to the gravitational focusing effect of the Sun. A better variable to compare results among experiments employing different targets is the minimum speed $v_\text{min}$ a dark matter particle must have to impart a recoil energy $E_\text{R}$ to a target nucleus. It is widely believed that the modulation expressed as a function of $v_\text{min}$ is common to all experiments, irrespective of the dark matter distribution. We point out that the annual modulation as a function of $v_\text{min}$, and in particular the times at which the rate is maximum and minimum, could be very different depending on the detector material. This would be an indication of a scattering cross section with non-factorizable velocity and target material dependence. Observing an annual modulation with at least two different target elements would be necessary to identify this type of cross section.

We investigate the allowed range of reheating temperature values in light of the Planck 2015 results and the recent joint analysis of Cosmic Microwave Background (CMB) data from the BICEP2/Keck Array and Planck experiments, using monomial and binomial inflationary potentials. While the well studied $\phi^2$ inflationary potential is no longer favored by current CMB data, as well as $\phi^p$ with $p>2$, a $\phi^1$ potential and canonical reheating ($w_{re}=0$) provide a good fit to the CMB measurements. In this last case, we find that the Planck 2015 $68\%$ confidence limit upper bound on the spectral index, $n_s$, implies an upper bound on the reheating temperature of $T_{re}\lesssim 6\times 10^{10}\,{\rm GeV}$, and excludes instantaneous reheating. The low reheating temperatures allowed by this model open the possiblity that dark matter could be produced during the reheating period instead of when the Universe is radiation dominated, which could lead to very different predictions for the relic density and momentum distribution of WIMPs, sterile neutrinos, and axions. We also study binomial inflationary potentials and show the effects of a small departure from a $\phi^1$ potential. We find that as a subdominant $\phi^2$ term in the potential increases, first instantaneous reheating becomes allowed, and then the lowest possible reheating temperature of $T_{re}=4\,{\rm MeV}$ is excluded by the Planck 2015 $68\%$ confidence limit.

We study the effects of kination and scalar-tensor pre-Big Bang Nucleosynthesis cosmologies on the non-resonant production of sterile neutrinos. We show that if the peak of the production rate of sterile neutrinos occurs during the non-standard cosmological phase, the relic number density of sterile neutrinos could be reduced with respect to the number expected in the standard cosmology. Consequently, current bounds on active-sterile neutrino mixing derived from the relic energy density of sterile neutrinos could be greatly relaxed. In particular, we show that the sterile neutrinos which could explain the anomalies found in short-baseline neutrino experiments are compatible with recent joint Planck upper limits on their contribution to the energy density of the Universe in a scalar-tensor or a low-reheating temperature pre-Big Bang Nucleosynthesis cosmology.

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 extend the halo-independent method of Fox, Liu, and Weiner to include energy resolution and efficiency with arbitrary energy dependence, making it more suitable for experiments to use in presenting their results. Then we compare measurements and upper limits on the direct detection of low mass ($\sim10$ GeV) weakly interacting massive particles with spin-independent interactions, including the upper limit on the annual modulation amplitude from the CDMS collaboration. We find that isospin-symmetric couplings are severely constrained both by XENON100 and CDMS bounds, and that isospin-violating couplings are still possible at the lowest energies, while the tension of the higher energy CoGeNT bins with the CDMS modulation constraint remains. We find the CRESST II signal is not compatible with the modulation signals of DAMA and CoGeNT.

We review fits to "light WIMPs" since the region was first mentioned relative to the DAMA collaboration data in 2003 to the present, analyzing the compatibility of potential signals and bounds in this region. We include dark halo independent data comparisons.

The motion of the Earth around the Sun causes an annual change in the magnitude and direction of the arrival velocity of dark matter particles on Earth, in a way analogous to aberration of stellar light. In directional detectors, aberration of weakly interacting massive particles (WIMPs) modulates the pattern of nuclear recoil directions in a way that depends on the orbital velocity of the Earth and the local galactic distribution of WIMP velocities. Knowing the former, WIMP aberration can give information on the latter, besides being a curious way of confirming the revolution of the Earth and the extraterrestrial provenance of WIMPs. While observing the full aberration pattern requires extremely large exposures, we claim that the annual variation of the mean recoil direction or of the event counts over specific solid angles may be detectable with moderately large exposures. For example, integrated counts over Galactic hemispheres separated by planes perpendicular to Earth's orbit would modulate annually, resulting in Galactic Hemisphere Annual Modulations (GHAM) with amplitudes larger than the usual non-directional annual modulation.

We discuss a novel dark matter signature relevant for directional detection of Weakly Interacting Massive Particles (WIMPs). For heavy enough WIMPs and low enough recoil energies, the maximum of the recoil rate is not in the direction of the average WIMP arrival direction but in a ring around it at an angular radius that increases with the WIMP mass and can approach 90 degrees at very low energies. The ring is easier to detect for smaller WIMP velocity dispersion and larger average WIMP velocities relative to the detector. In principle the ring could be used as an additional indication of the WIMP mass range.

The channeling of the recoiling nucleus in crystalline detectors after a WIMP collision would produce a larger scintillation or ionization signal in direct detection experiments than otherwise expected. I present estimates of channeling fractions obtained using analytic models developed from the 1960's onwards to describe channeling and blocking effects. We find the fractions to be too small to affect the fits to potential WIMP candidates. I also examine the possibility of detecting a daily modulation of the dark matter signal due to channeling.

The surface detector array of the Pierre Auger Observatory consists of 1600 water-Cherenkov detectors, for the study of extensive air showers (EAS) generated by ultra-high-energy cosmic rays. We describe the trigger hierarchy, from the identification of candidate showers at the level of a single detector, amongst a large background (mainly random single cosmic ray muons), up to the selection of real events and the rejection of random coincidences. Such trigger makes the surface detector array fully efficient for the detection of EAS with energy above $3\times 10^{18}$ eV, for all zenith angles between 0$^\circ$ and 60$^\circ$, independently of the position of the impact point and of the mass of the primary particle. In these range of energies and angles, the exposure of the surface array can be determined purely on the basis of the geometrical acceptance.

The dip model assumes that the ultra-high energy cosmic rays (UHECRs) above 10$^{18}$ eV consist exclusively of protons and is consistent with the spectrum and composition measure by HiRes. Here we present the range of cosmogenic neutrino fluxes in the dip-model which are compatible with a recent determination of the extragalactic very high energy (VHE) gamma-ray diffuse background derived from 2.5 years of Fermi/LAT data. We show that the largest fluxes predicted in the dip model would be detectable by IceCube in about 10 years of observation and are within the reach of a few years of observation with the ARA project. In the incomplete UHECR model in which protons are assumed to dominate only above 10$^{19}$ eV, the cosmogenic neutrino fluxes could be a factor of 2 or 3 larger. Any fraction of heavier nuclei in the UHECR at these energies would reduce the maximum cosmogenic neutrino fluxes. We also consider here special evolution models in which the UHECR sources are assumed to have the same evolution of either the star formation rate (SFR), or the gamma-ray burst (GRB) rate, or the active galactic nuclei (AGN) rate in the Universe and found that the last two are disfavored (and in the dip model rejected) by the new VHE gamma-ray background.

After a brief introduction to dark matter in general and to WIMPs as candidates, we review recent results of direct dark matter searches. We concentrate on older and more recent hints pointing to light WIMP's with mass below 10 GeV.

The channeling of the ion recoiling after a collision with a WIMP in direct dark matter crystalline detectors produces a larger scintillation or ionization signal than otherwise expected. Channeling is a directional effect which depends on the velocity distribution of WIMPs in the dark halo of our Galaxy and could lead to a daily modulation of the signal. Here we compute upper bounds to the expected amplitude of daily modulation due to channeling using channeling fractions that we obtained with analytic models in prior work. After developing the general formalism, we examine the possibility of finding a daily modulation due to channeling in the data already collected by the DAMA/NaI and DAMA/LIBRA experiments. We find that even the largest daily modulation amplitudes (of the order of 10% in some instances) would not be observable for WIMPs in the standard halo in the 13 years of data taken by the DAMA collaboration. For these to be observable the DAMA total rate should be 1/40 of what it is or the total DAMA exposure should be 40 times larger. The daily modulation due to channeling will be difficult to measure in future experiments. We find it could be observed for light WIMPs in solid Ne, assuming no background.

The channeling of the ion recoiling after a collision with a WIMP changes the ionization signal in direct detection experiments, producing a larger scintillation or ionization signal than otherwise expected. We give estimates of the fraction of channeled recoiling ions in solid Xe, Ar and Ne crystals using analytic models produced since the 1960's and 70's to describe channeling and blocking effects.

The channeling of the recoiling nucleus in crystalline detectors after a WIMP collision would produce a larger scintillation or ionization signal in direct detection experiments than otherwise expected. I present estimates of the importance of this effect in NaI, Si and Ge crystals, using analytic models developed from the 1960's onwards to describe channeling and blocking effects.

The channeling of the ion recoiling after a collision with a WIMP changes the ionization signal in direct detection experiments, producing a larger signal than otherwise expected. We give estimates of the fraction of channeled recoiling ions in NaI (Tl) crystals using analytic models produced since the 1960's and 70's to describe channeling and blocking effects. We find that the channeling fraction of recoiling lattice nuclei is smaller than that of ions that are injected into the crystal and that it is strongly temperature dependent.

The channeling of the ion recoiling after a collision with a WIMP changes the ionization signal in direct detection experiments, producing a larger signal scintillation or ionization than otherwise expected. We give estimates of the fraction of channeled recoiling ions in CsI crystals using analytic models produced since the 1960's and 70's to describe channeling and blocking effects.

The channeling of the ion recoiling after a collision with a WIMP changes the ionization signal in direct detection experiments, producing a larger signal than otherwise expected. We give estimates of the fraction of channeled recoiling ions in Si and Ge crystals using analytic models produced since the 1960's and 70's to describe channeling and blocking effects. We used data obtained to avoid channeling in the implantation of dopants in Si crystals to test our models.

Data collected by the Pierre Auger Observatory through 31 August 2007 showed evidence for anisotropy in the arrival directions of cosmic rays above the Greisen-Zatsepin-Kuz'min energy threshold, \nobreak{$6\times 10^{19}$eV}. The anisotropy was measured by the fraction of arrival directions that are less than $3.1^\circ$ from the position of an active galactic nucleus within 75 Mpc (using the V\'eron-Cetty and V\'eron $12^{\rm th}$ catalog). An updated measurement of this fraction is reported here using the arrival directions of cosmic rays recorded above the same energy threshold through 31 December 2009. The number of arrival directions has increased from 27 to 69, allowing a more precise measurement. The correlating fraction is $(38^{+7}_{-6})%$, compared with $21%$ expected for isotropic cosmic rays. This is down from the early estimate of $(69^{+11}_{-13})%$. The enlarged set of arrival directions is examined also in relation to other populations of nearby extragalactic objects: galaxies in the 2 Microns All Sky Survey and active galactic nuclei detected in hard X-rays by the Swift Burst Alert Telescope. A celestial region around the position of the radiogalaxy Cen A has the largest excess of arrival directions relative to isotropic expectations. The 2-point autocorrelation function is shown for the enlarged set of arrival directions and compared to the isotropic expectation.

We review the mechanism of production of dark matter particles in the early Universe, both in standard and non-standard pre-Big Bang Nucleosynthesis cosmologies. We concentrate mostly on the production of WIMPs.

The search for dark matter is a very wide and active field of research. Many potential hints of dark matter have appeared recently which led to a burst of theoretical activity and model building. I necessarily concentrate here only in some aspects of it. I review here some recent hints and some of the ways in which they could be explained.

We consider the compatibility of DAMA/LIBRA, CoGeNT, XENON10 and XENON100 results for spin-independent (SI) dark matter Weakly Interacting Massive Particles (WIMPs), particularly at low masses (~ 10 GeV), assuming a standard dark matter halo. The XENON bounds depend on the scintillation efficiency factor Leff for which there is considerable uncertainty. Thus we consider various extrapolations for Leff at low energy. With the Leff measurements we consider, XENON100 results are found to be insensitive to the low energy extrapolation. We find the strongest bounds are from XENON10, rather than XENON100, due to the lower energy threshold. For reasonable choices of Leff and for the case of SI elastic scattering, XENON10 is incompatible with the DAMA/LIBRA 3$\sigma$ region and severely constrains the 7-12 GeV WIMP mass region of interest published by the CoGeNT collaboration.

We calculate contained and upward muon flux and contained shower event rates from neutrino interactions, when neutrinos are produced from annihilation of the dark matter in the Galactic Center. We consider model-independent direct neutrino production and secondary neutrino production from the decay of taus, W bosons and bottom quarks produced in the annihilation of dark matter. We illustrate how muon flux from dark matter annihilation has a very different shape than the muon flux from atmospheric neutrinos. We also discuss the dependence of the muon fluxes on the dark matter density profile and on the dark matter mass and of the total muon rates on the detector threshold. We consider both the upward muon flux, when muons are created in the rock below the detector, and the contained flux when muons are created in the (ice) detector. We also calculate the event rates for showers from neutrino interactions in the detector and show that the signal dominates over the background for $150 {\rm GeV} <m_\chi < 1$ TeV for $E_{sh}^{th} = 100$ GeV.

In usual particle models, sterile neutrinos can account for the dark matter of the Universe only if they have masses in the keV range and are warm dark matter. Stringent cosmological and astrophysical bounds, in particular imposed by X-ray observations, apply to them. We point out that in a particular variation of the inert doublet model, sterile neutrinos can account for the dark matter in the Universe and may be either cold or warm dark matter candidates, even for masses much larger than the keV range. These Inert-Sterile neutrinos, produced non-thermally in the early Universe, would be stable and have very small couplings to Standard Model particles, rendering very difficult their detection in either direct or indirect dark matter searches. They could be, in principle, revealed in colliders by discovering other particles in the model.

The Pierre Auger Observatory is a hybrid detector for ultra-high energy cosmic rays. It combines a surface array to measure secondary particles at ground level together with a fluorescence detector to measure the development of air showers in the atmosphere above the array. The fluorescence detector comprises 24 large telescopes specialized for measuring the nitrogen fluorescence caused by charged particles of cosmic ray air showers. In this paper we describe the components of the fluorescence detector including its optical system, the design of the camera, the electronics, and the systems for relative and absolute calibration. We also discuss the operation and the monitoring of the detector. Finally, we evaluate the detector performance and precision of shower reconstructions.

Studies of the cosmic ray energy spectrum at the highest energies with the Pierre Auger Observatory.

Technical reports on operations and features of the Pierre Auger Observatory, including ongoing and planned enhancements and the status of the future northern hemisphere portion of the Observatory. Contributions to the 31st International Cosmic Ray Conference, Lodz, Poland, July 2009.

Studies of the correlations of ultra-high energy cosmic ray directions with extra-Galactic objects, of general anisotropy, of photons and neutrinos, and of other astrophysical effects, with the Pierre Auger Observatory. Contributions to the 31st ICRC, Lodz, Poland, July 2009.

Reports on the atmospheric monitoring, calibration, and other operating systems of the Pierre Auger Observatory. Contributions to the 31st International Cosmic Ray Conference, Lodz, Poland, July 2009.

Studies of the composition of the highest energy cosmic rays with the Pierre Auger Observatory, including examination of hadronic physics effects on the structure of extensive air showers. Submissions to the 31st ICRC, Lodz, Poland (July 2009).

From direct observations of the longitudinal development of ultra-high energy air showers performed with the Pierre Auger Observatory, upper limits of 3.8%, 2.4%, 3.5% and 11.7% (at 95% c.l.) are obtained on the fraction of cosmic-ray photons above 2, 3, 5 and 10 EeV (1 EeV = 10^18 eV) respectively. These are the first experimental limits on ultra-high energy photons at energies below 10 EeV. The results complement previous constraints on top-down models from array data and they reduce systematic uncertainties in the interpretation of shower data in terms of primary flux, nuclear composition and proton-air cross-section.

I review here some of the physics we are learning and expect to learn in the near future through the observation of cosmic rays. The study of cosmic rays involves a combination of data from accelerators, ground arrays, atmospheric fluorescence detectors and balloon and satellite experiments. I will discuss the data of the Pierre Auger Observatory, PAMELA, ATIC and FST among other experiments.

The DAMA/NaI and DAMA/LIBRA annual modulation data, which may be interpreted as a signal for the existence of weakly interacting dark matter (WIMPs) in our galactic halo, are examined in light of null results from other experiments. We use the energy spectrum of the combined DAMA modulation data given in 36 bins, and include the effect of channeling. Several statistical tools are implemented in our study: likelihood ratio with a global fit and with raster scans in the WIMP mass and goodness-of-fit (g.o.f.). These approaches allow us to differentiate between the preferred (global best fit) and allowed (g.o.f.) parameter regions. It is hard to find WIMP masses and couplings consistent with all existing data sets. For spin-independent (SI) interactions, the best fit DAMA regions are ruled out to the 3$\sigma$ C.L., even with channeling taken into account. However, for WIMP masses of ~8 GeV some parameters outside these regions still yield a moderately reasonable fit to the DAMA data and are compatible with all 90% C.L. upper limits from negative searches, when channeling is included. For spin-dependent (SD) interactions with proton-only couplings, a range of masses below 10 GeV is compatible with DAMA and other experiments, with and without channeling, when SuperK indirect detection constraints are included; without the SuperK constraints, masses as high as ~20 GeV are compatible. For SD neutron-only couplings we find no parameters compatible with all the experiments. Mixed SD couplings are examined: e.g. ~8 GeV mass WIMPs with a_n = +/- a_p are found to be consistent with all experiments. In short, there are surviving regions at low mass for both SI and SD interactions; if indirect detection limits are relaxed, some SD proton-only couplings at high masses also survive.

Weakly interacting massive particles (WIMPs) are one of very few probes of cosmology before Big Bang nucleosynthesis (BBN). We point out that in scenarios in which the Universe evolves in a non-standard manner during and after WIMP kinetic decoupling, the horizon mass scale at decoupling can be smaller and the dark matter WIMPs can be colder than in standard cosmology. This would lead to much smaller first objects in hierarchical structure formation. In low reheating temperature scenarios the effect may be large enough as to noticeably enhance indirect detection signals in GLAST and other detectors, by up to two orders of magnitude.

It is commonly assumed that the cosmological and astrophysical bounds on the mixings of sterile with active neutrinos are much more stringent than those obtained from laboratory measurements. We point out that in scenarios with a very low reheating temperature T_RH << 100 MeV at the end of (the last episode of) inflation or entropy creation, the abundance of sterile neutrinos becomes largely suppressed with respect to that obtained within the standard framework. Thus, in this case cosmological bounds become much less stringent than usually assumed, allowing sterile neutrinos to be ``visible'' in future experiments. Here, we concentrate on massive (mostly sterile) neutrinos heavier than 1 MeV.

We fit the recently published Pierre Auger ultra-high energy cosmic ray spectrum assuming that either nucleons or nuclei are emitted at the sources. We consider the simplified cases of pure proton, or pure oxygen, or pure iron injection. We perform an exhaustive scan in the source evolution factor, the spectral index, the maximum energy of the source spectrum Z E_{max}, and the minimum distance to the sources. We show that the Pierre Auger spectrum agrees with any of the source compositions we assumed. For iron, in particular, there are two distinct solutions with high and low E_{max} (e.g. 6.4 10^{20} eV and 2 10^{19} eV) respectively which could be distinguished by either a large fraction or the near absence of proton primaries at the highest energies. We raise the possibility that an iron dominated injected flux may be in line with the latest composition measurement from the Pierre Auger Observatory where a hint of heavy element dominance is seen.

We calculate the flux of "GZK-photons", namely the flux of Ultra High Energy Cosmic Rays (UHECR) consisting of photons produced by extragalactic nucleons through the resonant photoproduction of pions, the so called GZK effect. We We calculate the flux of "GZK-photons", namely the flux of Ultra High Energy Cosmic Rays (UHECR) consisting of photons produced by extragalactic nucleons through the resonant photoproduction of pions, the so called GZK effect. We show that, for primary nucleons, the GZK photon fraction of the total UHECR flux is between $10^{-4}$ and $10^{-2}$ above $10^{19}$ eV and up to the order of 0.1 above $10^{20}$ eV. The GZK photon flux depends on the assumed UHECR spectrum, slope of the nucleon flux at the source, distribution of sources and intervening backgrounds. Detection of this photon flux would open the way for UHECR gamma-ray astronomy. Detection of a larger photon flux would imply the emission of photons at the source or new physics. We compare the photon fractions expected for GZK photons and the minimal predicted by Top-Down models. We find that the photon fraction above $10^{19}$ eV is a crucial test for Top-Down models.

In a recently proposed model the cosmic rays spectrum at energies above EeV can be fitted with a minimal number of unknown parameters assuming that the extragalactic cosmic rays are only protons with a power law source spectrum. Within this minimal model, after fitting the observed HiRes spectrum with four parameters (proton injection spectrum power law index, maximum energy, minimum distance to sources and evolution parameter) we compute the flux of ultra-high energy photons due to photon-pion production and e+e- pair production by protons for several radio background models and a range of average extragalactic magnetic fields.

We calculate the flux of "GZK-photons", namely the flux of photons produced by extragalactic nucleons through the resonant photoproduction of pions, the so called GZK effect. This flux depends on the UHECR spectrum on Earth, of the spectrum of nucleons emitted at the sources, which we characterize by its slope and maximum energy, on the distribution of sources and on the intervening cosmological backgrounds, in particular the magnetic field and radio backgrounds. For the first time we calculate the GZK photons produced by nuclei. We calculate the possible range of the GZK photon fraction of the total UHECR flux for the AGASA and the HiRes spectra. We find that for nucleons produced at the sources it could be as large as a few % and as low as 10^{-4} above 10 EeV. For nuclei produced at the sources the maximum photon fraction is a factor of 2 to 3 times smaller above 10 EeV but the minimum could be much smaller than for nucleons. We also comment on cosmogenic neutrino fluxes.

In a recently proposed model the cosmic rays spectrum at energies above 10^{18} eV can be fitted with a minimal number of unknown parameters assuming that the extragalactic cosmic rays are only protons with a power law source spectrum ~E^{-alpha} and alpha~2.6. Within this minimal model, after fitting the observed HiRes spectrum with four parameters (proton injection spectrum power law index and maximum energy, minimum distance to sources and evolution parameter) we compute the flux of ultra-high energy photons due to photon-pion production, the GZK photons, for several radio background models and average extragalactic magnetic fields with amplitude between 10^{-11} G and 10^{-9} G. We find the photon fraction to be between 10^{-4} and 10^{-3} in cosmic rays at energies above 10^{19} eV. These small fluxes could only be detected in future experiments like Auger North plus South and EUSO.

We show that Majoron emission from a hot nascent neutron star can be anisotropic in the presence of a strong magnetic field. If Majorons carry a non-negligible fraction of the supernova energy, the resulting recoil velocity of a neutron star can explain the observed velocities of pulsars.

Neutrinos from far-away sources annihilating at the Z resonance on relic neutrinos may give origin to the extreme-energy cosmic rays (EECR). If ``Z-bursts'' are responsible for the EECR events, then we show that the non-observation of cosmic ray events at energies above 2 x 10^20 eV by the AGASA Collaboration implies a lower bound around 0.3 eV on the relic neutrino mass. Since this mass exceeds the mass-squared differences inferred from oscillation physics, the bound in fact applies to all three neutrino masses. Together with the upper bound provided by comparisons of the CMB anisotropy with large-scale structure, this bound leaves only a small interval for neutrino masses around 0.3 eV, if Z-bursts are to explain the existing EECR events.

We present here a scenario, based on a low reheating temperature T_R << 100 MeV at the end of (the last episode of) inflation, in which the coupling of sterile neutrinos to active neutrinos can be as large as experimental bounds permit (thus making this neutrino ``visible'' in future experiments). In previous models this coupling was forced to be very small to prevent a cosmological overabundance of sterile neutrinos. Here the abundance depends on how low the reheating temperature is. For example, the sterile neutrino required by the LSND result does not have any cosmological problem within our scenario.

We show that direct dark matter detection experiments can distinguish between pointlike and non-pointlike dark-matter candidates. The shape of the nuclear recoil energy spectrum from pointlike dark-matter particles, e.g., neutralinos, is determined by the velocity distribution of dark matter in the galactic halo and by nuclear form factors. In contrast, typical cross sections of non-pointlike dark matter, for example, Q-balls, have a new form factor, which decreases rapidly with the recoil energy. Therefore, a signal from non-pointlike dark matter is expected to peak near the experimental threshold and to fall off rapidly at higher energies. Although the width of the signal is practically independent of the dark matter velocity dispersion, its height is expected to exhibit an annual modulation due to the changes in the dark matter flux.

Neutrinos from far away sources annihilating at the Z resonance on relic neutrinos may give origin to the ultrahigh-energy cosmic rays. Here we present predictions of this mechanism with relic neutrinos lighter than 1 eV, which do not gravitationally cluster. We show that not only the super GZK events, but the "ankle" and all events above it can be accounted for. Most primaries above the ankle are predicted to be nucleons up to $10^{20.0}$ eV and photons at higher energies. We also find an accumulation at the GZK cutoff energy, a hint of which can be seen in the data.

Relic neutrinos with mass in the range indicated by Super-Kamiokande results if neutrino masses are hierarchial (about 0.07 eV) are many times deemed too light to be cosmologically relevant. Here we remark that these neutrinos may significantly contribute to the dark matter of the Universe (with a large lepton asymmetry $L$) and that their existence might be revealed by the spectrum of ultra high energy cosmic rays (maybe even in the absence of a large $L$).

Decays of superheavy relic particles may produce extremely energetic neutrinos. Their annihilations on the relic neutrinos can be the origin of the cosmic rays with energies beyond the Greisen-Zatsepin-Kuzmin cutoff. The red shift acts as a cosmological filter selecting the sources at some particular value z_e, for which the present neutrino energy is close to the Z pole of the annihilation cross section. We predict no directional correlation of the ultrahigh-energy cosmic rays with the galactic halo. At the same time, there can be some directional correlations in the data, reflecting the distribution of matter at red shift z=z_e. Both of these features are manifest in the existing data. Our scenario is consistent with the neutrino mass reported by Super-Kamiokande and requires no lepton asymmetry or clustering of the background neutrinos.

Relic neutrinos with mass 0.07 (+0.02/-0.04) eV, in the range consistent with Super-Kamiokande data, can explain the cosmic rays with energies in excess of the Greisen-Zatsepin-Kuzmin cutoff. The spectrum of ultra-high energy cosmic rays produced in this fashion has some distinctive features that may help identify their origin. Our mechanism does not require but is consistent with a neutrino density high enough to be a new kind of hot dark matter.

A supersymmetric inflationary stage dominated by an $F$-term has the problem that the flatness of the potential is spoiled by supergravity corrections, that is the slow-roll parameter $\eta$ gets contributions of order unity. We show that in $F$-term inflationary models based on strings there is natural way of obtaining small values of $\eta$. This happens in models of hybrid inflation based on orbifold constructions, in which a modulus $T$ field is responsible for the large value of the potential during inflation, and a second field $\phi$ with appropriate modular weight is responsible for the roll-over. We illustrate the mechanism with a model in which the inflaton potential is provided by gaugino condensation, leading to succesful inflation.

We consider the effects of a massive, unstable neutrino on the evolution of large--scale structure and anisotropies in the cosmic microwave background. Comparison with large--scale structure data allows us to rule out a wide range of masses and lifetimes for such neutrinos. We also define a range of masses and lifetimes which delay matter--radiation equality and improve the agreement with the data of Cold Dark Matter models with critical density.