The configuration of the three neutrino masses can take two forms, known as the normal and inverted hierarchies. We compute the Bayesian evidence associated with these two hierarchies. Previous studies found a mild preference for the normal hierarchy, and this was driven by the asymmetric manner in which cosmological data has confined the available parameter space. Here we identify the presence of a second asymmetry, which is imposed by data from neutrino oscillations. By combining constraints on the squared-mass splittings with the limit on the sum of neutrino masses of $\Sigma m_\nu < 0.13$ eV, and using a minimally informative prior on the masses, we infer odds of 42:1 in favour of the normal hierarchy, which is classified as "strong" in the Jeffreys' scale. We explore how these odds may evolve in light of higher precision cosmological data, and discuss the implications of this finding with regards to the nature of neutrinos. Finally the individual masses are inferred to be $m_1 = 3.80^{+26.2}_{-3.73} \, \text{meV}, m_2 = 8.8^{+18}_{-1.2} \, \text{meV}, m_3 = 50.4^{+5.8}_{-1.2} \, \text{meV}$ ($95\%$ credible intervals).

Cosmological observables show a dependence with the neutrino mass, which is partially degenerate with parameters of extended models of gravity. We study and explore this degeneracy in Horndeski generalized scalar-tensor theories of gravity. Using forecasted cosmic microwave background and galaxy power spectrum datasets, we find that a single parameter in the linear regime of the effective theory dominates the correlation with the total neutrino mass. For any given mass, a particular value of this parameter approximately cancels the power suppression due to the neutrino mass at a given redshift. The extent of the cancellation of this degeneracy depends on the cosmological large-scale structure data used at different redshifts. We constrain the parameters and functions of the effective gravity theory and determine the influence of gravity on the determination of the neutrino mass from present and future surveys.

We critically examine recent claims of a high solar metallicity by von Steiger \& Zurbuchen (2016) based on \textit{in situ} measurements of the solar wind, rather than the standard spectroscopically-inferred abundances (Asplund et al. 2009). We test the claim by Vagnozzi et al. (2016) that a composition based on the solar wind enables one to construct a standard solar model in agreement with helioseismological observations and thus solve the decades-old solar modelling problem. We show that, although some helioseismological observables are improved compared to models computed with spectroscopic abundances, most are in fact worse. The high abundance of refractory elements leads to an overproduction of neutrinos, with a predicted $^8$B flux that is nearly twice its observed value, and $^7$Be and CNO fluxes that are experimentally ruled out at high confidence. A combined likelihood analysis shows that models using the vSZ16 abundances fare worse than AGSS09 despite a higher metallicity. We also present astrophysical and spectroscopic arguments showing the vSZ16 composition to be an implausible representation of the solar interior, identifying the first ionisation potential effect in the outer solar atmosphere and wind as the likely culprit.

We demonstrate that a scalar field is unable to reverse its direction of motion while continuously exchanging energy with another fluid. If the rate of transfer is modulated by the scalar's acceleration, the field can undergo a rapid process of freezing, despite being displaced from the local minimum of its potential. This enables dark energy to form from any potential, regardless of its shape. The field's equation of state mimicks that of a cosmological constant. We present a physically motivated realisation in the form of a derivative neutrino-majoron coupling. Coherent motions, which form only once the neutrinos become non-relativistic, could be responsible for instigating the freezing process. This would provide a natural resolution to the dark energy coincidence problem, while avoiding the dynamical instabilities associated with mass-varying neutrino models. Finally we discuss possible means by which this model could be experimentally verified.

Recent constrains on the sum of neutrino masses inferred by analyzing cosmological data, show that detecting a non-zero neutrino mass is within reach of forthcoming cosmological surveys, implying a direct determination of the absolute neutrino mass scale. The measurement relies on constraining the shape of the matter power spectrum below the neutrino free streaming scale: massive neutrinos erase power at these scales. Detection of a lack of small-scale power, however, could also be due to a host of other effects. It is therefore of paramount importance to validate neutrinos as the source of power suppression at small scales. We show that, independent on hierarchy, neutrinos always show a footprint on large, linear scales; the exact location and properties can be related to the measured power suppression (an astrophysical measurement) and atmospheric neutrinos mass splitting (a neutrino oscillation experiment measurement). This feature can not be easily mimicked by systematic uncertainties or modifications in the cosmological model. The measurement of such a feature, up to 1% relative change in the power spectrum, is a smoking gun for confirming the determination of the absolute neutrino mass scale from cosmological observations. It also demonstrates the synergy of astrophysics and particle physics experiments.

