2 THE VCV OPTICAL CATALOGUE

A fundamental part of constructing the HMQ was to audit and incorporate the best previous comprehensive catalogue of quasars, which is of course the Véron-Cetty & Véron Catalogue of Quasars and Active Nuclei 13th edition (Véron-Cetty & Véron Reference Véron-Cetty and Véron2010). However, VCV included astrometry and identifications which were not always optically based, for legacy reasons. In Flesch (Reference Flesch2013), I promulgated the need to convert that to a fully optical presentation to enable easy inclusion into dynamic databases like NEDFootnote ⁴ and SIMBADFootnote ⁵, and in service of that I presented the true position or status of 449 VCV objects which hitherto had no obvious optical counterpart. I have continued that task to where the VCV catalogue is now fully mapped, line for line, into the optical J2000. This ‘VCV Optical catalogue’Footnote ⁶ identifies and resolves all anomalies including those reported in Flesch (Reference Flesch2013). Following is a summary of those issues, and lists of objects supplementary to those presented in Flesch (Reference Flesch2013).

VCV presented 167 567 QSOs/AGN of which 160 271 were subclassed as type 1 (or had no subclass), 6 203 as type 2, 926 as LINER galaxies, and 167 as starburst galaxies; also 1 374 BL Lac objects were listed, thus 168 941 in total. The HMQ optical mapping, with the benefit of SDSS data through DR12, resolves these as 156 322 type 1, 1 142 BL Lac objects, 4 662 type 2, 747 LINERs, 5 024 galaxies, 267 stars, 667 objects of unknown classification (usually from SDSS-DR7 or later SDSS releases which dropped them), 8 satellite streaks, 2 asteroids (i.e., 3-colour strobes on SDSS), 2 star spikes, 61 duplicates of other VCV objects, and 37 objects which are lost or never existed. It behoves me to explain these reclassifications, so first I list out the most involved discovery papers (the individual objects involved are listed in the VCV Optical catalogue with these citations):

1. (A) 4 113 objects were re-classified as stars or non-AGN galaxies by papers subsequent to VCV, usually SDSS releases:

   1. 1. 3 424 by SDSS-DR8 Aihara et al. (Reference Aihara2011) due to re-processing of the spectra via an improved algorithm.

   2. 2. 291 by SDSS-DR12Q Pâris et al. (Reference Pâris2015) which reclassified many objects via visual inspection of the spectra.

   3. 3. 267 by SDSS-DR9Q Pâris et al. (Reference Pâris2012) which reclassified many objects via visual inspection of the spectra.

   4. 4. 42 by SDSS-DR10Q Pâris et al. (Reference Pâris2014) which reclassified many objects via visual inspection of the spectra.

   5. 5. 39 by SDSS-DR9 Ahn et al. (Reference Ahn2012) via an improved automated pipeline.

   6. 6. 25 by SDSS-DR10 Ahn et al. (Reference Ahn2014) via an improved automated pipeline.

   7. 7. 14 by Kleinman et al. (Reference Kleinman2013) which identified white dwarfs from the SDSS-DR7 data.

   8. 8. 5 by Ge et al. (Reference Ge, Hu, Wang, Bai and Zhang2012) which identified double-peaked narrow emission line galaxies (NELGs).

   9. 9. 6 other objects from 4 papers.

