Many works have concentrated on the observable signatures of the dark matter being an ultralight axion-like particle (ALP). We concentrate on a particularly dramatic signature in the late-time cosmological matter power spectrum that occurs if the symmetry breaking that establishes the ALP happens after inflation -- white-noise density fluctuations that dominate at small scales over the adiabatic fluctuations from inflation. These fluctuations alter the early history of nonlinear structure formation. We find that for symmetry breaking scales of $f_A \sim 10^{13}-10^{15}$GeV, which requires a high effective maximum temperature after inflation, ALP dark matter with particle mass of $m_A \sim 10^{-13}-10^{-20}$eV could significantly change the number of high-redshift dwarf galaxies, the reionization history, and the Ly$\alpha$ forest. We consider all three observables. We find that the Ly$\alpha$ forest is the most constraining of current observables, excluding $f_A \gtrsim 10^{15}$GeV ($m_A \lesssim 10^{-17}$eV) in the simplest model for the ALP and considerably lower values in models coupled to a hidden strongly interacting sector ($f_A \gtrsim 10^{13}$GeV and $m_A \lesssim 10^{-13}$eV). Observations that constrain the extremely high-redshift tail of reionization may disfavor similar levels of isocurvature fluctuations as the forest. Future $z\sim 20-30$ 21cm observations have the potential to improve these constraints further using that the supersonic motions of the isocurvature-enhanced abundance of $\sim10^4M_\odot$ halos would shock heat the baryons, sourcing large BAO features.

We explore a potential new probe of reionization: the cross-correlation of high-redshift galaxies with maps of the thermal Sunyaev-Zel'dovich (tSZ) effect. We consider two types of high redshift galaxies: Lyman break galaxies (LBGs) and Lyman-$\alpha$ emitters (LAEs). LBGs and LAEs will be detected in large numbers at high redshift ($z \approx$ 4 to 7) by ongoing and future surveys. We consider a LBG sample like that expected from the The Rubin Observatory Legacy Survey of Space and Time (LSST), and a selection of LAEs modelled after the Subaru SILVERRUSH program, but covering a much larger sky fraction. The tSZ effect is sensitive to a line-of-sight integral of the ionized gas pressure, and can be measured across large patches of the sky using multi-frequency CMB surveys. We consider forecast tSZ maps from CMB Stage 4 and more futuristic observations. Using a suite of hydrodynamical simulations, we show that the high-redshift galaxies are correlated with the tSZ signal from reionization, with cross-power amplitude of order $10^{-15}$ at $\ell \sim 1000$. The cross-spectra between LBGs/LAEs with tSZ maps contain information about the reionization history of the Universe, such as the distribution of bubble sizes, and could in principle be used to directly measure the timing of reionization. The amplitude of the signal is small, however, and its detectability is hindered by low-redshift contributions to tSZ maps and by instrumental noise. If the the low-redshift contribution to the observed tSZ signal is suppressed by masking of massive halos, a combination of overlapping futuristic CMB and galaxy surveys could potentially probe this signal.

We examine the impact of the baryon-dark matter relative velocity on intergalactic small-scale structure and the 21 cm signal during the epoch of reionization (EoR). Streaming velocities reduced the clumping of the intergalactic medium (IGM) on mass scales of $\sim 10^4 - 10^8$ M$_{\odot}$. This effect produced a distinct baryon acoustic oscillation (BAO) feature in the EoR 21 cm power spectrum at wave numbers $k\sim 0.1$ h/Mpc, near which forthcoming surveys will be most sensitive. In contrast to the highly uncertain impact of streaming velocities on star formation, the effect on clumping can be calculated more precisely because it is set mainly by cosmology and straightforward gas dynamics. We quantify the latter using fully coupled radiation-hydrodynamics simulations that capture the Jeans scale of the pre-EoR gas. The clumping factor of ionized gas is reduced by 5-10\% in regions that had typical streaming velocities of 30 km/s at recombination. The suppression peaks $\approx 5$ Myr after a region is reionized, but washes out within 200 Myr due to pressure smoothing of the gas. Using these results, we model the corresponding impact on the EoR 21 cm power spectrum and find that the BAO feature is most likely to appear early in reionization ($\approx$ 10 \% ionization).During this phase, the signal may appear at the 1 \% (5 \%) level at $k \sim 0.1 (0.06)$ h/\Mpc with an amplitude that varies by less than a factor of 10 across a range of reionization histories. We also provide a simple model for the signal originating from streaming velocity's impact on ionizing sources, which can vary by 4 orders of magnitude depending on highly uncertain source properties. We find that the clumping-driven signal is likely to dominate the source-driven one unless population III star formation in halos of masses $10^6 - 10^8$ M$_{\odot}$ was efficient enough to drive the first 10\% of reionization.

