{'fraction_non_alphanumeric': 0.031186825667234526, 'fraction_numerical': 0.027541169789892107, 'mean_word_length': 4.40588163064833, 'pattern_counts': {'":': 0, '<': 4, '<?xml version=': 0, '>': 5, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
The Quasar Outflow Contribution to AGN Feedback: VLT Measurements of SDSS J0318-0600

We present high spectral resolution VLT observations of the BAL quasar SDSS J0318-0600. This high quality data set allows us to extract accurate ionic column densities and determine an electron number density of n_e=10^3.3 +/- 0.2 cm^-3 for the main outflow absorption component. The heavily reddened spectrum of SDSS J0318-0600 requires purely silicate dust with a reddening curve characteristic of predominately large grains, from which we estimate the bolometric luminosity. We carry out photoionization modeling to determine the total column density, ionization parameter and distance of the gas and find that the photionization models suggest abundances greater than solar. Due to the uncertainty in the location of the dust extinction, we arrive at two viable distances for the main ouflow component from the central source, 6 and 18 kpc, where we consider the 6 kpc location as somewhat more physically plausable. Assuming the canonical global covering of 20% for the outflow and a distance of 6 kpc, our analysis yields a mass flux of 120 M_sun yr^-1 and a kinetic luminosity that is ~0.1% of the bolometric luminosity of the object. Should the dust be part of the outflow, then these values are ~4x larger. The large mass flux and kinetic luminosity make this outflow a significant contributor to AGN feedback processes.


Introduction
Quasar outflows, seen as ultraviolet blueshifted absorption features in their spectra, may contribute to a variety of AGN feedback mechanisms (Elvis 2006). Such feedback processes are invoked to explain the: evolution of the super massive black hole (SMBH) (Silk & Rees 1998;Blandford & Begelman 1999); evolution of the host galaxy (Di , shape of the observed quasar luminosity function (Hopkins et al. 2005a;Springel et al. 2005;Menci et al. 2006) and the enrichment of the intra-cluster and intergalactic media (Moll et al. 2007).
The contribution of quasar outflows to AGN feedback depends mainly on their mass flux (Ṁ ) and kinetic luminosity (Ė k ), both of which depend linearly on the the total hydrogen column density (N H ) of the outflow and its characteristic distance (R) from the central source (see Equations 9 and 11 in § 6). Until recently, only loosely constrained estimates have been available for N H and R in the same outflow. Morris et al. (1986) and Tripp et al. (1996) were only able to find lower limits for the distance to outflows in GC 1556+335 and HS 1946+7658, respectively. Wampler et al. (1995), analyzing the spectrum of QSO 0059−2735, constrained the distance of the outflow to a range spanning over an order of magnitude. de Kool et al. (2001) determined the distance to the outflow observed in FIRST J1044+3656 to within a factor of 100, and had a similar uncertainty in the inferred N H (see their table 3), which led to an uncertainty inṀ andĖ k of approximately one order of magnitude. In the cases of FIRST J1214+2803 & FIRST J0840+3633, de Kool et al. (2002a,b) were able to constraiṅ M andĖ k to a similar accuracy. For the quasar 3C 191, Hamann et al. (2001) estimateḋ M andĖ k to within a factor of a few. However, beside their internal uncertainties, all the works cited above suffered from unquantified systematic errors due to the use of inadequate absorption models that led to unreliable measurements of column densities in the observed troughs, which are crucial for determining almost every physical aspect of the outflows: the ionization equilibrium and abundances, number density, distance, mass flux, and kinetic luminosity.
In order to determine accurateṀ andĖ k in quasar outflows we launched a multistage research program: We first developed analysis techniques that allow us to determine reliable trough column densities (Arav et al. 1999a(Arav et al. ,b, 2002(Arav et al. , 2005Scott et al. 2004b;Gabel et al. 2005). This was followed by identifying outflow targets in which we can find diagnostics for both N H and R. We searched through ∼50,000 QSOs with z>0.25 and r<19.3 in the SDSS Data Release 6 (Adelman- McCarthy et al. 2008), ∼100 objects were found to exhibit intrinsic outflow absorption systems containing troughs from excited or metastable levels of one or more of the following ions: He i, Si ii, Fe ii, Fe iii or Ni ii. The key spectral features to obtain distances and total column densities for outflows are a) absorption troughs of excited and metastable levels from which we can determine the electron density, and b) the availability of troughs from many ions, including lines from multiple ions of the same element. We then began observing these objects with a combination of high spectral resolution, large spectral coverage and high signal-to-noise ratios that in aggregate allows for the use of the analysis techniques mentioned above. Arav et al. (2008) and Korista et al. (2008) used this analysis method on VLT observations of QSO 2359-1241, deriving N H and hydrogen number density determinations with accuracies better than 25% and 20%, respectively.
The object we analyze in this paper, SDSS J0318−0600 was found in our search through the SDSS data base (although first discovered by Hall et al. 2002). We chose this object for a high resolution follow up observation, due to its high apparent brightness and its suitability for VLT observations. SDSS J0318−0600 is a luminous quasar (log L Bol ≈47.7 ergs s −1 , see §5; z=1.9668; Hall et al. 2002;Schneider et al. 2003) that shows a myriad of outflow components ranging from a high ionization broad absorption line system to several narrow absorption line (NAL) features spanning a velocity range from −2800 to −7500 km s −1 , with respect to the quasar's rest frame. In a survey by Hall et al. (2002), it was classified as an iron low-ionization broad absorption line system (FeLoBAL) and listed as an unusual quasar due to the BAL located at ∼27,000 km s −1 and the system of NALs, which appears as a BAL in higher ionization lines such as C iv. They noted that the Fe ii UV1, UV2 and UV3 multiplets were detected in absorption and that the UV38 multiplet appeared to be present as well. It was also noted that SDSS J0318−0600 has a heavily reddened spectrum. Prochter et al. (2006) listed this object as a strong Mg ii absorber with a rest frame absorption width >1.4 A, in a search of Mg ii absorbing objects. Herbert-Fort et al. (2006) showed that the object's absorption was strong, particularly in Si ii λ 1808, and labeled the absorber as a damped Lyman alpha system (DLA). However, in our VLT data, we detect absorption from excited ionic energy levels that conclusively classify the system as an outflow (Hamann et al. 1997).
The outline for this paper is as follows, in § 2, we present the VLT eschelle observations and their reduction. For § 3, we discuss our methods of determining and measuring individual kinematic components. We then concentrate on a physical analysis of the main component of the outflow: in section 4, we present our photoionization analysis of the spectrum and determine the hydrogen column density (N H ), electron density (n e ), distance (R) and temperature (T) for the main absorption component. We derive the mass flux and kinetic luminosity in section 5, and in § 6 we discuss our results. During October of 2006, we obtained 8 echelle spectra, each with an exposure time of 2850 sec, of the object SDSS J0318−0600 (z=1.967, Schneider et al. 2007, RA 03:18:56.7, Dec −06:00:37.4) with the VLT. The observations were carried out by European Southern Observatory (ESO) staff and we present the observational parameters in Table 1. The spectra were taken using the UV-Visual Echelle Spectrograph (UVES) with a slit width of 1.0 arcsec, and a central wavelength position of 5800 and 3460Å for the red (Grating/Filter=CD1/HER5) and blue (Grating/Filter=CD3/SHP700) detectors respectively. A slit width of 1.0 arcsec gives a resolving power of ∼40000. The data were read out in a 2×2 pixel-binned fashion, in high gain mode.
The data were extracted and calibrated using the standard UVES data reduction pipeline, based on ECHELLE routines, which is part of the reduction package MIDAS (Ballester et al. 2000). This calibration process can introduce quasi-periodic oscillations in amplitude in the final spectral product (see Arav et al. 2008, on their analysis of QSO 2359-1241), which can create up to 2-3% systematic errors in continuum placement. The observations were corrected to heliocentric and subject to a cosmic ray removal procedure (σ > 3.5). We re-binned the data to 7.0 km s −1 , approximately one resolution element, and normalized the spectrum using the methods described in Arav et al. (2008).

Detection & Measurements
The spectrum of SDSS J0318−0600 contains a plethora of absorption features available for potential diagnostics. In Figure 1, we show a sample of the VLT spectrum with line identifications, and present the full spectrum with identifications in our online figure. These identifications come from matching the Mg ii doublet lines and plotting location identifiers for other ionic lines. We also provide an online table for all identified lines, their wavelengths, oscillator strengths and references. We show in Figure 1 a portion of the full online figure that contains readily visible lines from both the ground as well as excited levels of Fe ii (lower energy levels up to 1873 cm −1 ).

Identification of Kinematic Components
We identify each of the individual kinematic absorption components in the spectrum using the Al ii λ1671 absorption line profile, which is a singlet transition with no nearby strong absorption features. In Figure 2, we show that there are approximately 11 separate kinematic components (labeled a − k), with component e clearly a blend of multiple insep- The redshift of the quasar is in black located at the top right. The log(gf ) values of the identified transitions are shown above the spectrum, where g is the statistical weight and f is the oscillator strength. The size of the vertical blue line provides the value, which can be read from the y-axis scale at the top right. Each identifier lists the ion, wavelength, and lower energy level (in 1000s of cm −1 ). The spectrum error is plotted in solid black below the spectrum (≈0.05 in normalized flux). arable components. Using Al ii as a template for the absorbers we search other ionic species for each of the 11 components.
Next, we identify the broad, strong i component in as many lines as possible and select unblended troughs from which we are able to measure ionic column densities. Figure 1 shows the UV 38 Fe ii multiplet, which exhibits several absorption features from component i. We note that the absorption feature at ∼4727Å is the k component of the Fe ii λ1608 resonance line. While this identification scheme is sufficient to locate the majority of the lines that show an i component, we also use the Al ii as a template (see §3.3), to search for the other 10 components in the Fe ii UV 1, 2 & 3 multiplet absorption complexes.
-7 - Fig. 2.-The main outflow troughs in SDSS J0318−0600, ordered by wavelength. We identify the 11 kinematic components found in Al ii from high to low velocities, which we label a through k. In the plots of Fe ii λ2383, λ2396, and λ2600 there are several contaminate lines present due to the telluric A and B bands as well as blending from other Fe ii troughs (identified in online figure). Each doublet also shows blending between components (e.g. C iv, Mg ii).  a Reflects the statistical errors from line fitting. This does not reflect the reliability of the measurement methods, which have increasing systematic error in order of increasing label number (see section 3.2).
b Method used to determine column density, described by the respective label number; 1: Power Law fitting, 2: Velocity Dependent Covering Factor fitting, 3: Apparent optical depth The Fe ii multiplets UV 1, 2 & 3 are the best candidates to identify Fe ii lines from the other 10 components to measure, because they contain the Fe ii transitions with the strongest oscillator strengths. However, this task proves to be non-trivial in SDSS J0318-0600. Identifying measurable lines from a saturated and heavily blended region (see Figure  2) is difficult, especially when considering the large number of strong resonance and excited state lines. Our strategy is to make what line identifications for all 11 components that we can using our Al ii template and search for other lines where the blending is not so prevalent. We provide the full list of detected ions, transition wavelengths, their oscillator strengths and lower level energies in the online table.
We find that for the lowest and highest velocity components (k and a), enough measurements or limit determinations are available that we can effectively determine the electron density (n H ), the ionization parameter (U H ), and the total hydrogen column density (N H , see §4). Unfortunately, blending prevents us from identifying enough measurable troughs for components b − h for proper analysis. We opt for the limiting case by taking the sum of these latter systems (see §3.2).

Column Density Determinations
With a list of measurable absorption troughs, we extract ionic column densities. We begin the process by measuring lines from the strongest component, (i), where we employ three primary methods for determining the column densities. We describe three additional methods for column density determination that are less reliable in section 3.3. We provide the measured values of ionic column densities as well as the method applied in Table 2.
Our first method of measuring column densities is to model a set of observed Fe ii resonance lines (e.g., Figure 3) using an inhomogenous absorption model, where we approximate the gas distribution across our line of sight to the background source with a power law (see de Kool et al. 2002c;Arav et al. 2008, especially Fig. 4, for a full description). The optical depth at a given velocity within the outflow is described by: where x is the spatial dimension in the plane of the sky (simplified to one dimension with no loss of generality, see Arav et al. 2005), a is the power law distribution index, and τ max (v) is the highest value of τ at a given velocity. We simultaneously solve for both τ max and a for each velocity resolution element of the resonance Fe ii lines and apply this distribution (a) to the excited Fe ii lines. We then apply the same distribution dictated by the Fe ii lines to all ions because this is the simplest assumption that imposes the strongest constraints.
We convert the modeled τ v to column density via (Savage & Sembach 1991): where λ is the wavelength of the line inÅ, f is the oscillator strength, and τ v is the optical depth in velocity space in km s −1 , which is integrated over the spatial dimension x. Using Equation 2, we compute the column density for each velocity element and integrate over the range in velocity, which we justify later in this section. In Table 2, we quote the power law estimation of the column densities (when available) as this is the most accurate and most plausibly physical explanation of non-black saturated absorption at large distances de Kool et al. 2002c). In Table 2 and throughout this paper, we will refer to this as method 1.
Method 2 is to model the Fe ii resonance lines with a velocity dependent line of sight (LOS) covering factor. The LOS covering factor in this case is a geometric step-function distribution (see Figure 4 in Arav et al. 2008). In other words, a fraction of the extended source is completely obscured by the outflow, while the rest of the source is uncovered. We assume that each radial velocity element has an independent LOS covering factor, which we again determine using the Fe ii resonance lines. We solve a fitting equation for τ v (Arav et al. 2005): where I is the normalized intensity within the trough and C(v) is the velocity dependent LOS covering factor. First, we solve Equation (3) for C(v) and τ (v) using at least two resonance lines from Fe ii. We then calculate the velocity dependent column density using τ (v) in Equation (2) and multiply the column density at each velocity by C(v) to obtain the average column density across the line of sight. Finally, we compute the total ionic column density by integrating across the velocity range of the trough. This method is comparable in accuracy to method 1 , and the column density results are approximately the same as found by method 1.
The last of the three methods, method 3, is simply the apparent optical depth method (AOD, Savage & Sembach 1991). This method assumes that all photons traveling to the observer pass through the same amount of gas at a given velocity, which gives the equation for normalized intensity: I(v)=e −τ (v) . We again convert to column density via Equation (2) and integrate to find the total ionic column density. This method does not account for partial LOS covering and saturation effects in line centers, which can underestimate column densities by up to an order of magnitude. However, we find that for SDSS J0318-066, in most cases a LOS covering factor of 1 (i.e. the AOD approximation) yields similar N ion values as those derived by the power law or partial covering method, a point we will discuss later in this section.
For the upcoming photoionization modeling, we adopt the ionic column densities as determined by method 1 (power law method) with some exceptions to this due to line blending that are noted in Table 2. We now present our analysis for the individual components.

Measurements of Component i
We begin our column density determination by fitting the Fe ii resonance lines from the strongest component, (i). In Figure 3, we show the fits for three Fe ii resonance lines using the power law method. We select these particular line transitions because they provide us with a range of oscillator strengths and are free of blending or telluric contamination. Using these lines, we derive velocity dependent column densities that we integrate for the total ionic column density.
Unfortunately, all strong Fe ii lines from the lower energy level 385 cm −1 are heavily blended in the UV 1, 2 or 3 multiplets with lines from both the resonance and other excited lines (e.g., Fe ii* λ2396, shown in Figure 2), which renders column density measurements based on these transitions to have large systematic errors. However, we find relatively strong lines outside of the UV 1, 2 or 3 multiplets for the 385 cm −1 energy state that are useful such as the Fe ii line λ1622 from the UV 38 multiplet. We also find two weaker lines from this energy state at λ1619 and λ2281 that are unblended. The oscillator strength values for these three lines are taken from Fuhr & Wiese (2006). Of the three lines, both the λ1618 and λ1622 are reported to have accurate oscillator strengths (rated as "A" 2 by Fuhr & Wiese), while the weaker λ2281 has a reasonably reliable oscillator strength (rated "B"). However, we find that using the oscillator strengths reported by Fuhr & Wiese to be problematic. While the two weaker lines are well fit with these oscillator strengths, the trough depths relative to λ1622 are unphysical. For example, scaling the λ1619 trough in apparent optical depth by the ratio of oscillator strengths (∼1.8) should be equal to the apparent optical depth in λ1622 (or larger in the case of partial covering), however the scaled optical depth of λ1619 is significantly smaller than that of the λ1622 trough. This means that at least one of the oscillator strengths could be less reliable than previously believed. Because of this issue, we quote an ionic column density for the 385 cm −1 energy level that is the mean of the fit to the λ1622 line alone and the fit to the two weaker lines with an uncertainty equal to the range in column density. The data are presented in 7 km s −1 bins with a solid line (statistical observational error bars for these data are also shown), while the model fit is shown by a dotted line. Each line is color coordinated with the wavelength located at the lower left of the plot.
The other excited state lines from Fe ii (e.g., arising from 668, 863, and 977 cm −1 ) are also subject to heavy blending in the multiplets, but unlike the resonance and 385 excited state lines there is not sufficient column density to provide deep troughs with high signal-tonoise for the weaker lines. Thus, the column density measurements from these energy levels are significantly more uncertain. The exceptions to this are lines from the lower energy level state 1873 cm −1 ( 4 F) and 7955 cm −1 ( 4 D).
We find 5 separate absorption troughs of Fe ii arising from the level at 1873 cm −1 , which are shown in Figure 4. We successfully fit these using the oscillator strengths from Kurucz & Bell (1995) to determine the column density for this level. However, the compilation of recommended oscillator strengths of Fe ii from Fuhr & Wiese (2006) yields values lower by ∼ 0.3 dex for two of the transitions, λ1637.4 and λ1702.0. For these transitions the recommended values are rated "C", while no oscillator strengths are given for the other two observed transitions. In addition to the difference in the absolute values, there is also a difference in the relative values, and in fact our fits to the absorption troughs in the spectrum favors the relative values of Kurucz & Bell (1995). Thus, we opt to use the relative oscillator strengths from Kurucz & Bell (1995) to fit all four troughs, but the level column density derived this way is increased by 0.3 dex, as required by the absolute oscillator strengths of

Fuhr & Wiese.
We find that the results from all three models (Power law, covering factor, and AOD) are similar for the Fe ii column density measurements of component i, which suggests that the LOS covering factor of the system is close to 1. Further evidence for a LOS covering factor≈1 case is given in Figure 5, where we show that we can use the Al ii template to fit the Al iii doublet. We scale two Al ii templates (one for Al iii λ1855 and one for Al iii λ1863) in τ space for each individual kinematic component for the blue and the red line of the doublet by the expected ratio of 2:1 and coadd the templates to generate a fit to the data. There exist only small deviations, and the overall fit to the doublet is good. The most obvious and important deviation is in the fit to the λ1863 i component, which is due to some saturation in the core of the Al ii template. This demonstrates that AOD is an adequate approximation and the LOS covering factor for all components is approximately 1.0. This same procedure gives similar results when applied to the Mg ii doublet as well.

The Physical Connection Between Material Seen in Different Troughs
Ultimately, we want to establish the physical characteristics of the outflow based on the column density determinations described above. In order to do so we use the integrated The blue and red tickmarks and letters indicate the component positions of the Al iii doublet members, Al iii λλ1855,1863 respectively. The dotted blue and red lines are the individual Al ii template fits for each member of the doublet (corresponding to the blue and red tickmarks). Each individual component is scaled in τ and the resulting templates match the expected doublet ratio of 2:1 for each of the components. We coadd, in τ space, the two Al ii templates and present the resulting fit in green. The fit agrees with the data with the exception of the λ1863 component i, which is likely due to saturation in the λ1855 counterpart.
column densities as input to photioinization models, which in turn determine the total hydrogen column density and number density in the outflow. This procedure depends on two assumptions: 1) that the physical conditions do not change as a function of velocity across a given trough. This assumption allow for the use of integrated column densities instead of working with each velocity resolution element separately (see Arav et al 2007); 2) That the troughs of different ions correspond to the same outflowing gas, which is a necessary condition for using these measurements in the same photoionization models. Both assumptions can be checked using the data at hand.
The first assumption can be tested by checking the ratio measured column densities amongst different lower levels of the same ion vary as a function of velocity. In figure 6 we plot the ratios of column densities between those arising from excited levels to those from the ground state level of Fe ii and Si ii. Within the uncertainties, the ratios across the troughs are velocity independent, showing that the electron number density is constant across the flow (see § 4 for elaboration). This result gives validity to our first assumption and justify the use of integrated column densities for component i. We also tested this for Ni ii/Ni ii* and combinations between Ni ii, Si ii and Fe ii resonance and excited lines and find that the ratio is constant in all cases. This constant ratio is also seen in the results for the outflow system in QSO 2359−1241 as described by Korista et al. (2008).
In order to validate the second assumption we need to check the column density ratio across troughs from different ions and especially compare low and high ionization troughs. To do so we use the τ (v) profile of Si ii λ1527 and attempt to fit the profiles of other troughs while using a single scale factor for each line. The results are shown in figure 7. It is evident that the scaled τ (v) profile of Si ii λ1527 fits the troughs of Si iv, Al ii, Al iii troughs and Mg ii remarkably well, lending strong support for the assumption that we are measuring the same outflowing gas. For C iv, which is the highest ionization line observed, the fit for the center of the trough is less conclusive since the trough is highly saturated at its core. However, we point out the the fit to the low velocity wing (-3950 to -3850 km s −1 ) is very good, which would be a coincident unless the C iv trough can be fitted with the same τ (v) profile of Si ii λ1527. In summary, the data support the two main assumption that we need in order to use all integrated column densities in our photoionization models. Si ii λ1527 data has been converted to τ (v), where τ (v)=-ln(I(v)), and scaled by a single factor for each of the other lines in order to produce the fits (red histograms) for these troughs (data shown in black). Some discrepancy between the Si ii λ1527 line and Al iii λ1863 occurs near 4000 km s −1 , which is due to contamination from the blue k component. These agreements show that the ratio of column densities between each ion and Si ii is constant across the profiles of the different troughs. i and the Al ii template for all other components. We add two templates in a similar fashion as for Figure 5. These fits, however, are scaled by one constant for the blue template and scaled by 2:1 for the red template. We also include the Si ii λ1527 resonance line and λ1533 excited line in the fit for C iv, which affect velocities greater than ∼ −5500. These 4 fits clearly demonstrate that for each ion, all 11 components are present, and that no unexplained components appear. We note that in the C iv fit there is evidence for a component at −3600 km s −1 with a large optical depth whereas the optical depth of this component Al ii appears to be significantly lower, thus the template fit is poor. The fits would improve by scaling each component by a different factor in optical depth.
To demonstrate the tie between the high and low ionization troughs for the other outflow components, we show in Figure 8 that the combination of the Si ii composite template (described in section 3.3 for component i) and the Al ii template (for all other components) can be used to recreate the doublet troughs of Mg ii, Al iii, Si iv and C iv. Similar to Figure 5, we add two templates together in τ (v). We scale the entire blue template (all 11 components) by one factor and scale the respective red template by 2:1 in τ (v). This method gives the fewest possible constraints for the fit. We find that all 11 kinematic components are present in both low and high ionization troughs, and that there are no unexplained kinematic components in the higher ionization state lines, which is indication that this is the same outflowing gas for both the high and low ionization troughs. We also note that these fits could be better if we scaled each individual component for the blue member and kept the 2:1 ratio to the red, as we did for the Al iii fit in Figure 5.

