A population of Optically Quiescent Quasars from WISE and SDSS

The growth of active galactic nuclei (AGN) occurs under some form of obscuration in a large fraction of the population. The difficulty in constraining this population leads to high uncertainties in cosmic X-ray background and galaxy evolution models. Using an SDSS-WISE cross-match, we target infrared luminous AGN ($W1-W2$>0.8, and monochromatic rest-frame luminosity above $\lambda L_{\lambda}$(12$\mu$m) $\approx$ 3 $\times$ 10$^{44}$ erg s$^{-1}$), but with passive galaxy-like optical spectra (Optically Quiescent Quasars; OQQs). We find 47 objects that show no significant [O III]$\lambda$5007 emission, a typically strong AGN optical emission line. As a comparison sample, we examine SDSS-selected Type 2 quasars (QSO2s), which show a significant [O III]$\lambda$5007 line by definition. We find a 1:16 ratio of OQQs compared to QSO2s, suggesting that the OQQ duty cycle is likely much shorter than that of QSO2s (though selection biases are not fully quantified). We consider observed properties in comparison with other galaxy types, and examine them for consistency with theories on their intrinsic nature: chiefly (a) a high covering factor for surrounding obscuring matter, preventing the detection of high-ionisation emission lines - `cocooned AGN'; or (b) ionised gas being absent on the kpc scales of the Narrow Line Region (NLR), perhaps due to a `switching on' or `young' AGN. OQQs do not obviously fit the standard paradigm for merger-driven AGN and host galaxy evolution, implying we may be missing part of the flow of AGN evolution.


INTRODUCTION
The census of active galactic nuclei (AGN) is, at present, highly incomplete.Dusty gas that feeds supermassive black hole growth can obscure the nucleus, resulting in an attenuation of AGN signatures along the line-of-sight (l.o.s.).This is not a minor effect: the majority of AGN are affected by obscuration.There are many lines of evidence pointing to this.Studies of the cosmic X-ray background radiation require obscuration by neutral gas with column densities ( H ) exceeding 10 22 cm −2 with between 2-5 as many times obscured AGN as unobscured ones (e.g.Setti & Woltjer 1989;Fabian & Iwasawa 1999;Gandhi & Fabian 2003;Gilli et al. 2007;Treister et al. 2009;Ueda et al. 2014;Ananna et al. 2019).In the optical, obscuration can successfully explain the variety of observed AGN classes.Optical Type 1 AGN show broad emission ★ E-mail: claire.l.greenwell@durham.ac.uk lines with widths of several thousand km s −1 arising on scales of ∼ 30  0.7 5100 light-days, where  5100 is the continuum rest-frame luminosity at 5100 Å in units of 10 44 erg s −1 (e.g.Kaspi et al. 2000).Observing these close nuclear scales requires an extinction-free l.o.s.. Type 1 and Type 2 AGN show narrower emission lines with widths of a few hundred km s −1 , which arise on scales of tens to thousands of pc.This difference in appearance can naturally be explained by an anisotropic distribution of dust which obscures the close-in broad lines, but not the larger scale narrower ones.In this way, the zoo of AGN classes can be unified.This model can also explain observed polarization fractions of broad emission lines (e.g.Antonucci & Miller 1985;Peterson 1997).In the radio, the radioloud 1 subset of obscured quasars ('radio galaxies') was the first significant population of powerful and heavily absorbed AGN to be followed up in detail (e.g.McCarthy 1993;Miley & De Breuck 2008;Toba et al. 2019); many X-ray studies have now shown them to be strongly obscured, on average (e.g.Gandhi et al. 2006;Tozzi et al. 2009;Wilkes et al. 2013).
The column density and geometry of obscuring matter are expected to naturally evolve as AGN and their host galaxies grow, and many models posit that the bulk of supermassive black hole growth occurs in highly obscured phases (e.g.Fabian 1999;Di Matteo et al. 2005;Hopkins et al. 2006).AGN are also known to appear to change classification over time, from Type 1 to 2 and vice versa (Yang et al. 2018).These objects, referred to as 'Changing Look' AGN, show changes in emission line and continuum flux over timescales up to a few years.The physical mechanisms behind these changes are not well understood: the two main theories are variation in the line-ofsight obscuration (e.g. a clumpy torus; Elitzur 2012), a change in the accretion rate (e.g.Sheng et al. 2017), or thermal changes in the inner accretion disk (e.g.Stern et al. 2018).Any unusual AGN appearance must be considered in this context -lack of emission lines could be a transitional state of changing obscuration, or a change in intrinsic line production.Objects in the process of 'switching on' are a rare, brief chapter in the growth of AGN, an important but not well understood period.Sources with extreme covering factors approaching unity could also probe a unique phase in AGN evolution, indicating either strong growth rates with plenty of available circumnuclear matter for accretion or perhaps sky covering as a result of merger-driven turbulence.Several studies have linked galaxy mergers with higher rates of obscured AGN in MIR selected samples (e.g.Glikman et al. 2012;Weston et al. 2017;Satyapal et al. 2017).Although Seyfert galaxies with high covering factors approaching unity have been inferred in detailed individual studies (Ramos Almeida et al. 2009), the fraction of highly covered AGN at high power is expected to be small (e.g.Toba et al. 2014;Stalevski et al. 2016, although the dependence on luminosity is weak and not without counter-evidence -e.g.Netzer et al. 2016).Theoretically, accretion from large scales is not expected to be isotropic and is likely to be mediated via disks, warps and other instabilities (e.g.Hönig 2019).
The infrared regime is particularly effective for studies of AGN dust covering factors, and for studies of AGN with absorbed optical signatures.This is because dust serves as a bolometer, absorbing and reprocessing the AGN power to the mid-infrared (MIR), providing a probe of the obscuring material in emission, as opposed to the absorption pathway provided in optical and X-ray studies.High angular resolution multiwavelength studies suggest that the MIR emission is effectively (to within a factor of a few) isotropic (Gandhi et al. 2009;Levenson et al. 2009;Asmus et al. 2015;Stalevski et al. 2016).Covering fractions derived from MIR AGN number counts and modelling of individual spectral energy distributions may also be a function of luminosity (Maiolino et al. 2007;Alonso-Herrero et al. 2011;Assef et al. 2013;Toba et al. 2013Toba et al. , 2014;;Toba et al. 2021), though this remains controversial (Roseboom et al. 2013;Lawrence & Elvis 2010;Assef et al. 2015).
One limitation of most works on covering fraction to-date is the requirement for the presence of AGN emission lines (typically forbidden lines from the Narrow Line Region; NLR) in the optical or near-infrared, used to confirm the presence of an AGN and/or redshift identification.If AGN emission lines are observed, this implies that some intrinsic power must be escaping the AGN environment and the source cannot be fully covered.The presence of dust in AGN NLRs has been known for some time (e.g., Netzer & Laor 1993;Haas et al. 2005;Netzer et al. 2006) and can attenuate forbidden line fluxes, but selecting fully covered AGN is difficult, given the obvious observational biases against identifying such a population.Thus, there have only been a limited number of studies of AGN in the high covering factor regime (e.g., Gandhi et al. 2002;Imanishi et al. 2007).The [O iii]5007 line is one of the strongest optical AGN emission lines and can be used as a proxy for bolometric luminosity (e.g., Heckman et al. 2004).It arises in the NLR, likely as a result of photoionization by AGN radiation.This places the line origin beyond the putative classical torus of unification models, so it should not be affected by l.o.s.nuclear reddening (although there is some effect, e.g. on the bolometric correction; Lamastra et al. 2009).Narrow line emission, and in particular [O iii], is assumed to be present in the majority of AGN and is often used for selection; e.g., BPT diagrams (Kewley et al. 2006) and Type 2 quasars (Zakamska et al. 2003).Sources lacking this line are unusual among AGN catalogues, and by selecting based on this absence we can target objects either (a) with unstable large scale obscuration preventing transmission of [O iii], or (b) in which the physical conditions result in no line formation in the first place; for example if the AGN is recently switched on, and the ionising radiation has not yet reached NLR scales.
Recent large surveys now enable studies searching for elusive AGN subtypes to be carried out.Here, we go beyond previous works in selecting candidate mid-infrared (MIR) quasars that show no optical signatures.For this, we use the latest all-sky MIR survey by the Wide-field Infrared Survey Explorer (WISE; Wright et al. 2010) mission, combined with optical spectroscopy from the Sloan Digital Sky Survey (SDSS) to select MIR-luminous sources with AGNlike MIR colours but with early-type galaxy optical spectra.In this way, we avoid star-forming galaxies and associated degeneracies in separating such systems from AGN (e.g.Trouille & Barger 2010).At lower luminosity, spectral dilution of AGN emission by the host galaxy cannot be neglected (e.g.Moran et al. 2002;Comastri et al. 2002;Cocchia et al. 2007;Civano et al. 2007;Caccianiga et al. 2007).At high luminosity, it becomes increasingly difficult for the host galaxy to dilute the AGN, so MIR quasar studies offer a clean probe of nuclear activity.
The structure of this paper is as follows.We describe the data used and our sample selection in Sections 2 and 3. Results follow in Section 4. A detailed discussion of implications, caveats, and comparisons with other source classes can be found in Section 6.An Appendix includes information on individual sources, the optical spectra and broadband SEDs.We assume a flat cosmology with  0 = 67.4km s −1 Mpc −1 and Ω Λ = 0.685 (Planck Collaboration et al. 2018).

WISE
The WISE satellite has carried out a highly-sensitive all-sky survey in four bands (1, 2, 3 and 4, centred on wavelengths of 3.4, 4.6, 12 and 22 µm, respectively).The effective angular resolution corresponds to a Gaussian with full-width-at-half-maximum 6 ′′ in 1 − 3, and 12 ′′ in 4.The SDSS fifteenth Data Release (Aguado 2019) includes pre-calculated astrometric cross matches with the earlier WISE data release (WISE-AllSky 2 ).The AllWISE release 3 includes data from both the cryogenic and post-cryogenic phases of the mission, and therefore contains better quality data, with better photometric sensivitity, particularly in the 1 and 2 bands.The cross matching with SDSS and WISE-AllSky was used, and the WISE magnitudes updated with AllWISE values.
We choose to use the AllWISE data rather than the more recent CatWISE (Marocco et al. 2021) because the 3 and 4 measurements are vital to the selection process, and these two bands were only in operation in the earlier cryogenic years of WISE.Although CatWISE may have provided more accurate photometry, (a) we are only looking at bright sources, which should be present in the earlier catalogs, and (b) in case of variation in source output it is best to use measurements from all four bands taken in the same time period.

