Showing 1–50 of 68 results for author: Carswell, R F

Varying alpha, blinding, and bias in existing measurements

Authors: Chung-Chi Lee, John K. Webb, Robert F. Carswell, Vladimir A. Dzuba, Victor V. Flambaum, Dinko Milaković

Abstract: The high resolution spectrograph ESPRESSO on the VLT allows measurements of fundamental constants at unprecedented precision and hence enables tests for spacetime variations predicted by some theories. In a series of recent papers, we developed optimal analysis procedures that both exposes and eliminates the subjectivity and bias in previous quasar absorption system measurements. In this paper we… ▽ More The high resolution spectrograph ESPRESSO on the VLT allows measurements of fundamental constants at unprecedented precision and hence enables tests for spacetime variations predicted by some theories. In a series of recent papers, we developed optimal analysis procedures that both exposes and eliminates the subjectivity and bias in previous quasar absorption system measurements. In this paper we analyse the ESPRESSO spectrum of the absorption system at z_{abs}=1.15 towards the quasar HE0515-4414. Our goal here is not to provide a new unbiased measurement of fine structure constant, alpha, in this system (that will be done separately). Rather, it is to carefully examine the impact of blinding procedures applied in the recent analysis of the same data by Murphy (2022) and prior to that, in several other analyses. To do this we use supercomputer Monte Carlo AI calculations to generate a large number of independently constructed models of the absorption complex. Each model is obtained using AI-VPFIT, with alpha fixed until a "final" model is obtained, at which point alpha is then released as a free parameter for one final optimisation. The results show that the "measured" value of alpha is systematically biased towards the initially-fixed value i.e. this process produces meaningless measurements. The implication is straightforward: to avoid bias, all future measurements must include alpha as a free parameter from the beginning of the modelling process. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Comments: 17 pages, 7 figures, and 5 tables. Submitted to MNRAS 1 Dec 2022

Addendum: Precision in high resolution absorption line modelling, analytic Voigt derivatives, and optimisation methods

Authors: Chung-Chi Lee, John K. Webb, Robert F. Carswell

Abstract: The parent paper to this Addendum describes the optimisation theory on which VPFIT, a non-linear least-squares program for modelling absorption spectra, is based. In that paper, we show that Voigt function derivatives can be calculated analytically using Taylor series expansions and look-up tables, for the specific case of one column density parameter for each absorption component. However, in man… ▽ More The parent paper to this Addendum describes the optimisation theory on which VPFIT, a non-linear least-squares program for modelling absorption spectra, is based. In that paper, we show that Voigt function derivatives can be calculated analytically using Taylor series expansions and look-up tables, for the specific case of one column density parameter for each absorption component. However, in many situations, modelling requires more complex parameterisation, such as summed column densities over a whole absorption complex, or common pattern relative ion abundances. This Addendum provides those analytic derivatives. △ Less

Submitted 29 December, 2021; originally announced December 2021.

Comments: 4 pages, 2 figures. Submitted to MNRAS 23 Dec 2021, accepted 24 Dec 2021

Precision in high resolution absorption line modelling, analytic Voigt derivatives, and optimisation methods

Authors: John K. Webb, Robert F. Carswell, Chung-Chi Lee

Abstract: This paper describes the optimisation theory on which VPFIT, a non-linear least-squares program for modelling absorption spectra, is based. Particular attention is paid to precision. Voigt function derivatives have previously been calculated using numerical finite difference approximations. We show how these can instead be computed analytically using Taylor series expansions and look-up tables. We… ▽ More This paper describes the optimisation theory on which VPFIT, a non-linear least-squares program for modelling absorption spectra, is based. Particular attention is paid to precision. Voigt function derivatives have previously been calculated using numerical finite difference approximations. We show how these can instead be computed analytically using Taylor series expansions and look-up tables. We introduce a new optimisation method for an efficient descent path to the best-fit, combining the principles used in both the Gauss-Newton and Levenberg-Marquardt algorithms. A simple practical fix for ill-conditioning is described, a common problem when modelling quasar absorption systems. We also summarise how unbiased modelling depends on using an appropriate information criterion to guard against over- or under-fitting. The methods and the new implementations introduced in this paper are aimed at optimal usage of future data from facilities such as ESPRESSO/VLT and HIRES/ELT, particularly for the most demanding applications such as searches for spacetime variations in fundamental constants and attempts to detect cosmological redshift drift. △ Less

Submitted 4 October, 2021; v1 submitted 25 August, 2021; originally announced August 2021.

Comments: 15 pages, 7 figures, to appear in MNRAS

Non-uniqueness in quasar absorption models and implications for measurements of the fine structure constant

Authors: Chung-Chi Lee, John K. Webb, Dinko Milaković, Robert F. Carswell

