Search

Chemical evolution models predict a gradual build-up of $^{13}$C in the universe, based on empirical nuclear reaction rates and assumptions on the properties of stellar populations. However, old metal-poor stars within the Galaxy contain more $^{13}$C than is predicted, suggesting that further refinements to the models are necessary. Gas at high redshift provides important supplementary information at metallicities $-2\lesssim$ [Fe/H] $\lesssim-1$, for which there are only a few measurements in the Galaxy. We obtained new, high-quality, VLT/ESPRESSO observations of the QSO B1331+170 and used them to measure $^{12}$C/$^{13}$C in the damped Lyman-$\alpha$ system (DLA) at $z_{abs}=1.776$, with [Fe/H]=-1.27. AI-VPFIT, an Artificial Intelligence tool based on genetic algorithms and guided by a spectroscopic information criterion, was used to explore different possible kinematic structures of the carbon gas. Three hundred independent AI-VPFIT models of the absorption system were produced using pre-set $^{12}$C/$^{13}$C values, ranging from 4 to 500. Our results show that $^{12}$C/$^{13}$C$=28.5^{+51.5}_{-10.4}$, suggesting a possibility of $^{13}$C production at low metallicity., Comment: 11 pages, 6 figures. Three appendices. To appear in MNRAS

Searches for spacetime variations of fundamental constants have entered an era of unprecedented precision. New, high quality quasar spectra require increasingly refined analytic methods. In this article, a continuation in a series to establish robust and unbiased methodologies, we explore how convergence criteria in non-linear least squares optimisation impact on quasar absorption system measurements of the fine structure constant alpha. Given previous claims for high-precision constraints, we critically examine the veracity of a so-called ``blinding'' approach, in which alpha is fixed at the terrestrial value during the model building process, releasing it as a free parameter only after the ``final'' absorption system kinematic structure has been obtained. We show that this approach results in an extended flat canyon in chi squared-alpha space, such that convergence is unlikely to be reached, even after as many as 1000 iterations. The fix is straightforward: alpha must be treated as a free parameter from the earliest possible stages of absorption system model building. The implication of the results presented here is that all previous measurements that have used initially-fixed alpha should be reworked., Comment: 10 pages, 13 figures. Submitted for publication

We report unbiased AI measurements of the fine structure constant alpha in two proximate absorption regions in the spectrum of the quasar HE0515-4414. The data are high resolution, high signal to noise, and laser frequency comb calibrated, obtained using the ESPRESSO spectrograph on the VLT. The high quality of the data and proximity of the regions motivate a differential comparison, exploring the possibility of spatial variations of fundamental constants, as predicted in some theories. We show that if the magnesium isotopic relative abundances are terrestrial, the fine structure constants in these two systems differ at the 7-sigma level. A 3-sigma discrepancy between the two measurements persists even for the extreme non-terrestrial case of 100% ^{24}Mg, if shared by both systems. However, if Mg isotopic abundances take independent values in these two proximate systems, one terrestrial, the other with no heavy isotopes, both can be reconciled with a terrestrial alpha, and the discrepancy between the two measurements falls to 2-sigma. We discuss varying constant and varying isotope interpretations and resolutions to this conundrum for future high precision measurements., Comment: 6 pages, 3 figures

Just 10 recurrent novae (RNe) - which erupt repeatedly on timescales shorter than one century - are known in our Galaxy. The most extreme RN known (located in the Andromeda galaxy), M31N 2008-12a, undergoes a nova eruption every year, and is surrounded by a vast nova "super-remnant", 134 pc in extent. Simulations predict that all RNe should be surrounded by similar vast shells, but previous searches have failed to detect them. KT Eri has recently been suggested to be a RN, and we have used the Condor Array Telescope to image its environs through multiple narrowband filters. We report the existence of a large ($\sim$ 50 pc diameter), H$\,\alpha$-bright shell centered on KT Eri, exactly as predicted. This strongly supports the claim that KT Eri is the 11th Galactic recurrent nova, and only the second nova known to be surrounded by a super-remnant. SALT spectra of the super-remnant demonstrate that its velocity width is consistent with that of M31-2008-12a., Comment: 12 pages, 8 figures

