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.

Observations of the redshift z=7.085 quasar J1120+0641 have been used to search for variations of the fine structure constant, alpha, over the redshift range 5.5 to 7.1. Observations at z=7.1 probe the physics of the universe when it was only 0.8 billion years old. These are the most distant direct measurements of alpha to date and the first measurements made with a near-IR spectrograph. A new AI analysis method has been employed. Four measurements from the X-SHOOTER spectrograph on the European Southern Observatory's Very Large Telescope (VLT) directly constrain any changes in alpha relative to the value measured on Earth (alpha_0). The weighted mean strength of the electromagnetic force over this redshift range in this location in the universe is da/a = (alpha_z - alpha_0)/alpha_0 = (-2.18 +/- 7.27) X 10^{-5}, i.e. we find no evidence for a temporal change from the 4 new very high redshift measurements. When the 4 new measurements are combined with a large existing sample of lower redshift measurements, a new limit on possible spatial variation of da/a is marginally preferred over a no-variation model at the 3.7 sigma level.

For around 100 years, hydrogen spectral modelling has been based on Voigt profile fitting. The semi-classical Voigt profile is based on a 2-level atom approximation. Whilst the Voigt profile is excellent for many circumstances, the accuracy is insufficient for very high column density damped Lyman-$\alpha$ absorption systems. We have adapted the quantum-mechanical Kramers-Heisenberg model to include thermal broadening, producing a new profile, the KHT profile. Interactions involving multiple discrete atomic levels and continuum terms, not accounted for in the Voigt model, generate asymmetries in the Lyman line wings. If not modelled, this can lead to significant systematics in parameter estimation when modelling real data. There are important ramifications in particular for measurements of the primordial deuterium abundance. However, the KHT model is complicated. We therefore present a simplified formulation based on Taylor series expansions and look-up tables, quantifying the impact of the approximations made. The KHT profile has been implemented within the widely-used VPFIT code.

(abridged) Observing the signature of accretion from the intergalactic medium (IGM) onto galaxies at z~3 requires the detection of faint (L<<L*) galaxies embedded in a filamentary matrix of low-density, metal-poor gas coherent over hundreds of kpc. We study the gaseous environment of three Lyman$\alpha$ emitters (LAEs) at z=2.7-2.8, found to be aligned in projection with a background QSO over ~250 kpc along the slit of a long-slit spectrum. The lack of detection of the LAEs in deep continuum images and the low inferred Ly$\alpha$ luminosities show the LAEs to be intrinsically faint, low-mass galaxies (L<0.1 L*, M_star< 0.1 M*). An echelle spectrum of the QSO reveals strong Ly-alpha absorption within $\pm200$ km/s from the LAEs. Our absorption line analysis leads to HI column densities in the range of log N(HI) =16-18. Associated absorption from ionic metal species CIV and SiIV constrains the gas metallicities to ~0.01 solar if the gas is optically thin, and possibly as low as ~0.001 solar if the gas is optically thick, assuming photoionization equilibrium. While the inferred metallicities are at least a factor of ten lower than expected metallicities in the interstellar medium (ISM) of these LAEs, they are consistent with the observed chemical enrichment level in the IGM at the same epoch. Total metal abundances and kinematic arguments suggest that these faint galaxies have not been able to affect the properties of their surrounding gas. The projected spatial alignment of the LAEs, together with the kinematic quiescence and correspondence between the LAEs and absorbing gas in velocity space suggests that these observations probe a possible filamentary structure. Taken together with the blue-dominant Ly$\alpha$ emission line profile of one of the objects, the evidence suggests that the absorbing gas is part of an accretion stream of low-metallicity gas in the IGM.

White dwarf atmospheres are subjected to gravitational potentials around $10^5$ times larger than occur on Earth. They provide a unique environment in which to search for any possible variation in fundamental physics in the presence of strong gravitational fields. However, a sufficiently strong magnetic field will alter absorption line profiles and introduce additional uncertainties in measurements of the fine structure constant. Estimating the magnetic field strength is thus essential in this context. Here we model the absorption profiles of a large number of atomic transitions in the white dwarf photosphere, including first-order Zeeman effects in the line profiles, varying the magnetic field as a free parameter. We apply the method to a high signal-to-noise, high-resolution, far-ultraviolet HST/STIS spectrum of the white dwarf G191-B2B. The method yields a sensitive upper limit on its magnetic field of $B < 2300$ Gauss at the $3\sigma$ level. Using this upper limit we find that the potential impact of quadratic Zeeman shifts on measurements of the fine structure constant in G191-B2B is 4 orders of magnitude below laboratory wavelength uncertainties.

We present the first data release (DR1) of the UVES Spectral Quasar Absorption Database (SQUAD), comprising 467 fully reduced, continuum-fitted high-resolution quasar spectra from the Ultraviolet and Visual Echelle Spectrograph (UVES) on the European Southern Observatory's Very Large Telescope. The quasars have redshifts $z=0$-5, and a total exposure time of 10 million seconds provides continuum-to-noise ratios of 4-342 (median 20) per 2.5-km/s pixel at 5500 \AA. The SQUAD spectra are fully reproducible from the raw, archival UVES exposures with open-source software, including our UVES_popler tool for combining multiple extracted echelle exposures which we document here. All processing steps are completely transparent and can be improved upon or modified for specific applications. A primary goal of SQUAD is to enable statistical studies of large quasar and absorber samples, and we provide tools and basic information to assist three broad scientific uses: studies of damped Lyman-$\alpha$ systems (DLAs), absorption-line surveys and time-variable absorption lines. For example, we provide a catalogue of 155 DLAs whose Lyman-$\alpha$ lines are covered by the DR1 spectra, 18 of which are reported for the first time. The HI column densities of these new DLAs are measured from the DR1 spectra. DR1 is publicly available and includes all reduced data and information to reproduce the final spectra.

We present two measurements of the temperature-density relationship (TDR) of the intergalactic medium (IGM) in the redshift range $2.55 < z < 2.95$ using a sample of 13 high-quality quasar spectra and high resolution numerical simulations of the IGM. Our approach is based on fitting the neutral hydrogen column density $N_{HI}$ and the Doppler parameter $b$ of the absorption lines in the \mlya\ forest. The first measurement is obtained using a novel Bayesian scheme which takes into account the statistical correlations between the parameters characterising the lower cut-off of the $b-N_{HI}$ distribution and the power-law parameters $T_0$ and $\gamma$ describing the TDR. This approach yields $T_0/ 10^3\, {\rm K}=15.6 \pm 4.4 $ and $\gamma=1.45 \pm 0.17$ independent of the assumed pressure smoothing of the small scale density field. In order to explore the information contained in the overall $b-N_{HI}$ distribution rather than only the lower cut-off, we obtain a second measurement based on a similar Bayesian analysis of the median Doppler parameter for separate column-density ranges of the absorbers. In this case we obtain $T_0/ 10^3\, {\rm K}=14.6 \pm 3.7$ and $\gamma=1.37 \pm 0.17$ in good agreement with the first measurement. Our Bayesian analysis reveals strong anti-correlations between the inferred $T_0$ and $\gamma$ for both methods as well as an anti-correlation of the inferred $T_0$ and the pressure smoothing length for the second method, suggesting that the measurement accuracy can in the latter case be substantially increased if independent constraints on the smoothing are obtained. Our results are in good agreement with other recent measurements of the thermal state of the IGM probing similar (over-)density ranges.

The primordial deuterium abundance probes fundamental physics during the Big Bang Nucleosynthesis and can be used to infer cosmological parameters. Observationally, the abundance can be measured using absorbing clouds along the lines of sight to distant quasars. Observations of the quasar PKS1937--101 contain two absorbers for which the deuterium abundance has previously been determined. Here we focus on the higher redshift one at $z_{abs} = 3.572$. We present new observations with significantly increased signal-to-noise ratio which enable a far more precise and robust measurement of the deuterium to hydrogen column density ratio, resulting in D/H = $2.62\pm0.05\times10^{-5}$. This particular measurement is of interest because it is among the most precise assessments to date and it has been derived from the second lowest column-density absorber [N(HI) $=17.9\mathrm{cm}^{-2}$] that has so-far been utilised for deuterium abundance measurements. The majority of existing high-precision measurements were obtained from considerably higher column density systems [i.e. N(HI) $>19.4\mathrm{cm}^{-2}$]. This bodes well for future observations as low column density systems are more common.

