We present results from the FUNK experiment in the search for hidden-photon dark matter. Near the surface of a mirror, hidden photons may be converted into ordinary photons. These photons are emitted perpendicular to the surface and have an energy equal to the mass of the dark matter hidden photon. Our experimental setup consists of a large, spherical mirror with an area of more than 14 m$^2$, which concentrates the emitted photons into its central point. Using a detector sensitive to visible and near-UV photons, we can exclude a kinetic-mixing coupling of stronger than $\chi \approx 10^{-12}$ in the mass range of 2.5 to 7 eV, assuming hidden photons comprise all of the dark matter. The experimental setup and analysis used to obtain this limit are discussed in detail.

Dark matter comprised of axion-like particles (ALPs) generated by the realignment mechanism in the post-inflationary scenario leads to primordial isocurvature fluctuations. The power spectrum of these fluctuations is flat for small wave numbers, extending to scales accessible with cosmological surveys. We use the latest measurements of Cosmic Microwave Background (CMB) primary anisotropies together with CMB lensing, Baryonic Acoustic Oscillations (BAO) and Sunyaev Zel'dovich (SZ) cluster counts to measure the amplitude and tilt of the isocurvature component. We find preference for a white-noise isocurvature component in the CMB primary anisotropies; this conclusion is, however, weakened by current large-scale structure (LSS) data. Interpreting the result as a conservative upper limit on the isocurvature component, the combined bound on the ALP mass from all probes is $m_{a} \gtrsim 10^{-19}$ eV, with some dependence on how $m_{a}$ evolves with temperature. The expected sensitivity of cosmic shear and galaxy clustering from future LSS experiments and CMB lensing suggests improved bounds of $m_{a} \gtrsim 10^{-18}$-$10^{-13}$ eV, depending on scale cuts used to avoid non-linearities and the ALP mass-temperature dependence.

Axion--photon interactions can lead to an enhancement of the electromagnetic field by parametric resonance in the presence of a cold axion background, for modes with a frequency close to half the axion mass. In this paper, we study the role of the axion momentum dispersion as well as the effects of a background gravitational potential, which can detune the resonance due to gravitational redshift. We show, by analytical as well as numerical calculations, that the resonance leads to an exponential growth of the photon field only if (a) the axion momentum spread is smaller than the inverse resonance length, and (b) the gravitational detuning distance is longer than the resonance length. For realistic parameter values, both effects strongly suppress the resonance and prevent the exponential growth of the photon field. In particular, the redshift due to the gravitational potential of our galaxy prevents the resonance from developing for photons in the observable frequency range, even assuming that all the dark matter consists of a perfectly cold axion condensate. For axion clumps with masses below $\sim 10^{-13}\, M_\odot$, the momentum spread condition is more restrictive, whereas, for more massive clumps, the redshift condition dominates.

Next-generation dark matter direct detection experiments will explore several orders of magnitude in the dark matter--nucleus scattering cross section below current upper limits. In case a signal is discovered the immediate task will be to determine the dark matter mass and to study the underlying interactions. We develop a framework to determine the dark matter mass from signals in two experiments with different targets, independent of astrophysics. Our method relies on a distribution-free, nonparametric two-sample hypothesis test in velocity space, which neither requires binning of the data, nor any fitting of parametrisations of the velocity distribution. We apply our method to realistic configurations of xenon and argon detectors such as XENONnT/DARWIN and DarkSide, and estimate the precision with which the DM mass can be determined. Once the dark matter mass is identified, the ratio of coupling strengths to neutrons and protons can be constrained by using the same data. The test can be applied for event samples of order 20 events, but promising sensitivities require $\gtrsim 100$ events.

We assume that dark matter is comprised of axion-like particles (ALPs) generated by the realignment mechanism in the post-inflationary scenario. This leads to isocurvature fluctuations with an amplitude of order one for scales comparable to the horizon at the time when the ALP field starts oscillating. The power spectrum of these fluctuations is flat for small wave numbers, extending to scales relevant for cosmological observables. Denoting the relative isocurvature amplitude at $k_*$ = 0.05 Mpc$^{-1}$ by $f_{\rm iso}$, Planck observations of the cosmic microwave background (CMB) yield $f_{\rm iso}$ < 0.31 at the 2$\sigma$-level. This excludes the hypothesis of post-inflationary ALP dark matter with masses $m_{a}$ < 10$^{-20}-$10$^{-16}$ eV, where the range is due to details of the ALP mass-temperature dependence. Future CMB stage IV and 21-cm intensity mapping experiments may improve these limits by 1$-$2 orders of magnitude in $m_{a}$.

