How are Red and Blue Quasars Different? The Radio Properties

A non-negligible fraction of quasars are red at optical wavelengths, indicating (in the majority of cases) that the accretion disc is obscured by a column of dust which extinguishes the shorter-wavelength blue emission. In this paper, we summarize recent work by our group, where we find fundamental differences in the radio properties of SDSS optically-selected red quasars. We also present new analyses, using a consistent color-selected quasar parent sample matched to four radio surveys (FIRST, VLA Stripe 82, VLA COSMOS 3 GHz, and LoTSS DR1) across a frequency range 144 MHz-3 GHz and four orders of magnitude in radio flux. We show that red quasars have enhanced small-scale radio emission (~kpc) that peaks around the radio-quiet threshold (defined as the ratio of 1.4 GHz luminosity to 6 micron luminosity) across the four radio samples. Exploring the potential mechanisms behind this enhancement, we rule out star-formation and propose either small-scale synchrotron jets, frustrated jets, or dusty winds interacting with the interstellar medium; the latter two scenarios would provide a more direct connection between opacity (dust; gas) and the production of the radio emission. In our future study, using new multi-band uGMRT data, we aim to robustly distinguish between these scenarios.


Introduction
Quasi-stellar objects (QSOs), also known as quasars, are the most powerful class of Active Galactic Nuclei (AGN). Their extremely high bolometric luminosities (up to 10 47-48 erg s −1 ) are now known to be caused by accretion onto a supermassive black hole (SMBH; masses of 10 8 -10 9 M ) near the Eddington limit, which places them as some of the most luminous objects in the Universe. Due to the unobscured view of the SMBH accretion disc, which peaks in the ultra-violet (UV), the majority of Type 1 QSOs have very blue optical colors. However, there is a small but significant subset with redder optical-infrared colors (referred to as "red QSOs").
Although red QSOs have been well studied in the literature, there are still conflicting views on how red and blue QSOs are related [1][2][3][4][5][6][7][8][9][10][11][12]. Some studies suggest they are intrinsically the same objects, with a red QSO simply representing a blue QSO observed at a higher inclination to the putative dusty torus, while other studies suggest an evolutionary scenario whereby a red QSO is a short-lived phase in QSO evolution. We have recently found that Sloan Digital Sky Survey (SDSS) red QSOs show ∼2-3 times higher radio-detection fraction compared to typical blue QSOs, which cannot be explained via a simple orientation model [8][9][10][11]. This radio-detection enhancement is driven by red QSOs with compact radio morphologies (<5 ; <43 kpc at z = 1.5) and luminosities placing them around the radio-loud/radio-quiet threshold [8]. We confirmed these initial results with deeper and higher resolution data from the Very Large Array (VLA) Stripe 82, VLA COSMOS 3 GHz, and LOFAR surveys [9,10]. Pushing to lower values of radio-loudness, we found that the radio enhancement decreases to unity at extreme radio-quiet values. Analyzing a sample of radio-compact red and blue QSOs with 0. 2 resolution e-MERLIN data, we found that the radio enhancement occurs on nuclear-galactic scales [11]. Overall, these results are inconsistent with red QSO being blue QSOs observed at a higher inclination angle, whereby the line of sight intersects with the top of the torus. If this were the case, we would expect the opposite result; i.e., enhanced radio emission in blue QSOs due to Doppler boosting arising from a more face-on orientation of the radio jet. Therefore, there must be factors other than orientation causing this enhancement, e.g., an enhancement in small-scale jets or winds causing radio-emitting shocks in the interstellar medium (ISM), enhanced star-formation (SF), or differences in the accretion properties.
In this paper, we undertake, for the first time, self-consistent analyses using four different radio datasets with the same QSO parent quasar sample, to ensure there are no differences in the construction of our samples that could be driving the radio results 1 . In Section 2, we discuss the radio data and our self-consistent selection of the QSO color-selected samples used in this study. In Section 3.1, we use these data to explore the radio enhancement in red QSOs across four orders of magnitude in radio flux over radio frequencies spanning 144 MHz-3 GHz, and, in Section 3.2, we compare the radio morphological sizes of red and blue QSOs. In Section 4, we discuss the potential origin of the radio emission. Throughout our work, we adopt a flat Λ-cosmology with H 0 = 70 kms −1 Mpc −1 , Ω M = 0.3, and Ω Λ = 0.7.

