SDSS-IV/SPIDERS: A Catalogue of X-Ray Selected AGN Properties

SDSS-IV/SPIDERS: A Catalogue of X-Ray Selected AGN Properties

Spectral Properties and Black Hole Mass Estimates for SPIDERS SDSS DR14 Type 1 AGN
D. Coffey    M. Salvato    A. Merloni    Th. Boller    K. Nandra    T. Dwelly    J. Comparat    A. Schulze    A. Del Moro    D. P. Schneider
Accepted 8 April 2019
Key Words.:
quasars: emission lines – quasars: general – galaxies: active – catalogues – surveys
11institutetext: Max-Planck-Institut für extraterrestrische Physik, Giessenbachstraße 1, 85748 Garching, Germany
11email: coffeydg@mpe.mpg.de
22institutetext: National Astronomical Observatory of Japan, Mitaka, Tokyo 181-8588, Japan 33institutetext: Department of Astronomy and Astrophysics, and the Institute for Gravitation and the Cosmos, The Pennsylvania State University, 525 Davey Laboratory, University Park, PA 16802, USA

This work presents the catalogue of optical spectral properties for all X-ray selected SPIDERS active galactic nuclei (AGN) up to SDSS DR14. SPIDERS (SPectroscopic IDentification of eROSITA Sources) is an SDSS-IV programme that is currently conducting optical spectroscopy of the counterparts to the X-ray selected sources detected in the ROSAT all-sky survey and the XMM-Newton slew survey in the footprint of the Extended Baryon Oscillation Spectroscopic Survey (eBOSS). The SPIDERS DR14 sample is the largest sample of X-ray selected AGN with optical spectroscopic follow-up to date. The catalogue presented here is based on a clean sample of 7344 2RXS () and 1157 XMM-Newton slew survey () type 1 AGN with spectroscopic coverage of the H and/or MgII emission lines. Visual inspection results for each object in this sample are available from a combination of literature sources and the SPIDERS group, which provide both reliable redshifts and source classifications. The spectral regions around the H and MgII emission lines have been fit in order to measure both line and continuum properties, estimate bolometric luminosities, and provide black hole mass estimates using the single-epoch (or photoionisation) method. The use of both H and MgII allows black hole masses to be estimated up to z  2.5. A comparison is made between the spectral properties and black hole mass estimates derived from H and MgII using the subsample of objects which have coverage of both lines in their spectrum. These results have been made publicly available as an SDSS-IV DR14 value added catalogue.

1 Introduction

A crucial requirement for understanding AGN evolution and demographics is the ability to select a sample for study in a complete and unbiased way. X-ray emission has been frequently used for AGN selection, and can distinguish the high energy emission associated with mass accretion by a black hole (BH) from inactive galaxies and stars. Combining wide-area X-ray surveys with the ability to classify large numbers of objects spectroscopically via the Sloan Digital Sky Survey (SDSS, 2000AJ....120.1579Y; 2006AJ....131.2332G) provides a powerful tool for the study of AGN.

SPIDERS (SPectroscopic IDentification of eROSITA Sources; PIs Merloni and Nandra) is an SDSS-IV (2017AJ....154...28B) eBOSS (2016AJ....151...44D) subprogramme that is currently conducting optical spectroscopy of extragalactic X-ray detections in wide-area ROSAT and XMM-Newton surveys (2017MNRAS.469.1065D). Lying at the bright end of the X-ray source population, these sources will also be detected by eROSITA (2012arXiv1209.3114M; 2016SPIE.9905E..1KP). The current SDSS DR14 (2018ApJS..235...42A) SPIDERS sample is a powerful resource for the multiwavelength analysis of AGN. This work aims to capitalise on the wealth of information already available by providing detailed optical spectral measurements, as well as estimates of BH masses and Eddington ratios.

An accurate measurement of the central supermassive black hole (SMBH) mass is necessary for the study of AGN and their coevolution with their host galaxies. BH mass has been found to scale with a number of host galaxy spheroid properties; stellar velocity dispersion (the relation, e.g. 2000ApJ...539L..13G; 2001ApJ...547..140M; 2002ApJ...574..740T), stellar mass (e.g. 1998AJ....115.2285M), and luminosity (e.g. 1995ARA&A..33..581K). These correlations suggest a symbiotic evolution of SMBHs and their host galaxies.

Reverberation mapping (RM) has been used to measure the approximate radius of the broad-line region (BLR) in AGN (e.g. 1972ApJ...171..467B; 1982ApJ...261...35C; 1982ApJ...255..419B; 1993PASP..105..247P; 2015PASP..127...67B; 2015ApJS..216....4S). This technique involves measuring the time delay between variations in the continuum emission, which is expected to arise from the accretion disk, and the induced variations in the broad emission lines. It was found that different emission lines have different time delays, which is expected if the BLR is stratified, with lines of lower ionisation being emitted further from the central ionising source (e.g. 1986ApJ...305..175G). For example, the high ionisation line CIV1549 has a shorter time delay than H (2000ApJ...540L..13P).

The RM effort has also revealed a tight relationship between the continuum luminosity and the radius of the BLR (2000ApJ...533..631K; 2006ApJ...644..133B; 2009ApJ...697..160B). Therefore, by using the measured luminosity as a proxy for the BLR radius, and measuring the BLR line-of-sight velocity from the width of the broad emission lines, BH masses can be estimated from a single spectrum (2002ApJ...571..733V; 2002MNRAS.337..109M; 2006ApJ...641..689V; 2011ApJ...742...93A; 2012ApJ...753..125S; 2013BASI...41...61S). This approach is known as the single-epoch, or photoionisation, method.

Since H is the most widely studied RM emission line it is therefore considered to be the most reliable line to use for single-epoch mass estimation. In addition, AGN H emission lines typically exhibit a clear inflection point between the broad and narrow line components, making the virial full width at half maximum (FWHM) measurement relatively straightforward (see section 3.2.1). The MgII line width correlates well with that of H (see section 7.2), and therefore MgII has also been used for single-epoch mass estimation (e.g. 2002MNRAS.337..109M). For SDSS spectra, either H or MgII is visible in the redshift range 0  z  2.5.

At higher redshifts, the broad, high-ionisation line CIV is available. The CIV line width does not correlate strongly with that of low ionisation lines (e.g. 2005MNRAS.356.1029B; 2012MNRAS.427.3081T) and this, along with the presence of a large blueshifted component (e.g. 2002AJ....124....1R) makes it difficult to employ CIV for mass estimation. A number of calibrations have been developed which aim to improve the mass estimates derived from CIV (2012ApJ...759...44D; 2013MNRAS.434..848R; 2013ApJ...770...87P; 2017MNRAS.465.2120C), however, whether CIV can provide reliable mass estimates when compared with low ionisation lines is still a subject of debate (see 2018MNRAS.478.1929M).

This paper is organised as follows: the selection of a reliable subsample of sources to be used for optical spectral fitting is discussed in section 2. Section 3 describes the method used to fit the H and MgII emission line regions. The methods for estimating BH mass and bolometric luminosity are discussed in section 4. The X-ray flux measurements used in this work are discussed in section 5. Section 6 describes where the catalogue containing the results of this work can be accessed. A comparison between the UV and optical spectral fitting results is given in section 7. Section 8 provides a discussion of the sample properties, and finally, section 9 includes a discussion of the reliability and limitations of the fitting procedure. In order to facilitate a direct comparison with previous studies based on X-ray surveys, a concordance flat cosmology was adopted where =0.3, =1-, and H=70 km s Mpc.

2 Preparing the Input Catalogue

2.1 X-ray Data

The ROSAT sample used in this work is part of the second ROSAT all-sky survey (2RXS) catalogue (2016A&A...588A.103B), which has a limiting flux of , which corresponds to a luminosity of at z = 0.5. Compared to the first ROSAT data release (1999A&A...349..389V), the 2RXS catalogue is the result of an improved detection algorithm, which uses a more detailed background determination relative to the original ROSAT pipeline. A full visual inspection of the 2RXS catalogue has been performed, which provides a reliable estimate of its spurious source content (see 2016A&A...588A.103B). The first XMM-Newton slew survey catalogue release 1.6 (XMMSL1; 2008A&A...480..611S) was also used in this work. This catalogue includes observations made by the European Photon Imaging Camera (EPIC) pn detectors while slewing between targets, and has a limiting flux of in the soft band, which corresponds to a luminosity of at z=0.5.