When Peccei-Quinn (PQ) symmetry breaking happens after inflation, the axion field takes random values in causally disconnected regions. This leads to fluctuations of order one in the axion energy density around the QCD epoch. These over-densities eventually decouple from the Hubble expansion and form so-called miniclusters. We present a semi-analytical method to calculate the average axion energy density, as well as the power spectrum, from the re-alignment mechanism in this scenario. Furthermore, we develop a modified Press & Schechter approach, suitable to describe the collapse of non-linear density fluctuations during radiation domination, which is relevant for the formation of axion miniclusters. It allows us to calculate the double differential distribution of gravitationally collapsed miniclusters as a function of their mass and size. For instance, assuming a PQ scale of $10^{11}$ GeV, minicluster masses range from about $5 \times 10^{-16}$ to $3 \times 10^{-13}$ solar masses and have sizes from about $4\times 10^4$ to $7\times 10^5$ km at the time they start to collapse.

In the preprint arxiv:1703.03425 "strong evidence" for the normal neutrino mass ordering is claimed. The authors obtain Bayesian odds of 42:1 in favour of the normal ordering. Their conclusion is based on adopting a flat logarithmic prior for the three neutrino masses. Such an assumption favours a hierarchical spectrum for the masses, which is much easier to accommodate for the normal mass ordering, and hence their prior assumption makes the inverted ordering much less likely a priori. We argue that the claimed "evidence" for normal ordering is almost entirely driven by the adopted prior and not due to the data itself.

We present a comprehensive review of keV-scale sterile neutrino Dark Matter, collecting views and insights from all disciplines involved - cosmology, astrophysics, nuclear, and particle physics - in each case viewed from both theoretical and experimental/observational perspectives. After reviewing the role of active neutrinos in particle physics, astrophysics, and cosmology, we focus on sterile neutrinos in the context of the Dark Matter puzzle. Here, we first review the physics motivation for sterile neutrino Dark Matter, based on challenges and tensions in purely cold Dark Matter scenarios. We then round out the discussion by critically summarizing all known constraints on sterile neutrino Dark Matter arising from astrophysical observations, laboratory experiments, and theoretical considerations. In this context, we provide a balanced discourse on the possibly positive signal from X-ray observations. Another focus of the paper concerns the construction of particle physics models, aiming to explain how sterile neutrinos of keV-scale masses could arise in concrete settings beyond the Standard Model of elementary particle physics. The paper ends with an extensive review of current and future astrophysical and laboratory searches, highlighting new ideas and their experimental challenges, as well as future perspectives for the discovery of sterile neutrinos.

We propose a simple method to quantify a possible exclusion of the inverted neutrino mass ordering from cosmological bounds on the sum of the neutrino masses. The method is based on Bayesian inference and allows for a calculation of the posterior odds of normal versus inverted ordering. We apply the method for a specific set of current data from Planck CMB data and large-scale structure surveys, providing an upper bound on the sum of neutrino masses of 0.14 eV at 95% CL. With this analysis we obtain posterior odds for normal versus inverted ordering of about 2:1. If cosmological data is combined with data from oscillation experiments the odds reduce to about 3:2. For an exclusion of the inverted ordering from cosmology at more than 95% CL, an accuracy of better than 0.02 eV is needed for the sum. We demonstrate that such a value could be reached with planned observations of large scale structure by analysing artificial mock data for a EUCLID-like survey.