Materials and Methods
In the following sections, we briefly discuss the radio data used in this paper (Section 2.1), the optical data used to construct our QSO parent sample (Section 2.2), and our final color-selected and radio-detected samples (Section 2.3). Survey Area (deg 2 ) 10 100 1000 Figure 1. The 5σ sensitivity versus radio frequency for the four radio surveys considered in this paper. The survey area is indicated by the size of the marker. The dotted grey lines indicate the equivalent radio sensitivities of each of the four surveys at different frequencies for compact radio sources with an α = −0.7: for an equivalent 1.4 GHz luminosity, the S82, LoTSS DR1, and C3GHz survey sensitivities are ∼3, 10, and 38 times deeper than that of FIRST, respectively.

Radio Data
An illustration of the sensitivity, frequency and sky coverage for each radio survey used in this paper is shown in Figure 1. This highlights how some radio surveys cover a large fraction of the sky, but are lacking in sensitivity (e.g., FIRST), whereas others are extremely sensitive, but only cover a very small area of the sky (e.g., C3GHz).  [15]). S82 has a 1. 8 spatial resolution taken primarily in the A-configuration. At its median sensitivity of 52 µJy/bm, it is roughly three times deeper than FIRST.
In the published S82 catalogue, peak flux densities (F peak ) are derived by fitting an elliptical Gaussian model to the source. Matching to the FIRST catalogue recovers over 97% of the FIRST-detected QSOs in this region [14].
The COSMOS field has also been observed at 1.4 GHz over 2 deg 2 in an earlier VLA radio survey [17]. This survey has a resolution of 1. 5 × 1. 4 and a sensitivity of ∼10.5 µJy/bm (∼5 times deeper than S82). To calculate the L 1.4 GHz , the 1.4 GHz flux from this survey is used, when available, and otherwise scaled from the 3 GHz flux using a spectral index of α = −0.7.

LoTSS DR1
The LOFAR Two-meter Sky Survey (LoTSS) is a 120-168 MHz survey of northern sky at a resolution of 6 [18]. The first full-resolution data release (DR1 [19]) covers 424 deg 2 in the HETDEX spring field (RA = 10 h 45 m to 15 h 30 m, Dec = 45 • to 47 • ). At it's median sensitivity of 71 µJy/bm, LoTSS is equivalently ∼3.5 times deeper than S82 and ∼10 times deeper than FIRST at 1.4 GHz, assuming a spectral index of α = −0.7. To calculate the L 1.4 GHz , the 144 MHz fluxes are scaled to 1.4 GHz using a spectral index of α = −0.7.

Optical Data: The SDSS DR14 Quasar Catalogue
The SDSS DR14 Quasar Catalogue contains 526,356 spectroscopically-selected QSOs across a 9376 deg 2 region, out to redshifts around z = 7 [20]. The catalogue includes previous spectroscopically-confirmed QSOs from the SDSS-I and II Legacy surveys [21] with QSOs targeted by the Baryon Oscillation Spectroscopic Survey (BOSS [22]) in SDSS-III [23] and the extended Baryon Oscillation Spectroscopic Survey (eBOSS [24]) in SDSS-IV. In this paper, the SDSS five-band optical photometry (u * g * r * i * z * ) is also utilized, corrected by the associated band-dependent Galactic extinction estimates. The quasar catalogue provides spectroscopic redshifts based on different estimators [20]; in this work, we use the most robust of these estimates (listed as Z in the catalogue).