Figure 1: Soft X-ray luminosity versus spectroscopic redshift for the samples presented in this work and the following previously published X-ray selected samples; XMM-XXL (2016MNRAS.457..110M; 2016MNRAS.459.1602L), CDFS (2017ApJS..228....2L), STRIPE82X (2016ApJ...817..172L), COSMOS (2016ApJ...817...34M; 2016ApJ...830..100M; 2016ApJ...819...62C), AEGIS-X (2015ApJS..220...10N), and the Lockman Hole deep field (LHDF) (2008A&A...479..283B; 2012ApJS..198....1F). For each sample, the 0.5-2 keV luminosities are shown, except for the 2RXS sample, where the 0.1-2.4 keV luminosities are shown, and the XMMSL1 sample, where the 0.2-2 keV luminosities from 2008A&A...480..611S are shown. The detection limit for the 2RXS and XMMSL1 samples are shown by the solid and dashed grey lines respectively. The X-ray luminosities for the 2RXS sample are derived from the classical flux estimates described in section 5.1, however it is noted here that some low count rate 2RXS sources do not have flux estimates. For sources that were detected in both 2RXS and XMMSL1, only the XMMSL1 luminosities are shown. Sources classified as stars have not been included in this figure.

2.2 The SPIDERS Programme

The SPIDERS programme has been providing SDSS spectroscopic observations of 2RXS and XMMSL1 sources111The SPIDERS programme targets both point-like and extended X-ray sources. This work focuses on the counterparts to point-like X-ray detections, which are predominantly AGN, and therefore, the samples discussed in the subsequent paragraphs are derived from the SPIDERS-AGN programme. in the eBOSS footprint. Before the start of the eBOSS survey in 2014, the SPIDERS team compiled a sample of X-ray selected spectroscopic targets and submitted this sample for spectroscopic follow-up using the BOSS spectrograph as part of the eBOSS/SPIDERS subprogramme (see 2017MNRAS.469.1065D, for further details on the SPIDERS programme). As of the end of eBOSS in February 2019, the eBOSS/SPIDERS survey has covered a sky area of . The SDSS DR14 SPIDERS sample presented in this work covers an area of ( of the final eBOSS/SPIDERS area).

The spectroscopic completeness achieved by the SPIDERS survey as of SDSS DR14 in the eBOSS area is for the sample as a whole, considering only high-confidence X-ray detections (see section 2.5), and considering sources with high-confidence X-ray detections and optical counterparts with magnitudes in the nominal survey limits (). Outside the eBOSS area, the spectroscopic completeness of this sample is lower: for the sample as a whole, considering only high-confidence X-ray detections, and considering sources with high-confidence X-ray detections and optical counterparts with magnitudes in the nominal survey limits. The spectroscopic completeness of the SDSS DR16 SPIDERS sample inside and outside the eBOSS area is expected to be similar to that of the sample presented here.

In addition to those targeted during eBOSS/SPIDERS, a large number of 2RXS and XMMSL1 sources received spectra during the SDSS-I/II (2000-2008; 2000AJ....120.1579Y) and the SDSS-III (2011AJ....142...72E) BOSS (2009-2014; 2013AJ....145...10D) surveys. This paper includes spectra obtained by eBOSS/SPIDERS up to DR14 (2014-2016) as well as spectra from SDSS-I/II/III.

2.3 Identifying IR Counterparts

To identify SPIDERS spectroscopic targets, the Bayesian cross-matching algorithm “NWAY” (2018MNRAS.473.4937S) was used to select AllWISE (2013yCat.2328....0C) infrared (IR) counterparts for the 2RXS and XMMSL1 X-ray selected sources in the BOSS footprint. The AllWISE catalogue consists of data obtained during the two main survey phases of the Wide-field Infrared Survey Explorer mission (WISE; 2010AJ....140.1868W) which conducted an all-sky survey in the 3.4, 4.6, 12, and 22 bands (magnitudes in these bands are denoted [W1], [W2], [W3], and [W4] respectively). The matching process used the colour-magnitude priors [W2] and [W2-W1] (see 2017MNRAS.469.1065D) which, at the depth of the 2RXS and XMMSL1 surveys, can distinguish between the correct counterparts and chance associations. These colours would not be efficient if the 2RXS survey was much deeper (see 2018MNRAS.473.4937S, for a complete discussion). The resulting 2RXS and XMMSL1 catalogues with AllWISE counterparts contained 53455 and 4431 sources respectively. AllWISE positions were then matched to photometric counterparts, where available, in SDSS.

2.4 Comparison with Previous X-ray Surveys

Figure 1 displays the sources in the 2RXS and XMMSL1 samples which have spectroscopic redshifts and measurements of the soft X-ray flux. For comparison, a series of previously published X-ray selected samples that have optical spectroscopic redshifts are also shown. The large number of sources present in the 2RXS and XMMSL1 samples motivated the optical spectroscopic analysis discussed in the following sections.

Figure 2: Sequence of quality cuts applied to the 2RXS and XMMSL1 samples to produce the subsample used for spectral analysis. The starting points (2RXS-AllWISE and XMMSL1-AllWISE) are the full samples of 2RXS and XMMSL1 selected sources with AllWISE IR counterparts in the BOSS footprint (see section 2.3). The steps in grey are those that have been discussed in 2017MNRAS.469.1065D.

2.5 Selecting a Reliable Subsample

The selection of SPIDERS spectroscopic targets was discussed in detail by 2017MNRAS.469.1065D. This section summarises the selection steps discussed in detail by 2017MNRAS.469.1065D and describes the additional cuts made in this work to select a sample for spectral analysis. The sequence of selection criteria used and the resulting sample size are shown in figure 2.

2RXS sources with an X-ray detection likelihood (EXI_ML) 10 were excluded since these detections are considered highly uncertain with a spurious fraction (see 2016A&A...588A.103B). XMMSL1 sources with an X-ray detection likelihood (XMMSL_DET_ML_B0) 10 were also excluded. 2018MNRAS.473.4937S show the distribution of flux with detection likelihood for both samples. These cuts returned 23245/53455 2RXS and 3803/4431 XMMSL1 sources.

The following cuts, which were described in detail in 2017MNRAS.469.1065D, have also been applied to the sample:

  • For each X-ray source, 2018MNRAS.473.4937S give the probability, p_any, that a reliable counterpart exists among the possible AllWISE associations. Sources with p_any ¡ 0.01 were removed. This returned 23046/23245 2RXS and 3558/3803 XMMSL1 sources.

  • For each X-ray source, the most probable AllWISE counterpart was chosen by selecting sources with match_flag=1. This returned 20585/23046 2RXS and 3321/3558 XMMSL1 sources.

  • For each AllWISE counterpart, the brightest SDSS-DR13 photometric source within the AllWISE matching radius was selected using FLAG_SDSSv5b_best=1. This returned 19385/20585 2RXS and 3063/3321 XMMSL1 sources.

  • Cases where the AllWISE-SDSS separation exceeded 3 arcsec were removed. This returned 18575/19385 2RXS and 2893/3063 XMMSL1 sources.

The results of these constraints are displayed in grey in figure 2. As shown above, sources with match_flag=1 were targeted; however, for 14% of the 2RXS sample and 10% of the XMMSL1 sample, more than one counterpart was highly likely. This implies that either the counterpart association was not reliable, or that the X-ray detection was the result of emission from multiple sources. These sources were not included in the discussion of optical spectral properties as a function of X-ray properties in section 8.3. After selecting the brightest SDSS-DR13 (2017ApJS..233...25A) photometric source, there were 15 cases where two unique 2RXS sources were matched to the same AllWISE/SDSS counterpart and 3 cases where two unique XMMSL1 sources were matched to the same AllWISE/SDSS counterpart. These sources were also removed.

Of these samples with reliable SDSS photometric counterparts, 8777 2RXS and 1315 XMMSL1 sources have received spectra during SDSS-I/II/III while 1122 2RXS and 221 XMMSL1 sources have received spectra during the SPIDERS programme (including SEQUELS), resulting in a sample of 9899 2RXS and 1536 XMMSL1 sources with spectra as of DR14. The distribution of SDSS i band fiber2 magnitudes for this sample (showing the different spectroscopic programmes) is presented in figure 3. Due to targeting constraints (as discussed in section 2.2), the sample completeness is much lower outside of the nominal magnitude limits for the survey ( for eBOSS).

Figure 3: Distribution of i-band fiber magnitudes (fiber2Mag). The coloured curves represent all of the sources with spectra, and the survey from which the spectra were taken. The grey histogram displays the X-ray sources with a reliable SDSS photometric counterpart, including stars which cannot be targets for spectroscopy due to their brightness.

2.6 Source Classification

Visual inspection results for each object in this sample are available from a combination of literature sources (2007AJ....133..313A; 2010AJ....139..390P; 2010AJ....139.2360S; 2017A&A...597A..79P) and the SPIDERS group. The SPIDERS visual inspection (see 2017MNRAS.469.1065D, for further details) provides a visual confirmation of the SDSS pipeline redshift and object classification. The results of this inspection include a flag indicating the confidence of the redshift, “CONF_BEST”, which can take the values 3 (highly secure), 2 (uncertain), 1 (poor/unusable), 0 (insufficient data). A confirmation of the source classification was also added during the visual inspection, which uses the categories QSO, broad absorption line QSO (BALQSO), blazar, galaxy, star, and none. 2007AJ....133..313A provide the broad line AGN (BLAGN) and narrow line AGN (NLAGN) classifications, which are defined based on the presence or absence of broad (FWHM ¿ 1000) permitted emission lines.

