We present the photometric determination of the bright-end (L_Lya>10^43.5 erg/s) of the Lya luminosity function (LF) within four redshifts windows in the interval 2.2<z<3.3. Our work is based on the Javalambre Photometric Local Universe Survey (J-PLUS) first data-release, which provides multiple narrow-band measurements over ~1000 deg^2, with limiting magnitude r~22. The analysis of high-z Lya-emitting sources over such a wide area is unprecedented, and allows to select a total of ~14,500 hyper-bright (L_Lya>10^43.3 erg/s) Lya-emitting candidates. We test our selection with two spectroscopic follow-up programs at the GTC telescope, confirming ~89% of the targets as line-emitting sources, with ~64% being genuine z~2.2 QSOs. We extend the 2.2<z<3.3 Lya LF for the first time above L_Lya~10^44 erg/s and down to densities of ~10^-8 Mpc^-3. Our results unveil with high detail the Schechter exponential-decay of the brightest-end of the Lya LF, complementing the power-law component of previous LF determinations at 43.3<Log_10(L_Lya / [erg/s])<44. We measure Phi^*=(3.33+-0.19)x10^-6, Log(L^*)=44.65+-0.65 and alpha=-1.35+-0.84 as an average over the redshifts we probe. These values are significantly different than the typical Schechter parameters measured for the Lya LF of high-z star-forming LAEs. This suggests that z>2 AGN/QSOs (likely dominant in our samples) are described by a structurally different LF than z>2 star-forming LAEs, namely with L^*_QSOs ~ 100 L^*_LAEs and Phi^*_QSOs ~ 10^-3 Phi^*_LAEs. Finally, our method identifies very efficiently as high-z line-emitters sources without previous spectroscopic confirmation, currently classified as stars (~2000 objects in each redshift bin, on average). Assuming a large predominance of Lya-emitting AGN/QSOs in our samples, this supports the scenario by which these are the most abundant class of z>2 Lya emitters at L_Lya>10^43.3 erg/s.

We study the highly magnified arc SGAS J122651.3+215220 caused by a star-forming galaxy at $z_s=2.93$ crossing the lensing caustic cast by the galaxy cluster SDSS J1226+2152 ($z_l=0.43$), using Hubble Space Telescope observations. We report in the arc several asymmetric surface brightness features whose angular separations are a fraction of an arcsecond from the lensing critical curve and appear to be highly but unequally magnified image pairs of underlying compact sources, with one brightest pair having clear asymmetry consistently across four filters. One explanation of unequal magnification is microlensing by intracluster stars, which induces independent flux variations in the images of individual or groups of source stars in the lensed galaxy. For a second possibility, intracluster dark matter subhalos invisible to telescopes effectively perturb lensing magnifications near the critical curve and give rise to persistently unequal image pairs. Our modeling suggests, at least for the most prominent identified image pair, that the microlensing hypothesis is in tension with the absence of notable asymmetry variation over a six-year baseline, while subhalos of $\sim 10^6$--$10^8\,M_\odot$ anticipated from structure formation with Cold Dark Matter typically produce stationary and sizable asymmetries. We judge that observations at additional times and more precise lens models are necessary to stringently constrain temporal variability and robustly distinguish between the two explanations. The arc under this study is a scheduled target of a Director's Discretionary Early Release Science program of the James Webb Space Telescope, which will provide deep images and a high-resolution view with integral field spectroscopy.

Axions are a viable candidate for Cold Dark Matter (CDM) which should generically form minihalos of sub-planetary masses from white-noise isocurvature density fluctuations if the Peccei-Quinn phase transition occurs after inflation. Despite being denser than the larger halos formed out of adiabatic fluctuations from inflation, axion minihalos have surface densities much smaller than the critical value required for gravitational lensing to produce multiple images or high magnification, and hence are practically undetectable as lenses in isolation. However, their lensing effect can be enhanced when superposed near critical curves of other lenses. We propose a method to detect them through photometric monitoring of recently discovered caustic transiting stars behind cluster lenses, under extreme magnification factors $\mu \gtrsim 10^3$--$10^4$ as the lensed stars cross microlensing caustics induced by intracluster stars. For masses of the first gravitationally collapsed minihalos in the range $\sim 10^{-15}$--$10^{-8}\,h^{-1}\,M_\odot$, we show that axion minihalos in galaxy clusters should collectively produce subtle surface density fluctuations of amplitude $\sim 10^{-4}$--$10^{-3}$ on projected length scales of $\sim 10$--$10^4\,$AU, which imprint irregularities in the microlensing light curves of caustic transiting stars. We estimate that, inside a cluster halo and over the age of the Universe, most of these minihalos are likely to avoid dynamic disruption by encounters with stars or other minihalos.

We examine the maximum possible strength of the global 21-cm absorption dip on the Cosmic Background Radiation at high-redshift caused by the atomic intergalactic medium, when the Lyman-$\alpha$ coupling is maximum, assuming no exotic cooling mechanisms from interactions with dark matter. This maximum absorption is limited by three inevitable factors that need to be accounted for: $(a)$ heating by energy transferred from the Cosmic Background Radiation to the hydrogen atoms via 21-cm transitions, dubbed as 21-cm heating; $(b)$ Ly$\alpha$ heating by scatterings of Ly$\alpha$ photons from the first stars; $(c)$ the impact of the expected density fluctuations in the intergalactic gas in standard Cold Dark Matter theory, which reduces the mean 21-cm absorption signal. Inclusion of this third novel effect reduces the maximum global 21-cm absorption by $\sim 10\%$. Overall, the three effects studied here reduce the 21-cm global absorption by $\sim 20\%$ at $z \simeq 17$.

The Southern Photometric Local Universe Survey (S-PLUS) is imaging ~9300 deg^2 of the celestial sphere in twelve optical bands using a dedicated 0.8 m robotic telescope, the T80-South, at the Cerro Tololo Inter-American Observatory, Chile. The telescope is equipped with a 9.2k by 9.2k e2v detector with 10 um pixels, resulting in a field-of-view of 2 deg^2 with a plate scale of 0.55"/pixel. The survey consists of four main subfields, which include two non-contiguous fields at high Galactic latitudes (8000 deg^2 at |b| > 30 deg) and two areas of the Galactic plane and bulge (for an additional 1300 deg^2). S-PLUS uses the Javalambre 12-band magnitude system, which includes the 5 u, g, r, i, z broad-band filters and 7 narrow-band filters centered on prominent stellar spectral features: the Balmer jump/[OII], Ca H+K, H-delta, G-band, Mg b triplet, H-alpha, and the Ca triplet. S-PLUS delivers accurate photometric redshifts (delta_z/(1+z) = 0.02 or better) for galaxies with r < 20 AB mag and redshift < 0.5, thus producing a 3D map of the local Universe over a volume of more than 1 (Gpc/h)^3. The final S-PLUS catalogue will also enable the study of star formation and stellar populations in and around the Milky Way and nearby galaxies, as well as searches for quasars, variable sources, and low-metallicity stars. In this paper we introduce the main characteristics of the survey, illustrated with science verification data highlighting the unique capabilities of S-PLUS. We also present the first public data release of ~336 deg^2 of the Stripe-82 area, which is available at http://datalab.noao.edu/splus.

A method to search for tidal streams and to fit their orbits based on maximum likelihood is presented and applied to the Gaia data. Tests of the method are performed showing how a simulated stream produced by tidal stripping of a star cluster is recovered when added to a simulation of the Gaia catalogue. The method can be applied to search for streams associated with known progenitors or to do blind searches in a general catalogue. As the first example, we apply the method to the globular cluster M68 and detect its clear tidal stream stretching over the whole North Galactic hemisphere, and passing within 5 kpc of the Sun. This is one of the closest tidal streams to us detected so far, and is highly promising to provide new constraints on the Milky Way gravitational potential, for which we present preliminary fits finding a slightly oblate dark halo consistent with other observations. We identify the M68 tidal stream with the previously discovered Fj\"orm stream by Ibata et al. The tidal stream is confirmed to contain stars that are consistent with the HR-diagram of M68. We provide a list of 115 stars that are most likely to be stream members, and should be prime targets for follow-up spectroscopic studies.

We examine a caustic-straddling arc at $z=0.9397$ in the field of the galaxy cluster MACS J0416.1-2403 ($z=0.397$) using archival multiband HST images and show that its surface brightness exhibits anomalies that can be explained by a single highly magnified star undergoing microlensing. First, we show that the surface brightness pattern is not perfectly symmetric across the cluster critical curve, which is inconsistent with a locally smooth lens model; the location of the candidate star exhibits the most significant asymmetry. Second, our analysis indicates that the asymmetric feature has $\sim 30\%$ higher flux in the 2012 visits compared to the Frontier Fields program visits in 2014. Moreover, the variable asymmetric feature shows an anomalous color between the F814W and F105W filters in 2014. These anomalies are naturally explained by microlensing induced variability of a caustic-transiting blue supergiant in a star-forming region, with a mean magnification factor around $\mu \sim 200$. We extend this study to a statistical analysis of the whole arc image and find tentative evidence of the increased mismatch of the two images in the proximity of the critical line. Robust detection of one or multiple caustic-transiting stars in this arc will enable detailed follow-up studies that can shed light on the small-scale structure of the dark matter inside the cluster halo.