Measurements of the Other Components
The other two components with a sufficient number of measurable absorption lines are the k and a components. However, we are unable to utilize methods 1 or 2 due to a lack of unblended Fe ii resonance lines. For component k, only one Fe ii resonance line (Fe ii λ1608) is measurable and unblended. Thus, we measure the column densities based on a LOS covering factor of 1. We further confirm this assumption of a LOS covering factor of 1, as with the i component, by checking that the λ1608 line scales well to the unblended portion of Fe ii λ2383 according to oscillator strength. Component a has no measurable Fe ii resonance lines. However, according to the Al iii doublet ( Figure 5), we find that a LOS covering of ∼1 holds for all components and we are able to resort to the AOD method for our measurements. The results for components k and a are shown in Tables 3 and Table 4 respectively.
Due to the lack of measurable, unblended troughs for components b through h, we use the apparent optical depth method, and integrate the total column density across the combined blended troughs for components b through h, −7150 to −4600 km s −1 in velocity, for the Si iv, C iv, Mg ii, and Al iii doublets as well as the resonance and excited lines of C ii (λ1335, λ1336) and Si ii (λ1527, λ1533). We present these results in Table 5.

Estimates and Limits
For situations where blending between components in doublets is prevalent (i.e. C iv, Si iv, C ii, and C ii*) or where the trough is heavily saturated (i.e., C iv), our primary method of column density extraction fails. Thus, we resort to other, less accurate, methods to estimate the column densities, which are important when constraining photoionization models. We develop three additional methods in order to estimate column densities in blended regions, which we explain here in order of decreasing accuracy.
a Method used to determine column density, described by the respective label number; 1: Power Law fitting, 2: Velocity Dependent Covering Factor fitting, 3: Apparent optical depth fitting, 4: "wing" fitting, 5: "tip" fitting, 6: Upper limit based on the noise method b Number of lines used in fitting. See Online Table for specific transition information.
c Column limit based on the limiting ratio "tip" fit for Si ii*/Si ii  In the first of our estimation methods (4), we opt to employ a template to map out approximate profiles for the troughs. We utilize two templates for this process. The first template is a fit to the Al ii λ1671 line for all 11 components. However, due to moderate saturation in the core of the Al ii i component, we create a composite template using the Si ii resonance lines for the i component.
To create this template, we use both the Si ii λ1808 and λ1527 lines. The λ1527 line is similar in shape to the Al ii λ1671 line, and while the core is moderately saturated, it provides us with a good fit to the wings of the trough. In contrast, the λ1808 line has a significantly weaker oscillator strength and therefore the core of this trough is not saturated, however the wings are buried in the noise. We therefore create a composite template using the AOD column densities of the λ1527 for the wings (-3820 to -3990 km s −1 and -4390 to -4580 km s −1 ), and the AOD column densities of the λ1808 line for the core (-4390 and -3990 km s −1 ). We then convert the resultant Si ii column density profile to τ (v) template for the λ1527 line.
To measure column densities, we scale this composite τ (v) template to match the red wings of the Mg ii, Si iv, and C iv doublets as well as C ii*, which are apparently free of blending. We assume that the column density ratio is constant across the troughs, which is supported by Figure 7, scale the composite template for the i component to match and then integrate the total column (assuming a C los ≈1) from the template. We refer to this method as the "wing" fitting method (method 4). The error estimation is provided by using two τ scaling values, one underestimating and one overestimating the noise level in the trough's wing. We note that these measurements use only one line, the red doublet member, which is reflected in Tables 2 through 5. To further test the validity of this method, we used the composite template to fit the wing of the Si ii λ1527 line and we find that the template is able to reproduce the column density we find in our power law fitting of the λ1808 line to within 10%.
The next method we designate as the "tip" fitting method, or method 5, where we estimate the column densities in a similar way to method 4. For example, we estimate the column density of O i λ1302 for the i component in the region containing blended lines from Si ii λ1304 and Si ii* λ1309. We find that the composite template fit for the i component agrees well with the deep trough at ∼3809Å in velocity space (shown in online figure 1), which appears to be the tip of O i absorption. Thus, we use the composite Si ii template to approximate the minimum and maximum amount of column density that could come from O i. Due to the uncertainty in continuum placement and possible saturation, we consider this method less accurate than method 4, which is reflected in the larger error bars.
Our last column density estimation method (the noise method or method 6) defines upper limits on the column density of ions based on the statistical error of the data. Because the column density of He i is important for photoionization models (see Arav et al. 2001;Korista et al. 2008), we estimate an upper limit for the column density of the excited 2 3 s level of He i (He i* hereafter) in component i for He i* λ3189. To do this, we scale the Si ii template for the i component and compare it to the region where the line would be detected. We vary the scaling to the noise level, find the largest possible optical depth that could be buried without detection and integrate across the template to determine the column density. We also perform this analysis for Mg i λ2853, Fe i λ2484, and C i λ1560; and quote the measured limits in Table 2.
Component k presents its own problems. Unlike the i component, here we effectively measure column densities for Fe ii resonance and excited levels (E=0 cm −1 to E=668 cm −1 ), but we lack a Si ii resonance line to provide a ratio to the Si ii λ1309 excited line. Therefore, we apply method 5 with the Al ii template to find limits on the Si ii. Moreover, we use the "wing" fitting method (method 4) with the same Al ii template to estimate the column density for C ii and C ii* in component k.
In component a, the number of measurable lines is even fewer than for component k. The only clean, unblended density diagnostic is C ii to C ii*. As the critical density for this ratio is low (∼10 2 cm −3 , Blum & Pradhan 1992;Froese Fischer 1983), it is important to find another density check. Thus, using the Al ii template, we estimate an upper limit to the column density ratio of Si ii* to Si ii to be 0.49 for component a.

