This white paper summarizes the activities of the Brazilian community concerning dark matter physics and highlights the importance of financial support to Brazilian groups that are deeply involved in experimental endeavours. The flagships of the Brazilian dark matter program are the Cherenkov Telescope Array, DARKSIDE, SBN and LHC experiments, but we emphasize that smaller experiments such as DAMIC and CONNIE constitute important probes to dark sectors as well and should receive special attention. Small experimental projects showing the potential to probe new regions of parameter space of dark matter models are encouraged. On the theoretical and phenomenological side, some groups are devoted to astrophysical aspects such as the dark matter density profile while others explore the signature of dark matter models at colliders, direct and indirect detection experiments. In summary, the Brazilian dark matter community that was born not long ago has grown tremendously in the past years and now plays an important role in the hunt for a dark matter particle.

Flavour oscillations experiments are suggesting the existence of a sterile, $4$th neutrino's generation with a mass of an eV order. This would mean an additional relativistic degree of freedom in the cosmic inventory, in contradiction with recent results from the Planck satellite, that have confirmed the standard value $N_{eff} \approx 3$ for the effective number of relativistic species. On the other hand, the Planck best-fit for the Hubble-Lema\^itre parameter is in tension with the local value determined with the Hubble Space Telescope, and adjusting $N_{eff}$ is a possible way to overcome such a tension. In this paper we perform a joint analysis of three complementary cosmological distance rulers, namely the CMB acoustic scale measured by Planck, the BAO scale model-independently determined by Verde {\it et al.}, and luminosity distances measured with JLA and Pantheon SNe Ia surveys. Two Gaussian priors were imposed to the analysis, the local expansion rate measured by Riess {\it et al.}, and the baryon density parameter fixed from primordial nucleosynthesis by Cooke {\it et al.}. For the sake of generality, two different models are used in the tests, the standard $\Lambda$CDM model and a generalised Chaplygin gas. The best-fit gives $N_{eff} \approx 4$ in both models, with a Chaplygin gas parameter slightly negative, $\alpha \approx -0.04$. The standard value $N_{eff} \approx 3$ is ruled out with $\approx 3\sigma$.

We consider the possibility of solar neutrino decay as a sub-leading effect on their propagation between production and detection. Using current oscillation data, we set a new lower bound to the $\nu_2$ neutrino lifetime at $\tau_2\, /\, m_2 \geq 7.2 \times 10^{-4}\,\,\hbox{s}\,.\,\hbox{eV}^{-1}$ at $99\%\,$C.L.. Also, we show how seasonal variations in the solar neutrino data can give interesting additional information about neutrino lifetime.

The recent high-statistics high-energy atmospheric neutrino data collected by IceCube open a new window to probe new physics scenarios that are suppressed in lower energy neutrino experiments. In this paper we analyze the IceCube atmospheric neutrino data to constrain the Violation of Equivalence Principle (VEP) in the framework of three neutrinos with non-universal gravitational couplings. In this scenario the effect of VEP on neutrino oscillation probabilities can be parametrized by two parameters $\Delta \gamma_{21}\equiv \gamma_2-\gamma_1$ and $\Delta\gamma_{31}\equiv \gamma_3-\gamma_1$, where $\gamma_i$'s denote the coupling of neutrino mass eigenstates to gravitational field. By analyzing the latest muon-tracks data sets of IceCube-40 and IceCube-79, besides providing the 2D allowed regions in $(\phi\Delta\gamma_{21},\phi\Delta\gamma_{31})$ plane, we obtain the upper limits $|\phi\Delta\gamma_{21}| < 9.1\times 10^{-27}$ (at 90\% C.L.) which improves the previous limit by $\sim4$ orders of magnitude and $|\phi\Delta\gamma_{31}| \lesssim 6\times 10^{-27}$ (at 90\% C.L.) which improves the current limit by $\sim1$ order of magnitude. Also we discuss in detail and analytically the effect of VEP on neutrino oscillation probabilities.