2. (B) 1 845 objects were originally classified as stars or non-AGN galaxies by their discovery papers but taken up by VCV as ‘active galaxies’, or were reclassified as non-AGN prior to VCV’s final edition but VCV retained the original classification:

   1. 1. 473 from Schneider, Schmidt, & Gunn (Reference Schneider, Schmidt and Gunn1994) which classified these as emission-line galaxies (ELGs) of z < 0.45.

   2. 2. 430 from SDSS-DR7 Abazajian et al. (Reference Abazajian2009) which classified these as stars or non-AGN galaxies, confirmed by their exclusion from SDSS-DR7Q Schneider et al. (Reference Schneider2010). Many of these had been classified as AGN by earlier publications.

   3. 3. 261 from the 2dF Galaxy Redshift Survey (2dF-GRS: Colless et al. Reference Colless2001) which is a large galaxy catalogue without QSO classifications. VCV extracted high-redshift objects and others with emission lines showing activity, but visual inspection of the high-redshift objects reveal that they are often confused with closely co-positioned neighbours. Some of these objects were cited elsewhere and included by VCV on that basis, but subsequent SDSS classifications are as galaxies only. The optical appearance of almost all of these objects is that of galaxies, and only 16 show radio/X-ray association. I see no avenue by which to retain these objects as type I AGN, so they are dropped.

   4. 4. 144 from Iovino, Clowes, & Shaver (Reference Iovino, Clowes and Shaver1996) were annotated by them with ‘?’ signifying spectral overlap or possible star. Many of these are anomalously bright and none show radio/X-ray association, so they are all dropped. Some of the remaining 1 140 objects from that paper are quite bright or reddish, so caution in their use is indicated.

   5. 5. 88 from the 2QZ Croom et al. (Reference Croom, Smith, Boyle, Shanks, Miller, Outram and Loaring2004), which classified 85 of these as ‘continuum’ objects and 3 as NELGs. None show radio/X-ray association and all of the continuum objects were marked by VCV as ‘BL?’ showing that they were questionable, so they are dropped.

   6. 6. 64 from Mauch & Sadler (Reference Mauch and Sadler2007) who classified these as LINERs.

   7. 7. 47 from Schneider, Schmidt, & Gunn et al. (Reference Schneider, Schmidt and Gunn1999) who, as with their 1994 paper, classified these as ELGs of z < 0.45.

   8. 8. 45 from La Franca et al. (Reference La Franca, Lissandrini, Cristiani, Miller, Hawkins and McGillivray1999) which classified these as ELGs.

   9. 9. 25 from Barcons et al. (Reference Barcons2007) which classified these as NELGs.

   10. 10. 18 from Sandage & Bedke (Reference Sandage and Bedke1994) which classified these NGC galaxies as LINERs and starburst.

   11. 11. 10 from Gaston (Reference Gaston1983a) and 2 from Gaston (Reference Gaston1983b), both of which surveyed SA57 (i.e., ‘selected area 57’), but these 12 are without quasar colours or radio/X-ray association. Gaston (Reference Gaston1983b), sec III, qualified that the SA57 photometry had ‘much greater scatter’ and ‘nonuniformity’ compared with concomitant CFHT and KPNO photometry, and in sec IIai stated a forthcoming paper, ‘Hoag et al. 1983’, of these SA57 quasars which however never appeared. Edwards, Beauchemin, & Borra (Reference Edwards, Beauchemin and Borra1988) resurveyed SA57, confirming no Gaston object, and related that Hoag in fact ‘did not confirm many of the Gaston candidates’.

   12. 12. 238 other objects from 117 papers, usually from VCV’s own annotations that they are questionable, or starburst galaxies, etc.

A total of 61 duplications are found in VCV, of which 39 were presented in Flesch (Reference Flesch2013) and one in Cupani et al. (Reference Cupani, Cristiani, D’Odorico, Milvang-Jensen and Krogager2011). The remaining 21 newly-found duplications are listed in Table 2 which displays the VCV names of the wrongly-sited duplicate and correctly-sited master, the offset between them in arcseconds, and the final HMQ name which is usually the original form. One of these, #10 on the list, deserves a special mention in that VCV took the redshift of the master from SDSS-DR6, but that redshift was corrected by SDSS-DR7 to that given by the original discoverers.

Table 2. 21 more duplicates in VCV13.

^aVCV Table 1 listed QSOs, Table 3 listed AGN.

As a further supplement to Flesch (Reference Flesch2013), Table 3 lists 31 more moves of 30+ arcsec for VCV objects. 29 of these are newly presented, whilst two, #s 4 and 5, are moved from the position given in Flesch (Reference Flesch2013) due to better information. Also, four objects were wrongly moved in Flesch (Reference Flesch2013), being Q 0004-019, IRAS 12202+1646, and MS 22225+2114 which are returned to their VCV location, and IXO 32 which is moved to its finding chart location at 09 10 20.2 +07 05 55 (J2000). The right-hand column is a comment on the move, ‘OA’=original author.

Table 3. 31 more moves of 30+ arcsec in VCV13.

^aVCV Table 1 listed QSOs, Table 3 listed AGN.