At low densities the standard ionisation history of the intergalactic medium (IGM) predicts a decreasing temperature of the IGM with decreasing density once hydrogen (and helium) reionisation is complete. Heating the high-redshift, low-density IGM above the temperature expected from photo-heating is difficult, and previous claims of high/rising temperatures in low density regions of the Universe based on the probability density function (PDF) of the opacity in Lyman-$\alpha$ forest data at $2<z<4$ have been met with considerable scepticism, particularly since they appear to be in tension with other constraints on the temperature-density relation (TDR). We utilize here an ultra-high signal-to-noise spectrum of the QSO HE0940-1050 and a novel technique to study the low opacity part of the PDF. We show that there is indeed evidence (at 90% confidence level) that a significant volume fraction of the under-dense regions at $z \sim 3$ has temperatures as high or higher than those at densities comparable to the mean and above. We further demonstrate that this conclusion is nevertheless consistent with measurements of a slope of the TDR in over-dense regions that imply a decreasing temperature with decreasing density, as expected if photo-heating of ionised hydrogen is the dominant heating process. We briefly discuss implications of our findings for the need to invoke either spatial temperature fluctuations, as expected during helium reionization, or additional processes that heat a significant volume fraction of the low-density IGM.

The brightest southern quasar above redshift $z=1$, HE 0515$-$4414, with its strong intervening metal absorption-line system at $z_{abs}=1.1508$, provides a unique opportunity to precisely measure or limit relative variations in the fine-structure constant ($\Delta\alpha/\alpha$). A variation of just $\sim$3 parts per million (ppm) would produce detectable velocity shifts between its many strong metal transitions. Using new and archival observations from the Ultraviolet and Visual Echelle Spectrograph (UVES) we obtain an extremely high signal-to-noise ratio spectrum (peaking at S/N $\approx250$ pix$^{-1}$). This provides the most precise measurement of $\Delta\alpha/\alpha$ from a single absorption system to date, $\Delta\alpha/\alpha=-1.42\pm0.55_{\rm stat}\pm0.65_{\rm sys}$ ppm, comparable with the precision from previous, large samples of $\sim$150 absorbers. The largest systematic error in all (but one) previous similar measurements, including the large samples, was long-range distortions in the wavelength calibration. These would add a $\sim$2 ppm systematic error to our measurement and up to $\sim$10 ppm to other measurements using Mg and Fe transitions. However, we corrected the UVES spectra using well-calibrated spectra of the same quasar from the High Accuracy Radial velocity Planet Searcher (HARPS), leaving a residual 0.59 ppm systematic uncertainty, the largest contribution to our total systematic error. A similar approach, using short observations on future, well-calibrated spectrographs to correct existing, high S/N spectra, would efficiently enable a large sample of reliable $\Delta\alpha/\alpha$ measurements. The high S/N UVES spectrum also provides insights into analysis difficulties, detector artifacts and systematic errors likely to arise from 25-40-m telescopes.

In this work, we investigate the abundance and distribution of metals in the intergalactic medium (IGM) at $\langle z \rangle \simeq 2.8$ through the analysis of an ultra-high signal-to-noise ratio UVES spectrum of the quasar HE0940-1050. In the CIV forest, our deep spectrum is sensitive at $3\,\sigma$ to lines with column density down to $\log N_{\rm CIV} \simeq 11.4$ and in 60 per cent of the considered redshift range down to $\simeq11.1$. In our sample, all HI lines with $\log N_{\rm HI} \ge 14.8$ show an associated CIV absorption. In the range $14.0 \le \log N_{\rm HI} <14.8$, 43 per cent of HI lines has an associated CIV absorption. At $\log N_{\rm HI} < 14.0$, the detection rates drop to $<10$ per cent, possibly due to our sensitivity limits and not to an actual variation of the gas abundance properties. In the range $\log N_{\rm HI} \ge 14$, we observe a fraction of HI lines with detected CIV a factor of 2 larger than the fraction of HI lines lying in the circum-galactic medium (CGM) of relatively bright Lyman-break galaxies hosted by dark matter haloes with $\langle M\rangle \sim10^{12}$ M$_{\odot}$. The comparison of our results with the output of a grid of photoionization models and of two cosmological simulations implies that the volume filling factor of the IGM gas enriched to a metallicity $\log Z/Z_{\odot} \ge -3$ should be of the order of $\sim 10-13$ percent. In conclusion, our results favour a scenario in which metals are found also outside the CGM of bright star-forming galaxies, possibly due to pollution by lower mass objects and/or to an early enrichment by the first sources.

This paper presents a self-consistent model for the evolution of gas produced in the debris disc of $\beta$ Pictoris. Our model proposes that atomic carbon and oxygen are created from the photodissociation of CO, which is itself released from volatile-rich bodies in the debris disc due to grain-grain collisions or photodesorption. While the CO lasts less than one orbit, the atomic gas evolves by viscous spreading resulting in an accretion disc inside the parent belt and a decretion disc outside. The temperature, ionisation fraction and population levels of carbon and oxygen are followed with the photodissociation region model Cloudy, which is coupled to a dynamical viscous $\alpha$ model. We present new gas observations of $\beta$ Pic, of C I observed with APEX and O I observed with Herschel, and show that these along with published C II and CO observations can all be explained with this new model. Our model requires a viscosity $\alpha$ > 0.1, similar to that found in sufficiently ionised discs of other astronomical objects; we propose that the magnetorotational instability is at play in this highly ionised and dilute medium. This new model can be tested from its predictions for high resolution ALMA observations of C I. We also constrain the water content of the planetesimals in $\beta$ Pic. The scenario proposed here might be at play in all debris discs and this model could be used more generally on all discs with C, O or CO detections.

From the Voigt profile fitting analysis of 183 intervening CIV systems at 1.7 < z < 3.3 in 23 high-quality UVES/VLT and HIRES/Keck QSO spectra, we find that a majority of CIV systems (~75%) display a well-characterised scaling relation between integrated column densities of HI and CIV with a negligible redshift evolution, when column densities of all the HI and CIV components are integrated within a given (-150, +150) km/sec range centred at the CIV flux minimum. The integrated CIV column density N(CIV, sys) increases with N(HI, sys) at log N(HI, sys) = 14.0--15.5 and log N(CIV, sys) = 11.8--14.0, then becomes almost independent of N(HI, sys) at log N(HI, sys) > 16, with a large scatter: at log N(HI, sys) = 14--22, log N(CIV, sys) = C1 / (log(NHI, sys) + C2) + C3, with C1 = -1.90+0.55, C2 = -14.11+0.19 and C3 = 14.76+0.17, respectively. The steep (flat) part is dominated by SiIV-free (SiIV-enriched) CIV systems. Extrapolating the N(HI, sys)-N(CIV, sys) relation implies that most absorbers with log N(HI) < 14 are virtually CIV-free. The N(HI, sys)-N(CIV, sys) relation does not hold for individual components, clumps or the integration velocity range less than +-100 km/sec. It is expected if the line-of-sight extent of CIV is smaller than HI and N(CIV, sys) decreases more rapidly than N(HI, sys) at the larger impact parameter, regardless of the location of the HI+CIV gas in the IGM filaments or in the intervening galactic halos.

We present the Voigt profile analysis of 132 intervening CIV+CIII components associated with optically-thin HI absorbers at 2.1 < z < 3.4 in the 19 high-quality UVES/VLT and HIRES/Keck QSO spectra. For log N(CIV) = [11.7, 14.1], N(CIII) is proportional to N(CIV) with an exponent (1.42 +- 0.11) and < N(CIII)/N(CIV) > = 1.0 +- 0.3 with a negligible redshift evolution. For 54 CIV components tied (aligned) with HI at log N(HI) = [12.2, 16.0] and log N(CIV) = [11.8, 13.8], the gas temperature T_b estimated from absorption line widths is well-approximated to a Gaussian peaking at log T_b ~ 4.4 +- 0.3 for log T_b = [3.5, 5.5], with a negligible non-thermal contribution. For 32 of 54 tied HI+CIV pairs, also tied with CIII at log N(CIII) = [11.7, 13.8], we ran both photoionisation equilibrium (PIE) and non-PIE (using a fixed temperature T_b) Cloudy models for the Haardt-Madau QSO+galaxy 2012 UV background. We find evidence of bimodality in observed and derived physical properties. High-metallicity branch absorbers have a carbon abundance [C/H]_temp > -1.0, a line-of-sight length L_temp < 20 kpc, and a total (neutral and ionised) hydrogen volume density log n(H, temp) = [-4.5, -3.3] and and log T_b = [3.9, 4.5]. Low-metallicity branch absorbers have [C/H]_temp < -1.0, L_temp = [20, 480] kpc and log n(H, temp) = [-5.2, -4.3] and log T_b ~ 4.5. High-metallicity branch absorbers seem to be originated from extended disks, inner halos or outflowing gas of intervening galaxies, while low-metallicity absorbers are produced by galactic halos or the surrounding IGM filament.