Color-Selected and Radio-Detected Samples
To define our sample of color-selected QSOs we used the g * (4770 Å) and i * (7625 Å) band extinction-corrected photometry in the DR14 catalogue and followed a similar approach to our previous studies [8][9][10]. First, we restricted to a redshift range of 0.2 < z < 2.4 to ensure our color selection is not affected by the Lyman break. In order to control for the impact of redshift on the g * − i * color selections, we sorted the QSOs by redshift and constructed g * − i * distributions in contiguous bins of 1000 sources. Our red and control QSOs (hereafter, rQSOs and cQSOs) were then selected above the top 90th percentile and within the middle 50th percentile of the observed SDSS g * − i * distribution, respectively; therefore, the cQSOs represent typical QSOs in terms of their color. For our study, here, the color-selected QSO sample was then restricted to the various radio survey regions, using non-zero local root mean square (RMS) values either provided through a radio mosaic or an online matching server.
Due to the loss of flux at UV-optical wavelengths, dust-reddened optically-selected rQ-SOs will be biased towards higher luminosities compared to blue QSOs at similar redshifts. Therefore, in order to ensure that our results are not due to differences in the intrinsic luminosity distributions of rQSOs and cQSOs, we also matched the samples in both redshift and luminosity: for this, we use the rest-frame 6 µm luminosity (L 6µm ) since this provides an extinction-insensitive measure of the intrinsic luminosity (e.g., Reference [25][26][27][28]). To calculate L 6µm , we identified mid-infrared (MIR) counterparts from the the Wide-Field Infrared Survey Explorer (WISE [29]), an all-sky survey which provides photometry in four bands (3.4, 4.6, 12 and 22 µm). We used the NASA/IPAC query engine to match the SDSS DR14 QSOs to the All-Sky WISE Source Catalogue (ALL-WISE), adopting a 2. 7 search radius, which ensures a 95.5% positional certainty [30], and required a detection with a signal-to-noise ratio (SNR) of greater than 2 in the WISE W1, W2, and W3 bands in order to derive an accurate estimate of their L 6µm . For the final color-selected samples, we match in L 6µm and redshift, using a tolerance of 0.2 and 0.05 dex, respectively 2 , following a similar approach to our original work (Section 2.3.2 in Reference [8]). The g * − i * color versus redshift distribution for our L 6µm -z matched rQSO and cQSO parent samples are shown in Figure 2. The g * − i * distribution shows that, despite encompassing 50% of the parent QSO population, the cQSOs cover a narrow region of color space, compared to the broad swathe of color space encompassed by the top 10% rQSOs.
The resulting number of color-selected and radio-detected rQSOs and cQSOs for each of the four radio surveys explored in this work are displayed in Table 1; the radio-detection fraction and enhancement 3 in the rQSOs compared to the cQSOs are also shown. To quantify how many of the rQSOs and cQSOs in our sample are 'radio-quiet', we adopted the same "radio-loudness" parameter (R) as that first used in Reference [8], defined as the dimensionless quantity: We also used the same radio-loud/radio-quiet threshold of R = −4.2, equivalent to a mechanical-to-radiative power ratio of P mech,sync /P rad,L 6µm ≈ 0.001, which is broadly consistent with the classical radio-quiet/radio-loud threshold (often defined using a 5 GHz-to-2500 Å flux ratio), but is less susceptible to obscuration from dust (see Reference [8] for full details). Table 1. The columns from left to right display the: (1) radio frequency of each dataset, (2) radio survey area, (3) number of rQSOs and cQSOs in the parent samples within each radio survey region, (4) the number of radio-detected rQSOs and cQSOs (the radio-detection fractions for each color-selected sample is displayed in the brackets), and (5) the radio enhancement of rQSOs over cQSOs; the radio enhancement for the rQSOs is >1 in each sample. It should be noted that the samples S82 and C3GHz in italics are not matched in L 6µm due to lack of sources.

Results
In this section, we use the same color-selected parent QSO sample, matched in both redshift and 6µm luminosity, in order to apply self-consistent analyses of the radio properties of rQSOs, using radio data from four different radio surveys, spanning radio frequencies of 0.144-3 GHz and equivalent 1.4 GHz radio fluxes of ∼0.01-1000 mJy. Figure 3 displays the L 6µm versus L 1.4 GHz distributions for the four radio-detected samples used in our analyses. The differences in the distributions of the data is related to the depth and areal coverage of the radio surveys. The C3GHz has the most sensitive radio data, probing down to very low R values, but covers a small area and consequently few luminous radiosources are detected. By comparison, the other surveys cover large areas and detect QSOs across a broad range of the R plane, providing good coverage overall particularly from the sensitive LoTSS survey. Figure 4 displays the cumulative radio-detection fraction of our four samples as a function of the 1.4 GHz radio flux. Across four orders of magnitude in 1.4 GHz flux, the rQSOs have a factor 2-3 higher radio-detection fraction in all four radio samples.   Figure 1). The histograms display the 1.4 GHz distributions for the FIRST and LoTSS samples: there is a higher number of radiodetected rQSOs compared to cQSOs. Due to the small survey area, we do not match in L 6µm for the S82 and C3GHz samples.

