arxiv: v1 [astro-ph.ga] 25 Aug 2016

Transcription

1 Preprint 26 August 2016 Compiled using MNRAS LATEX style file v3.0 Narrow C iv absorption doublets on quasar spectra of the Baryon Oscillation Spectroscopic Survey arxiv: v1 [astro-ph.ga] 25 Aug 2016 Zhi-Fu Chen 1,2,3,4, Qiu-Sheng Gu 1,3,4,, Luwenjia Zhou 1,3,4, and Yan-Mei Chen 1,3,4 1 School of Astronomy and Space Science, Nanjing University, Nanjing , P. R. China 2 Department of Physics and Telecommunication Engineering, Baise University, Baise , China 3 Key Laboratory of Modern Astronomy and Astrophysics, Nanjing University, Nanjing , P. R. China 4 Collaborative Innovation Center of Modern Astronomy and Space Exploration, Nanjing, , P.R. China 26 August INTRODUCTION The circumgalactic medium (CGM) surrounding galaxies is crucial to comprehending the feedback from and gas accretion into central objects, as well as galaxy evolution. In recent years, there are a number of successful results of directly detecting emission lines arising from CGM by deep imaging or spectroscopy observations (Martin et al. 2014; Hennawi & Prochaska 2013; Cantalupo et al. 2014; Husband et al. 2015; Hennawi et al. 2015; Martin et al. 2015; Arrigoni Battaia et al. 2016). However, it is still difficult to systematically survey the emission of CGM due to the low gas density. Luckily, the CGM would leave marks (so called absorption features) on the spectra of the background objects, which is an efficient way to investigate the properties and kinematics of CGM. The gas, which infalls from the intergalatic medium (IGM) into galaxies to fuel star formation and is evaporated from galaxies to enrich the metal abundance of zhichenfu@126.com ABSTRACT In this paper, we extend our works of Papers I and II (Chen et al. 2014a,b), which are assigned to systematically survey C ivλλ1548,1551 narrow absorption lines (NALs) with z abs z em on quasar spectra of the Baryon Oscillation Spectroscopic Survey (BOSS), to collect Civ NALs with z abs z em from blue to red wings of Civλ1549 emission lines. Together with Papers I and II, we have collected a total number of 41,479 Civ NALs with z abs in surveyed spectral region redward of Ly α until red wing of Civλ1549 emission line. We find that the stronger Civ NALs tend to be the more saturated absorptions, and associatedsystems (z abs z em ) seem to have larger absorption strengths when compared to intervening ones (z abs z em ). The redshift density evolution behavior of absorbers (the number of absorbers per redshift path) is similar to the history of the cosmic star formation. When compared to the quasar-frame velocity (β) distribution of Mg ii absorbers, the β distribution of C iv absorbers is broader at β 0, shows longer extended tail, and exhibits a larger dispersion for environmental absorptions. In addition, for associated C iv absorbers, we find that low-luminosity quasars seem to exhibit smaller β and stronger absorptions when compared to high-luminosity quasars. Key words: galaxies: active quasars: general quasars: absorption lines IGM, can be well characterized by metal absorption lines (Springel & Hernquist 2003; Dekel et al. 2009; Richter 2012; Davé, Finlator, & Oppenheimer 2012; Rubin et al. 2012; Bouché et al. 2013; Ford et al. 2014; Borthakur et al. 2015; Kacprzak et al. 2015; Nielsen et al. 2015; Muzahid et al. 2015, e.g.,). The cold gas with T K of CGM has been explored by various transitions from neutral hydrogen to high-ionization absorption lines, such as the absorptions of Ly α (e.g., Kollmeier et al. 2006; Balashev et al. 2014), Ca iiλλ3934,3969 (e.g., Zhu & Ménard 2013; Murga et al. 2015), Mg iiλλ2798,2803 (e.g., Kacprzak et al. 2011; Farina et al. 2014; Huang et al. 2016), Siiiiλ1206 (e.g., Keeney et al. 2013), Si ivλλ1393,1402 (e.g., Songaila 2005; Cooksey et al. 2011), C ivλλ1548,1551 (e.g., Cen & Chisari 2011; Burchett et al. 2013), N vλλ1239,1242 (e.g., Fox et al. 2007, 2009; Fechner & Richter 2009), and O viλλ1031,1037 (e.g., Tripp et al. 2008, 2011; Savage et al. 2014). The properties of absorbers (e.g., ionization state, strength, kinematics) are related to halo mass of galaxy, locations and distances between the center of galaxy and absorber, colors, types, luminosities, morphologies, inclination angles, c 0000 The Authors