SDSS
Data from the SDSS data release 15 (DR15; Aguado 2019), including Baryon Oscillation Spectroscopic Survey (BOSS) spectra (Dawson et al. 2013), were used for cross-matching with the WISE catalog.
The relevant tables are: • SpecObjAll: the base SDSS table for spectroscopic observations, containing all measured spectra, including potentially bad or duplicate data.Contains 4,851,200 sources.
• SpecObj: the sub table of SpecObjAll, containing only the clean data, filtered for duplicates.Contains 4,311,571 sources.
• PhotoObjAll: the base SDSS table for all photometric observations.Contains 1,231,051,050 measurements.
• PhotoObj: the subset of PhotoObjAll containing only primary and secondary objects (i.e.not a family object, outside the chunk boundary, or within a hole4 ).Contains 794,328,715 measurements.
• wise_allsky: the WISE catalog from the WISE-AllSky data release.Contains 563,921,584 sources.
• wise_xmatch: a joined table that contains pointers to SDSS and WISE-AllSky measurements, from astrometric cross-matches between the two.Contains 495,003,196 cross-matches.

SAMPLE SELECTION
The initial SDSS to WISE crossmatch (table wise_xmatch) uses a 4 arcsecond distance threshold.The nominal WISE positional accuracy is actually significantly better than this5 .This is also seen in the distribution of associated counterpart distances shown in Fig. 1, which peaks at <0.2 ′′ .A few selection tests showed that an increasing fraction of ambiguous counterpart identification (e.g. two potential SDSS sources for a single WISE source) for distances of more than 1 arcsecond.In fact, the vast majority (86.4%) of closest associations lie at distances of less than 1 ′′ , so this is the threshold we selected to build our catalog.We have chosen to do cross-matching based on positional separation only rather than using a more statistical method (e.g., Nway, Salvato et al. 2018) due to the low incidence of multiple potential counterparts.The probability of each match being correct is high, based on low WISE positional uncertainties.
The initial set was based on the wise_xmatch catalog, with every spectral and photometric match joined (where one existed).The total number of objects in this table is 495,003,196 6 .From these, we selected sources using the following criteria: (i) W1 − W2 ≥ 0.8: This colour threshold has been shown by Stern et al. (2012) to be an efficient AGN selection criterion, yielding AGN samples with very high (∼95%) reliability and good (∼80%) completeness at high X-ray luminosities.AGN are expected to heat dust to temperatures of several hundred degrees K, approaching sublimation, resulting in an SED peaking at a few microns and producing a red 1 − 2 colour.Whereas the colour cut alone can be contaminated by cool brown dwarfs, dust-reddened stars and star-forming galaxies (e.g.Stern et al. 2012;Yan et al. 2013;Hainline et al. 2014;Assef et al. 2018), our additional luminosity selection of only the most powerful sources is expected to weed out most such contaminants, preferentially selecting AGN alone (see step (vi)).We will return to this point in Section 6.
(ii) WISE -SDSS cross match distance -required to be <1 ′′ as discussed above.
(iii) Reliable redshift: • z_ERR/z ≤ 0.01 -a fractional error of 1 % or less on the redshift, indicating an accurate value.This criterion was introduced because a good zWarning flag does not necessarily imply a robust redshift measurement.Note that the median fractional  uncertainty in our final sample is 0.00016, i.e. much smaller than the adopted threshold of 1 % which is designed to weed out the most obvious unreliable sources.
• zWarning flag (SDSS) = 0, indicating high confidence in the redshift.Possible causes for a warning flag include poor or inconclusive fits, insufficient wavelength coverage, or problems with the instrumentation.
(iv) Redshift limit: 0 <  ≤ 0.83 (for SDSS spectroscopy) or 0 <  ≤ 1.065 (for BOSS spectroscopy) -this ensures that all spectra cover the redshifted [O iii]5007 line.For SDSS, this is identical to the redshift cut adopted by Reyes et al. (2008) for selecting [O iii]luminous Type 2 quasars, thus allowing us to contrast this class of objects with our targets.We additionally include spectra from newer survey programmes, and as these use a wider wavelength range spectrometer, we can include some further objects.
(v) Reliable WISE data: • WISE cc_flags = '0000'.The four zeroes refer to the four WISE bands and indicate that the WISE data are not affected by any known artifacts that can cause confusion or contamination in any of the bands.
• WISE W?ph_qual ̸ = 'U' AND W?ph_qual ̸ = 'X' AND W?ph_qual ̸ = 'Z'.This ensures a detection, reliable photometry, and measurable uncertainty in all WISE bands (with a signal:noise always above 2 and in the majority of sources above 3), allowing a robust measurement of the MIR source monochromatic fluxes (see next point) as well as MIR SEDs.
(vii) SDSS classification: CLASS = 'GALAXY' and SUBCLASS = '' -the CLASS selection removes sources with detected emission lines characteristic of QSOs (stars already having been culled by the redshift criteria).The SUBCLASS selection further removes sources with emission lines characteristic of Seyfert and LINERs, as well as star-forming and starburst systems7 .
(viii) Best available spectrum: SciencePrimary=1.This indicates that the spectra are considered by SDSS to be the best available for each object.
(ix) No significant detection of [O iii]5007.The emission-LinesPort table contains detailed flux fitting of a large number of SDSS spectra, but has the disadvantage of not having been applied to many newer spectra.In order to include the maximum number of possible spectra, we measured the significance of the [O iii]5007 line flux directly, similarly to the method used to fit [O iii] lines when selecting for QSO2s (Reyes et al. 2008): Gaussian curves were fit to rest frame positions of [O iii]5007 in the spectral data from SDSS.The outcome of this trial was that a subset of the objects had no apparent [O iii] emission line, and therefore these will comprise the final sample.A more detailed explanation of this trial is presented in Appendix B1. (x) Continuum SNR limit: Estimate of SNR in the direct region around where an [O iii] line would be detected.This is to remove any objects where the noise in the region of interest is overwhelming the measurement, leading to artificially high upper limits for the line flux.A more detailed explanation of this trial is presented in Appendix B3.
(xi) Visual inspection of sources: Automatic classification is not perfect, and a number of sources with e.g.obviously wrong redshift, classification, or visible emission lines were rejected in this final step.
An overview of the selection steps can be seen in Table 2.

Sample Contaminants
Although the selection steps are designed conservatively to produce a clean sample, we must consider the possibility of contaminant objects; i.e. sources that pass all the same criteria but are not predominantly AGN powered.
• Blazars (as discussed further in Section 6.3.7) are highly luminous and are generally lacking strong emission lines.However, the shape of their optical spectrum is often very distinct, and any clear blazars were removed in step (xi).
• Luminous dusty star-forming galaxies may reach luminosities comparable to the OQQ threshold in rare cases.To assess this probability, we compare with the catalogue produced by Chang et al. (2015), selecting non-quiescent galaxies (using the threshold from Carnall et al. 2023) and processing them through the OQQ MIR selection.We find that only ∼0.01% would pass the colour and luminosity criteria.If we also include the SDSS spectroscopic requirement for a null classification based on emission lines, no sources remain.This is logical, as any star formation intense enough to produce MIR emission on this level is likely to show strong emission lines on a galactic scale.The amount and distribution of dust that would be required to reduce the fluxes of these lines to the low levels detected (or upper limits of non-detections) would be highly unrealistic.
• Compact Obscured Nuclei (CONs) (as discussed further in Section 6.3.1) may represent a more likely class of contaminants, although still low probability overall.These sources contain a dense, IR-bright core that may be powered by compact star formation or an AGN, but are usually found in LIRGs, ULIRGs, and disturbed systems which typically show strong optical emission lines.CONs are currently thought to be rare, but both their intrinsic power source and true number counts are unknown.We should note that they would only be considered contaminants if they are powered by starformation -if they contain a strong AGN and reside in quiescent galaxies they would be correctly included as OQQs.

Comparison Population
As a comparative and complementary sample, we use a combination of: • The SDSS-selected sample of Type 2 quasars from Reyes et al. (2008).These are sources that are high luminosity ( bol > 10 45 erg s -1 ), show a significant [O iii]5007 line, and are optically obscured (Reyes et al. 2008).
• A later, but similar, sample from Yuan et al. (2016) selecting QSO2s from both SDSS and BOSS spectroscopy.These objects provide a valuable counterpoint to our selection, which is also based on high luminosity, but with no detectable [O iii]5007 flux.We combine the two tables, removing any duplicates.A small minority of sources were selected by Yuan et al. (2016) as having the wrong redshift fit by SDSS, and their corrected redshift is used in the processing.These are flagged in the table.The same procedure as used for the OQQs is used to match these QSO2s with WISE and the same luminosity cuts as outlined in Section 3 are applied to the combined catalog to produce the final comparison set.The sample includes only objects at z<1.06, the range where [O iii] can be detected in BOSS.

General sample properties
The number of sources that cumulatively pass all selections is 47.We refer to this sample as 'Optically Quiescent Quasars' or OQQs8 .Information on the final list is in Table 1.A sample of their spectra are in Figure A1.The sources span a range in redshift of 0.33 ≤  ≤ 0.94 (see Figure 2).No sources show significant [O iii]5007 lines, by selection, but other low level emission lines are present in several cases.
Figure 2 (main panel) shows the distribution across 12 µm luminosity and redshift for the OQQ sample and the QSO2 comparison set. Figure 2 (bottom) shows the redshift distribution of both the OQQ sample and the QSO2 comparison sample, and Figure 2 (right) shows the distribution with  12 .The K-S statistic comparing the OQQ and QSO2 distributions is 0.073 for  12 and 0.209 for redshift, with p-values 0.874 for  12 and 0.007 for redshift.These results indicate that the samples are likely to come from populations with the same distribution in terms of luminosity, but with a different distribution in redshift.This is partly due to the combined selection of QSO2s from different studies, and the upgrade of the SDSS spectrometer -the original SDSS spectrometer would only allow measurements of [O iii] up to  ∼0.83, hence the jump in distribution at this point.This does not invalidate their use as counterpart objects, but may have implications for any detailed analysis of evolution over time.