Laming (1990) predicted that the narrow Balmer line core of the ~3000 km/s shock in the SN 1006 remnant would be significantly polarized due to electron and proton impact polarization. Here, based on deep spectrally resolved polarimetry obtained with the European Southern Observatory (ESO)'s Very Large Telescope (VLT), we report the discovery of polarized line emission of polarization degree approx 1.3 percent with position angle orthogonal to the SNR filament. Correcting for an unpolarized broad line component, the implied narrow line polarization is approx 2.0 percent, close to the predictions of Laming (1990). The predicted polarization is primarily sensitive to shock velocity and post-shock temperature equilibration. By measuring polarization for the SN1006 remnant, we validate and enable a new diagnostic that has important applications in a wide variety of astrophysical situations, such as shocks, intense radiation fields, high energy particle streams and conductive interfaces.

The primordial deuterium abundance is an important tracer of the fundamental physics taking place during Big Bang Nucleosynthesis. It can be determined from absorption features along the line of sight to distant quasars. The quasar PKS1937-1009 contains two absorptions systems that have been used to measure the primordial deuterium abundance, the lower redshift one being at z_abs = 3.256. New observations of this absorber are of a substantially higher signal-to-noise and thus permit a significantly more robust estimate of the primordial deuterium abundance, leading to a D/H ratio of 2.45+/-0.28 x10^-5. Whilst the precision of the new measurement presented here is below that obtained from the recent cosmological parameter measurements by Planck, our analysis illustrates how a statistical sample obtained using similarly high spectral signal-to-noise can make deuterium a competitive and complementary cosmological parameter estimator and provide an explanation for the scatter seen between some existing deuterium measurements.

We present the first large, blind and uniformly selected survey for molecular hydrogen (H_2) in damped Lyman alpha systems (DLAs) with moderate-to-high resolution spectra. 86 DLAs were searched for absorption in the many Lyman and Werner H_2 transitions, with approximately 79% completeness for H_2 column densities above N(H_2)=10^17.5 cm^-2 for an assumed Doppler broadening parameter b=2 km s^-1. Only a single strong H_2 absorber was found -- a system detected previously in VLT/UVES spectra. Given our distribution of N(H_2) upper limits, this ~1% detection rate is smaller than expected from previous surveys at 99.8% confidence. Assuming the N(H_2) distribution shape from previous surveys, our detection rate implies a covering factor of ~1% for N(H_2) >= 10^17.5 cm^-2 gas in DLAs (<6% at 95% confidence). We obtained new Magellan/MagE spectra for 53 DLAs; 8 km s^-1 -resolution spectra were available for 27 DLAs. MagE's moderate resolution (~71 km s^-1) yields weaker N(H_2) upper limits and makes them dependent on the assumed Doppler parameter. For example, half the (relevant) previous H_2 detections have N(H_2) >= 10^18.1 cm^-2, a factor of just 3 higher than our median upper limit. Nevertheless, several tests suggest our upper limits are accurate, and they would need to be increased by 1.8 dex to bring our detection rate within 95% confidence of previous surveys.

The nature of the interaction between low-excitation gas filaments at ~10^4 K, seen in optical line emission, and diffuse X-ray emitting coronal gas at ~10^7 K in the centers of galaxy clusters remains a puzzle. The presence of a strong, empirical correlation between the two gas phases is indicative of a fundamental relationship between them, though as yet of undetermined cause. The cooler filaments, originally thought to have condensed from the hot gas, could also arise from a merger or the disturbance of cool circumnuclear gas by nuclear activity. Here, we have searched for intrinsic line emission polarization in cool core galaxy clusters as a diagnostic of fundamental transport processes. Drawing on developments in solar astrophysics, direct energetic particle impact induced polarization holds the promise to definitively determine the role of collisional processes such as thermal conduction in the ISM physics of galaxy clusters, while providing insight into other highly anisotropic excitation mechanisms such as shocks, intense radiation fields and suprathermal particles. Under certain physical conditions, theoretical calculations predict of order ten percent polarization. Our observations of the filaments in four nearby cool core clusters place stringent upper limits (<0.1%) on the presence of emission line polarization, requiring that if thermal conduction is operative, the thermal gradients are not in the saturated regime. This limit is consistent with theoretical models of the thermal structure of filament interfaces.

Measuring rest-frame ultraviolet rotational transitions from the Lyman and Werner bands in absorption against a bright background continuum is one of the few ways to directly measure molecular hydrogen (H2). Here we report the detection of Lyman-Werner absorption from H2 at z=0.56 in a sub-damped Ly-alpha system with neutral hydrogen column density N(HI) = 10^(19.5 +/- 0.2) cm^-2. This is the first H2 system analysed at a redshift < 1.5 beyond the Milky Way halo. It has a surprisingly high molecular fraction: log f(H2) > -1.93 +/- 0.36 based on modelling the line profiles, with a robust model-independent lower limit of f(H2) > 10^-3. This is higher than f(H2) values seen along sightlines with similar N(HI) through the Milky Way disk and the Magellanic clouds. The metallicity of the absorber is 0.19 +0.21 -0.10 solar, with a dust-to-gas ratio < 0.36 times the value in the solar neighbourhood. Absorption from associated low-ionisation metal transitions such as OI and FeII is observed in addition to OVI. Using Cloudy models we show that there are three phases present; a ~100 K phase giving rise to H2, a ~10^4 K phase where most of the low-ionisation metal absorption is produced; and a hotter phase associated with OVI. Based on similarities to high velocity clouds in the Milky Way halo showing H2 and the presence of two nearby galaxy candidates with impact parameters of ~10 kpc, we suggest that the absorber may be produced by a tidally-stripped structure similar to the Magellanic Stream.

We have investigated the distribution and evolution of ~3100 intergalactic HI absorbers with HI column densities log N(HI) = [12.75, 17.0] at 1.9 < z < 3.2, using 18 high resolution, high S/N quasar spectra obtained from the ESO VLT/UVES archive. We used two sets of Voigt profile fitting analysis, one including all the available high-order Lyman lines to obtain reliable HI column densities of saturated lines, and another using only the Ly-alpha lines. There is no significant difference between the results from the two fits. Combining our results with literature data, the mean number density at 0 < z < 4 is not well described by a single power law and strongly suggests that its evolution slows down at z < 1.5 at the high and low column density ranges. We also divided our entire HI absorbers at 1.9 < z < 3.2 into two samples, the unenriched forest and the CIV-enriched forest, depending on whether HI lines are associated with CIV at log N(CIV) > 12.2 within a given velocity range. The entire HI column density distribution function (CDDF) can be described as the combination of these two well-characterised populations which overlap at log N(HI) ~ 15. At log N(HI) < 15, the unenriched forest dominates, showing a similar power-law distribution to the entire forest. The CIV-enriched forest dominates at log N(HI) > 15, with its distribution function proportional to N(HI)^(-1.45). However, it starts to flatten out at lower N(HI), since the enriched forest fraction decreases with decreasing N(HI). The deviation from the power law at log N(HI) = [14, 17] shown in the CDDF for the entire HI sample is a result of combining two different HI populations with a different CDDF shape. The total HI mass density relative to the critical density is Omega(HI) ~ 1.6 x 10^(-6) h^(-1), where the enriched forest accounts for ~40% of Omega(HI).

