Abstract: We revisit the joint constraints in the mixed hot dark matter scenario in which both thermally produced QCD axions and relic neutrinos are present. Upon recomputing the cosmological axion abundance via recent advances in the literature, we improve the state-of-the-art analyses and provide updated bounds on axion and neutrino masses. By avoiding approximate methods, such as the instantaneous decoupling approximation, and limitations due to the limited validity of the perturbative approach in QCD that forced to artificially divide the constraints from the axion-pion and the axion-gluon production channels, we find robust and self-consistent limits. We investigate the two most popular axion frameworks: KSVZ and DFSZ. From Big Bang Nucleosynthesis (BBN) light element abundances data we find for the KSVZ axion Delta N-eff < 0.31 and an axion mass bound m(a) < 0.53 eV (i.e., a bound on the axion decay constant f(a) > 1.07 x 10(7) GeV) both at 95% CL. These BBN bounds are improved to Delta N-eff < 0.14 and m(a) < 0.16 eV (f(a) > 3.56 x 10(7) GeV) if a prior on the baryon energy density from Cosmic Microwave Background (CMB) data is assumed. When instead considering cosmological observations from the CMB temperature, polarization and lensing from the Planck satellite combined with large scale structure data we find Delta N-eff < 0.23, m(a) < 0.28 eV (f(a) > 2.02 x 10(7) GeV) and Sigma m(nu) < 0.16 eV at 95% CL. This corresponds approximately to a factor of 5 improvement in the axion mass bound with respect to the existing limits. Very similar results are obtained for the DFSZ axion. We also forecast upcoming observations from future CMB and galaxy surveys, showing that they could reach percent level errors for m(a) similar to 1 eV.

Abstract: We study the impact of the cosmological parameters uncertainties on the measurements of primordial non-Gaussianity through the large-scale non-Gaussian halo bias effect. While this is not expected to be an issue for the standard Lambda CDM model, it may not be the case for more general models that modify the large-scale shape of the power spectrum. We consider the so-called local non-Gaussianity model, parametrized by the f(NL) non-Gaussianity parameter which is zero for a Gaussian case, and make forecasts on f(NL) from planned surveys, alone and combined with a Planck CMB prior. In particular, we consider EUCLID- and LSST-like surveys and forecast the correlations among f(NL) and the running of the spectral index alpha(s), the dark energy equation of state w, the effective sound speed of dark energy perturbations c(s)(2), the total mass of massive neutrinos M-nu = Sigma m(nu), and the number of extra relativistic degrees of freedom N-nu(rel). Neglecting CMB information on f(NL) and scales k > 0.03h/Mpc, we find that, if N-nu(rel) is assumed to be known, the uncertainty on cosmological parameters increases the error on f(NL) by 10 to 30% depending on the survey. Thus the f(NL) constraint is remarkable robust to cosmological model uncertainties. On the other hand, if N-nu(rel) is simultaneously constrained from the data, the f(NL) error increases by similar to 80%. Finally, future surveys which provide a large sample of galaxies or galaxy clusters over a volume comparable to the Hubble volume can measure primordial non-Gaussianity of the local form with a marginalized 1-sigma error of the order Delta f(NL) similar to 2 – 5, after combination with CMB priors for the remaining cosmological parameters. These results are competitive with CMB bispectrum constraints achievable with an ideal CMB experiment.

Abstract: Collapsing supermassive stars (M greater than or similar to 3 x 10(4) M-circle dot) at high redshifts can naturally provide seeds and explain the origin of the supermassive black holes observed in the centers of nearly all galaxies. During the collapse of supermassive stars, a burst of non-thermal neutrinos is generated with a luminosity that could greatly exceed that of a conventional core collapse supernova explosion. In this work, we investigate the extent to which the neutrinos produced in these explosions can be observed via coherent elastic neutrino-nucleus scattering (CEvNS). Large scale direct dark matter detection experiments provide particularly favorable targets. We find that upcoming O(100) tonne-scale experiments will be sensitive to the collapse of individual supermassive stars at distances as large as O(10) Mpc.

Abstract: Accurate power spectrum (or correlation function) covariance matrices are a crucial requirement for cosmological parameter estimation from large scale structure surveys. In order to minimize reliance on computationally expensive mock catalogs, it is important to have a solid analytic understanding of the different components that make up a covariance matrix. Considering the matter power spectrum covariance matrix, it has recently been found that there is a potentially dominant effect on mildly non-linear scales due to power in modes of size equal to and larger than the survey volume. This beat coupling effect has been derived analytically in perturbation theory and while it has been tested with simulations, some questions remain unanswered. Moreover, there is an additional effect of these large modes, which has so far not been included in analytic studies, namely the effect on the estimated average density which enters the power spectrum estimate. In this article, we work out analytic, perturbation theory based expressions including both the beat coupling and this local average effect and we show that while, when isolated, beat coupling indeed causes large excess covariance in agreement with the literature, in a realistic scenario this is compensated almost entirely by the local average effect, leaving only similar to 10% of the excess. We test our analytic expressions by comparison to a suite of large N-body simulations, using both full simulation boxes and subboxes thereof to study cases without beat coupling, with beat coupling and with both beat coupling and the local average effect. For the variances, we find excellent agreement with the analytic expressions for k < 0.2 hMpc(-1) at z = 0.5, while the correlation coefficients agree to beyond k = 0.4 hMpc(-1). As expected, the range of agreement increases towards higher redshift and decreases slightly towards z = 0. We finish by including the large-mode effects in a full covariance matrix description for arbitrary survey geometry and confirming its validity using simulations. This may be useful as a stepping stone towards building an actual galaxy (or other tracer's) power spectrum covariance matrix.

Abstract: We consider whether the non-Gaussian scale-dependent halo bias can be used not only to constrain the local form of non-Gaussianity but also to distinguish among different shapes. In particular, we ask whether it can constrain the behavior of the primordial three-point function in the squeezed limit where one of the momenta is much smaller than the other two. This is potentially interesting since the observation of a three-point function with a squeezed limit that does not go like the local nor equilateral templates would be a signal of non-trivial dynamics during inflation. To this end we use the quasi-single field inflation model of Chen & Wang [1, 2] as a representative two-parameter model, where one parameter governs the amplitude of non-Gaussianity and the other the shape. We also perform a model-independent analysis by parametrizing the scale-dependent bias as a power-law on large scales, where the power is to be constrained from observations. We find that proposed large-scale structure surveys (with characteristics similar to the dark energy task force stage IV surveys) have the potential to distinguish among the squeezed limit behavior of different bispectrum shapes for a wide range of fiducial model parameters. Thus the halo bias can help discriminate between different models of inflation.