We present an update of the determination of the solar neutrino fluxes from a global analysis of the solar and terrestrial neutrino data in the framework of three-neutrino mixing. Using a Bayesian analysis we reconstruct the posterior probability distribution function for the eight normalization parameters of the solar neutrino fluxes plus the relevant masses and mixing, with and without imposing the luminosity constraint. We then use these results to compare the description provided by different Standard Solar Models. Our results show that, at present, both models with low and high metallicity can describe the data with equivalent statistical agreement. We also argue that even with the present experimental precision the solar neutrino data have the potential to improve the accuracy of the solar model predictions.

The 511 keV electron-positron annihilation line, most recently characterized by the INTEGRAL/SPI experiment, is highly concentrated towards the Galactic centre. Its origin remains unknown despite decades of scrutiny. We propose a novel scenario in which known extragalactic positron sources such as radio jets of active galactic nuclei (AGN) fill the intergalactic medium with MeV e+e- pairs, which are then accreted into the Milky Way. We show that interpreting the diffuse cosmic radio background (CRB) as arising from radio sources with characteristics similar to the observed cores and radio lobes in powerful AGN jets suggests that the intergalactic positron-to-electron ratio could be as high as 10^{-5}, although this can be decreased if the CRB is not all produced by pairs and if not all positrons escape to the intergalactic medium. Assuming an accretion rate of one solar mass per year of matter into the Milky Way, a positron-to-electron ratio of only 10^{-7} is already enough to account for much of the 511 keV emission of the Galaxy. A simple spherical accretion model predicts an emission profile highly peaked in the central bulge, consistent with INTEGRAL observations. However, a realistic model of accretion with angular momentum would likely imply a more extended emission over the disk, with uncertainties depending on the magnetic field structure and turbulence in the galactic halo.

While standard solar model (SSM) predictions depend on approximately 20 input parameters, SSM neutrino flux predictions are strongly correlated with a single model output parameter, the core temperature $T_c$. Consequently, one can extract physics from solar neutrino flux measurements while minimizing the consequences of SSM uncertainties, by studying flux ratios with appropriate power-law weightings tuned to cancel this $T_c$ dependence. We re-examine an idea for constraining the primordial C+N content of the solar core from a ratio of CN-cycle $^{15}$O to pp-chain $^8$B neutrino fluxes, showing that nonnuclear SSM uncertainties in the ratio are small and effectively governed by a single parameter, the diffusion coefficient. We point out that measurements of both CN-I cycle neutrino branches -- $^{15}$O and $^{13}$N $\beta$-decay -- could in principle lead to separate determinations of the core C and N abundances, due to out-of-equilibrium CN-cycle burning in the cooler outer layers of the solar core. Finally, we show that the strategy of constructing "minimum uncertainty" neutrino flux ratios can also test other properties of the SSM. In particular, we demonstrate that a weighted ratio of $^7$Be and $^8$B fluxes constrains a product of S-factors to the same precision currently possible with laboratory data.

The study of solar neutrinos has given since ever a fundamental contribution both to astroparticle and to elementary particle physics, offering an ideal test of solar models and offering at the same time relevant indications on the fundamental interactions among particles. After reviewing the striking results of the last two decades, which were determinant to solve the long standing solar neutrino puzzle and refine the Standard Solar Model, we focus our attention on the more recent results in this field and on the experiments presently running or planned for the near future. The main focus at the moment is to improve the knowledge of the mass and mixing pattern and especially to study in detail the lowest energy part of the spectrum, which represents most of solar neutrino spectrum but is still a partially unexplored realm. We discuss this research project and the way in which present and future experiments could contribute to make the theoretical framemork more complete and stable, understanding the origin of some "anomalies" that seem to emerge from the data and contributing to answer some present questions, like the exact mechanism of the vacuum to matter transition and the solution of the so called solar metallicity problem.