We used the Condor Array Telescope to obtain deep imaging observations through the luminance filter of the entirety of the NGC 5866 Group, including a very extended region surrounding the galaxy NGC 5907 and its stellar stream. We find that the stellar stream consists of a single curved structure that stretches $220$ kpc from a brighter eastern stream to a fainter western stream that bends to the north and then curls back toward the galaxy. This result runs contrary to a previous claim of a second loop of the stellar stream but is consistent with another previous description of the overall morphology of the stream. We further find that: (1) an extension of the western stream appears to bifurcate near its apex, (2) there is an apparent gap of $\approx 6$ kpc in the western stream due east of the galaxy, (3) contrary to a previous claim, there is no evidence of the remnant of a progenitor galaxy within the eastern stream, although (4) there are many other possible progenitor galaxies, (5) there is another structure that, if it is at the distance of the galaxy, stretches 240 kpc and contains two very large, very low-surface-brightness "patches" of emission, one of which was noted previously and another of which was not. We note the number and variety of stellar streams in the vicinity of NGC 5907 and the apparent gap in the western stream, which may be indicative of a dark subhalo or satellite in the vicinity of the galaxy., Comment: 15 pages, 12 figures

The "Condor Array Telescope" or "Condor" is a high-performance "array telescope" comprised of six apochromatic refracting telescopes of objective diameter 180 mm, each equipped with a large-format, very low-read-noise ($\approx 1.2$ e$^-$), very rapid-read-time ($< 1$ s) CMOS camera. Condor is located at a very dark astronomical site in the southwest corner of New Mexico, at the Dark Sky New Mexico observatory near Animas, roughly midway between (and more than 150 km from either) Tucson and El Paso. Condor enjoys a wide field of view ($2.29 \times 1.53$ deg$^2$ or 3.50 deg$^2$), is optimized for measuring both point sources and extended, very low-surface-brightness features, and for broad-band images can operate at a cadence of 60 s (or even less) while remaining sky-noise limited with a duty cycle near 100\%. In its normal mode of operation, Condor obtains broad-band exposures of exposure time 60 s over dwell times spanning dozens or hundreds of hours. In this way, Condor builds up deep, sensitive images while simultaneously monitoring tens or hundreds of thousands of point sources per field at a cadence of 60 s. Condor is also equipped with diffraction gratings and with a set of He II 468.6 nm, [O III] 500.7 nm, He I 587.5 nm, H$\alpha$ 656.3 nm, [N II] 658.4 nm, and [S II] 671.6 nm narrow-band filters, allowing it to address a variety of broad- and narrow-band science issues. Given its unique capabilities, Condor can access regions of "astronomical discovery space" that have never before been studied. Here we introduce Condor and describe various aspects of its performance., Comment: 23 pages, 16 figures

This article reviews the two major recent developments that significantly improved cosmological measurements of fundamental constants derived from high resolution quasar spectroscopy. The first one is the deployment of astronomical Laser Frequency Combs on high resolution spectrographs and the second one is the development of spectral analysis tools based on Artificial Intelligence methods. The former all but eliminated the previously dominant source of instrumental uncertainty whereas the latter established optimal methods for measuring the fine structure constant ($\alpha$) in quasar absorption systems. The methods can be used on data collected by the new ESPRESSO spectrograph and the future ANDES spectrograph on the Extremely Large Telescope to produce unbiased $\Delta\alpha/\alpha$ measurements with unprecedented precision., Comment: 10 pages, part of the HACK100 conference proceedings

Randomly oriented flattened spheroids have been used to describe a broad range of astrophysical phenomena. Here we use this geometric approach to derive equations representing lines of sight through quasar absorption clouds to constrain cloud sizes., Comment: 10 pages, 1 figure

Searches for variations of fundamental constants require accurate measurement errors. There are several potential sources of errors and quantifying each one accurately is essential. This paper addresses one source of uncertainty relating to measuring the fine structure constant on white dwarf surfaces. Detailed modelling of photospheric absorption lines requires knowing the underlying spectral continuum level. Here we describe the development of a fully automated, objective, and reproducible continuum estimation method, based on fitting cubic splines to carefully selected data regions. Example fits to the Hubble Space Telescope spectrum of the white dwarf G191-B2B are given. We carry out measurements of the fine structure constant using two continuum models. The results show that continuum placement variations result in small systematic shifts in the centroids of narrow photospheric absorption lines which impact significantly on fine structure constant measurements. This effect must therefore be included in the overall error budget of future measurements. Our results also suggest that continuum placement variations should be investigated in other contexts, including fine structure constant measurements in stars other than white dwarfs, quasar absorption line measurements of the fine structure constant, and quasar measurements of cosmological redshift drift., Comment: 10 pages, 4 figures. 4 additional files provided as supplementary material. Submitted to MNRAS 1 Dec 2022

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., Comment: 17 pages, 7 figures, and 5 tables. Submitted to MNRAS 1 Dec 2022

