At present, cosmological observations set the most stringent bound on the neutrino mass scale. Within the standard cosmological model ($\Lambda$CDM), the Planck collaboration reports $\sum m_\nu < 0.12\,\text{eV}$ at 95% CL. This bound, taken at face value, excludes many neutrino mass models. However, unstable neutrinos, with lifetimes shorter than the age of the universe $\tau_\nu \lesssim t_U$, represent a particle physics avenue to relax this constraint. Motivated by this fact, we present a taxonomy of neutrino decay modes, categorizing them in terms of particle content and final decay products. Taking into account the relevant phenomenological bounds, our analysis shows that 2-body decaying neutrinos into BSM particles are a promising option to relax cosmological neutrino mass bounds. We then build a simple extension of the type I seesaw scenario by adding one sterile state $\nu_4$ and a Goldstone boson $\phi$, in which $\nu_i \to \nu_4 \, \phi$ decays can loosen the neutrino mass bounds up to $\sum m_\nu \sim 1\,\text{eV}$, without spoiling the light neutrino mass generation mechanism. Remarkably, this is possible for a large range of the right-handed neutrino masses, from the electroweak up to the GUT scale. We successfully implement this idea in the context of minimal neutrino mass models based on a $U(1)_{\mu-\tau}$ flavor symmetry, which are otherwise in tension with the current bound on $\sum m_\nu$.

It is commonly assumed that the energy density of the Universe was dominated by radiation between reheating after inflation and the onset of matter domination 54,000 years later. While the abundance of light elements indicates that the Universe was radiation dominated during Big Bang Nucleosynthesis (BBN), there is scant evidence that the Universe was radiation dominated prior to BBN. It is therefore possible that the cosmological history was more complicated, with deviations from the standard radiation domination during the earliest epochs. Indeed, several interesting proposals regarding various topics such as the generation of dark matter, matter-antimatter asymmetry, gravitational waves, primordial black holes, or microhalos during a nonstandard expansion phase have been recently made. In this paper, we review various possible causes and consequences of deviations from radiation domination in the early Universe - taking place either before or after BBN - and the constraints on them, as they have been discussed in the literature during the recent years.

Precision measurements of the number of effective relativistic neutrino species and the primordial element abundances require accurate theoretical predictions for early Universe observables in the Standard Model and beyond. Given the complexity of accurately modelling the thermal history of the early Universe, in this work, we extend a previous method presented by the author to obtain simple, fast and accurate early Universe thermodynamics. The method is based upon the approximation that all relevant species can be described by thermal equilibrium distribution functions characterized by a temperature and a chemical potential. We apply the method to neutrino decoupling in the Standard Model and find $N_{\rm eff}^{\rm SM} = 3.045$ -- a result in excellent agreement with previous state-of-the-art calculations. We apply the method to study the thermal history of the Universe in the presence of a very light ($1\,\text{eV}<m_\phi < 1\,\text{MeV}$) and weakly coupled ($\lambda \lesssim 10^{-9}$) neutrinophilic scalar. We find our results to be in excellent agreement with the solution to the exact Liouville equation. Finally, we release a code: NUDEC_BSM (available in both Mathematica and Python formats), with which neutrino decoupling can be accurately and efficiently solved in the Standard Model and beyond: https://github.com/MiguelEA/nudec_BSM .

Local measurements of the Hubble constant currently disagree with the high-precision value that is inferred from the CMB under the assumption of a $\Lambda$CDM cosmology. The significance of this tension clearly motivates studying extensions of the standard cosmological model capable of addressing this outstanding issue. Broadly speaking, models that have been successful in reducing the the tension between the CMB and local measurements (without introducing additional tension in other datasets) require an additional component of the energy density in the Universe at a time close to recombination. In this contribution, I will show that the Majoron -- a pseudo-Goldstone boson arising from the spontaneous breaking of a global lepton number symmetry and often associated with the neutrino mass mechanism -- can help to reduce the Hubble tension. Importantly, I will also show that current CMB observations can constrain neutrino-Majoron couplings as small as $10^{-13}$, which within the type-I seesaw mechanism correspond to scales of lepton number breaking as high as $\sim 1\,\text{TeV}$.