Systematic Uncertainty
We recognize that there exist systematic errors that could contribute further to the uncertainty in our measurements. For example, our placement of the continuum, defined in Section 2, could result in increased errors of the Si ii λ1808 and Si ii* λ1817 column densities by 5%. The amount of saturation in the core of Si iv is another systematic error that could be a potential source of error in column density, which at maximum, could lead to an error in ionization parameter of 0.1 dex. Both of these are negligible compared to the uncertainties due to the chemical abundance choice.

Spectral diagnostics and modeling
The observed spectrum of SDSS J0318−0600 exhibits distinct absorption from C ii and C iv, O i, Mg ii, Al ii and Al iii, Si ii and Si iv, Cr ii, Fe ii, Co ii, Ni ii, and Zn ii for component i. In addition, we measure upper limits to the column densities of He i*, Fe i, C i and Mg i. As already shown, a similar absorption trough structure is shared by the strongest transitions of the singly ionized species, and those of intermediate ionization (Al III and Si IV), as well as C IV. Additionally, component i likely has the largest column density of any component in SDSS J0318−0600, given that it has the largest number of absorption line identifications and is the deepest absorption feature for transitions of low oscillator strengths and/or low relative elemental abundance. Therefore, we analyze component i here and analysis of all other absorption components will be presented in a subsequent paper (Bautista et al. 2009, in prep).
One of the main objectives of modeling the observed ionic column densities is the determination of the total hydrogen column density of the outflow N H and its distance from the central source R, which are key in the determination of kinetic luminosity of the outflow (see Section 6). The former is obtained directly through photoionization modeling, but the latter is not. Instead, modeling yields the so called ionization parameter defined as: where n H is the total hydrogen density, c is the speed of light, Φ H is the surface flux of ionizing photons, Q H is the rate of hydrogen ionizing photons emitted by the central object, which is defined by: where the integration range covers the range of ionizing photons, L ν is the luminosity at each frequency and h is the Planck constant. From this definition, Q H can be estimated directly from the observed luminosity (λL λ of the object and some knowledge of its spectral energy distribution (SED), see Sections 4.2 & 4.3). Thus, the distance R can be obtained from the diagnosed value of U H whenever one has an independent determination of n H , which in a cloud with fully ionized hydrogen is roughly equal to the electron density, n e .
Thus, the first step in the analysis is to determine n e from its effect via collisions on the relative populations of excited levels observed in the spectrum. Later, we discuss the nature of the SED, which is to be used in the determination of Q and as input for photoionization modeling. After this, we will proceed to detailed photoionization modeling.