Luminous stars in background galaxies straddling the lensing caustic of a foreground galaxy cluster can be individually detected due to extreme magnification factors of $\sim 10^2$--$10^3$, as recently observed in deep HST images. We propose a direct method to probe the presence of dark matter subhalos in galaxy clusters by measuring the astrometric perturbation they induce on the image positions of magnified stars or bright clumps: lensing by subhalos breaks the symmetry of a smooth critical curve, traced by the midpoints of close image pairs. For the giant arc at $z = 0.725$ behind the lensing cluster Abell 370 at $z = 0.375$, a promising target for detecting image pairs of stars, we find that subhalos of masses in the range $10^6$--$10^8\,M_\odot$ with the abundance predicted in the cold dark matter theory should typically imprint astrometric distortions at the level of $20$--$80\,{\rm mas}$. We estimate that $\sim 10\,$hr integrations with JWST at $\sim 1$--$3\,\mu{\rm m}$ may uncover several magnified stars whose image doublets will reveal the subhalo-induced structures of the critical curve. This method can probe a dynamic range in the subhalo to cluster halo mass ratio $m/M \sim 10^{-7}$--$10^{-9}$, thereby placing new constraints on the nature of dark matter.

We measure the cosmological bias factor of DLAs from their cross-correlation with the Ly$\alpha$ forest absorption, as a function of the DLA metal strength, defined from an average of equivalent widths of the strongest detectable low-ionization metal lines. A clear increase of the bias factor with metal strength is detected, as expected from a relation of metallicity and velocity dispersion with host halo mass. The relation is stronger after the metal strength is corrected for the HI column density, to make it more related to metallicity instead of metal column density. After correcting for the effects of measurement errors of the metal strength parameter, we find that the bias factor of DLAs with the weakest metal lines is close to unity, consistent with an origin in dwarf galaxies with host halo masses $\sim10^{10} M_{\odot}$, whereas the most metal rich DLAs have a bias factor as large as $b_{\rm DLA} \sim 3$, indicative of massive galaxies or galaxy groups in host halos with masses $\sim 10^{12} M_{\odot}$. Our result confirms the physical origin of the relation of bias factors measured from cross-correlation studies to the host halos of the absorbers.

We present a new determination of the large-scale clustering of the CIV forest (i.e., the absorption due to all CIV absorbers) using its cross-correlation with quasars in the Sloan Digital Sky Survey (SDSS) Data Release 12 (DR12). We fit a linear bias model to the measured cross-correlation. We find that the transmission bias of the CIV forest, $b_{Fc}$, at a mean redshift of $z=2.3$, obeys the relation $(1+\beta_c)b_{F c} = -0.024 \pm 0.003$. Here, $\beta_{c}$ is the linear redshift space distortion parameter of the CIV absorption, which can only be poorly determined at $\beta_c=1.1\pm 0.6$ from our data. This transmission bias is related to the bias of CIV absorbers and their host halos through the effective mean optical depth of the CIV forest, $\bar\tau_c$. Estimating a value $\bar \tau_c(z) \simeq 0.01$ from previous studies of the CIV equivalent width distribution, our measurement implies a CIV absorber bias near unity, with a large error due to uncertainties in both $\beta_c$ and $\bar\tau_c$. This makes it compatible with the higher DLA bias $b_{\rm DLA}\simeq 2$ measured previously from the cross-correlation of DLAs and the Lyman-$\alpha$ forest. We discuss the implications of the CIV absorber bias for the mass distribution of their host halos. More accurate determinations of $\bar \tau_c(z)$ and $\beta_c$ are necessary to obtain a more robust measurement of this CIV absorber bias.

We present the first paper of the series Origin of Metals around Galaxies (OMG) aimed to study the origin of the metals observed in the circumgalactic and intergalactic media. In this work we extract and build the catalogs of metal absorbers that will be used in future analyses, and make our results publicly available to the community. We design a fully automatic algorithm to search for absorption metal-line doublets of the species CIV, NV, SiIV and MgII in high-resolution ($R\gtrsim30\,000$) quasar spectra without human intervention, and apply it to the high-resolution and signal-to-noise ratio spectra of 690 quasars, observed with the UVES and HIRES instruments. We obtain $5\,656$ CIV doublets, $7\,919$ doublets of MgII, $2\,258$ of SiIV, and 239 of NV, constituting the largest high-resolution metal-doublet samples to date, and estimate the dependence of their completeness and purity on various doublet parameters such as equivalent width and redshift, using real and artificial quasar spectra. The catalogs include doublets with rest-frame line equivalent widths down to a few ${\rm m\AA}$, all detected at a significance above 3$\sigma$, and covering the redshifts between $1<z \lesssim 5$, properties that make them useful for a wide range of chemical evolution studies.

We investigate the large-scale structure of Lyman-alpha emission intensity in the Universe at redshifts z=2-3.5 using cross-correlation techniques. Our Lya emission samples are spectra of BOSS Luminous Red Galaxies from Data Release 12 with the best fit model galaxies subtracted. We cross-correlate the residual flux in these spectra with BOSS quasars, and detect a positive signal on scales 1-15 Mpc/h. We identify and remove a source of contamination not previously accounted for, due to the effects of quasar clustering on cross-fibre light. Corrected, our quasar-Lya emission cross-correlation is 50 % lower than that seen by Croft et al. for DR10, but still significant. Because only 3% of space is within 15 Mpc/h of a quasar, the result does not fully explore the global large-scale structure of Lya emission. To do this, we cross-correlate with the Lya forest. We find no signal in this case. The 95% upper limit on the global Lya mean surface brightness from Lya emission-Lya forest cross-correlation is mu < 1.2x10^-22 erg/s/cm^2/A/arcsec^2 This null result rules out the scenario where the observed quasar-Lya emission cross-correlation is primarily due to the large scale structure of star forming galaxies, Taken in combination, our results suggest that Lya emitting galaxies contribute, but quasars dominate within 15 Mpc/h. A simple model for Lya emission from quasars based on hydrodynamic simulations reproduces both the observed forest-Lya emission and quasar-Lya emission signals. The latter is also consistent with extrapolation of observations of fluorescent emission from smaller scales r < 1 Mpc.

Spinning dust grains in front of the bright Galactic synchrotron background can produce a weak absorption signal that could affect measurements of high redshift 21 cm absorption. At frequencies near 80 MHz where the EDGES experiment has reported 21\,cm absorption at $z \approx 17$, absorption could be produced by interstellar nanoparticles with radii $a \approx 50\AA$ in the cold interstellar medium at temperature $T \approx 50$ K. Atmospheric aerosols could contribute additional absorption. The strength of the absorption depends on the abundance of such grains and on their dipole moments, which are uncertain. The breadth of the absorption spectrum of spinning dust limits its possible impact on measurement of a relatively narrow 21 cm absorption feature.

J-PLUS is an ongoing 12-band photometric optical survey, observing thousands of square degrees of the Northern hemisphere from the dedicated JAST/T80 telescope at the Observatorio Astrof\'isico de Javalambre. T80Cam is a 2 sq.deg field-of-view camera mounted on this 83cm-diameter telescope, and is equipped with a unique system of filters spanning the entire optical range. This filter system is a combination of broad, medium and narrow-band filters, optimally designed to extract the rest-frame spectral features (the 3700-4000\AA\ Balmer break region, H$\delta$, Ca H+K, the G-band, the Mgb and Ca triplets) that are key to both characterize stellar types and to deliver a low-resolution photo-spectrum for each pixel of the sky observed. With a typical depth of AB $\sim 21.25$ mag per band, this filter set thus allows for an indiscriminate and accurate characterization of the stellar population in our Galaxy, it provides an unprecedented 2D photo-spectral information for all resolved galaxies in the local universe, as well as accurate photo-z estimates ($\Delta\,z\sim 0.01-0.03$) for moderately bright (up to $r\sim 20$ mag) extragalactic sources. While some narrow band filters are designed for the study of particular emission features ([OII]/$\lambda$3727, H$\alpha$/$\lambda$6563) up to $z < 0.015$, they also provide well-defined windows for the analysis of other emission lines at higher redshifts. As a result, J-PLUS has the potential to contribute to a wide range of fields in Astrophysics, both in the nearby universe (Milky Way, 2D IFU-like studies, stellar populations of nearby and moderate redshift galaxies, clusters of galaxies) and at high redshifts (ELGs at $z\approx 0.77, 2.2$ and $4.4$, QSOs, etc). With this paper, we release $\sim 36$ sq.deg of J-PLUS data, containing about $1.5\times 10^5$ stars and $10^5$ galaxies at $r<21$ mag.