Also, Flesch (Reference Flesch2013) identified 30 VCV objects to de-list because they either are not quasars or have information of such poor quality that the object is lost beyond any mechanism of recovery. In a similar vein, I present 14 VCV objects in Table 4 which are not seen at all, and the methods of Flesch (Reference Flesch2013) find no candidate. These objects are obviously not included in the HMQ catalogue.

Table 4. 14 more objects to delete from VCV13.

^aVCV Table 1 listed QSOs, Table 3 listed AGN.

^bVmag from VCV, can be blank.

HMQ has sourced much data directly from large source catalogues like SDSS, 2QZ Croom et al. (Reference Croom, Smith, Boyle, Shanks, Miller, Outram and Loaring2004), etc., so the actual uptake from VCV consists of a residue from smaller papers which amount to 19 127 type 1 objects, 1 083 BL Lac objects, and 1 594 type 2 objects for a total of 21 794 legacy objects. Most of the citations for these legacy objects have been copied from the VCV citation list.

3 INCLUSION OF SDSS QUASARS

The majority of quasars in this catalogue are, of course, sourced from the ongoing SDSS releases which have dwarfed in size and complexity all previous quasar surveys; a full description of the SDSS project is beyond the scope of this paper, so the reader is assumed to have some acquaintance with it. There are two simplifications in the HMQ citations: the SDSS data releases 1 and 2 are subsumed into the DR3 release, consistent with their presentation in the SDSS data archives, and DR11 which was an internal release only is similarly subsumed into DR12 and so not considered further.

In use of the SDSS, it is not enough to simply take up all SDSS ‘quasars’ because they are not as uniformly reliable as that. The SDSS project uses an automated pipeline which matches spectra to standard templates by which any spectrum well-matched to a quasar template is securely classified as a quasar. But of course many spectra are ambiguous fits to templates, especially for faint objects with noisy spectra. To cater for this, the SDSS consortium periodically issues manually-checked quasar catalogues which are much more reliable than the automated ones, the most recent being SDSS-DR12Q Pâris et al. (Reference Pâris2015). I use those manual catalogues as a definitive source although they do have a few false objects to remove—indeed, each edition typically drops some objects from previous editions. Objects not accepted into those manual catalogues (i.e., ‘pipeline-only’ objects) may not be quasars or their spectra could be too low quality (i.e.,‘low significance’) to be credible, but often they simply fell on the wrong side of threshold cutoffs such as magnitude or luminosity. Inspection of the spectra of such pipeline-only quasars shows that indeed many are plainly recognisable as quasar spectra, and many are not. Therefore, I needed to find criteria to bulk-separate the high significance pipeline-only quasars from the low, in order to accept the good ones into this comprehensive catalogue. Fortunately this could be done as follows.

As a background, the SDSS data releases 1 to 7 comprise the ‘legacy’ period in which the classifications of quasars, galaxies, stars, and ‘unknowns’ were released as separate incremental tables, and the manual catalogues were quasars only. There was no means of bulk-grading the pipeline-only quasars in the legacy data. The current ‘SDSS3’ data releases 8 to 12 are made as single tables which include all classifications and historical data, and the manual catalogues include a ‘superset’ of all manually inspected classifications. Key enhancements of the SDSS3 releases are the provision of a ZWARNING field which flags problems with the spectra, and a SUBCLASS field which displays values like ‘BROADLINE’ or ‘STARFORMING’, but is often blank. For the ZWARNING field, I follow documented recommendations and accept only objects with ZWARNING of 0 (no warnings) or 16 (saturation) as well as 256 (poor astrometry for plughole) and some 128 (unplugged, but flag is documented to be wrong for about 200 objects which have strong signal). For the SUBCLASS field, I find from inspection of hundreds of spectra that the better quality spectra bear a populated SUBCLASS field, and the pipeline classifications and redshifts for those objects are of a good significance. In contrast, the spectra without a subclass are usually noisy and the pipeline information of a low significance—which undoubtedly explains the pipeline’s inability to discern a subclass. Accordingly, I take a populated SUBCLASS field to signal a high significance spectrum.