The majoron, a pseudo-Goldstone boson arising from the spontaneous breaking of global lepton number, is a generic feature of many models intended to explain the origin of the small neutrino masses. In this work, we investigate potential imprints in the Cosmic Microwave Background (CMB) arising from massive majorons, should they thermalize with neutrinos after Big Bang Nucleosynthesis via inverse neutrino decays. We show that Planck2018 measurements of the CMB are currently sensitive to neutrino-majoron couplings as small as $\lambda \sim 10^{-13}$, which if interpreted in the context of the type-I seesaw mechanism correspond to a lepton number symmetry breaking scale $v_L \sim \mathcal{O}(100) \, {\rm GeV}$. Additionally, we identify parameter space for which the majoron-neutrino interactions, collectively with an extra contribution to the effective number of relativistic species $N_{\rm eff}$, can ameliorate the outstanding $H_0$ tension.

Big Bang Nucleosynthesis (BBN) is very sensitive to the cosmological expansion rate. If the gravitational constant $G$ took a different value during the nucleosynthesis epoch than today, the primordial abundances of light elements would be affected. In this work, we improve the bounds on this variation using recent determinations of the primordial element abundances, updated nuclear and weak reaction rates and observations of the Cosmic Microwave Background (CMB). When combining the measured abundances and the baryon density from CMB observations by Planck, we find $G_\mathrm{BBN}/G_0 = 0.99^{+0.06}_{-0.05}$ at $2\sigma$ confidence level. If the variation of $G$ is linear in time, we find $\dot{G}/G_0 = 0.7^{+3.8}_{-4.3}\times 10^{-12} \, \mathrm{yr}^{-1}$, again at $2\sigma$. These bounds are significantly stronger than those from previous primordial nucleosynthesis studies, and are comparable and complementary to CMB, stellar, solar system, lunar laser ranging, pulsar timing and gravitational wave constraints.

Invisible neutrino decay modes are difficult to target at laboratory experiments, and current bounds on such decays from solar neutrino and neutrino oscillation experiments are somewhat weak. It has been known for some time that Cosmology can serve as a powerful probe of invisible neutrino decays. In this work, we show that in order for Big Bang Nucleosynthesis to be successful, the invisible neutrino decay lifetime is bounded to be $\tau_\nu > 10^{-3}\,\text{s}$ at 95\% CL. We revisit Cosmic Microwave Background constraints on invisible neutrino decays, and by using Planck2018 observations we find the following bound on the neutrino lifetime: $\tau_\nu > (1.3-0.3)\times 10^{9}\,\text{s} \, \left({m_\nu}/{ 0.05\,\text{eV} }\right)^3$ at $95\%$ CL. We show that this bound is robust to modifications of the cosmological model, in particular that it is independent of the presence of dark radiation. We find that lifetimes relevant for Supernova observations ($\tau_\nu \sim 10^{5}\,\text{s}\, \left({m_\nu}/{ 0.05\,\text{eV} }\right)^3$) are disfavoured at more than $5\,\sigma$ with respect to $\Lambda$CDM given the latest Planck CMB observations. Finally, we show that when including high-$\ell$ Planck polarization data, neutrino lifetimes $\tau_\nu = (2-16)\times 10^{9}\,\text{s} \, \left({m_\nu}/{ 0.05\,\text{eV} }\right)^3$ are mildly preferred -- with a 1-2 $\sigma$ significance -- over neutrinos being stable.