Continuum Properties
With regard to the infrared properties of the sources, the WISE colours are plotted in Fig. 3 on the canonical WISE colour-colour plane.The colours of all sources are consistent with being AGNdominated in the MIR according to the Stern et al. (2012) colour cut, by selection.In addition, most of the sources also lie within the WISE AGN wedge proposed by Mateos et al. (2012).
Fig. 4 compares the WISE and SDSS luminosities of the sources.For each source, we plot the monochromatic 12 µm and 5100 Å observed (i.e.not absorption-corrected) rest-frame luminosities (in   ).
The same figure also includes QSO2s: as already mentioned, QSO2s are likely to be a complementary sample to OQQs.All QSO2s with a WISE counterpart and MIR luminosity  12 > 3 × 10 44 erg s −1 are included (as described in Section 3.2), mimicking our selection strategy, and resulting in a sample of 1,990 QSO2s.The OQQ sample is reasonably well matched to the sample of QSO2s: at 5100 Å the average luminosity of the QSO2s is slightly dimmer, and the same (to a lesser extent) at 12 µm.In the following steps, we fit lines to the OQQ and QSO2 populations in order to estimate the reddening effect of the obscuring material on an average rather than individual source level.This line should not be taken as an indication of a relationship between the two properties -as we are only fitting upper limits it will be heavily effected by e.g. level of obscuration and/or intrinsic luminosity of the source.For this reason we have not made any special considerations in the fitting method for the fact that these are upper limits.This fit is based upon the ordinary least squares bisector method of Isobe et al. (1990).There is significant scatter of 0.11 dex and a Spearman rank correlation coefficient of 0.40 indicates that the relation is not very strong.Nevertheless, this serves as a guide for comparison to OQQs.

Emission Line Properties
The dotted line denotes an [O iii] deficit of ∼253 × from the above fit.This is the mean deficit of our OQQs relative to the expected intrinsic [O iii] luminosity if they are typical quasars similar to QSO2s.Since these are upper limits, this should not be considered a true fit to the data, but merely allows us to estimate the scale of obscuration present.The actual deficits could be much larger, and potentially show a different relationship to 12 µm luminosity.One can perform a similar comparison using the Hα emission line, as it is less affected by reddening.The result is shown in Fig. 6.
The corresponding relation for QSO2s is 5007 fluxes from the source papers (Reyes et al. 2008;Yuan et al. 2016) were used, but Hα6562 were drawn from the emissionLinesPort table.Since we are only using the QSO2 sample for comparison, we did not perform our own Hα line flux measurements.Above  ≈ 0.35, when redshifted Hα is shifted to ∼ > 9000 Å, the sky noise can often render lines difficult for the pipeline to measure correctly.Therefore, we restricted the QSO2 sample to below this redshift for extracting robust automated Hα fluxes.This results in few QSO2s for the above fit, all of which have a significant Hα detection without having been specifically selected on this line, as one would expect for a standard quasar spectrum.The scatter in this case is 0.10 dex with a Spearman rank correlation coefficient of 0.48.The OQQs have a median distribution (either a detection or a limit) consistent with an Hα deficit of ≈ 27.This is less than the deficit found from [O iii] partly because of the lower reddening affecting Hα, and partly because redshifted Hα typically lies in the reddest portions of the spectra with strong noise from sky emission lines resulting in less sensitive non-detection flux limits.Another possibility is that Hα has a larger fraction of its flux coming from star formation, while [O iii] is likely primarily from the AGN.

Spectral Energy Distribution (SED) fitting
We use agnfitter (Calistro Rivera et al. 2016) to fit the SEDs of a selection of OQQs and determine contributions to the observed emission from different components of the AGN/galaxy system.This method uses the multiwavelength data available to find the most likely weighted contributions to the overall SED, using Markov Chain Monte Carlo (MCMC) to fit from a large library of models.The most important parameter for our analysis is the contribution to the SED from the AGN obscuring material.We select the candidates with available long wavelength data (2MASS 9 , Herschel, in addition to all four WISE bands, present in all OQQs by selection), as this is key to constraining some components; this totals only 9 sources 10 .Figure 7 presents the results from four OQQs that have reasonably well-constrained results.After experimenting with different MCMC lengths, we selected one sufficient to achieve an auto-correlation time indicating convergence for the OQQ AGNheated dust component (see Foreman-Mackey et al. 2013 and Cal-9 1.235 µm, 1.662 µm, 2.159 µm. 10Herschel bands used were any available out of 250 µm, 350 µm, and/or 500 µm.Of these nine sources: two have 2MASS only; five have some combination of 2MASS and Herschel; two have Herschel only. istro Rivera et al. 2016 for details).Settings used were: 100 walkers; two burn in sets (10,000 steps each); set length 200,000.
In some cases the hot dust component (e.g. the top two plots in Figure 7) rises towards the blue end.This would not occur in a purely physical model component, but it is an empirical template based on combinations of AGN SEDs from Silva et al. (2004), and the range of values these can take can be seen in Figure 1 of Calistro Rivera et al. (2016) and Figure 1 of Silva et al. (2004); they are based on  H and assumed AGN type.As can be seen in these figures, the peak remains unchanged, varying only with normalisation.With further data into the blue and near-UV, we may be able to take this as some indication that the OQQ in question is lightly obscured, or at least has some accretion disc emission visible, but currently degeneracy with 'host galaxy' emission (solid orange line) makes this conclusion premature.Nevertheless, while the general results should be used with caution, the key outcome is that a high normalisation of the 'hot dust' component (see Figure 7, thick purple line) is required to reproduce the MIR part of the OQQ SEDs regardless of the tail shape (i.e. the peak of the hot dust component at ∼6-20 µm in the rest frame).We conclude it is very likely each of these objects contains an AGN bright enough to heat its obscuring dust to produce this component.agnfitter makes two separate estimates of star formation rate (SFR) -one in the optical, and one in the FIR.The FIR SFR provides an estimate of potential obscured star formation, but we note that very few candidates have reliable data in that region, so we cannot make any conclusions about the whole OQQ population.This estimate (and that of stellar mass) is obtained from host galaxy template parameters (for more detail, see Calistro Rivera et al. 2016).

4000 Å break
The 4000 Å break (D4000 from this point on), an indicator of a significant population of old stars, was calculated as described in Balogh et al. (1999), as the ratio between red (4000 Å -4100 Å) and blue (3850 Å -3950 Å) continuum fluxes.These ranges are slightly more restrictive than in previous works (e.g.Hamilton 1985), as this was found to improve on the repeatability of D4000 between separate measurements, and shows less sensitivity to reddening.A large 4000 Å break could be indicative of significant host galaxy light dilution.The values are found in Table 1 and shown in Figure 8: 44 out of 47 (94%) objects have a significant D4000.The median value is 1.41 (including only the significant breaks, defined as >3 standard deviations above unity).The QSO2 sample found only 811/1988 (41%) had a significant D4000.

X-RAY OBSERVATIONS
In the absence of AGN spectroscopic signatures in the optical, X-ray observations arguably provide the best evidence for the presence of an AGN independent of the MIR signatures.We checked the heasarc archive11 for pointed or serendipitous observations by Xray satellites, and discovered several OQQs observed during targeted observations of other nearby objects.We also obtained joint XMM-Newton and NuSTAR (Harrison et al. 2013) data for our prototypical OQQ (see Section 6), and these results are discussed in Section 5.1.Details of a serendipitous observation with sufficient counts for fitting are presented in Section 5.2.The XMM-Newton mission has one of the best combinations of effective area and angular resolution of current X-ray missions, and is therefore an ideal tool to help us analyse the X-ray properties of OQQs.Chandra, which has very low background noise, can also be valuable in detecting dim sources.We also examine Swift-XRT data -this is generally less sensitive, but in cases where no XMM-Newton or Chandra data is available Swift-XRT can still be useful to place limits on the X-ray emission.For detected sources, we extracted source and background spectra using standard ftools tasks and recommended reduction steps for each instrument, and analysed the data using the xspec package (Arnaud 1996).We then assessed the false detection probability based on counts in the source region relative to counts in a background region, and this derived either estimated fluxes or upper limits for each target (see e.g.Lansbury et al. 2017).Figures 11 and 12 show the results of analysis of data from XMM-Newton, Chandra, and Swift-XRT for OQQ candidates with serendipitous observations -in total 11 separate OQQs (See Table D2).These are discussed further in Section 5.3.

X-ray Observations of OQQ J0751+4028
We obtained targeted observations of our prototypical candidate with the aim of confirming its AGN nature and constraining its properties.Following standard procedures with xspec we modelled the source, finding that it is under-luminous at 2-10 keV (compared to the IR-predicted intrinsic luminosity), and appears to be lightly obscured ( H ∼10 22 cm −2 ).This work is summarised here: for details of the observations and procedures, see Greenwell et al. (2022).
We fit several models, starting with an absorbed power law (cabs*zwabs*pow), moving on to more physically realistic models (mytorus, Murphy & Yaqoob 2009;BNsphere, Brightman & Nandra 2011) and finally a more complex model involving a lightly obscured leak in an otherwise Compton thick sphere of obscuration.Using Bayesian X-ray Analysis (BXA; Buchner et al. 2014) we can compare the results and determine the model most likely to have produced the data.We find that BNsphere is favoureda spherical obscuration model that implies a 'cocooned' obscuring structure may be a plausible explanation for the X-ray appearance of OQQ J0751+4028.This result is shown in Figure 9 (bottom panel).The SDSS spectrum of OQQ J0751+4028 is shown in Figure 10 (bottom spectrum), and (by selection) shows typical properties for an OQQ (compared to examples from the general sample in Figure A1 and online).For more detail see Greenwell et al. (2021).
Figure 11 shows OQQ J0751+4028 in context with the serendipitous observations and the empirical relationship between 12 µm and intrinsic 2-10 keV luminosity.The unobscured rest frame 2-10 keV luminosity is 4.39 × 10 43 erg s −1 , approximately a factor of 6 lower than expected given the 12 µm luminosity, but sufficient to categorically demonstrate the presence of an AGN.The photon index is harder than usual for an AGN (Γ=1.32 +0.21  −0.19 ), and obscuration is approximately Compton-thin but present (log ( H / cm −2 ) = 21.96+0.16  −0.36 ).If Γ is fixed to a more usual AGN-like value of 1.9, the obscuration increased only slightly to log ( H / cm −2 ) = 22.3.