Thus, by use of the SUBCLASS and ZWARNING fields I am able to bulk-separate the pipeline-only quasars into high significance and low significance spectra, thus enabling inclusion of those high significance quasars into this catalogue. An additional benefit is that there is a high overlap of the SDSS3 data with the legacy SDSS data; thus all legacy pipeline-only quasars can be dropped unless qualifying by this SDSS3-based rule. However, pipeline-only quasars are also kept if confidently associated to radio/X-ray detections, even if the spectrum is of low significance. Also, pipeline-only redshifts above z = 4 are excluded as per advice from Adam Myers; upon further inspection of high-z pipeline-only spectra, I have lowered this cutoff to z = 3.65.

The total counts of SDSS-discovered type-1 objects included into HMQ are 369 249 manually checked quasars and 20 298 pipeline-only quasars. Also, 11 821 SDSS redshifts have been applied to quasars from earlier papers, the SDSS redshifts being favoured due to their precision and reliability. Specific breakdowns are available from the HMQ data in which each object is referenced with its discovery citation and redshift citation. Table 5 lists the discovery papers which contribute most quasars to the HMQ catalogue, in numerate order; there it is seen that the SDSS releases dominate the quasar portfolio today. A note on the Table 5 column ‘# of type 1’: this shows how many objects were originally presented by that paper, but SDSS astrometrically-based object names can change slightly between releases; the HMQ presents the current names of SDSS objects, so the discovery citation is that of the SDSS release which first presented the current name. My processing showed that DR12Q added 121 173 new type-1 quasars, so the 154 958 displayed in Table 5 includes 33 785 quasars with small name changes from earlier releases.

Table 5. Top 25 discovery papers for the HMQ.

^aDR3Q and DR3 here include the earlier DR1 and DR2 releases.

^bThe Principal Galaxy Catalogue (PGC) is not actually a discovery paper, but is used as a reference for names of AGN galaxies.

4 NEW QUASARS

114 new quasars with spectral confirmation are presented in this catalogue. They are 92 quasars from Cyril Hazard’s work in the 1980s, 4 quasars from E. Margaret Burbidge (Reference Burbidge2003), and 18 quasars from Dan Weedman on a re-inspection of Weedman (Reference Weedman1985). In detail:

1. 1. 92 quasars from Cyril Hazard, the prominent early quasar researcher who identified many quasars from objective prism plates obtained with the UK Schmidt telescope in the 1980s Hazard (Reference Hazard1980). Hazard’s modus operandi was to share his quasars out with other researchers rather than publish them himself. As a consequence, many remained unpublished. In June 2013, he sent word that he is preparing a comprehensive presentation of his works, and we shared some data. These 92 new quasars were identified from that data, and Cyril graciously assented to their inclusion in this work. The quasar names are from his worksheets which I have prefaced with ‘HAZ’ to clarify their provenance; however, two objects had the same worksheet identifier of ‘122’, so I have renamed one as simply ‘HAZARD’. 24 of these objects have been covered by the subsequent SDSS DR12 release with all confirmed as quasars—these are readily identified in the catalogue by the redshift citation. A sample is shown in Table 6 and the complete set can be extracted from the catalogue by selecting on names beginning with ‘HAZ’, or by the citation of Hazard (Reference Hazard2013). Full information on these quasars and a larger set of work is in preparation by Cyril Hazard as of his last communication in September 2013.

   Table 6. Sample quasars from Cyril Hazard.

   

2. 2. Four quasars from E. Margaret Burbidge, observed on the Keck-I LRIS on 2003 October 2 . Six targets were selected via X-ray association in the vicinity of NGC 891, and five found to be quasars; the sixth (J022302.5+421731) was missed because the slit was mislaid onto two adjacent stars. Burbidge elected not to publish due to the sparse results, but released the material to associates, myself included. The five quasars are listed in Table 7 with X-ray associations. The spectra were reduced by C. Gutiérrez using IRAFFootnote ⁷ and are displayed in Figure 2 which I have annotated with line identifications. The first quasar has since been re-surveyed by the Serendipitous Extragalactic X-ray Source Identification program (SEXSI: Eckart et al. Reference Eckart, Stern, Helfand, Harrison, Mao and Yost2006), so the remaining four are newly published here. I have named them as EMBUR 2 (etc.) after the discoverer, and their citation is Burbidge (Reference Burbidge2003).

   

   Figure 2. Spectra of the five Burbidge quasars, obtained with Keck-I LRIS using a 5600Å dichroic.