An ultra-light axion field with mass $\sim 10^{-22}\ {\rm eV}$, also known as wave or fuzzy dark matter, has been proposed as a component of the dark matter in the Universe. We study the evolution of the axion dark matter distribution in the central region of a halo, assuming the mass is dominated by this axion field, and that gravity is the only important interaction. We calculate the excited axion states in the spherical gravitational potential generated by the self-gravitating ground-state, also known as soliton. These excited states are similar to the states of the hydrogen atom with quantum numbers $(n,l,m)$, here designating oscillation modes of a classical wave. At fixed $n$, the modes with highest $l$ have the lowest energy because of the extended mass distribution generating the potential. We use an approximate analytical treatment to derive the distribution of mass in these states when a steady-state is reached by dynamical relaxation, and find that a corona with a mass density profile $\rho \propto r^{-5/3}$ should be set up around the central soliton, analogous to the Bahcall-Wolf cusp predicted for the stellar distribution around a central black hole. The central soliton accretes dark matter from the corona as dynamical relaxation proceeds and negative orbital energy flows out. This density profile should remain valid out to the radius where the mass in the corona is comparable to the mass of the central soliton; further than that, the gravitational potential depends on the initial distribution of dark matter and the relaxation time increases rapidly with radius.

Recent observations of lensed galaxies at cosmological distances have detected individual stars that are extremely magnified when crossing the caustics of lensing clusters. In idealized cluster lenses with smooth mass distributions, two images of a star of radius $R$ approaching a caustic brighten as $t^{-1/2}$ and reach a peak magnification $\sim 10^{6}\, (10\, R_{\odot}/R)^{1/2}$ before merging on the critical curve. We show that a mass fraction ($\kappa_\star \gtrsim \, 10^{-4.5}$) in microlenses inevitably disrupts the smooth caustic into a network of corrugated microcaustics, and produces light curves with numerous peaks. Using analytical calculations and numerical simulations, we derive the characteristic width of the network, caustic-crossing frequencies, and peak magnifications. For the lens parameters of a recent detection and a population of intracluster stars with $\kappa_\star \sim 0.01$, we find a source-plane width of $\sim 20 \, {\rm pc}$ for the caustic network, which spans $0.2 \, {\rm arcsec}$ on the image plane. A source star takes $\sim 2\times 10^4$ years to cross this width, with a total of $\sim 6 \times 10^4$ crossings, each one lasting for $\sim 5\,{\rm hr}\,(R/10\,R_\odot)$ with typical peak magnifications of $\sim 10^{4} \left( R/ 10\,R_\odot \right)^{-1/2}$. The exquisite sensitivity of caustic-crossing events to the granularity of the lens-mass distribution makes them ideal probes of dark matter components, such as compact halo objects and ultralight axion dark matter.

We present a measurement of baryon acoustic oscillations (BAO) in the cross-correlation of quasars with the Ly$\alpha$-forest flux-transmission at a mean redshift $z=2.40$. The measurement uses the complete SDSS-III data sample: 168,889 forests and 234,367 quasars from the SDSS Data Release DR12. In addition to the statistical improvement on our previous study using DR11, we have implemented numerous improvements at the analysis level allowing a more accurate measurement of this cross-correlation. We also developed the first simulations of the cross-correlation allowing us to test different aspects of our data analysis and to search for potential systematic errors in the determination of the BAO peak position. We measure the two ratios $D_{H}(z=2.40)/r_{d} = 9.01 \pm 0.36$ and $D_{M}(z=2.40)/r_{d} = 35.7 \pm 1.7$, where the errors include marginalization over the non-linear velocity of quasars and the metal - quasar cross-correlation contribution, among other effects. These results are within $1.8\sigma$ of the prediction of the flat-$\Lambda$CDM model describing the observed CMB anisotropies. We combine this study with the Ly$\alpha$-forest auto-correlation function [2017A&A...603A..12B], yielding $D_{H}(z=2.40)/r_{d} = 8.94 \pm 0.22$ and $D_{M}(z=2.40)/r_{d} = 36.6 \pm 1.2$, within $2.3\sigma$ of the same flat-$\Lambda$CDM model.

We present a measurement of the DLA mean bias from the cross-correlation of DLA and the Ly$\alpha$ forest, updating earlier results of Font-Ribera et al. 2012 with the final BOSS Data Release and an improved method to address continuum fitting corrections. Our cross-correlation is well fitted by linear theory with the standard $\Lambda CDM$ model, with a DLA bias of $b_{\rm DLA} = 1.99\pm 0.11$; a more conservative analysis, which removes DLA in the Ly$\beta$ forest and uses only the cross-correlation at $r> 10{\rm h^{-1}\,Mpc}$, yields $b_{\rm DLA} = 2.00\pm 0.19$. This assumes the cosmological model from \cite{Planck2015} and the Ly$\alpha$ forest bias factors of Bautista et al. 2017, and includes only statistical errors obtained from bootstrap analysis. The main systematic errors arise from possible impurities and selection effects in the DLA catalogue, and from uncertainties in the determination of the Ly$\alpha$ forest bias factors and a correction for effects of high column density absorbers. We find no dependence of the DLA bias on column density or redshift. The measured bias value corresponds to a host halo mass $\sim 4\cdot10^{11} {\rm M_{\odot}}$ if all DLA were hosted in halos of a similar mass. In a realistic model where host halos over a broad mass range have a DLA cross section $\Sigma(M_h) \propto M_h^{\alpha}$ down to $M_h > M_{\rm min} =10^{8.5} {\rm M_{\odot}}$, we find that $\alpha > 1$ is required to have $b_{\rm DLA}> 1.7$, implying a steeper relation or higher value of $M_{\rm min}$ than is generally predicted in numerical simulations of galaxy formation.

We study the mean absorption spectrum of the Damped Lyman alpha population at $z\sim 2.6$ by stacking normalized, rest-frame shifted spectra of $\sim 27\,000$ DLAs from the DR12 of BOSS/SDSS-III. We measure the equivalent widths of 50 individual metal absorption lines in 5 intervals of DLA hydrogen column density, 5 intervals of DLA redshift, and overall mean equivalent widths for an additional 13 absorption features from groups of strongly blended lines. The mean equivalent width of low-ionization lines increases with $N_{\rm HI}$, whereas for high-ionization lines the increase is much weaker. The mean metal line equivalent widths decrease by a factor $\sim 1.1-1.5$ from $z\sim2.1$ to $z \sim 3.5$, with small or no differences between low- and high-ionization species. We develop a theoretical model, inspired by the presence of multiple absorption components observed in high-resolution spectra, to infer mean metal column densities from the equivalent widths of partially saturated metal lines. We apply this model to 14 low-ionization species and to AlIII, SIII, SiIII, CIV, SiIV, NV and OVI. We use an approximate derivation for separating the equivalent width contributions of several lines to blended absorption features, and infer mean equivalent widths and column densities from lines of the additional species NI, ZnII, CII${}^{*}$, FeIII, and SIV. Several of these mean column densities of metal lines in DLAs are obtained for the first time; their values generally agree with measurements of individual DLAs from high-resolution, high signal-to-noise ratio spectra when they are available.

The fraction of galaxies with strong Ly$\alpha$ emission has been observed to decrease rapidly with redshift at $z \ge 6$, after a gradual increase at $z< 6$. This has been interpreted as a hint of the reionization of the intergalactic medium (IGM): the emitted Ly$\alpha$ photons would be scattered by an increasingly neutral IGM at $z>6$. We study this effect by modeling the ionization and Ly$\alpha$ radiative transfer in the infall region and the IGM around a Ly$\alpha$ emitting galaxy (LAE), for a spherical halo model with the mean density and radial velocity profiles in the standard $\Lambda$CDM cosmological scenario. We find that the expected fast increase of the ionizing background intensity toward the end of the reionization epoch implies a rapid evolution of halo infall regions from being self-shielded against the external ionizing background to being mostly ionized. Whereas self-shielded infall regions can scatter the Ly$\alpha$ photons over a much larger area than the commonly used apertures for observing LAEs, the same infalling gas is no longer optically thick to the Ly$\alpha$ emission line after it is ionized by the external background, making the Ly$\alpha$ emission more compact and brighter within the observed apertures. Based on this simple model, we show that the observed drop in the abundance of LAEs at $z>6$ does not imply a rapid increase with redshift of the fraction of the whole IGM volume that is atomic, but is accounted for by a rapid increase of the neutral fraction in the infall regions around galaxy host halos.

We use flux-transmission correlations in \Lya forests to measure the imprint of baryon acoustic oscillations (BAO). The study uses spectra of 157,783 quasars in the redshift range $2.1\le z \le 3.5$ from the Sloan Digital Sky Survey (SDSS) Data Release 12 (DR12). Besides the statistical improvements on our previous studies using SDSS DR9 and DR11, we have implemented numerous improvements in the analysis procedure, allowing us to construct a physical model of the correlation function and to investigate potential systematic errors in the determination of the BAO peak position. The Hubble distance, $\DHub=c/H(z)$, relative to the sound horizon is $\DHub(z=2.33)/r_d=9.07 \pm 0.31$. The best-determined combination of comoving angular-diameter distance, $\DM$, and the Hubble distance is found to be $\DHub^{0.7}\DM^{0.3}/r_d=13.94\pm0.35$. This value is $1.028\pm0.026$ times the prediction of the flat-\lcdm model consistent with the cosmic microwave background (CMB) anisotropy spectrum. The errors include marginalization over the effects of unidentified high-density absorption systems and fluctuations in ultraviolet ionizing radiation. Independently of the CMB measurements, the combination of our results and other BAO observations determine the open-\lcdm density parameters to be $\om=0.296 \pm 0.029$, $\ol=0.699 \pm 0.100$ and $\Omega_k = -0.002 \pm 0.119$.