New light states thermally coupled to the Standard Model plasma alter the expansion history of the Universe and impact the synthesis of the primordial light elements. In this work, we carry out an exhaustive and precise analysis of the implications of MeV-scale BSM particles in Big Bang Nucleosynthesis (BBN) and for Cosmic Microwave Background (CMB) observations. We find that, BBN observations set a lower bound on the thermal dark matter mass of $m_\chi > 0.4\,\text{MeV}$ at $2\sigma$. This bound is independent of the spin and number of internal degrees of freedom of the particle, of the annihilation being s-wave or p-wave, and of the annihilation final state. Furthermore, we show that current BBN plus CMB observations constrain purely electrophilic and neutrinophilic BSM species to have a mass, $m_\chi > 3.7\,\text{MeV}$ at $2\sigma$. We explore the reach of future BBN measurements and show that upcoming CMB missions should improve the bounds on light BSM thermal states to $m_\chi > (10-15)\,\text{MeV}$. Finally, we demonstrate that very light BSM species thermally coupled to the SM plasma are highly disfavoured by current cosmological observations.

Although weakly interacting massive particles (WIMPs) have long been among the most studied and theoretically attractive classes of candidates for the dark matter of our universe, the lack of their detection in direct detection and collider experiments has begun to dampen enthusiasm for this paradigm. In this study, we set out to appraise the status of the WIMP paradigm, focusing on the case of dark matter candidates that interact with the Standard Model through a new gauge boson. After considering a wide range of $Z'$ mediated dark matter models, we quantitatively evaluate the fraction of the parameter space that has been excluded by existing experiments, and that is projected to fall within the reach of future direct detection experiments. Despite the existence of stringent constraints, we find that a sizable fraction of this parameter space remains viable. More specifically, if the dark matter is a Majorana fermion, we find that an order one fraction of the parameter space is in many cases untested by current experiments. Future direct detection experiments with sensitivity near the irreducible neutrino floor will be able to test a significant fraction of the currently viable parameter space, providing considerable motivation for the next generation of direct detection experiments.

The number of effective relativistic neutrino species represents a fundamental probe of the thermal history of the early Universe, and as such of the Standard Model of Particle Physics. Traditional approaches to the neutrino decoupling are either very technical and computationally expensive, or assume that neutrinos decouple instantaneously. In this work, we aim to fill the gap between these two approaches by modeling the neutrino decoupling in terms of two simple coupled differential equations for the electromagnetic and neutrino sector temperatures, in which all the relevant interactions are taken into account and which allows for a straightforward implementation of BSM species. Upon including finite temperature QED corrections we reach an accuracy on $N_{\rm eff}$ in the SM of $0.01$. We illustrate the usefulness of this approach to the neutrino decoupling by considering, in a model independent manner, the impact of MeV thermal dark matter on $N_{\rm eff}$. We show that Planck rules out electrophilic and neutrinophilic thermal dark matter particles of $m< 3.0\,\text{MeV}$ at 95\% CL regardless of their spin, and of their annihilation being $s$-wave or $p$-wave. We point out that thermal dark matter particles with non-negligible interactions with both electrons and neutrinos are more elusive to CMB observations than purely electrophilic or neutrinophilic ones. In addition, assisted by the accuracy of our approach, we show that CMB Stage-IV experiments will generically test thermal dark matter particles with $m \lesssim 15\,\text{MeV}$. We make publicly available the codes developed for this study at https://github.com/MiguelEA/nudec_BSM .

We search Dark Energy Survey (DES) Year 3 imaging for galaxy-galaxy strong gravitational lenses using convolutional neural networks, extending previous work with new training sets and covering a wider range of redshifts and colors. We train two neural networks using images of simulated lenses, then use them to score postage stamp images of 7.9 million sources from the Dark Energy Survey chosen to have plausible lens colors based on simulations. We examine 1175 of the highest-scored candidates and identify 152 probable or definite lenses. Examining an additional 20,000 images with lower scores, we identify a further 247 probable or definite candidates. After including 86 candidates discovered in earlier searches using neural networks and 26 candidates discovered through visual inspection of blue-near-red objects in the DES catalog, we present a catalog of 511 lens candidates.