Abstract: High resolution spectra of quasar absorption systems provide the best constraints on temporal or spatial changes of fundamental constants in the early universe. An important systematic that has never before been quantified concerns model non-uniqueness. The absorption structure is generally complicated, comprising many blended lines. This characteristic means any given system can be fitted equally… ▽ More High resolution spectra of quasar absorption systems provide the best constraints on temporal or spatial changes of fundamental constants in the early universe. An important systematic that has never before been quantified concerns model non-uniqueness. The absorption structure is generally complicated, comprising many blended lines. This characteristic means any given system can be fitted equally well by many slightly different models, each having a different value of α, the fine structure constant. We use AI Monte Carlo modelling to quantify non-uniqueness. Extensive supercomputer calculations are reported, revealing new systematic effects that guide future analyses: (i) Whilst higher signal to noise and improved spectral resolution produces a smaller statistical uncertainty for α, model non-uniqueness adds a significant additional uncertainty. (ii) Non-uniqueness depends on the line broadening mechanism used. We show that modelling the spectral data using turbulent line broadening results in far greater non-uniqueness, hence this should no longer be done. Instead, for varying αstudies, it is important to use the more physically appropriate compound broadening. (iii) We have studied two absorption systems in detail. Generalising thus requires caution. Nevertheless, if non-uniqueness is present in all or most quasar absorption systems, it seems unavoidable that attempts to determine the existence (or non-existence) of spacetime variations of fundamental constants is best approached using a statistical sample. △ Less

Submitted 1 July, 2021; v1 submitted 23 February, 2021; originally announced February 2021.

Comments: 17 pages, 11 figures

The evolution of the low-density HI intergalactic medium from z=3.6 to 0: Data, transmitted flux and HI column density

Authors: Tae-Sun Kim, Bart P. Wakker, Fahad Nasir, Robert F. Carswell, Blair D. Savage, James S. Bolton, Andrew J. Fox, Matteo Viel, Martin G. Haehnelt, Jane C. Charlton, Ben E. Rosenwasser

Abstract: We present a new, uniform analysis of the HI transmitted flux (F) and HI column density (N(HI)) distribution in the low-density IGM as a function of redshift z for 0<z<3.6 using 55 HST/COS FUV (Delta(z)=7.2 at z<0.5), five HST/STIS+COS NUV (Delta(z)=1.3 at z~1) and 24 VLT/UVES and Keck/HIRES (Delta(z)=11.6 at 1.7<z<3.6) AGN spectra. We performed a consistent, uniform Voigt profile analysis to comb… ▽ More We present a new, uniform analysis of the HI transmitted flux (F) and HI column density (N(HI)) distribution in the low-density IGM as a function of redshift z for 0<z<3.6 using 55 HST/COS FUV (Delta(z)=7.2 at z<0.5), five HST/STIS+COS NUV (Delta(z)=1.3 at z~1) and 24 VLT/UVES and Keck/HIRES (Delta(z)=11.6 at 1.7<z<3.6) AGN spectra. We performed a consistent, uniform Voigt profile analysis to combine spectra taken with different instruments, to reduce systematics and to remove metal-line contamination. We confirm previously known conclusions on firmer quantitative grounds in particular by improving the measurements at z~1. Two flux statistics at 0<F<1, the mean HI flux and the flux probability distribution function (PDF), show that considerable evolution occurs from z=3.6 to z=1.5, after which it slows down to become effectively stable for z<0.5. However, there are large sightline variations. For the HI column density distribution function (CDDF, f proptional to N(HI)^(-beta)) at log (N(HI)/1cm^-2)=[13.5, 16.0], beta increases as z decreases from beta~1.60 at z~3.4 to beta~1.82 at z~0.1. The CDDF shape at lower redshifts can be reproduced by a small amount of clockwise rotation of a higher-z CDDF with a slightly larger CDDF normalisation. The absorption line number per z (dn/dz) shows a similar evolutionary break at z~1.5 as seen in the flux statistics. High-N(HI) absorbers evolve more rapidly than low-N(HI) absorbers to decrease in number or cross-section with time. The individual dn/dz shows a large scatter at a given z. The scatter increases toward lower z, possibly caused by a stronger clustering at lower z. △ Less

Submitted 10 December, 2020; originally announced December 2020.

Comments: 24 pages, 17 figures, 8 tables, accepted for publication to MNRAS, 5 supplementary online-only tables

Getting the model right; an information criterion for spectroscopy

Authors: John K. Webb, Chung-Chi Lee, Robert F. Carswell, Dinko Milaković

Abstract: Robust model-fitting to spectroscopic transitions is a requirement across many fields of science. The corrected Akaike and Bayesian information criteria (AICc and BIC) are most frequently used to select the optimal number of fitting parameters. In general, AICc modelling is thought to overfit (too many model parameters) and BIC underfits. For spectroscopic modelling, both AICc and BIC lack in two… ▽ More Robust model-fitting to spectroscopic transitions is a requirement across many fields of science. The corrected Akaike and Bayesian information criteria (AICc and BIC) are most frequently used to select the optimal number of fitting parameters. In general, AICc modelling is thought to overfit (too many model parameters) and BIC underfits. For spectroscopic modelling, both AICc and BIC lack in two important respects: (a) no penalty distinction is made according to line strength such that parameters of weak lines close to the detection threshold are treated with equal importance as strong lines and (b) no account is taken of the way in which spectral lines impact on narrow data regions. In this paper we introduce a new information criterion that addresses these shortcomings, the "Spectral Information Criterion" (SpIC). Spectral simulations are used to compare performances. The main findings are (i) SpIC clearly outperforms AICc for high signal to noise data, (ii) SpIC and AICc work equally well for lower signal to noise data, although SpIC achieves this with fewer parameters, and (iii) BIC does not perform well (for this application) and should be avoided. The new method should be of broader applicability (beyond spectroscopy), wherever different model parameters influence separated small ranges within a larger dataset and/or have widely varying sensitivities. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: 13 pages, 7 figures