Cosmic Probes of Fundamental Physics take two primary forms: Very high energy particles (cosmic rays, neutrinos, and gamma rays) and gravitational waves. Already today, these probes give access to fundamental physics not available by any other means, helping elucidate the underlying theory that completes the Standard Model. The last decade has witnessed a revolution of exciting discoveries such as the detection of high-energy neutrinos and gravitational waves. The scope for major developments in the next decades is dramatic, as we detail in this report., Comment: Report of theTopical Group on Cosmic Probes of Fundamental Physics, for the U.S. decadal Particle Physics Planning Exercise (Snowmass 2021)

In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $\sigma_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,\sigma$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $\Omega_m$, and the amplitude or rate of the growth of structure ($\sigma_8,f\sigma_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged], Comment: Contribution to Snowmass 2021. 224 pages, 27 figures. Accepted for publication in JHEAp

Recent advances in spectroscopic instrumentation and calibration methods dramatically improve the quality of quasar spectra. Supercomputer calculations show that, at high spectral resolution, procedures used in some previous analyses of spacetime variations of fundamental constants are likely to generate spurious measurements, biased systematically towards a null result. Developments in analysis methods are also summarised and a prescription given for the analysis of new and forthcoming data., Comment: 17 pages, 6 figures, accepted for publication in "The Universe". Contribution for the AlteCosmoFun'21 conference volume, 6-10 September 2022, Szczecin, Poland

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., Comment: 4 pages, 2 figures. Submitted to MNRAS 23 Dec 2021, accepted 24 Dec 2021

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., Comment: 15 pages, 7 figures, to appear in MNRAS

Modifications were made to the Sloan Digital Sky Survey's Baryonic Oscillations Spectroscopic Survey (SDSS/BOSS) optical fibres assigned to quasar targets in order to improve the signal-to-noise ratio in the Lyman-alpha forest. However, the penalty for these modifications is that quasars observed in this way require additional flux correction procedures in order to recover the correct spectral shapes. In this paper we describe such a procedure, based on the geometry of the problem and other observational parameters. Applying several correction methods to four SDSS quasars with multiple observations permits a detailed check on the relative performances of the different flux correction procedures. We contrast our method (which takes into account a wavelength dependent seeing profile) with the BOSS pipeline approach (which does not). Our results provide independent confirmation that the geometric approach employed in the SDSS pipeline works well, although with room for improvement. By separating the contributions from four effects we are able to quantify their relative importance. Most importantly, we demonstrate that wavelength dependence has a significant impact on the derived spectral shapes and thus should not be ignored., Comment: 24 pages, 17 figures, submitted to A&A

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 \alpha, 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 \alpha, 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 \alpha 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., Comment: 17 pages, 11 figures

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., Comment: 13 pages, 7 figures

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 $\alpha$ 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 $\alpha$ 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 $\alpha$ 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 $\alpha$ measurements. We constrain spatial fluctuations in $\alpha$ to be $\Delta\alpha/\alpha \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 $\Delta\alpha/\alpha=-0.27\pm2.41\times10^{-6}$, consistent with no variation., Comment: Published in MNRAS, 500 1 1. 10 pages, 7 figures. Online supplementary material is available upon request

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., Comment: Accepted by MNRAS, 1 April 2021

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., Comment: Accepted for publication in Science Advances