2 2 Zhi-Fu Chen et al. and star formation rates (e.g., Kacprzak et al. 2007, 2015; Gauthier & Chen 2011; Ford et al. 2014; Farina et al. 2014; Muzahid et al. 2015). Thus, absorption lines provide key insight into galaxy s CGMs, feedback from central objects, evolutions as well as intrinsic properties of host galaxies. Outflows and/or jets are fundamental components of active galactic nucleus (AGNs), which play a role in the galaxy evolution by the form of feedback. They carry angular momentum, energy and huge amounts of material away from production sites into host galaxies, CGMs and even IGMs to regulate their star formation rate, metallicities, and element distributions. The outflows are often detected via blueshifted absorption lines. Broad absorption lines (BALs) imprinted on quasar spectra are undoubtedly intrinsic to quasar systems and are likely to be formed in quasar outflows. Narrow absorption lines (NALs) with line widths of a few hundreds kms 1 are very common on quasar spectra. Some previous studies (e.g., Lü, Qin, & Gupta 2007; Vanden Berk et al. 2008; Shen & Ménard 2012; Pan & Chen 2013) have revealed that a part of NALs with z abs z em is also intrinsic to quasar systems. Variation analysis of NALs suggests that some of NALs are probably produced in quasar outflows (e.g., Wise et al. 2004; Narayanan et al. 2004; Chartas et al. 2009; Hamann et al. 2011; Chartas et al. 2012; Chen et al. 2013, 2015). And also, NALs are more common than BALs on quasar spectra especially for Mg ii, because i) low ionization BALs (LoBALs) are rare, and ii) a significant fraction of Mgii NALs are formed in intervening absorbers (i.e., not physically related to quasar themselves). Therefore, taking full advantage of NALs is crucial to understanding quasar outflows and their influences on surrounding environment. The Sloan Digital Sky Survey (SDSS, York et al. 2000) is a widely ambitious and influential survey in the astronomical history, which used a dedicated 2.5-meter telescope (Gunn et al. 2006) at Apache Point Observatory, New Mexico, to obtain optical spectra with coverage from λ = Å at a resolution of R 2000 (Abazajian et al. 2009) during the first (SDSS-I, ) and second (SDSS-II, ) stages. Some groups (e.g., Nestor, Turnshek, & Rao 2005; Quider et al. 2011; Zhu & Ménard 2013; Cooksey et al. 2013; Seyffert et al. 2013) have systematically survey NALs on quasar spectra of the SDSS-I/II. However, most of them only focused on intervening NALs (z abs z em), which limits the utilization of NALs to study quasar outflows and their influence on surrounding environment. The SDSS program continued with the Third Sloan Digital Sky Survey (SDSS-III) using the same 2.5-meter telescope with update spectrographs (Smee et al. 2013), which collected data from July 2008 to June 2014 and obtained spectra with wavelength range from 3600Å to 10400Å at a resolution of R = (Eisenstein et al. 2011). The Baryon Oscillation Spectroscopic Survey (BOSS, Dawson et al. 2013) is the main dark time of legacy survey of the SDSS-III. The first spectral data release of the BOSS (SDSS DR9) includes 87,822 quasars (DR9Q, Pâris et al. 2012). On quasar spectra of the DR9Q, our early systematical surveys (Chen et al. 2014a,b, here after Papers I and II) of Civ NALs are focused on intervening NALs, which are similar to the surveys mentioned in the last paragraph. In order to systematically collect Civ NALs with z abs z em, we extend our program of Papers I and II to spectral data redward of C ivλ1549 emission lines. These NALs will favour future studies on AGN feedback and its influence on central objects and surrounding environment. The paper is organized as follows. Section 2 describes quasar sample selection, definition of surveyed spectral region, spectral analysis and NAL selection. We calculate the redshift path and estimate the completeness of our absorption line detection in Section 3. The strength and saturability of absorptions are characterized in Section 4. Section 5 is discussions, which statistically describes redshift density evolution of absorbers, β (see Equation 8 for the definition of β) distributions, dependence of β and absorption strengths on radio detection and UV continuum luminosity. A short summary is presented in Section 6. Throughout this paper, we adopt cosmological parameters of Ω M = 0.3, Ω Λ = 0.7, and H 0 = 70 km s 1 Mpc 1. 2 QUASAR SAMPLE AND SPECTRAL ANALYSIS 2.1 Quasar sample and surveyed spectral region In this paper, we extend our previous search of Civ NALs released in Papers I and II to spectral data redward of C ivλ1549 emission lines. Therefore, the quasar sample and selection criteria is the combination of Papers I and II. We briefly summarize the criteria used to construct quasar sample and surveyed spectral region as follows: (i) Surveyed spectral region. Noisy data shortward of 3800Å at observed frame, and the contamination of series of absorptions, for example, doublets of [O i]λ1302 and Siiiλ1304, and Ly α forest, lead to difficulties in searching for Civ NALs. Thus, spectral region is limited to the data longward of the maximum value of 3800Å and 1310 (1+z em)å at observed frame. Papers I and II mainly paid attention to intervening NALs. Therefore, another limit is adopted to constrain the data blueward of 10,000 kms 1 of C ivλ1549 emission lines. (ii) Signal-to-noise ratio per pixel ( S/N ). The data accompanied with large noisy would markedly reduce the surveyed efficiency of NALs and may greatly lead to misidentification. In surveyed spectral region, the S/N of all quasar spectra in the DR9Q has a median value of about 4 pixel 1. The cut S/N 4 pixel 1 is accepted to guarantee robust absorption line measurements. In terms of above limits, Papers I and II selected 37,241 quasars with 1.54 z em 5.16 from the DR9Qthat contains 87,822 quasars. Using the spectra of these 37,241 quasars whose z em are displayed with black dash line in Figure 1, in this paper we survey Civ NALs with z abs z em in spectral datafrom blueward of 10,000 kms 1 until redward of 20,000 kms 1 of Civλ1549 emission lines. 2.2 Spectral analysis Quider et al. (2011) automatically searched for Mg ii doublets on quasar spectra of the SDSS DR4. Using the same parent sample of Quider et al. (2011), Qin et al. (2013) visually inspected Mgii doublets one by one, and detected much more Mgii doublets. Although the visual detection