Radio Enhancement in rQSOs
Splitting the radio-detected samples into contiguous R bins, we calculated the enhancement in the radio-detection fraction of the full color-selected sample for the rQSOs. This is calculated as the ratio of the radio-detection fraction of the rQSOs to the radio-detection fraction of the cQSOs in each sample (see our previous studies for more details [8][9][10]). Figure  5 displays the radio enhancement for the FIRST-(shaded orange), combined S82+C3GHz-(shaded blue), and LoTSS-(shaded magenta)-detected samples. All three of these samples peak around R ∼ −4.5 with a factor ∼3-6 excess in the radio enhancement for rQSOs, decreasing to around unity towards the extreme radio-loud and radio-quiet values; we will explore this decrease in Section 4.1. The spread in the shaded regions represent the error bars, showing both the peak and decrease in radio-detection enhancement is statistically significant for all the samples. Expanding on the work from Reference [8], who find an enhancement in the FIRST-detection fraction of SDSS DR7 rQSOs around the radio-quiet/radio-loud threshold, we confirm this result at a much higher significance with our DR14 sample. We also push to lower R values where we find a decrease in radio-detection enhancement, consistent with what has been seen previously with deeper radio samples [9,10], but now, for the first time, with the shallower FIRST data.  The shaded error region was calculated using the method described in Reference [31] and corresponds to 1σ binomial uncertainties.  is displayed as a dashed vertical line. The combined S82 and C3GHz data points are taken from Reference [9] and the LoTSS data points are taken from Reference [11]. The FIRST data points are calculated using our new self-consistent sample from this study. The shaded error region was calculated using the method described in Reference [31] and corresponds to 1σ binomial uncertainties: in all four samples, the peak and decrease in the radio-detection enhancement is statistically significant. The radio enhancement of the C3GHz rQSOs with radio emission dominated by SF is shown as a black star, and those with radio emission dominated by AGN processes are shown as green stars.

Radio Morphology
The majority of the radio-detected rQSOs that contribute to the radio-excess shown in Figure 5 have compact morphologies at the resolution of FIRST (5 ; ∼43 kpc at z = 1.5 [8]). However, these size constraints are just upper limits, and higher-resolution radio data has shown that the rQSOs are typically compact over galaxy scales (e.g., <16 kpc on the basis of 1. 8 resolution data in the S82 field [9]).
Most recently, we have obtained 0. 2 e-MERLIN radio data for 19 and 20 FIRST radiocompact rQSOs and cQSOs, finding that the majority remain unresolved (suggesting <2 kpc scales) but also revealing a statistically significant enhancement in the fraction of rQSOs that are resolved over the ∼2-10 kpc scales of the host galaxy [11]. A comparison of the radio image of two rQSOs at FIRST 5 resolution and e-MERLIN 0. 2 resolution is shown in the left of Figure 6: in the lower example, there is clear extended emission at e-MERLIN resolutions, while the upper example remains unresolved. Figure 7 displays the L 1.4 GHz versus projected size of the e-MERLIN rQSOs and cQSOs; resolved sources are shown as circles, and unresolved sources are shown as upper limits. The resolved and unresolved source displayed in the thumbnails of Figure 6 are outlined in black. This demonstrates that radio data sensitive to smaller spatial scales (e.g., VLBA observations at ∼0. 01) is required to directly measure the resolved sizes for the majority of the sources that contribute to the radio excess (see Section 4.2 for a more detailed discussion).  [11]; the top source remains unresolved at the e-MERLIN resolution but the bottom source shows ∼kpc-scale extended radio emission. (Right) radio flux versus frequency for the two sources displayed in the thumbnails, using 150 MHz data from TGSS (green circle), 1.4 GHz data from FIRST (magenta square), and 3 GHz data from VLASS (blue cross); a 7σ upper limit of 24.5 µJy is used for the TGSS-undetected source (top). The grey shaded regions demonstrate the four radio bands over which we have recently obtained uGMRT data for all of our e-MERLIN observed sources, allowing us to construct sensitive radio SEDs (PI: V. Fawcett). The dashed and dotted blue lines indicate the expected fluxes for a spectral slope of α = −0.5, −0.7, and −0.9, demonstrating the variety in radio spectral slopes within the FIRST-compact sources. The vertical lines show the frequencies expected for GPS/CSS-like sources of 0.03, 0.1, 0.5, and 2 kpc extents. This demonstrates that, for the unresolved e-MERLIN source, with the uGMRT data, we will be able to constrain the radio sizes to <2 kpc scales if it is a GPS/CSS-like source, below the resolution of the e-MERLIN data. Reference [11]. The resolved sources are displayed as circles and the unresolved sources are shown as upper limits. The two rQSOs displayed in Figure 6 are outlined in black: source 1153+5651 is resolved, and 1323+3948 is unresolved.