3. 3. 18 quasars with probable redshifts from Dan Weedman. These were first presented in Weedman (Reference Weedman1985) but without redshifts because each showed only a single emission line which could not be securely identified as Lyα or CIII, etc. Today we have SDSS-based photometric redshifts for these 18, primarily from Richards et al. (Reference Richards2009), which have enabled unambiguous identifications of each single emission line on the premise that the photometric redshift is near to the true redshift. These quasars have names of WEE 35, WEE 49, etc., from the original paper, and have a redshift citation of D. W. Weedman (2014, private communication) except for three re-surveyed by DR12 which confirm these as quasars.

Table 7. Quasars from E. Margaret Burbidge (Reference Burbidge2003).

5 CLASSIFICATION OF THE PRESENTED OBJECTS

There are five object classes reported in this catalogue which are shown with their counts in the first two columns of Table 8. Type 1 objects (class=‘Q’ and ‘A’) are the main thrust of this catalogue and only high confidence objects with secure redshifts are included; where the discovery paper flags objects as questionable, those are dropped. SDSS quasars are included when presented by the manual surveys such as DR12Q, and also good confidence pipeline-only quasars are accepted as detailed in Section 3. However, the SDSS pipeline uses a ‘BROADLINE’ subclassifier for broadline emitters, so emitters not flagged as ‘BROADLINE’ are either narrow-line emitters or the spectrum was too low confidence to make a judgement. Therefore, to include pipeline-only quasars as type-1 objects, I require them to be either of good confidence with the ‘BROADLINE’ subclass, or of low confidence (i.e., blank subclass) with secure radio/X-ray association.

Table 8. Counts of Radio/X-ray associations in the HMQ.

Type 2 objects (class=‘K’) are included for user reference at ≈ 90% completeness from the literature, and are generally taken as narrow-line emitters although legacy usage varied. Also added are a selection of objects from the large SDSS releases which did not classify type 2 objects but included emitters not subclassed as ‘BROADLINE’, many of which will be type-2 objects. I used rules to include them as follows: firstly, 1 264 good confidence pipeline-only quasars not flagged as ‘BROADLINE’ were taken. Also, 1 089 manual quasars with those same attributes were taken, because the manual catalogues included but did not classify type-2 objects, see Pâris et al. (Reference Pâris2012), introduction, last paragraph. Next, 16 842 SDSS galaxies subclassed as ‘AGN’ but not ‘BROADLINE’ were taken as narrow line Seyferts although some will be LINERs; similarly, 3 247 SDSS galaxies subclassed as ‘AGN BROADLINE’ were taken up as AGN (class=‘A’). These 20 089 galaxy sourced objects are flagged in the HMQ comment field with an ‘a’ to signal their provenance. In this way, the germane active objects from the SDSS releases have been included, although the HMQ class=‘K’ consequently contains unquantified contamination from broadline quasars and LINERs.

BL Lac objects (class=‘B’) have been quality checked and ones without radio/X-ray association which VCV annotated as ‘BL?’ have been dropped. Newly surveyed such objects are usually presented as ‘candidates’, so I have added those where confirmed by radio/Xray association. There is an online compendium of blazars, the BZCAT Massaro et al. (Reference Massaro, Giommi, Leto, Marchegiani, Maselli, Perri, Piranomonte and Sclavi2009) which has been used as a check and a reference, but not as a definitive guide.

The AGN (class=‘A’) are in principle broadline Seyfert galaxies which are optically dominated by their disk, with QSOs (class=‘Q’) being core-dominated objects with no disk seen. Of course the actual data shows intermediate objects throughout the full range, in line with the unified model. No attempt is made to subclass the AGN as S1.2, S1.5, etc. VCV used a cosmology-based luminosity threshold to separate the AGN and QSOs, and re-released their final edition with a changed cosmology. By contrast, the approach here is simply a practical one to separate the visually extended AGN from the unresolved QSOs, analyse those groups to get a luminosity based divider, and extend that divider into the set of fainter unresolved objects. If done correctly, the separator should trace a path from the bright QSO-AGN divider to some natural break in the bulk fainter data profile as is seen in Figure 3. With some testing I have identified a simple luminosity equation

