the central engines of powerful type i active galactic nuclei ( agns ) contain gas that is moving at extreme velocities , apparent from their broad emission lines ( bels ) and , in the case of broad absorption line ( bal ) quasars , in blueshifted absorption troughs .
these spectral features arise from compact gas flows close to the accretion disk , apparent from the presence of high ionization lines and size scales from reverberation mapping @xcite .
these compact flows are probably the origin of larger - scale outflows ( e.g. chelouche & netzer 2005 ; crenshaw & kraemer 2007 ; storchi - bergmann 2010 ) and result in feedback that is believed to be important on many scales , regulating supermassive black hole ( smbh ) fueling , quenching star formation , and so driving galaxy - smbh coevolution ( e.g. silk and rees 1998 ; fabian 1999 ; sommerville et al . 2008 ; zubovas & king 2012 ) as well as influencing the intergalactic and intracluster medium and affecting galaxy and large - scale structure formation ( e.g. scannapieco & oh 2004 ; rafferty , mcnamara & nulsen 2006 ; sommerville et al . 2008 ; fabjan et al . 2010 ; wang 2015 ) . yet despite the importance of these processes , the connection between compact and large - scale outflows is not understood . in the case of the broad emission line region ( belr ) , even the basic kinematics are unclear ; independent reverberation mapping results support outflow , inflow , and orbital motion @xcite , and microlensing results have been consistent with both virial belrs @xcite and with more complex , multi - component flows @xcite .
bal quasars offer perhaps the clearest evidence of nuclear outflows in agns .
approximately 15% of optically selected quasars exhibit high velocity , blueshifted absorption features @xcite .
the outflows responsible for these features may be common to most quasars if this 15% represents a special line of sight to outflows with a similar covering factor @xcite , or a particular phase in the lifecycle of quasars .
an accurate measurement of the absorption profile of a bal quasar would go a long way to revealing the kinematics of these outflows , however this measurement is complicated by the fact that their intrinsic absorption profile can not be recovered without knowing the shape of the underlying bel that it absorbs .
these bel profiles can vary dramatically between quasars , and in the case of bal quasars may be strongly affected by a biased viewing angle .
the gravitationally lensed bal quasar h1413 + 117 ( the cloverleaf ; magain et al .
1988 ) offers a unique opportunity to simultaneously study both bel and bal flows in the same source .
the four lensed images of this @xmath1 quasar exhibit deep bal features in high - ionization lines , as well as low - ionization absorption in aliii , making it a lobal quasar .
significant microlensing has been observed in h1413 + 117 over the past 15 years @xcite . in microlensing ,
substructure within the foreground lensing galaxy causes changes in the magnification of individual lensed images .
the degree of microlensing depends on size of lensed region , and in the case of h1413 + 117 this has resulted enhanced magnification in the continuum in lensed image d compared to the broad emission lines in this component .    given a model of the microlensing magnification of the lens system
, these differential magnifications can be used to derive size constraints of emission regions .
this approach has been used to contrain accretion disk size and temperature gradient in several lensed quasars ( e.g. kochanek et al .
2004 ; bate et al . 2008 ; floyd , bate & webster 2009 ; blackburne et al .
2011 ; jimnez - vicente et al . 2014
; blackburne et al .
2015 ) , including in h1413 + 117 by @xcite ( hereafter h10 )
. it has also provided constraints on belr size @xcite and kinematics @xcite in a number of quasars , however no such constraints have been derived for the belr of h1413 + 117 .
another novel use of microlensing takes advantage of differential magnifications between lensed images to algebraically decompose the spectral contributions of different physical regions within the quasar @xcite .
@xcite ( hereafter a90 ) and h10 apply this technique to h1413 + 117 , shedding light on the true shape of the broad absorption profile , and in the case of h10 decomposing into microlensed and non - microlensed components which reveal substructure in each component and may point to a two - component outflow comprising a polar outflow and equatorial disk - wind .
in this paper we present new gemini north gmos integral field spectroscopy of h1413 + 117 , and based on the strong microlensing observed in lensed image d we derive new statistical size constraints on the belr .
we also present a new spectral decomposition , deriving belr and continuum components , as well as the intrinsic bal profiles for nv , siiv , civ , and aliii .    in section [ obs ]
we describe the observations and extraction of lensed image spectra . in section [ sigs ]
we analyse the microlensing signatures in this new epoch of data . in section [ microlensingbehavior ]
we use lensing simulations to explore the size dependence of the observed microlensing signatures and deriving size constraints for the belr . in section [ algebra ]
we decompose the spectra in terms of the broad emission line region and continuum , emphasizing the physical interpretation of the decomposed regions , and we derive the intrinsic bal absorption profile .
we discuss the interpretation of these decomposed components in the context of disk - wind outflow models in section [ disc ] . in section [ conc ]
we summarize our conclusions .
117.,width=302 ]    h1413@xmath2117 was observed on the night of april 5th , 2006 with the gmos integral field unit ( ifu ) @xcite on the gemini north telescope ( program i d : gn-2006a - q-15 ) .
3@xmath320 minute exposures were taken in one - slit mode using the b600 grating with a spectral range of 3961  to 6845  and a resolution of @xmath4 . at the redshift of h1413@xmath2117 ,
@xmath1 , this covers the ly@xmath5 through ciii ] broad emission lines .
seeing conditions were good , with a psf fwhm of 08 .
the data was reduced using the gemini iraf package v1.10 under iraf 2.14.1 .
we used calibration files taken as part of the standard gemini facility calibration unit ( gcal ) suite .
mostly , the standard gmos ifu reduction procedure was followed for overscan subtraction and trimming , bias subtraction , flat fielding and response correction , arc calibration , fiber extraction , sky subtraction , flux calibration , cube reconstruction with atmospheric dispersion correction .
additionally , we corrected the gcal flat and the twilight flat for scattered light by fitting a surface to residual light detected between the fiber bundles .
that correction was not necessary for the science observations , the signal being too weak to generate any significant scattered light .
the cosmic rays were corrected before extraction with la - cosmic @xcite , followed , after extraction , by a floor and ceiling masking of remaining extreme hits . the calibrated cubes for each science observation were aligned and combined with a customized version of the pyfu software package .
figure [ 2dimage ] shows the reduced data cube collapsed into a 2-d image .
all lensed components are distinct and well - resolved .
the lensed images are well resolved in these observations , allowing for a good extraction of their spectra .
we perform this extraction using the method described in @xcite .
we refer the reader to that paper for a detailed description , however in short : relative lensed image fluxes at each wavelength bin are determined by least - squares fitting of reconstructed point spread functions ( psfs ) .
a separate psf is determined for each 100 - wide spectral bin .
each given psf is initially estimated using carefully - masked fragments of the ifu frame ( stacked over 100 ) , and are iteratively improved with a series of least - squares minimizations of each frame .
figure  [ spec ] shows the extracted spectra of the four quasar images .
we see distinct broad emission lines for ly@xmath5 , nv , oi , cii , siiv , civ , and ciii ] as well as prominent iron emisison features at @xmath01600 - 1800 .
emission from heii and aliii is also apparent , blended with iron and ciii ] respectively .
clear redshifted broad absorption troughs are seen for ly@xmath5/nv , siiv , civ , and aliii .
gravitational lensing is achromatic , and so in the absence of microlensing , differential extinction , or time - dependent spectral variation , the spectra of lensed images of a multiply - imaged quasar are identical modulo a scaling factor .
microlensing is typified by a difference in the magnification of one or more lensed images compared to the time - averaged magnification for that lensed image ( which we will refer to as the macrolens magnification , or just macro - magnification ) .
_ differential microlensing _ is characterized by magnifications that vary between different emission regions , resulting in a change in the shape of the spectrum of the affected lensed image .
differential microlensing is sensitive to emission region size and projected location , is frequently seen between bels and continua , and also often shows wavelength dependence within these features .
spectral variability coupled with different light travel times can result in spectral differences between lensed images . in the case of h1413 + 117
the time delay is potentially significant , estimated , for example at 23.4@xmath64 days between images a and d @xcite .
this variability can potentially lead to differences in four key observables between lensed images : bal and bel profile shapes , continuum - to - bel flux ratio , and the spectral index of the continuum .
while it is difficult to extricate the differences due to intrinsic variability from those due to microlensing , we can comment on the effect of this on our analysis .
variation in broad absorption and emission line profile shapes would be very difficult to discriminate from microlensing .
fortunately the effect is probably minimal .
as we show in section [ fluxratios ] the bel profiles are identical between lensed images , precluding both differential microlensing and intrinsic variability in the belr spectrum .
it is not so easy to compare bal profiles , however we note that , while this profile has been observed to vary , this variation occurs over months and years , not days ( h10 ) .
as we will see in section [ microlensingbehavior ] we estimate the size of the belr ( and hence also the bal wind ) at @xmath7100 light days , and so their intrinsic variation timescales should be significantly longer than the time delay .
variation in the continuum - to - bel ratio between lensed images can not be ruled out , and we do see significant variation in this ratio , which we attribute to microlensing because it is consistent with the microlensing behavior of the lens in recent years .
in section [ algebra ] we use this difference in continuum - to - bel ratios to decompose the spectrum into intrinsic components .
it s important to note that for the purpose of the decomposition of spectral components , it does not matter why this ratio varies between images  whether by microlensing or intrinsic variability  the decomposition is still valid as long as the bels do nt show wavelength - dependent variations .
variation in the slope of the continuum could also be due to intrinsic variability , and we do indeed observe a difference in slope ( sect .
[ fluxratios ] ) , however we do nt investigate this in detail in this paper .
differential extinction between lensed images will not affect the continuum - to - bel ratios as our measurement of the continuum strength is made at the relevant bel wavelength range .
nor will it significantly affect broad line profiles , which appear identical between lensed images anyway .
differential extinction only affects the continuum slope and ratios between bels , which are not a focus of his work .          to inspect differential microlensing properties of h1413 + 117 we first obtain an estimate of the unabsorbed continuum spectrum for each lensed image by performing a linear fit to regions of the spectrum least affected by emission or absorption ( fig .
[ spec ] shows fits for images a & d ) .
prominent iron emission features force us to choose small windows at rest frame 14201460 , 16851692 , and 17951805 . while the continuum is unlikely to be linear over the full range of this spectrum
, we find that the fit gives a very plausible estimate , particularly redward of the nv line .
subtraction of the fitted continua allow us to measure the average flux ratios across each bel for each lensed image with respect to image a. these ratios and the average flux ratios for the continuum at each line center are given in table  [ fluxratiotable ] .
for the purpose of inspection of the microlensing properties , we find that changing the continuum level and slope by allowing variation of @xmath65% ( the width of the pale gray region surrounding fits in fig .
[ spec ] ) does not significantly affect results .
lcccccccc & & & & & + & line & continuum & & line & continuum & & line & continuum + nv & 0.78 @xmath6 0.02 & 0.92 @xmath6 0.02 & & 0.93 @xmath6 0.03 & 0.81 @xmath6 0.02 & & 0.45 @xmath6 0.03 & 0.77 @xmath6 0.03 + siiv line & 0.82 @xmath6 0.02 & 0.91 @xmath6 0.01 & & 0.93 @xmath6 0.02 & 0.79 @xmath6 0.01 & & 0.38 @xmath6 0.02 & 0.75 @xmath6 0.02 + civ & 0.81 @xmath6 0.01 & 0.91 @xmath6 0.01 & & 0.85 @xmath6 0.01 & 0.77 @xmath6 0.01 & & 0.41 @xmath6 0.01 & 0.73 @xmath6 0.01 + ciii ] & 0.84 @xmath6 0.02 & 0.88 @xmath6 0.02 & & 0.87 @xmath6 0.03 & 0.71 @xmath6 0.02 & & 0.38 @xmath6 0.03 & 0.68 @xmath6 0.02 + fe complex & & & & & + mean bel & & & & & + mean cont . & & & & & +    figure  [ linesratios ] ( upper panels ) compares the spectra of all lensed components around the broad emission lines . to highlight the effect of differential microlensing
we scale the spectra of lensed images b , c , & d so that the average continuum emission at each line ( from the continuum fits ) matches that of image a. we see significant differential microlensing between broad emission lines and continuum in lensed image d versus all other images , with the emission lines at lower magnification in d relative to its continuum .
other lensed images show weak differential magnification between lines and continuum . due to an unknown lens redshift
the lens model is not well established and so the macro - magnifications are not well known .
for this reason we do nt know which , if any of a , b , or c are free of microlensing .
nonetheless , we can conclude that image d is undergoing significant differential microlensing between belr and accretion disk .
these results are consistent with observations over the past decade ( angonin et al .
1990 , stensen et al .
1997 , monier et al .
1998 , anguita et al .
2008 , h10 ) , which show strong anomalous line - continuum flux ratios in image d compared to a , b , and c.    our observations show wavelength - dependent differential microlensing of the continuum in the d / a and c / a ratios .
the d / a ratio tends to more extreme values ( deviating from the average bel ratio ) towards shorter continuum wavelength , as was also observed by h10 .
it is tempting to interpret this as a simple inverse correlation between emission region size and microlensing intensity , however we note that the observed c / a ratios have the opposite behaviour , tending towards the average c / a bel ratio towards _
longer _ wavelengths . in section [ microlensingbehavior ]
we show that we can not expect a monatonic trend of magnification ratio with size across multiple orders of magnitude in emission region size .
proper analysis of differential microlensing across the continuum requires careful lensing simulation , however we do not perform that analysis in this work .    to investigate the microlensing state of the bels we determine the wavelength - dependent flux ratios for the bel seds for lensed images a and d. these are just the ratios of the spectra of images a and d after subtracting the continuum estimates described above . in figure  [ linesratios ]
( lower panels ) we can see that the flux ratios for the nv , siiv , civ , and ciii ] broad emission lines are flat over the wavelength range redward of bal absorption .
ly@xmath5 is omitted due to its overlapping the nv absorption trough , which makes continuum subtraction impossible .
this flatness indicates that there is no differential microlensing across the belr velocity structure redward of the bal absorption , and very little difference between these emission species .
the mean and individual bel flux ratios are given in table  [ fluxratiotable ] .
error intervals include standard error of the mean and uncertainty in the continuum fits .
variation of the level of continuum subtraction up to its estimated uncertainty at the location of each emission line does not change the observed flatness of these flux ratios .
in @xcite we analyzed microlensing signatures in the civ and ciii ] bels of q2237 + 0305 , the einstein cross . in that case
there was clear differential microlensing across the velocity structure of these broad lines , with magnification ratios varying by @xmath025 - 30% between line center and wings .
the bel magnification ratio approached the highly anomalous continuum magnification ratio at the highest velocities .
our simple model indicated that this trend was consistent with the broad line flow being a keplerian field , and not supporting a pure radiation - driven outflow without a substantial gravitational component .
this differential microlensing signature sharply contrasts the wavelength - independent bel magnification ratios seen in h1413 + 117 ( sect .
[ sigs ] ) .
this wavelength independence suggests that the belr in h1413 + 117 is large enough to average over fluctuations in the magnification surface .
given the very low redshift of the lensing galaxy in q2237 + 0305 ( @xmath8 ) and its commensurately large fluctuation scale - length , differential microlensing of its belr is not surprising . to explore the implications of this difference in more detail
we compare the expected microlensing behavior of these two quasars as a function of source size .
although the lensing galaxy of h1413 + 117 does not have a confirmed redshift , it is possible to explore plausible microlensing scenarios using estimates from literature : z@xmath9 = 0.94 based on the probability distribution from lens parameters @xcite and supported by the z@xmath10 0.9 photometric redshift of overdensities in the field @xcite ; and z@xmath9 = 1.4 @xcite , that is also comparable to the z @xmath0 1.5 absorption systems in the quasar spectrum @xcite . for the h1413 + 117 lens macro - model we use the convergence and shear estimates of @xcite ( @xmath11 case ) , which are nearly identical to those of @xcite . we assume an 80% smooth matter percentage , consistent with the values estimated in other lenses where the line of sight does not pass through the galactic core @xcite . for q2237 + 0305
we use the macro - model of @xcite , and trott & wayth ( priv .
comm . ) . in this case
lensed images are clustered around the star - dominated galactic core and microlensing is consistent with a 0% smooth matter percentage @xcite , which we employ here .
we obtained simulated maps of microlensing magnifications from the gerlumph database @xcite based on these models .
we simulate the expected microlensing behavior of these lenses , using monte carlo simulations .
emission region surface brightness distributions are simulated with gaussian profiles whose half - light radii represent emission region sizes . note that except in the case of very strong microlensing events , half - light radius rather than the detailed surface brightness distribution governs microlensing properties @xcite .
a pair of random locations on gaussian - convolved magnification maps ( one map for each of a pair of lensed images ) yields a possible magnification ratio between those lensed images . by varying half - light radius and location we
can then explore microlensing behavior as a function of emission region size .
we determine d / a ratios for h1413 + 117 at 50,000 location pairs and b / a ratios for q2237 + 0305 at 25,000 location pairs .
we explore sizes from @xmath12  m to @xmath13  m , constrained by the resolution and size of the gerlumph maps .        figure  [ ratiovssize ] ( top 3 panels ) shows magnification ratio versus emission region size for a random selection of 26 of these locations .
magnification ratios are normalized to unity at the macro - model magnification ratios ( the large - scale average ratios ) .
these figures , which we will refer to as _ convergence plots _ , illustrate two important points :    \1 ) as expected , magnification ratios will converge on the macro - model ratios above a certain emission region size ( the _ convergence size _ ) .
we see that for h1413 + 117 ( top 2 panels ) the d / a convergence size is @xmath14  m ( @xmath020 light days ) for both the @xmath15 and @xmath16 cases . if the wavelength independence that we observe in the bel magnification ratios ( sect .
[ fluxratios ] ) is due to such convergence , then this gives us an approximate lower limit on the size of observed belrs . in the case of q2237 + 0305 ( panel 3 )
the convergence of the b / a magnification ratio does not occur until significantly greater than @xmath17  m , beyond the limit of our simulated magnification maps , and so the observation of differential microlensing across the full velocity width of its belr is not surprising .
\2 ) we expect tracks on the convergence plot to be flat for sizes larger than the convergence size , however we also see many flat tracks for very small emission regions .
this independence of magnification ratio with size often occurs in cases where the magnification ratio deviates significantly from the macro - model value ( i.e. where there is strong microlensing ) . for h1413 + 117 ,
most simulated tracks are relatively flat below @xmath18  m. there is an intermediate range over which all tracks converge on the macro - model ratios  the _ convergence range_. the convergence range spans @xmath18  m to @xmath14  m for h1413 + 117 , and in it we see strong size dependence of magnification ratios , and here the only flat tracks are those already close to the macro - model ratio .
the interpretation of this behavior is straightforward : the small - scale flatness seen the convergence plots occurs when a small emission region falls in relatively featureless regions between strong magnification gradients in the caustic network . as size increases into the convergence range ,
our larger source is beginning to overlap high or low magnification structures in the caustic network , resulting in size - dependent , differential microlensing . when the emission region is large enough to sample many such structures , tracks again flatten at the macro - model magnification ratio .
the magnification ratios of the civ and ciii ] bels in h1413 + 117 vary by less than 5% across the unabsorbed component of their velocity structure ( sect .
[ fluxratios ] ) , and so must reside on a `` flat '' region of a convergence plot track .
we can place more formal constraints on the possible sizes of these emission regions by determining the number of tracks that show similar flatness over a given size range .
we conservatively assume that the broad line gas flow spans a size range of 0.5 dex across its entire velocity structure .
in fact the belr probably spans a larger range , but choosing a smaller range yields more tracks within our flatness constraint , and so produces weaker ( more conservative ) size constraints .
figure  [ ratiovssize ] ( bottom ) shows the percentage of tracks consistent with @xmath195% variation in the magnification ratio over 0.5 dex in emission region size as a function of the maximum size of the emission region , for both h1413 + 117 ( @xmath15 ) and q2237 + 0305
.    in the case of q2237 + 0305 figure [ ratiovssize ] shows a reasonable likelihood of observing size - independent microlensing over this full range of emission region sizes . for h1413 + 117 ,
size - independent microlensing is unlikely ( @xmath20 5% ) for belrs spanning 0.5 dex with maximum size in the range @xmath18  m to @xmath21  m ( @xmath20 1% for @xmath22  m to @xmath23
m ) . at larger sizes ,
size independence will become more probable , and extrapolating our results we estimate that size independence will become likely ( @xmath750% ) at belr sizes of 2 - 3@xmath24  m , or around 100 light days .
it is important to note that h10 also observe this size - independent microlensing of the belr at a different epoch , greatly decreasing the probability that the belr is in the @xmath18  m to @xmath21  m range . for a reasonable probability of observing size independent magnification ratios in the belr of h1413 + 117 in two independent epochs
, the belr probably needs to be larger than @xmath17  m ( @xmath040 light days ) , and any significantly luminous component likely arises well above the convergence size of @xmath14  m ( @xmath020 light days ) .
these constraints are for the @xmath15 case ; larger values of @xmath25 push us to larger belr sizes .
we can also consider the possibility that the entire belr is emitting at sizes smaller than @xmath18  m ( @xmath04 light days ) .
we think this unlikely for two related reasons : 1 ) our observations and those of h10 show wavelength - dependent differential microlensing in the continuum d / a ratio , which precludes flat tracks at small size needed for size independence in the bel magnification ratios ; 2 ) for a high likelihood of size independence at small scales at two epochs , the emission region must be closer to a few times @xmath26  m , which is approaching the scale of the accretion disk ( see sect .
[ sizes ] ) , and yet the d / a ratios indicate very strong differential microlensing between continuum and bels at multiple epochs ( see sect . [ fluxratios ] ) .    from the reverberation mapping belr size - luminosity relation for luminous quasars of @xcite , h1413 + 117 ( with its @xmath27  f@xmath28(1350)@xmath29  erg  s@xmath30 )
is expected to have a civ belr in the range of 2590 light days . from their microlensing study of of 16 lensed quasars (
not including h1413 + 117 ) , @xcite find a belr size range for high - ionization lines in their luminous subsample of @xmath31 light days .
both of these results support our estimate of a belr that is signficantly larger than 20 light days , and very unlikely to be @xmath204 light days .
we estimate the black hole mass in h1413 + 117 at @xmath32  m@xmath33 based on the civ width and using the @xcite m@xmath34fwhm relation .
we measure the civ fwhm of 3,400  km / s based on the unabsorbed , red side .
note that the intrinsic civ bel is either symmetric or shows a blue excess , as discussed in section [ belrcomp ] , and so the correct fwhm and corresponding black hole mass may be a larger .    in order to obtain an upper limit estimate on the size of the accretion disk in h1413 + 117 we assume a shakura - sunyaev thin disk @xcite with a @xmath35  m@xmath33 non - rotating black hole accreting at the eddington limit , and observed at 150  nm .
this yields a source half - light radius of @xmath36  m , shown in figure  [ ratiovssize ] . for our black hole mass estimate of @xmath32  m@xmath33 and the same accretion parameters the 150  nm source size is @xmath37  m. based on differential microlensing across the continuum , h10 estimate a continuum source size of @xmath38  m.
the accretion disk for both quasars is much smaller than the convergence range . for h1413
+ 113 both size - dependent and size - independent magnification ratios are reasonably common at small radii , and it s not surprising that the continuum shows some differential magnification in both our observations and those of h10 .
this is also true for q2237 + 0305 , however size dependence occurs less frequently at these small size scales , in @xmath020% of cases .
the d / a magnification ratio in the h1413 + 117 continuum tends towards the average d / a bel ratio towards longer continuum wavelengths ( although never getting close to the bel ratio ) , and so it s tempting to interpret this as the accretion disk becoming large enough to start averaging over causic structures and hence approach the macro - model magnifications .
however figure [ ratiovssize ] shows that even a large accretion disk is much smaller than the convergence range . in this regime ,
magnfication ratio can have a size dependence that either approaches the macro - model value or moves away from the macro - model value , even if all tracks ultimately converge on the same value .
in short , we can not interpret a trend towards the macro - model magnifications as a reliable proxy for increasing emission region size _ except _ in the convergence range .
it is therefore not surprising that the c / a continuum ratios tend in the opposite direction to the d / a ratios , with shorter wavelengths approaching the c / a bel ratio .
proper interpretation of the differential microlensing in the continuum requires higher resolution magnification maps and a careful statistical analysis , which we leave to a later study .
we can consider the spectra of the lensed images in h1413 + 117 to be linear combinations of intrinsic spectral components from different emission regions , each subject to absorption and lensing that may vary between lensed images and emission regions .
if the relative strengths of the emission components vary between one lensed image and another due to differential microlensing or intrinsic variability , then with some assumptions it can be possible to solve for these components .
@xcite first apply this approach to h1413 + 117 to derive an intrinsic bal absorption profile for siiv , civ , and aliii .
they describe the spectrum as composed of physical components : a continuum + bel .
h10 explore the approach in more detail , applying the method of @xcite to h1413 + 117 to extract ly@xmath5+nv , siiv , civ bal troughs . in the @xmath39 method ,
the quasar spectrum is explicitly described in terms of a component that is only macrolensed and one that is both micro- and macrolensed .
both approaches lead to similar results , and this is because their assumptions are similar : that the continuum is microlensed but the belr is not ; as a result the micro+macrolensed component of h10 is essentially the absorbed , lensed continuum .    here
we derive a new spectral decomposition of h1413 + 117 .
our approach is mathematically similar to both angonin et al . and h10 , and we discuss their connection in [ h10comparison ] .
we choose to formulate the decomposition in terms of physical emission regions rather than by their degree of microlensing .
we do this because : 1 ) our observation of almost constant magnification ratios across the emission line velocity structures ( see sect .
[ sigs ] ) ensures a straightforward extraction of the bel component , and 2 ) we include in our formulation a separate bal absorption term for both the continuum and belr , which we feel adds clarity to the interpretation of the derived components , and 3 ) any source of variation  microlensing or intrinsic variability coupled with time delay  enables a decomposition , and yet the possible presence of the latter means that results can not be easily interpreted as microlensed / non - microlensed components .
we describe the spectrum * * _ f_**@xmath40 of lensed image x as a linear combination of two components : 1 ) @xmath41 : the intrinsic continuum spectrum * * _ f_**@xmath42 multiplied by the bal absorption profile affecting the continuum * * _ a_**@xmath42 and the ( potentially wavelength - dependent ) magnification * * _ m_**@xmath43 affecting the continuum in image x , and 2 ) @xmath44 : a broad emission line spectrum * * _ f_**@xmath9 with its own , potentially different absorption profile * * _ a_**@xmath9 and lensing magnification * * _ m_**@xmath45 .
the * _ f _ * and * _ a _ * components are intrinsic to the quasar and identical for all lensed images ( modulo quasar variability , but see below )
. the lines of sight to different lensed images do not diverge significantly within the bal wind ( see sect .
[ lineofsight ] ) and so all lensed images see the same absorption from this outflow
. however separate absorption terms are required becaues the accretion disk and belr are viewed through different regions of the bal wind .
the magnification factor is different for each lensed image and in the presence of microlensing this factor may differ between the continuum and bels , and each of these may have wavelength dependence .
we have neglected to include an explicit term for reddening or to account for intrinsic variability in the quasar .
both of these effects cause wavelength - dependent multiplicative changes to component spectra , and so in terms of the decomposition we can consider these to be folded into the magnification terms ; the origin of these multiplicative effects does not affect the decomposition . for lensed images 1 and 2 we have :    @xmath46    @xmath47    the only known quantities in these equations
are the observed spectra , * * _ f_**@xmath48 and * * _ f_**@xmath49 .
although we do nt know the intrinsic magnification factors of each component , it is possible to measure their ratios . in the absence of differential absorption or reddening
these are just the wavelength - dependent flux ratios for each component .
we set @xmath50 and @xmath51 = * * _ m_**@xmath52**_f_**@xmath53 : the magnification ratios in continuum and emission line .
we can then solve for @xmath54 and @xmath55 for one lensed ]    @xmath56    we can think of these equations as representing a scaled subtraction of the spectrum of one lensed image ( * * _ f_**@xmath49 ) from another ( * * _ f_**@xmath48 ) , with the scaling factor @xmath57 giving complete subtraction of one spectral component . in eq .
[ cont ] , for example , * * _ f_**@xmath49 is scaled so that its emission line flux matches that of * * _ f_**@xmath48 .
as long as the line - to - continuum ratio is different between * * _ f_**@xmath48 and * * _ f_**@xmath49 then * * _ f_**@xmath58  * * _ f_**@xmath49 eliminates the bel spectrum leaving a representation of the absorbed continuum spectrum with a slope somewhere between that of lensed images 1 and 2 .
the scaling factor @xmath59 then recovers the normalization and slope of the continuum spectrum seen in image 1 .
a similar interpretation of eq . [ line ] eliminates the continuum component to give the lensed , absorbed bel spectrum for image 1 .    for our decomposition of h1413 + 117
we define * * _ f_**@xmath48 to be the spectrum of image d and * * _ f_**@xmath49 to be that of image a. we take @xmath60 to be the wavelength - dependent continuum ratio from the fits described in section  [ sigs ] . as seen in sect.[sigs ] , the bel magnification ratios are not wavelength dependent over each line and vary relatively little between siiv , civ , and ciii ] . to reconstruct spectral components over the full wavelength range of the spectra we take @xmath61 to be a constant : @xmath62 , the mean continuum - subtracted broad line flux ratio from table  [ fluxratiotable ] . for individual bel decompositions
we use their individual magnification ratios ( see below ) .
figure  [ decomp ] shows the resulting decomposition over the full spectrum .
the upper panel shows @xmath63 : the lensed , absorbed quasar continuum ( eq .  [ cont ] ) , with the continuum magnification factor of image d. we overplot the continuum fit for the image d spectrum .
the lower panel shows @xmath64 : the lensed , absorbed bel spectrum ( eq .
[ line ] ) , again as magnified in image d. we overplot the observed spectrum of image d minus the continuum fit and also an iron emission spectrum from the templates of @xcite , scaled to roughly match the iron features in the spectrum .
figure  [ civdecomp ] shows a nearly identical decomposition as fig .
[ decomp ] over the range of the civ line , with the only difference being that we take @xmath65 , the continuum - subtracted flux ratio over the red wing of the civ line .
[ civdecomp ] also shows the range of uncertainty for these decompositions ( red lines ) : the maximum and minimum decomposed bel and continuum components obtained by allowing the continuum flux to vary by @xmath66% ( see sect .
[ sigs ] ) and @xmath67 to vary over its uncertainty range .
it can be seen that the general shape of the decomposed components does not change significantly over the range of possible decompositions .
lens image d is experiencing by far the strongest microlensing between the continuum source and belr .
however we do see some differences in broad line to continuum flux ratios in a , b and c. figure  [ comparedecomp ] shows these same decompositions over the siiv and civ lines for image pairs a & d and for b & c. in these cases we take the a / d and b / c emission line flux ratios to be the average for the siiv and civ lines for each image pair .
we see very consistent results , albeit with significantly more noise in the b & c decomposition due to the smaller differences between these spectra .        if we divide @xmath68 , the lensed , _ absorbed _ continuum spectrum , by our estimate of the lensed , unabsorbed continuum * *
_ m_**@xmath69  * * _ f_**@xmath42 we recover an estimate of the intrinsic absorption profile affecting the continuum . figure [ pcygni ] shows this for the civ , siiv , nv , and aliii absorption troughs .    .
velocity zero - point is set at line center or at lower wavelength line center for doublets .
vertical dotted lines show line centers .
letters are designations for absorption trough features discussed in section [ features ] .
[ pcygni],width=453 ]      the derived spectral components are similar to those found by h10 , and the derived absorption profile is similar to that of @xcite .
we discuss briefly similarities and differences in these works .
decompose in terms of macrolensed - only ( their @xmath70 ) and microlensed + macrolensed ( @xmath71 ) components , and use the combined spectra of lensed images a and b ( @xmath72 ) as well as that of image d ( @xmath73 ) .
the spectral subtraction scaling factors used by h10 are @xmath74 , the magnification ratio of the macrolensed - only component , and @xmath75 , the magnification ratio of the microlensed + macrolensed component .
@xmath75 is taken as the continuum ratio in regions unaffected by broad line flux , and so is similar to our @xmath60 .
they measure a single value of @xmath75 for each emission line from the average continuum bounding that line .
this accounts for changes in the continuum magnification ratio from one line to the next , but not for variation across the wavelength range of the lines , however the latter is of order 1 - 2% and so is not a large effect .
@xmath74 is determined under the assumption that at least part of each broad emission line is not microlensed , so that within the wavelength range of each line there is a region of the reconstructed micro+macrolensed component ( @xmath71 ) that is equivalent to the continuum spectrum .
this can be interpreted as subtraction of @xmath72 from @xmath73 , with @xmath72 scaled to eliminate at least part of the emission line ( the part of the line closest to the continuum magnification ) . in practice , because the emission line magnifications are constant across their wavelength range ( see sect .
[ sigs ] ) , this should mean almost complete subtraction of emission lines and produces a result that is very similar to our equation 2 .
@xcite ( a90 ) describe their decomposition in less detail than h10 , however their approach is mathematically similar to both ours and that of h10 .
they express the spectrum of lensed image d to be that of the ( purportedly ) unmicrolensed a+b+c composite plus a microlensed , absorbed continuum component . after deriving a continuum spectrum from the a+b+c composite they fit parameters to yield the absorption profile .
this appears to be equivalent to performing a scaled subtraction of a+b+c from d to eliminate the bel component and at the same time scaling out the wavelength - dependent continuum , reddening , lensing magnifications .    in summary , although a90 and h10 formulate decompositions differently and make some different assumptions , all three approaches are mathematically similar and yield similar decompositions .
it is worth noting that in our formulation the interpretation of the decomposition does not depend on their being any macrolensed - only component , only that the continuum - belr ratio differs between the two lensed images . while both h10 and a90 assume that the belr is not microlensed , any such wavelength - independent microlensing or variability - related differences in the belr do not invalidate the decompositions .
we compare the details of these decomposition results in section  [ disc ] .
the continuum spectral component ( fig .
[ decomp ] & [ civdecomp ] , upper panels ) and the absorption profile that we derive from it ( fig .
[ pcygni ] ) provide a picture of the bal absorption along the line - of - sight column to the accretion disk , and the shape of these profiles allows us to draw specific conclusions about the structure and kinematics of the outflow .
we focus on interpretations in the context of rotating disk - wind models , in which the material responsible for both bal and bel features arises from the accretion disk and is driven out by radiation and possibly gas pressure @xcite . in the case of civ
the continuum is experiencing essentially no absorption until @xmath01,500  km / s blueward of the emission line central wavelength , at which point transmission drops sharply over @xmath0600  km / s to opacity @xmath76 .
this near - black region extends to @xmath07,500  km / s before transmission ramps roughly linearly to the end of the absorption feature at @xmath015,000  km / s blueward of line center .
two narrow emission peaks are seen in this latter region . in the case of siiv
, absorption also has a sharp onset blueward of the siiv @xmath771394 line center , in this case at @xmath03000  km / s .
the deep region of its absorption trough ends much more sharply than in the case of civ at @xmath07,500  km / s .
a single clear emission peak is observed in the deep part of the trough .
a weaker absorption trough appears redward of the deeper trough onset and may result from siiv @xmath771403 . in any case , all absorption is blueward of the siiv doublet . as with civ ,
the nv trough begins at a blueshifted @xmath01,500  km / s , while the aliii trough appears at @xmath04,500 - 5,000  km / s , though may show some absorption from @xmath02,000  km / s .
of particular interest is the sharp , blueshifted onset of this absorption , which is especially striking in civ .
this sharpness is not ambiguous , and is apparent over the full range of uncertainty in the decomposition ( fig .
[ civdecomp ] , lower ) . while bals with highly blueshifted absorption onsets are sometimes observed ( so - called detached troughs , blueshifted at several thousand km / s ; e.g. turnshek 1987 ) , our decomposition of h1413 + 117
provides the first unambiguous description of a sudden , blueshifted absorption onset that is still within the velocity range of the bel .
the steepness and blueshift of this onset tells us two things : 1 ) the outflow has already accelerated to some degree before either it enters our sightline or it enters the region of efficient absorption within our sightline ; and 2 ) the kinematics of absorbing gas in our sightline is strongly dominated by outflow , not orbital motion .
@xcite intepret highly detached absorption troughs as resulting from sightlines that graze the upper layer of an equitorial outflow , hence missing the vertical base of their wind model where the flow has a small radial component .
@xcite ( h13 ) apply radiative transfer code to their disk - wind model to predict civ bal profiles , also finding that these same grazing sightlines produce a very sharp , blueshifted onset . in these equitorial disk - winds we expect any deeper line of sight
those passing near the launch region of the wind  to show absorption at low velocities . not only
do they encounter launching gas prior to significant outward acceleration , but they also encounter gas with very high rotational velocities in the vicinity of the accretion disk .
this latter effect is important , as shown by h13 ; at launch , the bal flow shares the keplerian velocity of the accretion disk and so is dominated by rotational motion . if the inner disk is occluded by this gas along our line of sight , and if this gas is significantly absorbing in the vicinity of the accretion disk , then some fraction of the continuum will be absorbed by the component of the orbiting flow that is receding along our sightline , producing partial absorption of the continuum at lower and/or redshifted projected velocities .
on the other hand , even large rotational motion far from the accretion disk will only have a transverse velocity component within our line - of - sight column to the disk .
this argument against a rotational contribution to the line - of - sight velocity is supported by @xcite , who show that , for a spherically expanding outflow , a sharp absorption onset is only expected in cases with small rotational components .
any significant rotation results in a gradual absorption onset .
this same grazing sightline in the h13 model also appears to reproduce the ramped blue edge of the civ absorption trough that we observe in h1413 + 117 .
the models of @xcite show that a ramped increase in transmission at the blue edge indicates a large variation in both velocity and density in the absorbing region along the line of sight , which appears to occur when this line of sight is close to tracing a streamline in the flow .
this is consistent with a line of sight along the upper surface of the wind .
the interpretation of @xcite contradicts this .
they classify sdss bal quasars according to the appearance of siiv and aliii troughs alongside the ubiquitous civ trough , and interpret those containing all three ( h1413 + 117 does , making it a lobal ) as resulting from lines of sight on the underside of the outflow , intersecting the base of the wind .
the argument of the grazing sightline does not explain the different velocity shifts we see in the absorption onset for different line species , which should be the same if caused by our line of sight exiting the outflow .
while ion species with the highest ionization potentials , civ ( 64.5  ev ) and nv ( 97.9  ev ) , have similar @xmath01,500  km / s absorption onsets , the onset for siiv ( 45.1  ev ) is @xmath03,000  km / s and the deeper part of the aliii ( 29.4  ev ) trough starts at @xmath04,500 - 5,000  km / s ( though may show weaker absorption from @xmath02,000  km / s ) .
this increase in absorption onset with diminishing ionization potential may suggest that for siiv and aliii we are in fact seeing very sudden jumps in transition rate as ionizing radiation and/or gas density reach optimal levels as a function of radius and hence outflow velocity .
for example , this may be the radius where it is first possible for `` locally optimized clouds '' @xcite to form due to the diminishing ionization parameter .
a combination of the above may be needed to explain the onset of all absorption troughs ; we propose that for civ and the higher ionization nv we are seeing our line of sight exit the outflow near its launch point , while for lower ionization lines we are seeing the onset of favorable absorption conditions further out in the flow . in the ionization / radiative transfer calculations of h13
, high civ fractions are expected at the highest ionization parameters in their luminous benchmark quasar ( m@xmath78  m@xmath33 , eddington ratio = 0.2 ) , and so it is expected that civ absorption should be seen at its inner - most radii .
figures [ decomp ] & [ civdecomp ] ( lower ) shows a plausible bel spectrum , however it should be remembered that this component includes the bal absorption of this broad line emission .
it is not possible to extricate the true absorption profile @xmath79 from the line profiles @xmath80 . in the case of the civ line ( fig .
[ civdecomp ] , lower ) we see a region of near - zero emission in the rest frame 1515 - 1525  range ( @xmath045006400 km / s ) where we might otherwise expect to see broad - line emission  either civ or feii , given the presence of one or both surrounding this region , and the significant iron emission in this wavelength range of the @xcite template .
we take this as evidence of absorption of broad - line emission by the bal wind .
if this black region is due to absorption then we are seeing very high column densities _ and _ very high covering factors for absorption of the civ and/or feii emission at these velocities .
this is very interesting due to the presumably much larger extent of the source of this broad - line emission compared to the continuum  the line - emitting gas is completely enveloped by absorbing gas across the entire belr at these projected velocities .
further evidence of absorption of bel light by the bal wind comes from the significant absorption of ly@xmath5 by the nv trough in this source ( see also turnshek 1987 ) .
while the symmetry of the civ emission in this component resembles that of an unabsorbed bel , it should be remembered that civ very frequently shows a strong blue excess which may be absorbed in this case .
the siiv line shows a sharp drop in flux in its blue wing that also suggests absorption of bel light .
the iron emission template of @xcite shown in these figures provides a good match to several features in the emission line spectrum , and an imperfect match to others .
the complexity and potential variation in iron emission makes a perfect match unlikely , however the scale and shape of these features match the template adequately . as well as natural variation , it is also possible that we are seeing differential microlensing acrosss the iron emission spectrum .
the strong presence of iron features in component @xmath81 implies that some of this emission arises on size scales similar to other bels .
the fact that the iron template provides a reasonable match to the iron emission across the full wavelength range of the belr component ( fig .
[ decomp ] ) gives us confidence in our choice of continuum , and that we avoided most of the iron emission ; if not then we would expect a significant and wavelength - dependent mismatch between the decomposed bel spectrum and the template .
low - velocity - width emission peaks are seen within the region of the civ bal absorption trough , marked in figure [ civdecomp ] : in the continuum component centered at blueshifts of @xmath08,000  km / s ( _ a _ ) ) and @xmath011,000  km / s ( _ b _ ) and in the bel component at @xmath04,000  km / s ( _ c _ ) , @xmath07,500  km / s ( _ d _ ) , and @xmath013,000  km / s ( _ e _ ) .
an emission peak is also observed blueward of siiv @xmath771393.76 at @xmath04,500  km / s ( which we designate feature _ f _ , fig .
[ pcygni ] ) .
h10 observe a similar emission peak in their decomposition ; it appears in all epochs in @xmath06,5009,500  km / s blueward of civ in their @xmath82 ( bel ) component ( from the decomposition of lensed image d with respect to a & b ) .
this appears to correspond to our feature _
d_. also notable is a feature in their @xmath83 ( continuum ) component from the decomposition of image c which corresponds to our feature _ a_. a similar peak was also observed in the civ absorption profile derived by a90 at @xmath05,500  km / s . in all of these previous decompositions ,
signal quality makes it is difficult to rule out the presence of any of the other features that we note in the current work .
the features observed in the civ and siiv troughs of the continuum component ( _ a _ , _ b _ & _ f _ ) indicate regions of reduced absorption in the bal outflow . this may be due to radial fluctuations in overall outflow density and/or relative ion species density due to variations in ionization parameter .
this interpretation reflects that of h10 , who propose a strongly absorbing cloud embedded within the outflow but with lower velocity width that results in a double - peaked civ line . in the context of a sightline grazing the upper layer of a fast , equitorial outflow , @xcite
show that , in the hydrodynamic model of @xcite , `` blobs '' of failed outflow can cause clumpiness in this layer as they fall from the over - ionized high - inclination regions .
interpretation of the features in the bel component is more complicated .
d _ and _ e _ may result from clumpy absorption of a strong blue excess in civ and/or feii emission ( which may or may not be absorbed ) . as discussed in section [ belrcomp ] , the near - black region between features _
c _ and _ d _ does seem to indicate absorption of bel light by the outflow . whether this is also true of the dark region between features _
d _ and _ e _ is unclear , as feature _ e _ may be isolated feii emission .
c _ is probably iron emission .
it is clear that features _
a _ and _ d _ correspond to the same peak in the original spectrum , and so to the same region of diminished absorption .
b _ in the continuum component has no obvious counterpart in the bel component .
this may be due to intrinsically low bel flux in this region or a difference in the bel vs. continuum absorption profiles .
the differences observed in the bal profiles between lens images can not be due to differences in the line of sight .
the bal absorption profile affecting the continuum is determined by the region of the bal flow sampled by the column along the line of sight to the accretion disk , and so different lensed images will only experience different absorption if their line - of - sight columns diverge significantly within this outflow . the image separation between images a and d is 118 , which corresponds to a divergence angle of the a and d lines of sight from the source of 066 and 20 for @xmath15 and 1.4 respectively . assuming @xmath16 and a conservatively small accretion disk with half - light radius @xmath26 m ( ss disk with @xmath84 and @xmath85 ) , lines of sight have diverged by one half - light radius at @xmath0100 light years .
this is probably much larger than the strongly - absorbing region of the bal flow .
choosing a smaller @xmath25 gives an even smaller divergence . on the @xmath201000 light day size scale of a large belr @xcite , line - of - sight columns to the accretion disk
still have @xmath86% overlap .
line - of - sight columns to the belr , assuming 100 light day cross - section , diverge completely at @xmath028,000 ly , and so may be subject to different line - of - sight absorption features both within the host galaxy and in the igm .
this is also true of the continuum , however the continuum will be subject to differences in absorption features from a greater depth within the host galaxy .
these line - of - sight differences do account for the many differences in narrow absorption features observed in our spectra .
@xcite show that these features arise both in the igm and potentially the lensing galaxy , but probably not the host galaxy which has a higher redshift than observed features .
we do nt investigate the absorption features in detail in this paper , however we note that the @xmath87 feii lines ( @xmath88 2344.2 , 2374.5 , 2382.8 ) follow the same behavior in our spectra as observed by @xcite ; it is present in lensed images a and b and absent in c and d. several other features exhibit similar behavior , and we will investigate these in a separate study .
we have presented new spectroscopy of the gravitationally lensed broad absorption line quasar h1413 + 117 using gemini north s gmos ifu .
extracted spectra of the four lensed images reveal that differential microlensing strongly affects the relative magnification of the broad emission lines compared to the continuum in lensed image d , with the continuum in d appearing anomalously bright compared to all other lensed images . lensed images a , b , and c also show small differences in their continuum - to - broad line flux ratios relative to each other .
subtracting a fit to the continuum emission outside the broad emission and absorption regions allows us to determine the wavelength - dependent d / a flux ratios of the broad emission lines .
these flux ratios are flat , indicating no differential microlensing across the velocity structure of the belr .
monte carlo simulations were performed to analyze the dependence of the flatness of broad line flux ratios on the the size range spanned by the belr velocity structure .
it was found that belrs with an upper size limit in the @xmath89@xmath17  m range , and spanning 0.5 dex across their velocity structure are very unlikely to produce the observed flat flux ratios .
smaller belrs are also unlikely due to the strong differential magnification observed between broad lines and continuum .
therefore , it is inferred that the belr in h1413 + 117 is significantly larger than @xmath17  m , and produces no significant emission from regions smaller than @xmath14  m , or 20 light days .
we take advantage of the differential magnification between lensed images a and d to algebraically decompose the quasar spectrum into two components : the absorbed broad line spectrum and the absorbed continuum spectrum , and we use the latter to derive intrinsic absorption profiles .
our decomposition leads us to these conclusions :    @xmath90 a sharp , blue - shifted onset of absorption is observed in all prominent bal features in the continuum component . in the case of civ ,
nv , and ly@xmath5 this onset is at @xmath01,500  km / s , which may represent the inner edge of a radial outflow . in the context of typical disk - wind models this would place the onset of bal absorption near its launch point , and so in the vicinity of the accretion disk . in the case of siiv and aliii , which have a larger absorption onsets , we may be seeing increases in absorption efficiency of lower ionization species with distance from the ionizing source .
this sharp onset also indicates that the line - of - sight velocity is dominated by outflow , not rotational motion .
@xmath90 overall the derived absorption profile for civ resembles that of an accelerating outflow that is either spherical or viewed relatively close to a streamline , albeit with a @xmath01,500km / s initial velocity .
@xmath90 the civ and/or feii belr appears to be strongly absorbed in the range @xmath04,5006,400  km / s .
it is unclear whether emission features blueward of this trough are civ or feii , however the complete absorption of either of these features suggests a very high covering factor for the absorbing material at large radii .
combined with the apparent strong absorption near the outflow launch point , this suggests that gas responsible for bal features maintains a high - covering - factor over at least two orders of magnitude in size .
angonin m .- c .
, vanderriest c. , remy m. , surdej j. , 1990 , a&a , 233 , l5 allington - smith , j. , graham , m. , content , r. , dodsworth , g. , davies , r. , miller , b. w. , jorgensen , i. , hook , i. , crampton , d. , murowinski , r. , 2002 , pasp , 114 , 892 anguita , t. , faure , c. , yonehara , a. , wambsganss , j. , kneib , j .-
p . , covone , g. , alloin , d. , 2008 , a&a , 481 , 615 baldwin , j. , ferland , g. , korista k. , verner , d. , 1995 , apj , 455 , l119 bate , n. f. , floyd , d. j. e. , webster , r. l. , wyithe , j. s. b. , 2008 , mnras , 391 , 1955 bate , n. f. , floyd , d. j. e. , webster , r. l. , wyithe , j. s. b. , 2011 , apj , 731 , 71 bentz , m. c. , walsh , j. l. , barth , a. j. , et al . , 2009 , apj , 705 , 199 blackburne , j. a. , pooley , d. , rappaport , s. , schechter , p. l. , 2011 , apj , 729 , 34 blackburne , j. a. , kochanek , c. s. , chen , b. , dai , x. , chartas , g. , 2014 , apj , 798 , 95 blandford , r. d. , mckee , c. f. , 1982 , apj , 255 , 419b borguet , b. , hutsemkers , d. , 2010 , a&a , 515 , 22 chartas , g. , kochanek , c. s. , dai , x. , poindexter , s. , garmire , g. , 2009 , apj , 693 , 174 chelouche , d. , netzer , h. , 2005 , apj , 625 , 95c crenshaw , d. m. , kraemer , s. b. , 2007 , apj , 659 , 250 dai , x. , kochanek , c. s. , chartas , g. , kozlowski , s. , morgan , c. w. , garmire , g. , agol , e. , 2010 , apj , 709 , 278 denney , k. d. , peterson , b. m. , pogge , r. w. , et al . , 2009 , apj , 704 , l80 doroshenko , v. t. , sergeev , s. g. , klimanov , s. a. , pronik , v. i. , efimov , y. s. , 2012 , mnras , 426 , 416 eigenbrod , a. , courbin , f. , sluse , d. , meylan , g. , agol , e. , 2008 a&a , 480 , 647 elvis , m. , 2000 , apj , 545 , 63 fabian , a. , 1999 , mnras , 308 , l39 fabjan , d. , borgani , s. , tornatore , l. , saro , a. , murante , g. , dolag , k. , 2010 , mnras , 401 , 1670 faure , c. , alloin , d. , kneib , j-.p .
, courbin , f. , 2004 , a&a 428 , 741 filiz ak , n. , brandt , w. n. , hall , p. b. , et al . , 2014 , apj , 791 , 88 floyd , d. j. e. , bate , n. f. , webster , r. l. , 2009 , mnras , 398 , 233 gallagher , s. c. , everett , j. e. , 2007 , _ the central engine of active galactic nuclei _ , asp conference series , vol .
l. c. ho and j .- m .
wang gallagher , s. c. , hines , d. c. , blaylock , myra , priddey , r. s. , brandt , w. n. , egami , e. e. , 2007 , apj , 556 , 619 giustini , m. , proga , d. , 2012 , apj , 758 , 70 goicoechea ,
l. j. , shalyapin , v. n. 2010 , apj , 708 , 1001 guerras , e. , mediavilla , e. , jimenez - vicente , j. , kochanek , c. s. , muoz , j. a. , falco , e. , motta , v. , 2013 , apj , 764 , 160 hall , p. b. , brandt , w. n. , petitjean , p. , et al . , 2013 , mnras , 434 , 222 higginbottom , n. , knigge , c. , long , k. s. , sim , s. a. , matthews , j. h. , 2013 , mnras , 436 , 1390 hook , i. , jorgensen , i. , allington - smith , j. r. , davies , r. l. , metcalfe , n. , murowinski , r. g. , crampton , d. , 2004 , pasp , 116 , 425      jimnez - vicente , j. , mediavilla , e. , kochanek , c. s. , muoz , j. a. , motta , v. , falco , e. , mosquera , a. m. , 2014 , apj , 783 , 47 kaspi , s. , smith , p. s. ; netzer , h. ; maoz , d. , jannuzi , b. t. , giveon , u. , 2000 , apj , 533 , 631 kaspi , s. , brandt , w. n. , maoz , d. , netzer , h. , schneider , d. p. , shemmer , o. , 2007 , apj , 659 , 997 keeton , c. r. , burles , s. , schechter , p. l. , wambsganss , j. , 2006 , apj , 639 , 1 kneib , j .-
p . , alloin , d. , mellier , y. , et al . , 1998 ,
a&a , 329 , 827 kochanek , c. s. , 2004 , apj , 605 , 58 lamy , h. , hutsemkers , d. , 2004 , a&a , 427 , 107 magain , p. , surdej , j. , swings , j .-
p . , borgeest , u. , kayser , r. 1988 , nature , 334 , 325 mediavilla , e. , muoz , j. a. , falco , e. , motta , v. , guerras , e. , canovas , h. , jean , c. , oscoz , a. , mosquera , a. m. , 2009 , apj , 706 , 1451 monier , e. m. , turnshek , d. a. , lupie , o. l. , 1999 , apj , 496 , 177 mortonson , j. m. , schecter , p. l. , & wambsganss , j. , 2005 , apj , 628 , 594 mosquera , a. m. & kochanek , c. s. , 2011 , apj , 738 , 96 murray , n. , chiang , j. , 1995 , apj , 454 , l105 murray , n. , chiang , j. , grossman , s. a. , voit , g. m. , 1995 , apj , 451 , 498 monier , e. m. , turnshek , d. a. , lupie , o. l. 1998 , apj , 496 , 177 odowd , m. , bate , n. f. , webster , r. l. , wayth , r. , labrie , k. , 2011 , mnras , 415 , 1985 ostensen , r. , remy , m. , lindblad , p. o. , refsdal , s. , et al .
, 1997 , a&as , 126 , 393 peterson , b. m. , 1993 , pasp , 105 , 247 pooley , d. , rappaport , s. , blackburne , j. , schechter , p. l. , schwab , j. , wambsganss , j. , 2009 , apj , 697 , 1892 proga , d , stone , j. m. , kallman , t. r. , 2000 , apj , 543 , 686 proga , d , kallman , t. r. , 2004 , apj , 616 , 688 rafferty , d.a . , mcnamara , b.r . , & nulsen , p.e.j . , 2008 ,
apj , 687 , 899 reichard , t. a. , richards , g. t. , hall , r. b. , et al .
, 2003 , aj , 126 , 2594 richards , g. t. , keeton , c. r. , pindor , b. , et al . , 2004 , apj , 610 , 679 scannapieco , e. , oh , s. p. , 2004 ,
apj , 608 , 62 shakura , n. i. , sunyaev , r. a. , 1973 , a&a , 24 , 337 silk , j. , rees , m. , 1993 , a&a , 331 , l1 sim , s.a .
, long , k. s. , miller , l. , turner , t. j. , 2008 , mnras , 388 , 611 sluse , d. , claeskens , j .- f . , hutsemkers , d. , surdej , j. , 2007 , a&a , 468 , 885 sluse , d. , schmidt , r. , courbin , f. , hutsemkers , d. , meylan , g. , eigenbrod , a. , anguita , t. , agol , e. , wambsganss , j. , 2011 , a&a , 528 . 100 storchi - bergmann , t. , 2010 , _ proceedings of the iau symposium no .
267 _ , eds .
peterson , r.s .
somerville and t. storchi - bergmann somerville , r. s. , hopkins , p. f. , cox , t. j. , robertson , b. e. , hernquist , l. , 2008 , mnras , 391 , 481 surdej , j. , fraipont - caro , d. , gosset , d. , refsdal , s. , remy , m. , 1993 , proc .
31st lige international colloq .
( ligel univ .
lige ) , eds .
j. surdej , d. fraipont - caro , e. gosset , s. refsdal , m. remy trott c. m. , treu t. , koopmans l. v. e. , webster r. l. , 2010 , mnras , 401 , 1540 turnshek , d. a. , 1987 , in _ qso absorption lines : probing the universe _ , ed .
j. c. blades , d. turnshek , & c. norman , ( cambridge : cambridge university press ) , 17 van dokkum , p. g. , 2001 , pasp 113 , 1420 vernardos , g. , fluke , c. j. , bate , n. , f. , croton , d. , 2014 , apj , 211 , 16 vestergaard , m. & wilkes , b. j. , 2001 , apjss , 134 , 1 vestergaard , m. , peterson , b. m. , 2006 , apj , 641 , 689 wang , j. , 2015 , newa , 37 , 15 wayth , r. b. , odowd , m. j. , webster , r. l. , 2005 , mnras , 359 , 561 weymann , r. j. , morris , s. l. , foltz , c. b. , hewett , p. c. , 1991 , apj , 373 , 23 witt , h. , j. , mao , s. , schechter , p. , l. , 1995 , apj , 443 , 18 zubovas , k. , king , a. r. , 2012 , mnras , 426 , 2751