Astrocombs are ideal spectrograph calibrators whose limiting precision can be derived using a second, independent, astrocomb system. We therefore analyse data from two astrocombs (one 18 GHz and one 25 GHz) used simultaneously on the HARPS spectrograph at the European Southern Observatory. The first aim of this paper is to quantify the wavelength repeatability achieved by a particular astrocomb. The second aim is to measure wavelength calibration consistency between independent astrocombs, that is to place limits or measure any possible zero-point offsets. We present three main findings, each with important implications for exoplanet detection, varying fundamental constant and redshift drift measurements. Firstly, wavelength calibration procedures are important: using multiple segmented polynomials within one echelle order results in significantly better wavelength calibration compared to using a single higher-order polynomial. Segmented polynomials should be used in all applications aimed at precise spectral line position measurements. Secondly, we found that changing astrocombs causes significant zero-point offsets ($\approx 60{\rm cms}^{-1}$ in our raw data) which were removed. Thirdly, astrocombs achieve a precision of $\lesssim 4{\rm cms}^{-1}$ in a single exposure ($\approx 10\% $ above the measured photon-limited precision) and $1 {\rm cms}^{-1}$ when time-averaged over a few hours, confirming previous results. Astrocombs therefore provide the technological requirements necessary for detecting Earth-Sun analogues, measuring variations of fundamental constants and the redshift drift., Comment: 16 pages, 17 figures, accepted for publication in MNRAS

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

We recently presented a new "artificial intelligence" method for the analysis of high-resolution absorption spectra (Bainbridge and Webb, Mon. Not. R. Astron. Soc. 2017, 468,1639-1670). This new method unifies three established numerical methods: a genetic algorithm (GVPFIT); non-linear least-squares optimisation with parameter constraints (VPFIT); and Bayesian Model Averaging (BMA). In this work, we investigate the performance of GVPFIT and BMA over a broad range of velocity structures using synthetic spectra. We found that this new method recovers the velocity structures of the absorption systems and accurately estimates variation in the fine structure constant. Studies such as this one are required to evaluate this new method before it can be applied to the analysis of large sets of absorption spectra. This is the first time that a sample of synthetic spectra has been utilised to investigate the analysis of absorption spectra. Probing the variation of nature's fundamental constants (such as the fine structure constant), through the analysis of absorption spectra, is one of the most direct ways of testing the universality of physical laws. This "artificial intelligence" method provides a way to avoid the main limiting factor, i.e., human interaction, in the analysis of absorption spectra., Comment: 9 pages, 5 figures, published on 5 April 2017 in Universe

Hot white dwarf stars are the ideal probe for a relationship between the fine-structure constant and strong gravitational fields, providing us with an opportunity for a direct observational test. We study a sample of hot white dwarf stars, combining far-UV spectroscopic observations, atomic physics, atmospheric modelling and fundamental physics, in the search for variation in the fine structure constant. This variation manifests as shifts in the observed wavelengths of absorption lines, such as quadruply ionized iron (FeV) and quadruply ionized nickel (NiV), when compared to laboratory wavelengths. Berengut et al. (Phys. Rev. Lett. 2013, 111, 010801) demonstrated the validity of such an analysis using high-resolution Hubble Space Telescope (HST)/Space Telescope Imaging Spectrograph (STIS) spectra of G191-B2B. We have made three important improvements by: (a) using three new independent sets of laboratory wavelengths, (b) analysing a sample of objects, and (c) improving the methodology by incorporating robust techniques from previous studies towards quasars (the Many Multiplet method). A successful detection would be the first direct measurement of a gravitational field effect on a bare constant of nature. Here we describe our approach and present preliminary results from nine objects using both FeV and NiV., Comment: 6 pages, 1 figure, published on 30th March 2017 in Universe as part of VARCOSMOFUN'16 proceedings

Spectra observed with the Ultraviolet and Visual Echelle Spectrograph (UVES) on the European Southern Observatory's VLT exhibit long-range wavelength distortions. These distortions impose a systematic error on high-precision measurements of the fine-structure constant, $\alpha$, derived from intervening quasar absorption systems. If the distortion is modelled using a model that is too simplistic, the resulting bias in $\Delta\alpha/\alpha$ away from the true value can be larger than the statistical uncertainty on the $\alpha$ measurement. If the effect is ignored altogether, the same is true. If the effect is modelled properly, accounting for the way in which final spectra are generally formed from the co-addition of exposures made at several different instrumental settings, the effect can be accurately removed and the correct $\Delta\alpha/\alpha$ recovered., Comment: 8 pages, 7 figures, 3 tables