Electron density determination
Under excitation equilibrium conditions the column densities of lowly excited metastable levels can be used as diagnostics of electron density, assuming the atomic parameters and excitation processes responsible for excitation are sufficiently well understood. We have measured column densities for the 0 cm −1 ground state and low-lying metastable excited levels within Si ii, Fe ii, and Ni ii.
The simplest and probably most reliable of the electron density diagnostics is that from the ratio of the column densities of the Si ii levels, as shown in Figure 9. Here, the ratio of the excited 3s 2 3p 2 P o 3/2 at 287 cm −1 to the 3s 2 3p 2 P o 1/2 ground level is calculated from three different sets of effective collision strengths, Dufton & Kingston (1991), Tayal (2008), and Bautista et al. (2009, in prep.). We measure this column density ratio from the pair of absorption lines at 1808.0Å (resonant) and 1816.9Å (excited state). As described above, we have a good deal of confidence in their measured integrated optical depths, especially their ratio, and we will thus place the greatest weight on this electron density diagnostic. We caution, however, that as the weakest amongst such pairs of transitions commonly measured in Si ii, there are moderate uncertainties in their oscillator strengths, even in their relative values (about 0.1 dex). Two other such pairs of transitions that lie within our spectrum, centered on 1307Å and 1531Å, have smaller uncertainties in their oscillator strengths but are saturated in our spectrum and are not useful for direct measurements. These uncertainties dominate over the statistical measurement errors in the column density fits to the observed troughs at 1808.0Å and 1816.9Å. Figure 9 illustrates the electron density dependence of this level specific column density ratio, wherein the solid and dotted horizontal lines represent our measurement of this ratio and an estimated range of its uncertainty, the latter dominated by the uncertainties in the transition oscillator strengths. We derive an electron density of log(n e ) = 3.3 ± 0.2 (cm −3 ).
The various observed transitions arising from low-lying metastable levels of the Fe ii ion allow for independent electron density diagnostics. However, constraining the electron density from the relative populations within these metastable levels proved more difficult. We present these in Fig. 10 for column densities of the 385 cm −1 , 1873 cm −1 , and 7955 cm −1 levels. These are the only levels for which we have sufficient confidence in their measured column densities (see Table 2) to justify their utility as electron density diagnostics. In blue are the predictions from cloudy (v06.02, Ferland et al. 1998) for photoionization models presented in Korista et al. (2008) for the main outflow system in QSO J2359−1241 (solar abundances, fixed ionization parameter and total column density for which no hydrogen ionization front formed, and n e ≈ n H within the Fe + zone). Relevant mainly for the ratios involving the 7955 cm −1 level, the electron density weighted average temperature within the Fe + zone of these models was approximately 11,000 K. In red are the predictions of the collisional equilibrium model Fe ii atom of Verner et al. (1999) for a fixed electron temperature (10,000 K) at each electron density. This difference in adopted temperature has no significant effect on the populations of the 385 cm −1 and 1873 cm −1 levels, while the population of the 7955 cm −1 level is roughly linearly dependent on electron temperature for the conditions of interest (e.g., ∼ 10% differences are found for the two temperatures quoted). The horizontal solid and dashed lines are the measured values and estimated uncertainties, respectively. The small differences in the predictions pertaining to the 385 cm −1 level and the large differences pertaining to the 1873 cm −1 level are primarily due to differences in the radiative transition rates (A-values) for forbidden transitions used in the two sets of compu-  Bautista & Pradhan (1998, red), by full Cloudy photoionization models dominated by collisional excitation Verner et al. (1999, blue), and a model that accounts for C iv Bowen fluorescence (green). The measured column density ratios are depicted by horizontal solid lines, while the corresponding errors are indicated by dashed lines. We note that the red curve for the N (Fe ii* 7955)/N (Fe ii 0) ratio is located underneath the green line. The electron density predicted by the Si ii lines (log(n e ) = 3.3 ± 0.2 cm −3 , see Figure 9) and that from the 385 cm −1 level in Fe ii (log(n e ) ≈ 3.9 cm −3 ) differ somewhat in value. However, considering the moderate uncertainties in the relevant atomic data of the Si ii and Fe ii 385 cm −1 transitions available to us as diagnostics of electron density, and the larger uncertainties in the measurements in the Fe ii 385 cm −1 troughs, we do not consider these differences to be significant.
In contrast, the predicted electron density based on the relative population of the Fe ii a 4 F 9/2 level at 1873 cm −1 , being an order of magnitude smaller, does present a possible conundrum. Among all of the column density measurements of transitions from excited states of Fe ii we believe this one to be the most secure. However, the atomic data (A-values and collision strengths) involved in the computations of level populations of the a 4 F term carry greater uncertainties than those pertaining to the levels of the a 6 D (ground) term. This is because the a 4 F and a 6 D terms belong to different electron configurations, i.e. 3d 7 and 3d 6 4s respectively, and the energy difference between these terms is comparable in magnitude to the relativistic corrections to the Hamiltonian. Thus, the computation of accurate transition rates among levels of these terms requires fine tuning of large atomic configuration expansions (see Bautista 2008 for a more detailed discussion).
We also consider the possibility of some selective excitation process that might lead to a significantly lower population in Fe ii a 4 F 9/2 than would be inferred by collisions alone. Neither photoexcitation by the continuum radiation field nor Lyα fluorescence can explain the observed discrepancy. Another mechanism, known to occur in spectra of circumstellar envelopes of symbiotic stars and in the environs of gamma-ray bursts is Bowen fluorescence of Fe ii by C iv emission at 1548Å (e.g. Johansson 1983;Eriksson et al. 2008;Vreeswijk et al. 2007;Hartman & Johansson 2000;Michalitsianos et al. 1992). Pumping occurs as the 1548.19Å resonance line of C iv closely coincides with the 1548.20Å transition from the a 4 F 9/2 metastable level to the y 4 H o 11/2 level of Fe ii. Subsequently, this high level cascades down to the ground stage of the ion, so the population of the a 4 F 9/2 level is re-distributed among other low lying metastable levels (see also Verner et al. 2000). Fully self-consistent modeling of this fluorescence process requires a detailed treatment of radiative transfer of C iv and Fe ii lines within the cloud, which is beyond the scope of the present paper. Nevertheless, we can estimate the relative level populations of the observed levels in our spectra under the fluorescence scenario by artificially including its effect in the Bautista & Pradhan (1998) collisional equilibrium Fe ii model. We find that solving this discrepancy requires a photoexcitation rate in the 1548.20Å transition that is at least ten times greater than what is available through photoabsorption of radiation from the adopted dereddened SED (Section 4.2). We show the resulting effects on the level populations as the green curves in Fig. 10 (i.e., the red curve becomes the green one). Although, this fluorescence mechanism is able to correct the population of the a 4 F 9/2 level to agree with the observed column density at the same electron density given by Si ii, the pumping rate required seems too large even when considering the inclusion of C iv emission from the broad emission line region. Further investigation of this effect is left for future work. After considering all of these issues pertaining to the various uncertainties, we adopt the Si ii-determined electron density of log(n e ) = 3.3 ± 0.2 (cm −3 ).
Finally, we measure column densities from the ground 3d 9 2 D 5/2 level and the 3d 8 4s 4 F 9/2 metastable level at 8394 cm −1 for the Ni ii ion. We find that the population of the metastable level exceeds the predictions of pure collisional excitation conditions for any value of the electron density. This is to be expected because Ni ii is very efficiently excited by continuum radiation (Bautista et al. 1996;Lucy 1995). As a consequence Ni ii cannot be used as density diagnostic here, but instead serves as a check for the flux of photons below the hydrogen ionization threshold from the central source. As such, Ni ii can be used as a direct distance indicator to the outflow (see Section 4.4).