Discussion
In this paper, we have reviewed our recent research that identified fundamental differences in the radio properties between rQSOs and cQSOs. We also undertook a series of new analyses using the same color-selected parent QSO sample, which we then matched to four different radio surveys at various depths and frequencies. This allowed us to perform self-consistent analyses based on the same QSO sample across the different radio datasets. We confirmed our previous radio results, finding an enhancement in the radio-detection fraction of rQSOs compared to cQSOs that arises in radio compact systems around the radio-quiet threshold [8][9][10][11]. This small-scale enhanced radio emission found in red QSOs could be driven by SF, AGN processes, or differences in the accretion properties. In a recent study (Fawcett et al. (submitted)), we used broad-wavelength Very Large Telescope (VLT) X-shooter [32] spectra to explore the accretion properties of rQSOs and cQSOs, finding no significant differences between the rQSOs and cQSOs. Consequently, it appears unlikely that the radio differences between the rQSOs and cQSOs in our SDSS studies are driven by differences in the accretion properties. Therefore, in the following sections, we restrict our discussion to SF and AGN processes as drivers of the enhanced radio emission in rQSOs.

Origin of the Radio Emission: Star-Formation
It has been shown in previous studies that SF dominates the radio emission in the majority of classically radio-quiet quasars, although many can still be AGN/jet-powered [33][34][35]. In our previous study, we analyzed the spectral energy distributions (SEDs) for a sample of rQSOs and cQSOs and found no differences in the average star-formation rate properties [12]; this will be robustly tested in our future ALMA study. Therefore, it is not obvious that SF should be driving the differences we see in the radio properties of red and blue QSOs. In order to further test this, in an earlier study, we utilized the rest-frame (8-1000 µm) far-infrared (FIR) data available in the COSMOS field to provide constraints on SF [9]. We defined the origin of the radio emission as dominated by SF if the measured radio luminosity was within a factor of 3 of the radio-FIR relationship for star-forming galaxies, and AGN-dominated otherwise. Using these definitions, we then split the lower radio-loudness bins covered by our C3GHz sample into sources with radio emission either dominated by SF or AGN (black and green stars, respectively; see Figure 5). From Figure 5, it is clear that, at R < −5, SF-dominated sources have a radio enhancement consistent with unity (albeit low source statistics), expected if there are no differences in the SF properties of rQSOs and cQSOs. We also found the AGNdominated sources, although very uncertain, are consistent with being enhanced at all values of R, suggesting that the decrease in the radio enhancement of rQSOs at R < −5 is due to the increased contribution of SF to the radio emission. Therefore, the radio enhancement in rQSOs is likely driven by AGN mechanisms.

