text:
# Size estimates for intervening C iv absorbers from high resolution spectroscopy of APM 0827+5255.

## 1 Introduction

The advent of 10-m class telescopes in the last decade has provided us with an unprecedented view of the high-redshift universe. In particular, HIRES, the high-resolution echelle spectrograph on the Keck I telescope in Hawaii, has made it possible to obtain observations with resolution as high as $`6`$ km s-1. As a result, a number of detailed studies of absorption lines in the spectra of background QSOs have been carried out (for reviews see Rauch 1998; Hamann & Ferland 1999).

However, this improvement in spectral resolution cannot be easily accompanied by an improvement in spatial resolution: 2-dimensional information about the absorbers can normally only be obtained by comparing separate spectra of multiple quasar images, which are due to either gravitationally lensed single QSOs or actual QSO pairs which have a small angular separation on the plane of the sky. Because lensed QSOs have image separations up to a few arcseconds, they are useful for probing scales $`<100`$ kpc. QSO pairs have image separations of at least a few arcminutes, thus allowing one to probe scales of a few hundred kpc. From coincident absorption in more than one line of sight (LOS), minimum size (coherence length) estimates can be obtained immediately. Estimates of most probable sizes require a statistical ‘hit and miss’ approach (e.g. McGill 1990).

Such studies have led to size estimates for Lyman $`\alpha `$ absorbers ranging from a few kpc (Foltz et al. 1984; McGuill 1990; Bechtold & Yee 1995) to a few hundred kpc (Dinshaw et al. 1995; Smette et al. 1995; Dinshaw et al. 1998) to $`1`$ Mpc (Fang et al. 1996; Dinshaw et al. 1997; Young et al. 2001). Such large ‘sizes’ suggest that the typical, low column density ($`<10^{14.5}`$cm-2) Lyman $`\alpha `$ forest line does not originate in discrete galactic or protogalactic clouds. The emerging picture is consistent with results from hydrodynamical numerical simulations in a Cold Dark Matter (CDM) structure formation scenario according to which the Lyman $`\alpha `$ forest is a manifestation of a space filling photoionized intergalactic medium (IGM). This is the most important baryon reservoir in the Universe at high redshift (see e.g. Efstathiou, Schaye & Theuns 2000 for a review).

The C iv doublet ($`\lambda \lambda 1548.20,1550.78`$ Å) is the most commonly observed heavy-element absorption signature redward of the Lyman $`\alpha `$ emission line in spectra of high-redshift QSOs. The oscillator strengths for the 1548.20 and 1550.78 Å components are in a $`2:1`$ ratio, which translates to the same optical depth ratio for optically thin ($`\tau 1`$) lines.

Line parameters are consistent with a photoionized gas (Rauch et al. 1996). The C iv absorbers appear to be highly clustered (Sargent et al. 1988; Petitjean & Bergeron 1994; Loh, Quashnock & Stein 2001) and have been thought to originate in gas clouds orbiting in galactic haloes (Sargent et al. 1979; Steidel 1993; Petitjean & Bergeron 1994; Mo & Miralda-Escudé 1996). Hydrodynamical simulations also support a picture in which high-redshift C iv arises in protogalactic clumps (PGCs), i.e., small, future constituents of large galaxies, forming through gravitational instability in a hierarchical cosmogony (Haehnelt, Steinmetz & Rauch 1996; Rauch, Haehnelt & Steinmetz 1997).

Size estimates for heavy element absorbers range from a few kpc (Young et al. 1981; Smette et al. 1992; Monier, Turnshek & Lupie 1998) to a few tens of kpc (Crotts et al. 1994; Smette et al. 1995; Lopez, Hagen & Reimers 2000) or even hundreds of kpc (Shaver & Robertson 1983; Petitjean et al. 1998). A recent study using three gravitationally lensed QSOs (Rauch, Sargent & Barlow 2001) obtains size estimates on the order of kiloparsecs. It is also apparent that sizes are absorber-type dependent, with individual clouds for low-ionization absorbers having dimensions on the order of a few $`\times 10`$ pc as opposed to at least a few $`\times 100`$ pc for high-ionization absorbers (Rauch, Sargent & Barlow 1999; Rauch et al. 2001; Rigby, Charlton & Churchill 2002).

With an $`R`$-band magnitude of 15.2 (Irwin et al. 1998) APM 0827+5255 is one of the intrinsically most luminous QSOs known (bolometric luminosity $`<10^{15}L_{}`$, Egami et al. 2000; Ibata et al. 1999). From molecular CO emission lines the estimated quasar redshift is $`z=3.911`$ (Downes et al. 1999). It is thus well suited for QSO absorption line studies. Additionally, the object is gravitationally lensed (Ledoux et al. 1998; Ibata et al. 1999; Egami et al. 2000) but the redshift of the lens is not known. As the object is lensed into three images, it represents the first confirmed odd-image lens system (Lewis et al. 2002).

For this object it is not possible to obtain separate spectra for different LOS with ground-based observations (image separation $`0.38`$ arcsec). However, if only one LOS is covered, there is an excess in the continuum intensity that gets through to the observer. If the absorbing column is high enough, the observed optical depths of the components of a heavy element doublet are not in the expected ratio, providing a signature for this effect. We use the term Anomalous Doublet Ratio (ADR) for such doublets. In such a case, in order to obtain a satisfactory fit, one has to adjust the spectral zero level locally. A number of Mg ii systems show this effect in APM 0827+5255. By means of equivalent width modelling Petitjean et al. (2000) have obtained size estimates of $`1`$ kpc for Mg ii cloudlets. From separate near-infrared spectra for each LOS Kobayashi et al. (2002) obtain upper limit estimates of $`200`$ pc for Mg ii clouds in a damped Lyman $`\alpha `$ system.