We compare a sample of five high-resolution, high S/N Ly$\alpha$ forest spectra of bright $6<z \lesssim 6.5$ QSOs aimed at spectrally resolving the last remaining transmission spikes at $z>5$ with those obtained from mock absorption spectra from the Sherwood and Sherwood-Relics suites of hydrodynamical simulations of the intergalactic medium (IGM). We use a profile fitting procedure for the inverted transmitted flux, $1-F$, similar to the widely used Voigt profile fitting of the transmitted flux $F$ at lower redshifts, to characterise the transmission spikes that probe predominately underdense regions of the IGM. We are able to reproduce the width and height distributions of the transmission spikes, both with optically thin simulations of the post-reionization Universe using a homogeneous UV background and full radiative transfer simulations of a late reionization model. We find that the width of the fitted components of the simulated transmission spikes is very sensitive to the instantaneous temperature of the reionized IGM. The internal structures of the spikes are more prominant in low temeperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution ($\leq $ 8 km/s) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of $T_0= (11000 \pm 1600,10500\pm 2100,12000 \pm 2200)$ K at $z= (5.4,5.6,5.8)$. This is weakly dependent on the slope of the temperature-density relation, which is favoured to be moderately steeper than isothermal. In the inhomogeneous, late reionization, full radiative transfer simulations where islands of neutral hydrogen persist to $z\sim5.3$, the width distribution of the observed transmission spikes is consistent with the range of $T_0$ caused by spatial fluctuations in the temperature-density relation.

In this paper we present an efficient method for constraining both thermal and non-thermal Dark Matter (DM) scenarios with the Lyman-$\alpha$ forest, based on a simple and flexible parametrisation capable to reproduce the small scale clustering signal of a large set of non-cold DM (nCDM) models. We extract new limits on the fundamental DM properties, through an extensive analysis of the high resolution, high redshift data obtained by the MIKE/HIRES spectrographs. By using a large suite of hydrodynamical simulations, we determine constraints on both astrophysical, cosmological, and nCDM parameters by performing a full Monte Carlo Markov Chain (MCMC) analysis. We obtain a marginalised upper limit on the largest possible scale at which a power suppression induced by nearly any nCDM scenario can occur, i.e. $\alpha<0.03~{\rm{Mpc}}/h$ (2$\sigma$ C.L.). We explicitly describe how to test several of the most viable nCDM scenarios without the need to run any specific numerical simulations, due to the novel parametrisation proposed, and due to a new scheme that interpolates between the cosmological models explored. The shape of the linear matter power spectrum for standard thermal warm DM models appear to be in mild tension ($\sim 2\sigma$ C.L.) with the data, compared to non-thermal scenarios. We show that a DM fluid composed by both a warm (thermal) and a cold component is also in tension with the Lyman-$\alpha$ forest, at least for large $\alpha$ values. This is the first study that allows to probe the linear small scale shape of the DM power spectrum for a large set of nCDM models.

Fuzzy Dark Matter (FDM) represents an alternative and intriguing description of the standard Cold Dark Matter (CDM) fluid, able to explain the lack of direct detection of dark matter particles in the GeV sector and to alleviate small scales tensions in the cosmic large-scale structure formation. Cosmological simulations of FDM models in the literature were performed either with very expensive high-resolution grid-based simulations of individual haloes or through N-body simulations encompassing larger cosmic volumes but resorting on significant approximations in the FDM non-linear dynamics to reduce their computational cost. With the use of the new N-body cosmological hydrodynamical code AX-GADGET, we are now able not only to overcome such numerical problems, but also to combine a fully consistent treatment of FDM dynamics with the presence of gas particles and baryonic physical processes, in order to quantify the FDM impact on specific astrophysical observables. In particular, in this paper we perform and analyse several hydrodynamical simulations in order to constrain the FDM mass by quantifying the impact of FDM on Lyman-$\alpha$ forest observations, as obtained for the first time in the literature in a N-body setup without approximating the FDM dynamics. We also study the statistical properties of haloes, exploiting the large available sample, to extract information on how FDM affects the abundance, the shape, and density profiles of dark matter haloes.

