most active galactic nuclei ( agns ) show prominent fe ii  emission lines in their spectra .
the lines appear in several broad bands that represent thousands of individual transitions .
the stronger bands cover the ranges of 40005400   ( hereafter optical fe ii  lines ) , 28003500   and 20002600   ( hereafter uv fe ii  lines ; e.g. , @xcite and references therein ) .
the strongest fe ii  lines ( relative to h@xmath0 ) are observed in narrow - line seyfert 1 galaxies ( nls1s , see , e.g. , @xcite ) .
such objects are characterized by : 1 ) narrow ( @xmath1 2000 @xmath2 ) broad emission lines ; 2 ) weak [ o iii ]  line ; 3 ) and steep hard x - ray spectrum . these objects are found usually at the extreme end of the so - called eigenvector 1 sequence @xcite indicating high eddington ratios and several other properties that are not fully understood .    some nls1s show evidence for super - eddington accretion , with @xmath3  @xmath4 1 .
we refer to these objects as super - eddington accreting massive black holes ( seambhs ) .
our earlier work @xcite show that such objects are potentially a new kind of standard candle for cosmology . to test this idea , and to study in more detail the physical properties of extreme nls1s
, we have initiated a large reverberation mapping ( rm ) campaign to measure , accurately , the black hole ( bh ) mass of such sources .
the initial results of this campaign have been published in ( * ? ? ?
* hereafter paper i ) and ( * ? ? ?
* paper ii ) .
the present paper , the third in this series , is dedicated to the study of fe ii  emission lines in our sample of seambh candidates .
a fourth paper , submitted in parallel to this one , presents new data on h@xmath0  time lags and bh mass in five extreme seambhs and focuses on the modification of the @xmath5@xmath6  relationship in agns ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) in the presence of such sources .
the excitation mechanism of the fe ii  emission in agns has been discussed in numerous publications ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
most of these studies suggest an origin in the broad - line region ( blr ) but the line intensities calculated so far are in poor agreement with most observations .
other suggestions connect these lines to the outer part of the central accretion disk ( e.g. , * ? ? ?
* ) but the agreement with observations is still poor . in particular , none of the existing models can explain , satisfactorily , the relative intensity of the uv and optical fe ii  lines and the observed fe ii  spectrum in the range 20002600  .
line variability is an important tool in such studies since it can indicate , given a measured time lag relative to the continuum variations , the location of the fe ii - emitting gas .
such variations have been detected in a number of sources ( e.g. , @xcite ) , but robust lag measurements were not obtained .
this situation has changed , recently .
@xcite revisited the data of pg 1700@xmath7518 from @xcite and measured the light curve of the optical fe ii  lines .
they obtained a significant time lag relative to the 5100   continuum albeit with a very large uncertainty . @xcite and
@xcite adopted the multivariate correlation function ( mcf ) scheme of @xcite designed for photometric rm .
@xcite studied the nls1 sdssj113913.91 + 335551.1 using rm and found the fe ii  time lag is consistent with the one of the h@xmath0  line . @xcite
measured optical fe ii  light curves , and time lags , for three of the objects in @xcite , one of which is the object studies by @xcite .
they also find somewhat less significant fe ii  time lag for three other sources .
@xcite presents a tentative fe ii  size - luminosity relation and suggest that the fe ii  emission - region size is comparable to that of the h@xmath0line .
@xcite used spectroscopic measurements of optical fe ii  lines and found lags that are 1.5 and 1.9 times longer than the corresponding h@xmath0 .
all these studies confirm the photoionization origin of the fe ii  lines .
the success of the @xcite campaign is due both to the detailed and frequent spectroscopic observations and a novel method they used to correct for the host galaxy contribution to the optical spectra .
this paper presents the results obtained , so far , to measure the time lags of the optical fe ii  lines ( hereafter fe ii  lines ) in our seambh campaign .
we have attempted to detect such lags in 10 of the sources and were able to obtain statistically significant results in nine of them . in six of the newly measured sources , the fe ii  lag is entirely consistent with the h@xmath0  lag and in two others it is considerably longer .
section [ sec - obs ] gives a brief review of the observations and data reduction .
section [ sec - method ] describes our spectral fitting method , with emphasis on host galaxy and narrow line subtraction which we find to be crucial to the analysis .
more details are given in appendix [ sec - host ] .
section [ sec - result ] presents fe iilight curves and their analysis for both fe ii  and h@xmath0 , and compares the results of the new h@xmath0  lags to those presented in paper i and paper ii . in section [ sec - discussion ] , we plot the size - luminosity relation for fe iiand compare the lag and intensity of fe ii  with h@xmath0 . the implications and some additional interpretation
are also discussed .
section [ sec - summary ] gives a summary of the new results .
the details of the seambh campaign , including the observations , data reduction and analysis , were presented in paper i and paper ii . for completeness , we summarize the more important points below and discuss in detail the new method of galaxy and narrow line subtraction .      ten nls1s identified as seambh candidates were observed , spectroscopically and photometrically , between october 2012 and june 2013 .
objects names and coordinates are listed in table 1 of paper ii .
h@xmath0  time lags for three of the sources
( mrk 335 , mrk 142 , iras f12397@xmath73333 ) are presented in paper i and for five additional sources ( mrk 1044 , mrk 382 , mcg @xmath706@xmath826@xmath8012 , mrk 486 , mrk 493 ) in paper ii .
like many other nls1s ( * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* e.g. ) , all the sources in our sample show strong fe ii  emission lines and high @xmath3 .
this sample is , therefore , different from most agns in the local universe and none of the results presented below should be compared with earlier studies , like @xcite that address the general population properties .
the spectra were obtained using the yunnan faint object spectrograph and camera ( yfosc ) , mounted on the lijiang 2.4 m telescope at yunnan observatory of the chinese academy of sciences .
a longslit with projected width of 2@xmath95 was oriented to take the spectra of the object and a nearby non - varying comparison star simultaneously following , e.g. , @xcite and @xcite .
the comparison star is then used as a standard for flux calibration . for the lijiang 2.4 m telescope ,
the rotator is accurate and the tracking is stable .
thus the object and the comparison star were kept within the slit during the typical 30 min .
the distance between the object and the comparison star along the direction of slit width keeps less than 1 pixel ( 0@xmath9283 ) .
grism 14 was used and yielded spectra covering the wavelength range of 38007200   with a dispersion of 1.8   pixel@xmath10 .
the final spectral resolution , obtained by comparing the width of the [ o iii ]  emission line with the one measured from the sloan digital sky survey ( sdss ; @xcite ) spectrum of the same object , is roughly 500 @xmath2 .
all spectra are extracted in a uniform , large aperture of 8@xmath95 to minimize light losses .
the flux calibration is based on the comparison of the object flux to the comparison star flux in the 2@xmath95@xmath118@xmath95 apertures
. however , the host galaxy flux calibration is different because of the non - stellar image  the galaxy is extended and resolved .
this can result in apparent flux variations due to variable seeing and mis - centering .
this requires a different calibration procedure that has a large impact on the various light curves .
appendix [ sec - host ] gives details of this procedure .
the traditional method to measure the flux of the continuum and the broad emission lines , in most rm studies , is simple integration ( e.g. , * ? ? ?
a straight line is set between two line - free windows to define the agn continuum , and the flux of the emission line is measured by simple integration above the line .
this method works well for single , strong emission lines , e.g. , h@xmath0 , because there are no ( or only weak ) other emission lines in this range and the wavelength window is narrow enough for the continuum to be approximated by a straight line . in the case of fe ii , neither condition is satisfied .
fe ii  emission consists of thousands of lines that form a pseudo - continuum .
the most prominent features in the optical band are two bumps between 4500 to 5500  : one between h@xmath0  and h@xmath12 , and the other on the red side of [ o iii]@xmath135007 . even for quasars ,
whose host galaxies are extremely faint relative to the agn , it is hard to find `` pure '' continuum windows , and it is inappropriate to assume a straight line over such a wide wavelength range of more than 1000  . the strong host galaxy contribution in most of our objects , and the contamination by other emission lines ( e.g. , he ii  and coronal lines ) , make the situation much worse . appendix [ sec - feintg ] presents the fe ii  light curves measured by the traditional integration method .
all light curves have large scatter , and only less than a half show rough structures .
a simultaneous fitting including the fe ii  emission and all the other spectral components is necessary .
template - fitting is a widely used method to measure fe ii  emission in single - epoch spectra of agns ( see , e.g. , @xcite for a brief overview ) .
input - output simulations show that template - fitting is a reliable measurement of fe ii  emission with equivalent width ( ew ) @xmath4 25   in quasars @xcite .
it is also common to include a host galaxy component into the fitting of type i agns with strong host contribution ( e.g. , * ? ? ?
* ; * ? ? ?
the method has been adopted , recently , to measure the light curves in a few reverberation mapping studies , and proved to be successful .
@xcite reanalyzed the spectra of pg 1700@xmath7518 , using data from @xcite .
they detected a lag of @xmath14 days for fe ii  emission .
@xcite performed a very careful spectral decomposition including a power - law continuum , a host galaxy component , fe iitemplate , and other emission lines ( h@xmath0 , [ o iii ] , he ii , and he i ) . using this technique
they were able to obtain statistically significant reverberation lags of the fe ii  lines in two seyfert 1 galaxies , ngc 4593 and mrk 1511 .
the spectral fitting in the present paper follows the algorithm of @xcite , where several spectral components are fitted simultaneously by minimizing the @xmath15  via levenberg - marquardt method .
all the light curves are obtained directly from the results of the fitting . before fitting ,
we correct the calibrated spectra for galactic extinction and redshift .
we use the @xmath16-dependent galactic extinction law given by @xcite and @xcite and @xmath16 is assumed to be 3.1 . the redshift @xmath17 and @xmath18-band extinction
are taken from the nasa / ipac extragalactic database and @xcite , as listed in papers i and ii .
the left panel of figure [ fig - spec382 ] shows the fitting to the spectrum of mrk 382 taken at julian date ( jd ) 2456298 , after galactic extinction and redshift correction .
the fitting is performed in the rest frame wavelength range 41506280  , except two narrow windows around h@xmath12  and he i@xmath135876 .
the two lines are not blended with the major part of fe iiemission , and their study is beyond the scope of the present paper . we keep them out of the fitting to avoid introducing too many unnecessary parameters .
the observed spectra in the fitting windows are plotted in green , while those left out of the fitting are in black .
the fitting includes the following components : ( 1 ) a single power law , ( 2 ) fe ii  emission , ( 3 ) host galaxy .
these three components , which when combined together forms the pseudo - continuum , are plotted in blue .
( 4 ) h@xmath0  emission line plotted in magenta . ( 5 )
broad he ii  @xmath134686 emission line plotted in cyan .
( 6 ) narrow emission lines plotted in orange , including [ o iii ]  @xmath194959 , 5007 , he ii@xmath134686 , he i  @xmath134471 , and several coronal lines .
the summed model is plotted in red .
the bottom panel shows the residual spectrum .
note that , although he i  @xmath135876 is not in the fitting window , its profile is well recovered after removing the na i  @xmath195890 , 5896 ( na
d ) absorption lines from the host galaxy .    limited by the signal - to - noise ratio ( s / n ) for the individual - night spectra , there are degeneracies between several pairs of spectral components , including : ( a ) the agn power - law continuum and the host galaxy ; ( b ) fe iiemission and the broad he ii  line ; ( c ) fe ii  emission and the coronal lines .
thus , we first fit the high - s / n _ mean spectrum _ with all parameters set free .
we then fit each individual - night spectrum with the values of some parameters fixed to those obtained from the fitting of the mean spectrum .
the details of each spectral component , and the fitting parameters , are described in the following subsections .
a single power law is used to describe the featureless agn continuum .
it has two parameters : the flux density at 5100   ( @xmath20 ) , and the spectral index ( @xmath21 , defined as @xmath22 ) .
there is some degeneracy between the power - law continuum and the host galaxy .
a larger @xmath21 or a higher - flux galaxy ( @xmath23 ) both make the total spectrum redder ( note that the unphysical change in the color of the observed spectrum caused by weather or differential atmospheric refraction is avoided in our observation , as our flux calibration by the comparison star provides differential spectrophotometry in each wavelength bin of the spectrum ; see section [ sec - reduction ] ) .
for our observations , we find that the _ relative flux _ of the galaxy component is the main reason for the change in the total - spectrum slope .
figure [ fig - spec382 ] illustrates this point .
the right panel shows the spectrum of mrk 382 taken at jd 2456299 , just one day after the date of the spectrum shown in the left panel .
apparently , the flux at 5100   is @xmath2430% higher than the day before , and the color is redder .
however , the difference between the @xmath18-band magnitudes of the two nights is only @xmath240.01 mag .
we have therefore adopted an approach based on the assumption that the absolute flux of the host must be constant but the relative flux inside the slit can vary depending on the observing conditions .
this is illustrated in detail in appendix [ sec - host ] .
the conclusion about the agn continuum is in line with several , but not all studies discussing the relation between @xmath21 and luminosity in agns .
most recently , @xcite investigated this issue using the reverberation mapping data of the 17 quasars from @xcite .
they found no strong dependence of @xmath21 on the variability of the luminosity .
thus , in the fitting of each individual - night spectrum , we fix the value of @xmath21 to agree with the best fit to the mean spectrum , and let @xmath20  and @xmath23  free . comparing the left and right panels of figure [ fig - spec382 ] , only the galaxy component changes . in this way ,
the resultant @xmath20  light curve matches the @xmath18-band light curve well ( see appendix [ sec - host ] ) .
our fitting is consistent with no changes in spectral index as a function of luminosity .
there are several optical fe ii  templates available in the literature . among them , two templates , from @xcite and @xcite , are most widely used , both constructed from the spectrum of the nls1 galaxy i zw 1 .
@xcite find that the template from @xcite yields inconsistent he i  emission lines ; significant broad he i  @xmath134922 and @xmath135016 emission lines are needed while he i  @xmath134471 is constrained to have zero flux by the fitting .
we compared these two templates and found that the one from @xcite gives better fitting as judged by the smaller reduced @xmath15 .
no broad he
i  @xmath134922 and @xmath135016 emission lines are needed , which is consistent with the zero flux of he i@xmath134471 .
introducing the two broad he i  lines will add more parameters to the fit and will introduce additional degeneracy between he i@xmath134922 and h@xmath0 .
given this , we chose the fe ii  template from @xcite in this paper .
the fe ii  template is convolved with a gaussian function to be scaled , broadened , and shifted ( see * ? ?
* for details ) .
three parameters , the flux ( @xmath25 ) , the full width at half - maximum ( @xmath26 ) , and shift ( @xmath27 ) , are used for fe ii  emission in the fitting .
we let all the three parameters to be free in the fitting of each individual - night spectrum .
note that , @xmath25  is defined as the flux of the integrated fe ii  emission between 4434 and 4684   from the best - fit fe ii  model , following @xcite .
the spectral resolution of our spectra is rather low ( @xmath24500 @xmath2 ) , and in addition our objects are selected to have narrow h@xmath0  emission lines .
thus , it is hard to decompose the h@xmath0  emission line and remove the narrow component that comes from the narrow - line region .
considering that the contribution of the narrow component is weak and supposed to be constant over reverberation timescales @xcite , we treat the entire h@xmath0  emission line as one component .
the entire h@xmath0  profile is modeled by a gauss - hermite function @xcite .
all the five parameters of the gauss - hermite function are set free in the fitting of each individual - night spectrum .
the flux of the line ( @xmath28 ) is calculated from the best - fit model .      besides the strong [ o iii ]  @xmath194959 , 5007
, there are many other narrow emission lines in the wavelength range of our fit spectrum @xcite .
we identify narrow he ii  @xmath134686 , he i@xmath134471 , and several high - ionization forbidden coronal lines . in some objects ,
these narrow lines are so strong that their contamination is non - negligible to the fe ii  measurement .
figure [ fig - spec335 ] shows an individual - night spectrum of mrk 335 and our fit .
the strong coronal lines include in this case [ fe vii ]  @xmath135158 , [ fe vi ]  @xmath135176 , [ n i]@xmath135199 , [ ca v ]  @xmath135309 , [ fe vii ]  @xmath135721 , and [ fe vii]@xmath136086 .
adding these lines significantly improve the fit .
the intensities of the narrow lines differ from one object to the next . for each object , we first identify narrow lines in the mean spectrum by testing the goodness of fit . only lines that are identified in the mean spectrum are included in the fitting of individual - night spectra .
each narrow emission line is modeled by a gaussian .
these lines are coming from the narrow - line region and hence do not vary on the campaign time scale .
thus , we constrain them to have the same velocity width and shift with those of the [ o iii ]  @xmath135007 line .
the flux ratio relative to [ o iii ]  @xmath135007 is also kept constant , as given by the best fit to the mean spectrum ( the flux ratio of [ o iii]@xmath134959 to @xmath135007 is fixed to the theoretical value of 1/3 ) , in order to avoid the degeneracy between these lines and fe ii .
only the flux , velocity width , and shift of [ o iii ]  @xmath135007 are set free .
the broad he ii  @xmath134686 emission line is strong in the spectrum of some objects in our sample .
figure [ fig - spec335 ] shows an example .
the line , plotted in cyan , is much broader than h@xmath0 , and strongly blueshifted .
such a line profile is common for high - ionization lines in nls1s and consistent with that of the ultraviolet ( uv ) he ii  @xmath131640 line in sdss spectra @xcite .
it is even more prominent in the root - mean - square ( rms ) spectra of our objects , indicating large variability similar to the results in previous he ii  reverberation mapping studies ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the broad he ii  line is heavily blended with fe ii  emission .
the high - ionization lines in agn are often asymmetric @xcite , but there is no other high - ionization line in the wavelength range of our spectra to constrain the profile of he ii  line .
we first fit the line profile in the mean spectrum using a single gaussian .
we then kept the width and the shift in the individual - night spectra but leave the line intensity as a free parameter .
the intensity variation of the he ii  lines in our sample will be discussed in a forthcoming publication .
here we only note that the procedure described above seems to be consistent with the observations , and , in general , the lag of this line is much shorter than the h@xmath0  lag .
following @xcite , we use single simple stellar population models from @xcite as templates for the galaxy component . for most objects in our sample , the instantaneous - burst model with an age of 11 gyr and solar metallicity ( @xmath29 ) provided sufficiently good fit to the mean spectrum and consistent flux ratio compared with the measured _ hubble space telescope _ ( _ hst _ ) images ( see papers i and ii ) . in a few cases ( mrk 335 , mrk 142 , and mrk 42 )
, this template gives a flux that is much larger than that derived from the _ hst _ image .
the template with 11 gyr and @xmath30 provides better results and was adopted in these cases .    like the fe ii  template
, the galaxy template is also convolved with a gaussian to be scaled , broadened , and shifted .
as described previously , and detailed in appendix [ sec - host ] , the relative flux @xmath23  is free to vary in the fitting of individual - night spectra .
table [ tab - fitfix ] lists the values of the parameters fixed in the fitting of individual - night spectra for all the objects in our sample .
columns ( 5 ) to ( 12 ) list the narrow emission lines included in the fits .
as explained , their velocity widths and shifts are constrained to be the same with those of the [ o iii ]  @xmath135007 line while the relative intensity ratios are kept as measured from the mean spectrum .
the last column gives the galaxy template used . in the table
, a blank entry means that spectral component is not included for the object .
lccccccccccccc mrk 335 & @xmath81.45 & 7082 & @xmath8620 & & 0.042 & 0.115 & 0.043 & 0.032 & & 0.028 & 0.057 & 0.079 & 11gyr_z05 + mrk 1044 & @xmath82.01 & 5161 & @xmath8942 & & 0.216 & & 0.117 & 0.219 & & & 0.205 & 0.225 & 11gyr_z02 + iras 04416@xmath71215 & free & & & & & & & 0.095 & 0.042 & & & & + mrk 382 & @xmath81.96 & 5797 & @xmath81322 & & 0.028 & 0.074 & & & & 0.021 & 0.033 & 0.021 & 11gyr_z02 + mrk 142 & @xmath82.11 & 4824 & @xmath8757 & & 0.041 & 0.131 & 0.055 & 0.116 & & & 0.082 & 0.095 & 11gyr_z05 + mcg @xmath706@xmath826@xmath8012 & @xmath80.84 & 5191 & @xmath8927 & & 0.063 & 0.096 & & & & & 0.046 & 0.061 & 11gyr_z02 + iras f12397@xmath73333 & @xmath82.28 & 5506 & @xmath8334 & & & 0.034 & 0.013 & & 0.018 & & 0.008 & 0.010 & 11gyr_z02 + mrk 42 & @xmath80.70 & 1472 & @xmath855 & & 0.069 & 0.093 & & & & & & & 11gyr_z05 + mrk 486 & @xmath80.76 & 4877 & @xmath8724 & & & 0.171 & 0.027 & & & 0.067 & 0.050 & 0.140 & 11gyr_z02 + mrk 493 & @xmath80.91 & 4054 & @xmath81580 & & 0.132 & 0.261 & & & & & & & 11gyr_z02    figure [ fig - spec ] shows examples of fittings to individual - night spectra for the eight objects not presented so far .
the notations and colors are as the same as those in figure [ fig - spec382 ] for mrk 382 .
notes on two objects with unusual treatments are as follows .
+     +    * iras 04416@xmath71215 . * this object has the highest redshift ( 0.089 ) in our sample . from the images taken for slit centering , the point spread function of the object is as broad as that of the comparison star , meaning that there is a negligible galaxy contribution . for this object ,
the comparison star is fainter than the agn , especially in the red part .
we found that the shape of the spectrum of this object is not as well calibrated by the comparison star as in the other objects .
thus , we fit in a relatively narrow wavelength window ( 44305550  ) , but covering the two bumps of fe ii  emission on either side of h@xmath0 .
we also let the spectral index @xmath21 of the agn power - law continuum free to vary in individual - night spectra to compensate for the difference in spectral shape caused by the calibration .
as shown in figure [ fig - spec ] , the spectrum can be fitted well by including only a power law , fe iiemission , h@xmath0 , [ o iii ] , and a few narrow lines .
* iras f12397@xmath73333 . *
as discussed in paper i , the spectrum of this agn is probably affected by host galaxy extinction .
the observed spectra ( see figure 4 of paper i for an example ) have a very red color , and are badly fitted by our spectral model .
we perform an extinction correction after de - redshifting , assuming the galactic extinction law and @xmath31 = 1.44 mag estimated from the balmer decrement ( see details in paper i ) .
this dereddened spectrum is fit well with a typical spectral index as shown in figure [ fig - spec ] .      for each individual - night spectrum ,
the fluxes , fwhms , and velocity shifts ( with respect to [ o iii ] ) of fe ii  and h@xmath0  are calculated from the best - fit model .
the means of these properties are considered to be the measurements of the line profile , except for the width of h@xmath0  ( @xmath32 ; see below ) .
the standard deviations on those means are used as the uncertainties .
@xmath32  is underestimated in the fitting , as we do not include the narrow component of h@xmath0  in our model ( see section [ sec - fithb ] ) .
thus , we use a different method to estimate this width .
we add a gaussian to our model to represent the narrow h@xmath0  component .
the velocity width and shift are constrained to be the same as those of [ o iii ] .
we then fit the mean spectra twice , once assuming the flux of the narrow h@xmath0  line is 10% of the flux of [ o iii ]  @xmath135007 and once assuming a flux ratio of 0.2 .
the width of the broad h@xmath0  obtained with flux ratio set to 0.1 is adopted as the @xmath32 .
the uncertainty is obtained from the fit with a flux ratio set to 0.2 and the original fit that assumed no narrow h@xmath0 .
this method essentially resembles the one used in papers i and ii , while the spectral model here is more sophisticated .
table [ tab - profile ] lists the measurements of the 10 objects .
the instrumental broadening ( fwhm @xmath24 500 @xmath2 , section [ sec - fithb ] ) has been taken into account in the listed fwhms .
note that the h@xmath0  of iras 04416@xmath71215  has large velocity shift with respect to [ o iii ]  @xmath135007 , and its profile shows no feature of a narrow component with the same shift of [ o iii ] . adding a narrow component
makes the width of the broad h@xmath0  even narrower .
thus , the @xmath32  of iras 04416@xmath71215 , as other properties , was obtained from the measurements of individual - night spectra .
lr@  @xmath33  lr@  @xmath33  lr@  @xmath33  lcr@  @xmath33  lr@  @xmath33  lr@  @xmath33  l mrk 335 & 253 & 9 & 1947 & 143 & 49 & 37 & & 661 & 22 & 2096 & 170 & @xmath81 & 16 + mrk 1044 & 369 & 11 & 866 & 25 & 76 & 10 & & 380 & 16 & 1178 & 22 & 12 & 12 + iras 04416@xmath71215 & 292 & 9 & 1313 & 50 & 496 & 21 & & 149 & 4 & 1522 & 44 & 241 & 29 + mrk 382 & 27 & 4 & 1326 & 234 & 13 & 115 & & 39 & 2 & 1462 & 296 & @xmath844 & 35 + mrk 142 & 87 & 5 & 1512 & 69 & @xmath825 & 37 & & 78 & 5 & 1588 & 58 & @xmath8101 & 28 + mcg @xmath706@xmath826@xmath8012 & 41 & 4 & 1155 & 70 & @xmath821 & 41 & & 41 & 4 & 1334 & 80 & @xmath829 & 23 + iras f12397@xmath73333 & 550 & 28 & 1748 & 78 & 36 & 40 & & 405 & 18 & 1802 & 560 & @xmath850 & 26 + mrk 42 & 67 & 3 & 787 & 16 & 105 & 19 & & 57 & 2 & 802 & 18 & 87 & 19 + mrk 486 & 186 & 5 & 1790 & 88 & 95 & 33 & & 346 & 12 & 1942 & 67 & @xmath846 & 9 + mrk 493 & 102 & 3 & 780 & 9 & 171 & 9 & & 92 & 3 & 778 & 12 & 126 & 13      our fitting successfully reduces the scatter in the light curves due to the influence of the host galaxy contamination ; thus , @xmath20  represents the real agn continuum much better than the simply integrated 5100   flux ( @xmath34 ) , as shown in appendix [ sec - host ] .
the light curves of the emission lines are also generated directly from the best - fit values of the corresponding parameters obtained from the fits of individual - night spectra .
the errors of the fluxes given by the fitting are not large enough to account for the scatter in the fluxes of successive nights , an additional systematic error is estimated for each light curve ( as in paper i ) .
this systematic error is added in quadrature to the fitting error for the calculations of variability amplitudes and time lags below .
note that our treatment is different from that of @xcite , who measured the light curve of h@xmath0  by integrating the continuum - subtracted spectra .
the left columns of figure [ fig - lc ] and figure [ fig - lc42 ] show the light curves of @xmath20 , @xmath28 , and @xmath25  for nine objects with reliable lag measurements , and the remaining one ( mrk 42 ) , respectively . in this section
, we will calculate the variability amplitudes and the reverberation lags for the h@xmath0  and fe ii  emission lines .
we then compare the h@xmath0  time lags with those presented in papers i and ii .      +          +            ( middle ) , and fe ii  ( bottom ) for mrk 42 .
note the similarity of the light curves of h@xmath0  and fe ii .
, scaledwidth=40.0% ]      we use the quantity @xmath35  defined in @xcite to represent the variability amplitude .
this quantity is an estimate of the intrinsic variability over the errors .
the uncertainties are calculated following @xcite .
the results for fe ii  and h@xmath0  are listed in table [ tab - var ] . for the variability amplitude ratio of fe ii  to h@xmath0 ,
the range is from about 0.6 ( mrk 486 ) to 1.2 ( mrk 42 ) , except for mrk 382 which has the largest host galaxy contamination . on average
( except for mrk 382 ) , the fe iivalue is about 10% smaller , consistent with previous results ( e.g. , * ? ?
? * ; * ? ? ?
lr@  @xmath33  lr@  @xmath33  l mrk 335 & 3.1 & 0.3 & 3.0 & 0.3 + mrk 1044 & 2.6 & 0.3 & 3.7 & 0.4 + iras 04416@xmath71215 & 2.1 & 0.3 & 2.0 & 0.3 + mrk 382 & 11.5 & 1.2 & 4.1 & 0.4 + mrk 142 & 5.5 & 0.5 & 6.6 & 0.5 + mcg @xmath706@xmath826@xmath8012 & 8.1 & 1.2 & 9.2 & 1.2 + iras f12397@xmath73333 & 4.4 & 0.6 & 4.1 & 0.5 + mrk 42 & 3.6 & 0.7 & 2.9 & 0.4 + mrk 486 & 2.0 & 0.5 & 3.4 & 0.4 + mrk 493 & 2.1 & 0.7 & 3.1 & 0.5      the time lags between the agn continuum variations ( @xmath20 ) and the emission lines ( @xmath25  and @xmath28 ) are measured from the cross - correlation functions ( ccfs ) for the relevant light curves .
we use the interpolation cross - correlation function ( iccf ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) method to calculate the ccf , and adopt the centroid of the ccf , above 80% of the peak value ( @xmath36 ) as the time lag @xcite .
the uncertainty in the time lag measurement is estimated from the cross - correlation centroid distribution ( cccd ) given by random subset selection / flux randomization monte carlo realizations @xcite .
all objects in our sample , except for mrk 42 , have reliable time lag measurements for both fe ii  and h@xmath0 .
the right columns of figure [ fig - lc ] show the results of the ccf analysis . for each object ,
the top - right panel shows the autocorrelation function ( acf ) of the @xmath20  light curve , which is shown in the top - left panel .
the two lower panels in the right column show the ccfs ( in black ) for h@xmath0  ( middle - right ) and fe ii  ( bottom - right ) with respect to @xmath20  light curves .
the blue histograms are the corresponding cccds .
most cccds have a rather symmetric profile , except those for iras 04416@xmath71215 .
table [ tab - lag ] lists the time lags of fe ii  ( @xmath37 ) and h@xmath0  ( @xmath38 ) , their uncertainties , and the corresponding @xmath36 , for the nine objects with reliable lag measurements .
the listed time lags and uncertainties are in the rest frame after time - dilation correction .
both time lags for iras 04416@xmath71215  have uncertainties of highly unequal upper and lower limits , as the sequence of their asymmetric cccds . for other objects ,
the time lags are well determined by the single - peak ccfs with high @xmath36  and symmetric , narrow cccds .
lccccc mrk 335 & 0.48 & @xmath39 & & 0.70 & @xmath40 + mrk 1044 & 0.48 & @xmath41 & & 0.62 & @xmath42 + iras 04416@xmath71215 & 0.61 & @xmath43 & & 0.60 & @xmath44 + mrk 382 & 0.40 & @xmath45 & & 0.47 & @xmath46 + mrk 142 & 0.78 & @xmath47 & & 0.84 & @xmath48 + mcg @xmath706@xmath826@xmath8012 & 0.86 & @xmath49 & & 0.92 & @xmath50 + iras f12397@xmath73333 & 0.63 & @xmath51 & & 0.71 & @xmath52 + mrk 486 & 0.77 & @xmath53 & & 0.79 & @xmath54 + mrk 493 & 0.63 & @xmath55 & & 0.81 & @xmath56    the detection rate of fe ii  time lag in our sample is extremely high ( 9 in 10 ) , compared with previous reverberation - mapping experiments .
this is mainly attributed to the common feature of strong fe ii  emission in the seambh sample , the high - cadence observation , and sophisticated spectral fitting .
the failure to obtain reliable @xmath37  and @xmath38  in mrk 42 is caused by the @xmath20light curve , which has large scatter and no clear structure ( see figure [ fig - lc42 ] ) .
however , the light curves of @xmath25  and @xmath28  show hints for a similar structure which , unfortunately , is not enough to establish a time lag .    the clear reverberation of the fe ii  emission in response to the continuum supports the hypothesis that the fe ii  emission originates from photoionized gas in the blr .
the high detection rate indicates that this is prevalent in nls1 , which are agns with high @xmath3 .
it is possible that the origin of fe iiemission in agns with low @xmath3  is different , but there is little evidence to support this claim .
in fact , @xcite has already detected the reverberation of fe ii  in two broad - line seyfert 1 galaxies with h@xmath0  widths of fwhm @xmath57 4000 @xmath2  and much lower @xmath3 .
comparing the time lags of fe ii  and h@xmath0 , along with their velocity widths , provides some hints about the location and geometries of the two emission regions .
we will discuss these ideas in section [ sec - cmpfehb ] below .      in paper
i and ii , the reported time lags of h@xmath0  are obtained from the light curves measured by the more traditional integration method , without taking into account the contamination by the host galaxy and the narrow emission lines .
figure [ fig - hbsumcmp ] shows a comparison between the previously reported @xmath38  based on the integration method and the new results given by the fitting method reported here .
the dashed diagonal line denotes the 1:1 ratio and the region between the two dotted lines gives the 0.1 dex deviation from the line . for seven out of nine objects ,
the differences between @xmath38given by the two methods are less than 0.1 dex .
the two exceptions are iras 04416@xmath71215and mrk 1044 . in paper ii we show that using the integration method , we can not obtain a significant @xmath38  for iras 04416@xmath71215 . the lag for mrk 1044 reported in that paper
is about a factor 2 shorter than obtained here , with a very large uncertainty ( @xmath58 ) .
this value is consistent with the new one presented here within the uncertainties .
in addition , the uncertainty on @xmath38 in mrk 493 given by the integration method is very large with a negative lower limit ( @xmath59 ) .
the fitting method yields much smaller uncertainty with basically the same @xmath38 .
this is a direct result of the much smoother h@xmath0  light curve of mrk 493 given by the fitting method ( see figure [ fig - lc ] , and also figure 2 of paper ii ) .
thus , the time lags of h@xmath0given by the two methods are consistent with each other and @xmath38  is better defined by the fitting procedure in some of the cases .
time lags obtained by simple integration ( papers i and ii ) with those obtained by fitting ( this work ) .
the differences are typically less than 0.1 dex ( dotted lines ) except for iras 04416@xmath71215  ( not plotted because its @xmath38  can not be defined by integration ) and mrk 1044 ( very large uncertainty of the integration method ) .
mrk 493 also has very large uncertainty when measured by the integration method but the measured @xmath38 are very similar . , scaledwidth=45.0% ]
a simple theoretical expectation based on photoionization is that @xmath60 with @xmath61 .
the @xmath5@xmath6  relation for h@xmath0  has been compared with such predictions in several earlier publications ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* and references therein ) .
regarding fe ii  lines , @xcite provide a tentative @xmath5@xmath6  relation for the first time , from a small inhomogeneous sample of six agns . among them , two objects are reported by @xcite using the standard spectral fitting method , four others come from @xcite and @xcite using the mcf scheme of @xcite . here
, we revisit this issue by more than doubling the size of the sample by adding our nine newly measured sources . for a proper comparison with the @xcite results we used our best observed fluxes and uncertainty and a cosmology with @xmath62 mpc@xmath10 , @xmath63 , and @xmath64 .
figure [ fig - ferl ] shows @xmath37  vs. the 5100   luminosity for the nine objects in our sample ( black dots ) alongside the earlier results ( kindly provided by d. chelouche ) .
the green squares are those obtained by the mcf scheme ( including three with insignificant results , see figure 4 of @xcite for details ) .
the blue triangles are the two objects from @xcite for which we used the centroid of the published ccfs for fe ii  and the @xmath18-band light curves . using the fitexy method @xcite to fit all the 18 objects we find : @xmath65 which is plotted in a solid line .
this line deviates considerably from the relation presented in @xcite ( shown in a dotted line ) , but close to that of @xcite for h@xmath0  ( dashed line ) .
this result is expected since most of the objects in our sample show roughly equal @xmath37  and @xmath38(see table [ tab - lag ] ) , except for three objects ( mrk 335 , mrk382 , and mrk486 ) , while the previous studies have @xmath37  greater than @xmath38 .
our sample by itself shows no correlation at all between @xmath37  and the 5100 luminosity .
this is not surprising given the relatively small range in luminosity considered by us and the fact that for this sample , the correlation between @xmath38  and the continuum luminosity is also very weak .
we suggest that this lack of correlation is related also to the large @xmath3  of the objects in our sample , a topic which we discuss in great detail in paper iv ( du et al . ,
submitted ) .
@xmath6  relation .
the nine objects in our sample are plotted as black dots , the three objects with their names labeled have relatively large differences between @xmath38  and @xmath37 .
measurements from previous studies are plotted as blue triangles and green squares ( see text for details ) .
the solid line is the fit to the entire data set .
the dotted and dashed lines are the @xmath5@xmath6  relations in @xcite for fe ii  and @xcite for h@xmath0 , respectively .
, scaledwidth=45.0% ]      as shown in table [ tab - lag ] , @xmath37  is roughly equal to @xmath38  ( counting the uncertainties in both lags ) in six of our objects , longer in two ( mrk 335 and mrk 382 ) and somewhat shorter in one ( mrk 486 ) . a common feature of the three objects with unequal lags is the relatively low fe ii / h@xmath0  intensity ratio .
figure [ fig - rfelag ] shows the plot of @xmath37/@xmath38  vs. @xmath25/@xmath28  ( defined as @xmath66 ) .
the three objects with different time lags are labeled by their names and locate in the region of @xmath66  @xmath1 1 .
the six other objects show @xmath66  @xmath4 1 .
we also mark the approximate positions of the two objects measured by @xcite .
defined here is about 3/4 that in @xcite . ]
based on eye estimates of the published light curves shown in their paper , both objects have @xmath66  @xmath24 0.75 .
their positions in figure [ fig - rfelag ] are marked around 0.75 ( and horizontally shifted for clarity ) .
they both have different time lags than h@xmath0  and @xmath66  @xmath1 1 , consistent with the results of our sample .    vs. intensity ratios of fe ii  to h@xmath0 .
objects in our sample are shown as black dots .
the three objects labeled by names have different time lags of fe ii  and h@xmath0 , and all have low intensity ratio of fe ii  to h@xmath0 .
blue triangles show the approximate positions of the two objects in @xcite ( intensity ratios are estimated by eye ) .
the two points are slightly shifted horizontally for clarity ( see text for details . ) , scaledwidth=45.0% ]    we note , again , that our sample is highly biased toward agns with high @xmath3and thus high @xmath66 . in an unbiased agn sample , the fraction of sources with @xmath66  @xmath1 1 is much larger ( see , e.g. , figure 7 of @xcite and figure 1 of @xcite ) .
@xcite analysed the properties of a large number of low - redshift agns .
they found that the velocity width of fe ii  is systematically narrower ( @xmath67 ) than that of h@xmath0 , which may suggest , given keplerian velocities , that in those sources @xmath37  is longer than @xmath38 .
our present sample contains only 10 sources which is too few to test , systematically , any of these ideas .
we provide new fe ii  measurements for 10 nls1s and report statistically significant time lag measurements for nine of these sources .
all the observed nls1s are suspected to be seambhs .
this more than doubles the number of agn showing measurable fe ii  time lags .
our time lag measurements are based on high - cadence , high - s / n measurements at the yunnan observatory , and on a new , sophisticated , fitting analysis that takes into account the uncertainties caused by the apparent flux change of the host galaxy and several narrow emission lines .
we demonstrate , by a careful comparison with earlier measurements of @xmath38 , that the method can considerably improve the accuracy of the time lag measurement .
the main findings reported in this paper can be summarized as follows :    1 .
all 10 objects presented here show fe ii  variations with an amplitude of a few to ten percent . on average ,
this is about 10% smaller than the variability of the h@xmath0  line . 2 .
reliable fe ii  reverberation time lags with respect to the agn continuum are detected in nine objects , confirming the suggestion that the fe ii  emission originates from photoionized gas .
3 .   combining the new reverberation mapping results with those in previous work , shows a clear radius
luminosity relationship for fe ii  that is similar to the one for h@xmath0
. however , our sample by itself shows no such correlation due to the large intrinsic scatter over a small luminosity range .
the difference in the time lags of fe ii  and h@xmath0  depends on the intensity ratio of fe ii  to h@xmath0  ( @xmath66 ) .
the time lag of fe ii  is roughly equal to that of h@xmath0  in all the six objects with @xmath66  @xmath57 1 .
the fe ii  time lag is longer in mrk 335 and mrk 382 , those objects with @xmath66  @xmath1 1 , and shorter in mrk 486 with @xmath66  @xmath1 1 .
we acknowledge the support of the staff of the lijiang 2.4 m telescope .
funding for the telescope has been provided by cas and the people s government of yunnan province .
this research is supported by the strategic priority research program
the emergence of cosmological structures of the chinese academy of sciences , grant no .
xdb09000000 , by the nsfc through nsfc-11173023 , -11133006 , -11233003 , -11473002 , and by israel - china isf - nsfc grant 83/13 .
our flux calibration method provides the _ relative _ flux of the object with respect to that of the local comparison star . for two point sources kept in a line parallel to the slit , the fractions of light loss due to seeing , differential atmospheric refraction , and mis - centering are identical .
thus this method gives accurate relative flux calibration for the spectral components of the agn , including the featureless power law ( @xmath20 ) and the broad emission lines .
this is not the case for the extended host galaxies whose flux relative to the comparison star inside the slit can vary much more due to seeing variations and mis - centering .
in addition , different parts of the galaxy may be observed each time . as a result ,
the derived flux of the host galaxy is @xmath68 , where @xmath69 is the absolute flux ( which is constant ) , and @xmath70 is a factor accounting for all the effects of varying observing conditions in the flux calibration procedure .
@xmath70 could change from one exposure to the next , introducing an additional uncertainty into the integrated 5100 flux @xmath71 .    in papers i and ii , we estimated the flux of the host galaxies in the spectral extraction aperture using archival _ hst _
images for eight objects in our sample .
the resultant relative fluxes , ( @xmath23/@xmath20 ) , range from @xmath240.2 in mrk 335 to @xmath57 1 in mrk 382 .
thus , the apparent change of flux of the host galaxy is not negligible for those objects with strong host galaxy contribution . for mrk 382 ,
the measurement of the agn continuum was badly affected by this uncertainty , which forced us to use the @xmath18-band photometry instead . the new procedure adopted here ( see the main text ) enables us to solve for @xmath72 for each individual - night spectrum .
the detailed fitting of the various spectral components considerably reduces the noise in the various line and continuum light curves . in particular , for mrk 382 , we can now recover the 5100   agn
continuum variations to a much better accuracy , which is evident by the good agreement with the @xmath18-band photometry .
the big improvements due to the spectral fitting procedure are shown in figure [ fig - lc382host ] , which provides more information on the process for mrk 382 .
the top panel shows the light curve of @xmath34 , which is measured as in papers i and ii , except for the inclusion of galactic extinction correction and de - redshifting . the scatter in this light curve is very large , and hence it was replaced in paper ii by the @xmath18-band light curve .
the upper - middle and lower - middle panels show the results of the spectral fitting for @xmath23  and @xmath20 , respectively .
the light curve of @xmath20  is basically identical to that of the @xmath18-band ( bottom panel ) .
clearly , the unrealistic high @xmath34  around jd 2456300 was caused by the large contribution of the host galaxy caused by seeing and guiding fluctuations .
on other jds , the host contamination is relatively small , but still strong enough to affect the variability of @xmath20 , e.g. , the peak around jd 2456250 .    ,
top panel ) , fitted host galaxy ( @xmath72 described in appendix [ sec - host ] , upper - middle panel ) , fitted agn power - law continuum ( @xmath20 , lower - middle panel ) , and @xmath18-band magnitude ( bottom panel , see also paper ii ) .
note that @xmath34  is heavily contaminated by the host whose apparent flux is varying by @xmath70 , and the @xmath20  recovers , accurately , the intrinsic agn variability as judged by a comparison with the @xmath18-band observations . , scaledwidth=40.0% ]
we present the fe ii  light curves measured by the traditional integration method for a comparison with that by the present fitting scheme in this appendix .
the flux of the fe ii  bump on the red side of [ o iii ]  @xmath135007 is calculated by integrating in 51155465   the flux above the straight line set by two continuum windows 50855115 and 54655495   ( the three bands are in the rest frame ) .
figure [ fig - feintg ] shows these light curves . comparing with those given by the fitting method ( in figures [ fig - lc ] and [ fig - lc42 ] ) , all the light curves by integration show larger scatter . for four objects ( mrk 1044 , mrk 142 , mcg @xmath706@xmath826@xmath8012 , and iras f12397@xmath73333 ) , rough structures can be recognized in their light curves , and the derived fe iitime lags are consistent with those given by the fitting but with larger uncertainties .
the remaining objects fail to have reliable measurement of time lag , for their scattered light curves without clear structures .
we also tried to measure the flux of the fe ii  bump on the blue side of h@xmath0  by choosing an integration interval avoiding the broad he ii  line .
similarly , light curves with poor quality are obtained .
baldwin , j. a. , ferland , g. j. , korista , k. t. , hamann , f. , & lacluyz , a. 2004 , , 615 , 610 barth , a. j. , pancoast , a. , bennert , v. n. , et al .
2013 , , 769 , 128 barth , a. j. , pancoast , a. , thorman , s. j. , et al .
2011 , , 743 , l4 bentz , m. c. , denney , k. d. , grier , c. j. , et al .
2013 , , 767 , 149 bentz , m. c. , peterson , b. m. , netzer , h. , pogge , r. w. , & vestergaard , m. 2009 , , 697 , 160 bentz , m. c. , walsh , j. l. , barth , a. j. , et al .
2010 , , 716 , 993 bian , w .- h .
, huang , k. , hu , c. , et al . 2010 , , 718 , 460 boksenberg , a. , & netzer , h. 1977 , , 212 , 37 boller , t. , brandt , w. n. , & fink , h. 1996 , , 305 , 53 boroson , t. a. 2002 , , 565 , 78 boroson , t. a. , & green , r. f. 1992 , , 80 , 109 bruzual , g. , & charlot , s. 2003 , , 344 , 1000 cardelli , j. a. , clayton , g. c. & mathis , j. s. 1989 , , 345 , 245 chelouche , d. , rafter , s. e. , cotlier , g. i. , kaspi , s. , & barth , a. j. 2014 , , 783 , l34 chelouche , d. , & zucker , s. 2013 , , 769 , 124 collin - souffrin , s. , joly , m. , dumont , s. , & heidmann , n. 1980 , , 83 , 190 collin , s. , kawaguchi , t. , peterson , b. m. , & vestergaard , m. 2006 , , 456 , 75 du , p. , hu , c. , lu , k .- x . , et al .
2014 , , 782 , 45  ( paper i ) edelson , r. , turner , t. j. , pounds , k. , et al .
2002 , , 568 , 610 ferland , g. j. , hu , c. , wang , j .-
, et al . 2009 , , 707 , l82 gaskell , c. m. , & sparke , l. s. 1986 , , 305 , 175 gaskell , c. m. , & peterson , b. m. 1987 , , 65 , 1 grier , c. j. , peterson , b. m. , pogge , r. w. , et al . 2012 , , 744 , l4 ho , l. c. , & kim , m. 2009 , , 184 , 398 hu , c. , wang , j .-
, ho , l. c. , et al .
2008 , , 687 , 78 hu , c. , wang , j .- m . , ho , l. c. , et al .
2012 , , 760 , 126 joly , m. 1987 , , 184 , 33 kaspi , s. , maoz , d. , netzer , h. , et al .
2005 , , 629 , 61 kaspi , s. , smith , p. s. , netzer , h. , et al .
2000 , , 533 , 631 kollatschny , w. , bischoff , k. , & dietrich , m. 2000 , , 361 , 901 koratkar , a. p. , & gaskell , c. m. 1991 , , 75 , 719 kuehn , c. a. , baldwin , j. a. , peterson , b. m. , & korista , k. t. 2008 , , 673 , 69 maoz , d. , & netzer , h. 1989 , , 236 , 21 maoz , d. , netzer , h. , leibowitz , e. , et al .
1990 , , 351 , 75 maoz , d. , netzer , h. , peterson , b. m. , et al .
1993 , , 404 , 576 netzer , h. 1990 , in active galactic nuclei , ed .
r. d. blandford , h. netzer , & l. woltjer ( berlin : springer ) , 137 netzer , h. , & wills , b. j. 1983 , , 275 , 445 odonnell , j. e. 1994 , , 422 , 158 onken , c. a. , ferrarese , l. , merritt , d. , et al .
2004 , , 615 , 645 osterbrock , d. e. , & pogge , r. w. 1985 , , 297 , 166 peterson , b. m. , denney , k. d. , de rosa , g. , et al .
2013 , , 779 , 109 peterson , b. m. , ferrarese , l. , gilbert , k. m. , et al .
2004 , , 613 , 682 peterson , b. m. , wanders , i. , horne , k. , et al .
1998 , , 110 , 660 press , w. h. , teukolsky , s. a. , vetterling , w. t. , & flannery , b. p. 1992
, numerical recipes in fortran ( 2nd ed .
; cambridge : cambridge univ .
press ) rafter , s. e. , kaspi , s. , chelouche , d. , et al .
2013 , , 773 , 24 richards , g. t. , kruczek , n. e. , gallagher , s. c. , et al .
2011 , , 141 , 167 rodrguez - pascual , p. m. , alloin , d. , clavel , j. , et al .
1997 , , 110 , 9 schlafly , e. f. , & finkbeiner , d. p. 2011
, , 737 , 103 shapovalova , a. i. , popovic , l. c. , burenkov , a. n. , et al . 2012 , , 202 , 10 shen , y. , & ho , l. c. 2014 , , 513 , 210 sigut , t. a. a. , & pradhan , a. k. 1998 , , 499 , l139 sulentic , j. w. , marziani , p. , & dultzin - hacyan , d. 2000 , , 38 , 521 vanden berk , d. e. , richards , g. t. , bauer , a. , et al .
2001 , , 122 , 549 van der marel , r. p. , & franx , m. 1993 , , 407 , 525 vron - cetty , m .- p . , joly , m. , & vron , p. 2004
, , 417 , 515 vestergaard , m. , & peterson , b. m. 2005 , , 625 , 688 wang , j. , wei , j. y. , & he , x. t. 2005 , , 436 , 417 wang , j .-
m . , du , p. , hu , c. , et al .
2014 , , 793 , 108  ( paper ii ) wang , j .-
du , p. , valls - gabaud , d. , hu , c. , & netzer , h. 2013 , , 110 , 81301 wang , j .- m . , watarai , k .- y . , & mineshige , s. 2004 , , 607 , l107 white , r. j. , & peterson , b. m. 1994 , , 106 , 879 wills , b. j. , netzer , h. , & wills , d. 1985 , , 288 , 94 woo , j .- h . ,
schulze , a. , park , d. , et al .
2013 , , 772 , 49 york , d. g. , adelman , j. , anderson , j. e. , et al .
2000 , , 120 , 1579 zhang , x .- g .
2013 , , 435 , 2141 zhou , h. , wang , t. , yuan , w. , et al .
2006 , , 166 , 128