A reliable comparison of different dark matter (DM) searches requires models that satisfy certain consistency requirements like gauge invariance and perturbative unitarity. As a well-motivated example, we study two-mediator DM (2MDM). The model is based on a spontaneously broken $U(1)'$ gauge symmetry and contains a Majorana DM particle as well as two $s$-channel mediators, one vector (the $Z'$) and one scalar (the dark Higgs). We perform a global scan over the parameters of the model assuming that the DM relic density is obtained by thermal freeze-out in the early Universe and imposing a large set of constraints: direct and indirect DM searches, monojet, dijet and dilepton searches at colliders, Higgs observables, electroweak precision tests and perturbative unitarity. We conclude that thermal DM is only allowed either close to an $s$-channel resonance or if at least one mediator is lighter than the DM particle. In these cases a thermal DM abundance can be obtained although DM couplings to the Standard Model are tiny. Interestingly, we find that vector-mediated DM-nucleon scattering leads to relevant constraints despite the velocity-suppressed cross section, and that indirect detection can be important if DM annihilations into both mediators are kinematically allowed.

We consider how QCD axions produced by the misalignment mechanism could form galactic dark matter halos. We recall that stationary, gravitationally stable axion field configurations have the size of an asteroid with masses of order $10^{-13} $ solar masses (because gradient pressure is insufficient to support a larger object). We call such field configurations "drops". We explore whether rotating drops could be larger, and find that their mass could increase by a factor ~ 10. Remarkably this mass is comparable to the mass of miniclusters generated from misalignment axions in the scenario where the axion is born after inflation. We speculate that misalignment axions today are in the form of drops, contributing to dark matter like a distribution of asteroids (and not as a coherent oscillating background field). We consider some observational signatures of the drops, which seem consistent with a galactic halo made of axion dark matter.

We show that simplified models used to describe the interactions of dark matter with Standard Model particles do not in general respect gauge invariance and that perturbative unitarity may be violated in large regions of the parameter space. The modifications necessary to cure these inconsistencies may imply a much richer phenomenology and lead to stringent constraints on the model. We illustrate these observations by considering the simplified model of a fermionic dark matter particle and a vector mediator. Imposing gauge invariance then leads to strong constraints from dilepton resonance searches and electroweak precision tests. Furthermore, the new states required to restore perturbative unitarity can mix with Standard Model states and mediate interactions between the dark and the visible sector, leading to new experimental signatures such as invisible Higgs decays. The resulting constraints are typically stronger than the 'classic' constraints on DM simplified models such as monojet searches and make it difficult to avoid thermal overproduction of dark matter.

We summarize the determination of some neutrino properties from the global analysis of solar, atmospheric, reactor, and accelerator neutrino data in the framework of three-neutrino mixing as well as in some extended scenarios such as the mixing with eV-scale sterile neutrinos invoked for the interpretation of the short baseline anomalies, and the presence of non-standard neutrino interactions.

We show that a positive signal in a dark matter (DM) direct detection experiment can be used to place a lower bound on the DM capture rate in the Sun, independent of the DM halo. For a given particle physics model and DM mass we obtain a lower bound on the capture rate independent of the local DM density, velocity distribution, galactic escape velocity, as well as the scattering cross section. We illustrate this lower bound on the capture rate by assuming that upcoming direct detection experiments will soon obtain a significant signal. When comparing the lower bound on the capture rate with limits on the high-energy neutrino flux from the Sun from neutrino telescopes, we can place upper limits on the branching fraction of DM annihilation channels leading to neutrinos. With current data from IceCube and Super-Kamiokande non-trivial limits can be obtained for spin-dependent interactions and direct annihilations into neutrinos. In some cases also annihilations into $\tau\tau$ or $b\bar b$ start getting constrained. For spin-independent interactions current constraints are weak, but they may become interesting for data from future neutrino telescopes.

In Dark Matter direct detection we are facing the situation of some experiments reporting positive signals which are in conflict with limits from other experiments. Such conclusions are subject to large uncertainties introduced by the poorly known local Dark Matter distribution. We present a method to calculate an upper bound on the joint probability of obtaining the outcome of two potentially conflicting experiments under the assumption that the Dark Matter hypothesis is correct, but completely independent of assumptions about the Dark Matter distribution. In this way we can quantify the compatibility of two experiments in an astrophysics independent way. We illustrate our method by testing the compatibility of the hints reported by DAMA and CDMS-Si with the limits from the LUX and SuperCDMS experiments. The method does not require Monte Carlo simulations but is mostly based on using Poisson statistics. In order to deal with signals of few events we introduce the so-called "signal length" to take into account energy information. The signal length method provides a simple way to calculate the probability to obtain a given experimental outcome under a specified Dark Matter and background hypothesis.