We present a method to minimize, or even cancel out, the nuisance parameters affecting a measurement. Our approach is general and can be applied to any experiment or observation. We compare it with the bayesian technique used to deal with nuisance parameters: marginalization, and show how the method compares and improves by avoiding biases. We illustrate the method with several examples taken from the astrophysics and cosmology world: baryonic acoustic oscillations, cosmic clocks, Supernova Type Ia luminosity distance, neutrino oscillations and dark matter detection. By applying the method we recover some known results but also find some interesting new ones. For baryonic acoustic oscillation (BAO) experiments we show how to combine radial and angular BAO measurements in order to completely eliminate the dependence on the sound horizon at radiation drag. In the case of exploiting SN1a as standard candles we show how the uncertainty in the luminosity distance by a second parameter modeled as a metallicity dependence can be eliminated or greatly reduced. When using cosmic clocks to measure the expansion rate of the universe, we demonstrate how a particular combination of observables nearly removes the metallicity dependence of the galaxy on determining differential ages, thus removing the age-metallicity degeneracy in stellar populations. We hope that these findings will be useful in future surveys to obtain robust constraints on the dark energy equation of state.

We generate new standard solar models using newly analyzed nuclear fusion cross sections and present results for helioseismic quantities and solar neutrino fluxes. We discuss the status of the solar abundance problem and investigate whether nonstandard solar models with accretion from the protoplanetary disk might alleviate the problem. We examine a broad range of possibilities, analyzing both metal-enriched and metal-depleted accretion models and exploring three scenarios for the timing of the accretion. Only partial solutions are found: one can bring either the depth of the convective zone or the surface helium abundance into agreement with helioseismic results, but not both simultaneously. In addition, detailed results for solar neutrino fluxes show that neutrinos are a competitive source of information about the solar core and can help constrain possible accretion histories of the Sun. Finally, we briefly discuss how measurements of solar neutrinos from the CN-cycle could shed light on the interaction between the early Sun and its protoplanetary disk.

We explore how the Peccei-Quinn (PQ) axion parameter space can be constrained by the frequency-dependence dimming of radiation from astrophysical objects. To do so we perform accurate calculations of photon-axion conversion in the presence of a variable magnetic field. We propose several tests where the PQ axion parameter space can be explored with current and future astronomical surveys: the observed spectra of isolated neutron stars, occultations of background objects by white dwarfs and neutron stars, the light curves of eclipsing binaries containing a white dwarf. We find that the lack of dimming of the light curve of a detached eclipsing white dwarf binary recently observed, leads to relevant constraints on the photon-axion conversion. Current surveys designed for Earth-like planet searches are well matched to strengthen and improve the constraints on the PQ axion using astrophysical objects radiation dimming.

Cosmological probes are steadily reducing the total neutrino mass window, resulting in constraints on the neutrino-mass degeneracy as the most significant outcome. In this work we explore the discovery potential of cosmological probes to constrain the neutrino hierarchy, and point out some subtleties that could yield spurious claims of detection. This has an important implication for next generation of double beta decay experiments, that will be able to achieve a positive signal in the case of degenerate or inverted hierarchy of Majorana neutrinos. We find that cosmological experiments that nearly cover the whole sky could in principle distinguish the neutrino hierarchy by yielding 'substantial' evidence for one scenario over the another, via precise measurements of the shape of the matter power spectrum from large scale structure and weak gravitational lensing.

A class of well-behaved modified gravity models with long enough matter domination epoch and a late-time accelerated expansion is confronted with SNIa, CMB, SDSS, BAO and H(z) galaxy ages data, as well as current measurements of the linear growth of structure. We show that the combination of geometrical probes and growth data exploited here allows to rule out f(R) gravity models, in particular, the logarithmic of curvature model. We also apply solar system tests to the models in agreement with the cosmological data. We find that the exponential of the inverse of the curvature model satisfies all the observational tests considered and we derive the allowed range of parameters. Current data still allows for small deviations of Einstein gravity. Future, high precision growth data, in combination with expansion history data, will be able to distinguish tiny modifications of standard gravity from the LambdaCDM model.