In these proceedings we describe the WEAVE-QSO survey, which will observe around 400,000 high redshift quasars starting in 2018. This survey is part of a broader WEAVE survey to be conducted at the 4.2m William Herschel Telescope. We will focus on chiefly on the science goals, but will also briefly summarise the target selection methods anticipated and the expected survey plan. Understanding the apparent acceleration in the expansion of the Universe is one of the key scientific challenges of our time. Many experiments have been proposed to study this expansion, using a variety of techniques. Here we describe a survey that can measure this acceleration and therefore help elucidate the nature of dark energy: a survey of the Lyman-alpha forest (and quasar absorption in general) in spectra towards z>2 quasars (QSOs). Further constraints on neutrino masses and warm dark matter are also anticipated. The same data will also shed light on galaxy formation via study of the properties of inflowing/outflowing gas associated with nearby galaxies and in a cosmic web context. Gas properties are sensitive to density, temperature, UV radiation, metallicity and abundance pattern, and so constraint galaxy formation in a variety of ways. WEAVE-QSO will study absorbers with a dynamic range spanning more than 8 orders of magnitude in column density, their thermal broadening, and a host of elements and ionization species. A core principal of the WEAVE-QSO survey is the targeting of QSOs with near 100% efficiency principally through use of the J-PAS (r < 23.2) and Gaia (r < 20) data.

In recent years, the Lyman-$\alpha$ absorption observed in the spectra of high-redshift quasars has been used as a tracer of large-scale structure by means of the three-dimensional Lyman-$\alpha$ forest auto-correlation function at redshift $z\simeq 2.3$, but the need to fit the quasar continuum in every absorption spectrum introduces a broadband distortion that is difficult to correct and causes a systematic error for measuring any broadband properties. We describe a $k$-space model for this broadband distortion based on a multiplicative correction to the power spectrum of the transmitted flux fraction that suppresses power on scales corresponding to the typical length of a Lyman-$\alpha$ forest spectrum. Implementing the distortion model in fits for the baryon acoustic oscillation (BAO) peak position in the Lyman-$\alpha$ forest auto-correlation, we find that the fitting method recovers the input values of the linear bias parameter $b_{F}$ and the redshift-space distortion parameter $\beta_{F}$ for mock data sets with a systematic error of less than 0.5\%. Applied to the auto-correlation measured for BOSS Data Release 11, our method improves on the previous treatment of broadband distortions in BAO fitting by providing a better fit to the data using fewer parameters and reducing the statistical errors on $\beta_{F}$ and the combination $b_{F}(1+\beta_{F})$ by more than a factor of seven. The measured values at redshift $z=2.3$ are $\beta_{F}=1.39^{+0.11\ +0.24\ +0.38}_{-0.10\ -0.19\ -0.28}$ and $b_{F}(1+\beta_{F})=-0.374^{+0.007\ +0.013\ +0.020}_{-0.007\ -0.014\ -0.022}$ (1$\sigma$, 2$\sigma$ and 3$\sigma$ statistical errors). Our fitting software and the input files needed to reproduce our main results are publicly available.

We describe mock data-sets generated to simulate the high-redshift quasar sample in Data Release 11 (DR11) of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). The mock spectra contain Ly{\alpha} forest correlations useful for studying the 3D correlation function including Baryon Acoustic Oscillations (BAO). They also include astrophysical effects such as quasar continuum diversity and high-density absorbers, instrumental effects such as noise and spectral resolution, as well as imperfections introduced by the SDSS pipeline treatment of the raw data. The Ly{\alpha} forest BAO analysis of the BOSS collaboration, described in Delubac et al. 2014, has used these mock data-sets to develop and cross-check analysis procedures prior to performing the BAO analysis on real data, and for continued systematic cross checks. Tests presented here show that the simulations reproduce sufficiently well important characteristics of real spectra. These mock data-sets will be made available together with the data at the time of the Data Release 11.

The Lyman alpha forest power spectrum has been measured on large scales by the BOSS survey in SDSS-III at $z\sim 2.3$, has been shown to agree well with linear theory predictions, and has provided the first measurement of Baryon Acoustic Oscillations at this redshift. However, the power at small scales, affected by non-linearities, has not been well examined so far. We present results from a variety of hydrodynamic simulations to predict the redshift space non-linear power spectrum of the Lyman Alpha transmission for several models, testing the dependence on resolution and box size. A new fitting formula is introduced to facilitate the comparison of our simulation results with observations and other simulations. The non-linear power spectrum has a generic shape determined by a transition scale from linear to non-linear anisotropy, and a Jeans scale below which the power drops rapidly. In addition, we predict the two linear bias factors of the Lyman Alpha forest and provide a better physical interpretation of their values and redshift evolution. The dependence of these bias factors and the non-linear power on the amplitude and slope of the primordial fluctuations power spectrum, the temperature-density relation of the intergalactic medium, and the mean Lyman Alpha transmission, as well as the redshift evolution, is investigated and discussed in detail. A preliminary comparison to the observations shows that the predicted redshift distortion parameter is in good agreement with the recent determination of Blomqvist et al., but the density bias factor is lower than observed. We make all our results publicly available in the form of tables of the non-linear power spectrum that is directly obtained from all our simulations, and parameters of our fitting formula.

We derive constraints on cosmological parameters and tests of dark energy models from the combination of baryon acoustic oscillation (BAO) measurements with cosmic microwave background (CMB) and Type Ia supernova (SN) data. We take advantage of high-precision BAO measurements from galaxy clustering and the Ly-alpha forest (LyaF) in the BOSS survey of SDSS-III. BAO data alone yield a high confidence detection of dark energy, and in combination with the CMB angular acoustic scale they further imply a nearly flat universe. Combining BAO and SN data into an "inverse distance ladder" yields a 1.7% measurement of $H_0=67.3 \pm1.1$ km/s/Mpc. This measurement assumes standard pre-recombination physics but is insensitive to assumptions about dark energy or space curvature, so agreement with CMB-based estimates that assume a flat LCDM cosmology is an important corroboration of this minimal cosmological model. For open LCDM, our BAO+SN+CMB combination yields $\Omega_m=0.301 \pm 0.008$ and curvature $\Omega_k=-0.003 \pm 0.003$. When we allow more general forms of evolving dark energy, the BAO+SN+CMB parameter constraints remain consistent with flat LCDM. While the overall $\chi^2$ of model fits is satisfactory, the LyaF BAO measurements are in moderate (2-2.5 sigma) tension with model predictions. Models with early dark energy that tracks the dominant energy component at high redshifts remain consistent with our constraints. Expansion history alone yields an upper limit of 0.56 eV on the summed mass of neutrino species, improving to 0.26 eV if we include Planck CMB lensing. Standard dark energy models constrained by our data predict a level of matter clustering that is high compared to most, but not all, observational estimates. (Abridged)

The third generation of the Sloan Digital Sky Survey (SDSS-III) took data from 2008 to 2014 using the original SDSS wide-field imager, the original and an upgraded multi-object fiber-fed optical spectrograph, a new near-infrared high-resolution spectrograph, and a novel optical interferometer. All the data from SDSS-III are now made public. In particular, this paper describes Data Release 11 (DR11) including all data acquired through 2013 July, and Data Release 12 (DR12) adding data acquired through 2014 July (including all data included in previous data releases), marking the end of SDSS-III observing. Relative to our previous public release (DR10), DR12 adds one million new spectra of galaxies and quasars from the Baryon Oscillation Spectroscopic Survey (BOSS) over an additional 3000 sq. deg of sky, more than triples the number of H-band spectra of stars as part of the Apache Point Observatory (APO) Galactic Evolution Experiment (APOGEE), and includes repeated accurate radial velocity measurements of 5500 stars from the Multi-Object APO Radial Velocity Exoplanet Large-area Survey (MARVELS). The APOGEE outputs now include measured abundances of 15 different elements for each star. In total, SDSS-III added 2350 sq. deg of ugriz imaging; 155,520 spectra of 138,099 stars as part of the Sloan Exploration of Galactic Understanding and Evolution 2 (SEGUE-2) survey; 2,497,484 BOSS spectra of 1,372,737 galaxies, 294,512 quasars, and 247,216 stars over 9376 sq. deg; 618,080 APOGEE spectra of 156,593 stars; and 197,040 MARVELS spectra of 5,513 stars. Since its first light in 1998, SDSS has imaged over 1/3 of the Celestial sphere in five bands and obtained over five million astronomical spectra.