The main goal of this work is to analyse the type 1 AGN in the SPIDERS sample, and therefore only sources that have been classified via their optical spectra as either “BLAGN” or “QSO” were selected for spectroscopic analysis. This returned 7805/9899 2RXS and 1192/1536 XMMSL1 sources. Since the categories “BLAGN” and “QSO” are based on different classification criteria, there will be some overlap between the two sets of sources. Therefore, no distinction will be made between the two categories; instead, both sets of objects will be considered type 1 AGN in this work.

2.7 Contamination from Starburst Galaxies

Although our sample probes luminosity ranges typically associated with AGN emission, starburst galaxies are also powerful X-ray sources and may be present as contaminants in our AGN sample. The X-ray emission from starburst galaxies is expected to originate from a number of energetic phenomena including supernova explosions and X-ray binaries (e.g. 2002A&A...382..843P). Therefore, the X-ray emission from starburst galaxies can be expected to be correlated with the star formation rate (SFR). Using their sample of luminous infrared galaxies, and a sample of nearby galaxies from 2003A&A...399...39R, 2011A&A...535A..93P found that the total SFR is related to the soft X-ray luminosity as follows:

(1)

2015A&A...579A...2I, figure 3, show the specific SFR for the COSMOS (2007ApJS..172....1S) and GOODS (2004ApJ...600L..93G) surveys for a series of redshift bins in the range . The peak of the redshift distribution of the 2RXS/XMMSL1 samples presented in this work is . Therefore, assuming that the COSMOS/GOODS sample in the redshift bin 0.2-0.4 is a good representative of the 2RXS/XMMSL1 samples, the upper limit on the SFR that can be expected for galaxies in our sample is . According to equation 1, this corresponds to a soft X-ray luminosity of , which is below the lower range probed by our samples , see figure 1).

2.8 Redshift Constraints

Using the “CONF_BEST” flag, sources with uncertain redshift or spectral classification (identified during the visual inspection of the sample) were also removed. This process resulted in a sample of 7795/7805 2RXS sources and 1190/1192 XMMSL1 sources. In the spectral fitting procedure (described in section 3), the H and MgII lines were fit independently. Sources with H and MgII present in their optical spectrum were selected using the following logic:

Different redshift ranges have been used because the BOSS spectrograph has a larger wavelength coverage than the SDSS spectrograph. In some cases, parts of the fitting region will have been redshifted out of the SDSS/BOSS spectrograph range (2013AJ....146...32S), and therefore will not be fit. However, the redshift limits where chosen so that both samples contain the broad lines used for estimating BH mass. Sources with a median signal-to-noise ratio (S/N) less than or equal to five per resolution element were excluded from the spectral analysis since for these sources the broad line decomposition and resulting BH mass estimates may be unreliable (see 2009ApJ...692..246D; 2011ApJS..194...45S).

Table 1 lists the numbers of sources with spectral coverage of either H or MgII, while figure 4 shows the redshift distribution of these sources. There are 711 cases where the same optical counterpart was detected by both 2RXS and XMMSL1. The final combined (2RXS and XMMSL1) sample for spectral analysis contains 7790 unique type 1 sources.

MgII H H and MgII H or MgII
2RXS 3310 6268 2234 7344
XMMSL1 314 1070 227 1157
Total 3473 6654 2337 7790
Table 1: The coverage of the H and MgII emission lines in the two samples used in this work. There are 711 sources which were detected in both the 2RXS and XMMSL1 surveys. The “total” row lists the total number of unique sources obtained from combining the 2RXS and XMMSL1 samples.
Figure 4: Redshift distribution of the sample of type 1 AGN with coverage of and/or MgII.

3 Spectral Analysis

A series of scripts have been written to perform spectral fits using the MPFIT least-squares curve fitting routine (2009ASPC..411..251M). Each spectrum was corrected for Milky Way extinction using the extinction curve from 1989ApJ...345..245C, and the dust map from 1998ApJ...500..525S, with an R=3.1. No attempt has been made to estimate and correct for the intrinsic (host) extinction of each source222 Also note that extinction laws (e.g. Calzetti and Prevot) are based on samples of nearby SB and irregular type galaxies. Due to the lack of nearby passive galaxies, an extinction law for these galaxy types is not yet available.. Measured line widths were corrected for the resolution of the  SDSS/BOSS spectrographs. The H and MgII emission line regions were fit independently using similar methods described in the following sections333For each model parameter, the 1-sigma uncertainties from MPFIT were adopted..

3.1 Iron Emission Template

AGN typically exhibit FeII emission consisting of a large number of individual lines across the optical and UV regions of the spectrum. These lines appear to be blended, probably due to the motion of the gas from which they are emitted, and the magnitude of this broadening varies significantly from source to source. The presence of FeII emission in the optical and UV portions of the spectrum can be a significant complication when attempting to accurately measure line profiles. Therefore it is crucial that the model used to derive line widths for BH mass measurements also accounts for the nearby FeII emission.

Figure 5 shows the two FeII templates used in this work; the 2001ApJS..134....1V and 1992ApJS...80..109B templates used for the UV and optical regions of the spectrum, respectively. Both of these templates have been derived from the narrow line Seyfert 1 galaxy I Zwicky 1 which, due to its bright FeII emission and narrow emission lines, is an ideal candidate for generating the FeII template. In order to model the observed blending of the FeII emission, the templates were convolved with a Gaussian whose width was included as a free parameter in the fitting procedure.

Figure 5: Upper panel: The 2001ApJS..134....1V FeII template which was used when fitting the MgII emission line region. Lower panel: The 1992ApJS...80..109B FeII template which was used when fitting the H emission line region. In both cases, the original template is shown (blue) along with the same template convolved with a Gaussian with FWHM = 4000 km s (red). The vertical dashed lines correspond to the position of MgII and H at 2800 and 4861.
Figure 6: Examples of model fits to the H (left panel; plate=1159, MJD=52669, fiber=470) and MgII (right panel; plate=423, MJD=51821, fiber=250) spectral regions. The model components are colour-coded as follows; power law (orange), iron emission (violet), broad lines (blue), narrow lines (yellow), [OIII] shifted wings (red), and the total model (red). The panels beneath the spectra show the data/model ratio.

3.2 H

The region from 4420-5500  was fit for each spectrum. The continuum model consisted of a power law, a galaxy template, and the 1992ApJS...80..109B FeII emission template. The FeII template was convolved with a Gaussian while fitting, and the width of this Gaussian, along with the normalisation of the template were included as free parameters in the fit (see section 3.1). Previous spectral analyses of AGN spectra have assumed an early-type galaxy component in the model (2017MNRAS.472.4051C). Following this method, we use an early-type SDSS galaxy template444template number 24 on
http://classic.sdss.org/dr5/algorithms/spectemplates/
in the fit, and the normalisation of this template as well as the normalisation and slope of the power law were also free parameters. The use of a single, early-type galaxy template is an approximation, however, it is considered to be justified since AGN are typically found to reside in bulge dominated galaxies (e.g. 2005ApJ...627L..97G; 2007ApJ...660L..19P), and the spectroscopic fiber collects emission mostly from the bulge (which is characterised by an old stellar population) and the active nucleus.

The [OIII]4959 and [OIII]5007 narrow lines were each fit with two Gaussians, one used to fit the narrow core, and an additional Gaussian to account for the presence of blue-shifted wings which are often detected (2005AJ....130..381B). A single Gaussian was used to fit the HeII4686 emission line. To avoid overfitting the H line, the fitting process was run four times, with one, two, three, and four555Three broad Gaussians are used in addition to a single narrow component to account for the three distinct broad components that are expected to be present (see section 8.1) in at least some sources (2010MNRAS.409.1033M). Gaussian components used to fit the H line. For each fit, the velocity width and peak wavelength of one of the Gaussian components was fixed to that of [OIII]4959 and [OIII]5007 in order to aid the identification of the narrow H component. The normalisation ratio of the [OIII]4959 and [OIII]5007 lines was fixed to the expected value of 1:3 (e.g. 2000MNRAS.312..813S). The best-fit model was then selected using the Bayesian information criterion (BIC, 1978AnSta...6..461S), which can be written as

where n is the number of data points, k is the number of model parameters, and is the chi-square of the fit. The preferred model is that with the lowest BIC. An example of a fit to the H spectral region is shown in the left panel of figure 6.

3.2.1 Broad Line Decomposition