In this third paper in a series on the nature of extended, asymmetric Lyman alpha emitters at z ~ 3 we report the discovery, in an ultra-deep, blind, spectroscopic long-slit survey, of a Lyman alpha emitting halo around a QSO at redshift 3.045. The QSO is a previously known, obscured AGN. The halo appears extended along the direction of the slit and exhibits two faint patches separated by 17 proper kpc in projection from the QSO. Comparison of the 2-dimensional spectrum with archival HST ACS images shows that these patches coincide spatially with emission from a peculiar, dumbbell-shaped, faint galaxy. The assumptions that the Lyman alpha emission patches are originating in the galaxy and that the galaxy is physically related to the QSO are at variance with photometric estimates of the galaxy redshift. We show, however, that a population of very young stars at the redshift of the QSO may fit the existing rest frame broad band UV photometry of the galaxy. If this scenario is correct then the symmetry of the galaxy in continuum and Lyman alpha emission, the extension of the QSO's Lyman alpha emission in its direction, and the likely presence of a young stellar population in close proximity to a (short-lived) AGN suggest that this may be an example of AGN feedback triggering external star formation in high redshift galaxies.

The physical relationship between low-excitation gas filaments at ~10^4 K, seen in optical line emission, and diffuse X-ray emitting coronal gas at ~10^7 K in the centers of many galaxy clusters is not understood. It is unclear whether the ~10^4 K filaments have cooled and condensed from the ambient hot (~10^7 K) medium or have some other origin such as the infall of cold gas in a merger, or the disturbance of an internal cool reservoir of gas by nuclear activity. Observations of gas at intermediate temperatures (~10^5-10^6 K) can potentially reveal whether the central massive galaxies are gaining cool gas through condensation or losing it through conductive evaporation and hence identify plausible scenarios for transport processes in galaxy cluster gas. Here we present spectroscopic detection of ~10^5 K gas spatially associated with the H-alpha filaments in a central cluster galaxy, M87 in the Virgo Cluster. The measured emission-line fluxes from triply ionized carbon (CIV 1549 A) and singly ionized helium (HeII 1640 A) are consistent with a model in which thermal conduction determines the interaction between hot and cold phases.