The Spectral Energy Distribution
The SED incident on the outflowing gas has important consequences for the ionization and thermal structures within the gas. In what follows, we provide motivation for our choice of two standard SEDs to use in the photoionization models. While AGN are routinely observed in the optical/UV and X-rays (see Elvis et al. 1994), the critical EUV spectral region (∼1 Ry to 0.1 keV, rest frame) is largely inaccessible except near its low and high energy extremes. The present working model is that a thermal accretion disk spectrum dominates the energy contained within the continuum spectrum feature known as the "UV bump", sometimes called the "big blue bump" (Shields 1978;Malkan 1983). Such a feature is expected to peak (in units of νF ν ) at energies of a few Rydberg, depending on the mass of the supermassive black hole and the mass accretion rate (Shakura & Syunyaev 1973). Mathews & Ferland (1987) had this in mind when they constructed a simplified piece-wise power-law representation of a typical QSO SED, basing the relative strength and roll-over energy of the UV-FUV bump on measurements of the strengths of the He ii broad emission lines (photon counting) and overall energy balance within the broad emission line gas (e.g., comparing the strength of Lyα 1216 to the strengths of strong cooling lines such as C iv λ1549), and joining the turn-down onto a power-law representation of the X-ray spectrum for typical X-ray bright quasars (samples of which at that time were dominated by radio loud QSOs). This spectrum, hereafter referred to as MF87, has since been the standard SED in modeling photoionized gas associated with AGN. We adopt the MF87 SED as one of two standards for our photoionization models.
By the late 1990s, improved UV and X-ray observations of AGN led several investigators to reconsider the importance of accretion disk emission to the FUV spectra of AGN. In a sample of luminous QSOs with HST spectra, Zheng et al. (1997) found surprisingly soft continuum SEDs at wavelengths shortward of ∼1000Å. Shortly thereafter, Laor et al. (1997) reported that the soft X-ray spectra of quasars appeared to be pointing (in νF ν space) at the soft FUV slopes found by Zheng et al. and proposed that a large fraction of the UV bump inferred by Mathews & Ferland (1987) and standard accretion disk models was missing (see also Elvis et al. 1994). This led to investigations of cooler accretion disk models with thermal spectra that are Comptonized by a ∼ 10 8 K plasma, which result in power-law like instead of exponential like high energy tails in the spectrum (e.g. Blaes et al. 2001). Telfer et al. (2002) followed up on the work of Zheng et al. with a larger sample of objects with HST spectra, which improved the statistics at the shorter wavelengths. They found similarly but not quite as soft FUV slopes, and concluded for radio quiet quasars "that it is plausible to represent the entire typical ionizing continuum from ∼ 10 eV to ∼ 2 keV by a single power law" (α ≈ −1.6). With this in mind Hamann et al. (2001Hamann et al. ( , 2002 constructed a simple piece-wise power law radio quiet quasar SED to model intrinsic AGN absorption systems and broad emission lines. We adopt a slightly modified version of their SED as our second standard: from 1 micron to 1 Ry, α = −0.6; from 1 Ry to 2 keV, α = −1.6; from 2 keV to 50 keV, α = −0.9 and a high energy break beyond. This SED has an α ox ≈ −1.43, very nearly equal to that of the MF87 SED (α ox ≈ −1.40), while effectively excluding the presence of a UV bump peaking at FUV energies 5 . We refer to this as the "UV-soft SED", and we show it along with the MF87 SED in Figure 11.
With regards to these observational studies of quasar SEDs, the reader should bear in mind that significant uncertainties remain in the corrections in the UV for Galactic extinction (primarily in the adopted value of R V ), in the corrections for atomic absorption within the intervening IGM (mainly the Lyα forest), as well as in the quasar host galaxy extinction. No attempt was made to correct for the last of these in any of the above works. Thus from an observational point of view, the existence and magnitude of the 1000Å break and the apparent roll-down toward the observed X-rays, as suggested in Laor et al. (1997), remains unclear. Scott et al. (2004a), in a sample of FUSE spectra of luminous Seyfert1 galaxies and low redshift quasars, reported no FUV spectral break down to ∼600Å. However, some of the differences between this finding and that in the more luminous AGN sample of Telfer et al. (2002) may be intrinsic, due to differences in the accretion disk physics. We refer the reader to Shang et al. (2005) and Bonning et al. (2007) for recent detailed comparisons between observations and modern accretion disk model predictions, as well as critical evaluations of the current uncertainties. We consider it plausible that actual quasar FUV/EUV SEDs fall somewhere between our two adopted standards.

Reddening of the SED
In the previous section, we constructed a plausible UV-soft SED for SDSS J0318−0600. However, it is known that the UV spectrum of SDSS J0318−0600 is strongly affected by extinction due to dust. Hall et al. (2002) found that if the Small Magellanic Cloud reddening curve, typical for QSOs, were adopted the extinction correction E(B − V ) would be between ∼0.1 and ∼0.4 depending on whether the 2000-3000Å or 1250-2000Å rest frame wavelength ranges were fit to the SDSS composite QSO spectrum. Because the location of the bulk of the reddening of the SED is unknown, we must consider a scenario in which the SED incident on the outflow has been reddened by dust grains. Figure 12 shows a portion of the SED centered around the rest frame energy range of our observations. Here we show the unattenuated UV-soft SED described in the previous section and the FIRST Bright Quasar Survey composite spectrum of Brotherton et al. (2001). The observed spectrum of SDSS J0318−0600, also shown in the figure, is clearly affected by extinction. Five reddened versions of the UV-soft SED are also shown here; two of them are a result of the standard extinction curve of the Small Magellanic Cloud (Hutchings 1982) with E(B − V )=0.1 and 0.4, and additional curves generated from cloudy models of pure silicate grains with a distribution of the form These models were simulated with cloudy as the transmitted spectrum of a cloud with a Galactic dust to gas ratio, log n H (cm −3 ) = 4.0, log N H (cm −2 ) = 21.8 and log U H = +2.0. This high ionization parameter was chosen in order to minimize atomic features in the resulting spectrum and to allow grain extinction to dominate the attenuation. A cloud of this number density and ionization parameter would be located ∼ 30 pc from the central source roughly with the expected realm of the putative clumpy, dusty torus. The dust temperature in this cloud has a maximum value of 400K, which is below the sublimation temperature of silicate grains (∼ 1400 K for silicates). See Pier & Voit (1995); Nenkova et al. (2008) for a more extensive discussion of the properties of illuminated dusty tori associated with AGN. No graphite grains are included in the models because they lead to the 2175Å peak, which is absent in the spectrum. On the other hand, a large absorption peak from silicates at approximately 690Å (∼1.3 Ry) does correlate with the observed spectrum. The depth of such a peak is dependent on the size distribution of grains (i.e. in the value of α). A problem arises in determining the values of α and total extinction that are proportional to C, because from the UV spectrum alone no unique solution for both values can be found simultaneously. Fortunately, SDSS J0318−0600 was also observed by the Two Micron All Sky Survey (2MASS). The photometric measurements in the J, H, and K bands (rest frame 4162, 5602, and 7277Å respectively) are marked in Fig. 12. It is important to note here that the flux in the J band may be affected by Hδ and other emission lines in the rest frame of the quasar. With this large wavelength range we are able to determine the magnitude of extinction, while the value of α is determined by fitting the slope of the SDSS spectrum.
We find E(B − V ) = 0.15 for an extinction curve with a ratio of selective to total extinction The extinction curve of SDSS J0318−0600 with R V = 4.4 is quite different from that of the SMC with R V ∼ 2 and the diffuse Galactic Interstellar Medium (ISM) with R V = 3.1. The steepness of the present curve is more similar to that of dense molecular clouds and star forming regions with R V ≈ 5. These large values of R V are characteristic of a distribution of grains that favors larger grains than in typical ISM. This would result if either large grains are allowed to accrete in a dense environment shielded from UV radiation or small grains are efficiently destroyed in a violent event, e.g. shocks. Such may be the environment of the dusty clouds composing the clumpy tori of AGN. Once the extinction that affects the observed spectrum has been determined we are in the position to compute the bolometric luminosity of the QSO. To do that we adopt the UV-soft SED shown in Figure 12 of the previous section and integrate over the whole energy space. This gives L Bol = 5.0 × 10 47 ergs s −1 , assuming cosmological parameters of H 0 =70 km s −1 Mpc −1 , Ω m =0.3, and Ω Λ =0.7.
With the above normalization and our adopted (UV-soft) SED we find a hydrogen ionizing photon rate of Q = 3.6 × 10 57 photons s −1 . However, it is unclear at this time where most of the extinction occurs with respect to the outflow. If the grain extinction of the SED occurs interior to the outflow, Q should be determined from our preferred dust reddened model (α = 3.4 and E(B − V ) = 0.15). For this case, the effective Q = 2.1 × 10 56 photons s −1 . We will later use Q in conjunction with the derived values of gas number density and ionization parameter to estimate the distance of the outflow from the central source.

Photoionization modeling
We start by assuming constant total hydrogen density "clouds" of solar abundances from Lodders (2003). The optimization of all our models is based on the column densities of Si ii and Si iv. This is because they have constrained column density measurements from the spectrum, expand the whole range of ionization stages observed, and they come from the same chemical element, thus their relative strengths are independent of the chemical composition of the cloud.
In trying to reproduce the measured ionic column densities we have constructed a number of photoionization models exploring the physical conditions and our assumptions. For a qualitative description of the ionization structure of the cloud see Korista et al. (2008). It is important in all models that the average electron density predicted by the model coincide with that diagnosed from the Si ii lines in §4.1 within 0.1 dex. For every model, we compute the electron density averaged over the Si ii trough, as: This value is presented in Table 6 together with the adopted n H .
A summary of the results of different models is presented in Table 6. Here the first three column densities indicate the species observed, their measured column densities and respective uncertainties. The subsequent column densities present the differences between observations and the results of various models as described below. Here we also present the optimal U H and N H determined for every model.
We first try to model the cloud using solar abundances and two possible SEDs: our UVsoft SED described above and the UV-soft SED reddened as in the previous section. While the measured silicon column densities can be matched in both cases, this only happens when going very deep into the ionization front where n e and temperature have started to drop. Therefore, we see that in order to reproduce a log(< n e >) = 3.3 in the Si ii region the particle density of the cloud must be log(< n H >) = 3.8. Another concern about these models is their predictions for O i column densities, which come to be about two orders of magnitude greater than allowed by observations. These results are shown as models A and B in Table 6 for the unreddened and reddened SEDs respectively. Because O i is tightly driven by H i through charge exchange, the observed column density of O i is a direct indicator of the depth of the ionization front with respect to the total column density of the cloud. The same argument is true for the column density of Fe ii. The problems with this model are illustrated in Figure 13. Despite the fact that our measurement of O i is poor the discrepancy between this and models with solar composition is significant. The current measured column density is 1.9×10 15 cm −2 with an uncertainty of a factor of 2−3. On the other hand, models with solar composition overestimate the column density by 1.9 dex, i.e. 1.5×10 17 cm −2 . This would be the largest column density of all ions in the spectrum, about 5 times that of C iv. Despite the oscillator strength of O i being smaller than that of C iv, a column density as predicted would yield a very strong and completely saturated trough in the spectrum. This is clearly not the case. The fact that the present model goes too deep into the ionization front to reach the high observed ratio of Si ii to Si iv column densities indicates that the abundance of Si relative to O must be higher than solar.
Similar to Si, the Fe/O and Ni/O abundance ratios must also be supersolar because the ionic fractions of Fe ii, Ni ii, and O i are tied to neutral hydrogen through charge exchange. Thus, the Fe ii/O i and Ni ii/O i ratios are nearly independent of the physical condition of the cloud and can only be varied by changing the abundances of these elements (Bautista & Pradhan 1998;Bautista et al. 1996). Furthermore, we find that for solar abundances of all other elements, the minimum enhancement of Si, Fe, and Ni in gas phase that satisfies the are We note that this abundances pattern is similar to that expected in high metalicity bulges of galaxies. Ballero et al. (2008) studied the evolution of chemical abundances of spiral bulges hosting Seyfert nuclei by modeling the efficiency of star formation. They predict the abundances of various elements for Z/Z ⊙ = 4.23, 6.11, and 7.22. We present these abundances, which we use for photoionization modeling in Table 7. We present these elemental abundances for the first and last of the metallicities, and use them in the photoionization modeling to investigate the effects of a realistic span in gas metallicity. Ballero et al. (2008) do not predict abundances for Al, which is present in our spectra, thus we simply scale this linearly with Z to [Al/H]=+0.50 and 0.86 at Z/Z ⊙ = 4.23 and 7.22, respectively. We use the [He/H] vs. Z relation from those predicted in the chemical evolution models of (Hamann & Ferland 1993). The remaining elements through Zn untabulated by Ballero et al. were assigned abundances as estimated by Z. Model A: UV-soft SED and solar abundances; Model B: Reddened UV-soft SED and solar abundances; Model C: UV-soft SED and Ballero et al. (2008) abundances for Z/Z ⊙ = 4.2, see Table 7 Model D: Reddened UV-soft SED and Ballero et al. (2008) abundances for Z/Z ⊙ = 4.2, see Table 7 Model E: UV-soft SED and Ballero et al. (2008) abundances for Z/Z ⊙ = 7.2, see Table 7 Model F : Reddened UV-soft SED and Ballero et al. (2008) abundances for Z/Z ⊙ = 7.2, see Table 7 a Determined from the population ratio of Si ii excited (287 cm −1 ) to ground levels 2 abundances, also with a reddened SED and log U H =-3.13, (right panels). It is found that when solar abundances are adopted the measured Si column densities require a total column density that extends beyond the ionization front, which then results in overestimated O i and Fe ii column densities The results of these models are given as models C through F in Table 6. It can be seen that these abundances yield better agreement with respect to the measured column densities than when solar abundances are used. Moreover, the abundance enhancements for Z/Z ⊙ = 4.23 still yield overestimated column densities for O i by one to two sigma, depending on the SED used. On the other hand, the models with abundances for Z/Z ⊙ = 7.2 agree well in terms of the O i predictions if the reddened UV-soft SED is adopted, but underestimate the O i column density when the non-extinguished UV-soft SED is used.
Given our arguments for metalicities greater than solar in the outflow and the fact that the solar models are unable to match the n H , the true metalicity in QSO J0318−0600 is likely to be several times solar. Regardless of this we are able to well constrain U H and N H . Models C and E constrain log(U H ) to −2.75 and −2.85 and log(N H ) to 20.16 and 19.86 (cm −2 ), respectively, for the unreddened SED. While, models D and F constrain log(U H ) to −3.02 and −3.13 and log(N H ) to 20.07 and 19.73 (cm −2 ), respectively, for the reddened SED.
One significant discrepancy between all models and observations is in the carbon column densities. As indicated in Table 7, the [Si/C] ratio is rather large in these two elemental abundance models, and a relatively larger carbon abundance would reduce the discrepancies. Ballero et al. point out that their prediction for C is likely uncertain as it depends sensitively on the mass of the bulge. However, it is important to point out that our measurements for carbon are problematic, as we only have a lower limit for the C ii column density and large portions of the C iv troughs are highly saturated.
A relatively small uncertainty in U H and N H arises from the adopted standard SED. For example, should we adopt the MF87 SED rather than the UV-soft SED constructed in Section 4.2, our results for U H and N H would be affected at a ∼0.1 level. However, this effect is very small compared to the uncertainties resulting from the dust extinction of the SED and the adopted chemical abundances.
Another discrepancy in the models stems from the predictions for Al. Most of our models underestimate the column densities of Al ii and Al iii. The predicted column density ratios of Al ii to Al iii differ from observations by approximately a factor of four. Again, it could be that the true Al abundance is higher than assumed, which is not surprising as there are only very rough approximations on how this element scales with metalicity. But, beyond the absolute column densities there may also be differences in the predicted ionization structure of Al. The other problem is that the predicted column densities of Mg i always exceed the upper levels derived from our spectra. It is clear that the ionization structure of neutrals that form by recombination across the ionization front is not well described by our models. The cause for this problem is a combination of inaccuracies in the atomic parameters for these species and limitations of photoionization modeling codes to treat neutral ions with low ionization potential (<10.2 eV). The exception to this is He i*, which is consistent with the models predictions. This problem is intrinsic to all modeling codes that treat radiative transfer through the escape probability approximation, as neutral species with low ionization potential are greatly affected by Lyman α radiation. Thus, a very accurate treatment of radiative transfer is needed to confidently model these neutral species, which is beyond the scope of this paper.
3.6. The Ni ii lines as independent distance indicators Lucy (1995) demonstrated that Ni ii could be efficiently photo-excited by continuum radiation under nebular conditions. Bautista et al. (1996) constructed a detailed 142-level model of the Ni ii system for spectral analysis, which was updated with improved atomic data by Bautista (2004). Figure 14 shows the level populations ratio between the excited level at 8394 cm −1 and the ground level. The upper panel shows the ratio vs. n e under pure collisional excitation conditions. It is noted that the uncertainty in the measured ratio column densities is dominated by that in the oscillator strengths, particularly for troughs from the 8394 cm −1 . For these transitions we adopted oscillator strengths from Kurucz and Bell (1995), which in the case of the resonant transitions were discrepant by up to a factor of two from experimental determinations of Fedchak and Lawler (1999). Thus, we adopt a nominal 30% uncertainty in the observed ratios. Further theoretical and experimental work on the oscillator strengths of Ni ii is warranted.
From Figure 14 one sees that the observed ratio is consistent with the predictions of pure collisional excitation for a wide range in n e , that includes the conditions of the absorbed as diagnosed from Si ii. Hence, it should be possible to set an upper limit to the continuum flux that also contributes to the excitation of the 8394 cm −1 level. In the lower panel of the figure we plot the ratio against the flux of ionizing photons Φ for the same SED used in our photoionization models. It is important to notice that only photons less energetic than 1 Ry participate in photoexcitation of the ion and these are unaffected by the main atomic opacities of the cloud. It is found that the N(Ni ii 8394)/N(Ni ii 0) ratio is strongly enhanced by continuum fluorescence and this drives the ratio well above the Boltzman equilibrium populations. Under these conditions, induced photoexcitation of highly excited odd parity levels is balanced by rapid spontaneous radiative decay. For much more intense radiation fields, induced photoexcitation exceeds spontaneous decay and the level population ratio reaches a plateau. It can be seen that the observed column density ratio sets upper limits of log 10 (Φ/c) ≈ +1.5 or +0.4 depending on whether one uses a reddened or unreddened SED. Thus, from a total particle density of log(n H ) = 3.3 one derives a value for the ionization parameter log 10 (U H ) ≤ −1.8 or ≤ −2.9, for reddened and unreddened SEDs respectively. Both of these values are consistent with that derived through photoionization modeling. It is unfortunate that more stringent constrains could not be derived due to the large uncertainties in the oscillator strengths for the absorption lines in the spectrum.

CalculatingṀ andĖ k for the Main Outflow Component
In the previous section, we determined the ionization parameter, number density and total column column density of the main outflow component, (i). Here we calculate the mass and energetic charateristics of the outflow. First, we determine the mass contained within the wind. For simplicity, we assume that the outflow is in a form of a partial thin spherical shell at a distance R from the central source moving radially at constant speed v. The total mass contained in such an outflow is: where N H is the total column density of hydrogen, m p is the mass of a proton, µ=1.4 is the mean molecular weight of the plasma per proton, and Ω is the fraction of the shell occupied by the outflow. Therefore, assuming dN H /dt=0, the mass flux follows as the temporal derivative: the momentum flux is given by: and the kinetic luminosity is given by: For component i: v = −4200 km s −1 and table 6 gives the N H for each of our photoionization models. To calculate the distance R we use equation (4), where n e is determined in § 4.1, U H is taken from table 6 and Q H is given in § 4.4, log Q H (photons s −1 ) as 56.3 and 57.6 photons s −1 for the attenuated and unattenuated SEDs, respectively. In Table 8 we list the derived values for R,Ė k andṀ for each photoionization model.

The Global Covering Fraction
As shown above, we determined N H and R for the photoionization models and therefore each model is associated with a reliable estimates of the kinetic luminosity per unit solid angle. In order to obtain estimates for the total kinetic luminosity we need to address the remaining variable in equation (1), Ω, the fraction of the full sphere occupied by the outflow. There is no direct way to obtain the Ω of a given outflow from its spectrum as we only see the material along the line of sight. Statistically, C iv BALs are seen in 20% of all quasars (Hewett & Foltz 2003). This can be interpreted as: a) all quasars have BAL winds covering 20% of the solid angle, b) 20% of all quasars have winds covering the full solid angle, c) all quasars have a phase of full solid angle wind lasting 20% of their duty cycle. A combination of these 3 options is also possible. It is important to note that with all else equal, options a) and c) yield the same integrated mass flux and kinetic luminosity over the life time of an individual quasar, and that all 3 options yield similar integrated quantities for a sample of quasars.
However, we cannot measure distances to the majority of C iv outflows, but only to those outflows that show troughs from excited and metastable levels. For most of our objects (including SDSS J0318-0600), the spectra cover such transitions only from singly ionized species (e.g., Fe ii, Si ii, C ii, and Ni ii). Outflows with troughs from these ions are rare (Trump et al. 2006). Our comprehensive search of the 50,000 brightest SDSS quasars yielded 100 such outflows. We found that only ∼ 1% of the objects that show a C iv outflow show troughs from these singly ionized species. We therefore have two limits to explore with regards to the Ω of the outflow. First, outflows that exhibit troughs from these 4 singly ionized species excited states may represent all regular C iv outflows viewed through special conditions that allow for the formation of singly ionized troughs. This scenario yields Ω ≃ 0.2 as is the case for the frequency of C iv BALQSOs among all quasars. The second limit is Ω ≃ 0.002, which is the minimal covering fraction needed to explain the 100 occurrences in 50,000 objects.
Several arguments suggest that the first scenario is closer to the actual physical situation. Since outflows with troughs from singly ionized species are rare, we may simply be observationally biased by this low ionization selection effect as we lack excited and metastable level troughs from C iv or species with similar ionization. One physically motivated picture is that we are looking at a normal C iv outflow through a line of sight passing through semiopaque material at the edge of the putative AGN torus. Such a scenario explains why we see the singly ionized species (lowering of the ionization parameter); why most of the objects we measure have strong dust reddening (torus dust) and why these systems are so rare (a solid angle that has sufficient torus material to lower the ionization parameter significantly but not enough to make it a type II quasar). Aoki et al. (2009, in preparation) give support to the above scenario. A quasar showing strong Fe ii absorption troughs and strong dust reddning in 2002, showed much less dust extinction and very little Fe ii absorption in 2008. Thus, the amount of dust along our line of sight changed drastically between two epochs seperated by ∼3.5 years in the rest frame of the source. The interpretation is that the dust source has vacated the line-of-sight, therefore increasing the ionization parameter of the outflow and reducing the column densitiy of Fe ii. A 3.5 years time scale favors changes in a near-by obscuring material (i.e., the putative torus) over material at a kpc scale and therefore support both of our assertions: 1) that the dusty material seen in the spectra of these outflows is indeed much closer to the source than the outflows themselves; 2) that these excited-state singly ionized absorption outflows are naturally explained by a special line-of-sight grazing the putative torus.