(Abridged) We detect the large-scale structure of Lya emission in the Universe at redshifts z=2-3.5 by measuring the cross-correlation of Lya surface brightness with quasars in SDSS/BOSS. We use a million spectra targeting Luminous Red Galaxies at z<0.8, after subtracting a best fit model galaxy spectrum from each one, as an estimate of the high-redshift Lya surface brightness. The quasar-Lya emission cross-correlation we detect has a shape consistent with a LambdaCDM model with Omega_M =0.30^+0.10-0.07. The predicted amplitude of this cross-correlation is proportional to the product of the mean Lya surface brightness, <mu_alpha>, the amplitude of mass fluctuations, and the quasar and Lya emission bias factors. Using known values, we infer <mu_alpha>(b_alpha/3) = (3.9 +/- 0.9) x 10^-21 erg/s cm^-2 A^-1 arcsec^-2, where b_alpha is the Lya emission bias factor. If the dominant sources of Lya emission are star forming galaxies, we infer rho_SFR = (0.28 +/- 0.07) (3/b_alpha) /yr/Mpc^3 at z=2-3.5. For b_alpha=3, this value is a factor of 21-35 above previous estimates from individually detected Lya emitters, although consistent with the total rho_SFR derived from dust-corrected, continuum UV surveys. 97% of the Lya emission in the Universe at these redshifts is therefore undetected in previous surveys of Lya emitters. Our measurement is much greater than seen from stacking analyses of faint halos surrounding previously detected Lya emitters, but we speculate that it arises from similar Lya halos surrounding all luminous star-forming galaxies. We also detect redshift space anisotropy of the quasar-Lya emission cross-correlation, finding evidence at the 3.0 sigma level that it is radially elongated, consistent with distortions caused by radiative-transfer effects (Zheng et al. (2011)). Our measurements represent the first application of the intensity mapping technique to optical observations.

We present a measurement of the cross-correlation of MgII absorption and massive galaxies, using the DR11 main galaxy sample of the Baryon Oscillation Spectroscopic Survey of SDSS-III (CMASS galaxies), and the DR7 quasar spectra of SDSS-II. The cross-correlation is measured by stacking quasar absorption spectra shifted to the redshift of galaxies that are within a certain impact parameter bin of the quasar, after dividing by a quasar continuum model. This results in an average MgII equivalent width as a function of impact parameter from a galaxy, ranging from 50 kpc to more than 10 Mpc in proper units, which includes all MgII absorbers. We show that special care needs to be taken to use an unbiased quasar continuum estimator, to avoid systematic errors in the measurement of the mean stacked MgII equivalent width. The measured cross-correlation follows the expected shape of the galaxy correlation function, although measurement errors are large. We use the cross-correlation amplitude to derive the bias factor of MgII absorbers, finding bMgII = 2.33 \pm? 0.19, where the error accounts only for the statistical uncertainty in measuring the mean equivalent width. This bias factor is larger than that obtained in previous studies and may be affected by modeling uncertainties that we discuss, but if correct it suggests that MgII absorbers at redshift z \simeq 0:5 are spatially distributed on large scales similarly to the CMASS galaxies in BOSS. Keywords: galaxies: haloes, galaxies: formation, quasars: absorption lines, large-scale structure of universe

We report a detection of the baryon acoustic oscillation (BAO) feature in the flux-correlation function of the Ly{\alpha} forest of high-redshift quasars with a statistical significance of five standard deviations. The study uses 137,562 quasars in the redshift range $2.1\le z \le 3.5$ from the Data Release 11 (DR11) of the Baryon Oscillation Spectroscopic Survey (BOSS) of SDSS-III. This sample contains three times the number of quasars used in previous studies. The measured position of the BAO peak determines the angular distance, $D_A(z=2.34)$ and expansion rate, $H(z=2.34)$, both on a scale set by the sound horizon at the drag epoch, $r_d$. We find $D_A/r_d=11.28\pm0.65(1\sigma)^{+2.8}_{-1.2}(2\sigma)$ and $D_H/r_d=9.18\pm0.28(1\sigma)\pm0.6(2\sigma)$ where $D_H=c/H$. The optimal combination, $\sim D_H^{0.7}D_A^{0.3}/r_d$ is determined with a precision of $\sim2\%$. For the value $r_d=147.4~{\rm Mpc}$, consistent with the CMB power spectrum measured by Planck, we find $D_A(z=2.34)=1662\pm96(1\sigma)~{\rm Mpc}$ and $H(z=2.34)=222\pm7(1\sigma)~{\rm km\,s^{-1}Mpc^{-1}}$. Tests with mock catalogs and variations of our analysis procedure have revealed no systematic uncertainties comparable to our statistical errors. Our results agree with the previously reported BAO measurement at the same redshift using the quasar-Ly{\alpha} forest cross-correlation. The auto-correlation and cross-correlation approaches are complementary because of the quite different impact of redshift-space distortion on the two measurements. The combined constraints from the two correlation functions imply values of $D_A/r_d$ and $D_H/r_d$ that are, respectively, 7% low and 7% high compared to the predictions of a flat $\Lambda$CDM cosmological model with the best-fit Planck parameters. With our estimated statistical errors, the significance of this discrepancy is $\approx 2.5\sigma$.

We measure the large-scale cross-correlation of quasars with the Lyman alpha forest absorption, using over 164,000 quasars from Data Release 11 of the SDSS-III Baryon Oscillation Spectroscopic Survey. We extend the previous study of roughly 60,000 quasars from Data Release 9 to larger separations, allowing a measurement of the Baryonic Acoustic Oscillation (BAO) scale along the line of sight $c/(H(z=2.36) ~ r_s) = 9.0 \pm 0.3$ and across the line of sight $D_A(z=2.36) / ~ r_s = 10.8 \pm 0.4$, consistent with CMB and other BAO data. Using the best fit value of the sound horizon from Planck data ($r_s=147.49 Mpc$), we can translate these results to a measurement of the Hubble parameter of $H(z=2.36) = 226 \pm 8 km/s / Mpc$ and of the angular diameter distance of $D_A(z=2.36) = 1590 \pm 60 Mpc$. The measured cross-correlation function and an update of the code to fit the BAO scale (baofit) are made publicly available.

An analytical framework is presented to understand the effects of a fluctuating intensity of the cosmic ionising background on the correlations of the Ly{\alpha} forest transmission fraction measured in quasar spectra. In the absence of intensity fluctuations, the Ly{\alpha} power spectrum should have the expected cold dark matter power spectrum with redshift distortions in the linear regime, with a bias factor b_{\delta} and a redshift distortion parameter {\beta} that depend on redshift but are independent of scale. The intensity fluctuations introduce a scale dependence in both b_{\delta} and {\beta}, but keeping their product b_{\delta}{\beta} fixed. Observations of the Ly{\alpha} correlations and cross-correlations with radiation sources like those being done at present in the BOSS survey of SDSS-III (Busca et al. 2013; Slosar et al. 2013; Font-Ribera et al. 2014) have the potential to measure this scale dependence, which reflects the biasing properties of the sources and absorbers of the ionising background. We also compute a second term affecting the Ly{\alpha} spectrum, due to shot noise in the sources of radiation. This term is very large if luminous quasars are assumed to produce the ionising background and to emit isotropically with a constant luminosity, but should be reduced by a contribution from galaxies, and by the finite lifetime and anisotropic emission of quasars.

The Sloan Digital Sky Survey (SDSS) has been in operation since 2000 April. This paper presents the tenth public data release (DR10) from its current incarnation, SDSS-III. This data release includes the first spectroscopic data from the Apache Point Observatory Galaxy Evolution Experiment (APOGEE), along with spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS) taken through 2012 July. The APOGEE instrument is a near-infrared R~22,500 300-fiber spectrograph covering 1.514--1.696 microns. The APOGEE survey is studying the chemical abundances and radial velocities of roughly 100,000 red giant star candidates in the bulge, bar, disk, and halo of the Milky Way. DR10 includes 178,397 spectra of 57,454 stars, each typically observed three or more times, from APOGEE. Derived quantities from these spectra (radial velocities, effective temperatures, surface gravities, and metallicities) are also included.DR10 also roughly doubles the number of BOSS spectra over those included in the ninth data release. DR10 includes a total of 1,507,954 BOSS spectra, comprising 927,844 galaxy spectra; 182,009 quasar spectra; and 159,327 stellar spectra, selected over 6373.2 square degrees.

We present the Data Release 10 Quasar (DR10Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the first 2.5 years of the survey and that are confirmed as quasars via visual inspection of the spectra. The catalog also includes known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. The catalog contains 166,583 quasars (74,454 are new discoveries since SDSS-DR9) detected over 6,373 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra. The number of quasars with $z>2.15$ (117,668) is $\sim$5 times greater than the number of $z>2.15$ quasars known prior to BOSS. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII, MgII). The catalog identifies 16,461 broad absorption line quasars and gives their characteristics. For each object, the catalog presents five-band (u, g, r, i, z) CCD-based photometry with typical accuracy of 0.03 mag and information on the optical morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra cover the wavelength region 3,600-10,500\AA\ at a spectral resolution in the range 1,300$<$R$<$2,500; the spectra can be retrieved from the SDSS Catalog Archive Server. We also provide a supplemental list of an additional 2,376 quasars that have been identified among the galaxy targets of the SDSS-III/BOSS.