Baryon Acoustic Oscillations (BAO) in the radial direction offer a method to directly measure the Universe expansion history, and to set limits to space curvature when combined to the angular BAO signal. In addition to spectroscopic surveys, radial BAO might be measured from accurate enough photometric redshifts obtained with narrow-band filters. We explore the requirements for a photometric survey using Luminous Red Galaxies (LRG) to competitively measure the radial BAO signal and discuss the possible systematic errors of this approach. If LRG were a highly homogeneous population, we show that the photo-z accuracy would not substantially improve by increasing the number of filters beyond $\sim 10$, except for a small fraction of the sources detected at high signal-to-noise, and broad-band filters would suffice to achieve the target $\sigma_z = 0.003 (1+z)$ for measuring radial BAO. Using the LRG spectra obtained from SDSS, we find that the spectral variability of LRG substantially worsens the achievable photometric redshift errors, and that the optimal system consists of $\sim$ 30 filters of width $\Delta \lambda / \lambda \sim 0.02$. A $S/N > 20$ is generally necessary at the filters on the red side of the $H\alpha$ break to reach the target photometric accuracy. We estimate that a 5-year survey in a dedicated telescope with etendue in excess of 60 ${\rm m}^2 {\rm deg}^2$ would be necessary to obtain a high enough density of galaxies to measure radial BAO with sufficiently low shot noise up to $z= 0.85$. We conclude that spectroscopic surveys have a superior performance than photometric ones for measuring BAO in the radial direction.

Standard solar models (SSM) are facing nowadays a new puzzle: the solar composition problem. New determinations of solar metal abundances lead SSM calculations to conflict with helioseismological measurements, showing discrepancies that extend from the convection zone to the solar core and can not be easily assigned to deficiencies in the modelling of the solar convection zone. We present updated solar neutrino fluxes and uncertainties for two SSM with high (old) and low (new) solar metallicity determinations. The uncertainties in iron and carbon abundances are the largest contribution to the uncertainties of the solar neutrino fluxes. The uncertainty on the ^14N+p -> ^15O+g rate is the largest of the non-composition uncertainties to the CNO neutrino fluxes. We propose an independent method to help identify which SSM is the correct one. Present neutrino data can not distinguish the solar neutrino predictions of both models but ongoing measurements can help to solve the puzzle.

In this work we study the phenomenological consequences of the dependence of mass varying neutrinos on the neutrino density in the Sun, which we precisely compute in each point along the neutrino trajectory. We find that a generic characteristic of these scenarios is that they establish a connection between the effective Delta m^2 in the Sun and the absolute neutrino mass scale. This does not lead to any new allowed region in the oscillation parameter space. On the contrary, due to this effect, the description of solar neutrino data worsens for large absolute mass. As a consequence a lower bound on the level of degeneracy can be derived from the combined analysis of the solar and KamLAND data. In particular this implies that the analysis favours normal over inverted mass orderings. These results, in combination with a positive independent determination of the absolute neutrino mass, can be used as a test of these scenarios together with a precise determination of the energy dependence of the survival probability of solar neutrinos, in particular for low energies.

We compare, using a three neutrino analysis, the allowed neutrino oscillation parameters and solar neutrino fluxes determined by the experimental data available Before and After Neutrino 2004. New data available after Neutrino2004 include refined KamLAND and gallium measurements. We use six different approaches to analyzing the KamLAND data. We present detailed results using all the available neutrino and anti-neutrino data for Delta m^2_{12}, tan^2 theta_{12}, sin^2 theta_{13}, and sin^2 eta (sterile fraction). Using the same complete data sets, we also present Before and After determinations of all the solar neutrino fluxes, which are treated as free parameters, an upper limit to the luminosity fraction associated with CNO neutrinos, and the predicted rate for a 7Be solar neutrino experiment. The 1 sigma (3 sigma) allowed range of Delta m^2_{21} = (8.2 +- 0.3) (^+1.0_-0.8)times 10^{-5} eV^2 is decreased by a factor of 1.7 (5), but the allowed ranges of all other neutrino oscillation parameters and neutrino fluxes are not significantly changed. Maximal mixing is disfavored at 5.8 sigma and the bound on the mixing angle theta_{13} is slightly improved to sin^2 theta_{13}<0.048 at 3 sigma. The predicted rate in a 7Be neutrino-electron scattering experiment is (0.665 +-0.015) of the rate implied by the BP04 solar model in the absence of neutrino oscillations. The corresponding predictions for p-p and pep experiments are, respectively, 0.707 {+0.011}{-0.013} and 0.644 {+0.011}{-0.013}. We derive upper limits to CPT violation in the weak sector by comparing reactor anti-neutrino oscillation parameters with neutrino oscillation parameters. We also show that the recent data disfavor at 91 % CL a proposed non-standard interaction description of solar neutrino oscillations.