In this paper, we explore in detail the cosmological implications of an abelian $L_\mu-L_\tau$ gauge extension of the Standard Model featuring a light and weakly coupled $Z'$. Such a scenario is motivated by the longstanding $\sim \, 4 \sigma$ discrepancy between the measured and predicted values of the muon's anomalous magnetic moment, $(g-2)_\mu$, as well as the tension between late and early time determinations of the Hubble constant. If sufficiently light, the $Z'$ population will decay to neutrinos, increasing the overall energy density of radiation and altering the expansion history of the early universe. We identify two distinct regions of parameter space in this model in which the Hubble tension can be significantly relaxed. The first of these is the previously identified region in which a $\sim \, 10-20$ MeV $Z'$ reaches equilibrium in the early universe and then decays, heating the neutrino population and delaying the process of neutrino decoupling. For a coupling of $g_{\mu-\tau} \simeq (3-8) \times 10^{-4}$, such a particle can also explain the observed $(g-2)_{\mu}$ anomaly. In the second region, the $Z'$ is very light ($m_{Z'} \sim 1\,\text{eV}$ to $\text{MeV}$) and very weakly coupled ($g_{\mu-\tau} \sim 10^{-13}$ to $10^{-9}$). In this case, the $Z'$ population is produced through freeze-in, and decays to neutrinos after neutrino decoupling. Across large regions of parameter space, we predict a contribution to the energy density of radiation that can appreciably relax the reported Hubble tension, $\Delta N_{\rm eff} \simeq 0.2$.

We present a new mechanism of Baryogenesis and dark matter production in which both the dark matter relic abundance and the baryon asymmetry arise from neutral $B$ meson oscillations and subsequent decays. This set-up is testable at hadron colliders and $B$-factories. In the early Universe, decays of a long lived particle produce $B$ mesons and anti-mesons out of thermal equilibrium. These mesons/anti-mesons then undergo CP violating oscillations before quickly decaying into visible and dark sector particles. Dark matter will be charged under Baryon number so that the visible sector baryon asymmetry is produced without violating the total baryon number of the Universe. The produced baryon asymmetry will be directly related to the leptonic charge asymmetry in neutral $B$ decays; an experimental observable. Dark matter is stabilized by an unbroken discrete symmetry, and proton decay is simply evaded by kinematics. We will illustrate this mechanism with a model that is unconstrained by di-nucleon decay, does not require a high reheat temperature, and would have unique experimental signals -- a positive leptonic asymmetry in $B$ meson decays, a new decay of $B$ mesons into a baryon and missing energy, and a new decay of $b$-flavored baryons into mesons and missing energy. These three observables are testable at current and upcoming collider experiments, allowing for a distinct probe of this mechanism.