3 Narrow C iv absorption doublets 3 number of pixels over ±3σ, here σ is given by the Gaussian fitting of a certain candidate. (iv) Significant level of absorption line (here after S.L.). On the normalized spectrum, in terms of the largest depth (S abs ) of the Gaussian fitting relative to unity and the mean value (σ S) of the normalized flux uncertainty around the absorption feature over ±3σ, the S.L. of each absorption line candidate is obtained by Figure 1. Redshisft distributions. The black dash line is the z em of quasars used to search for Civ NALs, and the black solid one is the z abs of all Civ NALs from the current paper as well as Papers I and II. requires a huge amount of time, we follow the same method of Qin et al. (2013), which is also used in Papers I and II, to survey Civ NALs with z abs z em. We simply summarize the main steps used to survey Civ NALs in followings. (i) Considering the underlying continuum and broad emission lines, a combination of cubic spline and Gaussian functions is accepted to fit a pseudo-continuum for each spectrum in iterative model (e.g., Nestor, Turnshek, & Rao 2005; Quider et al. 2011). As an example, we overplot the pseudo-continuum fitting on the quasar spectrum with green line in Figure 2. The normalized spectrum is the quasar spectrum divided by the pseudo-continuum fitting, which is displayed with black line in the lower panel of Figure 2. The red and blue lines displayed in Figure 2 are ±1σ and -2σ flux uncertainty levels divided by the pseudo-continuum fittings, respectively. We disregard absorption features above 2σ and continual absorption troughs (BALs) with line widths larger than 2000 kms 1 at depths greater than ten percent below the pseudo-continuum fitting. (ii) C iv candidates are surveyed on the normalized spectrum. We invokea pair of Gaussian functions to fit each candidate, and visually check one by one. NALs often show a line width of less than a few hundred kms 1. Line widths of mini-bals are between those of NALs and BALs. We conservatively reject the absorption candidates with FWHM 800 kms 1, which are possible mini-bals and beyond the scope of a series of works of Civ NALs. (iii) Equivalent width at rest frame (W r) and corresponding uncertainty (σ w). The W r of each candidate is derived by integrating the Gaussian fitting. Based on the Gaussian fitting, the σ w is evaluated via N P 2 (λ i λ 0)σf 2 i λ 2 i i=1 σ w =, (1) (1+z) N P 2 (λ i λ 0) i=1 where λ i, P(λ i λ 0), σ fi, and N are the wavelength interval spanned by the pixel, line profile centered at λ 0, flux uncertainty normalized by pseudo-continuum fitting, and the S.L. = 1 S abs σ S, (2) (v) C iv NALs are selected from the candidates with S.L. 2.0 for both lines of each doublets. We get rid of absorption lines with W r< 0.2 Å, no matter what values of S.L. are, as we did in Papers I and II. A caveat here is that, the minimum equivalent width limit of absorption line is not a criterion to search for absorption candidates for our program. In other words, our sample completeness is determined by the S.L. rather than the equivalent width (see Section 3). From 37,241 quasar spectra, we have detected a total number of 41,479 Civ NALs with z abs , of which 23,336 doublets with z abs z em are released in Papers I and II, and 18,143 doublets with z abs z em are newly detected in this paper. We distribute z abs of these 41,479 Civ NALs with back solid line in Figure 1, and list the properties of each absorption line in Table 1. Here after, all analyses about Civ NALs are referred to the whole Civ NAL sample that includes doublets released in Papers I and II, and newly detected in this paper. 3 COMPLETENESS AND REDSHIFT PATH We run Monte Carlo simulation to test the efficiency of detection rate of absorbers with different signal-to-noise ratio. Using the spectral data of the BOSS quasar, we simulate 1000 C iv NALs with different S.L.. Then, based on the spectral resolution of the BOSS, we randomly insert simulated C iv doublets onto 400 quasar spectra, where the detected Civ systems have been marked out, in the surveyed spectral region of Civ NALs. Using the method and criteria described in Section 2, we collect the simulated doublets from the 400 quasar spectra. At a given S.L., we define the completeness of detection as the ratio of the number of detected doublets to that of simulated doublets. The result is plotted in Figure 3, which suggests that the detectability tends to be complete in range of S/N λ1548 > 4. In surveyed spectral region, the total redshift path covered by surveys as a function of signal-to-noise ratio per pixel (S/N λ1548 ) under the pseudo-continuum fittings is derived by Z(S/N λ1548 ) = N spec i=1 z max i z min i g i(s/n λ1548,z)dz, (3) where, the sum is over all quasar spectra used to survey absorption lines, the integral boundaries ([zi min,zi max ]) are constrained by surveyed spectral region, the sensitivity function g i(s/n λ1548,z) = 1 if detection threshold S/N lim S/N λ1548, otherwise g i(s/n λ1548,z) = 0. In the left panel of Figure 4, we display the cumulative redshift path as a

4 4 Zhi-Fu Chen et al. Figure 2. Spectrum of quasar SDSS J Upper panel: the pseudo-continuum fitting (green) is overplotted on the quasar spectrum. Lower panel: the quasar spectrum divided by the pseudo-continuum fitting (normalized spectrum). The red and blue lines are the ±1σ and 2σ flux uncertainty levels normalized by the pseudo-continuum fitting, respectively. The detected C iv NALs with S.L. 2 and W r 0.2 Å are marked with purple lines. Table 1. Catalog of C iv narrow absorption line systems SDSS NAME PLATEID MJD FIBERID z em z abs W λ1548 r σ λ1548 W r W λ1551 r Å Å Å Å σ λ1551 W r S.L. λ1548 S.L. λ1551 β Note: β = v/c = ((1 +z em) 2 (1 +z abs ) 2 )/((1 +z em) 2 +(1 +z abs ) 2 ). function of S/N λ1548. The right panel of Figure 4 shows the total redshift path as a function of redshift. The pinnacles at z 2.6 and 2.8, and the obvious features at z > 4.2 are likely due to poor subtractions of skylines, for example, [Oi] and Na D, and OH lines. The spectral regions with strong emission lines often have higher signal-to-noise ratio than those without. Thus, for our identification procedure, the detection rate of C iv doublet per redshift path could be higher for systems with z abs z em, which are located on or near to Civ emission lines, than those with z abs z em. This should be noted in the following statistics and discussion. 4 EQUIVALENT WIDTH DISTRIBUTION AND SATURABILITY The equivalent width of absorption line is an important parameter to reflect absorption strength and elemental abundance of absorber. In our Civ NAL catalog, about 72% of absorbers have Wr λ Å, 29% have Wr λ1548 1Å, and 5% have Wr λ1548 2Å. Using the SDSS quasar spectra, Cooksey et al. (2013) have assembled a catalog of intervening Civ NALs with υ βc > 5000 kms 1. In order to compare our detection to that of Cooksey et al. (2013), we also use the criterion of υ βc = 5000 kms 1 to divide our C iv NAL catalog into two subcatalogs of associated and intervening Civ NALs. Figure 5 shows the distributions of the W r, which is clear that the absorptions of associated systems tend to be stronger than those of intervening

