Cosmic voids offer an extraordinary opportunity to study the effects of massive neutrinos on cosmological scales. Because they are freely streaming, neutrinos can penetrate the interior of voids more easily than cold dark matter or baryons, which makes their relative contribution to the mass budget in voids much higher than elsewhere in the Universe. In simulations it has recently been shown how various characteristics of voids in the matter distribution are affected by neutrinos, such as their abundance, density profiles, dynamics, and clustering properties. However, the tracers used to identify voids in observations (e.g., galaxies or halos) are affected by neutrinos as well, and isolating the unique neutrino signatures inherent to voids becomes more difficult. In this paper we make use of the DEMNUni suite of simulations to investigate the clustering bias of voids in Fourier space as a function of their core density and compensation. We find a clear dependence on the sum of neutrino masses that remains significant even for void statistics extracted from halos. In particular, we observe that the amplitude of the linear void bias increases with neutrino mass for voids defined in dark matter, whereas this trend gets reversed and slightly attenuated when measuring the relative void-halo bias using voids identified in the halo distribution. Finally, we argue how the original behaviour can be restored when considering observations of the total matter distribution (e.g. via weak lensing), and comment on scale-dependent effects in the void bias that may provide additional information on neutrinos in the future.

Cosmic voids, the underdense regions in the universe, are particularly sensitive to diffuse density components such as cosmic neutrinos. This sensitivity is enhanced by the match between void sizes and the free-streaming scale of massive neutrinos. Using the massive neutrino simulations \texttt{MassiveNuS}, we investigate the effect of neutrino mass on dark matter halos as a function of environment. We find that the halo mass function depends strongly on neutrino mass and that this dependence is more pronounced in voids than in high-density environments. An observational program that measured the characteristic mass of the most massive halos in voids should be able to place novel constraints on the sum of the masses of neutrinos $\sum m_\nu$. The neutrino mass effect in the simulations is quite strong: In a 512$^3$ $h^{-3}$ Mpc$^3$ survey, the mean mass of the 1000 most massive halos in the void interiors is $(4.82 \pm 0.11) \times 10^{12} h^{-1}M_{\odot}$ for $\sum m_\nu = 0.6$ eV and $(8.21 \pm 0.13) \times 10^{12} h^{-1}M_{\odot}$ for $\sum m_\nu = 0.1$ eV. Subaru (SuMIRe), Euclid and WFIRST will have both spectroscopic and weak lensing surveys. Covering volumes at least 50 times larger than our simulations, they should be sensitive probes of neutrino mass through void substructure.

The Quijote simulations are a set of 43100 full N-body simulations spanning more than 7000 cosmological models in the $\{\Omega_{\rm m}, \Omega_{\rm b}, h, n_s, \sigma_8, M_\nu, w \}$ hyperplane. At a single redshift the simulations contain more than 8.5 trillions of particles over a combined volume of 43100 $(h^{-1}{\rm Gpc})^3$. Billions of dark matter halos and cosmic voids have been identified in the simulations, whose runs required more than 35 million core hours. The Quijote simulations have been designed for two main purposes: 1) to quantify the information content on cosmological observables, and 2) to provide enough data to train machine learning algorithms. In this paper we describe the simulations and show a few of their applications. We also release the Petabyte of data generated, comprising hundreds of thousands of simulation snapshots at multiple redshifts, halo and void catalogs, together with millions of summary statistics such as power spectra, bispectra, correlation functions, marked power spectra, and estimated probability density functions.

Do void statistics contain information beyond the tracer 2-point correlation function? Yes! As we vary the sum of the neutrino masses, we find void statistics contain information absent when using just tracer 2-point statistics. Massive neutrinos uniquely affect cosmic voids. We explore their impact on void clustering using both the DEMNUni and MassiveNuS simulations. For voids, neutrino effects depend on the observed void tracers. As the neutrino mass increases, the number of small voids traced by cold dark matter particles increases and the number of large voids decreases. Surprisingly, when massive, highly biased, halos are used as tracers, we find the opposite effect. The scale at which voids cluster, as well as the void correlation, is similarly sensitive to the sum of neutrino masses and the tracers. This scale dependent trend is not due to simulation volume or halo density. The interplay of these signatures in the void abundance and clustering leaves a distinct fingerprint that could be detected with observations and potentially help break degeneracies between different cosmological parameters. This paper paves the way to exploit cosmic voids in future surveys to constrain the mass of neutrinos.

Of the many proposed extensions to the $\Lambda$CDM paradigm, a model in which neutrinos self-interact until close to the epoch of matter-radiation equality has been shown to provide a good fit to current cosmic microwave background (CMB) data, while at the same time alleviating tensions with late-time measurements of the expansion rate and matter fluctuation amplitude. Interestingly, CMB fits to this model either pick out a specific large value of the neutrino interaction strength, or are consistent with the extremely weak neutrino interaction found in $\Lambda$CDM, resulting in a bimodal posterior distribution for the neutrino self-interaction cross section. In this paper, we explore why current cosmological data select this particular large neutrino self-interaction strength, and by consequence, disfavor intermediate values of the self-interaction cross section. We show how it is the $\ell \gtrsim 1000$ CMB temperature anisotropies, most recently measured by the Planck satellite, that produce this bimodality. We also establish that smaller scale temperature data, and improved polarization data measuring the temperature-polarization cross-correlation, will best constrain the neutrino self-interaction strength. We forecast that the upcoming Simons Observatory should be capable of distinguishing between the models.