Structure formation creates high temperature and density regions in the Universe that allow the conversion of matter into more stable states, with a corresponding emission of relativistic matter and radiation. An example of such a mechanism is the supernova event, that releases relativistic neutrinos corresponding to 99% of the binding energy of remnant neutron star. We take this phenomena as a starting point for an assumption that similar processes could occur in the dark sector, where structure formation would generate a late time conversion of cold dark matter into a relativistic form of dark matter. We performed a phenomenological study about the limits of this conversion, where we assumed a transition profile that is a generalized version of the neutrino production in supernovae events. With this assumption, we obtained an interesting modification for the constraint over the cold dark matter density. We show that when comparing with the standard \Lambda CDM cosmology, there is no preference for conversion, although the best fit is within 1\sigma\ from the standard model best fit. The methodology and the results obtained qualify this conversion hypothesis, from the large scale structure point of view, as a viable and interesting model to be tested in the future with small scale data, and mitigate discrepancies between observations at this scale and the pure cold dark matter model.

The relatively large value of neutrino mixing angle \theta_{13} set by recent measurements allows us to use solar neutrinos to set a limit on neutrino magnetic moment involving second and third families, \mu_{\mu\tau}. The existence of a random magnetic field in solar convective zone can produce a significant anti-neutrino flux when a non-vanishing neutrino magnetic moment is assumed. Even if we consider a vanishing neutrino magnetic moment involving the first family, electron anti-neutrinos are indirectly produced through the mixing between first and third families and non-vanishing \mu_{\mu\tau}. Using KamLAND limits on the solar flux of electron anti-neutrino, we set the limit \mu_{\mu\tau} < 0.5e-11 Bohr magneton for a reasonable assumption on the behavior of solar magnetic fields. This is the first time a limit on \mu_{\mu\tau} is established in the literature directly from neutrino interaction with magnetic fields, and, interestingly enough, is comparable with the limits on neutrino magnetic moment involving the first family and with the ones coming from modifications on electroweak cross section.

Recent results from the SNO, Super-Kamiokande and Borexino experiments do not show the expected upturn of the energy spectrum of events (the ratio $R \equiv N_{obs}/N_{SSM}$) at low energies. At the same time, cosmological observations testify for possible existence of additional relativistic degrees of freedom in the early Universe: $\Delta N_{eff} = 1 - 2$. These facts strengthen the case of very light sterile neutrino, $\nu_s$, with $\Delta m^2_{01} \sim (0.7 - 2) \cdot 10^{-5}$ eV$^2$, which mixes weakly with the active neutrinos. The $\nu_s$ mixing in the mass eigenstate $\nu_1$ characterized by $\sin^2 2\alpha \sim 10^{-3}$ can explain an absence of the upturn. The mixing of $\nu_s$ in the eigenstate $\nu_3$ with $\sin^2 \beta \sim 0.1$ leads to production of $\nu_s$ via oscillations in the Universe and to additional contribution $\Delta N_{eff} \approx 0.7 - 1$ before the big bang nucleosynthesis and later. Such a mixing can be tested in forthcoming experiments with the atmospheric neutrinos as well as in future accelerator long baseline experiments. It has substantial impact on conversion of the supernova neutrinos.

Solar neutrino physics enters a stage of precision measurements. In this connection we present a precise analytic description of the neutrino conversion in the context of LMA MSW solution of the solar neutrino problem. Using the adiabatic perturbation theory we derive an analytic formula for the $\nu_e$ survival probability which takes into account the non-adiabatic corrections and the regeneration effect inside the Earth. The probability is averaged over the neutrino production region. We find that the non-adiabatic corrections are of the order $10^{-9}-10^{-7}$. Using the formula for the Earth regeneration effect we discuss features of the zenith angle dependence of the $\nu_e$ flux. In particular, we show that effects of small structures at the surface of the Earth can be important.