Serendipitous X-ray observations of OQQ J1051+3241
One serendipitous OQQ had sufficient XMM-Newton counts across two observations (0781410101 and 0781410201; both only observed the source with the PN detector) to be analysed further with xspec.We fit this source (OQQ J1051+3241) with several models, as was done with the targeted observations of OQQ J0751+4028 (see Section 5.1 and Greenwell et al. 2022).We find the best fit model to be borus02 (Baloković et al. 2018), at fixed angle of inclination assuming viewing through the torus (cos  = 0.05).We do not find a significant Fe K emission line, and allowing the angle of inclination to vary does not favour any angle strongly.This object is more obscured and has a higher Γ than J0751: log ( H / cm −2 ) = 23.26±0.06,Γ=2.05±0.2(see Figure 9, top panel).For all sources (with the exception of our targeted observation; Section 5.1), we have assumed the same photon index as for this source, and used  WebPIMMS12 to estimate the luminosity of each source, given the net counts (or upper limit) observed.We find that the detected OQQs are underluminous compared with expectations.If the non-detected OQQs are equally underluminous, they may still be detectable with sufficiently long observations.The SDSS spectrum of OQQ J1051+3241 is shown in Figure 10 (top spectrum).In this case it shows a slight rise towards the  blue end, implying that this is not as strictly 'typical' of an OQQ as OQQ J0751+4028.

Discussion of X-ray results
Figure 11 shows the comparison between X-ray and 12 µm luminosity, including the predicted relationship between intrinsic X-ray and nuclear 12 µm luminosity from Asmus et al. (2015) and between X-ray and 6 µm luminosity (Stern 2015); for the latter, 6 µm luminosities were converted to 12 µm luminosities using a relationship derived from the QSO template of Hao et al. (2007).The OQQ X-ray luminosities are slightly lower than the relations would predict, but this may be due to the difference between observed and intrinsic IR luminosities -if the OQQs are heavily obscured (whether from the circum-nuclear material or larger scale obscuration), their intrinsic luminosities will be higher.OQQ IR luminosities refer to the total galaxy emission, not the nuclear emission alone, as with current data there is no reliable way to extract or convert this for individual objects.Figure 11 shows where obscuration-reduced MIR-predicted X-ray luminosities (Asmus et al. 2015;Stern 2015) would fall.This would imply that most OQQs lie somewhere between Compton thick and unobscured.However, if they are intrinsically similar to OQQ J0751+4028 (Greenwell et al. 2022) the obscuration depth is low for an obscured AGN ( H ∼10 22 cm −2 ), and thus the intrinsic values for the OQQs may be only slightly higher.Figure 12 shows the relationship between [O iii] and X-ray luminosity found by Lamastra et al. (2009), with the QSO2s falling where expected.Upper limits on [O iii] luminosity place the OQQs significantly below expectations, reinforcing the disconnect between observed [O iii] and intrinsic properties.D1).We also show [O iii] predicted values for intrinsic luminosities (solid line; Lamastra et al. 2009).

DISCUSSION
The combination of WISE all-sky MIR photometry and SDSS optical spectroscopy have enabled us to identify an interesting class of sources showing all the characteristics of luminous AGN in the MIR (based upon their colours and luminosities) but no obvious signatures of the corresponding AGN activity in the optical.
Multiple physical scenarios can explain these observations, principally based on whether the optical AGN signatures are obscured, or intrinsically absent.In the former, we can consider the OQQ in terms of the spatial scale of its obscuration.As opposed to the putative donut shape envisaged by the zeroth order AGN unified torus scheme, full or near full sky covering by dusty 'cocoons' could easily explain the optical absence of emission signatures by not allowing an extended NLR to form or be ionized in the first place.Such cocoons would still reprocess the intrinsic AGN emission to the MIR, and the high WISE luminosities (along with bright hot dust component, where SED fitting is available) of these sources strongly suggest that the underlying power sources are luminous quasars, as does the MIR comparison with SDSS-selected QSO2s.For the latter scenario, AGN emission lines are not currently produced by the OQQ.We can envisage this as a 'young' AGN -recently switched on, with the AGN radiation not yet ionising the NLR, and consequently no narrow emission lines.It is still likely that this is an obscured AGN, i.e. that it is being viewed through a dusty torus, which would hide any broad emission lines (produced closer to the supermassive black hole (SMBH), and therefore earlier in the switching on process).As in the 'cocooned' picture, reprocessed intrinsic emission can still produce the observed IR properties.
We previously identified a prototypical example of this class, described in Greenwell et al. (2021) and Greenwell et al. (2022).We found that it indicated a new subclass of AGN that does not easily fit into current models: • The spectral shape of the optical continuum is that of a typical galaxy, including showing a strong 4000 Å break.
• Strong upper limits on emission lines that may indicate the presence of an AGN -there is no [O iii] by selection.Others may be present at a very low level.
• The WISE IR colours are firmly within the AGN wedge (Mateos et al. 2012), a stricter criterion than that required for the main sample selection.
• The 12 µm luminosity is very high, consistent with the quasar regime: there is a high chance of a hidden AGN being present, and it is not likely that this comes from a different mechanism.
• X-ray luminosity is luminous enough to clearly indicate the presence of an AGN, but under-luminous compared to expectations from 12 µm luminosity (see 5.1).
• An upper limit on radio flux is available from FIRST (Becker et al. 1995): <0.695 mJy/beam.