We present new measurements of the free-streaming of warm dark matter (WDM) from Lyman-$\alpha$ flux-power spectra. We use data from the medium resolution, intermediate redshift XQ-100 sample observed with the X-shooter spectrograph ($z=3 - 4.2$) and the high-resolution, high-redshift sample used in Viel et al. (2013) obtained with the HIRES/MIKE spectrographs ($z=4.2 - 5.4$). Based on further improved modelling of the dependence of the Lyman-$\alpha$ flux-power spectrum on the free-streaming of dark matter, cosmological parameters, as well as the thermal history of the intergalactic medium (IGM) with hydrodynamical simulations, we obtain the following limits, expressed as the equivalent mass of thermal relic WDM particles. The XQ-100 flux power spectrum alone gives a lower limit of 1.4 keV, the re-analysis of the HIRES/MIKE sample gives 4.1 keV while the combined analysis gives our best and significantly strengthened lower limit of 5.3 keV (all 2$\sigma$ C.L.). The further improvement in the joint analysis is partly due to the fact that the two data sets have different degeneracies between astrophysical and cosmological parameters that are broken when the data sets are combined, and more importantly on chosen priors on the thermal evolution. These results all assume that the temperature evolution of the IGM can be modelled as a power law in redshift. Allowing for a non-smooth evolution of the temperature of the IGM with sudden temperature changes of up to 5000K reduces the lower limit for the combined analysis to 3.5 keV. A WDM with smaller thermal relic masses would require, however, a sudden temperature jump of $5000\,K$ or more in the narrow redshift interval $z=4.6-4.8$, in disagreement with observations of the thermal history based on high-resolution resolution Lyman-$\alpha$ forest data and expectations for photo-heating and cooling in the low density IGM at these redshifts.

We present a semi-analytic model for the Lyman-$\alpha$ forest that is inspired by the Halo Model. This model is built on the absorption line decomposition of the forest. Flux correlations are decomposed into those within each absorption line (the 1-absorber term) and those between separate lines (the 2-absorber term), treating the lines as biased tracers of the underlying matter fluctuations. While the nonlinear exponential mapping between optical depth and flux requires an infinite series of moments to calculate any statistic, we show that this series can be resumed (truncating at the desired order in the linear matter overdensity). We focus on the $z=2-3$ line-of-sight power spectrum. Our model finds that 1-absorber term dominates the power on all scales, with most of its contribution coming from HI columns of $10^{14}-10^{15}\;\mathrm{cm^{-2}}$, while the smaller 2-absorber contribution comes from lower columns that trace overdensities of a few. The prominence of the 1-absorber correlations indicates that the line-of-sight power spectrum is shaped principally by the lines' number densities and their absorption profiles, with correlations between lines contributing to a lesser extent. We present intuitive formulae for the effective optical depth as well as the large-scale limits of 1-absorber and 2-absorber terms, which simplify to integrals over the HI column density distribution with different equivalent-width weightings. With minimalist models for the bias of absorption systems and their peculiar velocity broadening, our model predicts values for the density bias and velocity gradient bias that are consistent with those found in simulations.

We investigate constraints on scalar dark matter (DM) by analyzing the Lyman-alpha forest, which probes structure formation at medium and small scales, and also by studying its cosmological consequences at high and low redshift. For scalar DM that constitutes more than 30% of the total DM density, we obtain a lower limit m >~ 10^{-21} eV for the mass of scalar DM. This implies an upper limit on the initial field displacement (or the decay constant for an axion-like field) of phi <~ 10^{16} GeV. We also derive limits on the energy scale of cosmic inflation and establish an upper bound on the tensor-to-scalar ratio of r < 10^{-3} in the presence of scalar DM. Furthermore, we show that there is very little room for ultralight scalar DM to solve the "small-scale crisis" of cold DM without spoiling the Lyman-alpha forest results. The constraints presented in this paper can be used for testing generic theories that contain light scalar fields.

We present constraints on the masses of extremely light bosons dubbed fuzzy dark matter from Lyman-$\alpha$ forest data. Extremely light bosons with a De Broglie wavelength of $\sim 1$ kpc have been suggested as dark matter candidates that may resolve some of the current small scale problems of the cold dark matter model. For the first time we use hydrodynamical simulations to model the Lyman-$\alpha$ flux power spectrum in these models and compare with the observed flux power spectrum from two different data sets: the XQ-100 and HIRES/MIKE quasar spectra samples. After marginalization over nuisance and physical parameters and with conservative assumptions for the thermal history of the IGM that allow for jumps in the temperature of up to $5000\rm\,K$, XQ-100 provides a lower limit of 7.1$\times 10^{-22}$ eV, HIRES/MIKE returns a stronger limit of 14.3$\times 10^{-22}$ eV, while the combination of both data sets results in a limit of 20 $\times 10^{-22}$ eV (2$\sigma$ C.L.). The limits for the analysis of the combined data sets increases to 37.5$\times 10^{-22}$ eV (2$\sigma$ C.L.) when a smoother thermal history is assumed where the temperature of the IGM evolves as a power-law in redshift. Light boson masses in the range $1-10 \times10^{-22}$ eV are ruled out at high significance by our analysis, casting strong doubts that FDM helps solve the "small scale crisis" of the cold dark matter models.