(abridged) We present a new analysis of a large sample of quasar absorption-line spectra obtained using UVES (the Ultraviolet and Visual Echelle Spectrograph) on the VLT (Very Large Telescope) in Chile. In the VLT sample (154 absorbers), we find evidence that alpha increases with increasing cosmological distance from Earth. However, as previously shown, the Keck sample (141 absorbers) provided evidence for a smaller alpha in the distant absorption clouds. Upon combining the samples an apparent variation of alpha across the sky emerges which is well represented by an angular dipole model pointing in the direction RA=(17.3 +/- 1.0) hr, dec. = (-61 +/- 10) deg, with amplitude (0.97 +0.22/-0.20) x 10^(-5). The dipole model is required at the 4.1 sigma statistical significance level over a simple monopole model where alpha is the same across the sky (but possibly different to the current laboratory value). The data sets reveal a number of remarkable consistencies: various data cuts are consistent and there is consistency in the overlap region of the Keck and VLT samples. Assuming a dipole-only (i.e. no-monopole) model whose amplitude grows proportionally with `lookback-time distance' (r=ct, where t is the lookback time), the amplitude is (1.1 +/- 0.2) x 10^(-6) GLyr^(-1) and the model is significant at the 4.2 sigma confidence level over the null model [Delta alpha]/alpha = 0). We apply robustness checks and demonstrate that the dipole effect does not originate from a small subset of the absorbers or spectra. We present an analysis of systematic effects, and are unable to identify any single systematic effect which can emulate the observed variation in alpha.

By comparing the widths of absorption lines from OI, SiII and FeII in the redshift z=2.076 single-component damped Lyman alpha absorption system in the spectrum of Q2206-199 we establish that these absorption lines arise in Warm Neutral Medium gas at ~12000 +/- 3000K. This is consistent with thermal equilibrium model estimates of ~ 8000K for the Warm Neutral Medium in galaxies, but not with the presence of a significant cold component. It is also consistent with, but not required by, the absence of CII* fine structure absorption in this system. Some possible implications concerning abundance estimates in narrow-line WNM absorbers are discussed.

We present measurements of carbon, oxygen, silicon, and iron in quasar absorption systems existing when the universe was roughly one billion years old. We measure column densities in nine low-ionization systems at 4.7 < z < 6.3 using Keck, Magellan, and VLT optical and near-infrared spectra with moderate to high resolution. The column density ratios among C II, O I, Si II, and Fe II are nearly identical to sub-DLAs and metal-poor ([M/H] < -1) DLAs at lower redshifts, with no significant evolution over 2 < z < 6. The estimated intrinsic scatter in the ratio of any two elements is also small, with a typical r.m.s. deviation of <0.1 dex. These facts suggest that dust depletion and ionization effects are minimal in our z > 4.7 systems, as in the lower-redshift DLAs, and that the column density ratios are close to the intrinsic relative element abundances. The abundances in our z > 4.7 systems are therefore likely to represent the typical integrated yields from stellar populations within the first gigayear of cosmic history. Due to the time limit imposed by the age of the universe at these redshifts, our measurements thus place direct constraints on the metal production of massive stars, including iron yields of prompt supernovae. The lack of redshift evolution further suggests that the metal inventories of most metal-poor absorption systems at z > 2 are also dominated by massive stars, with minimal contributions from delayed Type Ia supernovae or AGB winds. The relative abundances in our systems broadly agree with those in very metal-poor, non-carbon-enhanced Galactic halo stars. This is consistent with the picture in which present-day metal-poor stars were potentially formed as early as one billion years after the Big Bang.

We previously reported Keck telescope observations suggesting a smaller value of the fine structure constant, alpha, at high redshift. New Very Large Telescope (VLT) data, probing a different direction in the universe, shows an inverse evolution; alpha increases at high redshift. Although the pattern could be due to as yet undetected systematic effects, with the systematics as presently understood the combined dataset fits a spatial dipole, significant at the 4.2-sigma level, in the direction right ascension 17.5 +/- 0.9 hours, declination -58 +/- 9 degrees. The independent VLT and Keck samples give consistent dipole directions and amplitudes, as do high and low redshift samples. A search for systematics, using observations duplicated at both telescopes, reveals none so far which emulate this result.

[ABRIDGED] We examine the velocity structure in the gas associated with \ion{H}{1} in the damped Ly$\alpha$ absorption system at redshift $z=1.7764$ towards the QSO $1331+170$ using 21cm data, optical and STIS spectra. We find at least two, and possibly three, components showing \ion{C}{1} lines. One of these has Doppler parameter $b=0.55${\kms}, corresponding to a kinetic temperature of 220K if the broadening is thermal. We re-examine the H$_2$ analysis undertaken by \citet{Cui05} using the neutral carbon velocity structure, and find a model which is, consistent with a mixture of collisional and background radiation excitation of the observed H$_2$ rotational levels. For singly ionized heavy elements we find eight components covering a velocity range of $\sim 110$ {\kms}. The \ion{H}{1} structure is expected to follow some combination of the singly ionized and neutral gas, but the 21cm absorption profile is considerably different. This may be because of the different extent and brightness distributions of the radio and optical background sources, and so the spin temperature derived by comparing the Ly$\alpha$ and 21cm line strengths has little physical meaning. The neutral and singly ionized heavy element line profiles also show significant differences, and so the dominant components in each appear to be physically distinct. Attempts to use the range of atomic masses to separate thermal and turbulent components of their Doppler widths were not generally successful. The velocity structure in all ionization stages up to $+3$, apart from the neutral heavy elements, is sufficiently complex that it is difficult to separate out the corresponding velocity components for different ionization levels and determine their column densities.

We report two detections of deuterated molecular hydrogen (HD) in QSO absorption-line systems at $z > 2$. Toward J2123-0500, we find $N$(HD) $= 13.84 \pm 0.2$ for a sub-DLA with metallicity $\simeq 0.5Z_{\odot}$ and $N$(H$_2$) = $17.64 \pm 0.15$ at $z = 2.0594$. Toward FJ0812+32, we find $N$(HD) $= 15.38 \pm 0.3$ for a solar-metallicity DLA with $N$(H$_2$) = $19.88 \pm 0.2$ at $z = 2.6265$. These systems have ratios of HD to H$_2$ above that observed in dense clouds within the Milky Way disk and apparently consistent with a simple conversion from the cosmological ratio of D/H. These ratios are not readily explained by any available model of HD chemistry and there are no obvious trends with metallicity or molecular content. Taken together, these two systems and the two published $z > 2$ HD-bearing DLAs indicate that HD is either less effectively dissociated or more efficiently produced in high-redshift interstellar gas, even at low molecular fraction and/or solar metallicity. It is puzzling that such diverse systems should show such consistent HD/H$_2$ ratios. Without clear knowledge of all the aspects of HD chemistry that may help determine the ratio HD/H$_2$, we conclude that these systems are potentially more revealing of gas chemistry than of D/H itself and that it is premature to use such systems to constrain D/H at high-redshift.

We present an analysis of the Ly-a forest toward 3C 273 from the Space Telescope Imaging Spectrograph at ~7 km/s resolution, along with re-processed data from the Far Ultraviolet Spectroscopic Explorer. The high UV flux of 3C 273 allows us to probe the weak, low z absorbers. The main sample consists of 21 HI absorbers that we could discriminate to a sensitivity of log NHI~ 12.5. The redshift density for absorbers with 13.1<log NHI<14.0 is ~1.5 sigma below the mean for other lines of sight; for log NHI >= 12.5, it is consistent with numerical model predictions. The Doppler parameter distribution is consistent with other low z samples. We find no evidence for a break in the column density power-law distribution to log NHI=12.3. A broad Ly-a absorber (BLA) is within Delta v =< 50 km/s and 1.3 local frame Mpc of two ~0.5L* galaxies, with an OVI absorber ~700 km/s away, similarly close to three galaxies and indicating overdense environments. We detect clustering on the Delta v<1000 km/s scale at 3.4 sigma significance for log NHI >= 12.6, consistent with the level predicted from hydrodynamical simulations, and indication for a Ly-a forest void at 0.09<z<0.12. We find at least two components for the z=0.0053 Virgo absorber, but the total NHI column is not significantly changed.

Molecular transitions recently discovered at redshift z_abs=2.059 toward the bright background quasar J2123-0050 are analysed to limit cosmological variation in the proton-to-electron mass ratio, mu=m_p/m_e. Observed with the Keck telescope, the optical spectrum has the highest resolving power and largest number (86) of H_2 transitions in such analyses so far. Also, (7) HD transitions are used for the first time to constrain mu-variation. These factors, and an analysis employing the fewest possible free parameters, strongly constrain mu's relative deviation from the current laboratory value: dmu/mu =(+5.6+/-5.5_stat+/-2.7_sys)x10^{-6}. This is the first Keck result to complement recent constraints from three systems at z_abs>2.5 observed with the Very Large Telescope.

Molecular transitions recently discovered at redshift z_abs=2.059 toward the bright background quasar J2123-0050 are analysed to limit cosmological variation in the proton-to-electron mass ratio, mu=m_p/m_e. Observed with the Keck telescope, the optical echelle spectrum has the highest resolving power and largest number (86) of H_2 transitions in such analyses so far. Also, (seven) HD transitions are used for the first time to constrain mu-variation. These factors, and an analysis employing the fewest possible free parameters, strongly constrain mu's relative deviation from the current laboratory value: dmu/mu =(+5.6+/-5.5_stat+/-2.9_sys)x10^{-6}, indicating an insignificantly larger mu in the absorber. This is the first Keck result to complement recent null constraints from three systems at z_abs>2.5 observed with the Very Large Telescope. The main possible systematic errors stem from wavelength calibration uncertainties. In particular, distortions in the wavelength solution on echelle order scales are estimated to contribute approximately half the total systematic error component, but our estimate is model dependent and may therefore under or overestimate the real effect, if present. To assist future mu-variation analyses of this kind, and other astrophysical studies of H_2 in general, we provide a compilation of the most precise laboratory wavelengths and calculated parameters important for absorption-line work with H_2 transitions redwards of the hydrogen Lyman limit.

Gas at intermediate temperature between the hot X-ray emitting coronal gas in galaxies at the centers of galaxy clusters, and the much cooler optical line emitting filaments, yields information on transport processes and plausible scenarios for the relationship between X-ray cool cores and other galactic phenomena such as mergers or the onset of an active galactic nucleus. Hitherto, detection of intermediate temperature gas has proven elusive. Here, we present FUV imaging of the "low excitation" emission filaments of M87 and show strong evidence for the presence of CIV 1549 A emission which arises in gas at temperature ~10^5K co-located with Halpha+[NII] emission from cooler ~10^4K gas. We infer that the hot and cool phases are in thermal communication, and show that quantitatively the emission strength is consistent with thermal conduction, which in turn may account for many of the observed characteristics of cool core galaxy clusters.

We present the first direct evidence for cold gas in a high redshift DLA galaxy. We measured several multiplets of weak neutral carbon (CI) transitions in order to perform a curve of growth analysis. A delta chi-squared test constrains the best fit Doppler parameter, b = 0.33_{-0.04}^{+0.05} km/s, and logN(CI) = 13.30 +- 0.2 cm^-2. This Doppler parameter constrains the kinetic temperature of the gas to T <= 78 K (T <= 115 K, 2 sigma). We used the associated CI fine structure lines to constrain the volume density of the gas, n(HI) ~ 40 - 200 cm^-3 (2 sigma), resulting in a lower limit on the cloud size of approximately 0.1 - 1 parsec. While it is difficult to determine the metallicity of the cold component, the absence of Cr II indicates that the cold cloud suffers a high level of dust depletion. Additionally, the large amount of Lyman and Werner-band molecular hydrogen absorption (log N(H2)_{total} = 19.88 cm^-2, f_{H_2} >= 0.06) with an excitation temperature of T_{ex} = 46 K as determined by the rotational J = 0 and J = 1 states, is consistent with the presence of cold gas. We propose that this cloud may be gravitationally confined and may represent a transition gas-phase from primarily neutral atomic gas, to a colder, denser molecular phase that will eventually host star formation.

We present a strong constraint on variation of the proton-to-electron mass ratio, mu, over cosmological time scales using molecular hydrogen transitions in optical quasar spectra. Using high quality spectra of quasars Q0405-443, Q0347-383 and Q0528-250, variation in mu relative to the present day value is limited to (Delta mu)/mu=(2.6 +/- 3.0) x 10^(-6). We reduce systematic errors compared to previous works by substantially improving the spectral wavelength calibration method and by fitting absorption profiles to the forest of hydrogen Lyman alpha transitions surrounding each H_2 transition. Our results are consistent with no variation, and inconsistent with a previous approx 4 sigma detection of mu variation involving Q0405-443 and Q0347-383. If the results of this work and those suggesting that alpha may be varying are both correct, then this would tend to disfavour certain grand unification models.

We compare the improved measurement of the Lya forest flux probability distribution at 1.7<z<3.2 presented by Kim et al. (2007) to a large set of hydrodynamical simulations of the Lya forest with different cosmological parameters and thermal histories. The simulations are in good agreement with the observational data if the temperature-density relation for the low density intergalactic medium (IGM), T=T_0 Delta^{gamma-1}, is either close to isothermal or inverted (gamma<1). Our results suggest that the voids in the IGM may be significantly hotter and the thermal state of the low density IGM may be substantially more complex than is usually assumed at these redshifts. We discuss radiative transfer effects which alter the spectral shape of ionising radiation during the epoch of HeII reionisation as a possible physical mechanism for achieving an inverted temperature-density relation at z~3.

CODEX and ESPRESSO are concepts for ultra-stable, high-resolution spectrographs at the E-ELT and VLT, respectively. Both instruments are well motivated by distinct sets of science drivers. However, ESPRESSO will also be a stepping stone towards CODEX both in a scientific as well as in a technical sense. Here we discuss this role of ESPRESSO with respect to one of the most exciting CODEX science cases, i.e. the dynamical determination of the cosmic expansion history.

We have performed an extensive Voigt profile analysis of the neutral hydrogen (HI) and metal absorption present in a sample of 18 high resolution, high signal-to-noise QSO spectra observed with VLT/UVES. We use this analysis to separate the metal contribution from the HI absorption and present an improved measurement of the flux probability distribution function (PDF) due to HI absorption alone at <z> = 2.07, 2.52, and 2.94. The flux PDF is sensitive to the continuum fit in the normalised flux range 0.8 < F < 1.0 and to metal absorption at 0.2 < F < 0.8. Our new measurements of the flux PDF due to HI absorption alone are systematically lower at 0.2 < F < 0.8 by up to 30% compared to the widely used measurement of McDonald et al. (2000), based on a significantly smaller sample of Keck/HIRES data. This discrepancy is probably due to a combination of our improved removal of the metal absorption and cosmic variance, since variations in the flux PDF between different lines-of-sight are large. The HI effective optical depth at 1.7 < z < 4 is best fit with a single power law, (0.0023 +-0.0007) (1+z)^(3.65 +- 0.21), in good agreement with previous measurements from comparable data. As also found previously, the effect of noise on the flux distribution is not significant in high resolution, high signal-to-noise data.

Recent observations suggest the incidence of strong intervening MgII absorption systems along the line-of-sight to gamma ray burst (GRB) afterglows is significantly higher than expected from analogous quasar sightlines. One possible explanation is a geometric effect, arising because MgII absorbers only partially cover the quasar continuum regions, in which case MgII absorbers must be considerably smaller than previous estimates. We investigate the production of abnormal absorption profiles by partial coverage and conclude that the lack of any known anomalous profiles in observed systems, whilst constraining, cannot on its own rule out patchy MgII absorbers. In a separate test, we look for differences in the distribution function of MgII equivalent widths over quasar continuum regions and CIII] emission lines. We show that these anomalies should be observable in any scenario where MgII absorbers are very small, but they are not present in the data. We conclude that models invoking small MgII cloudlets to explain the excess of absorbers seen towards GRBs are ruled out.

We carried out a survey for high-metallicity C IV absorbers at redshift z ~ 2.3 in the spectra of 9 high-quality quasar spectra. Using a novel analysis technique, based on detections of C IV lines and automatically determined upper limits on the column densities of H I, C III, N V, and O VI, we find a large (dN/dz > 7) population of photo-ionized, compact (R ~ 10^2 pc), metal-rich (Z >~ Z_solar) C IV clouds with moderate densities (n_H ~ 10^{-3.5} cm^{-3}), properties that we show are robust with respect to uncertainties in the ionization model. In particular, local sources of ionizing radiation, overabundance of oxygen, departures from ionization equilibrium, and collisional ionization would all imply more compact clouds. The clouds are too small to be self-gravitating and pressure confinement is only consistent under special conditions. We argue that the clouds are, in any case, likely to be short-lived and we demonstrate that this implies that the clouds could easily have been responsible for the transport of all metals that end up in the intergalactic medium (IGM). When the clouds reach pressure equilibrium with the general, photo-ionized IGM, the metals will still be concentrated in small high-metallicity patches, but they will look like ordinary, low-metallicity absorbers. We conclude that intergalactic metals are poorly mixed on small scales and that nearly all of the IGM, and thus the universe, may therefore be of primordial composition.

The low density hydrogen and helium in the IGM probed by QSO absorption lines is sensitive to the amplitude and spectral shape of the metagalactic UV background. We use realistic HI and HeII Ly-alpha forest spectra, constructed from state-of-the-art hydrodynamical simulations of a Lambda-CDM Universe, to confirm the reliability of using line profile fitting techniques to infer the ratio of the metagalactic HI and HeII ionization rates. We further show that the large spatial variations and the anti-correlation with HI absorber density observed in the ratio of the measured HeII to HI column densities can be explained in a model where the HI ionization rate is dominated by the combined UV emission from young star forming galaxies and QSOs and the HeII ionization rate is dominated by emission from QSOs only. In such a model the large fluctuations in the column density ratio are due to the small number of QSOs expected to contribute at any given point to the HeII ionization rate. A significant contribution to UV emission at the HeII photoelectric edge from hot gas in galaxies and galaxy groups would decrease the expected fluctuations in the column density ratio. Consequently, this model appears difficult to reconcile with the large increase in HeII opacity fluctuations towards higher redshift. Our results further strengthen previous suggestions that observed HeII Ly-alpha forest spectra at z~2-3.5 probe the tail end of the reionization of HeII by QSOs.

We present STIS 7km/s resolution data of the Ly-a forest toward PKS0405-123 (z=0.574). We make two samples. The strong sample (60 systems) has column density logNHI>13.3 over 0.002<z<0.423, and the weak sample (44) has logNHI>13.1 at 0.020<z<0.234; 7 absorbers show metals, all with OVI, but often offset in velocity from HI. The strong sample Doppler parameters show <b>=47+-22km/s; for the weak, <b>=44+-21km/s. Line blending and s/n effects likely inflate b. The redshift density dN/dz is consistent with previous, lower resolution results for logNHI>14.0. For 13.1<logNHI<14.0, we find a dN/dz overdensity of ~0.2-0.3 dex at 0.127<z<0.234, which we believe arises from cosmic variance. We find Ly-a clustering over Delta v<250km/s for logNHI>13.3, consistent with a numerical model. There is a void in the strong sample at 0.0320<z<0.0814. We detect radial velocity correlations for Delta v<250km/s between Ly-a absorbers and 39 galaxies in a 10' field. The correlation grows with minimum NHI, peaking for logNHI ~> 13.5-14.0. The strength is similar to that of the galaxy-galaxy correlation for our sample, implying that such Ly-a systems have mass log M/Msun= 11.3+1.0-0.6. Including lower column density systems in the sample shows correlations only out to Delta v <125km/s, as would be expected for smaller density perturbations. Galaxy counts and local HI column density appear correlated, peaking at ~5000km/s. We also give column densities for Galactic species. (Abridged)

The flux power spectra of the Lyman-alpha forest from a sample of 27 QSOs taken with the high resolution echelle spectrograph UVES on VLT are presented. We find a similar fluctuation amplitude at the peak of the ``3D'' flux power spectrum at k ~ 0.03 (km/sec)^(-1) as the study by Croft et al. (2002), in the same redshift range. The amplitude of the flux power spectrum increases with decreasing redshift if corrected for the increase in the mean flux level as expected if the evolution of the flux power spectrum is sensitive to the gravitational growth of matter density fluctuations. This is in agreement with the findings of McDonald et al. (2000) at larger redshift. The logarithmic slope of the "3D" flux power spectrum, P_F(k), at large scales k < 0.03 (km/sec)^(-1), is 1.4 +- 0.3, i.e. 0.3 shallower than that found by Croft et al. (2002) but consistent within the errors.

Spectroscopic observations of distant quasars have resulted in the detection of molecular hydrogen in intervening damped Lyman-alpha absorption clouds (DLAs). We use observations compiled from different experimental groups to show that the molecular hydrogen abundance exhibits a dramatic increase over a cosmological time period corresponding to 13% to 24% of the age of the universe. We also tentatively show that the heavy element abundances in the same gas clouds exhibit a faster and more well-defined cosmological evolution compared to the general DLA population over the same time baseline. We argue that this latter point is unsurprising, because the general DLA population arises in a wide variety of galaxy types and environments, and thus a spans broad range of ISM gas-phases and abundances at the same cosmic time. DLAs exhibiting H2 absorption may therefore circumvent this problem, efficiently identifying a narrower class of objects, and provide a more sensitive probe of cosmological chemical evolution.

We demonstrate that the lyman alpha forest flux power spectrum of ``randomised'' QSO absorption spectra is comparable in shape and amplitude to the flux power spectrum of the original observed spectra. In the randomised spectra a random shift in wavelength has been added to the observed absorption lines as identified and fitted with VPFIT. At 0.03 s/km <k<0.1 s/km the ``3D'' power spectrum of the randomised flux agrees with that of observed spectra within the errors. At larger scales it is still greater than 50 % of that of the observed spectra. At smaller scales the flux power spectrum is dominated by metal lines. Lines of increasing column density contribute to the ``3D'' flux power spectrum at increasingly larger scales. Lines with 13< log (N_HI/cm^{-2})<15 dominate at the peak of the ``3D'' power spectrum while strong absorbers with log (N_HI/cm^{-2})>15 dominate at large scales, k< 0.03 s/km. We further show that a fraction of up to 20 % of the mean flux decrement is contributed by strong absorbers. Analysis of the flux power spectrum which use numerical simulations with too few strong absorption systems calibrated with the observed mean flux are likely to underestimate the inferred rms fluctuation amplitude and the slope of the initial dark matter matter power spectrum. In a combined analysis with other data which constrains the dark matter power spectrum on larger scales this can result in a spurious detection of a running spectral index.

We report a further analysis of the ratio of deuterium to hydrogen (D/H) using HST spectra of the z=0.701 Lyman limit system towards the QSO PG1718+481. Initial analyses of this absorber found it gave a high D/H value, 1.8 - 3.1 \times 10^{-4} (Webb et al. 1998), inconsistent with several higher redshift measurements. It is thus important to critically examine this measurement. By analysing the velocity widths of the DI, HI and metal lines present in this system, Kirkman et al. (2001) report that the additional absorption in the blue wing of the lya line can not be DI, with a confidence level of 98%. Here we present a more detailed analysis, taking into account possible wavelength shifts between the three sets of HST spectra used in the analysis. We find that the constraints on this system are not as strong as those claimed by Kirkman et al. The discrepancy between the parameters of the blue wing absorption and the parameters expected for DI is marginally worse than 1 sigma. Tytler et al.(1999) commented on the first analysis of Webb et al.(1997,1998), reporting the presence of a contaminating lower redshift Lyman limit system, with log[N(HI)] = 16.7 at z=0.602, which biases the N(HI) estimate for the main system. Here we show that this absorber actually has log[N(HI)] < 14.6 and does not impact on the estimate of N(HI) in the system of interest at z = 0.701. The purpose of the present paper is to highlight important aspects of the analysis which were not explored in previous studies, and hence help refine the methods used in future analyses of D/H in quasar spectra.

A new analysis of Keck/HIRES observations of the broad absorption line QSO APM 0827+5255 indicates that a number of intervening C IV absorbers give rise to absorption lines for which the observed optical depths for 1548, 1550 A doublet components are not in the expected 2:1 ratio. To compensate for the effect, a local adjustment of the zero-level is required. We model this effect as coverage of one line of sight to this gravitationally lensed QSO and perform a set of simulations to select a sample of lines for which our model provides an explanation for the effect. We use lines in this sample to obtain estimates for minimum C IV absorber sizes from total coverage and the separations of the lines of sight for a range of lens redshifts, z_l, and two cosmologies. We also obtain best estimates for overall sizes from a statistical `hit and miss' approach. For z_l=0.7 our results set a lower limit to sizes of C IV absorbers of ~0.3 h_72^{-1} kpc (~0.5 h_72^{-1} kpc) for Omega_M=1, Omega_Lambda=0 (Omega_M=0.3, Omega_Lambda=0.7), in agreement with other results from similar work but are limited by sample size and the uncertainty in z_l. Our method can be used to detect lensed QSOs and to probe absorber sizes when separate spectra cannot be obtained for each line of sight.