We explore the constraints and phenomenology of possibly the simplest scenario that could account at the same time for the active neutrino masses and the dark matter in the Universe within a gauged $U(1)_{B-L}$ symmetry, namely right-handed neutrino dark matter. We find that null searches from lepton and hadron colliders require dark matter with a mass below 900 GeV to annihilate through a resonance. Additionally, the very strong constraints from high-energy dilepton searches fully exclude the model for $ 150 \, \text{GeV} < m_{Z'} < 3 \, \text{TeV}$. We further explore the phenomenology in the high mass region (i.e. masses $\gtrsim \mathcal{O}(1) \, \text{TeV}$) and highlight theoretical arguments, related to the appearance of a Landau pole or an instability of the scalar potential, disfavoring large portions of this parameter space. Collectively, these considerations illustrate that a minimal extension of the Standard Model via a local $U(1)_{B-L}$ symmetry with a viable thermal dark matter candidate is difficult to achieve without fine-tuning. We conclude by discussing possible extensions of the model that relieve tension with collider constraints by reducing the gauge coupling required to produce the correct relic abundance.

The elastic scattering between dark matter particles and radiation represents an attractive possibility to solve a number of discrepancies between observations and standard cold dark matter predictions, as the induced collisional damping would imply a suppression of small-scale structures. We consider this scenario and confront it with measurements of the ionization history of the Universe at several redshifts and with recent estimates of the counts of Milky Way satellite galaxies. We derive a conservative upper bound on the dark matter-photon elastic scattering cross section of $\sigma_{\gamma \rm{DM}} < 8 \times 10^{-10} \, \sigma_T \, \left(m_{\rm DM}/{\rm GeV}\right)$ at $95\%$~CL, about one order of magnitude tighter than previous {constraints from satellite number counts}. Due to the strong degeneracies with astrophysical parameters, the bound on the dark matter-photon scattering cross section derived here is driven by the estimate of the number of Milky Way satellite galaxies. Finally, we also argue that future 21~cm probes could help in disentangling among possible non-cold dark matter candidates, such as interacting and warm dark matter scenarios. Let us emphasize that bounds of similar magnitude to the ones obtained here could be also derived for models with dark matter-neutrino interactions and would be as constraining as the tightest limits on such scenarios.

We use present cosmological observations and forecasts of future experiments to illustrate the power of large-scale structure (LSS) surveys in probing dark matter (DM) microphysics and unveiling potential deviations from the standard $\Lambda$CDM scenario. To quantify this statement, we focus on an extension of $\Lambda$CDM with DM-neutrino scattering, which leaves a distinctive imprint on the angular and matter power spectra. After finding that future CMB experiments (such as COrE+) will not significantly improve the constraints set by the Planck satellite, we show that the next generation of galaxy clustering surveys (such as DESI) could play a leading role in constraining alternative cosmologies and even have the potential to make a discovery. Typically we find that DESI would be an order of magnitude more sensitive to DM interactions than Planck, thus probing effects that until now have only been accessible via $N$-body simulations.

Stringent constraints from direct detection experiments and the Large Hadron Collider motivate us to consider models in which the dark matter does not directly couple to the Standard Model, but that instead annihilates into hidden sector particles which ultimately decay through small couplings to the Standard Model. We calculate the gamma-ray emission generated within the context of several such hidden sector models, including those in which the hidden sector couples to the Standard Model through the vector portal (kinetic mixing with Standard Model hypercharge), through the Higgs portal (mixing with the Standard Model Higgs boson), or both. In each case, we identify broad regions of parameter space in which the observed spectrum and intensity of the Galactic Center gamma-ray excess can easily be accommodated, while providing an acceptable thermal relic abundance and remaining consistent with all current constraints. We also point out that cosmic-ray antiproton measurements could potentially discriminate some hidden sector models from more conventional dark matter scenarios.

We analyze a simple extension of the Standard Model (SM) with a dark sector composed of a scalar and a fermion, both singlets under the SM gauge group but charged under a dark sector symmetry group. Sterile neutrinos, which are singlets under both groups, mediate the interactions between the dark sector and the SM particles, and generate masses for the active neutrinos via the seesaw mechanism. We explore the parameter space region where the observed Dark Matter relic abundance is determined by the annihilation into sterile neutrinos, both for fermion and scalar Dark Matter particles. The scalar Dark Matter case provides an interesting alternative to the usual Higgs portal scenario. We also study the constraints from direct Dark Matter searches and the prospects for indirect detection via sterile neutrino decays to leptons, which may be able to rule out Dark Matter masses below and around 100 GeV.

Utilizing an exhaustive set of simplified models, we revisit dark matter scenarios potentially capable of generating the observed Galactic Center gamma-ray excess, updating constraints from the LUX and PandaX-II experiments, as well as from the LHC and other colliders. We identify a variety of pseudoscalar mediated models that remain consistent with all constraints. In contrast, dark matter candidates which annihilate through a spin-1 mediator are ruled out by direct detection constraints unless the mass of the mediator is near an annihilation resonance, or the mediator has a purely vector coupling to the dark matter and a purely axial coupling to Standard Model fermions. All scenarios in which the dark matter annihilates through $t$-channel processes are now ruled out by a combination of the constraints from LUX/PandaX-II and the LHC.

We explore an extended cosmological scenario where the dark matter is an admixture of cold and additional non-cold species. The mass and temperature of the non-cold dark matter particles are extracted from a number of cosmological measurements. Among others, we consider tomographic weak lensing data and Milky Way dwarf satellite galaxy counts. We also study the potential of these scenarios in alleviating the existing tensions between local measurements and Cosmic Microwave Background (CMB) estimates of the $S_8$ parameter, with $S_8=\sigma_8\sqrt{\Omega_m}$, and of the Hubble constant $H_0$. In principle, a sub-dominant, non-cold dark matter particle with a mass $m_X\sim$~keV, could achieve the goals above. However, the preferred ranges for its temperature and its mass are different when extracted from weak lensing observations and from Milky Way dwarf satellite galaxy counts, since these two measurements require suppressions of the matter power spectrum at different scales. Therefore, solving simultaneously the CMB-weak lensing tensions and the small scale crisis in the standard cold dark matter picture via only one non-cold dark matter component seems to be challenging.

In this paper we explore the possibility that the sterile neutrino and Dark Matter sectors in the Universe have a common origin. We study the consequences of this assumption in the simple case of coupling the dark sector to the Standard Model via a global $U(1)_{B-L}$, broken down spontaneously by a dark scalar. This dark scalar provides masses to the dark fermions and communicates with the Higgs via a Higgs portal coupling. We find an interesting interplay between Dark Matter annihilation to dark scalars - the CP-even that mixes with the Higgs and the CP-odd which becomes a Goldstone boson, the Majoron - and heavy neutrinos, as well as collider probes via the coupling to the Higgs. Moreover, Dark Matter annihilation into sterile neutrinos and its subsequent decay to gauge bosons and quarks, charged leptons or neutrinos lead to indirect detection signatures which are close to current bounds on the gamma ray flux from the galactic center and dwarf galaxies.

In recent years, direct detection, indirect detection, and collider experiments have placed increasingly stringent constraints on particle dark matter, exploring much of the parameter space associated with the WIMP paradigm. In this paper, we focus on the subset of WIMP models in which the dark matter annihilates in the early universe through couplings to either the Standard Model $Z$ or the Standard Model Higgs boson. Considering fermionic, scalar, and vector dark matter candidates within a model-independent context, we find that the overwhelming majority of these dark matter candidates are already ruled out by existing experiments. In the case of $Z$ mediated dark matter, the only scenarios that are not currently excluded are those in which the dark matter is a fermion with an axial coupling and with a mass either within a few GeV of the $Z$ resonance ($m_{\rm DM} \simeq m_Z/2$) or greater than 200 GeV, or with a vector coupling and with $m_{\rm DM} > 6$ TeV. Several Higgs mediated scenarios are currently viable if the mass of the dark matter is near the Higgs pole ($m_{\rm DM} \simeq m_H/2$). Otherwise, the only scenarios that are not excluded are those in which the dark matter is a scalar (vector) heavier than 400 GeV (1160 GeV) with a Higgs portal coupling, or a fermion with a pseudoscalar (CP violating) coupling to the Standard Model Higgs boson. With the exception of dark matter with a purely pseudoscalar coupling to the Higgs, it is anticipated that planned direct detection experiments will probe nearly the entire range of models considered in this study.

The value of the tensor-to-scalar ratio $r$ in the region allowed by the latest $Planck$ 2015 measurements can be associated to a large variety of inflationary models. We discuss here the potential of future Cosmic Microwave Background cosmological observations in disentangling among the possible theoretical scenarios allowed by our analyses of current $Planck$ temperature and polarization data. Rather than focusing only on $r$, we focus as well on the running of the primordial power spectrum, $\alpha_s$ and the running of thereof, $\beta_s$. Our Fisher matrix method benefits from a detailed and realistic appraisal of the expected foregrounds. Future cosmological probes, as the COrE mission, may be able to reach an unprecedented accuracy in the extraction of $\beta_s$ and rule out the most favoured inflationary models.