Recent observations by
Further support for this scenario comes from the fact that all our objects show a C iv BAL, which is morphologically indistinguishable from C iv BALs in the ubiquitous high ionization BALQSOs. We therefore conclude that it is more likely that we are seeing normal Ω ≃ 0.2 outflows through a rare line of sight.
A direct resolution of this issue necessitates imaging of the winds. This is a difficult task due to the high redshifts of the objects and the presence of a very bright quasar at the center of the image. Nonetheless, we hope to pursue this line of research in the future. Another approach is to observe at wavelengths shorwards of Lyα where there are a few excited state troughs from higher ionization species. In particular S iv/S iv*λλ1062.66,1072.97. This ion is similar to C iv in ionization structure and indeed S iv troughs are much more prevalent than those of the singly ionized species. Measuring N H and R for these objects will therefore eliminate the uncertainty connected with the rarity of outflows showing troughs from singly ionized species, and will yield robust measurements forṀ andĖ k . We are currently pursuing this project, which is challenging due to the heavy Lyα forest contamination.

The Super Massive Black Hole and Host Galaxy
Given the bolometric luminosity and assuming that the super massive black hole (SMBH) is accreting near the Eddington limit we find that the amount of mass being accreted is: where we assumed an accretion efficiency (ǫ r ) of 0.1 (Salpeter 1964) and c is the speed of light. This large accretion rate (the outcome of a very luminous quasar) suggest that the mass flux values we derive (see Table 8) are not unreasonable. We expect such an outflow to have roughly an order of magnitude higher mass flux than the accretion rate onto the SMBH since the original outflow entrained material on its journey to several kpc distance. Simple momentum conservation requires that an outflow which started at ∼20,000-30,000 km s −1 to have 5-7 times more mass flux by the time it slowed to 4200 km s −1 . Therefore, an outflow outflow such as this could have began withṀ outf low ∼Ṁ acc .
Next, assuming Eddington accretion, we can derive the SMBH mass and infer a mass for the host galaxy bulge. The black hole mass derived from Eddington accretion is (Peterson 1997): where σ T is the Thompson scattering cross-section for an electron, G is the gravitational constant, M BH is the mass of the SMBH, and m p is the mass of the proton. Thus, the spheroidal bulge component of the host galaxy (assumed to be approximately the mass of the host galaxy) has a mass M bulge =3.3×10 12 M ⊙ based on a SMBH to bulge ratio of 0.0012 (McLure & Dunlop 2002). This mass can be compared with the total mass ejected by the outflow over the ∼ 10 8 year duty cycle of the quasar, which is a few times 10 10 M ⊙ .