We consider the effect of the Sun's gravitational potential on the local phase space distribution of dark matter particles, focusing on its implication for the annual modulation signal in direct detection experiments. We perform a fit to the modulation signal observed in DAMA/LIBRA and show that the allowed region shrinks if Solar gravitational focusing (GF) is included compared to the one without GF. Furthermore, we consider a possible signal in a generic future direct detection experiment, irrespective of the DAMA/LIBRA signal. Even for scattering cross sections close to the current bound and a large exposure of a xenon target with 270 ton yr it will be hard to establish the presence of GF from data. In the region of dark matter masses below 40 GeV an annual modulation signal can be established for our assumed experimental setup, however GF is negligible for low masses. In the high mass region, where GF is more important, the significance of annual modulation itself is very low. We obtain similar results for lighter targets such as Ge and Ar. We comment also on inelastic scattering, noting that GF becomes somewhat more important for exothermic scattering compared to the elastic case.

Neutrino oscillations involving eV-scale neutrino mass states are investigated in the context of global neutrino oscillation data including short and long-baseline accelerator, reactor, and radioactive source experiments, as well as atmospheric and solar neutrinos. We consider sterile neutrino mass schemes involving one or two mass-squared differences at the eV^2 scale denoted by 3+1, 3+2, and 1+3+1. We discuss the hints for eV-scale neutrinos from nu_e disappearance (reactor and Gallium anomalies) and nu_mu->nu_e appearance (LSND and MiniBooNE) searches, and we present constraints on sterile neutrino mixing from nu_mu and neutral-current disappearance data. An explanation of all hints in terms of oscillations suffers from severe tension between appearance and disappearance data. The best compatibility is obtained in the 1+3+1 scheme with a p-value of 0.2% and exceedingly worse compatibilities in the 3+1 and 3+2 schemes.

Dark matter N-body simulations suggest that the velocity distribution of dark matter is anisotropic. In this work we employ a mass model for the Milky Way whose parameters are determined from a fit to kinematical data. Then we adopt an ansatz for the dark matter phase space distribution which allows to construct self-consistent halo models which feature a degree of anisotropy as a function of the radius such as suggested by the simulations. The resulting velocity distributions are then used for an analysis of current data from dark matter direct detection experiments. We find that velocity distributions which are radially biased at large galactocentric distances (up to the virial radius) lead to an increased high velocity tail of the local dark matter distribution. This affects the interpretation of data from direct detection experiments, especially for dark matter masses around 10 GeV, since in this region the high velocity tail is sampled. We find that the allowed regions in the dark matter mass-cross section plane as indicated by possible hints for a dark matter signal reported by several experiments as well as conflicting exclusion limits from other experiments shift in a similar way when the halo model is varied. Hence, it is not possible to improve the consistency of the data by referring to anisotropic halo models of the type considered in this work.

We present halo-independent methods to analyze the results of dark matter direct detection experiments assuming inelastic scattering. We focus on the annual modulation signal reported by DAMA/LIBRA and present three different halo-independent tests. First, we compare it to the upper limit on the unmodulated rate from XENON100 using (a) the trivial requirement that the amplitude of the annual modulation has to be smaller than the bound on the unmodulated rate, and (b) a bound on the annual modulation amplitude based on an expansion in the Earth's velocity. The third test uses the special predictions of the signal shape for inelastic scattering and allows for an internal consistency check of the data without referring to any astrophysics. We conclude that a strong conflict between DAMA/LIBRA and XENON100 in the framework of spin-independent inelastic scattering can be established independently of the local properties of the dark matter halo.

In the context of supersymmetric models in which small Dirac neutrino masses are generated by supersymmetry breaking, a mainly right-handed (RH) mixed sneutrino can be an excellent cold dark matter (DM) candidate. We perform a global analysis of the Minimal Supersymmetric Standard Model (MSSM)+RH neutrino parameter space by means of Markov Chain Monte Carlo scans. We include all relevant constraints from collider and dark matter searches, paying particular attention to nuclear and astrophysical uncertainties. Two distinct cases can satisfy all constraints: heavy sneutrino DM with mass of order 100 GeV, as well as light sneutrino DM with mass of about 3-6 GeV. We discuss the implications for direct and indirect dark matter searches, as well as for SUSY and Higgs searches at the LHC for both, the light and the heavy sneutrino dark matter case. The light sneutrino case is excluded by the 125-126 GeV Higgs signal.