A new and automated method is presented for the analysis of high-resolution absorption spectra. Three established numerical methods are unified into one "artificial intelligence" process: a genetic algorithm (GVPFIT); non-linear least-squares with parameter constraints (VPFIT); and Bayesian Model Averaging (BMA). The method has broad application but here we apply it specifically to the problem of measuring the fine structure constant at high redshift. For this we need objectivity and reproducibility. GVPFIT is also motivated by the importance of obtaining a large statistical sample of measurements of $\Delta\alpha/\alpha$. Interactive analyses are both time consuming and complex and automation makes obtaining a large sample feasible. In contrast to previous methodologies, we use BMA to derive results using a large set of models and show that this procedure is more robust than a human picking a single preferred model since BMA avoids the systematic uncertainties associated with model choice. Numerical simulations provide stringent tests of the whole process and we show using both real and simulated spectra that the unified automated fitting procedure out-performs a human interactive analysis. The method should be invaluable in the context of future instrumentation like ESPRESSO on the VLT and indeed future ELTs. We apply the method to the $z_{abs} = 1.8389$ absorber towards the $z_{em} = 2.145$ quasar J110325-264515. The derived constraint of $\Delta\alpha/\alpha = 3.3 \pm 2.9 \times 10^{-6}$ is consistent with no variation and also consistent with the tentative spatial variation reported in Webb et al (2011) and King et al (2012)., Comment: 34 pages, 10 figures, 16 tables. Submitted to MNRAS

We present an analysis of 23 absorption systems along the lines of sight towards 18 quasars in the redshift range of $0.4 \leq z_{abs} \leq 2.3$ observed on the Very Large Telescope (VLT) using the Ultraviolet and Visual Echelle Spectrograph (UVES). Considering both statistical and systematic error contributions we find a robust estimate of the weighted mean deviation of the fine-structure constant from its current, laboratory value of $\Delta\alpha/\alpha=\left(0.22\pm0.23\right)\times10^{-5}$, consistent with the dipole variation reported in Webb et al. and King et al. This paper also examines modelling methodologies and systematic effects. In particular we focus on the consequences of fitting quasar absorption systems with too few absorbing components and of selectively fitting only the stronger components in an absorption complex. We show that using insufficient continuum regions around an absorption complex causes a significant increase in the scatter of a sample of $\Delta\alpha/\alpha$ measurements, thus unnecessarily reducing the overall precision. We further show that fitting absorption systems with too few velocity components also results in a significant increase in the scatter of $\Delta\alpha/\alpha$ measurements, and in addition causes $\Delta\alpha/\alpha$ error estimates to be systematically underestimated. These results thus identify some of the potential pitfalls in analysis techniques and provide a guide for future analyses., Comment: 14 pages, 7 figures. Published by the Monthly Notices of the Royal Astronomical Society. v2: Updated reference Dumont & Webb (2016)