Optical Reddening
(Non)detection of emission lines can be used to place constraints on the optical depth of the theoretical cocoons, if present (see Section 6.2.2).We assume that the QSO2s have low extinction to their NLRs (an approximation), and that therefore the relationship between 12 µm and [O iii] luminosity (Hα luminosity) (equations 1(2) above) can be used to predict the approximate intrinsic flux of a cocooned source at the same 12 µm luminosity (assuming that all OQQs had [O iii] signals at the upper limit of detection.The extinction  V can be calculated using this predicted flux, the actual measured flux, and the standard Milky Way extinction law of Cardelli et al. (1989) with an  V = 3.1.The resultant values of the extinction  V as derived from the two lines are listed in Table 1.Since QSO2s are also likely to suffer mild NLR extinction (Reyes et al. 2008), these should be regarded as lower limits on the true  V .The  V limits measured from the [O iii] non-detections range over 3.4-10.5mag, with a median value of 4.69.
What would be the corresponding column density of gas for this extinction?A factor of 100 deficit in [O iii] does not need extreme columns.Our median extinction  V = 4.69 mag is equivalent to a neutral gas column of  ⊥ H ∼ 7 × 10 22 cm −2 along the galaxy axis, assuming the gas-to-dust ratios found in AGN environments by Maiolino et al. (2001).The parameter  ⊥ H is used to denote the fact that these gas columns are relevant for the NLR, which is likely to probe angles out of the direct l.o.s., thus making the distinction from the  los H measured from X-rays in the previous section.Further speculations on the gas columns or spatial distributions with only two detected objects would be premature.The  ⊥ H values are lower limits inferred from the limits on the optical extinction, so it is unexpected that the  los H would be lower.Figure 13 shows the  ⊥ H values inferred from [O iii] extinction.In a small number of cases, this can also be calculated using Hα (see Table 1).Where a Hα-inferred value is available, the predicted  H is consistently lower, and as some of these emission lines are measured rather than upper limits, are likely to be closer to the true value.

Intrinsic Nature
It is difficult at this stage to make definitive statements about the true intrinsic nature of OQQs, but we can speculate on some likely possibilities.A schematic comparison between QSO2s and theories on OQQs can be found in Figure 14.In Section 6.3 we discuss how these theories might relate OQQs to various other object classes.

Young AGN
A form of changing look AGN (e.g.Yang et al. 2018), objects belonging to this class have been recently activated such that they heat the torus to produce strong IR signals.However, radiation has not extended far enough to ionise the outer NLR, where the [O iii] we normally see in QSO2s is produced (middle diagram in Figure 14).For example, Gezari et al. (2017) present an example of a 'switching-on' AGN that goes from the appearance of a LINER to a Type-I quasar over less than 1 year.No [O iii] is found in the 'before' or 'after' spectra, due to [O iii] following longer term AGN activity than other narrow emission lines, such as Hα, which increase strongly.If OQQs are an obscured analogue to this object, it follows that we are unable to detect other emission lines.They must be observed through the torus, as otherwise we would still observe the broad emission lines produced closer in.If this is the case, and we expect broad emission lines to be present, then detection of these in spectral ranges less subject to obscuration could provide evidence for this explanation (e.g.NIR measurements; Onori et al. 2017a).An analogous and opposite population would be those described in Saade et al. (2022) which have low X-ray luminosity when compared with their strong [O iii] emission, with one possible explanation being a 'switching-off' AGN.They find an example with NLR and nuclear MIR emission but low X-ray luminosity and conclude that this object has recently become inactive.Discovery of this object among their relatively small parent sample implies that AGN caught 'mid-change' may not be extremely rare; we might expect 'switching-on' AGN to be found in similar numbers.

'Cocooned' AGN
This type of OQQ would be more mature than a young AGN, and has similar intrinsic features to a QSO2, but with a 'cocoon' of obscuring material preventing either observation of the NLR (right diagram in Figure 14), or its formation entirely.This 'cocoon' is likely to be unstable, and as such these objects may only have a short lifetime, suggested by the small number counts found when compared to QSO2s.At  ≈ 0.5 for most of our sample, a physical size of 1 kpc corresponds to an angular size of 0.16 ′′ , so any circumnuclear dust on sub-kpc scales would certainly be unresolved in most groundbased imaging.If the cocoon is comparable in radius to the dust sublimation radius, it may be even smaller -pc-scale -and thus unresolvable even with instruments such as ALMA.However, high resolution imaging might be expected to show very late stage mergers or other disturbances in the centres of these galaxies, as the source of the 'cocoon', despite the 'cocoon' itself being unobservable.We searched public archives such as MAST13 and the Virtual Observatory14 for any high angular resolution imaging of our targets, but these fields have not been observed with HST, nor with any other large ground-based observatory.In SDSS imaging (as well as WISE), the targets do not show any special characteristics, preferred morphologies or orientations, nor obvious disturbances.This appears to support small physical scales for any putative dust cocoons in our sample.However, it must be cautioned that SDSS imaging is seeing-limited and not very deep, so for fainter sources, signatures of interaction may be missed.

Comparison with other object classes
Although the idea of optically quiescent AGN is not new, there are no close counterparts of our selected sample, at least in the local universe.In this section, we briefly compare and contrast other related samples with the OQQs.

ULIRGs
Our selected targets have total MIR powers within or close to the ULIRG regime.ULIRGs tend to be hosted in chaotic appearing systems, with little clear structure and sometimes lack important optical emission lines, but are not usually completely devoid of strong emission signatures.For example, Mrk 231 has much stronger Hα than [O iii]5007 emission (Lacy et al. 1982;Malkan et al. 2017).Several local ULIRGs are also considered as being 'buried' AGN obscured along all lines of sight (e.g.Imanishi et al. 2001;Oyabu et al. 2011).These are typically associated with sources that have LINER-like optical spectra with clear emission lines (Imanishi et al. 2007).
As discussed in § 4.3, 15 of the 28 sources for which redshifted Hα lies within the SDSS spectral range show significant Hα emission.These sources may be similar to ULIRGs such as Mrk 231, though the comparison with SDSS-selected QSO2s implies that Hα also suffers some extinction, as already discussed.The remaining objects clearly differ from local ULIRGs in their lack of Hα.This may be explained by correspondingly higher NLR dust covering factors in the 'cocooned AGN' scenario.It could also be a sign that the dust distribution in the galaxy is less centrally located, additionally obscuring Hα from sources other than AGN ionisation.
Compact Obscured Nuclei (CONs) are found in LIRGs and ULIRGs (e.g.Falstad et al. 2021).They are dense, IR bright cores found at the centre of an optically normal galaxy.The strength of the IR emission indicates a heavily obscured nucleus, with either a compact starburst or AGN powering the luminosity.One intensively studied example is NGC 4418 (Costagliola et al. 2013;Asmus et al. 2014;Ohyama et al. 2019), which appears optically as an Sa type galaxy, and is part of an interacting pair.Extreme visual obscuration is demonstrated by a very deep silicate absorption feature in the MIR, and high resolution imagery shows a small, optically thick core.Ohyama et al. (2019) find that a dusty wind and compact starburst activity are present in the core, but that an additional energy source in the core is still required to explain the observations; it is still currently unknown whether this is an AGN or an extremely compact starburst.In a study with NuSTAR, Ricci et al. (2021) do not detect it and find that a SF model fits the data better, implying either no AGN or extremely thick column densities.
CONs are not very well understood, and there are not many examples.Falstad et al. (2021) search LIRGs and ULIRGs for compact sources of HCN vibrational emission, which are only found in extreme environments.They find that 40% of ULIRGs and 20% of LIRGs contain a CON, noting that their sample size is not large, and that the selection in the FIR does introduce some biases.This indicates that CONs are relatively common in (U)LIRGs; however, (U)LIRGs themselves are not extremely common.
OQQs are also bright in the IR, but are at larger distances than these prototypical CONs tend to be.Without high resolution imaging, it is impossible to definitively say whether the OQQs are equally compact.Unlike NGC 4418, OQQs appear to be small red galaxies with no obvious signs of interaction.

Weak Line and Lineless Quasars
At higher redshifts, several studies have been made of 'Weak Line Quasars' (WLQs; e.g.Meusinger & Balafkan 2014).Similarly to OQQs, these appear as normal quasars in terms of continuum emission, but show weak or absent emission lines.Examination of WLQs ( ⪅ 2.5; OQQs:  ⪅ 1.0) with Chandra snapshots and multiwavelength analysis from Wu et al. (2012) and Luo et al. (2015) find spectra that rise in the optical blue end which are generally lacking in stellar absorption lines, in contrast to the majority of the OQQs.They tend to have weak but detectable X-ray emission.At higher redshift, Shemmer et al. (2010) examine the region around the H and [O iii] lines of two WLQs at  ≈ 3.5, finding extremely weak H and no detectable [O iii].Wu et al. (2012) show results consistent with a weak, or 'anaemic' BLR producing weaker emission lines, and Luo et al. (2015) suggest that a 'puffed up' inner accretion disc could be preventing ionisation of the BLR.WLQs are not selected or restricted in terms of MIR colour or luminosity, which is a key factor in our OQQ selection.Where possible, spectroscopic measurements of more emission lines, along with targeted X-ray observations, could shed light on the link between low redshift OQQs and higher redshift WLQs.
Laor & Davis (2011) suggest a cold accretion disc could produce 'Lineless Quasars' -if the accretion rate of the AGN is low, then no ionisation can occur, and no emission lines will be observed.They argue that in the case of a cold disc, the ionising component is the X-ray power law alone.These may cross over with WLQs if the X-ray luminosity is high enough, but will appear lineless if   ≪   .Considering the underluminous X-ray results of OQQ J0751+4028 (see Section 5.1) this may be one reasonable explanation for an intrinsic lack of lines, if that is the case.The optical-NUV-FUV continuum would be expected to be blue in shape, and at this time the wavelength coverage of OQQs does not extend far enough to further investigate this possibility.

Low X-ray Scattering Fraction AGN
An interesting population of AGN selected from hard X-ray surveys has been identified by Ueda et al. (2007).These sources show apparently low X-ray scattering fractions with scattering fractions at least an order of magnitude lower than observed in typical Seyferts, and weak [O iii] as expected from a geometrically thick torus (Ueda et al. 2015).Interstellar dust reddening is also likely to play a role in depressing the scattered flux, as evidenced from the fact that a significant fraction appear to be hosted in edge-on host galaxies (Hönig et al. 2014).Gupta et al. (2021) present a detailed study on 386 Swift-BAT AGN, looking at the relationship between scattering fraction and various AGN properties.They find a significant correlation between scattering fraction and ratio of [O iii] to X-ray luminosity -i.e.low soft X-ray scattering fraction correlates with weak [O iii].In the context of OQQs, it follows that if the [O iii] emission is intrinsically not present (as in the 'Young' AGN theory), there may also be less scattering of soft X-rays.However, if the obscuring material in the 'cocoon' is sufficiently thick to absorb soft X-rays, we may not be able to distinguish between the theories in this way.In the same paper they find that low scattering fraction AGN are also likely to have higher intrinsic column densities, compounding the problem of selecting obscured AGN via optical methods.
Since our sample is MIR-selected, it may be better suited to search for Compton-thick counterparts to low X-ray scattering AGN.Detailed comparison of the OQQs to this sample will require more X-ray observations than are available so far.

X-ray Bright Optically Normal Galaxies (XBONGs)
The class of objects known as XBONGs are X-ray selected AGN that are hosted in optically-normal galaxies lacking (optical) AGN signatures.This lack has been attributed to extinction along all lines of sight, or due to dilution by host galaxy light in a substantial fraction.However, as opposed to our targets, XBONGs tend to generally possess lower Seyfert-like luminosities (Moran et al. 2002;Comastri et al. 2002;Civano et al. 2007;Caccianiga et al. 2007;Smith et al. 2014).In such cases, it can be easier for the host galaxy to dominate over the AGN and dilute AGN optical signatures.This becomes much harder for the quasar-like luminosities for our targets, and we do not believe dilution to be important for the sample presented herein.
In order to demonstrate this, in Fig. 15 we compare the OQQ broadband SEDs with the low resolution broadband AGN template compiled by Assef et al. (2010).The SEDs have been normalised at 12 µm for comparison.All wavelengths referred to here are in rest-frame.The median optical-to-MIR flux ratio for our sample is 0.31 when measured at 5500 Å and 12 µm.The template AGN SED clearly lies above most of our sample in terms of median optical normalised fluxes, with an optical-to-MIR ratio of 0.40.This implies that host galaxy dilution alone is unlikely to be a major effect, if our sample harbours typical central AGN.Not all XBONG classifications are considered to be a result of dilution, and XBONGs that are good candidates to be fully covered (Maiolino et al. 2003) may constitute the lower luminosity end of our selected sample.

High Intrinsic Accretion Rate AGN
An anti-correlation between [O iii] and the relative intensity of Fe II has been shown in analyses of emission-line and accretion properties of AGN (Eigenvector 1; e.g.Boroson & Green 1992;Shen & Ho 2014).This relationship can be explained by differences in Eddington ratio -an intrinsically high accretion rate could be associated with a low [O iii] luminosity.To assess the plausibility of the scenario for statistically significant numbers of the OQQ population, accurate measurements of the black hole (BH) masses would be needed; while this is out of the scope of this paper, analysis of NIR broad lines that may be less obscured (e.g.NIR measurements; Onori et al. 2017b) could provide the necessary information.Results from OQQ J0751+4028 do not support this (Greenwell et al. 2021(Greenwell et al. , 2022)), as it is not luminous enough to imply a sufficiently high Eddington ratio.

Other IR-selected Quasars
Using the WISE-WISE colour plane, as done here for OQQs, Hviding et al. (2018) focus on a different region of the plane, choosing targets more likely to be heavily obscured.They also make no selection based on emission lines, and the majority of their sources appear more optically typical than the OQQs.However, while there is no direct overlap between these samples, they are likely to be close cousins with unusually heavy obscuration instead of unusually distributed.Further investigation of the targets with insufficient spectral lines or features to assign a spectroscopic redshift leads Hviding et al. (2018) to hypothesise that a subset of these may be powerful AGN with emission lines attenuated by thick line of sight obscuration.
Extremely red quasars (ERQs) are quasars that show an unusually high IR to optical ratio (e.g.Banerji et al. 2015;Glikman et al. 2022).Ross et al. (2015) select these objects based on a comparison between SDSS and WISE magnitudes, producing a sample of objects up to  ∼4.An interesting subset of these, at  ∼2-3, show extremely broad emission lines which are not easily explained by current models.Some of these ERQs could be a Type 1 analogue to 'cocooned' OQQs -i.e.sources with a high covering factor (not the complete 4 steradians) and possibly obscuring material on large scales, but that retain some (unobscured but possibly reddened) line-of-sight to their nucleus; therefore presenting with broad AGN emission lines.Spectropolarimetry results from Type 2 ERQs (Alexandroff et al. 2018) indicate that the visible emission lines are due to reflection from a dusty disc wind rather than a direct line-of-sight.A potential interpretation is that these are similar to 'cocooned' OQQs, but with a small opening that allows the opportunity for reflected emission to be visible.Furthermore, an attempt at OQQ selection in large surveys based on photometry alone or low-quality spectroscopic data may include a number of ERQs, therefore it is important to understand their physical differences and how they might be distinguished.
Another interesting population found with extremely high infrared luminosities are Dust Obscured Galaxies (DOGs).These objects are selected based on an excess of 24 µm (from Spitzer; Dey et al. 2008) or 22 µm (from WISE; Toba et al. 2015) emission compared to optical, and are thus selected to be underluminous in the rest-frame UV.They often show signs of intense star formation or obscured AGN activity (Dey et al. 2008;Toba et al. 2017), with MIR colours indicating AGN presence in a large number of them, particularly at higher luminosities (Toba & Nagao 2016;Toba et al. 2017).The high luminosity along with dust obscuration seem to place them into the high-growth post-merger phase of the galaxy merger scenario for AGN-host galaxy growth (Hopkins et al. 2008).Their selection process results in a subset of DOGs that lack broad emission lines, and thus may be closer to OQQs than Type 1 ERQs.However, evidence of star formation and the presence of strong narrow lines shows that they are not an identical population.Spectropolarimetry results may suggest an evolutionary connection between ERQs, Hot DOGs, and other reddened quasars (Assef et al. 2022).Glikman et al. (2012) select candidate red quasars based on 2MASS (IR) and FIRST (radio) data.A detailed examination of the radio properties of OQQs is outside the scope of this paper, but we searched the FIRST catalogue (Becker et al. 1995) and found 7 FIRST counterparts to OQQ candidates, of which four are brighter than 10mJy at 1.4 GHz.The OQQs are therefore not consistently radio bright or dim, and unlikely to be similar, population-wise, to radio-selected objects.

Blazars
One potential source of contamination in this sample is likely to be blazars -the chance viewing of AGN head on to a relativistic jet (Landt et al. 2004), producing a bright, variable source, strong in radio emission.These are known to be very bright objects, with some exhibiting very flat spectra.As such, there may be a number that pass the criteria specified in Section 3, but that are not relevant objects for this paper.It would be useful therefore to be able to estimate the level of such contaminants in our sample, and if possible to remove some in order to improve the reliability of the sample.D'Abrusco et al. ( 2014) created a catalogue of -ray emitting blazars, assembled from WISE sources.Comparison of overlapping sources from this catalogue (WIBRaLS) with our OQQ sample, and with respect to various parameters, could help identify the unwanted objects.Massaro et al. (2012) identified the WISE Gamma-ray Strip (WGS), which uses IR colour-colour diagrams to distinguish possible blazars.They divide these into two sub-categories -BZBs (BL Lac objects, with featureless optical spectra -galactic absorption lines or weak emission lines only), and BZQs (flat spectrum radio quasars, where the optical spectrum should show broad emission lines).According to these definitions, the most likely category for any contaminant objects will be BZB.Massaro et al. (2012) quantify the likelihood of a source being within the WGS, based on its location on three WISE colour-colour diagrams, and the uncertainties on those colours.Figure 16 shows the application of this method to our objects.This is not considered a sufficiently reliable classification to remove any sources from the OQQ list, merely to assess the level of contamination from BZB sources.
Furthermore, given the flat spectrum characteristic of blazars, they tend to have weak values of D4000. Figure 8 shows that the OOQs tend to show significant breaks, and this can be seen in the individual optical spectra (Figure A).Hence a visual check of the spectra substantially reduces the likelihood of blazar contamination (Section B2).
One property common to all blazar types is that they exhibit high, irregular variability across the whole spectrum.As such, this may be a useful way to distinguish likely blazars from OQQs.
After the completion of the main WISE mission, the NEOWISE project (Mainzer et al. 2011) continues to scan the sky, doing so approximately every six months (depending on sky location) with the W1 and W2 bands (W3 and W4 cannot operate without cryogenic cooling, and this was exhausted during the main phase).The primary science aim of this phase is the study of near Earth objects.However as the single exposure images and source magnitudes are available from the NEOWISE database 15 , it provides a measure of the multiyear timescale infra-red variability.
To approximately quantify the level of variability in each source, the data from each WISE measurement epoch was averaged, then these combined points fit to a time-invariant line. 2 was used to assess whether this invariant state could reasonably explain the observed variation in WISE detections 16 .If it could not, then intrinsic uncertainty was added to the points.The size of this intrinsic 15 http://wise2.ipac.caltech.edu/docs/release/neowise/ 16The W3 (12 µm) band magnitudes were chosen as most representative of the region of interest.A more detailed investigation would examine the effects with luminosities rather than magnitudes.uncertainty was used as an estimate of the variability of the source, and Figure 17 shows the results of this. Figure 18 shows the fitting with the OQQ showing highest variability, OQQ J1208+1159.This source is flagged as a WGS BZB (see Figure 16), and is detected in FIRST (Becker et al. 1995) as a compact source with integrated flux 67.31 mJy.This combined with the MIR variability suggests that it may be a blazar.
The number of objects found that are above the threshold for being within the WGS BZB is 17 out of 47 OQQ.The number of objects with matches in the WIBRaLS table is 8 out of 47 OQQsthese are all also above the BZB threshold.The number of objects with high variability is currently difficult to quantify without a comparison threshold.The QSO2 variability is slightly lower on average (mean 0.031/median 0.020 for OQQ, mean 0.026, median 0.017, for QSO2s) and includes no very high variability objects.
If we take the lower limit of high variability as 0.1, then there are two OQQs above this threshold, both of which are likely BZBs.Of the 47 OQQs below the threshold, 36% are likely BZBs.In comparison, 57 QSO2s (out of 1,990; 2.9%) pass the BZB threshold overall.An examination of the variability of a large sample of non-AGN and various types of AGN would be illuminating, but outside the scope of this work.Due to the high uncertainties on these methods, it is difficult to fully quantify the number of BL Lac objects that could be interfering in the results here: anything from 3% to 25% may be possible.

Evolutionary Context
It is important to consider the cosmological situations where we might find OQQs.Many works have examined the life cycle of galaxies and AGN during mergers (e.g.Di Matteo et al. 2005;Springel et al. 2005;Hopkins et al. 2006).As the two galaxies interact, the interstellar medium from both comes together, and the resultant enhanced density of matter leads to increased star formation, and (as the bulk of the merging galaxy reaches the nucleus of its pair) we may see increased AGN activity.This brightening starts as the AGN is still enshrouded in dust and gas from the merger, but as the emission increases in intensity, energy and matter feedback starts to clear out the obscuring matter.Our OQQ sample objects, for the most part, do not show signs of high SFR as may be expected if they lie in this part of the sequence.Thus, OQQs do not obviously fit the standard paradigm for AGN and host galaxy evolution.Springel et al. (2005) show that in a modelled merger of two disc galaxies, including AGN feedback and BH growth, the peak of SFR precedes the black hole accretion rate (BHAR) peak by a small amount of time.SFR begins to drop before the BHAR, so there may be a short period of time where the SFR is well below peak, but BHAR is high.As BHAR is proportional to AGN luminosity, the OQQ objects might be objects in this short transitional phase.We find fewer OQQs compared to QSO2s, suggesting that the OQQ duty cycle is shorter than the duty cycle of Type 2 quasars.Interestingly, the difference in median redshift between these two samples (e.g., Fig. 2) implies that there may be genuine evolutionary differences between their populations, with OQQs peaking earlier than QSO2s.However, incompleteness issues could impact both source samples at the upper redshift end, and the number of nearby OQQs may be artificially low simply due to easier detectability of faint [O iii].
The two galaxies are coalescing at this point so would not necessarily be clearly identifiable as advanced mergers, especially with the images currently available.Additional high resolution imagery could be valuable in identifying any late stage mergers among the sample.For the OQQs to be objects existing in this brief period of time, we may expect strong outflows to be visible, if the SFR is being quenched by AGN activity.These will manifest as blueshifted absorption lines, or broadened and asymmetric emission lines at the systemic velocity of the source.With the current spectra this is impossible to either confirm or rule out, but measurements of select likely candidates with more sensitive instruments could provide more insight into these scenarios.
Simulations by Blecha et al. (2018) suggest that MIR selection of AGN in mergers may miss AGN, particularly low apparent luminosity and very late stage mergers.This work selects for high IR luminosity objects in order to produce a more reliable candidate list; lower luminosity OQQs may be missed from MIR selection due to (a) selection being less complete in general for less bright objects, (b) as shown in Blecha et al. if OQQs are late stage mergers.
The absence of obvious merger signatures may instead point towards high mass ratio mergers, with the smaller galaxy not massive enough to significantly disrupt the morphology of the larger primary, nor massive enough to trigger substantial star formation throughout the body of the primary.However, an efficient mechanism to channel the gas from the donor galaxy to the central SMBH would be required in this case.Alternatively, whether secular processes such as efficient gas flows can trigger AGN activity without simultaneous star formation remains unclear.
Alternatively, the 'Young AGN' scenario presents a different evolutionary perspective.If the triggering event is recent enough, the AGN could be switching on, mimicking the phase described by Schawinski et al. (2015) as 'Optically Elusive AGN'.They show that there could be a period of ∼10 4 years in which the AGN is visible in X-rays but has not yet photoionized the extended area that will later become the NLR.As the dusty torus (the source of reprocessed IR emission, and the IR colour that we select for in this work) is closer to the central engine than the NLR, we expect that sign of AGN life to appear prior to any narrow optical emission lines.Schawinski et al. (2015) calculated the full duty cycle of an AGN as ∼10 5 years, and the optically elusive phase as ∼10 4 years.If we expect the IR bright OQQ phase to be some fraction of the optically elusive phase, this puts the OQQ lifetime at less than ∼10% of the 'normal AGN' duty cycle.The viewing angle of OQQs and QSO2s is likely to be similar -both viewed through the torus, obscuring any broad emission lines -so assuming they have similar covering factors, we would expect to see them in the same proportions as the full range of AGN types: approximately 10:1.We find ∼1000 QSO2s and 47 OQQs: approximately 20:1; encouragingly close, but somewhat fewer OQQs than expected, especially if we consider that some OQQs may be 'cocooned' rather than 'young'.As stated above, however, the 10:1 ratio is a lower limit, based on the full 'optically elusive' timescale.If every OQQ is a 'young AGN' our results are consistent with the IR bright stage being approximately two thirds of the 'optically elusive' phase.These are merely crude estimates that do not fully account for selection and modelling biases, much of which remain unknown.
An investigation of the morphologies of OQQs could provide some interesting insight into their relationship with their host galaxies; for example, Rigby et al. (2006) found from an examination of the morphologies of X-ray-selected optically-dull AGN that such sources exhibited a wider range of axis ratios (major-minor axis, assuming an elliptical shape) as compared to optically-active AGN which showed only very round axis ratios.They conclude that extranuclear dust within the host galaxy is responsible for an opticallydull nature.SDSS photometry is severely limited for this purpose and is unlikely to be able to accurately resolve the size or axis ratio of the host galaxy except in unusual cases.Future surveys, such as Euclid (Euclid Collaboration et al. 2022), will greatly improve on this.

SUMMARY
By selecting MIR-luminous sources with red WISE colours and a marked absence of optical emission lines, we have presented a new sample of quasars that are candidates for being enshrouded within dust cocoons.The overall sample shows broadband continuum properties very similar to those of SDSS-selected QSO2s, but with an [O iii]5007 emission line suppression factor of ∼ >250 (see Section 4.3).We show that host galaxy dilution is unlikely to be a major effect in most cases.
The idea of fully enshrouded AGN is far from new.Various models suggest that the bulk of AGN growth occurs in highly obscured phases (e.g.Fabian & Iwasawa 1999;Hopkins et al. 2006).The shape of the cosmic X-ray background spectrum also requires similar numbers of Compton-thick ( H > 10 24 cm −2 ) obscured AGN compared to unobscured.Ananna et al. (2019) calculate the required proportion at  = 1.0 as 56% ± 7% (50% ± 9% at  = 0.1).A number of dedicated legacy surveys with NuSTAR that are selected in the infrared are finding CT fractions ≳30% within ≲200 Mpc (e.g., the NuSTAR Local AGN NH Distribution Survey -NuLANDS, Boorman et al. in prep., and the [Ne V] survey; Annuar et al. 2020).In any evolutionary paradigm, this would naturally include fully obscured sources in the stochastic competition between gas feeding from large scales and AGN feedback, especially if there is a delay between the two competing processes (noting that NuLANDS currently includes optical classification as part of the sample definition; the preliminary infrared selection should still include such fully-covered objects).It may also imply the existence of obscured objects that do not occur as part of the merger-driven growth phase, but instead represent a separate path.
The receding torus model also suggests that the covering factor of AGN decreases with luminosity (e.g.Lawrence 1991) so fully covered AGN may be more likely to occur at low AGN luminosities.This issue is far from settled, however, with conflicting results from a variety of studies in different bands (e.g.Glikman et al. 2007;Lawrence & Elvis 2010;Toba et al. 2014;Assef et al. 2015).An element of bias can be introduced by requiring that optical emission lines be present for secure source redshift and type identificationany sample based on this selection criteria may be missing objects with some emission lines blocked, either completely or to a level below measurable.The strength of our selection criteria is to instead rely upon the reprocessed thermal spectrum (i.e.MIR colours) and luminosities, allowing us to search for objects in which potential 4 or near 4 covering completely extinguishes the optical lines.
• We show that unknown AGN with unusual properties can be selected with a combination of multi-wavelength selection techniques, finding 47 'Optically Quiescent' Quasars (OQQs; strong MIR objects with no [O iii] emission) -see Section 3.
• Comparisons with Type 2 quasars (QSO2s; similarly MIR bright, but selected for strong narrow [O iii]) show similar underlying properties in e.g.redshift, colour, and luminosity distributions, but fewer numbers of OQQs -see Section 3.2.
• The exception to this is a strong deficit between the detected [O iii] emission from QSO2s and the upper limits placed on the OQQ lines -see Section 4.3.
• We consider what physical explanations may produce the observed properties (Section 6.2), principally: -'Cocooned' AGN -an accreting black hole entirely shrouded in obscuring material.
-'Young' AGN -a recently switched on AGN, resulting in contrasting signals originating at different distances from the black hole.
-Intrinsically weak lined AGN -some properties of the accretion disc can cause ionisation of narrow lines to be less prominent.
• Comparisons with other populations of AGN show that OQQs do not fit into any of these groups, but may be similar to some, depending on their intrinsic nature -e.g., Weak Line Quasars may be similar to OQQs that are caused by unusual accretion activity that results in low ionisation; Compact Obscured Nuclei may be similar to 'cocooned' AGN -see Section 6.3.contributions directly funded by ESA Member States and NASA.This research made use of data from the NuSTAR mission, a project led by the California Institute of Technology, managed by the Jet Propulsion Laboratory, and funded by NASA.This research has made use of the NuSTAR Data Analysis Software (NuSTARDAS) jointly developed by the ASI Science Data Center(ASDC, Italy)and the California Institute of Technology(USA).The NASA/IPAC Infrared Science Archive (IRSA) operated by JPL under contract with NASA, and SIMBAD operated at CDS, Strasbourg, France, were the main databases queried.The USNO Image and Catalogue Archive is operated by the U.S. Naval Observatory.This research has made use of data obtained through the High Energy Astrophysics Science Archive Research Center Online Service, provided by the NASA/Goddard Space Flight Center.This research made use of Astropy 17 , a community-developed core Python package for Astronomy (Astropy Collaboration et al. 2013, 2018).Many thanks also to Taiki Kawamura for his helpful comments and suggestions.

B2 Visual check
All spectra that passed the automatic check were then inspected by eye, and discarded or flagged based on several criteria.These cuts were done with the aim of making a highly reliable rather than complete sample, so that we are as certain as possible about the relevant nature of each passing candidate.A summary of these results is shown in Table B2.
(i) Visible [O iii] -the automatic check sometimes fails to find a significant [O iii] signal, despite there being a small but clearly apparent peak.We visually examine each spectrum, and based on the region around [O iii] and the width and height of the visible peak, we discard or keep objects based on whether the [O iii] seems truly absent.
(ii) Noisy data -we discard any candidates where the noise in the region of the [O iii] line is visibly worse than the flux there, in order to remove any where the [O iii] is likely to be lost in the noise.The S/N is calculated in detail in a later step (Section B3) (iii) Incorrect class -candidates were selected based on SDSS spectroscopic classifications, but rarely an object is misclassified, and we discard any objects where the class appears unreasonable.
(iv) Uncertain redshift -candidates were already selected based on SDSS spectroscopic redshift showing no warning, but this automatic fit fails occasionally.We discard any objects where it is not clear how the redshift has been assigned.By design, the OQQ spectra show no [O iii], and as shown in Figure 6, Hα is also generally low, where measurable.However, some optical spectra show other indicators of AGN presence and/or other interesting features.
• Blue continuum shapes: although any sources with extremely blue spectra were removed in Section 3, some remain with a blue upturn below the 4000 Å break (for example, sources shown in Figure B3, top), which may indicate some contribution to the continuum from the accretion disc.This is a small number of sources and are unlikely to represent a major part of the OQQ population, but the combination of continuum shape and lack of emission lines may make these interesting sources.Alternatively, this blue emission could be due to recent star formation, especially if seen in combination with other signatures such as [O ii].
• Broad emission lines: generally emission lines are low level or narrow for the majority of the sample, but there are some exceptions.J0948+0958 in Figure B3 (middle) shows broad Hα emission, possibly implying a source with greater velocity dispersion than expected from the galaxy, and likely AGN-related.Hα is present in several more sources but generally less broad.It is outside of the clean spectroscopic range for a large proportion of the OQQs, making judgement of population properties in this case harder.
• SF vs. dead galaxy indications: some SED fitting shows a contribution to MIR luminosities from SF (see Section 4.4), but only a small fraction of OQQs have sufficient data to draw reliable conclusions about their properties from these fits.Nevertheless, it implies that not all of the OQQs are dead galaxies with very little star formation, and observed properties of some spectra also lead to this conclusion.Narrow emission lines such as [O ii] may indicate a younger stellar population and thus recent or ongoing star formation, and some sources (e.g.sources in Figure B3, (bottom)) may show this.A strong 4000 Å break and Ca H and K absorption lines could instead imply an older population.However, as discussed in Section 4.5, this may not be conclusive.
We can conclude from inspection of the spectra (along with other properties) that it is unlikely that OQQs represent a single intrinsic type.Rather, it is more likely that they represent multiple classes that, due to obscuration or accretion properties, appear similar in terms of OQQ selection (see Section 6.2).Future work is needed to effectively and conclusively distinguish between intrinsic properties.

B3 Continuum SNR check
To eliminate sources where the spectrum in the region of [O iii] is significantly noisy -i.e. that there may be [O iii] present at a low level -we estimate the noise as follows (summary of results in Table B2 as 'Poor S/N'): (i) Cut the spectrum to two regions on either side of the expected position of [O iii] (vertical bands in Fig. B1, bottom right panel), avoiding any other potential emission or absorption lines.
(ii) Estimate the noise in these regions as the standard deviation on the flux data.
(iii) Estimate the signal as the mean on the flux data.(iv) Discard any candidate OQQ with the ratio of these >2.

B4 Duplicate Spectra
For the fitting and selection process, and the properties in the main sections of the paper, we use only the spectra assigned by SDSS as  SciencePrimary (a flag that indicates that this spectrum is considered the best available for this source).However, for a minority of sources there are multiple spectra available, so we compare these to see if there is any significant difference in flux, class, redshift, or emission line flux.Most of these duplicate spectra are either well matched, or are candidates that are dismissed for high levels of noise.We discuss the exceptions to this below.
B4.1 J0002: 0.592036, -0.4320923Two available spectra; generally consistent in appearance and flux, but one has been assigned the SDSS subclass 'Star-forming' (our selection process requires subclass 'null').This may be due to a change in emission line height.We choose to keep this candidate based on its otherwise promising appearance..2 J0231: 37.99002, 0.648491 Fifteen available spectra; generally consistent in appearance and flux, but again one has been assigned the SDSS subclass 'Starforming'.We also choose to keep this candidate based on its otherwise promising appearance; as the potentially star-forming galaxy is only one of the available spectra, and the rest are 'null', as required.B4.3 J0911: 137.96392, 29.824825Three available spectra.All with consistent SDSS class and redshift, but clear changes in flux between spectra, particularly noting a change in shape in one.The primary spectra is the most recent, so we have retained this object in the sample, but note that it may be variable.

APPENDIX C: COLOUR IMAGES
This section presents PanSTARRS colour images of a sub-sample of OQQs in Figure C1.The full sample is available online.2).Labelled in green ('[O iii] not detected') show the objects that pass this test and move on to visual checks, and those labelled in red ('[O iii] detected') are discarded.Grey sections of the spectra show regions removed for this test to avoid other possible emission lines biasing the fit.