Discussion
In this section we compare the results from this object to other outflows found in the literature; discuss the location of the dust in SDSS J0318−0600, which is a major cause for uncertainty in the derived mass flux and kinetic luminosity; elaborate on the relationship between the dynamical time scale of the outflow and the duty cycle of the quasar; and conclude by describing the possible role of these outflows in AGN feedback phenomena In Table 9 we compare all the analyses in the literature that can be used to obtain mass flux and kinetic luminosity for outflows in individual quasars. We also note that with the exception of SDSS J0318−0600, every absorber shows evidence for inhomogeneous line of sight (LOS) covering. For outflow absorption, we refer to homogeneity as the spatial distribution of gas across our line of sight to the quasar (see section 3.2), where a homogeneous flow would uniformly cover the background emission source (LOS covering factor of 1). Initially, this may be attributed is to the large distance of the absorber, however, the outflow in 3C 191 shows evidence of inhomogeneous covering at a greater distance than we find component i of SDSS J0318−0600.   Table 5 in Wampler et al. (1995) d Distance derived from Fe ii fluorescence and no photoionization modeling was performed for this object SDSS J0318-0600 is also the most heavily reddened in the sample, which leaves us with the question of the location of the obscuring dust. Based on the value from Schlegel et al. (1998), we know that the Milky Way Galactic reddening is minimal, and thus the reddening is intrinsic to the quasar. There exist three possible locations of the dust with respect to the absorber: interior, dispersed within, and exterior. The two locations for the dust we consider here due to their potential impact on the modeling of the outflow are that the dust is interior to the absorber or the dust is associated with the absorber. If the dust is interior to the absorber, a plausible source of dust would be the edge of the putative obscuring torus of the AGN unified model (Antonucci & Miller 1985). This is supported by the object discussed by Aoki et al. (2009, in preperation) and furthermore would provide the rare sight line needed to match the frequency of detections for low ionization BAL outflows that show iron absorption (FeLoBAL), as discussed in section 5.1. By contrast, galactic wind imaging in nearby galaxies shows that dust co-exists inside the winds and may be associated with the outflows (e.g., M82, Ichikawa et al. 1994). Elvis et al. (2002) proposed the possibility of grains forming within AGN outflows, as conditions allowed. We note though that if the outflowing gas contains the observed dust, the amount of dust depletion must be substantially smaller than in the ISM or SMC in order to produce the observe ionic column densities. Therefore, while we cannot strongly exclude the possibility of a dusty outflow, locating the dust at the edge of the putative torus is physically more plausible and explains several observed features in the FeLoBAL population: the heavy reddening of FeLoBALs, the low reddening experienced by C iv BAL quasar; and the higher prevalence of high ionization BALSQSOs compared with FeLoBALs.
The dynamical time for component i of the outflow is t d = R/v ≈ 10 6 yr (assuming R = 6 kpc). Therefore, if we are seeing a representative stage of the outflow, over its ∼10 8 yr duty cycle (Hopkins et al. 2005b) the quasar must eject roughly hundred times more mass than is seen in component i. Obvious candidates for more ejected mass are the other 10 components of the outflow. However, Bautista et al. (2009, in prep) show that these components collectively do not contribute appreciably to the total mass output for SDSS J0318−0600. Earlier episodes of mass ejection may not be detectable in absorption. Let us assume that component i is continuing at a constant velocity and no longer entraining any significant mass. Within 10 7 yr, its distance will increase by a factor of ten and therefore its optical depth will drop by a factor of a hundred, making it undetectable in absorption. To complete the picture we note that in addition to the FeLoBAL system, SDSS J0318−0600 has a high ionization BAL at ∼27000 km s −1 . This outflow component is probably much closer and therefore represent a very young mass ejection episode that in time may evolve to something resembling component i. Similar situations exist in other objects (e.g. SDSS J0838+2955; Moe et al. 2009).
Given these caveats the kinetic luminosity, such as the one we have measured in SDSS J0318−0600, provides ∼(0.1-0.6)Ω 0.2 % of the bolometric luminosity as feedback. We note that the measured mass of component i represents only the cold phase (T=10 4 K) of the outflow. The hot phase of the gas (or X-ray, T=10 6 K) has been shown in Seyfert galaxies (e.g., NGC 3783; Netzer et al. 2003) to contribute up to 50× more column density than the cold phase, therefore it is possible that the gas in the hot phase of the outflow in SDSS J0318−0600 could contribute significantly to the mass flux and kinetic luminosity. For comparison, AGN feedback models of galaxy growth and evolution require a mechanical energy input of roughly 5% of the bolometric luminosity of the quasar. For example, Di  use this ratio to recreate the observed M BH -σ relation and Hopkins et al. (2005a) show in their models that the same value is required to 'uncover' a shrouded quasar after a merger. Thus, the absorption outflow we observe in SDSS J0318−0600 can be a significant contributer for AGN feedback processes.
Beyond the effects of feedback on the host galaxy, if SDSS J0318−0600 is representative of the type of quasars found in clusters, then we also see that an outflow with the kinetic luminosity of SDSS J0318−0600 would contribute energy to the gas that forms the ICM. The models of Scannapieco & Oh (2004) show that a 5% feedback ratio would suffice to heat the ICM during formation of the cluster and provide the well known L x -T relationship (Kaiser 1986;Arnaud & Evrard 1999) Also, related to clusters are X-ray cavities in the ICM, which coincide with relativistic radio jets (e.g. NGC1275 in the Perseus cluster, Fabian et al. 2006). We suggest that the sub-relativistic outflows could be related to the bubbles or perhaps are the driving mechanisms for bubble inflation. This follows as the inflation of a bubble is a simpler process provided a wide angle wind (Sternberg et al. 2007). Rafferty et al. (2006) shows that the power needed to inflate these bubbles is ∼10 45 ergs s −1 , which is similar to the kinetic luminosity we measure in SDSS J0318−0600 of ∼1-3×10 45 ergs s −1 . These similarities in fraction of the bolometric luminosity and power output suggest that quasar outflows play significant roles in influencing their host galaxy and neighborhoods through feedback.
2. Using the Al II λ1671 trough, we divide the outflow into 11 kinematic components ranging in velocity from -7400 to -2800 km s −1 . We utilize these components to identify the corresponding absorption troughs in several lines from other ions.
3. For three of these components, we measure many ionic column densities and develop methods to estimate ionic column densities in blended and saturated regions. The Al iii doublet data show that all 11 components follow the expected doublet ratio of 2:1 in optical depth and therefore the covering of the emission source is homogeneous. Further tests for components i and k with Fe ii resonance lines demonstrate the same result.
4. The measured column densities of Si ii and Si ii*, allow us to determine a number density of 10 3.3 cm −3 for component i, which is the main and innermost component.