A new era of fine structure constant measurements at high redshift

Authors: Dinko Milaković, Chung-Chi Lee, Robert F. Carswell, John K. Webb, Paolo Molaro, Luca Pasquini

Abstract: New observations of the quasar HE0515$-$4414 have been made using the HARPS spectrograph on the ESO 3.6m telescope, aided by the Laser Frequency Comb (LFC). We present three important advances for $α$ measurements in quasar absorption spectra from these observations. Firstly, the data have been wavelength calibrated using LFC and ThAr methods. The LFC wavelength calibration residuals are six times… ▽ More New observations of the quasar HE0515$-$4414 have been made using the HARPS spectrograph on the ESO 3.6m telescope, aided by the Laser Frequency Comb (LFC). We present three important advances for $α$ measurements in quasar absorption spectra from these observations. Firstly, the data have been wavelength calibrated using LFC and ThAr methods. The LFC wavelength calibration residuals are six times smaller than when using the standard ThAr calibration. We give a direct comparison between $α$ measurements made using the two methods. Secondly, spectral modelling was performed using Artificial Intelligence (fully automated, all human bias eliminated), including a temperature parameter for each absorption component. Thirdly, in contrast to previous work, additional model parameters were assigned to measure $α$ for each individual absorption component. The increase in statistical uncertainty from the larger number of model parameters is small and the method allows a substantial advantage; outliers that would otherwise contribute a significant systematic, possibly corrupting the entire measurement, are identified and removed, permitting a more robust overall result. The $z_{abs} = 1.15$ absorption system along the HE0515$-$4414 sightline yields 40 new $α$ measurements. We constrain spatial fluctuations in $α$ to be $Δα/α\leq 9 \times 10^{-5}$ on scales $\approx 20\;{\rm km\,s}^{-1}$, corresponding to $\approx25\;$kpc if the $z_{abs} = 1.15$ system arises in a $1\;$Mpc cluster. Collectively, the 40 measurements yield $Δα/α=-0.27\pm2.41\times10^{-6}$, consistent with no variation. △ Less

Submitted 5 March, 2021; v1 submitted 24 August, 2020; originally announced August 2020.

Comments: Published in MNRAS, 500 1 1. 10 pages, 7 figures. Online supplementary material is available upon request

Artificial intelligence and quasar absorption system modelling; application to fundamental constants at high redshift

Authors: Chung-Chi Lee, John K. Webb, R. F. Carswell, Dinko Milaković

Abstract: Exploring the possibility that fundamental constants of Nature might vary temporally or spatially constitutes one of the key science drivers for the European Southern Observatory's ESPRESSO spectrograph on the VLT and for the HIRES spectrograph on the ELT. High-resolution spectra of quasar absorption systems permit accurate measurements of fundamental constants out to high redshifts. The quality o… ▽ More Exploring the possibility that fundamental constants of Nature might vary temporally or spatially constitutes one of the key science drivers for the European Southern Observatory's ESPRESSO spectrograph on the VLT and for the HIRES spectrograph on the ELT. High-resolution spectra of quasar absorption systems permit accurate measurements of fundamental constants out to high redshifts. The quality of new data demands completely objective and reproducible methods. We have developed a new fully automated Artificial Intelligence-based method capable of deriving optimal models of even the most complex absorption systems known. The AI structure is built around VPFIT, a well-developed and extensively-tested non-linear least-squares code. The new method forms a sophisticated parallelised system, eliminating human decision-making and hence bias. Here we describe the workings of such a system and apply it to synthetic spectra, in doing so establishing methods of importance for future analyses of VLT and ELT data. The results show that modelling line broadening for high-redshift absorption components should include both thermal and turbulent components. Failing to do so means it is easy to derive the wrong model and hence incorrect parameter estimates. We also argue that model non-uniqueness can be significant, such that it is not feasible to expect to derive an unambiguous estimate of the fine structure constant alpha from one or a small number of measurements. No matter how optimal the modelling method, it is a fundamental requirement to use a large sample of measurements to meaningfully constrain temporal or spatial alpha variation. △ Less

Submitted 2 April, 2021; v1 submitted 6 August, 2020; originally announced August 2020.

Comments: Accepted by MNRAS, 1 April 2021

Measuring the fine structure constant on a white dwarf surface; a detailed analysis of Fe V absorption in G191-B2B

Authors: J. Hu, J. K. Webb, T. R. Ayres, M. B. Bainbridge, J. D. Barrow, M. A. Barstow, J. C. Berengut, R. F. Carswell, V. Dumont, V. Dzuba, V. V. Flambaum, C. C. Lee, N. Reindl, S. P. Preval, W. -Ü. L. Tchang-Brillet