\begin{equation*} magnitude = 22 + 5 \times \text{log}_{10}( z \times (1 + z )), \end{equation*}

which traces the required path. However, type 1 objects of z < 0.1 (as a rough cut) are all taken as AGN Seyferts, consistent with the absence of nearby quasars. Thus, I need to adjust the separator to terminate at z = 0.1. Doing so gives the following:

\begin{equation*} {magnitude = 22 + 5 \times \text{log}_{10}( z \times (0.8 + z ) - 0.09)}, \end{equation*}

as the separator which is bounded by Vmag = 17.0 at z = 0.1. However, such bright objects can be classified by their psf. The separator is needed for the set of objects which are psf-unresolved just because they are so faint; these are SDSS objects with g ⩾ 22.0 (‘survey-quality threshold’, Pâris et al. Reference Pâris2014, Section 2.3) and Vmag ⩾ 21.0 for other objects. The separator is also needed for brighter objects which have no psf class, e.g., USNO-B photometry. This luminosity-based separator equation is that of a static Euclidean manifold with time dilation, so is non-physical in a scale-invariant universe, but fits the data well. Figure 3 shows the separator at x=0 as a local minimum in the data profile of the faint objects. The QSOs to the left (brighter) of the separator comprise the bulk of the objects; the twin peaks are an artefact of SDSS BOSSFootnote ⁸ selection. The AGN to the right (fainter) of the separator have a Gaussian-looking profile which, if real, implies some underlying physical distinction between the QSOs and AGN.

Figure 3. Vmag offsets from luminosity separator for faint (see text) or no-psf type-1 objects, binned by 0.1 mag. Chart produced with TOPCAT (Taylor Reference Taylor2005).

In summary of type 1 objects, extended psf morphology signals an AGN, with SDSS morphology used for SDSS photometry. Galaxies from the Principal Galaxy Catalogue Paturel et al. (Reference Paturel, Petit, Prugniel, Theureau, Rousseau, Brouty, Dubois and Cambrésy2003) are taken as extended, as well as any type-1 of z < 0.1. Brighter unresolved objects (as defined above) are classified as QSO, and fainter unresolved objects, and those without psf class, are classified as QSO or AGN using the above separator equation; the overall outcome is ≈ 94% QSOs.

6 OPTICAL PHOTOMETRY

The HMQ sources optical photometry from a background data pool as is documented in Flesch & Hardcastle (Reference Flesch and Hardcastle2004), but some points are worth highlighting. The optical data is presented as red and blue psf and magnitudes, and are preferentially sourced from Cambridge Automatic Plate Measuring machine (APM) McMahon & Irwin (Reference McMahon, Irwin, MacGillivray and Thomson1992) and United States Naval Observatory Monet et al. (Reference Monet2003) USNO-B data. In particular, their data from the first-epoch National Geographic-Palomar Observatory Sky Survey (POSS-I) are preferred because POSS-I O is violet 4 100Å and so well-offset from POSS-I E 6 500Å and both plates were always taken on the same night, thus the red–blue colour is correct even for variable objects. Other surveys’ photometry are used where POSS-I data is incomplete, and SDSS photometry is used where both red and blue magnitudes are not already available from the APM or USNO-B. The HMQ ReadMe has additional information about the photometry and how it is flagged in the data.

7 RADIO AND X-RAY ASSOCIATIONS

HMQ displays radio/X-ray associations just by giving the radio/X-ray identifier(s) without elaboration. Beneath this simple presentation are the calculations detailed in Flesch & Hardcastle (Reference Flesch and Hardcastle2004) which quantify causality, i.e., the confidence that the association is true. In brief, the method is that areal densities of optical profiles around radio/X-ray detections are compared to background averages to obtain the causality percentage. These calculations are done against anonymous optical data, and catalogues like the Atlas of Radio/X-ray Associations Flesch (Reference Flesch2010) which report associations to anonymous optical objects do need to display those calculated probabilities. But when that anonymous object turns out to be a quasar, which are commonly radio/X-ray emitters, then that is a strong confirmation that the association is correct. Thus, since HMQ has quasars only, this very simple presentation style can be used.