We have developed two independent methods to measure the one-dimensional power spectrum of the transmitted flux in the Lyman-$\alpha$ forest. The first method is based on a Fourier transform, and the second on a maximum likelihood estimator. The two methods are independent and have different systematic uncertainties. The determination of the noise level in the data spectra was subject to a novel treatment, because of its significant impact on the derived power spectrum. We applied the two methods to 13,821 quasar spectra from SDSS-III/BOSS DR9 selected from a larger sample of over 60,000 spectra on the basis of their high quality, large signal-to-noise ratio, and good spectral resolution. The power spectra measured using either approach are in good agreement over all twelve redshift bins from $<z> = 2.2$ to $<z> = 4.4$, and scales from 0.001 $\rm(km/s)^{-1}$ to $0.02 \rm(km/s)^{-1}$. We determine the methodological and instrumental systematic uncertainties of our measurements. We provide a preliminary cosmological interpretation of our measurements using available hydrodynamical simulations. The improvement in precision over previously published results from SDSS is a factor 2--3 for constraints on relevant cosmological parameters. For a $\Lambda$CDM model and using a constraint on $H_0$ that encompasses measurements based on the local distance ladder and on CMB anisotropies, we infer $\sigma_8 =0.83\pm0.03$ and $n_s= 0.97\pm0.02$ based on \ion{H}{i} absorption in the range $2.1<z<3.7$.

We measure the large-scale cross-correlation of quasars with the Lyman \alpha\ forest absorption in redshift space, using ~ 60000 quasar spectra from Data Release 9 (DR9) of the Baryon Oscillation Spectroscopic Survey (BOSS). The cross-correlation is detected over a wide range of scales, up to comoving separations r of 80 Mpc/h. For r > 15 Mpc/h, we show that the cross-correlation is well fitted by the linear theory prediction for the mean overdensity around a quasar host halo in the standard \Lambda CDM model, with the redshift distortions indicative of gravitational evolution detected at high confidence. Using previous determinations of the Lyman \alpha\ forest bias factor obtained from the Lyman \alpha\ autocorrelation, we infer the quasar bias factor to be b_q = 3.64^+0.13_-0.15 at a mean redshift z=2.38, in agreement with previous measurements from the quasar auto-correlation. We also obtain a new estimate of the Lyman \alpha\ forest redshift distortion factor, \beta_F = 1.1 +/- 0.15, slightly larger than but consistent with the previous measurement from the Lyman \alpha\ forest autocorrelation. The simple linear model we use fails at separations r < 15 Mpc/h, and we show that this may reasonably be due to the enhanced ionization due to radiation from the quasars. We also provide the expected correction that the mass overdensity around the quasar implies for measurements of the ionizing radiation background from the line-of-sight proximity effect.

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 report a detection of the baryon acoustic oscillation (BAO) feature in the three-dimensional correlation function of the transmitted flux fraction in the \Lya forest of high-redshift quasars. The study uses 48,640 quasars in the redshift range $2.1\le z \le 3.5$ from the Baryon Oscillation Spectroscopic Survey (BOSS) of the third generation of the Sloan Digital Sky Survey (SDSS-III). At a mean redshift $z=2.3$, we measure the monopole and quadrupole components of the correlation function for separations in the range $20\hMpc<r<200\hMpc$. A peak in the correlation function is seen at a separation equal to $(1.01\pm0.03)$ times the distance expected for the BAO peak within a concordance $\Lambda$CDM cosmology. This first detection of the BAO peak at high redshift, when the universe was strongly matter dominated, results in constraints on the angular diameter distance $\da$ and the expansion rate $H$ at $z=2.3$ that, combined with priors on $H_0$ and the baryon density, require the existence of dark energy. Combined with constraints derived from Cosmic Microwave Background (CMB) observations, this result implies $H(z=2.3)=(224\pm8){\rm km\,s^{-1}Mpc^{-1}}$, indicating that the time derivative of the cosmological scale parameter $\dot{a}=H(z=2.3)/(1+z)$ is significantly greater than that measured with BAO at $z\sim0.5$. This demonstrates that the expansion was decelerating in the range $0.7<z<2.3$, as expected from the matter domination during this epoch. Combined with measurements of $H_0$, one sees the pattern of deceleration followed by acceleration characteristic of a dark-energy dominated universe.

We describe fitting methods developed to analyze fluctuations in the Lyman-{\alpha} forest and measure the parameters of baryon acoustic oscillations (BAO). We apply our methods to BOSS Data Release 9. Our method is based on models of the three-dimensional correlation function in physical coordinate space, and includes the effects of redshift-space distortions, anisotropic non-linear broadening, and broadband distortions. We allow for independent scale factors along and perpendicular to the line of sight to minimize the dependence on our assumed fiducial cosmology and to obtain separate measurements of the BAO angular and relative velocity scales. Our fitting software and the input files needed to reproduce our main BOSS Data Release 9 results are publicly available.

We present the first measurement of the large-scale cross-correlation of Lyman alpha forest absorption and Damped Lyman alpha systems (DLA), using the 9th Data Release of the Baryon Oscillation Spectroscopic Survey (BOSS). The cross-correlation is clearly detected on scales up to 40 Mpc/h and is well fitted by the linear theory prediction of the standard Cold Dark Matter model of structure formation with the expected redshift distortions, confirming its origin in the gravitational evolution of structure. The amplitude of the DLA-Lyman alpha cross-correlation depends on only one free parameter, the bias factor of the DLA systems, once the Lyman alpha forest bias factors are known from independent Lyman alpha forest correlation measurements. We measure the DLA bias factor to be b_D = (2.17 +/- 0.20) beta_F^{0.22}, where the Lyman alpha forest redshift distortion parameter beta_F is expected to be above unity. This bias factor implies a typical host halo mass for DLAs that is much larger than expected in present DLA models, and is reproduced if the DLA cross section scales with halo mass as M_h^alpha, with alpha= 1.1 +/- 0.1 for beta_F=1. Matching the observed DLA bias factor and rate of incidence requires that atomic gas remains extended in massive halos over larger areas than predicted in present simulations of galaxy formation, with typical DLA proper sizes larger than 20 kpc in host halos of masses ~ 10^12 solar masses. We infer that typical galaxies at z ~ 2 to 3 are surrounded by systems of atomic clouds that are much more extended than the luminous parts of galaxies and contain ~ 10% of the baryons in the host halo.

The Baryon Oscillation Spectroscopic Survey (BOSS) is designed to measure the scale of baryon acoustic oscillations (BAO) in the clustering of matter over a larger volume than the combined efforts of all previous spectroscopic surveys of large scale structure. BOSS uses 1.5 million luminous galaxies as faint as i=19.9 over 10,000 square degrees to measure BAO to redshifts z<0.7. Observations of neutral hydrogen in the Lyman alpha forest in more than 150,000 quasar spectra (g<22) will constrain BAO over the redshift range 2.15<z<3.5. Early results from BOSS include the first detection of the large-scale three-dimensional clustering of the Lyman alpha forest and a strong detection from the Data Release 9 data set of the BAO in the clustering of massive galaxies at an effective redshift z = 0.57. We project that BOSS will yield measurements of the angular diameter distance D_A to an accuracy of 1.0% at redshifts z=0.3 and z=0.57 and measurements of H(z) to 1.8% and 1.7% at the same redshifts. Forecasts for Lyman alpha forest constraints predict a measurement of an overall dilation factor that scales the highly degenerate D_A(z) and H^{-1}(z) parameters to an accuracy of 1.9% at z~2.5 when the survey is complete. Here, we provide an overview of the selection of spectroscopic targets, planning of observations, and analysis of data and data quality of BOSS.

We present the Data Release 9 Quasar (DR9Q) catalog from the Baryon Oscillation Spectroscopic Survey (BOSS) of the Sloan Digital Sky Survey III. The catalog includes all BOSS objects that were targeted as quasar candidates during the survey, are spectrocopically confirmed as quasars via visual inspection, have luminosities Mi[z=2]<-20.5 (in a $\Lambda$CDM cosmology with H0 = 70 km/s/Mpc, $\Omega_{\rm M}$ = 0.3, and $\Omega_{\Lambda}$ = 0.7) and either display at least one emission line with full width at half maximum (FWHM) larger than 500 km/s or, if not, have interesting/complex absorption features. It includes as well, known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. This catalog contains 87,822 quasars (78,086 are new discoveries) detected over 3,275 deg$^{2}$ with robust identification and redshift measured by a combination of principal component eigenspectra newly derived from a training set of 8,632 spectra from SDSS-DR7. The number of quasars with $z>2.15$ (61,931) is ~2.8 times larger than the number of z>2.15 quasars previously known. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 7,533 broad absorption line quasars and gives their characteristics. For each object the catalog presents five-band (u,g,r,i,z) CCD-based photometry with typical accuracy of 0.03 mag, and information on the morphology and selection method. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys.