Abstract: The gravitational potential phi = GM/Rc^2 at the surface of the white dwarf G191-B2B is 10,000 times stronger than that at the Earth's surface. Numerous photospheric absorption features are detected, making this a suitable environment to test theories in which the fundamental constants depend on gravity. We have measured the fine structure constant, alpha, at the white dwarf surface, used a newly… ▽ More The gravitational potential phi = GM/Rc^2 at the surface of the white dwarf G191-B2B is 10,000 times stronger than that at the Earth's surface. Numerous photospheric absorption features are detected, making this a suitable environment to test theories in which the fundamental constants depend on gravity. We have measured the fine structure constant, alpha, at the white dwarf surface, used a newly calibrated Hubble Space Telescope STIS spectrum of G191-B2B, two new independent sets of laboratory Fe V wavelengths, and new atomic calculations of the sensitivity parameters that quantify Fe V wavelength dependency on alpha. The two results obtained are: dalpha/alpha = 6.36 +/- [0.33(stat) + 1.94(sys)] X 10^{-5} and dalpha/alpha = 4.21 +/- [0.47(stat) + 2.35(sys)] X 10^{-5}. The measurements hint that the fine structure constant increases slightly in the presence of strong gravitational fields. A comprehensive search for systematic errors is summarised, including possible effects from line misidentifications, line blending, stratification of the white dwarf atmosphere, the quadratic Zeeman effect and electric field effects, photospheric velocity flows, long-range wavelength distortions in the HST spectrum, and variations in the relative Fe isotopic abundances. None fully account for the observed deviation but the systematic uncertainties are heavily dominated by laboratory wavelength measurement precision. △ Less

Submitted 27 August, 2020; v1 submitted 21 July, 2020; originally announced July 2020.

Comments: Submitted to MNRAS 21 July 2020. Resubmitted after refereeing 27 August 2020

Four direct measurements of the fine-structure constant 13 billion years ago

Authors: Michael R. Wilczynska, John K. Webb, Matthew Bainbridge, Sarah E. I. Bosman, John D. Barrow, Robert F. Carswell, Mariusz P. Dabrowski, Vincent Dumont, Ana Catarina Leite, Chung-Chi Lee, Katarzyna Leszczynska, Jochen Liske, Konrad Marosek, Carlos J. A. P. Martins, Dinko Milakovic, Paolo Molaro, Luca Pasquini

Abstract: 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… ▽ More 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. △ Less

Submitted 17 March, 2020; originally announced March 2020.

Comments: Accepted for publication in Science Advances

Probing the thermal state of the intergalactic medium at $z>5$ with the transmission spikes in high-resolution Ly$α$ forest spectra

Authors: Prakash Gaikwad, Michael Rauch, Martin G. Haehnelt, Ewald Puchwein, James S. Bolton, Laura C. Keating, Girish Kulkarni, Vid Iršič, Eduardo Bañados, George D. Becker, Elisa Boera, Fakhri S. Zahedy, Hsiao-Wen Chen, Robert F. Carswell, Jonathan Chardin, Alberto Rorai

Abstract: We compare a sample of five high-resolution, high S/N Ly$α$ 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 tr… ▽ More We compare a sample of five high-resolution, high S/N Ly$α$ 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. △ Less

Submitted 31 March, 2020; v1 submitted 27 January, 2020; originally announced January 2020.

Comments: 20 pages (+10 pages appendices), 12 figures (+14 figures appendices); Accepted in MNRAS; Main results are summarized in Fig. 10, Fig. 12 and Table 3

Quantum mechanics at high redshift -- Modelling Damped Lyman-$α$ absorption systems

Authors: C. C. Lee, J. K. Webb, R. F. Carswell

Abstract: 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-$α$ absorption systems. We have adapted the quantum-mechanical Kramers-Heisenberg model to include th… ▽ More 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-$α$ 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. △ Less

Submitted 15 November, 2019; v1 submitted 7 October, 2019; originally announced October 2019.

Comments: 17 pages, 7 figures, 1 table

Probing IGM accretion onto faint Lyα emitters at z~2.8

Authors: Fakhri S. Zahedy, Michael Rauch, Hsiao-Wen Chen, Robert F. Carswell, Brian Stalder, Antony A. Stark