5 Narrow C iv absorption doublets 5 Table 2. The redshift number density with S/N λ z bin N abs Z N/ z Figure 3. Completeness as a function of S/N λ1548. The black line accounts for the systems with S/N λ1548 2, and the red dash line only considers the systems with S/N λ and W r 0.5Å. The detection tends to be complete when S/N λ ones (the black solid and dash lines). Figure 5 also indicates that the W r distribution of intervening systems included in our catalog is similar to that of Cooksey et al. (2013) (the red and black solid lines). These distributions show a break around Wr λ1548 = 0.5Å, which might suggest an uncomplete detection for systems with Wr λ1548 < 0.5Å. The weak systems often have small S.L., and the uncomplete detection suggested by the black solid line of Figure 3 might be mainly introduced by systems with Wr λ1548 < 0.5Å. When the Monte Carlo simulation for estimating the detected completeness (see Section 3) is limited to systems with Wr λ Å, we show the detected completeness as a function of S/N λ1548 with red dash line in Figure 3, which implies that the limit of Wr λ Å does not change the cut, namely S/N λ1548 > 4, of the detected completeness. The saturability of resonant absorption doublet can be obtained by the ratio of absorption strength(dr, Strömgren 1948). For C iv resonant doublet, on the theoretical hand, DR Wr λ1548 /Wr λ1551 can vary from 1 for full saturation to 2 for full unsaturation. The error of the DRcan be evaluated via σ DR DR = ( Wλ1548 r Wr λ1551 ) 2 +( Wλ1548 r (W λ1551 W λ1551 r r ) 2 ) 2 (4) Figure 6displays Wr λ1551 vs Wr λ1548. We findthatabout 73% of Civ NALs fall into the range of 1 DR 2, and only 7.7% of Civ NALs are beyond the range of 1 σ DR DR 2+σ DR. Meanwhile, at Wr λ1548 > 1Å, the DR monotonically declines with increasing Wr λ1548, indicating that the stronger absorptions correspond to the more saturated absorptions that follows the curve of growth. 5 DISCUSSION 5.1 The redshift density evolution of absorbers Systematical survey of NALs on spectra of large quasar sample is in favor of studies of the absorber evolution and the connection between absorbers and galaxies. The redshift density is the count of absorbers with specified limits and [1.4544,1.6544) [1.6544,1.8544) [1.8544,1.9544) [1.9544,2.0544) [2.0544,2.1544) [2.1544,2.2544) [2.2544,2.3544) [2.3544,2.4544) [2.4544,2.6544) [2.6544,2.9544) [2.9544,3.3544) [3.3544,3.8544) [3.8544,4.9224] Note N abs represents the number of detected Civ NALs in corresponding redshift intervals. normalized by redshift path length. In a specified redshift range, the redshift density of absorbers with a specified interval of S/N λ1548 is derived via N abs N z = 1 Z(S/N λ1548 i i ) where Z(S/N λ1548 ) is calculated by Equation 3, and the sum is over all C iv NALs. Here, we only propagate the uncertaintyof Z(S/N λ1548 )intotheuncertaintyof N/ z. Using the spectral flux uncertainty limited by the SDSS, we invoke the Monte Carlo simulation to resample the spectral flux, and then limit the uncertainty of Z(S/N λ1548 ) from Poisson statistics. Thus, the uncertainty of N/ z can be simply propagated by (5) ( N/ z) = N abs ( Z(S/N λ1548 )) Z(S/N λ1548 ) 2 (6) In several redshift ranges and at S/N λ1548 2, Table 2 lists the counts of Civ NALs, Z(S/N λ1548 ), and N/ z, and Figure 7 displays the distribution of N/ z as a function of redshift with red filled circles. We also show the N/ z of Cooksey et al. (2013) with green filled triangles in Figure 7, which are measured from quasar spectra of the SDSS- I/II. It is clear that the N/ z of ours and Cooksey et al. (2013) have parallel evolution profile. The specified values of N/ z depend on methodologies used to determine redshift path. Thus N/ z related to Z(S.L.) cannot be compared directly to N/ z related to Z(W λ r ). The redshift path of this work is limited by S.L., while that of Cooksey et al. (2013) is constrained by W λ r. Therefore, for comparison, the N/ z of Cooksey et al. (2013) plotted in Figure 7 has been scaled arbitrarily. The connection between the metal absorptions and the star formation rate densities has been widely explored in many previous works (e.g., Gauthier & Chen 2011; Ménard et al. 2011; Zhu & Ménard 2013; Huang et al. 2016). The cosmic star formation history peaks at z 2 and decreases by about one magnitude from z 1 to z 0

6 6 Zhi-Fu Chen et al. Figure 4. Left panel: cumulative redshift path versus S/N λ1548. Right panel: total redshift path versus redshift. The pinnacles at z 2.6 and 2.8, and the obvious features at z > 4.2 are likely due to poor subtractions of sky lines of [Oi] and Na D, and OH lines in the SDSS spectra. Figure 5. W r distributions of Civ NALs. The thick black solid line represents Civ NALs with υ βc > 5000 kms 1 (intervening) included in our catalog, the dash line represents those with υ βc 5000 kms 1 (associated) included in our catalog, and the thin red solid line represents C iv NALs of Cooksey et al. (2013). A break near 0.5Åis obvious. (see recent review Madau & Dickinson 2014). In order to establish the connection between N/ z and cosmic star formation histories, we scale the star formation rate densities ( ρ = dm /(dtdadx)) with dx/dz, where X is the comoving distance and A is the comoving area. The scaling ρ is plotted with open diamonds in Figure 7. The global evolution shapes of ρdx/dz and N/ z are very similar. The redshift number density N/ z 1 relates to the volume density of absorbers n vol and their physical crosssection σ phys. Therefore, the evolution of N/ z is possibly dominatedbytheproductofn vol σ phys.civ metalabsorbing 1 More than 60% C IV NALs show DR < 1+σ DR, which means that most of Civ NALs tend to be saturated absorptions. The column density derived from the curve of growth would intrinsically show large uncertainty for saturated NALs, which leads to large uncertainties in the estimation of mass density of Civ absorbing gas. Figure 6. Local data point densities with grayscale. The red lines represent DR = Wr λ1548 /Wr λ1551 = 1 and 2, which correspond the theoretical limits of full saturation and unsaturation, respectively. gas relates to galaxy (e.g., Chen, Lanzetta, & Webb 2001; Adelberger et al. 2005; Tinker & Chen 2008; Martin et al. 2010; Burchett et al. 2013; Bordoloi et al. 2014). Both the number of galaxies and their surrounding gas (halo) likely evolve over cosmic time, so evolutions in both n vol and σ phys can change the value of N/ z: N z n vol(z)σ phys (z). (7) Using the UV luminosity functions from Bouwens et al. (2007), Reddy & Steidel (2009) and Oesch et al. (2010), Cooksey et al. (2013) derived the UV-selected galaxy volume density (n vol ) by integrating the best fitting Schechter function, which increase by a factor of 2 to 3 in redshift range of 5 to 1.5 (see Figure 12 of Cooksey et al. (2013)). In similar redshift range, Figure 7 and Table 2 show that the N/ z increases by a factor of 3 to 4. In this case, the physical cross-section σ phys of absorbing cloud could modestly increase over cosmic time. This increase might be from a comprehensive evolution of absorber s UV ionizing background, metallicity, and geometric size.