We show how constraints on the time integrated event rate from a given dark matter (DM) direct detection experiment can be used to set a stringent constraint on the amplitude of the annual modulation signal in another experiment. The method requires only very mild assumptions about the properties of the local DM distribution: that it is temporally stable on the scale of months and spatially homogeneous on the ecliptic. We apply the method to the annual modulation signal in DAMA/LIBRA, which we compare to the bounds derived from the constraints on the time-averaged rates from XENON10, XENON100, CDMS and SIMPLE. Assuming a DM mass of 10 GeV, we show that a DM interpretation of the DAMA/LIBRA signal is excluded at 6.3sigma (4.6sigma) for isospin conserving (violating) spin-independent interactions, and at 4.9sigma for spin-dependent interactions on protons.

This white paper addresses the hypothesis of light sterile neutrinos based on recent anomalies observed in neutrino experiments and the latest astrophysical data.

We consider an electro-weak scale model for Dark Matter (DM) and radiative neutrino mass generation. Despite the leptophilic nature of DM with no direct couplings to quarks and gluons, scattering with nuclei is induced at the 1-loop level through photon exchange. Effectively, there are charge-charge, dipole-charge and dipole-dipole interactions. We investigate the parameter space consistent with constraints from neutrino masses and mixing, charged lepton-flavour violation, perturbativity, and the thermal production of the correct DM abundance, and calculate the expected event rate in DM direct detection experiments. We show that current data from XENON100 start to constrain certain regions of the allowed parameter space, whereas future data from XENON1T has the potential to significantly probe the model.

We derive constraints on the annual modulation signal in Dark Matter (DM) direct detection experiments in terms of the unmodulated event rate. A general bound independent of the details of DM distribution follows from the assumption that the motion of the earth around the sun is the only source of time variation. The bound is valid for a very general class of particle physics models and also holds in the presence of an unknown unmodulated background. More stringent bounds are obtained, if modest assumptions on symmetry properties of the DM halo are adopted. We illustrate the bounds by applying them to the annual modulation signals reported by the DAMA and CoGeNT experiments in the framework of spin-independent elastic scattering. While the DAMA signal satisfies our bounds, severe restrictions on the DM mass can be set for CoGeNT.

We present the results of a global neutrino oscillation data analysis within the three-flavour framework. We include latest results from the MINOS long-baseline experiment (including electron neutrino appearance as well as anti-neutrino data), updating all relevant solar (SK II+III), atmospheric (SK I+II+III) and reactor (KamLAND) data. Furthermore, we include a recent re-calculation of the anti-neutrino fluxes emitted from nuclear reactors. These results have important consequences for the analysis of reactor experiments and in particular for the status of the mixing angle $\theta_{13}$. In our recommended default analysis we find from the global fit that the hint for non-zero $\theta_{13}$ remains weak, at 1.8$\sigma$ for both neutrino mass hierarchy schemes. However, we discuss in detail the dependence of these results on assumptions concerning the reactor neutrino analysis.

New predictions for the anti-neutrino flux emitted by nuclear reactors suggest that reactor experiments may have measured a deficit in the anti-neutrino flux, which can be interpreted in terms of oscillations between the known active neutrinos and new sterile states. Motivated by this observation, we perform a re-analysis of global short-baseline neutrino oscillation data in a framework with one or two sterile neutrinos. While one sterile neutrino is still not sufficient to reconcile the signals suggested by reactor experiments and by the LSND and MiniBooNE experiments with null results from other searches, we find that, with the new reactor flux prediction, the global fit improves considerably when the existence of two sterile neutrinos is assumed.