We present the first results from an ongoing survey for Damped Lyman-alpha systems (DLAs) in the spectra of z>2 quasars observed in the course of the Baryon Oscillation Spectroscopic Survey (BOSS), which is part of the Sloan Digital Sky Survey (SDSS) III. Our full (non-statistical) sample, based on Data Release 9, comprises 12,081 systems with log N(HI)>=20, out of which 6,839 have log N(HI)>=20.3. This is the largest DLA sample ever compiled, superseding that from SDSS-II by a factor of seven. Using a statistical sub-sample and estimating systematics from realistic mock data, we probe the N(HI) distribution at <z> = 2.5. Contrary to what is generally believed, the distribution extends beyond 10^22 cm^-2 with a moderate slope of index\approx-3.5. This result matches surprisingly well the opacity-corrected distribution observed at z = 0. The cosmological mass density of neutral gas in DLAs is found to be Omega_g_DLA~10^-3, evolving only mildly over the past 12 billion years.

The Sloan Digital Sky Survey III (SDSS-III) presents the first spectroscopic data from the Baryon Oscillation Spectroscopic Survey (BOSS). This ninth data release (DR9) of the SDSS project includes 535,995 new galaxy spectra (median z=0.52), 102,100 new quasar spectra (median z=2.32), and 90,897 new stellar spectra, along with the data presented in previous data releases. These spectra were obtained with the new BOSS spectrograph and were taken between 2009 December and 2011 July. In addition, the stellar parameters pipeline, which determines radial velocities, surface temperatures, surface gravities, and metallicities of stars, has been updated and refined with improvements in temperature estimates for stars with T_eff<5000 K and in metallicity estimates for stars with [Fe/H]>-0.5. DR9 includes new stellar parameters for all stars presented in DR8, including stars from SDSS-I and II, as well as those observed as part of the SDSS-III Sloan Extension for Galactic Understanding and Exploration-2 (SEGUE-2). The astrometry error introduced in the DR8 imaging catalogs has been corrected in the DR9 data products. The next data release for SDSS-III will be in Summer 2013, which will present the first data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) along with another year of data from BOSS, followed by the final SDSS-III data release in December 2014.

We present a study of the effect of High Column Density (HCD) systems on the Lyman \alpha forest correlation function on large scales. We study the effect both numerically, by inserting HCD systems on mock spectra for a specific model, and analytically, in the context of two-point correlations and linear theory. We show that the presence of HCDs substantially contributes to the noise of the correlation function measurement, and systematically alters the measured redshift-space correlation function of the Lyman \alpha forest, increasing the value of the density bias factor and decreasing the redshift distortion parameter $\beta_\alpha$ of the Lyman \alpha forest. We provide simple formulae for corrections on these derived parameters, as a function of the mean effective optical depth and bias factor of the host halos of the HCDs, and discuss the conditions under which these expressions should be valid. In practice, precise corrections to the measured parameters of the Lyman \alpha forest correlation for the HCD effects are more complex than the simple analytical approximations we present, owing to non-linear effects of the damped wings of the HCD systems and the presence of three-point terms. However, we conclude that an accurate correction for these HCD effects can be obtained numerically and calibrated with observations of the HCD-Lyman \alpha cross-correlation. We also discuss an analogous formalism to treat and correct for the contaminating effect of metal lines overlapping the Lyman \alpha forest spectra.

Massive spectroscopic surveys of high-redshift quasars yield large numbers of correlated Lyman {\alpha} absorption spectra that can be used to measure large-scale structure. Simulations of these surveys are required to accurately interpret the measurements of correlations and correct for systematic errors. An efficient method to generate mock realizations of Lyman {\alpha} forest surveys is presented which generates a field over the lines of sight to the survey sources only, instead of having to generate it over the entire three-dimensional volume of the survey. The method can be calibrated to reproduce the power spectrum and one-point distribution function of the transmitted flux fraction, as well as the redshift evolution of these quantities, and is easily used for modeling any survey systematic effects. We present an example of how these mock surveys are applied to predict the measurement errors in a survey with similar parameters as the BOSS quasar survey in SDSS-III.

Using a sample of approximately 14,000 z>2.1 quasars observed in the first year of the Baryon Oscillation Spectroscopic Survey (BOSS), we measure the three-dimensional correlation function of absorption in the Lyman-alpha forest. The angle-averaged correlation function of transmitted flux (F = exp(-tau)) is securely detected out to comoving separations of 60 Mpc/h, the first detection of flux correlations across widely separated sightlines. A quadrupole distortion of the redshift-space correlation function by peculiar velocities, the signature of the gravitational instability origin of structure in the Lyman-alpha forest, is also detected at high significance. We obtain a good fit to the data assuming linear theory redshift-space distortion and linear bias of the transmitted flux, relative to the matter fluctuations of a standard LCDM cosmological model (inflationary cold dark matter with a cosmological constant). At 95% confidence, we find a linear bias parameter 0.16<b<0.24 and redshift-distortion parameter 0.44<beta<1.20, at central redshift z=2.25, with a well constrained combination b(1+\beta)=0.336 +/- 0.012. The errors on beta are asymmetric, with beta=0 excluded at over 5 sigma confidence level. The value of beta is somewhat low compared to theoretical predictions, and our tests on synthetic data suggest that it is depressed (relative to expectations for the Lyman-alpha forest alone) by the presence of high column density systems and metal line absorption. These results set the stage for cosmological parameter determinations from three-dimensional structure in the Lyman-alpha forest, including anticipated constraints on dark energy from baryon acoustic oscillations.

We study whether non-linear gravitational effects of relic neutrinos on the development of clustering and large-scale structure may be observable by weak gravitational lensing. We compute the density profile of relic massive neutrinos in a spherical model of a cluster of galaxies, for several neutrino mass schemes and cluster masses. Relic neutrinos add a small perturbation to the mass profile, making it more extended in the outer parts. In principle, this non-linear neutrino perturbation is detectable in an all-sky weak lensing survey such as EUCLID by averaging the shear profile of a large fraction of the visible massive clusters in the universe, or from its signature in the general weak lensing power spectrum or its cross-spectrum with galaxies. However, correctly modeling the distribution of mass in baryons and cold dark matter and suppressing any systematic errors to the accuracy required for detecting this neutrino perturbation is severely challenging.

The Sloan Digital Sky Survey (SDSS) started a new phase in August 2008, with new instrumentation and new surveys focused on Galactic structure and chemical evolution, measurements of the baryon oscillation feature in the clustering of galaxies and the quasar Ly alpha forest, and a radial velocity search for planets around ~8000 stars. This paper describes the first data release of SDSS-III (and the eighth counting from the beginning of the SDSS). The release includes five-band imaging of roughly 5200 deg^2 in the Southern Galactic Cap, bringing the total footprint of the SDSS imaging to 14,555 deg^2, or over a third of the Celestial Sphere. All the imaging data have been reprocessed with an improved sky-subtraction algorithm and a final, self-consistent photometric recalibration and flat-field determination. This release also includes all data from the second phase of the Sloan Extension for Galactic Understanding and Evolution (SEGUE-2), consisting of spectroscopy of approximately 118,000 stars at both high and low Galactic latitudes. All the more than half a million stellar spectra obtained with the SDSS spectrograph have been reprocessed through an improved stellar parameters pipeline, which has better determination of metallicity for high metallicity stars.

Although the powering mechanism for quasars is now widely recognized to be the accretion of matter in a geometrically thin disk, the transport of matter to the inner region of the disk where luminosity is emitted remains an unsolved question. Miralda-Escud\'e & Kollmeier (2005) proposed a model whereby quasars are fuelled when stars are captured by the accretion disk as they plunge through the gas. Such plunging stars can then be destroyed and deliver their mass to the accretion disk. Here we present the first detailed calculations for the capture of stars originating far from the accretion disk near the zone of influence of the central black hole. In particular we examine the effect of adding a perturbing mass to a fixed stellar cusp potential on bringing stars into the accretion disk where they can be captured. The work presented here will be discussed in detail in an upcoming publication Kennedy et al. (2010).

Measurements of the baryonic acoustic oscillation (BAO) peak in the redshift-space correlation function yield the angular diameter distance D_A(z) and the Hubble parameter H(z) as a function of redshift, constraining the properties of dark energy and space curvature. We discuss the perturbations introduced in the galaxy correlation function by gravitational lensing through the effect of magnification bias and its cross-correlation with the galaxy density. At the BAO scale, gravitational lensing adds a small and slowly varying component to the galaxy correlation function and does not change its shape significantly, through which the BAO peak is measured. The relative shift in the position of the BAO peak caused by gravitational lensing in the angle-averaged correlation function is 10^-4 at z=1, rising to 10^-3 at z=2.5. Lensing effects are stronger near the line-of-sight, however the relative peak shift increases only to 10^-3.3 and 10^-2.4 at z=1 and z=2.5, when the galaxy correlation is averaged within 5 degrees of the line-of-sight (containing only 0.4% of the galaxy pairs in a survey). Furthermore, the lensing contribution can be measured separately and subtracted from the observed correlation at the BAO scale.