We present the Lyman-$\alpha$ flux power spectrum measurements of the XQ-100 sample of quasar spectra obtained in the context of the European Southern Observatory Large Programme "Quasars and their absorption lines: a legacy survey of the high redshift universe with VLT/XSHOOTER". Using $100$ quasar spectra with medium resolution and signal-to-noise ratio we measure the power spectrum over a range of redshifts $z = 3 - 4.2$ and over a range of scales $k = 0.003 - 0.06\,\mathrm{s\,km^{-1}}$. The results agree well with the measurements of the one-dimensional power spectrum found in the literature. The data analysis used in this paper is based on the Fourier transform and has been tested on synthetic data. Systematic and statistical uncertainties of our measurements are estimated, with a total error (statistical and systematic) comparable to the one of the BOSS data in the overlapping range of scales, and smaller by more than $50\%$ for higher redshift bins ($z>3.6$) and small scales ($k > 0.01\,\mathrm{s\,km^{-1}}$). The XQ-100 data set has the unique feature of having signal-to-noise ratios and resolution intermediate between the two data sets that are typically used to perform cosmological studies, i.e. BOSS and high-resolution spectra (e.g. UVES/VLT or HIRES). More importantly, the measured flux power spectra span the high redshift regime which is usually more constraining for structure formation models.

We present the calculation of the Lyman-alpha (Lyman-$\alpha$) transmitted flux fluctuations with full relativistic corrections to the first order. Even though several studies exist on relativistic effects in galaxy clustering, this is the first study to extend the formalism to a different tracer of underlying matter at unique redshift range ($z=2-5$). Furthermore, we show a comprehensive application of our calculations of the Quasar-Lyman-$\alpha$ cross-correlation function. Our results indicate that the signal of relativistic effects is sizeable at Baryonic Acoustic Oscillation (BAO) scale mainly due to the large difference in density bias factors of our tracers. We construct an observable, the anti-symmetric part of the cross-correlation function, that is dominated by the relativistic signal and offers a new way to measure the relativistic terms at relatively small scales. The analysis shows that relativistic effects are important when considering cross-correlations between tracers with very different biases, and should be included in the data analysis of the current and future surveys. Moreover, the idea presented in this paper is highly complementary to other techniques and observable trying to isolate the effect of the relativistic corrections and thus test the validity of the theory of gravity beyond the Newtonian regime.

We present a comprehensive analysis of high resolution hydrodynamic simulations in terms of Lyman-alpha and Lyman-beta one dimensional flux power spectra ($P_{\alpha\alpha}$ and $P_{\beta\beta}$). In particular, we focus on the behaviour that the flux auto-power spectra and cross-power spectra ($P_{\alpha\beta}$) display when the intergalactic medium (IGM) thermal history is changed in a range of values that bracket a reference model, while cosmological parameters are kept fixed to best fit the cosmic microwave background data. We present empirical fits that describe at the sub-percent level the dependence of the power spectra on the thermal parameters. At the largest scales, the power spectra show a constant bias between each other that is set by the parameters describing the IGM thermal state. The cross-power spectrum has an oscillatory pattern and crosses zero at a scale which depends on $T_0$, the IGM temperature at the mean density, for reasonable values of the power-law index $\gamma$ of the IGM temperature-density relation ($T=T_0(1+\delta)^{\gamma-1}$). By performing a Fisher matrix analysis, we find that the power spectrum $P_{\beta\beta}$ is more sensitive to the thermal history than $P_{\alpha\alpha}$ alone, due to the fact that it probes denser regions than Lyman-alpha. When we combine the power and cross spectra the constraints on $\gamma$ can be improved by a factor $\sim 4$, while the constraints on $T_0$ improve by a factor of $\sim 2$. We address the role of signal-to-noise and resolution by mocking realistic observations and we conclude that the framework presented in this work can significantly improve the knowledge of the IGM thermal state, which will in turn guarantee better constraints on IGM-derived cosmological parameters.