The Homestake result is about ~ 2 \sigma lower than the Ar-production rate, Q_{Ar}, predicted by the LMA MSW solution of the solar neutrino problem. Also there is no apparent upturn of the energy spectrum (R \equiv N_{obs}/N_{SSM}) at low energies in SNO and Super-Kamiokande. Both these facts can be explained if a light, \Delta m^2_{01} ~ (0.2 - 2) \cdot 10^{-5} eV^2, sterile neutrino exists which mixes very weakly with active neutrinos: \sin^2 2\alpha ~ (10^{-5} - 10^{-3}). We perform both the analytical and numerical study of the conversion effects in the system of two active neutrinos with the LMA parameters and one weakly mixed sterile neutrino. The presence of sterile neutrino leads to a dip in the survival probability in the intermediate energy range E = (0.5 - 5) MeV thus suppressing the Be, or/and pep, CNO as well as B electron neutrino fluxes. Apart from diminishing Q_{Ar} it leads to decrease of the Ge-production rate and may lead to decrease of the BOREXINO signal and CC/NC ratio at SNO. Future studies of the solar neutrinos by SNO, SK, BOREXINO and KamLAND as well as by the new low energy experiments will allow us to check this possibility. We present a general analysis of modifications of the LMA energy profile due to mixing with new neutrino states.

We show that the non-standard neutrino interactions can play a role as sub-leading effect on the solar neutrino oscillations. We observe that very small flavor universality violations of order of 0.1-0.2 G_F is sufficient to induce two phenomena: suppression of the eletronic neutrino earth regeneration and a shift of the resonance layer in the sun. We obtain these phenomena even in the absence of any flavor changing interactions. We discuss their consequences and confront with a global analysis of solar+KamLAND results. We conclude that a new compatibility region in the Delta m^2 x tan^2\theta, which we call very low Large Mixing Angle region, is found for Delta m^2 ~ 1e-5 eV^2 and tan^2\theta= 0.45.

We have performed analysis of the solar neutrino data including results from the SNO salt phase as well as the combined analysis of the solar and the KamLAND results. The best fit values of neutrino parameters are Delta m^2 = 7.1e-5 eV^2, tan^2\theta = 0.40 with the boron flux f_B = 1.04. New SNO results strongly disfavor maximal mixing and the h-LMA region (Delta m^2 > 1e-4 eV^2) which is accepted now at the 3-sigma level. We find the 3-sigma upper bounds: Delta m^2 < 1.7e-4$ eV^2 and tan^2\theta < 0.64, and the lower bound Delta m^2 > 4.8e-5 eV^2. Non-zero 13-mixing does not change these results significantly. The present data determine quantitatively the physical picture of the solar neutrino conversion. At high energies relevant for SNO and Super-Kamiokande the deviation of the effective survival probability from the non-oscillatory value is about 10 - 14%. The oscillation effect contribution to this difference about 10% and the Earth regeneration is about 3 - 4%. At low energies (E < 1 MeV) the matter corrections to vacuum oscillation effect are below 5%. The predictions for the forthcoming measurements are given which include the spectral distortion and CC/NC ratio at SNO, the Day-Night asymmetry, the KamLAND spectrum and rate.

There are reasons to believe that mechanisms exist in the solar interior which lead to random density perturbations in the resonant region of the Large Mixing Angle solution to the solar neutrino problem. We find that, in the presence of these density perturbations, the best fit point in the (sin^2(2\theta), Delta_m^2) parameter space moves to smaller values, compared with the values obtained for the standard LMA solution. Combining solar data with KamLAND results, we find a new compatibility region, which we call VERY-LOW LMA, where sin^2(2\theta) ~ 0.6 and Delta_m^2~2e-5 eV^2, for random density fluctuations of order 5% < \xi< 8%. We argue that such values of density fluctuations are still allowed by helioseismological observations at small scales of order 10 - 1000 km deep inside the solar core.

The first KamLAND results are in a very good agreement with the predictions made on the basis of the solar neutrino data and the LMA realization of the MSW mechanism. We perform a combined analysis of the KamLAND (rate, spectrum) and the solar neutrino data with a free boron neutrino flux f_B. The best fit values of neutrino parameters are Delta m^2 = 7.3e-5 eV^2, tg^2 theta = 0.41 and f_B = 1.05 with the 1 sigma intervals: Delta m^2 = (6.2 - 8.4)e-5 eV^2, tg^2 theta = 0.33 - 0.54. We find the 3 sigma upper bounds: Delta m^2 < 2.8e-4 eV^2 and tg^2 theta < 0.84, and the lower bound Delta m^2 > 4e-5 eV^2. At 99% C.L. the KamLAND spectral result splits the LMA region into two parts with the preferred one at Delta m^2 < 1e-4 eV^2. The higher Delta m^2 region is accepted at about 2 sigma level. We show that effects of non-zero 13-mixing, sin^2 theta_{13} < 0.04, are small leading to a slight improvement of the fit in higher Delta m^2 region. In the best fit point we predict for SNO: CC/NC = 0.33 +0.05-0.03 and A_{DN}(SNO) = 2.8+-0.8 % (68% C.L.), and A_{DN}(SNO) < 9 % at the 3 sigma level. Further improvements in the determination of the parameters are discussed and implications of the solar neutrino and KamLAND results are considered.