The radio/X-ray associations were useful in making many identifications, e.g., the quasar ‘Q 1224-116’ is identified as the r = 18.3 b = 18.6 object at 12 27 25.0-11 33 39 (J2000) because of the X-ray association 3XMM J122725.1-113340, whereas previously it was identified to the nearby galaxy 2MASSi J1227247-113334; the original author did not distinguish between them. This theme was encountered in many places.

The ReadMe itemises the radio/X-ray surveys used, but of course they are all the largest ones like the NRAO VLA Sky Survey catalog (NVSS: Condon et al. Reference Condon, Cotton, Greisen, Yin, Perley, Taylor and Broderick1998), the Faint Images of the Radio Sky at Twenty-cm survey catalog (FIRST: White et al. Reference White, Becker, Helfand and Gregg1997), and the Sydney University Molonglo Sky Survey Murphy et al. (Reference Murphy, Mauch, Green, Hunstead, Piestrzynska, Kels and Sztajer2007) for the radio data, and Chandra, XMM-Newton, Swift, and ROSAT catalogues for the X-ray data. Optical field solutions are calculated from the raw source positions of all these catalogues as described in Flesch (Reference Flesch2010). Radio and X-ray associations are calculated independently of each other.

The radio/X-ray identifiers are displayed in four columns which are one for X-ray, one for core radio, and two for double radio lobes. Some association identifications are seen in the sample data of Table 1 and Table 8 shows the total HMQ count of association types by object class. However, radio lobes are uncommon ( < 1%), so when there are no lobes those two columns are instead used to display additional radio/X-ray identifiers from other surveys; this conveys a fuller sense of the data available. The ReadMe gives the rules used. The four columns are usually enough to display all available associations, but any overflow is not shown; this happens for 2 891 objects with a total of 4 241 X-ray associations dropped, representing 5.7% of all X-ray associations. An example is the blazar ‘ON 231’, seen in Table 1 with 4 associations but has additional X-ray associations 2RXP J122132.0+281403, 1WGA J1221.5+2813 and 1RXS J122131.1+281402. Such extra identifiers would be nice to show but it is not essential because they all refer to the same physical X-ray source seen on the sky, so are duplicates in that sense.

8 INCLUSION OF PHOTOMETRIC QUASARS

This catalogue includes 25 015 high-confidence (99%) SDSS photometric quasar candidates which also show radio and/or X-ray associations. These objects are well assured of being true quasars because two independent procedures, the photometry calculations and the radio/X-ray association calculations, separately identify the same objects as likely quasars, and the calibrated joint calculation given below quantifies that to the 99% likelihood required.

These objects are within the scope of this catalogue because the confidence in the spectroscopically confirmed quasars, i.e., that they are true quasars, are also rated at about 99%. This assessment is based on historical performance, and three examples convey the point. First, the VCV quasar catalogue, heretofore the authoritative standard, was found above to have ≈ 5 000 galaxies or stars amongst its 168 941 nominal active objects for an error rate of 3% and thus performance of 97%. Those have been removed but an unquantified residue remains from legacy surveys. The second example is that of the 20 298 SDSS pipeline-only quasars, documented above, which are included because the preponderance of them are true quasars, but some noisy spectra are present which conceal an error rate which I can estimate only informally at about 5%. Thirdly, even the high-quality SDSS manual surveys have a typical drop rate of 0.15% per each new edition (Schneider et al. Reference Schneider2010, Section 6), and in the case of DR10Q which presented 169 459 QSOs, the subsequent DR12Q release reclassified 1 655 of those into galaxies and dropped 25 entirely, for a total 0.99% drop rate. So, false positives slip through even under state-of-the-art conditions. Thus, over the set of all classified quasars in the literature, my assessment of an overall 1% error rate is a fair call.

Accordingly, the confirmed objects can be supplemented with photometric candidates which as a group also have 99% confidence of being true quasars. These photometric quasars are from the SDSS-based NBCKDE (i.e., Nonparametric Bayes classifier – kernel density estimate) Richards et al. (Reference Richards2009) and XDQSO (i.e., Extreme-deconvolution QSO targets) Bovy et al. (Reference Bovy2011) photometric quasar candidate catalogues. These catalogues give nominal estimates of pQSO (probability of being a true QSO) for each object which however are not calibrated against confirmed classified objects. I have performed that calibration against DR12Q classifications, and the details of that calibration and discussions of the NBCKDE and XDQSO attributes are in Appendix A. This yields a calibrated photometrically-based QSO likelihood for each candidate.