We assess how well next generation neutrinoless double beta decay and normal neutrino beta decay experiments can answer four fundamental questions. 1) If neutrinoless double beta decay searches do not detect a signal, and if the spectrum is known to be inverted hierarchy, can we conclude that neutrinos are Dirac particles? 2) If neutrinoless double beta decay searches are negative and a next generation ordinary beta decay experiment detects the neutrino mass scale, can we conclude that neutrinos are Dirac particles? 3) If neutrinoless double beta decay is observed with a large neutrino mass element, what is the total mass in neutrinos? 4) If neutrinoless double beta decay is observed but next generation beta decay searches for a neutrino mass only set a mass upper limit, can we establish whether the mass hierarchy is normal or inverted? We base our answers on the expected performance of next generation neutrinoless double beta decay experiments and on simulations of the accuracy of calculations of nuclear matrix elements.

We analyze all available solar and related reactor neutrino experiments, as well as simulated future 7Be, p-p, pep, and ^8B solar neutrino experiments. We treat all solar neutrino fluxes as free parameters subject to the condition that the total luminosity represented by the neutrinos equals the observed solar luminosity (the `luminosity constraint'). Existing experiments show that the p-p solar neutrino flux is 1.02 +- 0.02 (1 sigma) times the flux predicted by the BP00 standard solar model; the 7Be neutrino flux is 0.93^{+0.25}_{-0.63} the predicted flux; and the ^8B flux is 1.01 +- 0.04 the predicted flux. The neutrino oscillation parameters are: Delta m^2 = 7.3^{+0.4}_{-0.6}\times 10^{-5} eV^2 and tan^2 theta_{12} = 0.41 +- 0.04. We evaluate how accurate future experiments must be to determine more precisely neutrino oscillation parameters and solar neutrino fluxes, and to elucidate the transition from vacuum-dominated to matter-dominated oscillations at low energies. A future 7Be nu-e scattering experiment accurate to +- 10 % can reduce the uncertainty in the experimentally determined 7Be neutrino flux by a factor of four and the uncertainty in the p-p neutrino flux by a factor of 2.5 (to +- 0.8 %). A future p-p experiment must be accurate to better than +- 3 % to shrink the uncertainty in tan^2 theta_{12} by more than 15 %. The idea that the Sun shines because of nuclear fusion reactions can be tested accurately by comparing the observed photon luminosity of the Sun with the luminosity inferred from measurements of solar neutrino fluxes. Based upon quantitative analyses of present and simulated future experiments, we answer the question: Why perform low-energy solar neutrino experiments?

We use the recently reported KamLAND measurements on oscillations of reactor anti-neutrinos, together with the data of previously reported solar neutrino experiments, to show that: (1) the total 8B neutrino flux emitted by the Sun is 1.00(1.0 \pm 0.06) of the standard solar model (BP00) predicted flux, (2) the KamLAND measurements reduce the area of the globally allowed oscillation regions that must be explored in model fitting by six orders of magnitude in the Delta m^2-tan^2 theta plane, (3) LMA is now the unique oscillation solution to a CL of 4.7sigma, (4) maximal mixing is disfavored at 3.1 sigma, (5) active-sterile admixtures are constrained to sin^2 eta<0.13 at 1 sigma, (6) the observed ^8B flux that is in the form of sterile neutrinos is 0.00^{+0.09}_{-0.00} (1 sigma), of the standard solar model (BP00) predicted flux, and (7) non-standard solar models that were invented to completely avoid solar neutrino oscillations are excluded by KamLAND plus solar at 7.9 sigma . We also refine quantitative predictions for future 7Be and p-p solar neutrino experiments.

The 8B solar neutrino flux inferred from a global analysis of solar neutrino experiments is within 11% (1 sigma) of the predicted standard solar model value if only active neutrinos exist, but could be as large as 1.7 times the standard prediction if sterile neutrinos exist. We show that the total 8B neutrino flux (active plus sterile neutrinos) can be determined experimentally to about 10% (1 sigma) by combining charged current measurements made with the KamLAND reactor experiment and with the SNO CC solar neutrino experiment, provided the LMA neutrino oscillation solution is correct and the simulated performance of KamLAND is valid. Including also SNO NC data, the sterile component of the 8B neutrino flux can be measured by this method to an accuracy of about 12% (1 sigma) of the standard solar model flux. Combining Super-Kamiokande and KamLAND measurements and assuming the oscillations occur only among active neutrinos, the 8B neutrino flux can be measured to 6% (1 sigma); the total flux can be measured to an accuracy of about 9%. The total 7Be solar neutrino flux can be determined to an accuracy of about 28% (1 sigma) by combining measurements made with the KamLAND, SNO, and gallium neutrino experiments. One can determine the total 7Be neutrino flux to a one sigma accuracy of about 11% or better by comparing data from the KamLAND experiment and the BOREXINO solar neutrino experiment provided both detectors work as expected. The pp neutrino flux can be determined to about 15% using data from the gallium, KamLAND, BOREXINO, and SNO experiments.

With the goal of identifying signatures that select specific neutrino oscillation parameters, we test the robustness of global oscillation solutions that fit all the available solar and reactor experimental data. We use three global analysis strategies previously applied by different authors and also determine the sensitivity of the oscillation solutions to the critical nuclear fusion cross section, S_{17}(0), for the production of 8B. The favored solutions are LMA, LOW, and VAC in order of g.o.f. The neutral current to charged current ratio for SNO is predicted to be 3.5 +- 0.6 (1 sigma), which is separated from the no-oscillation value of 1.0 by much more than the expected experimental error. The predicted range of the day-night difference in charged current rates is (8.2 +- 5.2)% and is strongly correlated with the day-night effect for neutrino-electron scattering. A measurement by SNO of either a NC to CC ratio > 3.3 or a day-night difference > 10%, would favor a small region of the currently allowed LMA neutrino parameter space. The global oscillation solutions predict a 7Be neutrino-electron scattering rate in BOREXINO and KamLAND in the range 0.66 +- 0.04 of the BP00 standard solar model rate, a prediction which can be used to test both the solar model and the neutrino oscillation theory. Only the LOW solution predicts a large day-night effect(< 42%) in BOREXINO and KamLAND. For the KamLAND reactor experiment, the LMA solution predicts 0.44 of the standard model rate; we evaluate 1 sigma and 3 sigma uncertainties and the first and second moments of the energy spectrum.

For active and sterile neutrinos, we present the globally allowed solutions for two neutrino oscillations. We include the SNO CC measurement and all other relevant solar neutrino and reactor data. Five active neutrino oscillation solutions (LMA, LOW, SMA, VAC, and Just So2) are currently allowed at 3 sigma; three sterile neutrino solutions (Just So2, SMA, and VAC) are allowed at 3 sigma. The goodness of fit is satisfactory for all eight solutions. We also investigate the robustness of the allowed solutions by carrying out global analyses with and without: 1) imposing solar model constraints on the 8B neutrino flux, 2) including the Super-Kamiokande spectral energy distribution and day-night data, 3) including a continuous mixture of active and sterile neutrinos, 4) using an enhanced CC cross section for deuterium (due to radiative corrections), and 5) a optimistic, hypothetical reduction by a factor of three of the error of the SNO CC rate. For every analysis strategy used in this paper, the most favored solutions all involve large mixing angles: LMA, LOW, or VAC. The favored solutions are robust, but the presence at 3 sigma of individual sterile solutions and the active Just So2 solution is sensitive to the analysis assumptions.