7 Narrow C iv absorption doublets 7 Figure 7. The reshift density evolution of absorbers from two surveys: this work (Civ NALs, red filled circles), Cooksey et al. (2013)(C iv NALs, green filled triangles). The star formation rate densities ( ρ ) based on the data from Madau & Dickinson (2014) are displayed with black unfilled diamonds. The mappings between the left-axes of the triangles and the circles, and between the left and right y-axes are arbitrary. 5.2 The quasar-frame velocity distribution NALs are mainly arisen from: (1) material of intervening galaxies lying on quasar sightlines; and (2) material associated with quasar themselves, which can be related to quasars host galaxies, winds/outflows and surrounding environment. The former produces intervening NALs with z abs z em, and the latter brings about associated NALs with z abs z em. Associated NALs are alluring since they provide a powerful tool to explore AGN feedback. BALs are unambiguously associated with quasars. However, it is difficult to identify whether one NAL is related to quasar or not. The line variability and partial coverage analysis are frequently applied to diagnose the associated NALs. Only small sample of variable NALs are detected up to nowadays (e.g., Narayanan et al. 2004; Wise et al. 2004; Hacker et al. 2013; Chen et al. 2015). Therefore, the statistics is often utilized to decompose associated population from a large NAL catalog (e.g., Misawa et al. 2007; Nestor, Hamann, & Rodriguez Hidalgo 2008; Wild et al. 2008). If NALs are associated to quasars, then the differences between z abs and z em should be dominated by motions of absorbers relative to quasars, and the relative velocities can be derived by β υ c = (1+zem)2 (1+z abs ) 2 (1+z em) 2 +(1+z abs ) 2. (8) We note that z abs determined from NALs have high accuracy, while z em measured from broad emission lines can show an uncertainty from 100 kms 1 to 3000 kms 1 (e.g., Hewett & Wild 2010; Shen et al. 2011; Shen 2016). Using the near-infrared quasar spectra, Shen (2016) revealed that, after accounting for measurement errors, redshifts estimated from broad Mg iiλ2798 emission lines have an intrinsical uncertainty of about 200 kms 1 with respect to the redshifts determined by [O iii]λ5007 emission lines. This small uncertaintyimplies that Mgiiemission line is abetter indicator to estimate quasar systemic redshift relative tociii],civ, and Ly α emission lines (e.g., Hewett & Wild 2010; Shen et al. 2011). Chen, Gu, & Chen (2015) has systematically surveyed Mgii NALs on 43,260 BOSS quasar spectra of the DR9Q from redward of Ly α to red wing of Mg iiλ2798 emission lines, which results in a total number of 18,598 Mgii doublets with z abs The method used in Chen, Gu, & Chen (2015) to search for NALs is exactly the same as this work. In order to compare the β distribution of Civ NALs to that of Mgii NALs, we make use of the Mgii NALs included in Chen, Gu, & Chen (2015). We limit Civ NALs to quasars with z em 2.2 so that Mgiiλ2798 emission lines can be available by the BOSS spectra. Surveyed spectral region of Civ NALs is between Ly α and Civλ1549 emission lines, which constrains β < 0.16 for C iv NALs. And that of Mgii NALs is between Ly α and Mgiiλ2798 emission lines, which permits measurement of much higher β values for Mgii NALs. In order to make a consistent comparison, we limit both Civ and Mgii NALs with β < Taking into account of z em 2.2 and β < 0.15, there are 5922 Civ NALs and 3012 Mgii NALs. We display the β distributions of both Civ and Mgii NALs in Figure 8. In the calculation of β, we have made an assumption that all NALs are produced by quasar self-absorptions. In fact, the catalogs contain absorptions from (1) the first population: cosmologically intervening galaxies; (2) the second population: quasar host galaxies or galaxy clusters/groups that quasars reside in, which is often called as environmental absorptions; and (3) the third population: quasar outflows. Although ultra-fast outflows have been detected by many researchers (e.g., Hamann et al. 2011; Tombesi et al. 2011; Chen & Qin 2013; Gupta et al. 2013), the vast majority of absorbers with a large β belong to the first population. When absorbers fall into the second population, some absorbers move toward the observer and some away, and the β would show a normal distribution centred at β = 0. In addition, galaxies around quasars are clustered/grouped, so one expects an excess of associated absorbers with respect to cosmologically intervening absorbers. Absorbers of the third population move towards the observer, thus they have β > 0 and are expected to destroy the normal distribution of β of the second population absorptions. Figure 8 clearly shows that there are an abnormal excess in range of β < 0.02 and a flat tail in range of β > This suggests that the data in range of β > 0.02 would be dominated by absorptions of the first population, and those in range of β < 0.02 be mainly contributed from absorptions of the second and third populations. In theory, quasar outflows do not affect the β distribution of associated absorbers with β < 0. In Figure 8, absorbers with β > 0.06 are probably contributed from cosmologically intervening absorptions. We explore a combination of a Gaussian function accounting for environmental absorptions and a linear function accounting for intervening absorptions to fit the data with β < 0 and > This results in Gaussian components with σ = for C iv NALs and σ = for Mgii NALs. The fitting results are plotted with color lines in Figure 8. It is clear that a combination of a Gaussian and a linear functions can not well characterize the overall distributions, and the excesses located at β 0.01 are very remarkable with respect to