We have performed hydrodynamical simulations to investigate the effects of galactic winds on the high-redshift (z=3) universe. Strong winds suppress the formation of low-mass galaxies significantly, and the metals carried by them produce C IV absorption lines with properties in reasonable agreement with observations. The winds have little effect on the statistics of the H I absorption lines, because the hot gas bubbles blown by the winds fill only a small fraction of the volume and because they tend to escape into the voids, thereby leaving the filaments that produce these lines intact.

A search for O VI at redshifts corresponding to Ly-alpha lines in the z_em ~ 2.4 QSOs HE1122-1648 and HE2217-2818 reveals that a substantial fraction of those with H I column densities log N(HI) > 14 (cm^{-2}) are highly ionized and show some heavy element enrichment. If these two sight lines are typical, then the O VI systems contain a cosmologically significant fraction of the baryons and the metals in the universe. For most systems the temperatures derived from the line widths are too low for collisional ionization to be responsible for the O VI lines. Photoionization models with a substantial hard ultraviolet flux can reproduce the observations for densities that are in good agreement with a model assuming local, hydrostatic equilibrium and heavy element abundances in the range ~ 10^{-3} - 10^{-2} solar. Photoionization by a UV flux much softer than that predicted by Haardt & Madau (1996) for a background dominated by quasars can be ruled out. Finally, we find one system with a very low H I column density for which both photoionization and collisional ionization models yield a metallicity close to solar and a density that is inconsistent with gravitational confinement, unless the gas fraction is negligible.