New physics in the neutrino sector might be necessary to address anomalies between different neutrino oscillation experiments. Intriguingly, it also offers a possible solution to the discrepant cosmological measurements of $H_0$ and $\sigma_8$. We show here that delaying the onset of neutrino free-streaming until close to the epoch of matter-radiation equality can naturally accommodate a larger value for the Hubble constant $H_0=72.3 \pm 1.4$ km/s/Mpc and a lower value of the matter fluctuations $\sigma_8=0.786\pm 0.020$, while not degrading the fit to the cosmic microwave background (CMB) damping tail. We achieve this by introducing neutrino self-interactions in the presence of a non-vanishing sum of neutrino masses. This strongly interacting neutrino cosmology prefers $N_{\rm eff} = 4.02 \pm 0.29$, which has interesting implications for particle model-building and neutrino oscillation anomalies. We show that the absence of the neutrino free-streaming phase shift on the CMB can be compensated by shifting the value of other cosmological parameters, hence providing an important caveat to the detections made in the literature. Due to their impact on the evolution of the gravitational potential at early times, self-interacting neutrinos and their subsequent decoupling leave a rich structure on the matter power spectrum. In particular, we point out the existence of a novel localized feature appearing on scales entering the horizon at the onset of neutrino free-streaming. While the interacting neutrino cosmology provides a better global fit to current cosmological data, we find that traditional Bayesian analyses penalize the model as compared to the standard cosmological. Our analysis shows that it is possible to find radically different cosmological models that nonetheless provide excellent fits to the data, hence providing an impetus to thoroughly explore alternate cosmological scenarios.

Cosmic voids, the less dense patches of the Universe, are promising laboratories to extract cosmological information. Thanks to their unique low density character, voids are extremely sensitive to diffuse components such as neutrinos and dark energy, and represent ideal environments to study modifications of gravity, where the effects of such modifications are expected to be more prominent. Robust void-related observables, including for example redshift-space distortions (RSD) and weak lensing around voids, are a promising way to chase and test new physics. Cosmological analysis of the large-scale structure of the Universe predominantly relies on the high density regions. Current and upcoming surveys are designed to optimize the extraction of cosmological information from these zones, but leave voids under-exploited. A dense, large area spectroscopic survey with imaging capabilities is ideal to exploit the power of voids fully. Besides helping illuminate the nature of dark energy, modified gravity, and neutrinos, this survey will give access to a detailed map of under-dense regions, providing an unprecedented opportunity to observe and study a so far under-explored galaxy population.

The discovery of the electromagnetic counterpart to GW170817 severely constrains the tensor mode propagation speed, eliminating a large model space of Horndeski theory. We use the cosmic microwave background data from Planck and the joint analysis of the BICEP2/Keck Array and Planck, galaxy clustering data from the SDSS LRG survey, BOSS baryon acoustic oscillation data, and redshift space distortion measurements to place constraints on the remaining Horndeski parameters. We evolve the Horndeski parameters as power laws with both the amplitude and power law index free. We find a 95% CL upper bound on the present-day coefficient of the Hubble friction term in the cosmological propagation of gravitational waves is 2.38, whereas General Relativity gives 2 at all times. While an enhanced friction suppresses the amplitude of the reionization bump of the primordial B-mode power spectrum at $\ell < 10$, our result limits the suppression to be less than 0.8%. This constraint is primarily due to the scalar integrated Sachs-Wolfe effect in temperature fluctuations at low multipoles.

Some authors have proposed that electron energy distributions in H II regions and planetary nebulae may be significantly nonthermal, and kappa-distributions have been suggested as being appropriate. Here it is demonstrated that the electron energy distribution function is extremely close to a Maxwellian up to electron kinetic energies ~13 eV in HII regions, and up to ~16eV in planetary nebulae: kappa-distributions are inappropriate. The small departures from a Maxwellian have negligible effects on line ratios. When observed line ratios in H II regions deviate from models with a single electron temperature, it must arise from spatial variations in electron temperature, rather than local deviations from a Maxwellian.

We present an analysis of the complex gas hydrodynamics in the X-ray luminous galaxy cluster RXJ1347.5-1145 caught in the act of merging with a subcluster to its southeast using a combined $186$ ks Chandra exposure, $2.5$ times greater than previous analyses. The primary cluster hosts a sloshing cold front spiral traced by four surface brightness edges $5.^{\prime \prime}85^{+0.04}_{-0.03}$ west, $7.^{\prime \prime}10^{+0.07}_{-0.03}$ southeast, $11.^{\prime \prime}5^{+1.3}_{-1.2}$ east, and $16.^{\prime \prime}7^{+0.3}_{-0.5}$ northeast from the primary central dominant galaxy, suggesting the merger is in the plane of the sky. We measure temperature and density ratios across these edges, confirming they are sloshing cold fronts. We observe the eastern edge of the subcluster infall shock, confirming the observed subcluster is traveling from the southwest to the northeast in a clockwise orbit. We measure a shock density contrast of $1.38^{+0.16}_{-0.15}$ and infer a Mach number $1.25\pm0.08$ and a shock velocity of $2810^{+210}_{-240}$ km s$^{-1}$. Temperature and entropy maps show cool, low entropy gas trailing the subcluster in a southwestern tail, consistent with core shredding. Simulations suggest a perturber in the plane of the sky on a clockwise orbit would produce a sloshing spiral winding counterclockwise, opposite to that observed. The most compelling solution to this discrepancy is that the observed southeastern subcluster is on its first passage, shock heating gas during its clockwise infall, while the main cluster's clockwise cold front spiral formed from earlier encounters with a second perturber orbiting counterclockwise.