Possible hints for WIMP dark matter with mass around 10 GeV coming from the DAMA, CoGeNT, and maybe also CRESST experiments are presented, and confronted with constraints from CDMS and XENON data. Focusing on spin-independent (SI) WIMP--nucleus interactions, I elaborate on the difficulties to make the hints consistent with each other and to evade the constraints, mentioning energy scale uncertainties, quenching and light-yield factors, as well as uncertainties on halo properties. In the present situation it seems hard to reconcile all data within the SI framework, which suggests that if the experimental anomalies were indeed due to dark matter a more exotic mechanism (to be identified) had to be at work.

We perform a global fit to data from Dark Matter (DM) direct detection experiments, including the recent CDMS-II results. We discuss possible interpretations of the DAMA annual modulation signal in terms of spin-independent and spin-dependent DM-nucleus interactions, both for elastic and inelastic scattering. We find that in the spin-dependent inelastic scattering off protons a good fit to all data is obtained. We present a simple toy model realizing such a scenario. In all the remaining cases the DAMA allowed regions are disfavored by other experiments or suffer from severe fine tuning of DM parameters with respect to the galactic escape velocity. Finally, we also entertain the possibility that the two events observed in CDMS-II are an actual signal of elastic DM scattering, and we compare the resulting CDMS-II allowed regions to the exclusion limits from other experiments. In this arXiv version of the manuscript we also provide in appendix A the updated fits including recent CoGeNT results.

We review the present status of three-flavour neutrino oscillations, taking into account the latest available neutrino oscillation data presented at the Neutrino 2008 Conference. This includes the data released this summer by the MINOS collaboration, the data of the neutral current counter phase of the SNO solar neutrino experiment, as well as the latest KamLAND and Borexino data. We give the updated determinations of the leading 'solar' and 'atmospheric' oscillation parameters. We find from global data that the mixing angle $\theta_{13}$ is consistent with zero within $0.9\sigma$ and we derive an upper bound of $\sin^2\theta_{13} < 0.035 (0.056)$ at 90% CL (3$\sigma$).

We consider the hypothesis that Dark Matter (DM) has tree-level interactions only with leptons. Such a framework, where DM recoils against electrons bound in atoms, has been proposed as an explanation for the annually modulated scintillation signal in DAMA/LIBRA versus the absence of a signal for nuclear recoils in experiments like CDMS or XENON10. However, even in such a leptophilic DM scenario there are loop induced DM-hadron interactions, where the photons emitted from virtual leptons couple to the charge of a nucleus. Using a general effective field theory approach we show that, if such an interaction is induced at one or two-loop level, then DM-nucleus scattering dominates over DM-electron scattering. This is because the latter is suppressed by the bound state wave function. One obtains a situation similar to standard DM-nucleus scattering analyses with considerable tension between the results of DAMA and CDMS/XENON10. This conclusion does not apply in the case of pseudoscalar or axial vector coupling between DM and leptons, where the loop diagrams vanish. In this case the explanation of the DAMA signal in terms of DM-electron scattering is strongly disfavored by the spectral shape of the signal. Furthermore, if DM can annihilate into neutrinos or tau leptons, the required cross sections are excluded by many orders of magnitude using the Super-Kamiokande bound on neutrinos from DM annihilations in the Sun.

We propose a scenario where Dark Matter (DM) annihilates into an intermediate state which travels a distance $\lambda \equiv v/\Gamma$ on the order of galactic scales and then decays to Standard Model (SM) particles. The long lifetime disperses the production zone of the SM particles away from the galactic center and hence, relaxes constraints from gamma ray observations on canonical annihilation scenarios. We utilize this set up to explain the electron and positron excesses observed recently by PAMELA, ATIC, and FERMI. While an explanation in terms of usual DM annihilations seems to conflict with gamma ray observations, we show that within the proposed scenario, the PAMELA/ATIC/FERMI results are consistent with the gamma ray data. The distinction from decay scenarios is discsussed and we comment on the prospects for DM production at LHC. The typical decay length $\lambda \gtrsim 10$ kpc of the intermediate state can have its origin from a dimension six operator suppressed by a scale $\Lambda \sim 10^{13}$ GeV, which is roughly the seesaw scale for neutrino masses.

