Quasar absorption lines provide a precise test of whether the fine-structure constant, α, is the same in different places and through cosmological time. We present a new analysis of a large sample of ...quasar absorption-line spectra obtained using the Ultraviolet and Visual Echelle Spectrograph (UVES) on the Very Large Telescope (VLT) in Chile. We apply the many-multiplet method to derive values of Δα/α≡ (α
z
−α0)/α0 from 154 absorbers, and combine these values with 141 values from previous observations at the Keck Observatory in Hawaii. In the VLT sample, we find evidence that α increases with increasing cosmological distance from Earth. However, as previously shown, the Keck sample provided evidence for a smaller α in the distant absorption clouds. Upon combining the samples, an apparent variation of α across the sky emerges which is well represented by an angular dipole model pointing in the direction RA = 17.3 ± 1.0 h and Dec. =−61°± 10°, with amplitude
. The dipole model is required at the 4.1σ statistical significance level over a simple monopole model where α is the same across the sky (but possibly different from the current laboratory value). The data sets reveal remarkable consistencies: (i) the directions of dipoles fitted to the VLT and Keck samples separately agree; (ii) the directions of dipoles fitted to z < 1.6 and z > 1.6 cuts of the combined VLT+Keck samples agree; and (iii) in the equatorial region of the dipole, where both the Keck and VLT samples contribute a significant number of absorbers, there is no evidence for inconsistency between Keck and VLT. The amplitude of the dipole is clearly larger at higher redshift. 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) × 10−6 GLyr−1 and the model is significant at the 4.2σ confidence level over the null model (Δα/α≡ 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 α. To the best of our knowledge, this result is not in conflict with any other observational or experimental result.
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_{\rm H\, \small {I}}$ and the Doppler parameter b of the absorption lines in the Lyα forest. The first measurement is obtained using a novel Bayesian scheme that takes into account the statistical correlations between the parameters characterizing the lower cut-off of the $b\hbox{--}N_{\rm H\, \small {I}}$ distribution and the power-law parameters T0 and γ describing the TDR. This approach yields T0/103 K = 15.6 ± 4.4 and γ = 1.45 ± 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\hbox{--}N_{\rm H\, \small {I}}$ 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 T0/103 K = 14.6 ± 3.7 and γ = 1.37 ± 0.17 in good agreement with the first measurement. Our Bayesian analysis reveals strong anticorrelations between the inferred T0 and γ for both methods as well as an anticorrelation of the inferred T0 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.
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 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.
ABSTRACT
We compare a sample of five high-resolution, high S/N Ly α forest spectra of bright 6 < z < ∼6.5 QSOs aimed at spectrally resolving the last remaining transmission spikes at z > 5 with ...those obtained from mock absorption spectra from the Sherwoodand Sherwood–Relics simulation 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 characterize 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 prominent in low temperature models of the IGM. The width distribution of the observed transmission spikes, which require high spectral resolution (≤ 8 km s−1) to be resolved, is reproduced for optically thin simulations with a temperature at mean density of T0 = (11 000 ± 1600, 10 500 ± 2100, 12 000 ± 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 ∼ 5.3, the width distribution of the observed transmission spikes is consistent with the range of T0 caused by spatial fluctuations in the temperature–density relation.
ABSTRACT
New observations of the quasar HE0515−4414 have been made, aided by the Laser Frequency Comb (LFC), using the HARPS spectrograph on the ESO 3.6m telescope. We present three important ...advances for α measurements in quasar absorption spectra from these observations. First, 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 Δα/α ≤ 9 × 10−5 on scales $\approx \!\! {20}\, {\rm km\, s}^{-1}$, corresponding to $\approx 25\,$kpc if the $z$abs = 1.15 system arises in a 1Mpc cluster. Collectively, the 40 measurements yield Δα/α = −0.27 ± 2.41 × 10−6, consistent with no variation.
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 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 a narrow spectral line impacts only on a very small section of the overall data. 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 data set and/or have widely varying sensitivities.
ABSTRACT
This paper describes the optimization 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 optimization 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 summarize how unbiased modelling depends on using an appropriate information criterion to guard against overfitting or underfitting. 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 space–time variations in fundamental constants and attempts to detect cosmological redshift drift.
ABSTRACT
The parent paper to this Addendum describes the optimization 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 parameterization, such as summed column densities over a whole absorption complex, or common pattern relative ion abundances. This Addendum provides those analytic derivatives.
ABSTRACT
The high resolution spectrograph ESPRESSO on the VLT allows measurements of fundamental constants at unprecedented precision and hence enables tests for space–time variations predicted by ...some theories. In a series of recent papers, we developed optimal analysis procedures that expose and eliminate the subjectivity and bias in previous quasar absorption system measurements. In this paper, we analyse the ESPRESSO spectrum of the absorption system at zabs = 1.15 towards the quasar HE 0515-4414. Our goal here is not to provide a new unbiased measurement of Δα/α in this system (that will be done separately). Rather, it is to carefully examine the impact of blinding procedures applied in many previous measurements of the fine structure constant in quasar absorption systems. To do this, we emulate previous procedures, using 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 Δα/α fixed until a ‘final’ model for the absorption system is obtained, at which point Δα/α is then released as a free parameter for one final optimization. The results show that the value of Δα/α obtained in this way is systematically biased towards the initially fixed value i.e. this process produces measurements that are unrelated to the true value of Δα/α. The implication is straightforward: to avoid bias, all future measurements must include Δα/α as a free parameter from the beginning of the modelling process.