8 8 Zhi-Fu Chen et al. Figure 8. Relative velocity distributions of absorbers with β < 0.15 and z em 2.2. The y-axes are the counts of absorbers normalized by the total numbers of absorbers contained in corresponding subsamples. The red curves are the reduced χ 2 fittings of a combination of a Gaussian function centred at β = 0 and a linear function using the data with β < 0 and > The blue and orange dash lines are the Gaussian and linear components, respectively. Left panel: the KS test suggests a probability of P < 10 6 that the two populations draw from the same parent sample. Middle panel: Civ absorbers with σ = for the Gaussian component. Right panel: Mgii absorbers with σ = for the Gaussian component. fitting curves, which likely are the contribution of quasar outflows. Since Civ (64.49 ev) has much higher ionization potential than Mgii ( ev), we expect the associated Civ and Mg ii NALs to show different statistical properties. Recalling Figure 8, in range of β < 0.02, the β distribution of Mgiiabsorbers ismorecompactthanthesame valuesofciv absorbers. KS test derives a probability of P < 10 6 that theβ ofciv andmgiiabsorbers drawfrom thesame parent sample. The multi-component fittings overplotted with red lines imply that the dispersion of C iv environmental absorptions (σ = ) is much larger than the same value of Mg ii environmental absorptions (σ = ). In addition, with respect to multi-component fittings, the excessive tail of Mgii absorbers reaches a value of β 0.02, but that of Civ absorbers can extend out to a higher value of β > This suggests that the higher ionization outflows have the larger maximum velocity. Using Civ and Mgii NALs detected from quasar spectra of the SDSS DR3, the model of Wild et al. (2008) predicts that the scales that the quasar photoevaporates environmental absorbers are related to ionization energies of ions. The radius of 300 kpc for survived Civ absorbers is much smaller than that of 800 kpc for survived Mg ii absorbers. This discrepancy is also clearly reflected by our multi-component fittings (see Figure 8): σ = for Civ absorbers, and σ = for Mgii ones. Using the emission redshifts of the SDSS pipeline, Nestor, Hamann, & Rodriguez Hidalgo (2008) obtained a dispersion with avalueof913 kms 1 for Civ environmental absorbers after deconvolving the redshift uncertainty. When based on redshifts determined by Mg ii emission lines, they found a value of about 450 kms 1. For our multi-component fitting, we obtain a dispersion with a value of 728 kms 1, which is slightly larger than that of 450 kms 1. Our limit of z em 2.2can reduce,onacertainextent,theuncertainties of emission redshifts. However, we can not guarantee that the significant Mg ii emission line was imprinted on each quasar spectrum. Thus, the uncertainties of emission redshifts could be larger than those determined by Mg ii emission lines. 5.3 Dependence on the radio detection The observational area of the Faint Images of the Radio Sky at Twenty centimeter survey (FIRST, Becker, White, & Helfand 1995) coincides with that of the SDSS. Over most of the area, the FIRST survey is limited to a detection threshold of 1 mjy. In order to investigate the dependence of associated absorber properties on radio detection, we match our parent quasar sample with the FIRST within a radius of 30. Figure 8 shows that some outflow absorbers can reach a high velocity. Here, we loosely constrain associated absorbers with β < Combining the limits of β < 0.03 and z em 2.2, there are 2741 and 198 Civ NALs imprinted on the spectra of FIRST undetected and detected quasars, respectively. In terms of the luminosity at 20 cm, the quasar with a luminosity > W/Hz is defined as radio-loud one, inversely radio-quiet one. According to this definition, the FIRST detected quasars of the 198 C iv absorbers belong to radio-loud population. Because the small sample of the FIRST detected quasars, we do not discuss the dependence of associated absorber properties on radio morphologies. We display β distributions in Figure 9. No significant differenceis shown betweenthe β ofthe FIRSTdetected and undetected NALs. This similar result is also exhibited in the W r distributions (see Figure 10). The β distributions do not coincide with those of Wild et al. (2008), who found that the β measured from radio-loud quasar spectra are more compact at β 0 than the same values measured from radio-quiet ones. There are only 69 FIRST detected quasars satisfied our criteria. This number is too small to allow us to further investigate the different β distributions between our results and Wild et al. (2008). 5.4 Dependence on the UV continuum luminosity The ionizing radiation from quasars produces feedback to surrounding gas, which can affect gaseous ionization state, density and kinematics. The dependence of this feedback on the quasar luminosity can be studied by absorbers surrounding quasars (e.g., Pan & Chen 2013; Johnson, Chen, & Mulchaey 2015). We exclude quasars that

9 Narrow C iv absorption doublets 9 Figure 9. Relative velocity distributions of absorbers with β < 0.03 and z em 2.2. The y-axis are the counts of absorbers normalized by the total numbers of absorbers contained in corresponding subsamples. The solid (dash) lines represent the quasars without (with) the FIRST detection. KS test yields a P = 0.03 Figure 11. Relative velocity distributions of absorbers with β < 0.03 and z em 2.2. The y-axis are the counts of absorbers normalized by the total numbers of absorbers contained in corresponding subsamples. The solid (dash) lines represent quasars with continuum luminosities less (greater) than corresponding median values. KS test yields a P < Figure 10. W r distributions of absorbers with β < 0.03 and z em 2.2. The y-axis are the counts of absorbers normalized by the total numbers of absorbers contained in corresponding subsamples. The solid (dash) lines represent the quasars without (with) the FIRST detection. KS test yields a P = 0.7 Figure 12. W r distributions of absorbers with β < 0.03 and z em 2.2. The y-axis are the counts of absorbers normalized by the total numbers of absorbers contained in corresponding subsamples. KS test yields a P < the flux data at 1350Å are not available by the BOSS spectra due to quasar redshifts and that are located at z em> 2.2, and limit absorbers with β < This results in a sample of 2866 Civ NALs. We divide the resulting sample at median value of the continuum luminosity at rest frame 1350Å (L 1350 = ergs 1 ) into low- and high-luminosity quasar samples, where continuum luminosities are measured from pseudo-continuum fitting. The β and W r distributions of both low- and highluminosity quasar samples are shown in Figures 11 and 12, respectively. KS test yields a P < for the β and a P < for the W r, which suggests that the ionizing radiation of quasars acts on both β and W r distributions. We also note that the low-luminosity quasars have stronger absorption strengths. The SDSS quasars are obtained by flux limited surveys, and consequently the high luminosity quasars can be ob- tained at high redshifts. Therefore, the differences of β and W r distributions between low- and high-luminosity quasars might be the evolution result of quasar surrounding gas. We find that both low- and high-luminosity quasars have similar median z em. Thus, the differences reflected by Figures 11 and 12 are not originated in the evolution of quasar surrounding gas. It is interesting that the high-luminosity quasars tend to show larger velocity (β) Civ NALs with smaller absorption strengths, which might be related to the radiative acceleration by quasar luminosity and ionization condition of CGMs around quasars. For high-luminosity quasars, the enormous power of outflows and the energetic radiation from accretion disc could produce stronger feedback to surrounding gas of quasars. In this case, surrounding gas of highluminosity quasars is possibly accelerated to higher velocity and photo-ionized to higher level.