Origin of the Radio Emission: AGN
From the e-MERLIN observations of the FIRST radio-compact QSOs, we found a statistically significant excess in the fraction of rQSOs with kpc-scale extended emission when compared to the cQSOs. However, the majority of both the rQSO and cQSO sources remained unresolved on <2 kpc scales. The mechanism driving this compact radio emission could be wind-driven shocks, small-scale synchrotron jets, or larger jets with synchrotron selfabsorption. A potential signature of the latter scenario is a turnover in the radio SED due to self-absorbed synchrotron radiation from a frustrated jet, similar to Compact Steep Spectrum (CSS) and Gigahertz-Peaked Spectrum (GPS) sources, which have well-defined peaks around 100 MHz and 1 GHz, respectively [36], or a flat spectrum due to a self-absorbed core [37]. Our future work analyzing 4-band 150 MHz-1.4 GHz upgraded Giant Meter Radio Telescope (uGMRT) data of our e-MERLIN sample (PI: V. Fawcett), in addition to the VLA Sky Survey (VLASS [38]) 3 GHz data, will be crucial for constructing sensitive 0.144-3 GHz radio SEDs to test the main mechanism driving the radio emission. We already have some insight into what we might expect to find with the uGMRT data from our previous study, where we constrained the radio spectral slopes of rQSOs and cQSOs that were compact in LoTSS at two frequencies: 144 MHz and 1.4 GHz from LoTSS DR1 and FIRST, respectively [10]. From this analysis, we found a broad range in spectral indices that peaked around α ∼ −0.7, consistent with synchrotron emission by normal jets or winds, but found no significant differences between the rQSOs and cQSOs. The right plot of Figure 6 displays the radio fluxes from the 150 MHz TIFR GMRT Sky Survey (TGSS; [39]), or an upper limit for the undetected unresolved source, FIRST, and VLASS for the two sources shown in the left thumbnails. The dashed and dotted blue lines indicate the expected fluxes for a spectral slope of α = −0.5, −0.7, and −0.9, and the grey shaded regions are the radio frequencies covered by our uGMRT data, where we could see a spectral turnover in the radio SEDs of our rQSOs and cQSOs.
In addition to providing an insight into the emission mechanism, the turnover frequency of a peak in the radio SED also varies inversely with the linear size of the radio emission for GPS/CSS-like sources, providing sub-kpc constraints on the scale of the radio emission without the need for high resolution imaging [40]. In Figure 6, the vertical lines show the frequency turnovers expected for GPS/CSS-like sources of 0.03, 0.1, 0.5, and 2 kpc extents, demonstrating that, with our uGMRT data, we will directly constrain the number of rQSOs and cQSOs that have frustrated jets on ∼0.1-2 kpc scales. If the turnover frequency occurs at <150 MHz, then, this would imply synchrotron self-absorption by a frustrated jet on scales >2 kpc, resolvable for the majority of our sources with e-MERLIN (e.g., the bottom panel of Figure 6). Therefore, we will be able to rule out this scenario for the unresolved e-MERLIN sources that display no spectral turnover in their radio SEDs.
For the sources that show no evidence in their radio SED for synchrotron self-absorption, the radio emission could be driven by either small-scale synchrotron jets or wind-driven shocks. If the enhanced radio emission is due to small-scale synchrotron jets, then, it is not immediately clear why there would be a close connection between the radio emission and optical color. Exploring the origin of the red colors, in Fawcett et al. (submitted), we use VLT/X-shooter spectra to explore the cause of the red colors and find that the reddening in rQSOs is fully consistent with moderate amounts of dust (A V ∼ 0.1-1 mags), as also inferred in our recent work [12]. To further demonstrate this result, Figure 8 displays DR14 SDSS rQSO (red) and cQSO (blue) median composites, in the redshift range 1.0 < z < 1.6 and luminosity range 45.3 < log L 6µm < 47.0 erg s −1 . A dust-reddened cQSO composite (A V ∼ 0.5 mags) is shown by the dotted black line, which describes the shape of the red composite remarkably well. A direct connection between the presence of dust and radio emission could be through shocks, whereby a dusty wind from the nucleus interacts with the ISM. We have previously found evidence for larger FWHMs in the broader wing component of the [O III]λ5007 in rQSOs compared to cQSOs at z < 0.8; evidence for stronger winds within the rQSO population [12]. In Alexander et al. (in prep), we also find a link between rQSOs and Low-ionization Broad Absorption Line QSOs (LoBALs; BALQSOs that display additional absorption in low-ionization species, such as Mg II and Al III, and are known to host powerful outflows, e.g., Reference [41]); LoBALs tend to have redder optical colors on average and have enhanced radio emission (also identified in other studies, e.g., Reference [42]). Therefore, although we cannot rule out small-scale synchrotron jets, wind-driven shocks interacting with the ISM can more directly associate the excess radio with the presence of dust and is consistent with our previous results. 3 3 × 10 3 4 × 10 3