Next, QSO likelihoods are needed for the radio/X-ray associations. The calculations described in Section 7 obtain the confidence of association percentage, i.e., causality, and the classified objects used in that analysis include a percentage of QSOs; thus, causality × QSO% = QSO likelihood. This likelihood is already calibrated against confirmed objects as part of the calculation. Double radio lobes are included via a heuristic algorithm, also calibrated.

Now the calibrated photometric based QSO likelihood P₁ and the radio/X-ray association based QSO likelihood P₂ are combined as (expressing the likelihoods as fractions)

\begin{equation*} {P_{\text{tot}} = 1 / (1 + ((1 - P_{1})(1 - P_{2}) / P_{1}P_{2} ))}, \end{equation*}

This equation requires a component to be above 50% for it to make a positive contribution to the outcome; therefore we consider only objects with both P₁ > 50% and P₂ > 50%. The goal is to identify the set of candidates with P _tot ≥ 99%. Table 9 shows the result, binned by rounded P_tot, and the three right-hand columns show the results for fellow candidates which have since been classified by SDSS-DR12Q (with ZWARNING=0 for best precision). Column 6 gives the actual QSO percentage of classified objects, and shows a performance break between the rows binned at 96% and 97%, so I provisionally accept the bins of 97% to 100% (i.e., the bottom 4 rows) as equivalent actual performers, and sum them at bottom. For that total set, the column 3 calculated likelihood is 98.977% which is acceptably close to the 99% criterion, and the column 6 actual QSO percentage of 99.175% meets the criterion. In this way, the 25 015 are selected.

Table 9. QSO likelihoods for the SDSS candidates, with confirmed percentages.

These 25 015 photometric objects are given class=‘q’ in this catalogue to identify them as objects without full spectra, but photometric redshifts are displayed. These are sourced as available from NBCKDE which provides redshifts for 14 202 NBCKDE-only and shared objects. The remaining 10 813 objects are XDQSO-only, but XDQSO has not published photometric redshifts. These I have calculated via a clustering method akin to rainflow analysis which is described in Appendix B. Note that photometric redshift calculations sometimes fail to achieve a result, for example when two contending redshift values are adjudged equally likely. This was the case for an additional 627 otherwise-qualifying candidates, representing a failure rate of 2.4%, which were accordingly excluded from the HMQ because redshifts are required. All photometric redshifts are rounded to 0.1 z in the catalogue to help convey that they are estimates. Sample objects are seen at top and bottom rows of Table 1.

9 MISCELLANEOUS NOTES

The HMQ includes 6 770 new spectroscopically confirmed QSOs (and 4 560 type-2) from 101 papers published since the VCV 13th edition, apart from the SDSS surveys. My thanks to those authors who sent me lists of their objects. Also, 180 QSOs from 67 early papers overlooked by VCV have also been added.

The biggest ongoing task throughout this project was de-duplication of objects across papers and surveys, while preserving true doublets. These are all resolved except maybe for two BL Lac objects (or is it one) which are SHBL J23435+3440 at J234333.8+344005 and BZBJ2343+3439 at J234333.5+343948. These objects are 14 arcsec apart, and both have redshift=0.336. The SHBL object has radio NVSS J234333.7+344000 and the BZB object has X-ray 1SXPS J234333.4+343949. The BZB object has white BL Lac colours while the SHBL object is quite red, but its discovery paper was clear on the identification. I suspect the SHBL object is the true one and the BZB object a duplicate, but possibly both are valid so both are in this catalogue.

The HMQ is a catalogue of quasars and not of gravitational lens images, so I have de-duplicated multiple lensed signatures into single identifications where I was aware of them; the optically brightest signature is the one retained.

I write NVSS radio detections in the form NVSS Jhhmmss.s+ddmmss which is the standard form used by radio/X-ray surveys, although the recommended NVSS naming convention is NVSS Jhhmmss+ddmmss. I prefer the standard form because when both FIRST and NVSS report the same signature, their astrometry agrees well, often to the arcsecond, so the standard form displays useful astrometric information. The recommended NVSS name can be recovered by simply removing the time decimal.