Source author record

Robert F. Carswell

Robert F. Carswell appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

17works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

17 published item(s)

preprint2025arXiv

Detecting the IGM metal enrichment with the 2-point correlation function of the flux. Application to the UVES deep spectrum

The distribution and the abundance of metals in the intergalactic medium (IGM) have strong implications on galaxy formation and evolution models. The ionic transitions of heavy elements in quasar spectra can probe both the mechanisms and the sources of chemical pollution, but high resolution, high signal-to-noise ratio (S/N) spectra are required, as the IGM absorbers can be too weak for direct detection. In this work, we investigate the IGM metallicity, focusing on the detection of the weak absorption lines. We exploited the cosmological tool of the two-point correlation function (TPCF) and applied it to the transmitted flux in the CIV forest of the ultra-high S/N UVES spectrum of the quasar HE0940-1050 (z~3). We also "deabsorbed" the strongest circum-galactic medium (CGM) systems to reveal the underlying IGM signal. For each test, we generated 1000 mock spectra with shuffled line positions to derive an estimate for the TPCF noise level. The TPCF shows a clear peak at the characteristic velocity separation of the CIV doublet, but when deabsorbing the CGM contribution (i.e. all metal lines and CIV lines associated with logN(HI)>14.0), the peak is no longer significant at 1$σ$, although 7 weak CIV systems remain visible. Even adding up to 135 additional weak mock CIV systems (logN(HI)<14.0) to the spectrum does not produce a significant CIV peak. Eventually, when we create a synthetic spectrum with gaussian distributed noise and same S/N as the complete spectrum, we remove the signal caused by the spectral intrinsic features and thus find a peak compatible with a metallicity of -3.80<[C/H]<-3.50. We conclude that the TPCF method is not sensitive to the presence of the weakest systems in the real spectrum, despite the extremely high S/N and high resolution of the data. However, the results of this statistical technique would possibly change when combining more than one line of sight.

preprint2021arXiv

A new era of fine structure constant measurements at high redshift

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.

preprint2021arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2015arXiv

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

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.

preprint2014arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

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

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.

preprint2013arXiv

Polarization Diagnostics for Cool Core Cluster Emission Lines

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.

preprint2013arXiv

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

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).

preprint2012arXiv

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

(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.

preprint2011arXiv

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

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.

preprint2010arXiv

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

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.

preprint2010arXiv

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

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.

preprint2009arXiv

Direct Evidence of Cold Gas in DLA 0812+32B

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.

preprint2008arXiv

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

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.