Abstract: (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$α$ emitters (LAEs) at z=2.7-2.8, found to be aligned in projection with a background QSO over ~250 kpc… ▽ More (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$α$ 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$α$ 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$α$ 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. △ Less

Submitted 22 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS. 13 pages and 6 figures + an Appendix

Constraining the magnetic field on white dwarf surfaces; Zeeman effects and fine structure constant variation

Authors: J. Hu, J. K. Webb, T. R. Ayres, M. B. Bainbridge, J. D. Barrow, M. A. Barstow, J. C. Berengut, R. F. Carswell, V. A. Dzuba, V. V. Flambaum, J. B. Holberg, C. C. Lee, S. P. Preval, N. Reindl, W. -Ü L. Tchang-Brillet

Abstract: 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 measurement… ▽ More 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σ$ 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. △ Less

Submitted 30 December, 2018; originally announced December 2018.

Comments: 10 pages, 5 figures

The UVES Spectral Quasar Absorption Database (SQUAD) Data Release 1: The first 10 million seconds

Authors: Michael T. Murphy, Glenn G. Kacprzak, Giulia A. D. Savorgnan, Robert F. Carswell

Abstract: 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-t… ▽ More 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 Å. 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-$α$ systems (DLAs), absorption-line surveys and time-variable absorption lines. For example, we provide a catalogue of 155 DLAs whose Lyman-$α$ 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. △ Less

Submitted 14 October, 2018; originally announced October 2018.

Comments: 21 pages, 18 figures. Accepted by MNRAS. All final quasar spectra, reduced contributing exposures, and supplementary material available via https://github.com/MTMurphy77/UVES_SQUAD_DR1

A new Measurement of the Intergalactic Temperature at $z\sim 2.55-2.95$

Authors: Alberto Rorai, Robert F. Carswell, Martin G. Haehnelt, George D. Becker, James S. Bolton, Michael T. Murphy

Abstract: 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\ for… ▽ More 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 $γ$ describing the TDR. This approach yields $T_0/ 10^3\, {\rm K}=15.6 \pm 4.4 $ and $γ=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 $γ=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 $γ$ 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. △ Less

Submitted 2 November, 2017; originally announced November 2017.

Comments: Accepted for publication in MNRAS

A precise deuterium abundance: Re-measurement of the z=3.572 absorption system towards the quasar PKS1937-101

Authors: S. Riemer-Sørensen, S. Kotuš, J. K. Webb, K. Ali, V. Dumont, M. T. Murphy, R. F. Carswell

Abstract: 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 f… ▽ More 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. △ Less

Submitted 20 March, 2017; originally announced March 2017.

Comments: 13 pages, 5 figures, 6 tables, accepted by MNRAS

Exploring the Thermal State of the Low-Density Intergalactic Medium at z=3 with an Ultra-High Signal-to-Noise QSO Spectrum

Authors: A. Rorai, G. D. Becker, M. G. Haehnelt, R. F. Carswell, J. S. Bolton, S. Cristiani, V. D'Odorico, G. Cupani, P. Barai, F. Calura, T. -S. Kim, E. Pomante, E. Tescari, M. Viel

Abstract: 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 b… ▽ More 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-$α$ 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. △ Less

Submitted 11 November, 2016; originally announced November 2016.

Comments: Accepted for publication on MNRAS

High-precision limit on variation in the fine-structure constant from a single quasar absorption system

Authors: Srđan M. Kotuš, Michael T. Murphy, Robert F. Carswell

Abstract: 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 ($Δα/α$). A variation of just $\sim$3 parts per million (ppm) would produce detectable velocity shifts between its many strong metal transition… ▽ More 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 ($Δα/α$). 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 $Δα/α$ from a single absorption system to date, $Δα/α=-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 $Δα/α$ 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. △ Less

Submitted 13 September, 2016; originally announced September 2016.

Comments: 27 pages, 10 figures. Accepted by MNRAS. Quasar spectra and absorption profile fits are available at https://doi.org/10.5281/zenodo.51715

Metals in the z~3 intergalactic medium: results from an ultra-high signal-to-noise ratio UVES quasar spectrum

Authors: V. D'Odorico, S. Cristiani, E. Pomante, R. F. Carswell, M. Viel, P. Barai, G. D. Becker, F. Calura, G. Cupani, F. Fontanot, M. G. Haehnelt, T-S. Kim, J. Miralda-Escude, A. Rorai, E. Tescari, E. Vanzella

Abstract: 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\,σ$ to lines with column density down to $\log N_{\rm CIV} \simeq 11.4$ and in 60 per cent of the conside… ▽ More 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\,σ$ 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. △ Less

Submitted 13 September, 2016; v1 submitted 22 August, 2016; originally announced August 2016.

Comments: 23 pages, 28 figures. Accepted by MNRAS on 2016 August 24. Received 2016 August 10; in original form 2016 April 15

Journal ref: 2016 MNRAS 463, 2690

Triple-ionised carbon associated with the low-density neutral hydrogen gas at 1.7 < z < 3.3: the integrated N(HI)-N(CIV) relation

Authors: Tae-Sun Kim, R. F. Carswell, C. Mongardi, A. M. Partl, J. P. Muecket, P. Barai, S. Criatiani

Abstract: 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… ▽ More 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. △ Less

Submitted 18 January, 2016; originally announced January 2016.

Comments: Accepted for publication on MNRAS, 26 pages, 20 figures, 4 tables. On-line materials are found in the submitted civ.tar.gz file: complete Table 2, complete Table 3, complete Table 4, velocity plots civ1.pdf, civ2.pdf, civ3.pdf, civ4.pdf and civ5.pdf

Evidence of bimodal physical properties of intervening, optically-thin CIII absorbers at z ~ 2.5

Authors: Tae-Sun Kim, Robert F. Carswell, Drake Ranquist

Abstract: 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 (ali… ▽ More 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. △ Less

Submitted 3 December, 2015; originally announced December 2015.

Comments: Accepted for publication in MNRAS. For the published paper: 29 pages, 24 figures, 7 tables. For online data, 53 velocity plots (c3_velplot.pdf), profile fitting line parameters (table2.dat, table2.pdf)

Discovery of polarized line emission in SN1006

Authors: W. B. Sparks, J. E. Pringle, R. F. Carswell, K. S. Long, M. Cracraft

Abstract: 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 appro… ▽ More 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. △ Less

Submitted 18 November, 2015; originally announced November 2015.

Comments: ApJ Letters, accepted

The Magellan Uniform Survey of Damped Lyman alpha Systems II: Paucity of Strong Molecular Hydrogen Absorption

Authors: Regina A. Jorgenson, Michael T. Murphy, Rodger Thompson, Robert F. Carswell

Abstract: 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.… ▽ More 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. △ Less

Submitted 3 July, 2014; originally announced July 2014.

Comments: 21 pages, 20 figures, Accepted by MNRAS

Polarization Diagnostics for Cool Core Cluster Emission Lines

Authors: William B. Sparks, James E. Pringle, Robert F. Carswell, G. Mark Voit, Megan Donahue, Misty Cracraft, Eileen T. Meyer, James H. Hough, Nadine Manset

Abstract: 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 filamen… ▽ More 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. △ Less

Submitted 12 November, 2013; originally announced November 2013.

Comments: 49 pages, 14 figures, accepted for publication in The Astrophysical Journal

The evolution of HI and CIV quasar absorption line systems at 1.9 < z < 3.2

Authors: Tae-Sun Kim, Adrian M. Partl, Robert F. Carswell, Volker Müller

Abstract: 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 satu… ▽ More 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). △ Less

Submitted 26 February, 2013; originally announced February 2013.

Comments: 26 pages, 20 figures, accepted for AA, in press

A z=3.045 Lyman alpha emitting halo hosting a QSO and a possible candidate for AGN-triggered star-formation

Authors: Michael Rauch, George D. Becker, Martin G. Haehnelt, Robert F. Carswell, Jean-Rene Gauthier

Abstract: 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 prope… ▽ More 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. △ Less

Submitted 11 February, 2013; originally announced February 2013.

Comments: MNRAS, in press

A high molecular fraction in a sub-damped absorber at z=0.56

Authors: Neil H. M. Crighton, Jill Bechtold, Robert F. Carswell, Romeel Davé, Craig B. Foltz, Buell T. Jannuzi, Simon L. Morris, John M. O'Meara, J. Xavier Prochaska, Joop Schaye, Nicolas Tejos

Abstract: 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 firs… ▽ More 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. △ Less

Submitted 24 April, 2013; v1 submitted 2 October, 2012; originally announced October 2012.

Comments: Final version accepted by MNRAS. Several minor changes compared to the first submitted version, and the conclusions remain unchanged. The data and code used are available here: https://github.com/nhmc/H2

Hundred Thousand Degree Gas in the Virgo Cluster of Galaxies

Authors: W. B. Sparks, J. E. Pringle, R. F. Carswell, M. Donahue, R. Martin, M. Voit, M. Cracraft, N. Manset, J. H. Hough

Abstract: 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 d… ▽ More 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. △ Less

Submitted 22 March, 2012; originally announced March 2012.

Comments: 10 pages, 2 figures; to appear in ApJL

Spatial variation in the fine-structure constant -- new results from VLT/UVES

Authors: Julian A. King, John K. Webb, Michael T. Murphy, Victor V. Flambaum, Robert F. Carswell, Matthew B. Bainbridge, Michael R. Wilczynska, F. Elliot Koch

Abstract: (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 evi… ▽ More (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. △ Less

Submitted 21 February, 2012; originally announced February 2012.

Comments: 47 pages, 35 figures. Accepted for publication by Monthly Notices of the Royal Astronomical Society. Please see http://astronomy.swin.edu.au/~mmurphy/pub.html for an ASCII version of table A1 and the full set of Voigt profile fits for appendix C

The kinetic temperature in a damped Lyman-alpha absorption system in Q2206-199 - an example of the warm neutral medium

Authors: R. F. Carswell, G. D. Becker, R. A. Jorgenson, M. T. Murphy, A. M. Wolfe

Abstract: 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 pr… ▽ More 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. △ Less

Submitted 14 February, 2012; originally announced February 2012.

Comments: 9 pages, 3 figures. MNRAS accepted

Iron and alpha-element Production in the First One Billion Years after the Big Bang

Authors: George D. Becker, Wallace L. W. Sargent, Michael Rauch, Robert F. Carswell

Abstract: 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 ide… ▽ More 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. △ Less

Submitted 21 November, 2011; originally announced November 2011.

Comments: Accepted for publication in ApJ

A cold component and the complex velocity structure of DLA1331+170

Authors: R F Carswell, R A Jorgenson, A M Wolfe, M T Murphy

Abstract: [ABRIDGED] We examine the velocity structure in the gas associated with \ion{H}{1} in the damped Ly$α$ 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 br… ▽ More [ABRIDGED] We examine the velocity structure in the gas associated with \ion{H}{1} in the damped Ly$α$ 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$α$ 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. △ Less

Submitted 14 October, 2010; originally announced October 2010.

Comments: 49 pages, 8 figures, 6 tables. Accepted for publication in MNRAS

Indications of a spatial variation of the fine structure constant

Authors: J. K. Webb, J. A. King, M. T. Murphy, V. V. Flambaum, R. F. Carswell, M. B. Bainbridge

Abstract: 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… ▽ More 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. △ Less

Submitted 31 October, 2011; v1 submitted 23 August, 2010; originally announced August 2010.

Comments: 5 pages, 5 figures, published on 31st October 2011 in Physical Review Letters

Journal ref: Phys. Rev. Lett., 107, 191101, 2011

Cosmological concordance or chemical coincidence? Deuterated molecular hydrogen abundances at high redshift

Authors: J. Tumlinson, A. L. Malec, R. F. Carswell, M. T. Murphy, R. Buning, N. Milutinovic, S. L. Ellison, J. X. Prochaska, R. A. Jorgenson, W. Ubachs, A. M. Wolfe

Abstract: 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… ▽ More 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. △ Less

Submitted 30 June, 2010; originally announced July 2010.

Comments: accepted for publication in ApJ Letters, 5 pages emulate ApJ text

The Low-Redshift Lyman Alpha Forest toward 3C 273

Authors: G. M. Williger, R. F. Carswell, R. J. Weymann, E. B. Jenkins, K. R. Sembach, T. M. Tripp, R. Davé, L. Haberzettl, S. R. Heap

Abstract: 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 den… ▽ More 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. △ Less

Submitted 17 February, 2010; originally announced February 2010.

Comments: 30 pages, 10 figures, MNRAS, in press

Journal ref: MNRAS, v405, 1736 (2010)

New limit on a varying proton-to-electron mass ratio from high-resolution optical quasar spectra

Authors: Adrian L. Malec, Ruth Buning, Michael T. Murphy, Nikola Milutinovic, S. L. Ellison, J. Xavier Prochaska, Lex Kaper, Jason Tumlinson, Robert F. Carswell, Wim Ubachs

Abstract: 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 f… ▽ More 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. △ Less

Submitted 22 January, 2010; originally announced January 2010.

Comments: 6 pages, 3 figures (3 EPS files). Contribution to Proc. IAU XXVIIth General Assembly, Joint Discussion 9, "Are the fundamental constants varying with time?". In P. Molaro, E. Vangioni-Flam, eds, Memorie della Societa Astronomica Italiana (MmSAIt), Vol. 80, pp. 882-887 (2009)

Keck telescope constraint on cosmological variation of the proton-to-electron mass ratio

Authors: Adrian L. Malec, Ruth Buning, Michael T. Murphy, Nikola Milutinovic, S. L. Ellison, J. Xavier Prochaska, Lex Kaper, Jason Tumlinson, Robert F. Carswell, Wim Ubachs

Abstract: 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 transition… ▽ More 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. △ Less

Submitted 22 January, 2010; originally announced January 2010.

Comments: 15 pages, 8 figures (8 EPS files), 3 tables. Accepted by MNRAS. ArXiv copy includes full version of Fig. 1 (additional 8 pages, 7 EPS files). Complete version of Table 1 available at http://astronomy.swin.edu.au/~mmurphy/pub.html

Direct Evidence of Cold Gas in DLA 0812+32B

Authors: Regina A. Jorgenson, Arthur M. Wolfe, J. Xavier Prochaska, Robert F. Carswell

Abstract: 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… ▽ More 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. △ Less

Submitted 24 August, 2009; originally announced August 2009.

Comments: Accepted for publication in ApJ, 5 pages, 8 figures

Stringent null constraint on cosmological evolution of the proton-to-electron mass ratio

Authors: Julian A. King, John K. Webb, Michael T. Murphy, Robert F. Carswell

Abstract: 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 p… ▽ More 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. △ Less

Submitted 17 November, 2008; v1 submitted 28 July, 2008; originally announced July 2008.

Comments: 12 pages. A shorter (4 page) version of this paper has been accepted for publication in Physical Review Letters. Changes in v2: replaced v1 with the version accepted for publication in Physical Review Letters

Journal ref: Phys.Rev.Lett.101:251304,2008

From Espresso to Codex

Authors: J. Liske, L. Pasquini, P. Bonifacio, F. Bouchy, R. F. Carswell, S. Cristiani, M. Dessauges, S. D'Odorico, V. D'Odorico, A. Grazian, R. Garcia-Lopez, M. Haehnelt, G. Israelian, C. Lovis, E. Martin, M. Mayor, P. Molaro, M. T. Murphy, F. Pepe, D. Queloz, R. Rebolo, S. Udry, E. Vanzella, M. Viel, T. Wiklind , et al. (2 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 13 February, 2008; originally announced February 2008.

Comments: To appear in the Proceedings of the Workshop "Science with the VLT in the ELT era", 8-12 October 2007, Garching, A. Moorwood, ed

Possible evidence for an inverted temperature-density relation in the intergalactic medium from the flux distribution of the Lyman-alpha forest

Authors: J. S. Bolton, M. Viel, T. -S. Kim, M. G. Haehnelt, R. F. Carswell

Abstract: 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 (I… ▽ More 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. △ Less

Submitted 19 February, 2008; v1 submitted 14 November, 2007; originally announced November 2007.

Comments: 16 pages, 6 figures, accepted for publication in MNRAS following minor revision. The accepted version includes an expanded discussion of the flux power spectrum

Journal ref: Mon.Not.Roy.Astron.Soc. 386 (2008) 1131-1144

An improved measurement of the flux distribution of the Ly-alpha forest in QSO absorption spectra: the effect of continuum fitting, metal contamination and noise properties

Authors: T. -S. Kim, J. S. Bolton, M. Viel, M. G. Haehnelt, R. F. Carswell

Abstract: 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… ▽ More 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. △ Less

Submitted 13 November, 2007; originally announced November 2007.

Comments: Accepted 2007 September 3 for publication on MNRAS, 21 pages, 13 figures and 6 tables. Data and high-resolution versions of Figs. 2, 3 and 4 and of the paper are available on http://adlibitum.oat.ts.astro.it/viel/PDF/

Journal ref: Mon.Not.Roy.Astron.Soc. 382 (2007) 1657-1674

A large population of metal-rich, compact, intergalactic C IV absorbers - Evidence for poor small-scale metal mixing

Authors: Joop Schaye, Robert F. Carswell, Tae-Sun Kim

Abstract: 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)… ▽ More 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. △ Less

Submitted 21 May, 2007; v1 submitted 26 January, 2007; originally announced January 2007.

Comments: Accepted for publication in MNRAS, 29 pages, 24 figures. Minor changes

Journal ref: Mon.Not.Roy.Astron.Soc.379:1169-1194,2007

Spatial fluctuations in the spectral shape of the UV background at 2<z<3 and the reionization of helium

Authors: James S. Bolton, Martin G. Haehnelt, Matteo Viel, Robert F. Carswell

Abstract: 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 m… ▽ More 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. △ Less

Submitted 29 November, 2005; v1 submitted 9 August, 2005; originally announced August 2005.

Comments: 14 pages, 7 figures, accepted for publication in MNRAS. Minor changes to submitted version, most notably an extra section (5.2) and figure (fig.6)

Journal ref: Mon.Not.Roy.Astron.Soc.366:1378-1390,2006

The Low-Redshift Lyman Alpha Forest toward PKS 0405-123

Authors: G. M. Williger, S. R. Heap, R. J. Weymann, R. Dave, E. Ellingson, R. F. Carswell, T. M. Tripp, E. B. Jenkins

Abstract: 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,… ▽ More 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) △ Less

Submitted 30 August, 2005; v1 submitted 29 May, 2005; originally announced May 2005.

Comments: 59 pages, 6 tables, 14 figures, 1 standalone, ApJ, in press; minor revisions/clarifications for acceptance, including galaxy list

Journal ref: Astrophys.J.636:631-653,2006

Cosmological Evolution of Heavy Element and Molecular Hydrogen Abundances

Authors: Stephen J. Curran, John K. Webb, Michael T. Murphy, Robert F. Carswell

Abstract: 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 tent… ▽ More 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. △ Less

Submitted 27 April, 2004; v1 submitted 15 November, 2003; originally announced November 2003.

Comments: 5 pages, 2 figures. Accepted by MNRAS Letters. v2: Added table summarizing H2-bearing DLA properties, added figure showing [Fe/H] vs. redshift, added more discussion

Journal ref: Mon.Not.Roy.Astron.Soc. 351 (2004) L24

The power spectrum of the flux distribution in the Lyman-alpha forest of a Large sample of UVES QSO Absorption Spectra (LUQAS)

Authors: T. -S. Kim, M. Viel, M. G. Haehnelt, R. F. Carswell, S. Cristiani

Abstract: 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 r… ▽ More 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. △ Less

Submitted 11 May, 2004; v1 submitted 7 August, 2003; originally announced August 2003.

Comments: 18 pages, 9 PS figures, 6 tables. Note that the k-values of the 1D flux power spectrum had been erroneously shifted by half a bin size (in log k) in the previous version. All the other results are unaffected. New tables can be found at http://www.ast.cam.ac.uk/~rtnigm/luqas.htm

Journal ref: Mon.Not.Roy.Astron.Soc.347:355,2004

The effect of (strong) discrete absorption systems on the lyman alpha forest flux power spectrum

Authors: M. Viel, M. G. Haehnelt, R. F. Carswell, T. -S. Kim

Abstract: 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 th… ▽ More 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. △ Less

Submitted 6 August, 2003; originally announced August 2003.

Comments: 6 pages, 3 figures, 1 table. Submitted to MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.349:L33,2004

The Deuterium Abundance in the z=0.7 absorber towards QSO PG1718+4807

Authors: Neil H. M. Crighton, John K. Webb, Robert F. Carswell, Ken M. Lanzetta

Abstract: 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… ▽ More 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. △ Less

Submitted 27 June, 2003; originally announced June 2003.

Comments: 11 pages, 9 figures. Accepted by MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.345:243,2003