We have used the ADR to select lines from intervening C iv systems which may be due to coverage of one LOS towards APM 0827+5255. The observations and data reduction are described in Section 2, the fitting procedure in Section 3 and our partial coverage model in Section 4. We have run a set of simulations to establish whether the effect can be genuine for each fitted line, and these are described in Section 5. We have then used a maximum likelihood analysis to obtain estimates for the most probable sizes of the C iv absorbers (Section 6). We have used the currently favoured cosmological model with matter density $`\mathrm{\Omega }_M=0.3`$, and vacuum energy density $`\mathrm{\Omega }_\mathrm{\Lambda }=0.7`$ (Perlmutter et al. 1999; Riess et al. 1998; Riess et al. 2001). To facilitate comparisons with previous work, where $`q_0=0.5(\mathrm{\Omega }_M=1,\mathrm{\Omega }_\mathrm{\Lambda }=0)`$ is used almost exclusively, we have also carried out the calculations for such a model. $`H_0=72h_{72}`$ km s-1 Mpc-1 (Freedman et al. 2001) throughout. We discuss our results in Section 7 and conclude in Section 8.

## 2 Observations and data reduction

The data used in this analysis were obtained through a programme of high resolution observations on the 10-m Keck I telescope in Hawaii in April and May 1998 using HIRES with the TeK 2048$`\times `$2048 CCD at the Nasmyth focus. Details of the observations can be found in Ellison et al. (1999b). The data have been made publicly available (Ellison et al. 1999a) but we did not use the version available on the internet because there are a number of discontinuities due to inappropriate merging of the orders (also noticed by Petitjean et al. 2000). Instead, we used the original data files. However, in two of the exposures there was a spurious systematic shift in wavelength of $`0.2`$ Å, due to incorrect dating and, hence, incorrect heliocentric velocity correction. After correcting for this problem, we used Tom Barlow’s MAuna Kea Echelle Extraction programme (see http://www2.keck.hawaii.edu:3636/realpublic/inst/hires/makeewww/index.html) to obtain a set of disjoint echelle orders, wavelength calibrated and rebinned onto a linear wavelength scale of 0.04 Å per pixel.

We have flux calibrated each echelle order individually by means of the available standard stars and standard IRAF routines to ensure that there were no discontinuities where they were joined. The orders were joined by means of specially written interactive IRAF scripts which ensured a smooth joining of adjacent and overlapping parts with weights proportional to the local signal-to-noise ratio (SNR). A continuous spectrum has thus been formed for which the continuum does not change abruptly. This continuum was fitted using cubic spline functions in spectral regions deemed free of absorption.

## 3 Absorption Line Fitting

Apart from the 23 C iv systems mentioned in Ellison et al. 1999b we have fitted another 5 C iv systems. We have detected one more system which is rather uncertain and we have not used it in this analysis. Fitted redshift values range from 2.376647 to 3.671271. The latter value corresponds to a radial velocity of $`1.5\times 10^4`$ km s-1 with respect to the QSO. This reasonably ensures that the corresponding absorbers are intervening.

The lines were fitted by means of the Voigt profile fitting programme vpfit (Webb 1987; Cooke 1994; http://www.ast.cam.ac.uk/$``$rfc/vpfit.html) which uses $`\chi ^2`$ minimization to produce lists of $`b=\sqrt{2}\sigma `$ (where $`\sigma `$ is the standard deviation in a Gaussian distribution of velocities for the absorbing atoms), column density, $`\mathrm{log}N`$, and redshift, $`z`$, for all lines fitted. Where the C iv lines of interest were blended with transitions at different redshifts all lines were fitted together to produce a self-consistent fit. The wavelengths and oscillator strengths used for the C iv doublet are from Griesmann & Kling (2000). In some cases the lines of interest were blended with atmospheric transitions which were identified by means of the standard stars and an atmopheric line database. In order to obtain a satisfactory fit in the region of interest, we fitted the atmospheric transitions with Voigt profiles labelled ‘??’, which were assigned the Lyman $`\alpha `$ rest wavelength and oscillator strength. This option was also used in a few cases where it has not been possible to identify absorption features with certainty. In all, we have fitted 100 C iv doublets in the 28 systems which fall within the redshift interval mentioned above.

## 4 Model

### 4.1 Motivation

For a number of C iv doublets the fit showed a tendency to be below the data points for the 1548 Å line and above the data points for the 1550 Å line. This was an indication that the optical depth 2:1 ratio of the model profile was not adequate for the particular doublet, and it was a manifestation of an ADR, as mentioned above (see Figure 1). The fitting parameters for such lines are shown in Table LABEL:tab:adj. We stress that, short of including many extremely narrow and unphysical lines ($`b`$ instrumental resolution), the only way to compensate for the effect is to adjust the zero level. In particular, the effect cannot be attributed to continuum uncertainties in the region of these lines. The zero level for a given line which shows an ADR is effectively ‘stretched’ by the addition of a term, whose value is determined iteratively during the fitting process. This does not affect the zero level of any other lines which do not show an ADR, regardless of how close in velocity space these may be to the adjusted line. This adjustment can thus selectively compensate for excess flux due to partial coverage from some, but not all, absorbing clouds.

### 4.2 Theoretical Considerations

We now show how an ADR follows from relatively straightforward partial coverage theory.

For a given true optical depth, $`\tau _{\mathrm{true}}\frac{N}{b}`$, when there is partial coverage, the observed intensity at a given wavelength is given by

$$I=(1f)I_0+fI_0e^{\tau _{\mathrm{true}}}I_0e^{\tau _{\mathrm{obs}}}$$
(1)
where $`f\frac{\mathrm{Flux}_{\mathrm{covered}}}{\mathrm{Flux}_{\mathrm{total}}}`$ is the covering factor, $`I_0`$ is the unabsorbed continuum intensity and $`\tau _{\mathrm{obs}}`$ is the optical depth that corresponds to the observed line profile. In what follows we shall use the label $`F`$ to indicate a covering factor value obtained by means of observed and known image fluxes.

It then follows that for a doublet

$$I_1=e^{\tau _1}=1f+fe^{\tau _{\mathrm{strong}}}$$
(2)
$$I_2=e^{\tau _2}=1f+fe^{\tau _{\mathrm{weak}}}$$
(3)
where $`\tau _{\mathrm{strong}}`$ ($`\tau _{\mathrm{weak}}`$) is the true optical depth of the strong (weak) component, $`\tau _1`$, $`\tau _2`$ are the observed (apparent) optical depths and $`I_1`$, $`I_2`$ are the observed residual intensities (normalized by the local continuum, $`I_0`$) at the same velocity for the strong and weak component of the doublet, respectively.

If the lines are optically thin, then from these exact expressions the relation

$$f=\frac{1+I_2^22I_2}{1+I_12I_2}$$
(4)
can be obtained which relates the covering factor to the observed residual intensities (Hamann & Ferland 1999). In what follows we shall use the label $`R`$ to indicate a covering factor value obtained by means of measured continuum normalized residual intensities. All lines for which an adjustment of the zero level has been performed with vpfit are listed in Table LABEL:tab:adj. The zero level adjustment, determined with vpfit, as a fraction of the continuum is given in column 7 and its $`1\sigma `$ error in column 8. The covering factor determined from Equation 4, $`f_R`$, for these lines is given in column 9 of the Table. In principle, since the zero level adjustment is introduced to compensate for an excess continuum flux, the value for this adjustment should be in agreement with $`1f_R`$. The last column in Table LABEL:tab:adj shows the value of $`1f_R`$ for all entries, $`f_R`$, in column 9. In other words, the last column can be thought of as giving an ‘expected’ adjustment based on the residual intensities. If one then compares corresponding entries in columns 7 (actual adjustment used) and 10 (expected adjustment from residual intensities), one may notice that these agree to within $`1\sigma `$ for line 17, to within $`>1\sigma `$ for lines 1, 4 and 15, and to within $`2\sigma `$ or more for the rest of the lines. We call this an R-type comparison (see below) and we shall return to this point later.

Further, from Equations 2 and 3 it is clear that the ratio of the observed optical depths, $`\tau _1/\tau _2`$, depends both on $`f`$ and on the true optical depth for each component, and thus, for given $`b`$, on the column density. This theoretical result is illustrated in Figure 2 for a line with $`b=10`$ km s-1 and different values of the column density and covering factor. It can be seen that for low column densities and/or high covering factors, this ratio is close to 2. This corresponds to the case where an ADR is not observable. However, for given $`f`$ and $`b`$ there is a threshold in column density, above which the ratio will deviate significantly from 2 and, depending on the SNR in the region of the line, it will not be possible to fit the line with a model Voigt profile (which assumes a ratio of 2) without obtaining large residuals. This effect is seen in Figure 1: relative to $`\tau _2`$, $`\tau _1`$ is less than what vpfit expects and, hence, the fit appears markedly below the 1548 Å line and above the 1550 Å line.

### 4.3 Simulating an ADR

This discrepancy is due to excess flux getting through at the velocity where a particular doublet is observed. If the true optical depth and the covering factor are known, the observed, continuum normalized intensity is given by Equations 2 or 3. Thus, although the true continuum normalized profile is $`e^\tau `$, the profile actually observed will be

$$\frac{1f}{f}+e^\tau .$$
(5)
This Equation is really the same as Equation 2 (or Equation 3) but with a different normalization. In practice this means that, if one wishes to simulate this effect for various values of the true optical depth and covering factor, one needs to add a constant $`\alpha \frac{1f}{f}`$ to a continuum normalized spectrum. For example, if there are two LOS of equal brightness and one is fully covered while the other is uncovered, then $`f=0.5`$, $`\alpha =1`$ and if one starts with a given, true absorption profile in a unit continuum, the observable effects of partial coverage, if any, will be reproduced if one adds a constant $`=1`$ to all data points. The different normalization of Equation 5 is convenient because excess flux can be expressed in terms of a single additive constant, $`\alpha `$.

Thus a choice needs to be made on an appropriate value for $`f`$ (or, equivalently, $`\alpha `$) since this is necessary for simulating an ADR. In particular, in the case of APM 0827+5255 the most plausible lensing model (Egami et al. 2000, Ibata et al. 1999), as well as recent observations with the Hubble Space Telescope (Lewis et al. 2002), indicate that the lensed system is made up of three images, A, B and C, for which the flux ratio is A : B : C $`=1:0.773:0.175`$. However, image C lies closer to A, so A is less likely to be the only one covered. Therefore, it is plausible to assume that in most cases there are, effectively, two images, A+C and B, either of which, or both, may be covered by intervening absorbers. Also, since C contributes less than 10% of the total flux, its presence (or absence) will make little difference to the ADR. There is also further, more quantitative evidence that such a partial coverage configuration is plausible. This was obtained as follows:

In Table LABEL:tab:covs we have tabulated all partial coverage combinations which, from a geometric point of view, are reasonably probable. From the known multiple image fluxes we calculated the covering factors, $`f_F`$, and the corresponding ‘expected zero level adjustments’, $`1f_F`$, for each partial coverage configuration. We then compared the latter with columns 7 and 8 (actual adjustment performed, $`\pm `$ error) of Table LABEL:tab:adj (we call this a comparison of ‘type F’). This comparison suggests that, among the simple partial coverage models in Table LABEL:tab:covs, the one with $`\alpha =1.520`$ ($`f0.4`$) provides the best physical basis for most of the adjustments presented in Table LABEL:tab:adj, except those whose index number is marked by an x. Alternatively, it is possible to compare the last column of Table LABEL:tab:adj (‘expected zero level adjustment’, $`1f_R`$, from the residual intensities) to columns 7 and 8 of the same Table (we call this a comparison of ‘type R’). However, for any given line, $`f_R`$ is likely to be less reliable than $`f_F`$ as it depends on the residual intensity measured, which is strongly affected by blending due to other nearby lines. We have thus not based our choice of model on an R-type comparison. Nevertheless, it is interesting to note that for line x15 which is only moderately blended, no simple model among those shown in Table LABEL:tab:covs comes close to being acceptable, on the basis of an F-comparison. However, this line is only moderately blended, and an R-comparison shows that corresponding values in columns 7 and 10 of Table LABEL:tab:adj are fairly close. There is thus little doubt about the reality of the effect, but, perhaps, there is a more complicated physical situation which requires more sophisticated modelling. On the other hand, model ‘A covered, B+C uncovered’ would be adequate for lines x28 and x29 (on the basis of our chosen F-comparison). In any case, we chose not to discard lines for which there was poor agreement between columns 7 and 10 of Table LABEL:tab:adj (the R-comparison failed) as they might provide useful insight after the simulations were complete.

To summarise, by performing Voigt profile fitting of intervening C iv lines, we have obtained a sample of 100 fitted $`b`$$`\mathrm{log}N`$ pairs. For 12 of these lines an ADR was present and an adjustment of the zero level has been applied in order that a good fit be obtained. We have chosen a set of 4 possible partial coverage configurations for this object and, using the known flux ratios of the three images, have obtained 4 values for $`1f_F`$, the ‘expected zero level adjustment’. One of these values agrees best with the greatest number of adjustments actually imposed on the 12 fitted lines. We thus chose the corresponding model ($`f0.4`$) as the best candidate for a physical explanation of the observed ADR.

## 5 Simulations

### 5.1 Need for Simulations

A reliable way of establishing that the chosen partial coverage model can (or cannot) give rise to a detectable ADR for a given $`b`$$`\mathrm{log}N`$–SNR combination is to perform a set of simulations, as explained in the next subsection.

Once the detectability, or not, of the ADR has been established, this can be compared to the actual situation in the observed spectrum. For each $`b`$$`\mathrm{log}N`$–SNR combination there are four possibilities:

1. ADR not detectable and not seen in observed spectrum,
2. ADR detectable and seen in observed spectrum,
3. ADR detectable but not seen in observed spectrum,
4. ADR not detectable but seen in observed spectrum.

Cases (ii) and (iii) are useful for estimating absorber sizes: In case (ii) partial coverage is established and, therefore, the transverse size of the absorber must be smaller than the proper separation of the LOS at the absorbing redshift. Conversely, in case (iii) total coverage is established and the transverse size of the absorber must be larger than the proper separation of the LOS.

Further, a maximum likelihood approach can be used to estimate most probable absorber sizes from the number of established partial and total coverage cases.

### 5.2 Method

Each of the 100 fitted $`b`$$`\mathrm{log}N`$ pairs mentioned above gives a single true $`\tau `$ value and can be used to produce an artificial continuum normalized absorption profile, $`e^\tau `$. If one adds a term $`\alpha =1.520`$ to each data point in this profile, one effectively introduces excess flux due to partial coverage from our chosen physical model (Equation 5).

Because the SNR varies considerably over the observed spectrum from $`25`$ to $`120`$, for each region containing a C iv doublet component an average SNR value was determined by inspection of the spectrum using the task splot in IRAF.

For each $`b`$$`\mathrm{log}N`$(C iv) pair two artificial spectra were produced. The first spectrum contained 120 C iv doublets with a profile $`e^\tau `$. The spacing of the lines was chosen empirically so that when subsequently the lines would be fitted with vpfit, blending complications would be absent, i.e. even the broadest lines were well apart from each other. The spectrum was resampled to the resolution of the observed spectrum (0.04 Å per pixel) and smoothed to the instrumental resolution of 6.6 km s-1. Gaussian noise was added with $`\sigma =\frac{1}{\mathrm{SNR}}`$ to obtain the desirable final SNR, as measured beforehand for each region.

The above process was repeated to produce another set of 120 simulations of the same line aiming to reproduce the effects of excess flux. This time gaussian noise with $`\sigma =\frac{1+\alpha }{\mathrm{SNR}}`$ was added to the unit continuum spectrum and both data and continuum were then increased by $`\alpha `$. This ensured that the desired SNR was obtained in this case as well.

Next, vpfit was used to fit the artificial spectra. Each of the 120 artificial doublets in the spectrum was fitted independently, so that 120 separate $`\chi ^2`$ values were obtained. The starting $`b`$$`\mathrm{log}N`$ values for the vpfit iterations were exactly the same as those used to make the simulated lines, so that one would expect to obtain a good fit, unless the effects of excess flux were significant.

vpfit was allowed to fit only a single Voigt profile to each simulated line. The option of automatically putting in additional narrow lines to improve poor fits was not used since such lines were known to be spurious by construction. The fitting region for each doublet component was centred on the component ‘redshift’ and extended over $`3\times \mathrm{FWHM}`$, where $`\mathrm{FWHM}=2\sqrt{\mathrm{ln}2}b`$ is the full width at half minimum. This was done to avoid including too much continuum which can make a poor fit less pronounced.

If for a given $`b`$$`\mathrm{log}N`$ pair the higher continuum introduced an observable ADR, then, on average, the second set of 120 fits should be markedly poorer than the first set. Otherwise, the two distributions of $`\chi ^2`$(reduced) obtained should not differ significantly. Figure 4 shows a single simulation of line 3 with the Voigt profile fit superposed. An ADR was reproduced at $`>2\sigma `$ level in the central region of the line. Note that a $`2:1`$ doublet ratio was recovered for this simulated line when the zero level was adjusted by $`0.62\pm 0.02`$, which agrees both with the original adjustment of $`0.54\pm 0.11`$ (Table LABEL:tab:adj) and with $`1f_F=0.603`$ for a simulated line with $`\alpha =1.520`$ (Table LABEL:tab:covs). Independently from this result from vpfit, the continuum normalized residual intensities gave a covering factor $`f_R=0.378`$ ($`1f_R=0.622`$), also in good agreement. Additionally, the original Voigt parameters of the line were recovered after the adjustment of the zero level ($`b=5.0\pm 0.2`$ km s-1, $`\mathrm{log}N=13.44\pm 0.05`$; compare with Table LABEL:tab:adj). This shows that this type of adjustment successfully compensates for partial coverage effects and that using fitted $`b`$$`\mathrm{log}N`$ values to simulate lines with an adjusted zero level is physically meaningful.

For each simulated $`b`$$`\mathrm{log}N`$ pair we performed a two-sided Kolmogorov–Smirnov (KS) test for the two cumulative distributions, $`C_1`$ and $`C_{1+\alpha }`$ of $`\chi ^2`$(reduced) obtained, in order to see if the ADR had been reproduced, in a statistically significant sense, by the simulations. If the significance level for the KS statistic, defined as $`D\mathrm{max}C_1(\chi ^2)C_{1+\alpha }(\chi ^2)`$, was less than 0.05, the null hypothesis that the two distributions were the same was discarded. This was interpreted as evidence that the second set of 120 fits was poorer and thus an ADR had been reproduced. The significance results for all such lines are given in Table LABEL:tab:sim0.4. We illustrate by plotting the cumulative distributions of the $`\chi ^2`$ values for fits to simulations of line 3 in Figure 3.

If the corresponding originally fitted line also showed an ADR \[case (ii) in Section 5.1\], the assumed partial covering model was taken to be the physical reason (a p-case. In work where separate spectra are available for each LOS the term ‘anticoincidence’ is commonly used for such cases). On the other hand, if the effect appeared only in a simulation \[case (iii) in Section 5.1\], then, as far as this model was concerned, there was no partial but total coverage (a t-case or ‘coincidence’ in work with separate spectra).

In Section 6 we explain how we used the number of t and p-cases to obtain an estimate for the absorber sizes for two distinct absorber geometries.

### 5.3 Results

According to the reasoning above, for $`\alpha =1.520`$ Table LABEL:tab:sim0.4 establishes the detectability of the ADR for the 10 simulated lines listed there, of which 6 did not show the effect in the observed spectrum (and were thus taken to be t-cases) and 4 also showed the effect in the observed spectrum (p-cases).

More generally, all fitted and simulated lines can be divided into the four groups introduced in Section 5.1:

1. For the majority of lines no ADR was detected either during Voigt profile fitting or in simulations. Such lines were not used further (but see discussion).
2. For lines 3, 4, 15 and 17 the zero level has been adjusted during fitting and the ADR was reproduced in simulations. These p lines are indicated by solid ticks in Figures 5 and 6. These lines belong to each of the two systems ($`z`$ 3.10 and $`z`$ 3.13) that show the effect and were used in this analysis (see Table LABEL:tab:adj).
3. For lines 7, 18, 25, 26, 45 and 96 an ADR was only produced in simulations with varying degrees of significance. Lines 7 and 18 are from the two systems at $`z`$ 3.10 and $`z`$ 3.13 for other lines of which partial coverage was suggested by the simulations. Lines 25 and 26 belong to the $`z`$ 3.17 system in which line x15 shows an ADR which cannot be modelled with $`\alpha =1.520`$. Lines 45 and 96 are from lower redshift systems. These t lines are indicated by dot-short-dashed ticks in Figures 5 to 9.
4. Finally, for lines 1, 2 ($`z`$ 3.10 system) and 13, 14, 16 ($`z`$ 3.13 system) an ADR was not reproduced in simulations although they have been fitted with a zero level adjustment. These lines are indicated by long-dashed ticks and numbers in Figures 5 and 6 and were not used further.

The existence of this last group of lines in a way provides an a posteriori justification of simulations. To check the effect that zero level adjustment of such lines had on the overall goodness-of-fit, we modified the fit for the $`z`$ 3.10 and 3.13 system as follows:

We used the last fitting results for each region as input to vpfit but excluded group (iv) lines from zero level adjustment. Otherwise, we did not introduce any new lines. This procedure had a negligible effect on the fit for the $`z`$ 3.10 system (reduced $`\chi ^2`$ changed from $`1.09`$ to $`1.10`$). The original adjustment of the zero level for lines 1 and 2 was not necessary and was probably a result of the empirical nature of this adjustment combined with the fact that these lines are close to the stronger line 3. In the case of the $`z`$ 3.13 system and lines 13, 14 and 16 the situation is similar since the reduced $`\chi ^2`$ changed from $`1.16`$ to $`1.17`$.

The simulations thus allowed us to distinguish for which lines an ADR may have a physical cause.

The choice of an upper limit for the significance of the KS statistic was somewhat arbitrary. We chose 5% in an effort not to discard too many lines, since few lines were available anyway. However, the $`\chi ^2`$ cumulative distribution plots for lines 25 and 26 strongly indicate that a lower value might be appropriate since $`\chi ^2`$ values for fits of simulated lines with excess flux did not show a clear trend to be higher. Since we were mainly interested in illustrating a technique, we formally accepted the 5% level and used these two lines in further analysis. If this caveat is taken into account, Table LABEL:tab:sim0.4 suggests that a threshold column density for an observable ADR in this model is given by $`\mathrm{log}N`$(C iv$`>13.1`$ with $`b10`$ km s-1. This explains the failure to reproduce lines in group (iv) above since all of these lines (except, perhaps, for line 2, which, however, has large $`\sigma _{\mathrm{log}N}`$) are below this threshold.

## 6 Estimating Absorber Sizes

### 6.1 Direct Estimates

Miminum size estimates can be obtained straightforwardly from lines which correspond to total coverage.

The proper separation at reshift $`z`$ of the light paths from the two images of a lensed QSO with emission redshift $`z_{\mathrm{em}}`$ due to a lens at redshift $`z_\mathrm{l}`$ is given by

$$S(z)=\frac{D(z_\mathrm{l},z_{\mathrm{con}}=0)D(z,z_{\mathrm{em}})}{D(z_\mathrm{l},z_{\mathrm{em}})}\delta \varphi ,$$
(6)
(Young et al. 1981) where $`\delta \varphi `$ is the observed angular separation of the images. The light paths converge at $`z_{\mathrm{con}}=0`$ and the angular diameter distances, $`D`$, are given by

$$D(z_1,z_2)=\frac{D_C(z_2)D_C(z_1)}{1+z_2}$$
(7)
where

$$D_C(z)=\frac{c}{H_0}_0^z\frac{dz^{}}{\sqrt{\mathrm{\Omega }_M(1+z^{})^3+\mathrm{\Omega }_\mathrm{\Lambda }}}$$
(8)
is the comoving distance for a flat universe ($`\mathrm{\Omega }_K=0`$; Peebles 1993; Hogg 2000). Equ. 6 holds for $`z>z_\mathrm{l}`$ as assumed here. Egami et al. (2000) argue that the lensing galaxy should be at $`z_\mathrm{l}3`$ based on the Einstein and core radii of their model. The damped Lyman $`\alpha `$ system at $`z2.974`$ may be compatible with this estimate. Additionally, Petitjean et al. (2000) suggest that strong Mg ii systems are prime lense candidates because the roughly equal brightnesses of images A and B suggest the LOS are traversing the central regions of the lensing object. We have thus calculated the separation for five values of $`z_\mathrm{l}`$ as shown in Tables LABEL:tab:sep and LABEL:tab:sep2. These values include an arbitrary value of 0.7, not corresponding to any known object, three values corresponding to observed Mg ii systems and one value from the damped Lyman $`\alpha `$ system mentioned. The observed angular separation of image B from the combined A+C image was taken to be $`\delta \varphi =0.369`$ arcsec. There is some ambiguity as to what one means by the location of the A+C image, and we have simply chosen the location of the flux-weighted centroid of the A+C image. Size estimates scale with $`\delta \varphi `$ (this still holds, approximately, in the statistical approach of the next section) and so the effect is of the order of no more than a few percent. Thus the results remain qualitatively unaffected.

### 6.2 Maximum Likelihood Analysis

We then applied a maximum likelihood method in order to estimate the most probable size of the absorbers, given the information that we have obtained through this work. We assumed two simple geometries for absorbers of uniform size.

The probability that a spherical cloud is intersected by both LOS, given that it is intersected by one is

$$\varphi _s(X)=\frac{2}{\pi }\left\{\mathrm{arccos}\left[X(z)\right]X(z)\sqrt{1X(z)^2}\right\}$$
(9)
(McGill 1990), for $`X[0,1]`$, and zero otherwise. Here $`X(z)S(z)/2R`$, where $`R`$ is the absorber radius. For randomly inclined circular disks the probability is

$`\varphi _d(X)={\displaystyle _{\pi /2}^{\pi /2}}{\displaystyle \frac{\mathrm{cos}\theta }{\pi }}`$ (10)
$`\times \left\{\mathrm{arccos}\left[{\displaystyle \frac{X(z)}{\mathrm{cos}\theta }}\right]{\displaystyle \frac{X(z)}{\mathrm{cos}\theta }}\sqrt{1{\displaystyle \frac{X(z)^2}{\mathrm{cos}^2\theta }}}\right\}d\theta `$
for $`X<\mathrm{cos}\theta `$ and zero otherwise (McGill 1990). Here $`\theta `$ is the angle of inclination of the disk relative to the lines of sight. We were actually interested in the case where there is absorption at either LOS, so that the probability that there is absorption at the other LOS is

$$\psi _{s,d}=\frac{\varphi _{s,d}}{2\varphi _{s,d}}$$
(11)
\[by analogy to Dinshaw et al. (1997) who apply this to two separate spectra\]. The probability of getting the estimated number of cases where there is coverage of one or both lines of sight is then given by the likelihood function

$$(R)=\underset{i}{}\psi _{s,d}\left[X(z_i)\right]\underset{j}{}\left\{1\psi _{s,d}\left[X(z_j)\right]\right\},$$
(12)
where $`i`$ labels cases of total coverage (t lines) and $`j`$ cases of single coverage (p lines).

The results of the likelihood analysis for different lensing redshifts are summarized in Figure 10 and Tables LABEL:tab:lik and LABEL:tab:lik2. The plots show the differential likelihood distributions for the two assumed geometries, $`(R)`$, which peak at the most probable radii for given $`z_\mathrm{l}`$. Also plotted are the cumulative distributions from which median $`R`$ values have been calculated together with 95% confidence intervals. For spherical absorbers in a non-zero $`\mathrm{\Lambda }`$ universe the most likely sizes ($`2R`$) are in the range from 1.74 to 9.26 $`h_{72}^1`$ kpc. Median sizes are in the range from 2.6 to 13.9 $`h_{72}^1`$ kpc. For randomly inclined disks the most likely sizes lie between 2.56 and 13.58 $`h_{72}^1`$ kpc. Median sizes are between 3.92 and 20.84 $`h_{72}^1`$ kpc. In a zero $`\mathrm{\Lambda }`$ universe these values are smaller by a factor of $`1.4`$.

## 7 Discussion

We have used a simple model to select candidate lines for which the ADR is due to coverage of one LOS. Our model only considered two LOS, the corresponding images to which, nevertheless, are not equally bright. As in previous work on the subject, only two simple geometries have been considered for the absorbing structures. These geometric assumptions are simple not only in terms of the assumed shapes but also in that these structures more realistically should have no sharp edge: even in the simplest scenario the column density has a profile which falls smoothly radially outwards. Thus the absorbers giving rise to lines x15, x28 and x29 may conceivably have the same column density profile as the rest but they may be positioned with respect to a LOS in such a way that they present a different covering factor. However, incorporation of such more realistic models is only possible when separate spectra for different LOS are available (e.g. D’ Odorico et al. 1998, Rauch et al. 2001).

The use of simulated spectra to reproduce an ADR is crucial because the observed ADR by itself does not provide enough information to establish the partial coverage effect. For example, group (iv) lines did not show an ADR when simulated. Also the magnitude of the adjustment does not necessarily constrain the model sufficiently. As explained, in theory one would expect the adjustment to agree with $`1f_R`$ calculated from the residual intensities (columns 7 and 10 in Table LABEL:tab:adj). In practice, this was not always the case. This was probably due to blending with nearby lines. Blending affects the observed residual intensities at the velocity position of the line of interest, leading to an inaccurate value for $`1f_R`$. It may also have an effect on the magnitude of the adjustment performed, leading to further discrepancy between the values in columns 7 and 10. For example, line 3 (adjustment = $`0.54\pm 11`$, $`1f_R=0.211566`$) was simulated with an ADR but lies 5.7 km s-1 from line 2 ($`b=12.8`$ km s-1, $`\mathrm{log}N=13.07`$) and 13.6 km s-1 from line 7 ($`b=18.6`$ km s-1, $`\mathrm{log}N=13.53`$). Although when simulated and fitted the line was in excellent agreement with the $`\alpha =1.520`$ model, for the observed and fitted line there was some discrepancy (though not within the error) between actual adjustment and $`1f_F`$, and a more serious discrepancy between actual adjustment and $`1f_R`$. Of course some degree of agreement is necessary, at least between $`1f_F`$ and actual adjustment. As explained before, this is how we chose the best value of $`\alpha `$ for our simulations.

On the other hand, line 17 does not suffer from as much blending, is stronger and the discrepancy is smaller.

More insight on the threshold Voigt parameters a line needs to have for partial coverage to lead to an observable ADR, and on the reality of the effect can be obtained by inspection of group (i) lines. Specifically, such lines should be below the estimated threshold mentioned above (Section 5.3), in that one or more of the following are true: The lines have low $`\mathrm{log}N`$, large $`b`$, or they lie in a region with low SNR. Parameters for all group (i) lines with $`\mathrm{log}N`$(C iv$`13.0`$ are shown in Table LABEL:tab:group1. It can be seen that usually lines are too broad for their column density. For line 93 the problem may be aggravated by a relatively low SNR.

Given the small size of the sample, these results cannot be considered representative of intervening C iv absorbers in general. Rather they should be considered as an illustration of what can be learned from a compilation of similar studies of gravitationally lensed QSOs. From Equations 9 to 12 it follows that these results depend on LOS separations for the sample lines and on lensing redshifts chosen. However, taken at face value they strongly suggest that C iv absorption arises in structures with sizes on the order of kiloparsecs. It can be seen from Tables LABEL:tab:lik and LABEL:tab:lik2 that, for the three $`z_\mathrm{l}`$ values 1.062, 1.1727 and 1.181, estimated sizes are of the order of $`2`$ kpc. On the other hand, Tables LABEL:tab:sep and LABEL:tab:sep2 suggest that the structures responsible for intervening C iv absorption cannot be smaller than $`0.2h_{72}^1`$ kpc for $`\mathrm{\Omega }_M=1`$, $`\mathrm{\Omega }_\mathrm{\Lambda }=0`$ ($`0.3h_{72}^1`$ kpc for $`\mathrm{\Omega }_M=0.3`$, $`\mathrm{\Omega }_\mathrm{\Lambda }=0.7`$).

Up to this point we have implicitly assumed that the estimated ‘sizes’ are those of coherent, discrete entities. In particular, the maximum likelihood method is based on an assumption about the geometric shape of such objects. However, whether one has separate spectra for different LOS or not, the same observations can also be due to an ensemble of correlated objects of smaller characteristic size. In this case total coverage provides an estimate for a coherence length of clustered absorbers. Alternatively, there might be a lack of discrete objects at any scale, with partial coverage being the effect of variations in an inhomogeneous density field permeating the IGM.

If total coverage is due to correlated discrete absorbers, then in practice we may be ‘measuring’ transverse sizes of clumpy systems. A way to investigate this issue is by observing systems which show several lines above the ADR detectability threshold and close in velocity space. If these proved to be p-type lines, we would probably be dealing with a clumpy system covering one LOS. We claim that Figure 5 argues against such a picture for the $`z`$ 3.10 system since line 7 (t) lies between lines 3 and 4 (p).

There are two main caveats concerning our results. First, the sample of lines is small and this is expected to affect the estimates themselves because it is not just the fraction of hits and misses that is important (Fang et al. 1996). The most probable size calculated reflects the LOS separation over a small redshift range. It is unclear whether the result would remain similar if a broader redshift range were covered. If such a sample were available, it would be possible to carry out size estimates in different redshift bins so as to look for any evolutionary trends in the estimated absorber sizes. Evolution of the UV background field would be expected to have an effect on absorber sizes. If these are predominantly photoionized, we might expect size estimates for C iv absorbers to become smaller at low $`z`$. Note however that apparent redshift evolution may be an artefact due to LOS separation (Fang et al. 1996; Dinshaw et al. 1998) and, therefore, difficult to confirm.

Second, the redshift of the lens is unknown; multiple lensing would introduce further complications. Thus these results only give an order of magnitude estimate of C iv absorber size.

It is worth noting that these results agree in order of magnitude, at least for the lower $`z_\mathrm{l}`$ values, with recent results from work on a larger sample with separate spectra for different LOS (Rauch et al. 2001). These authors found that there is very little difference between column densities in separate spectra for LOS separations below about $`300h_{50}^1`$ pc. On the other hand, column density differences reach 50% just below LOS separations of $`1h_{50}^1`$ kpc. From energy arguments these authors conclude that C iv clouds are the end products of ancient galactic outflows. There is also independent evidence for such a scenario. Theuns et al. (2002) showed that hydrodynamical simulations with feedback produce C iv absorption lines with properties in reasonable agreement with observations. They used the same Voigt profile fits to C iv lines which we have used in the present paper, together with fits to one more QSO to demonstrate this result. It is thus possible that the structures examined here represent clouds whose baryonic content has been transported, or at least metal-enriched, to large distances (up to $`1`$ Mpc) from the central regions of their parent galaxies.

Our results also agree with the results of Petitjean et al. (2000) for Mg ii systems. These authors used equivalent width models in which one LOS was saturated and one optically thin to estimate that $`\mathrm{log}N`$(Mg ii) for $`z<1.7`$ absorbers changes by at least an order of magnitude over less than $`1h_{75}^1`$ kpc. However, there is no indication that the two populations may belong to the same class, as there are no low ionization lines associated to any of our C iv lines which show an ADR. Although there are some low ionization lines associated to some of the other C iv lines in the full sample, these do not show multiplets, so that it is not possible to look for evidence of partial coverage. There is also no overlap in redshift between the two populations.

On the other hand, if our somewhat higher estimates for $`z_\mathrm{l}=2.974`$ are closer to reality, then these would support the picture according to which intervening C iv absorbers arise in PGCs (Haehnelt et al. 1996, Rauch et al. 1997).

Note that in total there are 8 (10) components which are either p or t, depending on whether we count lines 25 and 26. 4 out of these are p which is a fraction of 0.5 (0.4). Of course, this is a single case and its statistical significance cannot be estimated. However, even if, in a worst case scenario, we trust only one (observed) line, this result says that the probability of partial coverage for lines above the detectability limit is not negligible. This in turn suggests that for gravitationally lensed QSOs partial coverage by intervening structures is not rare and observing it is mainly an issue of detectability.

## 8 Conclusions

The spectra used in this work are the result of careful wavelength and flux calibration, and order merging. The observed ADRs are thus real, and we have shown that in many cases they can be simulated by means of a simple physical model. From LOS separations at different redshifts we obtained lower limit estimates of $`300h_{72}^1`$ pc for C iv absorbers. From a statistical approach we obtained most probable sizes of a few kiloparsecs, but, because the sample is small, this estimate is rather uncertain.

We have briefly discussed the possible nature of the absorbers based on our size estimates. Clearly, more work needs to be done before we can decide among different possibilities. In this respect it would be interesting to carry out a cross-correlation study between Lyman break galaxies and our C iv absorbers similar to the work by Adelberger et al. (2002) who found that the two classes may be related.

Ideally, one would like to have separate spectra for different LOS. However, given that lensed QSOs may have small image separations, it is difficult to find suitable candidates. The technique used here compensates for this problem as it allows probing of closer LOS. All that is required is the existence of two images with a known angular separation as well as an estimate for the redshift of the lens. As in the case studied here, the separation may be obtained in the course of different observations with different instruments (e.g. space based) or in different wavebands (e.g. radio) where the necessary resolution can be attained. Obtaining a reliable value for $`z_l`$ is more difficult as it depends on a combination of modelling and detection of the lensing mass, but if there are several candidates, several models can be run.

In principle the method could be extended to more than two images. However, with several images of different brightnesses the uncertainties in the zero level adjustment parameter would probably allow different configurations. It may be possible to add a probability argument for these as well, but adding this degree of complexity for what would be a small sample is not likely to yield worthwhile results.

Finally, it is worth noting that ADRs can be useful in a different way as well. In cases of spectra of QSOs which are not known to be lensed, the need for same, or similar, ADR corrections for several absorbers would be a strong indication that a QSO is lensed. The method illustrated in this paper would then have to be modified to include the angular separation as an unknown as well, with more speculative results.

## Acknowledgments

The observations were made at the W. M. Keck Observatory, which is operated as a scientific partnership between the California Institute of Technology and the University of California; it was made possible by the generous support of the W. M. Keck Foundation.