5.
We show that the column density ratio of Ni ii excited to resonance can be used to independently determine the flux of non-ionizing photons, which when combined with the assumed SED yields a direct indication of the distance. These values agree with that derived from Si ii.
6. We analyze the intrinsically reddened spectrum of SDSS J0318−0600, and find that it requires a different extinction curve from the known Small Magellanic Cloud (SMC) curve. We fit the spectrum with a reddened SED due to pure silicate dust and a power law grain size distribution with a power of α=3.4, which resulted in an R V =4.4. Based on our model of the intrinsically reddened SED, we estimate the bolometric luminosity (log L Bol = 47.7) and the number of ionizing photons illuminating the absorber with and without dust extinction (log Q = 56.3 or 57.6).
7. For component i, our photoionization models suggest supersolar abundances. The output total hydrogen column density and the so called ionization parameter for both cases with and without dust extinction are N H = 10 19.9 or 10 20.0 cm −1 and log U H = −3.1 or −2.7 respectively. 8. Given the gas density and ionization parameter, and an estimate of the source luminosity, we determine the distance to the outflow for attenuated and unattenuated SEDs of 5.9 or 17 kpc respectively, and determine the corresponding kinetic luminosities to be (0.7 or 2.5)×10 45 Ω 0.2 ergs s −1 , and the mass fluxes to be (120 or 450) Ω 0.2 M ⊙ yr −1 . We find that component i is the dominant source of energy and mass when compared to the other kinematic components, which are evaluated in Bautista et al. (2009, in prep).