The narrow H and [OIII] components are required to have widths . Any of the additional Gaussians used to fit MgII and H with FWHM  are considered “broad”. This threshold of 800 is taken from the approximate division between broad and narrow FWHM distributions in the lower panels of figure 12. The virial FWHM used for BH mass estimation is the FWHM of the line profile defined by the sum of these broad Gaussian components (see figure 7). A major challenge with using the single-epoch method for estimating BH mass is decomposing the broad and narrow components of the line in order to measure the virial FWHM. Figure 7 presents an example of the decomposition of a broad H line. In this case, the narrow H core can be easily distinguished and removed before measuring the virial FWHM. However, there are many cases where the broad and narrow components are blended, making it difficult to successfully identify the appropriate virial FWHM. There are also cases where there is a clear distinction between two broad line components that are shifted in wavelength relative to each other (known as “double-peaked emitters”). How one should interpret the single-epoch BH mass estimates for these unusual objects is uncertain (also see section 9.2).

Figure 7: An example of the decomposition of a typical AGN H line (plate=7276, MJD=57061, fiber=470). The horizontal dashed line represents the FWHM used for BH mass estimation. The vertical dashed line indicates the rest-frame wavelength of H. See section 3.2.1 for further details.

3.3 MgII

The region from 2450-3050 was fit for each spectrum. As in the case of the H fits, a power law, an early-type galaxy template (5 Gyr old elliptical galaxy; 1998ApJ...509..103S; 2007ApJ...663...81P), and the 2001ApJS..134....1V FeII emission template were used to fit the continuum. Again, the FeII template normalisation, and width of the Gaussian smoothing applied to the template, were included as free parameters in the fit. The MgII line is a doublet; however, due to the close spacing and virial broadening of the lines, it usually appears as a single broad component in AGN spectra. The narrow MgII line cores are usually not observed in AGN spectra, therefore the MgII profile was fit using three broad Gaussians. An example of a fit to the MgII spectral region is presented in the right panel of figure 6.

A B C Reference
MgII, L 1.816 0.584 1.712 2012ApJ...753..125S
H, L 0.91 0.5 2 2006ApJ...641..689V
H, L 0.895 0.52 2 2011ApJ...742...93A
Table 2: BH mass calibrations used in this work. A, B, and C are the calibration constants for single-epoch mass estimation (see equation 2).
Figure 8: Differences between the BH mass calibrations used in this work (see section 4 for further details). The vertical blue lines indicate the mean value of each distribution. The standard deviation, , of each distribution is also shown.

4 Bolometric Luminosity and BH Mass Estimation

Bolometric luminosities were estimated from the monochromatic luminosities using the bolometric corrections from 2006ApJS..166..470R; 2011ApJS..194...45S:

These bolometric corrections have been derived using mean AGN SEDs; however, 2006ApJS..166..470R note that using a bolometric correction resulting from a single mean SED can result in bolometric luminosities with inaccuracies up to 50.

Under the assumption that the BLR gas is virialised, the single-epoch method can be used to estimate BH mass as follows:

(2)

where is the monochromatic luminosity at wavelength , and FWHM is the full width at half maximum of the broad component of the emission line. A, B, and C are constants that are calibrated using RM results and vary depending on which line is used.

Over the years, many groups have provided calibrations of equation 2 for MgII and H. In this work, the calibrations from 2006ApJ...641..689V and 2011ApJ...742...93A are used for H. 2006ApJ...641..689V based their work on an updated study of the relationship (2005ApJ...629...61K; 2006ApJ...644..133B) and a reanalysis of the RM mass estimates (2004ApJ...613..682P) and therefore presented an improved mass calibration relative to previous studies. 2011ApJ...742...93A provide a mass calibration that is based on the relationship from 2009ApJ...705..199B. 2006ApJ...641..689V and 2011ApJ...742...93A both provide similar calibrations for single-epoch H mass estimation, as can be seen from the left panel of figure 8.

The 2012ApJ...753..125S calibration is used in this work for MgII. This calibration is based on a sample of 60 high-luminosity () quasars in the redshift range 1.5 - 2.2. 2012ApJ...753..125S use the 2006ApJ...641..689V mass estimates as a reference when determining their MgII calibration. The centre and right panels of figure 8 show the comparison between the 2012ApJ...753..125S MgII calibration and the calibrations from 2006ApJ...641..689V and 2011ApJ...742...93A. These calibrations agree reasonably well, with the standard deviation in both cases, which is likely due to the fact that these BH mass estimates were derived using two different emission lines. A list of the three BH mass calibrations used in this work is given in Table 2.

BH masses were computed for each of these calibrations and are included in the catalogue (see section A). BH masses were only estimated for sources with a detected broad line component (see section 3.2.1). These BH mass estimates were then used to estimate the Eddington luminosity and the Eddington ratio

where c is the speed of light, G is the gravitational constant, M is the BH mass, m is the proton mass, and is the Thomson scattering cross-section.

Figure 9: Comparison between the classical and Bayesian methods for estimating 2RXS fluxes. The deviation from a ratio of one at fainter fluxes results from the attempt to correct for the effect of the Eddington bias (see section 5.1 for further details).
Figure 10: Comparison between the results obtained from fitting the H and MgII regions of the spectrum. Upper left: 5100 and 2500 monochromatic luminosities. Upper right: Bolometric luminosities estimated from the 5100 and 3000 monochromatic luminosities. Lower left: Virial FWHM measured using H and MgII. Lower right: BH mass estimates derived from H (using the 2011ApJ...742...93A calibration) and MgII (using the 2012ApJ...753..125S calibration). In each panel, the solid black line is the unity line.

5 X-ray Flux Estimates

Since X-ray detections are available for all objects in this sample, X-ray flux estimates have also been included in the catalogue. XMMSL1 fluxes in the 0.2-12 keV range from 2008A&A...480..611S are included. 2008A&A...480..611S convert the XMMSL1 count rates to fluxes using a spectral model consisting of an absorbed power law with a photon index of 1.7 and N. The 2RXS fluxes were estimated using the method outlined below.

5.1 Estimating 2RXS X-Ray Fluxes

Many of the sources in the 2RXS sample have flux measurements close to the ROSAT flux limit (10). Therefore, when estimating fluxes for this sample, it was necessary to correct for the Eddington bias. This was done by adopting a Bayesian method to derive a probability distribution of fluxes based on the known distribution of AGN as a function of flux. Following 1991ApJ...374..344K, 2009ApJS..180..102L, and 2011MNRAS.414..992G, the probability of a source having flux f, given an observed number of counts C, is

(3)

where C is the total number of observed source and background counts, T is the mean expected total counts in the detection cell for a given flux, and is the prior, which is the distribution of AGN per unit X-ray flux interval. The exact expression for the prior was taken from 2008MNRAS.388.1205G, equation 1.

Source and background counts were taken from the 2RXS catalogue (2016A&A...588A.103B). A flux-count rate conversion factor, which was required to estimate T in equation 3, was derived using XSPEC (1996ASPC..101...17A) assuming a model consisting of a power law (with following 2017MNRAS.469.1065D) absorbed by the Milky Way column density. This method was used to estimate the flux in the full ROSAT band (0.1-2.4 keV) as well the monochromatic flux at 2 keV.

The fluxes resulting from the method described above with and without applying the prior (termed “Bayesian” and “classical”, respectively) are compared in Figure 9. The disagreement between the two flux estimates increases with decreasing flux, which is expected since, without the prior, the classical method fails to account for the Eddington bias. Low count rate sources in this sample would be assigned unrealistically low Bayesian fluxes. To avoid this, the flux was left as undetermined when the Bayesian flux estimate was more than a factor of ten smaller than the classical flux estimate.

6 Accessing the Data

The results from the spectral analysis discussed above, along with X-ray flux measurements and visual inspection results, have been made available in an SDSS DR14 value added catalogue which is available at http://www.sdss.org/dr14/data_access/value-added-catalogs/. Additionally, an extended version of the catalogue will be maintained at http://www.mpe.mpg.de/XraySurveys/SPIDERS/SPIDERS_AGN/ The column description for the catalogue is given in appendix A.

Figure 11: Bolometric luminosity versus redshift (left panel) and bolometric luminosity versus BH mass (right panel) for the sample presented in this work and the 2011ApJS..194...45S sample. Sources with -derived BH masses are shown in blue, and sources with MgII-derived BH masses are displayed in green.

7 Comparing the UV and Optical Fitting Results

A subsample of sources have spectral measurements available from both the MgII and H spectral regions. In order to test the consistency of the independent fits to these two regions, properties measured from both were compared.

7.1 Relation

A subsample of AGN whose spectra cover the rest-frame wavelengths 2500 and 5100 was selected using the following criteria:

Of this sample, 1718 sources had reliable measurements of both and . The relation was fit using the LINMIX (2007ApJ...665.1489K) package. LINMIX is a Bayesian linear regression algorithm that accounts for uncertainties in both dependent and independent variables, as well as non-detections. The upper left panel of figure 10 shows the distribution and the best-fit relation

(4)

with a regression intrinsic scatter of 0.0151. The comparison between the estimated bolometric luminosities derived from the 3000 and 5100 monochromatic fluxes is shown in the upper right panel of figure 10. Equation 4 can be used to estimate L from L, which allows low redshift sources to be included in the analysis discussed in section 8.3.

7.2 Comparing MgII and H FWHM Measurements