What can we learn from solar neutrino observations? Is there any solution to the solar neutrino anomaly which is favored by the present experimental panorama? After SNO results, is it possible to affirm that neutrinos have mass? In order to answer such questions we analyze the current available data from the solar neutrino experiments, including the recent SNO result, in view of many acceptable solutions to the solar neutrino problem based on different conversion mechanisms, for the first time, using the same statistical procedure. This allows us to do a direct comparison of the goodness of the fit among different solutions, from which we can discuss and conclude on the current status of each proposed dynamical mechanism. These solutions are based on different assumptions: (a) neutrino mass and mixing, (b) non-vanishing neutrino magnetic moment, (c) the existence of non-standard flavor-changing and non-universal neutrino interactions and (d) the tiny violation of the equivalence principle. We investigate the quality of the fit provided by each one of these solutions not only to the total rate measured by all the solar neutrino experiments but also to the recoil electron energy spectrum measured at different zenith angles by the Super-Kamiokande collaboration. We conclude that several non-standard neutrino flavor conversion mechanisms provide a very good fit to the experimental data which is comparable with (or even slightly better than) the most famous solution to the solar neutrino anomaly based on the neutrino oscillation induced by mass.

We present an updated global analysis of two-flavor MSW solutions to the solar neutrino problem in terms of conversions of nu_e into active or sterile neutrinos. We perform a fit to the full data set corresponding to the 825-day Super-Kamiokande data sample as well as to Chlorine, GALLEX and SAGE experiments. We use all measured total event rates as well as Super-Kamiokande data on the zenith angle dependence, energy spectrum and seasonal variation of the events. For conversions into active neutrinos we find that, although the data on the total event rates favours the Small Mixing Angle (SMA) solution, once the full data set is included both SMA and Large Mixing Angle (LMA) solutions give an equivalent fit to the data. The best-fit points for the combined analysis are Delta m^2=3.6 10^{-5} eV^2 sin^2(2theta)=0.79 with chi^2_min=35.4/30 d.o.f and Delta m^2=5.1~ 10^{-6} eV^2 and sin^2(2theta)=5.5 10^{-3} with chi^2_min=37.4/30 d.o.f. In contrast with the earlier 504-day study of Bahcall-Krastev-Smirnov our results indicate that the LMA solution is not only allowed, but slightly preferred. We also show that seasonal effects, although small, may still reach 8.5 % at the best-fit LMA solution without conflict with the hints of a day-night variation (6% is due to the eccentricity of the Earth's orbit). For conversions into sterile neutrinos only the SMA solution is allowed with best-fit point Delta m^2=5.0 10^{-6} eV^2 and sin^2(2theta)=3. 10^{-3} and chi^2_min=40.2/30 d.o.f. We also consider departures of the Standard Solar Model of Bahcall and Pinsonneault 1998 by allowing arbitrary 8^B and hep fluxes. Best fit is obtained for ^8B/^8B_SSM=0.61 (1.37) and hep/hep_SSM=12 (38) for the SMA (LMA) solution.

MSW solutions of the solar neutrino problem predict a seasonal dependence of the zenith angle distribution of the event rates, due to the non-zero latitude at the Super-Kamiokande site. We calculate this seasonal dependence and compare it with the expectations in the no-oscillation case as well as just-so scenario, in the light of the latest Super-Kamiokande 708-day data. The seasonal dependence can be sizeable in the large mixing angle MSW solution and would be correlated with the day-night effect. This may be used to discriminate between MSW and just-so scenarios and should be taken into account in refined fits of the data.