(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., Comment: 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

Molecular hydrogen transitions in quasar spectra can be used to constrain variation in the proton-to-electron mass ratio, $\mu\equiv m_p/m_e$, at high redshifts ($z\gtrsim 2$). We present here an analysis of a new spectrum of the quasar Q0528$-$250, obtained on VLT/UVES (the Ultraviolet and Visual Echelle Spectrograph, on the Very Large Telescope), and analyse the well-known H$_2$ absorber at $z=2.811$ in this spectrum. For the first time we detect HD (deuterated molecular hydrogen) in this system with a column density of $\log_{10}(N/\mathrm{cm^{-2}})=13.27 \pm 0.07$; HD is sensitive to variation in $\mu$, and so we include it in our analysis. Using 76 H$_2$ and 7 HD transitions we constrain variation in $\mu$ from the current laboratory value to be $\Delta\mu/\mu = (0.3\pm 3.2_\mathrm{stat} \pm 1.9_\mathrm{sys})\times 10^{-6}$, which is consistent with no cosmological variation in $\mu$, as well as with previous results from other H$_2$/HD absorbers. The main sources of systematic uncertainty relate to accurate wavelength calibration of the spectra and the re-dispersion of multiple telescope exposures onto the one pixel grid., Comment: 30 pages, 6 figures + 14 supplementary figures. Accepted for publication for MNRAS

The Keck telescope's HIRES spectrograph has previously provided evidence for a smaller fine-structure constant, alpha, compared to the current laboratory value, in a sample of 143 quasar absorption systems: da/a=(-0.57+/-0.11)x10^{-5}. This was based on a variety of metal-ion transitions which, if alpha varies, experience different relative velocity shifts. This result is yet to be robustly contradicted, or confirmed, by measurements on other telescopes and spectrographs; it remains crucial to do so. It is also important to consider new possible instrumental systematic effects which may explain the Keck/HIRES results. Griest et al. (2009, arXiv:0904.4725v1) recently identified distortions in the echelle order wavelength scales of HIRES with typical amplitudes +/-250m/s. Here we investigate the effect such distortions may have had on the Keck/HIRES varying alpha results. We demonstrate that they cause a random effect on da/a from absorber to absorber because the systems are at different redshifts, placing the relevant absorption lines at different positions in different echelle orders. The typical magnitude of the effect on da/a is ~0.4x10^{-5} per absorber which, compared to the median error on da/a in the sample, ~1.9x10^{-5}, is relatively small. Consequently, the weighted mean value changes by less than 0.05x10^{-5} if the corrections we calculate are applied. Nevertheless, we urge caution, particularly for analyses aiming to achieve high precision da/a measurements on individual systems or small samples, that a much more detailed understanding of such intra-order distortions and their dependence on observational parameters is important if they are to be avoided or modelled reliably. [Abridged], Comment: 9 pages, 2 figures, 2 tables. Invited contribution to Proc. IAU XXVIIth General Assembly, Joint Discussion 9, "Are the fundamental constants varying with time?". To appear in P. Molaro, E. Vangioni-Flam, eds, Memorie della Societa Astronomica Italiana (MmSAIt), Vol. 80. Complete version of Table 1 available at http://astronomy.swin.edu.au/~mmurphy/pub.html

Recent attempts to constrain cosmological variation in the fine structure constant, alpha, using quasar absorption lines have yielded two statistical samples which initially appear to be inconsistent. One of these samples was subsequently demonstrated to not pass consistency tests; it appears that the optimisation algorithm used to fit the model to the spectra failed. Nevertheless, the results of the other hinge on the robustness of the spectral fitting program VPFIT, which has been tested through simulation but not through direct exploration of the likelihood function. We present the application of Markov Chain Monte Carlo (MCMC) methods to this problem, and demonstrate that VPFIT produces similar values and uncertainties for (Delta alpha)/(alpha), the fractional change in the fine structure constant, as our MCMC algorithm, and thus that VPFIT is reliable., Comment: IAU 2009 JD9 conference proceedings. 6 pages, 2 figures. MmSAI, vol. 80 in press Paolo Molaro & Elisabeth Vangioni eds

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., Comment: 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

In their Letter [Phys. Rev. Lett. 92, 121302 (2004)] (also [Astron. Astrophys. 417, 853 (2004)]), Srianand et al. analysed optical spectra of heavy-element species in 23 absorption systems along background quasar sight-lines, reporting limits on relative variations in the fine-structure constant: da/a=(-0.06+/-0.06) x 10^{-5}. Here we demonstrate basic flaws in their analysis, using the same data and absorption profile fits, which led to spurious values of da/a and significantly underestimated uncertainties. We conclude that these data and fits offer no stringent test of previous evidence for a varying alpha. In their Reply (arXiv:0711.1742) to this Comment, Srianand et al. state or argue several points regarding their original analysis and our new analysis. We discuss these points here, dismissing all of them because they are demonstrably incorrect or because they rely on a flawed application of simple statistical arguments., Comment: 1+2 pages, 1 EPS figure. Page 1 accepted as PRL Comment on arXiv:astro-ph/0402177 . Further details available in arXiv:astro-ph/0612407 . v2: Added critical discussion of Reply from Srianand et al. (arXiv:0711.1742)

We critically review the current null results on a varying fine-structure constant, alpha, derived from VLT/UVES quasar absorption spectra, focusing primarily on the many-multiplet analysis of 23 absorbers from which Chand et al. (2004) reported a weighted mean relative variation of da/a=(-0.06+/-0.06) x 10^{-5}. Our analysis of the same reduced data, using the same fits to the absorption profiles, yields very different individual da/a values with uncertainties typically larger by a factor of ~3. We attribute the discrepancies to flawed parameter estimation techniques in the original analysis and demonstrate that the original da/a values were strongly biased towards zero. Were those flaws not present, the input data and spectra should have given a weighted mean of da/a=(-0.44+/-0.16) x 10^{-5}. Although this new value does reflect the input spectra and fits (unchanged from the original work -- only our analysis is different), we do not claim that it supports previous Keck/HIRES evidence for a varying alpha: there remains significant scatter in the individual da/a values which may stem from the overly simplistic profile fits in the original work. Allowing for such additional, unknown random errors by increasing the uncertainties on da/a to match the scatter provides a more conservative weighted mean, da/a=(-0.64+/-0.36) x 10^{-5}. We highlight similar problems in other current UVES constraints on varying alpha and argue that comparison with previous Keck/HIRES results is premature., Comment: 11 pages, 7 figures (10 EPS files). Accepted by MNRAS. v3: Added figure summarizing all current UVES constraints; added references and corrected typos

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., Comment: 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

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., Comment: 11 pages, 9 figures. Accepted by MNRAS

We present an analysis of galaxy and QSO absorption line pairs toward 24 QSOs at redshifts between z~0.2 and 1 in an effort to establish the relationship between galaxies and absorption lines in physical proximity to QSOs. We demonstrate the existence of a galaxy proximity effect, in that galaxies in the vicinities of QSOs do not show the same incidence and extent of gaseous envelopes as galaxies far from QSOs. We show that the galaxy proximity effect exists to galaxy-QSO velocity separations of ~ 3000 km/s, much larger than the size of a typical cluster (~ 1000 km/s), i.e. it is more comparable to the scale of the sphere of influence of QSO ionizing radiation rather than the scale of galaxy-QSO clustering. This indicates that the QSO ionizing radiation rather than some dynamical effect from the cluster environment is responsible for the galaxy proximity effect. We combine previous findings that (1) many or most Lya absorption lines arise in extended galaxy envelopes, and (2) galaxies cluster around QSOs to show that the magnitude of the Lya forest proximity effect is underestimated. Consequently, determinations of the UV ionizing background intensity using the proximity effect are likely overestimated. We use the galaxy-QSO cross-correlation function measured from our data to estimate the magnitude of this overestimate and find that it could be as high as a factor of 20 at z<1. This can have strong implications for models of the origin and evolution of the ionizing background, and may indicate that QSOs produce sufficient ionizing flux at all redshifts to account for the entire background radiation field., Comment: 26 pages, 3 figures, to appear in ApJ, October 20, 2001

(Abridged) We present results of a program to obtain and analyze HST WFPC2 images and ground-based images of galaxies identified in an imaging and spectroscopic survey of faint galaxies in fields of HST spectroscopic target QSOs. We confirm and improve the results presented by Lanzetta et al. (1995) and Chen et al. (1998) that (1) extended gaseous envelopes are a common and generic feature of galaxies of a wide range of luminosity and morphological type, (2) the extent of tenuous gas (N(H I) > 10^{14} cm^{-2} around galaxies scales with galaxy B-band luminosity as r\propto L_B^{0.39+/-0.09}, and (3) galaxy interactions do not play an important role in distributing tenuous gas around galaxies in most cases. We further demonstrate that (4) the gaseous extent of galaxies scales with galaxy K-band luminosity as r\propto L_K^{0.28+/-0.08}, and (5) tenuous gas around typical L_* galaxies is likely to be distributed in spherical halos of radius ~ 180 h^{-1} kpc of covering factor of nearly unity. Because galaxies of all morphological types possess extended gaseous halos and because the extent of tenuous gas around galaxies scales with galaxy K-band luminosity, we argue that galaxy mass---rather than recent star-formation activity---is likely to be the dominant factor that determines the extent of tenuous gas around galaxies. Nevertheless, applying the scaling relationship between the extent of Lya absorbing gas around galaxies and galaxy B-band luminosity, the results of our analysis also suggest that the number density evolution of Lya absorption systems may serve to constrain the evolution of the comoving galaxy B-band luminosity density (at least for the redshift interval between z~0 and z~1 that has been studied in our survey)., Comment: 46 pages, 7 figures, to appear in ApJ, October 1, 2001

(Abridged) We report the discovery of extended CIV gaseous envelopes around galaxies of a wide range of luminosity and morphological type. First, we show that CIV absorption systems are strongly clustered around galaxies on velocity scales of v < 250 km/s and impact parameter scales of rho < 100 h^{-1} kpc but not on larger velocity or impact parameter scales. Next, adopting measurements of galaxy properties presented in previous papers, we examine how properties of the CIV absorption systems depend on properties of the galaxies. On the basis of 14 galaxy and absorber pairs and 36 galaxies that do not produce corresponding CIV absorption lines to within sensitive upper limits, we find that: (1) Galaxies of a range of morphological type and luminosity appear to possess extended CIV gaseous envelopes of radius R ~ 100 h^{-1} kpc, with abrupt boundaries between the CIV absorbing and non-absorbing regions. (2) The extent of CIV-absorbing gas around galaxies scales with galaxy B-band luminosity as R \propto L_B^{0.5 +/- 0.1} but does not depend strongly on galaxy surface brightness, redshift, or morphological type. And (3) the covering factor of CIV clouds within ~ 100 h^{-1} kpc of galaxies is nearly unity, but there is a large scatter in the mean number of clouds encountered along the line of sight. The most significant implication of the study is that galaxies of a wide range of luminosity and morphological type are surrounded by chemically enriched gas that extends for at least ~ 100 h^{-1} kpc. We consider various scenarios that may have produced metals at large galactic distance and conclude that accreting satellites are most likely to be responsible for chemically enriched gas at large galactic distances to regular looking galaxies., Comment: 19 pages, 3 figures, to appear in ApJ, July 20 2001

We use the dynamics of a galaxy, set up initially at a constant proper distance from an observer, to derive and illustrate two counter-intuitive general relativistic results. Although the galaxy does gradually join the expansion of the universe (Hubble flow), it does not necessarily recede from us. In particular, in the currently favored cosmological model, which includes a cosmological constant, the galaxy recedes from the observer as it joins the Hubble flow, but in the previously favored cold dark matter model, the galaxy approaches, passes through the observer, and joins the Hubble flow on the opposite side of the sky. We show that this behavior is consistent with the general relativistic idea that space is expanding and is determined by the acceleration of the expansion of the universe -- not a force or drag associated with the expansion itself. We also show that objects at a constant proper distance will have a nonzero redshift; receding galaxies can be blueshifted and approaching galaxies can be redshifted., Comment: 8 pages including 6 figures, to appear in Am. J. Phys., 2003. Reference added in postscript

Although the extra-solar planets discovered so far are of the giant, gaseous, type, the increased sensitivity of future surveys will result in the discovery of lower mass planets. The detection of O2 in the atmosphere of a rocky extra-solar planet would be a potential indicator of a life. In this paper we address the specific issue of whether we would be able to detect the O2 A-band absorption feature in the atmosphere of a planet similar to the Earth, if it were in orbit around a nearby star. Our method is empirical, in that we use observations of the Earth's O2 A-band, with a simple geometric modification for a transiting extra-solar planet, allowing for limb-darkening of the host star. We simulate the spectrum of the host star with the superposed O2 A-band absorption of the transiting planet, assuming a spectral resolution of 7 km/s (typical of current echelle spectrographs), for a range of spectral signal-to-noise ratios. The main result is that we could reliably detect the O2 A-band of the transiting planet for host stars with radii 0.3 solar or less. However, using existing instrumentation and 8m telescopes, this requires target M stars with m(V) of approximately 10 or brighter for integration times of about 10 hours or less. The number of such stars over the sky is small. Larger aperture telescopes and/or improved instrumentation efficiency would enable surveys of M stars down to m(V) = 13 and greatly improve the chances of discovering life elsewhere., Comment: 14 pages, 6 figures, to appear in Publ. Astron. Soc. Austr

Here we describe results of our attempt to determine what types of galaxies are responsible for the Lyman-alpha forest absorption systems, based on our ongoing imaging and spectroscopic survey of faint galaxies in fields of HST spectroscopic target QSOs. Our primary conclusions are that the bulk of the Lyman-alpha forest arises in more or less normal galaxies (that span the normal range of luminosity and morphology) and that any ``unseen'' low surface brightness galaxies are unlikely to contribute significantly to the luminosity density of the universe., Comment: 8 pages, 4 figures, invited review for IAU Colloquium 171, The Low Surface Brightness Universe

We present new statistical study on the origin of Lyman-alpha absorption systems that are detected in QSO fields. The primary objective is to investigate environmental dependence of the relationship between galaxies and Lyman-alpha absorption systems. We find that galaxies exhibit no environmental dependence in giving rise to Lyman-alpha absorbers in their halos. Conversely, Lyman-alpha absorbers do not exhibit preference as to the environment where they reside., Comment: 4 pages, 1 figure, to appear in proceedings of the Xth Rencontres de Blois on "The Birth of Galaxies," LaTex style file included