10 10 Zhi-Fu Chen et al. 6 SUMMARY Using the BOSS quasar spectra (DR9Q), combining the two early works (Papers I and II), we have detected a total numberof41,479 Civ NALswith z abs inwavelength range between Ly α and C ivλ1549 emission lines. The detectability tend to be complete when S/N λ1548 > 4. Ours main results are as follows. (i) 92.3% of Civ NALs fall into the range of 1 σ DR DR 2+σ DR, and tend to be the more saturated absorptions with the stronger absorptions that follows the curve of growth. The associated systems show significant absorptions when compared to the intervening ones. (ii) The redshift density evolution shape of C iv absorbers is similar to the history of the cosmic star formation, which indicates that C iv absorbers may be connected to the star formation activities within galaxies. Thus, absorption lines imprinted on the spectra of background objects are possibly a powerful tool to investigate faint and even invisible galaxies. (iii) Comparing β distributions of C iv and Mg ii absorbers, we find clear distinctions. The β of Civ absorbers are morecompact toβ 0, andhavelonger extendedtail. In addition, the dispersion of C iv environmental absorptions is much larger than that of Mg ii environmental absorptions. The Civ has an ionization energy of ev, which is obviously different from that of ev of the Mgii. The clearly different ionization energies of Civ and Mgii may play an important role to the distinct β properties between the two kinds of absorbers. (iv) For associated C iv absorbers (β < 0.03), when compared with the low-luminosity quasars, high-luminosity quasars tend to have larger β and weaker absorptions. The line number density of absorbers depends not only on their volume density, but on physical cross-section and ionizing background of absorbing clouds. One of our future works will investigate the properties of absorbers that intersect both sightlines of a quasar pair (two close quasars on sky). This might provide a useful limit to the geometric cross-section of absorbing cloud. Deep imaging/spectroscopy observations to this kind of absorbers might also provide useful information about absorption line formation by comparing absorption strength and profile of absorption lines and characteristics of galaxy halo, such as gas distribution, star formation rate, orientations of disc and outflow, and so on. The properties of CGMs and outflows of quasars are crucial to comprehend the quasar feedback and evolution of host galaxy. Time variability and partial coverage analysis is two efficient ways to identify absorption lines of quasar outflow and environment. Therefore, high resolution observations of the quasars with variable absorption line could provide well opportunity to investigate the properties of quasar outflow and CGM. The authors are very grateful to the anonymous referee for his/her critical comments and instructive suggestions, which significantly strengthened the analyses in this work. This work was supported by the National Natural Science Foundation of China (NO , grants , and ), National Basic Research Program (973 program No. 2013CB834905), the Guangxi Natural Science Foundation (2015GXNSFBA139004), and the Guangxi University of Science and Technology research projects (No. KY2015YB289). Y.M.C. acknowledges support from the Opening Project of Key Laboratory of Computational Astrophysics, National Astronomical Observatories, Chinese Academy of Sciences. Funding for SDSS-III has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the U.S. Department of Energy Office of Science. The SDSS-III web site is SDSS-III is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS- III Collaboration including the University of Arizona, the Brazilian Participation Group, Brookhaven National Laboratory, Carnegie Mellon University, University of Florida, the French Participation Group, the German Participation Group, Harvard University, the Instituto de Astrofisica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, Max Planck Institute for Extraterrestrial Physics, New Mexico State University, New York University, Ohio State University, Pennsylvania State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University. REFERENCES Abazajian K. N., et al., 2009, ApJS, 182, Adelberger K. L., Shapley A. E., Steidel C. C., Pettini M., Erb D. K., Reddy N. A., 2005, ApJ, 629, 636 Balashev S. A., Klimenko V. V., Ivanchik A. V., Varshalovich D. A., Petitjean P., Noterdaeme P., 2014, MNRAS, 440, 225 Arrigoni Battaia F., Hennawi J. F., Cantalupo S., Prochaska J. X., 2016, arxiv, arxiv: Becker R. H., White R. L., Helfand D. J., 1995, ApJ, 450, 559 Bordoloi R., et al., 2014, ApJ, 796, 136 Borthakur S., et al., 2015, ApJ, 813, 46 Bouché N., Murphy M. T., Kacprzak G. G., Péroux C., Contini T., Martin C. L., Dessauges-Zavadsky M., 2013, Sci, 341, 50 Bouwens R. J., Illingworth G. D., Franx M., Ford H., 2007, ApJ, 670, 928 Burchett J. N., Tripp T. M., Werk J. K., Howk J. C., Prochaska J. X., Ford A. B., Davé R., 2013, ApJ, 779, L17 Cantalupo S., Arrigoni-Battaia F., Prochaska J. X., Hennawi J. F., Madau P., 2014, Natur, 506, 63 Cen R., Chisari N. E., 2011, ApJ, 731, 11 Chen H.-W., Lanzetta K. M., Webb J. K., 2001, ApJ, 556, 158 Chen Z.-F., Gu Q.-S., Chen Y.-M., 2015, ApJS, 221, 32 Chen Z.-F., Gu Q.-S., Chen Y.-M., Cao Y., 2015, MNRAS, 450, 3904 Chen Z.-F., Li M.-S., Huang W.-R., Pan C.-J., Li Y.-B., 2013, MNRAS, 434, 3275 Chen Z.-F., Qin Y.-P., 2013, ApJ, 777, 56 Chen Z.-F., Qin Y.-P., Pan C.-J., Huang W.-R., Qin M., Wu H.- N., 2014a, ApJS, 210, 7 (PaperI) Chen Z.-F.,QinY.-P.,QinM.,Pan C.-J.,PanD.-S.,2014b, ApJS, 215, 12 (PaperII) Chartas G., Eracleous M., Misawa T., Giustini G., Charlton J., 2012, ASPC, 460, 37