A subsample of AGN whose spectra cover the broad and MgII emission lines was selected using the following criteria

Of this sample, 2323 sources had FWHM measurements for both H and MgII. The lower left panel of figure 10 displays the virial FWHM measurements from H and MgII. The resulting best-fit relation, fit using LINMIX, is

with a regression intrinsic scatter of 0.005. This deviation from the one-to-one relation has also been observed by 2009ApJ...707.1334W, who reported a slope of , and 2012ApJ...753..125S, who found a slope of . The single-epoch BH mass relations (equation 2) account for the slope; when the correct BH mass calibration is used, the H and MgII virial FWHM measurements yield BH masses that are in close agreement (see the lower right panel of figure 10).

8 Sample Properties

Figure 11 presents the comparison between this sample and the full sample of SDSS DR7 AGN with optical spectral properties measured by 2011ApJS..194...45S in the bolometric luminosity-redshift and bolometric luminosity-BH mass planes. As discussed in section 1, -derived BH masses are used where available (shown in blue), while MgII-derived masses are used for the remaining higher-redshift sources (shown in green). The left panel of figure 11 shows that this sample populates the low-redshift, high-luminosity region of the parameter space, which is partially due to the high flux threshold of the X-ray selection. From the right panel of figure 11 it can be seen that the sample presented in this work appears to be well bounded by the Eddington limit at least up to .

8.1 H Line Components

Section 3.2 described how the H line profile was fit with either one, two, three, or four Gaussian components. Figure 12 displays the resulting distribution of FWHM measurements (the panels are split based on the number of Gaussian components required to fit the line). There is a clear peak in the distribution at low FWHM associated with the narrow H core typically measured at a few hundred . Above the distribution is bimodal (in the two lower panels) with a large number of sources showing evidence for the “very broad component” (VBC) of at FWHM also discussed in 2010MNRAS.409.1033M.

It has been suggested that the VBC is emitted from a distinct physical region, and is possibly the result of line emission from the accretion disk (e.g. 2009MNRAS.400..924B). If the VBC represents emission from the accretion disk, then a strong VBC may result in a bias towards a higher BH mass estimate, since the single-epoch method assumes a calibration that is based on the luminosity-BLR radius relation. However, since the kinematics and physical origin of the VBC remains uncertain, detected VBCs have not been excluded from the broad line profiles used to measure the virial FWHM in this analysis (as discussed in section 3.2.1).

Figure 12: Distribution of H Gaussian FWHM values. The panels are split based on the number of Gaussians required to fit the line. The coloured histograms each represent one of up to four possible Gaussians used to fit the H line. The grey histograms represent the sum of the individual coloured histograms.
Figure 13: FWHM of the broad component of H versus R. The left panel displays the sample described in section 8.2 (grey), and sources with a median spectral S/N20 (blue). The right panel presents the same subsample of high S/N sources, colour coded to indicate the trend in Eddington ratio across the distribution. The grey dashed line marks the division between population A () and population B () sources. The size of the typical uncertainties, multiplied by a factor of five, in both variables for the high-S/N subsample is also shown.

8.2 This Sample in the 4D Eigenvector 1 Context

The 4D Eigenvector 1 (4DE1) system (1992ApJS...80..109B; 2000ApJ...536L...5S; 2011BaltA..20..427S) aims to define a set of parameters that uniquely account for AGN diversity. Two main 4DE1 parameters are the FWHM of the broad component of H () and the strength of the FeII emission relative to that of H, defined as

where and are the fluxes of the FeII emission in the 4434-4684 range and broad line, respectively. A sample of 2098 sources with measurements of these parameters and reliable spectral fits () was selected. The left panel of figure 13 shows the distribution of this sample in the 4DE1 parameter space (grey). It is expected that a reliable measurement of the FeII component will be difficult for many of the lower S/N sources (see 2003ApJS..145..199M). For this reason, the subset of sources in figure 13 with a median S/N greater than or equal to 20 per resolution element is also shown (blue). The right panel of figure 13 presents the higher S/N sources, colour-coded as a function of Eddington ratio. The expected trend of increasing Eddington ratio towards smaller and larger R is observed for this sample of high-S/N sources. Typically, sources with both high and high are not observed. If these sources exist, they may be difficult to detect since strong FeII emission might conceal a faint broad component. The potential bias in the 4DE1 plane source distribution due to model limitations and spectral S/N is discussed in sections 9.1.2 and 9.1.3.

The grey dashed line in the right panel of figure 13 indicates the division between population A () and population B () sources in the 4DE1 context (see 2011BaltA..20..427S). Population A sources often possess Lorentzian broad line profiles, and it has been suggested that Gaussian fits to population A broad lines will result in an underestimation of the BH mass (see 2014AdSpR..54.1406S).

8.3 Relationship Between AGN X-ray and Optical Emission

Quasars exhibit a non-linear relationship between their X-ray and UV emission, usually represented by the parameter

where are the monochromatic luminosities and frequencies at 2 keV and , respectively (2003AJ....125.2876V; 2005AJ....130..387S; 2006AJ....131.2826S; 2007ApJ...665.1004J; 2008ApJS..176..355K; 2009ApJ...690..644G; 2009ApJS..183...17Y; 2010A&A...512A..34L). The parameter is considered to be a proxy for the relative contribution of the UV accretion disk emission and the X-ray emission from the surrounding corona to the total luminosity. In order to study this relationship, a sample of sources with measurements of the 2keV, 2500Å, and 5100Å  luminosities was selected. For lower redshift sources without spectral coverage of 2500Å, equation 4 was used to estimate the 2500Å  luminosity from the 5100Å  luminosity. Extended sources were removed in order to prevent additional scatter in the relationship due to the contribution of the host galaxy. This was done by requiring that the SDSS g band “stellarity”666 For a description of how cModelMag_g and psfMag_g are measured see https://www.sdss.org/dr12/algorithms/magnitudes/ (defined as S(g) = cModelMag_g- psfMag_g) lies between 0.1. This sample does not contain X-ray sources with more than one potential AllWISE counterpart and therefore avoids cases where the X-ray detection includes emission from more than one object. This selection process resulted in a sample of 4777 sources. Figure 14 shows the parameter versus the monochromatic luminosity at 2500. The relation was fit using LINMIX, which gave the following best-fit result

with a regression intrinsic scatter of 0.0034. This slope is consistent with previous results from the literature (e.g. 2008ApJS..176..355K).

Figure 14: versus UV luminosity. The dashed line is the best linear fit to the distribution. The size of the typical uncertainties in both variables is also shown.

9 Interpreting the Data and Limitations

In this section, the reliability and limitations of the sample will be discussed.

9.1 Measuring the FeII Emission

Distinguishing the FeII component from the continuum emission becomes more difficult when using low S/N spectra. In addition, for a given S/N, it may also be more difficult to detect FeII emission if the intrinsic broadening of the FeII lines is large, since broader, blended FeII emission lines are more likely to be fit by the model as continuum emission (see 2003ApJS..145..199M). Using simulated AGN spectra, 2003ApJS..145..199M estimate the minimum detectable optical FeII emission as a function of H width for different bins of S/N.

A poor fit to the FeII emission may affect the accuracy of the BH mass estimates, since FeII emission can influence measurements of both the broad line width (see section 9.1.1) and the continuum luminosity. FeII emission may also conceal a broad component thus biasing a source’s position in the 4DE1 plane (figure 13). These potential issues are tested in the following three sections.

9.1.1 Accuracy of the Broad Emission Line FWHM Measurements for Sources with FeII Continuum Emission

A poor fit to the FeII emission may affect the measurement of the broad emission line width. To quantify the magnitude of this effect, the fitting script (using four Gaussians to fit ) was run with and without the FeII template on a sample of randomly selected sources. The fit without the FeII template represents the most extreme case where the FeII emission is completely ignored by the model. Therefore, the change in line widths measured by these two models should be the upper limit on what can be expected for cases where the FeII fit is inadequate. Figure 15 shows that the line width dispersion induced by ignoring the presence of FeII emission is .

Figure 15: Comparison of the FWHM measurements derived using a model with and without an FeII template. The vertical dashed line shows the mean value of the distribution.

9.1.2 Model Limitations in Detecting Sources in the 4DE1 Plane

Sources with both large and large are typically not observed, however, this absence may be due to model limitations; at high the broad component may be concealed beneath the FeII emission, and therefore may not be detected. The experiment outlined in this section was carried out in order to determine whether the spectral fitting code used in this work would return accurate measurements for sources with high and values.

A parameter space defined by and was divided into a 1212 grid. 10 S/N bins between 5 and 50 (a representative range for the samples presented in this work) were selected for each point on the grid, and 10 spectra were simulated for each combination, resulting in 14400 simulated spectra. For the parameters that were fixed in this experiment, the interquartile mean of the best-fit values for the type 1 AGN in this sample were used. The line profile was modelled with one narrow and one broad Gaussian. The wavelength range was set to (as in section 3.2) and the logarithmic wavelength spacing777https://www.sdss.org/dr12/spectro/spectro_basics/ was set to be equal to that of SDSS spectra;