× 10
Wavelength (Å) 10 4 Normalised Flux (arbitrary units) rQSO cQSO cQSO +AV = 0.48 Figure 8. SDSS DR14 composites for our rQSO (red) and cQSO (blue) samples, in a redshift range 1.0 < z < 1.6 and luminosity range 45.3 < log L 6µm < 47.0 erg s −1 . The dotted black line shows the cQSO composite with a dust extinction of A V = 0.48 mags that is remarkably similar to the rQSO composite. This amount of dust extinction is consistent with that found from fitting the individual VLT/X-shooter spectra of a sample of rQSOs and cQSOs at z ∼ 1.5 (Fawcett et al. (submitted)) and in our SED analysis of SDSS rQSOs and cQSOs [12]. The major emission lines are indicated by a dashed vertical grey line.

Conclusions
We have undertaken new, self-consistent analyses using the same DR14 QSO parent sample, exploring the radio properties of a sample of 0.2 < z < 2.4 red and typical QSOs. We confirm previous results but for the first time utilizing four different radio datasets at various frequencies and depths: FIRST, S82, C3GHz, and LoTSS DR1. Overall, we find: • Red QSOs have enhanced compact 0.144-3 GHz radio emission that peaks around the radio-quiet threshold (see Figures 4-7): We find an enhancement in the radio-detection fraction of red QSOs across a 1.4 GHz radio flux range of 0.1-1000 mJy. Exploring the radio-loudness parameter, we show that the radio enhancement in all four samples, covering radio frequencies of 0.144-3 GHz, peaks around the radio-quiet threshold and decreases towards both the extreme radio-loud and radio-quiet ends. We confirm the results from our previous work exploring the FIRST properties of SDSS DR7 red QSOs, but now at a higher significance using a 10 times larger sample. We also push to lower radio-quiet values than previously explored, confirming the decrease in radio-detection enhancement found with the deeper radio data. We confirm that the majority of radiocompact red QSOs have radio emission on <2 kpc scales [11] and show that our future uGMRT study could provide sub-kpc constraints on the scale of the radio emission. • The enhanced radio emission in red QSOs is likely due to dusty winds or frustrated jets (see Section 4 and Figures 6 and 8): Comparing a red QSO SDSS composite to a moderately dust-reddened blue composite, we show that dust is the likely cause of the red colors. This result is also confirmed in our upcoming X-shooter study where we also do not find any significant differences in the accretion properties between red and blue QSOs (Fawcett et al. (submitted)). In a previous study, we used the FIR data in the COSMOS field to show that the enhanced radio emission in red QSOs is likely driven by AGN processes (see Figure 5), consistent with Reference [12]. A self-consistent scenario that explains these results, in addition to the radio enhancement, is that red QSOs reside in a more dust and gas rich environment, in which the radio emission is due to winddriven shocks or frustrated jets interacting with the ISM/circumnuclear environment.
Understanding the connection between red QSOs and the compact radio emission promises to provide a key insight on the relationship between red and blue QSOs: are red QSOs a transition stage between obscured and unobscured QSOs, or is the relationship more complex? We have found enhanced nuclear-galaxy scale radio emission in red QSOs compared to typical blue QSOs that could be driven by either frustrated jets, winds, or smallscale synchrotron jets. Due to the link with opacity/dust, emerging evidence from our studies suggests that the most plausible scenario is either dusty winds or frustrated jets interacting with the ISM/circumnuclear environment. With our future uGMRT+VLASS study, we aim to robustly distinguish between these two scenarios through the construction of sensitive 0.15-3 GHz radio SEDs.  Our previous studies had small differences in the selection of the QSO samples: e.g., different SDSS data releases and different sized control samples. 2 Due to a lack of sources, we are unable to match in luminosity for the S82 or C3GHz samples. 3 Defined as the radio-detection fraction of the rQSOs, divided by the radio-detection fraction of the cQSOs.