11 Narrow C iv absorption doublets 11 Chartas G., et al., 2009, NewAR, 53, 128 Dawson K. S., et al., 2013, AJ, 145, 10 Cooksey K. L., Kao M. M., Simcoe R. A., O Meara J. M., Prochaska J. X., 2013, ApJ, 763, 37 Cooksey K. L., Prochaska J. X., Thom C., Chen H.-W., 2011, ApJ, 729, 87 Davé R., Finlator K., Oppenheimer B. D., 2012, MNRAS, 421, 98 Dekel A., et al., 2009, Natur, 457, 451 Eisenstein D. J., et al., 2001, AJ, 122, 2267 Eisenstein D. J., et al., 2011, AJ, 142, 72 Farina E. P., Falomo R., Scarpa R., Decarli R., Treves A., Kotilainen J. K., 2014, MNRAS, 441, 886 Fechner C., Richter P., 2009, A&A, 496, 31 Ford A. B., Davé R., Oppenheimer B. D., Katz N., Kollmeier J. A., Thompson R., Weinberg D. H., 2014, MNRAS, 444, 1260 Fox A. J., Petitjean P., Ledoux C., Srianand R., 2007, A&A, 465, 171 Fox A. J., Prochaska J. X., Ledoux C., Petitjean P., Wolfe A. M., Srianand R., 2009, A&A, 503, 731 Gauthier J.-R., Chen H.-W., 2011, MNRAS, 418, 2730 Gupta A., Mathur S., Krongold Y., Nicastro F., 2013, ApJ, 772, 66 Gunn J. E., et al., 2006, AJ, 131, 2332 Hacker T. L., Brunner R. J., Lundgren B. F., York D. G., 2013, MNRAS, 434, 163 Hamann F., Kanekar N., Prochaska J. X., Murphy M. T., Ellison S., Malec A. L., Milutinovic N., Ubachs W., 2011, MNRAS, 410, 1957 Hennawi J. F., Prochaska J. X., 2013, ApJ, 766, 58 Hennawi J. F., Prochaska J. X., Cantalupo S., Arrigoni-Battaia F., 2015, Sci, 348, 779 Hewett P. C., Wild V., 2010, MNRAS, 405, 2302 Huang Y.-H., Chen H.-W., Johnson S. D., Weiner B. J., 2016, MNRAS, 455, 1713 Husband K., Bremer M. N., Stanway E. R., Lehnert M. D., 2015, MNRAS, 452, 2388 Johnson S. D., Chen H.-W., Mulchaey J. S., 2015, MNRAS, 452, 2553 Kacprzak G. G., Churchill C. W., Barton E. J., Cooke J., 2011, ApJ, 733, 105 Kacprzak G. G., Churchill C. W., Steidel C. C., Murphy M. T., Evans J. L., 2007, ApJ, 662, 909 Kacprzak G. G., Muzahid S., Churchill C. W., Nielsen N. M., Charlton J. C., 2015, ApJ, 815, 22 Keeney B. A., Stocke J. T., Rosenberg J. L., Danforth C. W., Ryan-Weber E. V., Shull J. M., Savage B. D., Green J. C., 2013, ApJ, 765, 27 Kollmeier J. A., Miralda-Escudé J., Cen R., Ostriker J. P., 2006, ApJ, 638, 52 Lü L.-Z., Qin Y.-P., Gupta A. C., 2007, ApJ, 669, 74 Ménard B., Wild V., Nestor D., Quider A., Zibetti S., Rao S., Turnshek D., 2011, MNRAS, 417, 801 Madau P., Dickinson M., 2014, ARA&A, 52, 415 Martin C. L., Scannapieco E., Ellison S. L., Hennawi J. F., Djorgovski S. G., Fournier A. P., 2010, ApJ, 721, 174 Martin D. C., Chang D., Matuszewski M., Morrissey P., Rahman S., Moore A., Steidel C. C., 2014, ApJ, 786, 106 Martin D. C., Matuszewski M., Morrissey P., Neill J. D., Moore A., Cantalupo S., Prochaska J. X., Chang D., 2015, Natur, 524, 192 Matejek M. S., Simcoe R. A., 2012, ApJ, 761, 112 Misawa T., Charlton J. C., Eracleous M., Ganguly R., Tytler D., Kirkman D., Suzuki N., Lubin D., 2007, ApJS, 171, 1 Murga M., Zhu G., Ménard B., Lan T.-W., 2015, MNRAS, 452, 511 Muzahid S., Kacprzak G. G., Churchill C. W., Charlton J. C., Nielsen N. M., Mathes N. L., Trujillo-Gomez S., 2015, ApJ, 811, 132 Narayanan D., Hamann F., Barlow T., Burbidge E. M., Cohen R. D., Junkkarinen V., Lyons R., 2004, ApJ, 601, 715 Nestor D., Hamann F., Rodriguez Hidalgo P., 2008, MNRAS, 386, 2055 Nestor D. B., Turnshek D. A., Rao S. M., 2005, ApJ, 628, 637 Nielsen N. M., Churchill C. W., Kacprzak G. G., Murphy M. T., Evans J. L., 2015, ApJ, 812, 83 Oesch P. A., et al., 2010, ApJ, 725, L150 Pâris I., et al., 2012, A&A, 548, A66 Pan C.-J., Chen Z.-F., 2013, JApA, 34, 357 Pérez-Ràfols I., Miralda-Escudé J., Lundgren B., Ge J., Petitjean P., Schneider D. P., York D. G., Weaver B. A., 2015, MNRAS, 447, 2784 Qin Y.-P., et al., 2013, PASJ, 65, 8 Quider A. M., Nestor D. B., Turnshek D. A., Rao S. M., Monier E. M., Weyant A. N., Busche J. R., 2011, AJ, 141, 137 Reddy N. A., Steidel C. C., 2009, ApJ, 692, 778 Richter P., 2012, ApJ, 750, 165 Rubin K. H. R., Prochaska J. X., Koo D. C., Phillips A. C., 2012, ApJ, 747, L26 Savage B. D., Kim T.-S., Wakker B. P., Keeney B., Shull J. M., Stocke J. T., Green J. C., 2014, ApJS, 212, 8 Seyffert E. N., Cooksey K. L., Simcoe R. A., O Meara J. M., Kao M. M., Prochaska J. X., 2013, ApJ, 779, 161 Shen Y., 2016, ApJ, 817, 55 Shen Y., et al., 2011, ApJS, 194, 45 Shen Y., Ménard B., 2012, ApJ, 748, 131 Smee S. A., et al., 2013, AJ, 146, 32 Songaila A., 2005, AJ, 130, 1996 Springel V., Hernquist L., 2003, MNRAS, 339, 289 Strömgren B., 1948, ApJ, 108, 242 Tinker J. L., Chen H.-W., 2008, ApJ, 679, Tombesi F., Cappi M., Reeves J. N., Palumbo G. G. C., Braito V., Dadina M., 2011, ApJ, 742, 44 Tripp T. M., et al., 2011, Sci, 334, 952 Tripp T. M., Sembach K. R., Bowen D. V., Savage B. D., Jenkins E. B., Lehner N., Richter P., 2008, ApJS, 177, Vanden Berk D., et al., 2008, ApJ, 679, Wild V., et al., 2008, MNRAS, 388, 227 Wise J. H., Eracleous M., Charlton J. C., Ganguly R., 2004, ApJ, 613, 129 York D. G., et al., 2000, AJ, 120, 1579 Zhu G., Ménard B., 2013, ApJ, 773, 16 Zhu G., Ménard B., 2013, ApJ, 770, 130