APPENDIX D: X-RAY DATA SOURCES
Here we summarise the source papers used for QSO2 X-ray information and observation IDs of OQQ data used in Section 5.    46227001,46227006,46227004,46227005,46227003

Figure 1 .
Figure 1.SDSS-WISE nearest counterpart cross-matching distance distribution.Numbers in orange indicate the fraction of sources in the accepted (left) and rejected (right) sections.

Figure 2 .
Figure 2. The 12 µm luminosity and redshift distribution of the OQQ sample and the QSO2s.

Fig. 5
Fig. 5 compares the 3 limits on the [O iii] line luminosity with the observed MIR power.The dashed line denotes the relation  [O III] = 4 × 10 −5  12 , which is the median ratio of the [O iii] limits to the 12 µm detections for the OQQs.The plot shows that our sources have a similar MIR luminosity distribution as QSO2s above the selection threshold of  12 > 3 × 10 44 erg s −1 .On the other hand, the  [O III] limits of OQQs clearly lie much below the corresponding typical line detections for QSO2s.

Figure 3 .
Figure 3.The WISE colours of the OQQ sample plotted on the standard colour-colour plane.The dashed line shows the Stern cut (Stern et al. 2012), which all candidates pass by selection.The solid lines are the AGN wedge as proposed by Mateos et al. (2012); most OQQ also pass according to this test.