Astronomical observations have shown that the expansion of the universe is at present accelerating, consistently with a constant negative pressure or tension. This is a major puzzle because we do not understand why this tension is so small compared to the Planck density; why, being so small, it is not exactly zero; and why it has precisely the required value to make the expansion start accelerating just at the epoch when we are observing the universe. The recently proposed conjecture by Afshordi that black holes create a gravitational aether owing to quantum gravity effects, which may be identified with this invisible tension, can solve this coincidence problem. The fact that the expansion of the universe is starting to accelerate at the epoch when we observe it is a necessity that is implied by our origin in a planet orbiting a star that formed when the age of the universe was of the same order as the lifetime of the star. This argument is unrelated to any anthropic reasoning.

We present predictions for the fluorescent Lyman-alpha emission signature arising from photoionized, optically thick structures in Smoothed Particle Hydrodynamic (SPH) cosmological simulations of a Lambda-CDM universe using a Monte Carlo Lyman-alpha radiative transfer code. We calculate the expected Lyman-alpha image and 2-dimensional spectra for gas exposed to a uniform ultraviolet ionizing background as well as gas exposed additionally to the photoionizing radiation from a local quasar, after correcting for the self-shielding of hydrogen. As a test of our numerical methods and for application to current observations, we examine simplified analytic structures that are uniformly or anisotropically illuminated. We compare these results with recent observations. We discuss future observing campaigns on large telescopes and realistic strategies for detecting fluorescence owing to the ambient metagalactic ionization and in regions close to bright quasars. While it will take hundreds of hours on the current generation of telescopes to detect fluorescence caused by the ultraviolet background (UVB) alone, our calculations suggest that of order ten sources of quasar-induced fluorescent Lyman-alpha emission should be detectable after a 10 hour exposure in a 10 arcmin^2 field around a bright quasar. These observations will help probe the physical conditions in the densest regions of the intergalactic medium as well as the temporal light curves and isotropy of quasar radiation.

We revisit the well known discrepancy between the observed number of Milky Way (MW) dwarf satellite companions and the predicted population of cold dark matter (CDM) sub-halos, in light of the dozen new low luminosity satellites found in SDSS imaging data and our recent calibration of the SDSS satellite detection efficiency, which implies a total population far larger than these dozen discoveries. We combine a dynamical model for the CDM sub-halo population with simple, physically motivated prescriptions for assigning stellar content to each sub-halo, then apply observational selection effects and compare to the current observational census. As expected, models in which the stellar mass is a constant fraction F(Omega_b/Omega_m) of the sub-halo mass M_sat at the time it becomes a satellite fail for any choice of F. However, previously advocated models that invoke suppression of gas accretion after reionization in halos with circular velocity v_c <~ 35 km/s can reproduce the observed satellite counts for -15 < M_V < 0, with F ~ 10^{-3}. Successful models also require strong suppression of star formation BEFORE reionization in halos with v_c <~ 10 km/s; models without pre-reionization suppression predict far too many satellites with -5 < M_V < 0. Our models also reproduce the observed stellar velocity dispersions ~ 5-10 km/s of the SDSS dwarfs given the observed sizes of their stellar distributions, and model satellites have M(<300 pc) ~ 10^7 M_sun as observed even though their present day total halo masses span more than two orders of magnitude. Our modeling shows that natural physical mechanisms acting within the CDM framework can quantitatively explain the properties of the MW satellite population as it is presently known, thus providing a convincing solution to the `missing satellite' problem.

A new model for the formation of Jovian planets is proposed. We consider planets forming at large distances from a protostar (>100 AU) through direct fragmentation of a gas cloud, by the same formation mechanism as wide stellar and brown dwarf binaries. We model the gravitational evolution of a system of these distant planets and a second population formed in a disk closer to the star. We compute the typical closest approach of these planets to the star (i.e., smallest pericenter) over the course of their evolution. When the planets reach a pericenter within a gaseous disk surroundig the star, dynamical friction from this disk slows down the planet at each plunge, causing its orbit to be gradually circularized and made coplanar with the disk. After the disk dissipates, a large fraction of these planets may be left at orbits small enough to be detected in present radial velocity surveys. A brief analytic derivation of the rate of orbital energy dissipation during these disk crossings is presented. Observational tests of this model are discussed.

We present predictions of the correlation between the Lyman-alpha forest absorption in quasar spectra and the mass within \sim 5 Mpc/h (comoving) of the line of sight, using fully hydrodynamic and hydro-PM numerical simulations of the cold dark matter model supported by present observations. The observed correlation based on galaxies and the Lya forest can be directly compared to our theoretical results, assuming that galaxies are linearly biased on large scales. Specifically, we predict the average value of the mass fluctuation, <delta_m>, conditioned to a fixed value of the Lya forest transmitted flux delta_F, after they have been smoothed over a 10 Mpc/h cube and line of sight interval, respectively. We find that <delta_m>/sigma_m as a function of delta_F/sigma_F has a slope of 0.6 at this smoothing scale, where sigma_m and sigma_F are the rms dispersions (this slope should decrease with the smoothing scale). We show that this value is largely insensitive to the cosmological model and other Lya forest parameters. Comparison of our predictions to observations should provide a fundamental test of our ideas on the nature of the Lya forest and the distribution of galaxies, and can yield a measurement of the bias factor of any type of galaxies that are observed in the vicinity of Lya forest lines of sight.

Lyman-alpha forest lines in QSO spectra have typical widths of 20-50 km/s. Low column density absorbers in cosmological simulations are large, diffuse structures, and the Hubble flow across the spatially extended absorber is usually the dominant contribution to the width of its associated absorption line. Thermal broadening is unimportant over most of the spectrum, and peculiar velocities tend to make absorption features narrower rather than broader. As a consequence of Hubble flow broadening, there is a close relation between local Lyman-alpha optical depth and local neutral hydrogen density, which is well approximated by the Gunn-Peterson formula. The physics that governs the unshocked intergalactic medium leads to a tight correlation between the neutral hydrogen density and the underlying gas and dark matter overdensity. For many purposes, it is simpler to regard a Lyman-alpha forest spectrum as a continuous, non-linear map of the density field rather than a collection of discrete lines. This continuous field view of the Lyman-alpha forest can be applied to measurement of the baryon density parameter, testing of cosmological models, and robust determination of the shape and amplitude of the primordial mass power spectrum.

We examine the correlational property of $\lya$ clouds in detail and compare it to that of mass and galaxies. We show that the correlation strength of $\lya$ clouds is somewhat weaker than that of the underlying matter, which in turn is weaker than that of galaxies (biased galaxy formation). On the scales probed, $10-300$km/s, for $\lya$ clouds, we find that higher density, higher optical depth, higher column density regions are more strongly clustered than lower density, lower optical depth, lower column density regions. Thus, a consistent picture seems to emerge: the correlation strength for a given set of objects is positively correlated with their characteristic global density and the differences among the correlations of galaxies, $\lya$ clouds and mass reflect the differences in density that each trace. Significant positive correlations with a strength of $0.1-1.0$ are found for $\lya$ clouds in the velocity range $50-300$km/s. Among the correlational measures examined, an optical depth correlation function (Equation 5) proposed here may serve as the best correlational measure. It reasonably faithfully represents the true correlation of the underlying matter, enabling a better indication of both matter correlation and the relationship between galaxies and $\lya$ clouds. Furthermore, it appears to be an alternative to the conventional line-line correlation function with the virtue that it does not require ambiguous post-observation fitting procedures such as those commonly employed in the conventional line-finding methods. Neither does it depend sensitively on the observational resolution (e.g., FWHM), insofaras the clouds are resolved Conveniently, it can be easily measured with the current observational sensitivity without being contaminated significantly by the presence of noise.

Observations of gravitational lensing indicate that the mass distribution in clusters of galaxies (where most of the mass is dark matter) is highly peaked towards the center, while X-ray observations imply that the gas is more extended. The additional fact that the gas is cooling in the center has often led one to expect that the gas temperature should be lower near the center, and therefore the gas should be more concentrated than the dark matter. We show that such expectation is not correct, and that the gas temperature must remain approximately constant within the cooling region in order to have consistency with the observed X-ray profiles and lenses. A multiphase cooling flow naturally produces an approximately constant temperature profile, and a more extended distribution for the gas compared to the mass. Cool phases are deposited at relatively large radius, while hot phases are adiabatically heated as they flow inwards and can keep the average temperature constant. Thus, cooling flows result in an {\it increase} of the central temperature, relative to a case where there is no cooling and the gas follows the mass distribution. The increased central temperatures caused by cooling flows give a characteristic core radius to the gas profiles, which is of order the cooling radius. This provides a natural explanation for the typical cores observed in X-ray clusters. It also brings into better agreement with observations the rate of cluster evolution expected in self-similar hierarchical models. We propose that clusters having core radii much larger than their cooling radii are in the process of merging and are not in dynamical equilibrium.