Combining a new, increased dataset of 8 QSOs covering the Ly-alpha forest at redshifts 1.5 < z < 3.6 from VLT/UVES observations with previously published results, we have investigated the properties of the Ly-alpha forest at 1.5 < z < 4. We have fitted Voigt profiles to the Ly-alpha absorption lines as in previous studies, and have, for two QSOs, fitted Ly-alpha at z ~ 2.2 and higher order of Lyman lines down to 3050 Angstrom simultaneously. For a given HI column density range, the Ly-alpha forest at 1.5 < z < 4 shows the monotonic evolution, in terms of the line number density and the HI opacity. For higher column density clouds at 10^{14} cm^(-2),there is a variation in the line number density from sightline to sightline at z < 2.5. This variation is stronger for higher column density systems, probably due to more gravitationally evolved structures at lower z. The correlation strength of the step optical depth correlation function shows the strong evolution from z ~ 3.3 to z ~ 2.1, although there is a large scatter along different sightlines. The analyses of the Ly-beta forest at z ~ 2.2 are, in general, in good agreement with those of the Ly-alpha forest.

The temperature of the low-density intergalactic medium is set by the balance between adiabatic cooling resulting from the expansion of the universe, and photo-heating by the UV-background. We have analysed the Lyman-alpha forest of eleven high-resolution quasar spectra using wavelets, and find strong evidence of a marked jump in the temperature at the mean density, T_0, of 60 per cent around a redshift z=3.3, which we attribute to reionization of HeliumII. The jump can be seen in all three of our spectra that straddle redshift 3.3, at a significance of 99 per cent. Below z=3.1, our results are consistent with a smooth cooling down of the universe, as expected when adiabatic expansion dominates over photo-heating by a UV-background from QSOs and galaxies. We find no evidence of thermal fluctuations on scales greater than 5000 km/s larger than 50 per cent, which could be detected by our method, suggesting that the IGM follows a reasonably well-defined temperature-density relation. We demonstrate that the mean wavelet amplitude scales inversely with T_0 and calibrate the relation with hydrodynamical simulations. We find T_0= 12000K at z>3.6. Such high temperature suggest that Hydrogen reionization occured relatively recent.

Spectra from the Keck HIRES instrument of the Lyman alpha forests in the lines of sight to the A and C components of the gravitationally lensed QSO Q1422+231 were used to investigate the structure of the intergalactic medium at mean redshift z = 3.3 on sub-kpc scales. We measured the cross-correlation amplitude between the two Lyman alpha forests for a mean transverse separation of 120 pc, and computed the RMS column density and velocity differences between individual absorption systems seen in both lines of sight. The RMS differences between the velocity centroids of the Lyman alpha forest lines were found to be less than about 400 m/s, for unsaturated absorption lines with column densities in the range 12<log N(HI)<14.13. The rate of energy transfer into the low density IGM on a typical scale of 100 pc seems to be lower by 3-4 orders of magnitude than the rate measured earlier for strong CIV metal absorption systems. The tight correlation between HI column density and baryonic density in the intergalactic medium was used to obtain an upper limit on the RMS fluctuations of the baryonic density field. The fraction of the absorption lines that are different across the lines of sight was used to determine the filling factor of the universe for gas which has suffered recent hydrodynamic disturbances. We thereby derived upper limits on the filling factor of galactic outflows at high redshift. Short-lived, short-range ancient winds are essentially unconstrained by this method but strong winds blowing for a substantial fraction of a Hubble time (at z = 3.3) appear to fill less than 20 percent of the volume of the universe.

The z ~ 0.701 absorption system towards the QSO PG1718+4807 is the only example of a QSO absorption system which might have a deuterium/hydrogen ratio approximately ten times the value found towards other QSOs. We have obtained new STIS spectra from the Hubble Space Telescope of the Lyman alpha and Lyman limit regions of the system. These spectra give the redshift and velocity dispersion of the neutral hydrogen which produces most of the observed absorption. The Lyman alpha line is too narrow to account for all of the observed absorption. It was previously known that extra absorption is needed on the blue side of the main H I near to the expected position of deuterium. The current data suggests with a 98% confidence level that the extra absorption is not deuterium. Some uncertainty persists because we have a low signal to noise ratio and the extra absorption - be it deuterium or hydrogen - is heavily blended with the Lyman alpha absorption from the main hydrogen component.