Figure 4 .
Figure 4. Comparison of WISE, SDSS and bolometric luminosities.OQQs are represented by outlined blue squares, and QSO2s by small green diamonds.Our luminosity threshold is shown by a horizontal dashed line.The dotted line represents where we would expect the points to lie if the IR and optical luminosities were equal.

Figure 5 .
Figure 5.A comparison of the 3 upper limits on the [O iii]5007 luminosity of the OQQ (large outlined blue square markers) and of the QSO2 sample (small green diamond markers), and their respective 12 µm luminosities.The dashed line shows the slope of the relationship between QSO2  12 and measured [O iii], scaled down by the average deficit.

Figure 6 .
Figure6.A comparison of the Hα6562 luminosity (or its upper limit -out of the sample of 62 OQQs, 17 have detected Hα emission lines), with OQQs (large blue square markers) and QSO2s (small green diamond markers), and their respective 12 µm luminosities (for those objects where Hα6562 is not redshifted out of the SDSS range).The dashed line shows the slope of the relationship between QSO2  12 and measured Hα, reduced by the average deficit.

Figure 7 .
Figure 7. agnfitter results for four candidate OQQs with long wavelength IR data.The hot dust component -that most indicative of an AGN -is highlighted (thick purple line).Errors are shown for this component and the total only.Photometric bands used in this plot (in addition to SDSS ugriz and WISE W1-4) were as follows: OQQ J0002-0025 all 2MASS and 250 µm; OQQ J0143+0151 all 2MASS and the three Herschel bands stated in footnote 10; OQQ J0231-0351 three Herschel bands; OQQ J1454+1440 350 µm only.