When measuring the one-dimensional power spectrum of the Ly$\alpha$ forest, it is common to measure the power spectrum in flux fluctuations red-ward of the Ly$\alpha$ emission of quasars and subtract this power from the measurements of the Ly$\alpha$ flux power spectrum. This removes excess power present in the Ly$\alpha$ forest which is believed to be dominated by metal absorption by the low-redshift metals uncorrelated with the neutral hydrogen aborbing in Ly$\alpha$. In this brief report we note that, assuming the contaminants are additive in optical depth, the correction contains a second order term. We estimate the magnitude of this term for two currently published measurements of the 1D Ly$\alpha$ flux power spectrum and show that it is negligible for the current generation of measurements. However, future measurements will have to take this into account when errorbars improve by a factor of two or more.

The Lyman-$\beta$ forest refers to a region in the spectra of distant quasars that lies between the rest-frame Lyman-$\beta$ and Lyman-$\gamma$ emissions. The forest in this region is dominated by a combination of absorption due to resonant Ly$\alpha$ and Ly$\beta$ scattering. When considering the 1D Ly$\beta$ forest in addition to the 1D Ly$\alpha$ forest, the full statistical description of the data requires four 1D power spectra: Ly$\alpha$ and Ly$\beta$ auto-power spectra and the Ly$\alpha$-Ly$\beta$ real and imaginary cross-power spectra. We describe how these can be measured using an optimal quadratic estimator that naturally disentangles Ly$\alpha$ and Ly$\beta$ contributions. Using a sample of approximately 60,000 quasar sight-lines from the BOSS Data Release 9, we make the measurement of the one-dimensional power spectrum of fluctuations due to the Ly$\beta$ resonant scattering. While we have not corrected our measurements for resolution damping of the power and other systematic effects carefully enough to use them for cosmological constraints, we can robustly conclude the following: i) Ly$\beta$ power spectrum and Ly$\alpha$-Ly$\beta$ cross spectra are detected with high statistical significance; ii) the cross-correlation coefficient is $\approx 1$ on large scales; iii) the Ly$\beta$ measurements are contaminated by the associated OVI absorption, which is analogous to the SiIII contamination of the Ly$\alpha$ forest. Measurements of the Ly$\beta$ forest will allow extension of the usable path-length for the Ly$\alpha$ measurements while allowing a better understanding of the physics of intergalactic medium and thus more robust cosmological constraints.

We use the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 9 (DR9) to detect and measure the position of the Baryonic Acoustic Oscillation (BAO) feature in the three-dimensional correlation function in the Lyman-alpha forest flux fluctuations at a redshift z=2.4. The feature is clearly detected at significance between 3 and 5 sigma (depending on the broadband model and method of error covariance matrix estimation) and is consistent with predictions of the standard LCDM model. We assess the biases in our method, stability of the error covariance matrix and possible systematic effects. We fit the resulting correlation function with several models that decouple the broadband and acoustic scale information. For an isotropic dilation factor, we measure 100x(alpha_iso-1) = -1.6 ^{+2.0+4.3+7.4}_{-2.0-4.1-6.8} (stat.) +/- 1.0 (syst.) (multiple statistical errors denote 1,2 and 3 sigma confidence limits) with respect to the acoustic scale in the fiducial cosmological model (flat LCDM with Omega_m=0.27, h=0.7). When fitting separately for the radial and transversal dilation factors we find marginalised constraints 100x(alpha_par-1) = -1.3 ^{+3.5+7.6 +12.3}_{-3.3-6.7-10.2} (stat.) +/- 2.0 (syst.) and 100x(alpha_perp-1) = -2.2 ^{+7.4+17}_{-7.1-15} +/- 3.0 (syst.). The dilation factor measurements are significantly correlated with cross-correlation coefficient of ~ -0.55. Errors become significantly non-Gaussian for deviations over 3 standard deviations from best fit value. Because of the data cuts and analysis method, these measurements give tighter constraints than a previous BAO analysis of the BOSS DR9 Lyman-alpha forest sample, providing an important consistency test of the standard cosmological model in a new redshift regime.

We study the compression of information present in the correlated perturbations to the luminosity distance in the low-redshift ($z<0.1$) supernovae Ia due to peculiar velocities of these supernovae. We demonstrate that the na\"{i}ve compression into angular velocity power spectrum does not work efficiently, due to thickness of the spherical shell over which the supernovae are measured. Instead, we show that measurements can be compressed into measurements of $f^2P(k)$, where $f$ is the logarithmic rate of growth of linear perturbations and $P(k)$ is their power spectrum. We develop an optimal quadratic estimator and show that it recovers all information for $\Lambda CDM$ models for surveys of $N\sim10,000$ or more supernovae. We explicitly demonstrate robustness with respect to the assumed fiducial model and the number of power spectrum bins. Using mock catalogues of SNe Ia we estimate that future low redshift surveys will be able to probe $\sigma_8$ to 6% accuracy with $10,000$ SNe Ia.