Recent numerical simulations have lead to a paradigm shift in our understanding of the intergalactic medium, and the loss of a physical justification for Voigt profile fitting of the Lyman-alpha forest. Many individual lines seen in simulated spectra have significant departures from the Voigt profile, yet could be well fitted by a blend of two or more such lines. We discuss the expected effect on the line profiles due to ongoing gravitational structure formation and Hubble expansion. We develop a method to detect departures from Voigt profiles of the absorption lines in a statistical way and apply this method to simulated Lyman-alpha forest spectra, confirming that the profiles seen do statistically differ from Voigt profiles.

The spatial distribution of Ly-alpha forest absorption systems towards ten QSOs has been analysed to search for large-scale structure over the redshift range 2.2 < z < 3.4. The QSOs form a closely spaced group on the sky and are concentrated within a 1 deg^2 field. We have employed a technique based on the first and second moments of the transmission probability density function which is capable of identifying and assessing the significance of regions of over- or underdense Ly-alpha absorption. We find evidence for large-scale structure in the distribution of Ly-alpha forest absorption at the > 99 per cent confidence level. In individual spectra we find overdense Ly-alpha absorption on scales of up to 1200 km s^-1. There is also strong evidence for correlated absorption across line of sight pairs separated by < 3 h^-1 proper Mpc (q_0 = 0.5). For larger separations the cross-correlation signal becomes progressively less significant.

(abridged) We present observations of the high-redshift QSO GB1759+7539 (z=3.05) obtained with HIRES on the Keck 10m telescope. The spectrum has a resolution of FWHM = 7 km/s, and a typical signal-to-noise ratio per 2 km/s pixel of ~25 in the Ly-alpha forest region, and ~60 longward of the Ly-alpha emission. The observed Ly-alpha forest systems have a mean redshift of z=2.7. There is marginal evidence of clustering along the line of sight over the velocity range 100 < v < 250 km/s. The 1-point and 2-point joint probability distributions of the transmitted flux for the Ly-alpha forest were calculated, and shown to be very insensitive to the heavy element contamination. We could find no evidence of Voigt profile departures due to infalling gas, as observed in the simulated forest spectra. Twelve heavy-element absorption systems were identified, including damped Lyman-alpha systems at z=2.62 and 2.91. The C, N, O, Al, Si, P, S, Mg, Fe, and Ni absorption features of these systems were studied, and the elemental abundances calculated for the weak unsaturated lines. The systems have metallicities of Z(2.62) ~1/20 Z(solar) and Z(2.91) \~1/45 Z(solar). Both systems appear to have a low dust content. They show an over-abundance of alpha-elements relative to Fe-peak elements, and an under-abundance of odd atomic number elements relative to even. Nitrogen was observed, and found to be under-abundant relative to oxygen, in line with the time delay model of primary nitrogen production. C II* was also seen, allowing us to determine an upper limit to the cosmic microwave background temperature at z=2.62 of T(CMB)<12.9K.

We report on high resolution observations (R=35000) of the Hubble Deep Field South QSO J2233-606 obtained with the University College London Echelle Spectrograph at the AAT. We present spectral data and an associated absorption line list for the wavelength region 3530 < lambda < 4390A. The data has a mean signal-to-noise ratio in the continuum of approximately 15 per 0.05A resolution element at lambda > 3700A.

We present a new HST/STIS spectrum of the z=3.18 quasar PKS 1935-692 and summarize the spectral features shortwards of 304A in the rest frame likely to be caused by foreground HeII Lyman-alpha absorption. In accord with previous results on two other quasars at similar redshifts, we demonstrate a correlation with the HI Lyman-alpha forest absorption, and show that much of the helium absorption is caused by a comparable quantity of more diffuse gas with Omega~0.01, that is not detected in HI. The helium ionization zone around the quasar is detected as well as a void seen in both HI and HeII. The properties of the absorption are in broad agreement with those of the other quasars and with models of the protogalactic gas distribution and ionization at this redshift.

We present a new technique for detecting structure on Mpc scales in the Lyman-alpha forest. The technique is easy to apply in practice since it does not involve absorption line fitting but is rather based on the statistics of the transmitted flux. It identifies and assesses the statistical significance of regions of over- or underdense Lyman-alpha absorption and is fairly insensitive to the quality of the data. Using extensive simulations we demonstrate that the new method is significantly more sensitive to the detection of large-scale structure in the Lyman-alpha forest than a traditional two-point correlation function analysis of fitted absorption lines.

Multidimensional cosmologies allow for variations of fundamental physical constants over the course of cosmological evolution, and different versions of the theories predict different time dependences. In particular, such variations could manifest themselves as changes of the proton-to-electron mass ratio \mu=m_p/m_e over the period of ~ 10^{10} years since the moment of formation of high-redshift QSO spectra. Here we analyze a new, high-resolution spectrum of the z=2.81080 molecular hydrogen absorption system toward the quasar PKS 0528-250 to derive a new observational constraint to the time-averaged variation rate of the proton-to-electron mass ratio. We find |\dot{\mu} / \mu| < 1.5 \times 10^{-14}/year, which is much tighter than previously measured limits.

Recent HST/GHRS observations of the z = 0.7010 absorber toward the QSO 1718+4807 (selected because of apparently ideal characteristics for measuring D/H) yield a D/H significantly higher D/H than some recent high-redshift measurements. Our analysis indicates D/H}= 1.8 - 2.5 x 10^{-4}. This may indicate a cosmological inhomogeneity in the deuterium abundance of at least a factor of ten.

Of the light elements, the primordial abundance of deuterium, (D/H)_p, provides the most sensitive diagnostic for the cosmological mass density parameter Omega_B. Recent high redshift (D/H) measurements are highly discrepant, although this may reflect observational uncertainties. The larger (D/H) values, which imply a low Omega_B and require the Universe to be dominated by non-baryonic matter (dynamical studies indicate a higher total density parameter), cause problems for galactic chemical evolution models since they have difficulty in reproducing the large decline down to the lower present-day (D/H). Conversely, low (D/H) values imply an Omega_B greater than derived from ^7Li and ^4He abundance measurements, and may require a deuterium abundance evolution that is too low to easily explain. Here we report the first measurement at intermediate redshift, where the observational difficulties are smaller, of a gas cloud with ideal characteristics for this experiment. Our analysis of the z = 0.7010 absorber toward 1718+4807 indicates (D/H) = 2.0 +/- 0.5 x 10^{-4} which is in the high range. This and other independent observations suggests there may be a cosmological inhomogeneity in (D/H)_p of at least a factor of ten.

We use new observations of very weak CIV absorption lines associated with high-redshift Lyman-alpha absorption systems to measure the high-redshift Lyman-alpha line two-point correlation function (TPCF). These very weak CIV absorption lines trace small-scale velocity structure that cannot be resolved by Lyman-alpha absorption lines. We find that (1) high-redshift Lyman-alpha absorption systems with N(HI) > 3.10^14 cm^{-2} are strongly clustered in redshift, (2) previous measurements of the Lyman-alpha line TPCF underestimated the actual clustering of the absorbers due to unresolved blending of overlapping velocity components, (3) the present observations are consistent with the hypothesis that clustering of Lyman-alpha absorption systems extends to lower column densities, but maybe with smaller amplitude in the correlation function, and (4) the observed clustering is broadly compatible with that expected for galaxies at z \sim 2-3. We interpret these results as suggesting that many or most Lyman-alpha absorbers may arise in galaxies even at high redshifts, and, therefore, that the Lyman-alpha forest probes processes of galaxy formation and evolution for redshifts \lesssim 5.

We present results from 12 km/s resolution echelle spectroscopy of the bright $z=4.694$ QSO BR1202-0725. A preliminary analysis shows that high metallicity narrow line absorption clouds are present up to the redshift of the quasar. A damped Ly$\alpha$ system with an HI column density of $3.1 \times 10^{20}$ cm$^{-2}$ at $z=4.383$ has an [O/H] ratio that is about 0.01 solar, while another absorption system with an HI column density of $5.0 \times 10^{16}$ cm$^{-2}$ at $z=4.672$, may have an O/H ratio that is twice solar. An upper limit, or the possible detection of deuterium in this absorption cloud gives a D/H ratio of about 1.5\,10$^{-4}$. Because this cloud is metal rich at least some of the cloud gas has been processed after the beginning of the Universe.