These spectra were fit using a version of the fitting script which used one narrow and one broad Gaussian component to fit . The minimum S/N required for the fitting script to return the correct and combinations was then determined. In order to consider an and combination detectable at a given S/N, at least 7/10 spectra were required to have best fit and values that agreed with the input values.

Figure 16 shows the detectable and combinations for each point on the grid, along with the minimum S/N required to detect that combination. It is clear from figure 16 that at the S/N levels available in this sample, a large region of the 4DE1 parameter space would not be detected.

Figure 16: Minimum S/N required to detect a range of and combinations. The blue grid indicates the range of the parameter space covered in the simulation described in section 9.1.2. Points on the grid that do not have a minimum S/N indicator represent and combinations that are not detectable even at the highest S/N used in this experiment. Sources detected at these and combinations are likely to be spurious (see figure 17).
Figure 17: Comparison between the measured and simulated and values for the highest and lowest S/N bins used in section 9.1.2. For clarity, each figure displays only one of the ten sets of spectra for that S/N. The values used to simulate the spectra (grey points) are connected to the corresponding best fit measurements (except for cases where either the broad or FeII components were not detected). For clarity, the figures do not show a small number of unphysically high measurements.

9.1.3 Bias in the and Distribution due to Model Limitations

Figure 17 shows the comparison between the simulated and measured and values for the highest and lowest S/N bins used in section 9.1.2. The left panel of figure 17 shows that at low S/N the results are clearly biased against high and values. At higher S/N (figure 17, right panel), the accuracy of the lower left quadrant measurements is significantly improved. However, even at S/N=45.5 (which is approximately the upper end of the S/N distribution of the samples presented in this work) the high - measurements deviate significantly from the corresponding “true” values. This may suggest that the L-shaped distribution of sources in the 4DE1 plane (e.g. figure 13) is at least in part due to model limitations.

9.2 Reliability of the Single-Epoch Method for Mass Estimation

Assuming that AGN broad emission lines are produced by gas whose motion is dominated by the gravitational potential of the central SMBH, the single-epoch method is expected to produce reliable mass estimates when compared to RM (see 2006ApJ...641..689V), with a systematic uncertainty of 0.3-0.4 dex. However, it is not clear how to measure the virial FWHM of lines that deviate from this norm.

The spectrum shown in the left panel of figure 18 is an example of a source which exhibits a double-peaked H line profile, where a clear inflection point is visible between two velocity-shifted broad line components. Double-peaked broad line profiles in AGN are expected to be the result of emission from the accretion disk (1988MNRAS.230..353P; 1989ApJ...339..742C; 1994ApJS...90....1E; 2003AJ....126.1720S; 2003ApJ...599..886E). 2007MNRAS.376.1335Z have found that single-epoch BH mass estimates obtained from double-peaked line profiles are significantly larger than BH mass estimates derived from stellar velocity dispersion measurements. 2007MNRAS.376.1335Z suggest that this discrepancy is the result of an overestimation of the BLR radius by the single-epoch mass calibrations for these objects. Therefore, the BH mass estimates provided in this work for sources which exhibit double-peaked broad emission lines should be treated with caution.

The right panel of figure 18 shows an example of narrow absorption in the UV portion of the spectrum caused by intervening absorbing material along the line of sight to the AGN. Sources identified during the visual inspection as having narrow absorption lines have been flagged in the catalogue (column 189; flag_abs). These sources were fit using the model described in section 3.3 with the absorption line regions masked. However, in many cases, the absorption features distort the broad MgII line, and therefore the resulting BH mass estimates may not be reliable.

Figure 18: Left panel: Example of a source with a double-peaked H line profile (plate=7512, MJD=56777, fiber=321). Right panel: Example of a source showing narrow UV absorption features which have been masked when fitting the model (plate=8188, MJD=57348, fiber=946). See section 9.2 for details.

10 Conclusions

This work presents a catalogue of spectral properties for all SPIDERS type 1 AGN up to SDSS DR14. Visual inspection results were used to select a reliable subsample for spectral analysis, and the spectral regions around H and MgII were fit with a multicomponent model. Using the single-epoch method, BH masses, bolometric luminosities, Eddington ratios, along with additional spectral parameters were measured. A catalogue containing these results has been made available as part of a set of SDSS DR14 value added catalogues. This catalogue also includes the results of a visual inspection of the sample, and is available at http://www.mpe.mpg.de/XraySurveys/SPIDERS/SPIDERS_AGN/.

11 Acknowledgements

DC has participated in the International Max Planck Research School on Astrophysics at the Ludwig Maximilians University Munich. DC also acknowledges financial support from the Max Planck Society. DC would also like to thank Riccardo Arcodia and Julien Wolf for many helpful discussions. The authors would like to thank Josephine Reisinger for her contribution to the visual inspection of sources in this sample. Finally, the authors would like to thank the anonymous referee for providing a thorough critique of the paper which greatly improved its content.

Funding for the Sloan Digital Sky Survey IV has been provided by the Alfred P. Sloan Foundation, the U.S. Department of Energy Office of Science, and the Participating Institutions. SDSS-IV acknowledges support and resources from the Center for High-Performance Computing at the University of Utah. The SDSS web site is www.sdss.org.

SDSS-IV is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS Collaboration including the Brazilian Participation Group, the Carnegie Institution for Science, Carnegie Mellon University, the Chilean Participation Group, the French Participation Group, Harvard-Smithsonian Center for Astrophysics, Instituto de Astrofísica de Canarias, The Johns Hopkins University, Kavli Institute for the Physics and Mathematics of the Universe (IPMU) / University of Tokyo, Lawrence Berkeley National Laboratory, Leibniz Institut für Astrophysik Potsdam (AIP), Max-Planck-Institut für Astronomie (MPIA Heidelberg), Max-Planck-Institut für Astrophysik (MPA Garching), Max-Planck-Institut für Extraterrestrische Physik (MPE), National Astronomical Observatories of China, New Mexico State University, New York University, University of Notre Dame, Observatário Nacional / MCTI, The Ohio State University, Pennsylvania State University, Shanghai Astronomical Observatory, United Kingdom Participation Group, Universidad Nacional Autónoma de México, University of Arizona, University of Colorado Boulder, University of Oxford, University of Portsmouth, University of Utah, University of Virginia, University of Washington, University of Wisconsin, Vanderbilt University, and Yale University.

Funding for SDSS-III has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, and the U.S. Department of Energy Office of Science. The SDSS-III web site is http://www.sdss3.org/.

SDSS-III is managed by the Astrophysical Research Consortium for the Participating Institutions of the SDSS-III Collaboration including the University of Arizona, the Brazilian Participation Group, Brookhaven National Laboratory, Carnegie Mellon University, University of Florida, the French Participation Group, the German Participation Group, Harvard University, the Instituto de Astrofísica de Canarias, the Michigan State/Notre Dame/JINA Participation Group, Johns Hopkins University, Lawrence Berkeley National Laboratory, Max Planck Institute for Astrophysics, Max Planck Institute for Extraterrestrial Physics, New Mexico State University, New York University, Ohio State University, Pennsylvania State University, University of Portsmouth, Princeton University, the Spanish Participation Group, University of Tokyo, University of Utah, Vanderbilt University, University of Virginia, University of Washington, and Yale University.

Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site is http://www.sdss.org/.

The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participating Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington.

This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, and NEOWISE, which is a project of the Jet Propulsion Laboratory/California Institute of Technology. WISE and NEOWISE are funded by the National Aeronautics and Space Administration.

Plot colours were in part based on www.ColorBrewer.org, by Cynthia A. Brewer, Penn State. The TOPCAT tool (2005ASPC..347...29T) was used during this work.

Appendix A Catalogue Column Description

