The density of dark matter near the Sun, ρ_{DM, ⊙}, is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way’s dark matter halo. Estimates to date have typically assumed that the Milky Way’s stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. We assess the impact of these assumptions on determinations of ρ_{DM, ⊙} by applying a free-form, steady-state, Jeans method to two different N-body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia–Sausage–Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans–Poisson equations on our ability to correctly extract ρ_{DM, ⊙} from the simulated data. We find that common approximations employed in the literature – axisymmetry and a locally flat rotation curve – can lead to significant systematic errors of up to a factor ∼1.5 in the recovered surface mass density ∼2 kpc above the disc plane, implying a fractional error on ρ_{DM, ⊙} of the order of unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate of ρ_{DM, ⊙}, consistent with the true value within our 95 per cent confidence intervals for realistic 20 per cent uncertainties on the baryonic surface density of the disc. Other terms – the axial tilt, 2nd Poisson and time-dependent terms – contribute less than 10 per cent to ρ_{DM, ⊙} (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis.

This paper presents the main features of a new and updated version of the program PArthENoPE, which the community has been using for many years for computing the abundances of light elements produced during Big Bang Nucleosynthesis. This is the third release of the PArthENoPE code, after the 2008 and the 2018 ones, and will be distributed from the code's website, http://parthenope.na.infn.it. Apart from minor changes, the main improvements in this new version include a revisited implementation of the nuclear rates for the most important reactions of deuterium destruction, 2H(p,γ)3He, 2H(d, n)3He and 2H(d, p)3H, and a re-designed GUI, which extends the functionality of the previous one. The new GUI, in particular, supersedes the previous tools for running over grids of parameters with a better management of parallel runs, and it offers a brand-new set of functions for plotting the results.

We present an updated global fit of neutrino oscillation data in the simplest three-neutrino framework. In the present study we include up-to-date analyses from a number of experiments. Concerning the atmospheric and solar sectors, besides the data considered previously, we give updated analyses of IceCube DeepCore and Sudbury Neutrino Observatory data, respectively. We have also included the latest electron antineutrino data collected by the Daya Bay and RENO reactor experiments, and the long-baseline T2K and NO nu A measurements, as reported in the Neutrino 2020 conference. All in all, these new analyses result in more accurate measurements of theta (13), theta (12), Delta m212 and Delta m312. The best fit value for the atmospheric angle theta (23) lies in the second octant, but first octant solutions remain allowed at similar to 2.4 sigma. Regarding CP violation measurements, the preferred value of delta we obtain is 1.08 pi (1.58 pi) for normal (inverted) neutrino mass ordering. The global analysis still prefers normal neutrino mass ordering with 2.5 sigma statistical significance. This preference is milder than the one found in previous global analyses. These new results should be regarded as robust due to the agreement found between our Bayesian and frequentist approaches. Taking into account only oscillation data, there is a weak/moderate preference for the normal neutrino mass ordering of 2.00 sigma. While adding neutrinoless double beta decay from the latest Gerda, CUORE and KamLAND-Zen results barely modifies this picture, cosmological measurements raise the preference to 2.68 sigma within a conservative approach. A more aggressive data set combination of cosmological observations leads to a similar preference for normal with respect to inverted mass ordering, namely 2.70 sigma. This very same cosmological data set provides 2 sigma upper limits on the total neutrino mass corresponding to Sigma m(nu)< 0.12 (0.15) eV in the normal (inverted) neutrino mass ordering scenario. The bounds on the neutrino mixing parameters and masses presented in this up-to-date global fit analysis include all currently available neutrino physics inputs.

We present a consistent framework to set limits on properties of light sterile neutrinos coupled to all three active neutrinos using a combination of the latest cosmological data and terrestrial measurements from oscillations, β-decay, and neutrinoless double-β-decay (0νββ) experiments. We directly constrain the full 3+1 active-sterile mixing matrix elements |Uα4|2, with α∈(e,μ,τ), and the mass-squared splitting Δm241≡m24−m21. We find that results for a 3+1 case differ from previously studied 1+1 scenarios where the sterile is coupled to only one of the neutrinos, which is largely explained by parameter space volume effects. Limits on the mass splitting and the mixing matrix elements are currently dominated by the cosmological datasets. The exact results are slightly prior dependent, but we reliably find all matrix elements to be constrained below |Uα4|2≲10−3. Short-baseline neutrino oscillation hints in favor of eV-scale sterile neutrinos are in serious tension with these bounds, irrespective of prior assumptions. We also translate the bounds from the cosmological analysis into constraints on the parameters probed by laboratory searches, such as mβ or mββ, the effective mass parameters probed by β-decay and 0νββ searches, respectively. When allowing for mixing with a light sterile neutrino, cosmology leads to upper bounds of mβ<0.09  eV and mββ<0.07  eV at 95% CL, more stringent than the limits from current laboratory experiments.

Neutrino non-standard interactions (NSI) with electrons are known to alter the picture of neutrino de coupling from the cosmic plasma. NSI modify both flavour oscillations through matter effects, and the annihilation and scattering between neutrinos and electrons and positrons in the thermal plasma. In view of the forthcoming cosmological observations, we perform a precision study of the impact of non universal and flavour-changing NSI on the effective number of neutrinos, N_eff. We present the variation of N_eff arising from the different NSI parameters and discuss the existing degeneracies among them, from cosmology alone and in relation to the current bounds from terrestrial experiments. Even though cosmology is generally less sensitive to NSI than these experiments, we find that future cosmological data would provide competitive and complementary constraints for some of the couplings and their combinations.

This report summarises progress made in estimating the local density of dark matter (ρ_DM,⊙), a quantity that is especially important for dark matter direct detection experiments. We outline and compare the most common methods to estimate ρ_DM,⊙ and the results from recent studies, including those that have benefited from the observations of the ESA/Gaia satellite. The result of most local analyses coincide within a range of , while a slightly lower range of is preferred by most global studies. In light of recent discoveries, we discuss the importance of going beyond the approximations of what we define as the ideal Galaxy (a steady-state Galaxy with axisymmetric shape and a mirror symmetry across the mid-plane) in order to improve the precision of ρ_DM,⊙ measurements. In particular, we review the growing evidence for local disequilibrium and broken symmetries in the present configuration of the Milky Way, as well as uncertainties associated with the galactic distribution of baryons. Finally, we comment on new ideas that have been proposed to further constrain the value of ρ_DM,⊙, most of which would benefit from Gaia's final data release.

We present in this work a new calculation of the standard-model benchmark value for the effective number of neutrinos, N-eff(SM), that quantifies the cosmological neutrinoto-photon energy densities. The calculation takes into account neutrino flavour oscillations, finite-temperature effects in the quantum electrodynamics plasma to O(e(3)), where e is the elementary electric charge, and a full evaluation of the neutrino-neutrino collision integral. We provide furthermore a detailed assessment of the uncertainties in the benchmark N(eff)(SM )value, through testing the value's dependence on (i) optional approximate modelling of the weak collision integrals, (ii) measurement errors in the physical parameters of the weak sector, and (iii) numerical convergence, particularly in relation to momentum discretisation. Our new, recommended standard-model benchmark is N-eff(SM) 3.0440 +/- 0.0002, where the nominal uncertainty is attributed predominantly to errors incurred in the numerical solution procedure (vertical bar delta N-eff vertical bar similar to 10(-4)), augmented by measurement errors in the solar mixing angle sin(2) theta(12) (vertical bar delta N-eff vertical bar similar to 10(-4)).

Aims. We infer the gravitational potential of the Galactic disk by analysing the phase-space densities of 120 stellar samples in 40 spatially separate sub-regions of the solar neighbourhood, using Gaia's second data release (DR2), in order to quantify spatially dependent systematic effects that bias this type of measurement.

Methods. The gravitational potential was inferred under the assumption of a steady state in the framework of a Bayesian hierarchical model. We performed a joint fit of our stellar tracers' three-dimensional velocity distribution, while fully accounting for the astrometric uncertainties of all stars as well as dust extinction, and we also masked angular areas of known open clusters. The inferred gravitational potential is compared, post-inference, to a model for the baryonic matter and halo dark matter components.

Results. We see an unexpected but clear trend for all 40 spatially separate sub-regions: Compared to the potential derived from the baryonic model, the inferred gravitational potential is significantly steeper close to the Galactic mid-plane (less than or similar to 60 pc), but flattens such that the two agree well at greater distances (similar to 400 pc). The inferred potential implies a total matter density distribution that is highly concentrated to the Galactic mid-plane and decays quickly with height. We see a dependence on the Galactic radius that is consistent with a disk scale length of a few kiloparsecs. Apart from this, there are discrepancies between stellar samples, implying spatially dependent systematic effects which are, at least in part, explained by substructures in the phase-space distributions.

Conclusions. In terms of the inferred matter density distribution, the very low matter density that is inferred at greater heights (greater than or similar to 300 pc) is inconsistent with the observed scale height and matter distribution of the stellar disk, which cannot be explained by a misunderstood density of cold gas or other hidden mass. Our interpretation is that these results must be biased by a time-varying phase-space structure, possibly a breathing mode, that is large enough to affect all stellar samples in the same manner.

We present a new method for inferring the gravitational potential of the Galactic disk, using the time-varying structure of a phase-space spiral in the (z, w)-plane (where z and w represent vertical position and vertical velocity). Our method of inference extracts information from the shape of the spiral and disregards the bulk density distribution that is usually used to perform dynamical mass measurements. In this manner, it is complementary to traditional methods that are based on the assumption of a steady state. Our method consists of fitting an analytical model for the phase-space spiral to data, where the spiral is seen as a perturbation of the stellar number density in the (z, w)-plane. We tested our method on one-dimensional simulations, which were initiated in a steady state and then perturbed by an external force similar to that of a passing satellite. We were able to retrieve the true gravitational potentials of the simulations with high accuracy. The gravitational potential at 400-500 parsec distances from the disk mid-plane was inferred with an error of only a few percent. This is the first paper of a series in which we plan to test and refine our method on more complex simulations, as well as apply our method to Gaia data.

Using the method that was developed in the first paper of this series, we measured the vertical gravitational potential of the Galactic disk from the time-varying structure of the phase-space spiral, using data from Gaia as well as supplementary radial velocity information from legacy spectroscopic surveys. For eleven independent data samples, we inferred gravitational potentials that were in good agreement, despite the data samples' varied and substantial selection e ffects. Using a model for the baryonic matter densities, we inferred a local halo dark matter density of 0.0085 +/- 0.0039 M(circle dot)pc(-3) = 0.32 +/- 0.15 GeV cm(-3). We were also able to place the most stringent constraint on the surface density of a thin dark disk with a scale height <= 50 pc, corresponding to an upper 95% confidence limit of roughly 5 M(circle dot)pc(-2) (compared to the previous limit of roughly 10 M(circle dot)pc(-2), given the same scale height). For the inferred halo dark matter density and thin dark disk surface density, the statistical uncertainties are dominated by the baryonic model, which potentially could also su ffer from a significant systematic error. With this level of precision, our method is highly competitive with traditional methods that rely on the assumption of a steady state. In a general sense, this illustrates that time-varying dynamical structures are not solely obstacles to dynamical mass measurements, but they can also be regarded as assets containing useful information.