We discuss the interpretation of the annual modulation signal seen in the DAMA experiment in terms of spin-independent elastic WIMP scattering. Taking into account channeling in the crystal as well as the spectral signature of the modulation signal we find that the low-mass WIMP region consistent with DAMA data is confined to WIMP masses close to $m_\chi \simeq 12$ GeV, in disagreement with the constraints from CDMS and XENON. We conclude that even if channeling is taken into account this interpretation of the DAMA modulation signal is disfavoured. There are no overlap regions in the parameter space at 90% CL and a consistency test gives the probability of $1.2\times 10^{-5}$. We study the robustness of this result with respect to variations of the WIMP velocity distribution in our galaxy, by changing various parameters of the distribution function, and by using the results of a realistic N-body dark matter simulation. We find that only by making rather extreme assumptions regarding halo properties can we obtain agreement between DAMA and CDMS/XENON.

We review the present status of global analyses of neutrino oscillations, taking into account the most recent neutrino data including the latest KamLAND and K2K updates presented at Neutrino2004, as well as state-of-the-art solar and atmospheric neutrino flux calculations. We give the two-neutrino solar + KamLAND results, as well as two-neutrino atmospheric + K2K oscillation regions, discussing in each case the robustness of the oscillation interpretation against departures from the Standard Solar Model and the possible existence of non-standard neutrino physics. Furthermore, we give the best fit values and allowed ranges of the three-flavour oscillation parameters from the current worlds' global neutrino data sample and discuss in detail the status of the small parameters $\alpha \equiv \Dms/\Dma$ as well as $\sin^2\theta_{13}$, which characterize the strength of CP violating effects in neutrino oscillations. We also update the degree of rejection of four-neutrino interpretations of the LSND anomaly in view of the most recent developments.

We propose an explanation of the LSND evidence for electron antineutrino appearance based on neutrino decay. We introduce a heavy neutrino, which is produced in pion and muon decays because of a small mixing with muon neutrinos, and then decays into a scalar particle and a light neutrino, predominantly of the electron type. We require values of $g m_4\sim$ few eV, $g$ being the neutrino--scalar coupling and $m_4$ the heavy neutrino mass, e.g. $m_4$ in the range from 1 keV to 1 MeV and $g \sim 10^{-6} - 10^{-3}$. Performing a fit to the LSND data as well as all relevant null-result experiments, we show that all data can be explained within this decay scenario. In the minimal version of the decay model, we predict a signal in the upcoming MiniBooNE experiment corresponding to a transition probability of the same order as seen in LSND. In addition, we show that extending our model to two nearly degenerate heavy neutrinos it is possible to introduce CP violation in the decay, which can lead to a suppression of the signal in MiniBooNE running in the neutrino mode. We briefly discuss signals in future neutrino oscillation experiments, we show that our scenario is compatible with bounds from laboratory experiments, and we comment on implications in astrophysics and cosmology.

We discuss a goodness-of-fit method which tests the compatibility between statistically independent data sets. The method gives sensible results even in cases where the chi^2-minima of the individual data sets are very low or when several parameters are fitted to a large number of data points. In particular, it avoids the problem that a possible disagreement between data sets becomes diluted by data points which are insensitive to the crucial parameters. A formal derivation of the probability distribution function for the proposed test statistic is given, based on standard theorems of statistics. The application of the method is illustrated on data from neutrino oscillation experiments, and its complementarity to the standard goodness-of-fit is discussed.

Prompted by recent solar and atmospheric data, we re-analyze the four-neutrino description of current global neutrino oscillation data, including the LSND evidence for oscillations. The higher degree of rejection for non-active solar and atmospheric oscillation solutions implied by the SNO neutral current result as well as by the latest 1489-day Super-K atmospheric neutrino data allows us to rule out (2+2) oscillation schemes proposed to reconcile LSND with the rest of current neutrino oscillation data. Using an improved goodness of fit (gof) method especially sensitive to the combination of data sets we obtain a gof of only 1.6 times 10^{-6} for (2+2) schemes. Further, we re-evaluate the status of (3+1) oscillations using two different analyses of the LSND data sample. We find that also (3+1) schemes are strongly disfavoured by the data. Depending on the LSND analysis we obtain a gof of 5.6 times 10^{-3} or 7.6 times 10^{-5}. This leads to the conclusion that all four-neutrino descriptions of the LSND anomaly, both in (2+2) as well as (3+1) realizations, are highly disfavoured. Our analysis brings the LSND hint to a more puzzling status.