Figure 8 .
Figure8.4000 Å break data.Blue squares are the OQQs, with error bars shown for those with D4000 significantly greater than 1. Green markers are the QSO2 comparison sample.The average D4000 is clearly higher for OQQs than QSOs, implying a difference in the host galaxy population.However, this may be a selection rather than intrinsic effect.

Figure 10 .
Figure 10.SDSS spectra of OQQ J1051+3241 and OQQ J0751+4028.J1051 shows a slight blue rise which may indicate that there is either some view of the accretion disc or recent star formation (see Appendix B2), whereas J0751 shows the continuum shape and lack of emission lines of typical OQQs.

Figure 11 .
Figure 11.Observed 2-10 keV luminosities of the OQQs (blue markers), plotted against their 12 µm luminosities.OQQ J0751+4028 is shown as a large cross marker.For comparison, available QSO2 luminosities are shown as diamonds (see TableD1).We also show MIR predicted values for intrinsic luminosities (dashed line;Stern 2015).If this relationship is scaled to the luminosity of J0751 (pale grey;Greenwell et al. 2022), the other detected OQQs are more consistent with this expectation.Upper limits may be better explained by a Compton thick scenario: dotted blue line shows the obscuration of the S15 relation by Compton thick material.
Figure 11.Observed 2-10 keV luminosities of the OQQs (blue markers), plotted against their 12 µm luminosities.OQQ J0751+4028 is shown as a large cross marker.For comparison, available QSO2 luminosities are shown as diamonds (see TableD1).We also show MIR predicted values for intrinsic luminosities (dashed line;Stern 2015).If this relationship is scaled to the luminosity of J0751 (pale grey;Greenwell et al. 2022), the other detected OQQs are more consistent with this expectation.Upper limits may be better explained by a Compton thick scenario: dotted blue line shows the obscuration of the S15 relation by Compton thick material.

Figure 13 .
Figure 13.Inferred  ⊥ H values from [O iii] (i.e.obscuring the NLR, distinct from the  los H related to X-ray obscuration).

Figure 14 .
Figure 14.Schematic showing (left) a typical QSO2, producing a spectrum with strong [O iii], (middle) an OQQ that is intrinsically 'young', with no [O iii] produced, and (right) a 'cocooned' OQQ, producing typical [O iii] emission that is obscured by a larger region of material, with the observer positioned to one side of the diagram.Inset boxes show a zoom in on the [O iii] region of the spectrum -a strong line in the QSO2, but not detected in the two OQQs.Components are as follows: black central SMBH, blue accretion disk, dark brown standard torus, green ionised gas, lighter brown 'cocoon'.Composition of dust/obscuring material may be different across various parts of these structures.

Figure 16 .
Figure 16.WISE colour-colour diagram (as in selection; see Figure 3) showing the categorisation of the OQQs.OQQs that pass in all three planes of Massaro et al. (2012) are outlined with yellow squares.

Figure 17 .
Figure17.Intrinsic spread added to force the WISE data to be consistent with an unvarying source -higher values indicate that the source is more variable, and less likely to be stable.

Figure 18 .
Figure 18.W2 variability of the OQQ showing significant variation (>0.1 mag).It appears to show random variability, although the large time steps in available WISE data make statistical assessment of the variability difficult.

Figure A1 .
Figure A1.SDSS spectra of the first eight OQQs by RA, smoothed to show the continuum shape of each spectrum.Dashed lines in dark red mark key features: the 4000Å break (see Section 4.5), [O iii], and Hα.Full sample available online.
) (iv) Fit function to data ( =  +  + ( 0 , , ℎ)) (v) Simulate an [O iii] line for parameters randomly selected from their uncertainty bounds (vi) Calculate line flux for each simulated line (vii) Look at lower bound on line flux from simulations -if less than zero, we can conclude that the line only has an upper limit (FigureB2shows spectra for a selection of [O iii] detected and non-detected candidates) (viii) Select subset to pass onto the next step.

Figure B1 .
Figure B1.Full spectrum, showing cut out with the region used in the fits.Purple + markers show the possible location of any spectral lines.Red star marker shows the location of [O iii].Vertical shaded regions show emission lines cut for the fit.Bottom right panel shows the region used for SNR calculations.Vertical shaded areas show the SNR estimate regions.Noise is estimated from these and compared to the continuum flux, as described in Section B3.

Figure B2 .
Figure B2.Candidate OQQs, randomly selected from the catalogue at Step 8 (see Table2).Labelled in green ('[O iii] not detected') show the objects that pass this test and move on to visual checks, and those labelled in red ('[O iii] detected') are discarded.Grey sections of the spectra show regions removed for this test to avoid other possible emission lines biasing the fit.

Figure B3 .
Figure B3.OQQ spectra showing selected spectra with the following properties: (top) spectra showing blue rises, (middle) spectra showing strong and/or broad Hα emission, (bottom) spectra showing narrow [O ii] emission.

Figure C1 .
Figure C1.PanSTARRS colour images of the first 8 OQQs by RA.Bars show approximate scale of objects based on redshift.Full sample available online.

Table 2 .
Process of cutting down objects.A full description of the steps is available in Section 3.

Table B1 .
Bounds imposed on fit parameters.Slope and intercept of background continua were allowed to vary freely, whereas parameters of the Gaussian fit to the emission line were restricted to physically likely values.Gaussian fit to find the presence of strong [O iii].However, we found that due to the low level of our target [O iii] lines, the second peak caused the fit to be confused by the noise.We therefore fit only one Gaussian.Next we assess the uncertainty on the line flux, with the aim of checking whether the lower bound is less than zero, and hence whether only an upper bound can be placed on the emission.The continuum can be ignored in this simulation, as it cancels out from the equation.Finally, we select OQQs from the candidate objects that pass all other tests (as outlined in Section 3). double

Table B2 .
Counts of objects removed or kept at all stages of checks outlined in this Appendix.

Table D2 .
In Section 5 we examined available serendipitous X-ray observations of OQQs, where available.The targeted observations obtained for J0751 are listed in Section 5.1.The remainder of instruments and observation IDs used are listed here.