The results from the spectral fitting code have been compiled into a single dataset published as an SDSS DR14 value added catalogue. A description of the columns included in the catalogue is given below. Catalogue entries that are either undetermined or not relevant for a given source are set to “-99”. Columns that provide information on the X-ray detections have been taken from 2008A&A...480..611S and 2016A&A...588A.103B.

  1. xray_detection: Flag indicating whether the X-ray source was detected in the 2RXS or XMMSL1 survey.

  2. name: IAU name of the X-ray detection.

  3. RA: Right ascension of the X-ray detection (J2000) [degrees].

  4. DEC: Declination of the X-ray detection (J2000) [degrees].

  5. ExiML: Existence likelihood of the 2RXS detection.

  6. ExpTime_2RXS: Exposure time (2RXS) [s].

  7. DET_ML: Detection likelihood of the XMMSL1 detection in the 0.2-12 keV range.

  8. ExpTime_XMMSL: Exposure time (XMMSL1) [s].

  9. f_2RXS: Flux in the 0.1-2.4 keV range (2RXS) [].

  10. errf_2RXS: Uncertainty in the flux in the 0.1-2.4 keV range (2RXS) [].

  11. fden_2RXS: Flux density at 2 keV (2RXS) [].

  12. errfden_2RXS: Uncertainty in the flux density at 2 keV (2RXS) [].

  13. l_2RXS: Luminosity in the 0.1-2.4 keV range (2RXS) [].

  14. errl_2RXS: Uncertainty in the luminosity in the 0.1-2.4 keV range (2RXS) [].

  15. l2keV_2RXS: Monochromatic luminosity at 2 keV (2RXS) [].

  16. errl2keV_2RXS: Uncertainty in the monochromatic luminosity at 2 keV (2RXS) [].

  17. f_XMMSL: Flux in the 0.2-12 keV range (XMMSL1; from 2008A&A...480..611S) [].

  18. errf_XMMSL: Uncertainty in the flux in the 0.2-12 keV range (XMMSL1; from 2008A&A...480..611S) [].

  19. l_XMMSL: Luminosity in the 0.2-12 keV range (XMMSL1; from 2008A&A...480..611S) [].

  20. errl_XMMSL: Uncertainty in the luminosity in the 0.2-12 keV range (XMMSL1; from 2008A&A...480..611S) [].

  21. Plate: SDSS plate number.

  22. MJD: MJD that the SDSS spectrum was taken.

  23. FiberID: SDSS fiber identification.

  24. DR14_RUN2D: Spectroscopic reprocessing number.

  25. DR14_PLUGRA: Right ascension of the drilled fiber position [degrees].

  26. DR14_PLUGDEC: Declination of the drilled fiber position [degrees].

  27. redshift: Source redshift based on the visual inspection results.

  28. CLASS_BEST: Source classification based on the visual inspection results.

  29. CONF_BEST: Visual inspection redshift and classification confidence flag.

  30. DR14_ZWARNING: Warning flag for SDSS spectra.

  31. DR14_SNMEDIANALL: Median S/N ratio per pixel of the spectrum.

  32. Instrument: Flag indicating which spectrograph was used (SDSS or BOSS) to measure the spectrum.

  33. norm1_mgII: Normalisation of the first Gaussian used to fit the MgII line [].

  34. errnorm1_mgII: Uncertainty in the normalisation of the first Gaussian used to fit the MgII line [].

  35. peak1_mgII: Wavelength of the peak of the first Gaussian used to fit the MgII line [].

  36. errpeak1_mgII: Uncertainty in the wavelength of the peak of the first Gaussian used to fit the MgII line [].

  37. width1_mgII: Width of the first Gaussian used to fit the MgII line [Å].

  38. errwidth1_mgII: Uncertainty in the width of the first Gaussian used to fit the MgII line [Å].

  39. fwhm1_mgII: FWHM of the first Gaussian used to fit the MgII line [].

  40. errfwhm1_mgII: Uncertainty in the FWHM of the first Gaussian used to fit the MgII line [].

  41. shift1_mgII: Wavelength shift of the peak of the first Gaussian used to fit the MgII line relative to the rest-frame wavelength [].

  42. norm2_mgII: Normalisation of the second Gaussian used to fit the MgII line [].

  43. errnorm2_mgII: Uncertainty in the normalisation of the second Gaussian used to fit the MgII line [].

  44. peak2_mgII: Wavelength of the peak of the second Gaussian used to fit the MgII line [].

  45. errpeak2_mgII: Uncertainty in the wavelength of the peak of the second Gaussian used to fit the MgII line [].

  46. width2_mgII: Width of the second Gaussian used to fit the MgII line [Å].

  47. errwidth2_mgII: Uncertainty in the width of the second Gaussian used to fit the MgII line [Å].

  48. fwhm2_mgII: FWHM of the second Gaussian used to fit the MgII line [].

  49. errfwhm2_mgII: Uncertainty in the FWHM of the second Gaussian used to fit the MgII line [].

  50. shift2_mgII: Wavelength shift of the peak of the second Gaussian used to fit the MgII line relative to the rest-frame wavelength [].

  51. norm3_mgII: Normalisation of the third Gaussian used to fit the MgII line [].

  52. errnorm3_mgII: Uncertainty in the normalisation of the third Gaussian used to fit the MgII line [].

  53. peak3_mgII: Wavelength of the peak of the third Gaussian used to fit the MgII line [].

  54. errpeak3_mgII: Uncertainty in the wavelength of the peak of the third Gaussian used to fit the MgII line [].

  55. width3_mgII: Width of the third Gaussian used to fit the MgII line [Å].

  56. errwidth3_mgII: Uncertainty in the width of the third Gaussian used to fit the MgII line [Å].

  57. fwhm3_mgII: FWHM of the third Gaussian used to fit the MgII line [].

  58. errfwhm3_mgII: Uncertainty in the FWHM of the third Gaussian used to fit the MgII line [].

  59. shift3_mgII: Wavelength shift of the peak of the third Gaussian used to fit the MgII line relative to the rest-frame wavelength [].

  60. norm_heII: Normalisation of the Gaussian used to fit the HeII line [].

  61. errnorm_heII: Uncertainty in the normalisation of the Gaussian used to fit the HeII line [].

  62. peak_heII: Wavelength of the peak of the Gaussian used to fit the HeII line [].

  63. errpeak_heII: Uncertainty in the wavelength of the peak of the Gaussian used to fit the HeII line [].

  64. width_heII: Width of the Gaussian used to fit the HeII line [Å].

  65. errwidth_heII: Uncertainty in the width of the Gaussian used to fit the HeII line [Å].

  66. fwhm_heII: FWHM of the Gaussian used to fit the HeII line [].

  67. errfwhm_heII: Uncertainty in the FWHM of the Gaussian used to fit the HeII line [].

  68. shift_heII: Wavelength shift of the peak of the Gaussian used to fit the HeII line relative to the rest-frame wavelength [].

  69. norm1_hb: Normalisation of the first Gaussian used to fit the line [].

  70. errnorm1_hb: Uncertainty in the normalisation of the first Gaussian used to fit the line [].

  71. peak1_hb: Wavelength of the peak of the first Gaussian used to fit the line [].

  72. errpeak1_hb: Uncertainty in the wavelength of the peak of the first Gaussian used to fit the line [].

  73. width1_hb: Width of the first Gaussian used to fit the line [Å].

  74. errwidth1_hb: Uncertainty in the width of the first Gaussian used to fit the line [Å].

  75. fwhm1_hb: FWHM of the first Gaussian used to fit the line [].

  76. errfwhm1_hb: Uncertainty in the FWHM of the first Gaussian used to fit the line [].

  77. shift1_hb: Wavelength shift of the peak of the first Gaussian used to fit the line relative to the rest-frame wavelength [].

  78. norm2_hb: Normalisation of the second Gaussian used to fit the line [].

  79. errnorm2_hb: Uncertainty in the normalisation of the second Gaussian used to fit the line [].

  80. peak2_hb: Wavelength of the peak of the second Gaussian used to fit the line [].

  81. errpeak2_hb: Uncertainty in the wavelength of the peak of the second Gaussian used to fit the line [].

  82. width2_hb: Width of the second Gaussian used to fit the line [Å].

  83. errwidth2_hb: Uncertainty in the width of the second Gaussian used to fit the line [Å].

  84. fwhm2_hb: FWHM of the second Gaussian used to fit the line [].

  85. errfwhm2_hb: Uncertainty in the FWHM of the second Gaussian used to fit the line [].

  86. shift2_hb: Wavelength shift of the peak of the second Gaussian used to fit the line relative to the rest-frame wavelength [].

  87. norm3_hb: Normalisation of the third Gaussian used to fit the line [].

  88. errnorm3_hb: Uncertainty in the normalisation of the third Gaussian used to fit the line [].

  89. peak3_hb: Wavelength of the peak of the third Gaussian used to fit the line [].

  90. errpeak3_hb: Uncertainty in the wavelength of the peak of the third Gaussian used to fit the line [].

  91. width3_hb: Width of the third Gaussian used to fit the line [Å].

  92. errwidth3_hb: Uncertainty in the width of the third Gaussian used to fit the line [Å].

  93. fwhm3_hb: FWHM of the third Gaussian used to fit the line [].

  94. errfwhm3_hb: Uncertainty in the FWHM of the third Gaussian used to fit the line [].

  95. shift3_hb: Wavelength shift of the peak of the third Gaussian used to fit the line relative to the rest-frame wavelength [].

  96. norm4_hb: Normalisation of the fourth Gaussian used to fit the line [].

  97. errnorm4_hb: Uncertainty in the normalisation of the fourth Gaussian used to fit the line [].

  98. peak4_hb: Wavelength of the peak of the fourth Gaussian used to fit the line [].

  99. errpeak4_hb: Uncertainty in the wavelength of the peak of the fourth Gaussian used to fit the line [].

  100. width4_hb: Width of the fourth Gaussian used to fit the line [Å].

  101. errwidth4_hb: Uncertainty in the width of the fourth Gaussian used to fit the line [Å].

  102. fwhm4_hb: FWHM of the fourth Gaussian used to fit the line [].

  103. errfwhm4_hb: Uncertainty in the FWHM of the fourth Gaussian used to fit the line [].

  104. shift4_hb: Wavelength shift of the peak of the fourth Gaussian used to fit the line relative to the rest-frame wavelength [].

  105. norm1_OIII4959: Normalisation of first Gaussian used to fit the line [].

  106. errnorm1_OIII4959: Uncertainty in the normalisation of first Gaussian used to fit the line [].

  107. peak1_OIII4959: Wavelength of the peak of the first Gaussian used to fit the line [].

  108. errpeak1_OIII4959: Uncertainty in the wavelength of the peak of the first Gaussian used to fit the line [].

  109. width1_OIII4959: Width of the first Gaussian used to fit the line [Å].

  110. errwidth1_OIII4959: Uncertainty in the width of the first Gaussian used to fit the line [Å].

  111. fwhm1_OIII4959: FWHM of the first Gaussian used to fit the line [].

  112. errfwhm1_OIII4959: Uncertainty in the FWHM of the first Gaussian used to fit the line [].

  113. shift1_OIII4959: Wavelength shift of the peak of the first Gaussian used to fit the line relative to the rest-frame wavelength [].

  114. norm2_OIII4959: Normalisation of second Gaussian used to fit the line [].

  115. errnorm2_OIII4959: Uncertainty in the normalisation of second Gaussian used to fit the line [].

  116. peak2_OIII4959: Wavelength of the peak of the second Gaussian used to fit the line [].

  117. errpeak2_OIII4959: Uncertainty in the wavelength of the peak of the second Gaussian used to fit the line [].

  118. width2_OIII4959: Width of the second Gaussian used to fit the line [Å].

  119. errwidth2_OIII4959: Uncertainty in the width of the second Gaussian used to fit the line [Å].

  120. fwhm2_OIII4959: FWHM of the second Gaussian used to fit the line [].

  121. errfwhm2_OIII4959: Uncertainty in the FWHM of the second Gaussian used to fit the line [].

  122. shift2_OIII4959: Wavelength shift of the peak of the second Gaussian used to fit the line relative to the rest-frame wavelength [].

  123. norm1_OIII5007: Normalisation of first Gaussian used to fit the line [].

  124. errnorm1_OIII5007: Uncertainty in the normalisation of first Gaussian used to fit the line [].

  125. peak1_OIII5007: Wavelength of the peak of the first Gaussian used to fit the line [].

  126. errpeak1_OIII5007: Uncertainty in the wavelength of the peak of the first Gaussian used to fit the line [].

  127. width1_OIII5007: Width of the first Gaussian used to fit the line [Å].

  128. errwidth1_OIII5007: Uncertainty in the width of the first Gaussian used to fit the line [Å].

  129. fwhm1_OIII5007: FWHM of the first Gaussian used to fit the line [].

  130. errfwhm1_OIII5007: Uncertainty in the FWHM of the first Gaussian used to fit the line [].

  131. shift1_OIII5007: Wavelength shift of the peak of the first Gaussian used to fit the line relative to the rest-frame wavelength [].

  132. norm2_OIII5007: Normalisation of second Gaussian used to fit the line [].

  133. errnorm2_OIII5007: Uncertainty in the normalisation of second Gaussian used to fit the line [].

  134. peak2_OIII5007: Wavelength of the peak of the second Gaussian used to fit the line [].

  135. errpeak2_OIII5007: Uncertainty in the wavelength of the peak of the second Gaussian used to fit the line [].

  136. width2_OIII5007: Width of the second Gaussian used to fit the line [Å].

  137. errwidth2_OIII5007: Uncertainty in the width of the second Gaussian used to fit the line [Å].

  138. fwhm2_OIII5007: FWHM of the second Gaussian used to fit the line [].

  139. errfwhm2_OIII5007: Uncertainty in the FWHM of the second Gaussian used to fit the line [].

  140. shift2_OIII5007: Wavelength shift of the peak of the second Gaussian used to fit the line relative to the rest-frame wavelength [].

  141. norm_pl1: Normalisation of the power law fit to the MgII continuum region [].

  142. errnorm_pl1: Uncertainty in the normalisation of the power law fit to the MgII continuum region [].

  143. slope_pl1: Slope of the power law fit to the MgII continuum region.

  144. errslope_pl1: Uncertainty in the slope of the power law fit to the MgII continuum region.

  145. norm_pl2: Normalisation of the power law fit to the H continuum region [].

  146. errnorm_pl2: Uncertainty in the normalisation of the power law fit to the H continuum region [].

  147. slope_pl2: Slope of the power law fit to the H continuum region.

  148. errslope_pl2: Uncertainty in the slope of the power law fit to the H continuum region.

  149. norm_gal1: Normalisation of the galaxy template used to fit the MgII continuum region [].

  150. errnorm_gal1: Uncertainty in the normalisation of the galaxy template used to fit the MgII continuum region [].

  151. norm_gal2: Normalisation of the galaxy template used to fit the H continuum region [].

  152. errnorm_gal2: Uncertainty in the normalisation of the galaxy template used to fit the H continuum region [].

  153. norm_feII1: Normalisation of the iron template used to fit the MgII continuum region [].

  154. errnorm_feII1: Uncertainty in the normalisation of the iron template used to fit the MgII continuum region [].

  155. norm_feII2: Normalisation of the iron template used to fit the H continuum region [].

  156. errnorm_feII2: Uncertainty in the normalisation of the iron template used to fit the H continuum region [].

  157. fwhm_feII1: FWHM of the Gaussian kernel convolved with the iron template used to fit the MgII continuum region [].

  158. errfwhm_feII1: Uncertainty in the FWHM of the Gaussian kernel convolved with the iron template used to fit the MgII continuum region [].

  159. fwhm_feII2: FWHM of the Gaussian kernel convolved with the iron template used to fit the H continuum region [].

  160. errfwhm_feII2: Uncertainty in the FWHM of the Gaussian kernel convolved with the iron template used to fit the H continuum region [].

  161. r_feII: Flux ratio of the 4434-4684 FeII emission to the broad component of H.

  162. OIII_Hbeta_ratio: Flux ratio of [OIII]5007 to H.

  163. virialfwhm_mgII: FWHM of the MgII broad line profile [].

  164. errvirialfwhm_mgII: Uncertainty in the FWHM of the MgII broad line profile [].

  165. virialfwhm_hb: FWHM of the broad line profile [].

  166. errvirialfwhm_hb: Uncertainty in the FWHM of the broad line profile [].

  167. mgII_chi: Reduced chi-squared () of the fit to the MgII region.

  168. hb_chi: Reduced chi-squared () of the fit to the region.

  169. l_2500: Monochromatic luminosity at 2500 [].

  170. errl_2500: Uncertainty in the monochromatic luminosity at 2500 [].

  171. l_3000: Monochromatic luminosity at 3000 [].

  172. errl_3000: Uncertainty in the monochromatic luminosity at 3000 [].

  173. l_5100: Monochromatic luminosity at 5100 [].

  174. errl_5100: Uncertainty in the monochromatic luminosity at 5100 [].

  175. l_bol1: Bolometric luminosity derived from the monochromatic luminosity at 3000 [].

  176. errl_bol1: Uncertainty in the bolometric luminosity derived from the monochromatic luminosity at 3000 [].

  177. l_bol2: Bolometric luminosity derived from the monochromatic luminosity at 5100 [].

  178. errl_bol2: Uncertainty in the bolometric luminosity derived from the monochromatic luminosity at 5100 [].

  179. logBHMVP_hb: BH mass derived from the H line using the 2006ApJ...641..689V calibration [].

  180. errlogBHMVP_hb: Uncertainty in the BH mass derived from the H line using the 2006ApJ...641..689V calibration [].

  181. logBHMA_hb: BH mass derived from the H line using the 2011ApJ...742...93A calibration [].

  182. errlogBHMA_hb: Uncertainty in the BH mass derived from the H line using the 2011ApJ...742...93A calibration [].

  183. logBHMS_mgII: BH mass derived from the MgII line using the 2012ApJ...753..125S calibration [].

  184. errlogBHMS_mgII: Uncertainty in the BH mass derived from the MgII line using the 2012ApJ...753..125S calibration [].

  185. l_edd1: Eddington luminosity based on the BH mass estimate derived using the 2012ApJ...753..125S calibration [].

  186. l_edd2: Eddington luminosity based on the BH mass estimate derived using the 2011ApJ...742...93A calibration [].

  187. edd_ratio1: Eddington ratio defined as l_bol1/l_edd1.

  188. edd_ratio2: Eddington ratio defined as l_bol2/l_edd2.

  189. flag_abs: Flag indicating whether or not strong absorption lines have been observed in the spectrum. flag_abs is set to either 0 (no absorption present) or 1 (absorption present).

References

Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
363442
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description