AGN Unified Model Research Articles

Unsupervised machine learning for the classification of astrophysical X-ray sources

ABSTRACT The automatic classification of X-ray detections is a necessary step in extracting astrophysical information from compiled catalogues of astrophysical sources. Classification is useful for the study of individual objects, statistics for population studies, as well as for anomaly detection, that is, the identification of new unexplored phenomena, including transients and spectrally extreme sources. Despite the importance of this task, classification remains challenging in X-ray astronomy due to the lack of optical counterparts and representative training sets. We develop an alternative methodology that employs an unsupervised machine learning approach to provide probabilistic classes to Chandra Source Catalog sources with a limited number of labelled sources, and without ancillary information from optical and infrared catalogues. We provide a catalogue of probabilistic classes for 8756 sources, comprising a total of 14 507 detections, and demonstrate the success of the method at identifying emission from young stellar objects, as well as distinguishing between small- and large-scale compact accretors with a significant level of confidence. We investigate the consistency between the distribution of features among classified objects and well-established astrophysical hypotheses such as the unified active galactic nucleus model. This provides interpretability to the probabilistic classifier. Code and tables are available publicly through GitHub. We provide a web playground for readers to explore our final classification at https://umlcaxs-playground.streamlit.app.

Comparing the host galaxy ages of X-ray selected AGN in COSMOS

We explore the properties of the host galaxies of X-ray selected AGN in the COSMOS field using the Chandra Legacy sample and the LEGA-C survey VLT optical spectra. Our main goal is to compare the relative ages of the host galaxies of the obscured and unobscured AGN by means of the calcium break Dn(4000) and the Hδ Balmer line. The host galaxy ages are examined in conjunction with other properties such as the galaxy stellar mass, and star-formation rate as well as the AGN Eddington ratio. Our sample consists of 50 unobscured or mildly obscured (NH < 1023 cm−2) and 23 heavily obscured AGN (NH > 1023 cm−2) in the redshift range z = 0.6 − 1. We take specific caution to create control samples in order to match the exact luminosity and redshift distributions for the obscured and unobscured AGN. The majority of unobscured AGN appear to live in young galaxies in contrast to the obscured AGN which appear to live in galaxies located between the young and old galaxy populations. This finding may be in contrast to those evolutionary AGN unification models which postulate that the AGN begin their life in a heavy obscuration phase. The host galaxies of the obscured AGN have significantly lower levels of specific star-formation. At the same time the obscured AGN have lower Eddington ratios indicating a link between the star-formation and the black hole accretion. We find that the distribution of the stellar masses of the host galaxies of obscured AGN is skewed towards higher stellar masses in agreement with previous findings. Our results on the relative age of obscured AGN are valid when we match our obscured and unobscured AGN samples according to the stellar mass of their host galaxies. All the above results become less conspicuous when a lower column density (log NH(cm−2) = 21.5 or 22) is used to separate the obscured and unobscured AGN populations.

A practicable estimation of opening angle of dust torus in Type-1.9 AGN with double-peaked broad Hα

ABSTRACT In this manuscript, an independent method is proposed to estimate opening angle of dust torus in AGN, through unique properties of Type-1.9 AGN with double-peaked broad Hα (Type-1.9 DPAGN) coming from central accretion disc. Type-1.9 AGN without broad Hβ can be expected by the commonly accepted unified model of AGN, considering central BLRs seriously obscured by dust torus with its upper boundary in the line of sight. For the unique Type-1.9 DPAGN, accretion disc originations of double-peaked broad Hα can be applied to determine the inclination angle of the central accretion disc, which is well accepted as substitute of the half opening angle of the central dust torus. Then, among low redshift Type-1.9 DPAGN in SDSS, SDSS J1607+3319 at redshift 0.063 is collected, and the half opening angle of the central dust torus is determined to be around 46 ± 4°, after considering the disfavoured BBH system to explain the double-peaked broad Hα through long-term none variabilities and disfavoured local physical conditions to explain the disappearance of broad Hβ through virial BH mass properties. The results indicate that more detailed studying on dust torus of AGN can be appropriately done through Type-1.9 DPAGN in the near future.

Coevolution of Broad Emission Lines and X-Ray Spectrum in Changing-look AGNs

Changing-look active galactic nuclei (CLAGNs) show the disappearance and reappearance of broad emission lines in a few years, which challenges the orientation-based AGN unification model. We reduce the X-ray data for five well-studied CLAGNs that show a strong change in broad emission lines in the past several decades. We find that the X-ray photon index, Γ, and the Eddington-scaled X-ray luminosity, L 2–10 keV/L Edd, normally follow negative and positive correlations when the Eddington ratio is lower and higher than a critical value of ∼10−3. We find that the CLAGNs observed with broad Hβ emission lines stay in the positive part of the Γ–L 2–10 keV/L Edd correlation, while the broad Hβ lines become weak or disappear in the anticorrelation part of the Γ–L 2–10 keV/L Edd correlation, which suggests that the evolution of the broad lines should be correlated with the evolution of the underlying accretion process. We further find that the CLAGNs are consistent with the other different types of AGNs in the L bol–L bol/L Edd correlation. These results support that the CLAGNs belong to a special stage of AGNs with a bolometric Eddington ratio ∼1%, where the broad emission lines are easily affected by the strong variation in ionization luminosity that is caused by the transition of accretion modes.

WISE View of Changing-look Active Galactic Nuclei: Evidence for a Transitional Stage of AGNs

Abstract The discovery of changing-look active galactic nuclei (CLAGNs) with a significant change in optical broad emission lines (optical CLAGNs) and/or strong variation of line-of-sight column densities (X-ray CLAGNs) challenges the orientation-based AGN unification model. We explore mid-infrared (mid-IR) properties for a sample of 57 optical CLAGNs and 11 X-ray CLAGNs based on the Wide-field Infrared Survey Explorer archive data. We find that Eddington-scaled mid-IR luminosities of both optical and X-ray CLAGNs stay just between those of low-luminosity AGNs and luminous QSOs. The average Eddington-scaled mid-IR luminosities for optical and X-ray CLAGNs are ∼0.4% and ∼0.5%, respectively, which roughly correspond to the bolometric luminosity of transition between a radiatively inefficient accretion flow and a Shakura–Sunyaev disk. We estimate the time lags of the variation in the mid-IR behind that in the optical band for 13 CLAGNs with strong mid-IR variability, where the tight correlation between the time lag and the bolometric luminosity (τ–L) for CLAGNs roughly follows that found in the luminous QSOs.

Giant Radio Quasars: Composite Optical Spectra

Abstract We present the composite optical spectrum for the largest sample of giant radio quasars (GRQs). They represent a rare subclass of radio quasars due to their large projected linear sizes of radio structures, which exceed 0.7 Mpc. To construct the composite spectrum, we combined the optical spectra of 216 GRQs from the Sloan Digital Sky Survey (SDSS). As a result, we obtained the composite spectrum covering the wavelength range from 1400 Å to 7000 Å. We calculated the power-law spectral slope for the GRQ’s composite, obtaining α λ = −1.25, and compared it with that of the smaller-sized radio quasars, as well as with the quasar composite spectrum obtained for a large sample of SDSS quasars. We obtained that the GRQ’s continuum is flatter (redder) than the continuum of comparison quasar samples. We also show that the continuum slope depends on core and total radio luminosity at 1.4 GHz, being steeper for higher radio luminosity bins. Moreover, we found that there is a flattening of the continuum with the increase in the projected linear size of the radio quasar. We show that α λ is orientation-dependent, being steeper for a higher radio core-to-lobe flux density ratio, which is consistent with AGN unified model predictions. For two GRQs, we fit the spectral energy distribution using the X-CIGALE code to compare the consistency of results obtained in the optical part of the electromagnetic spectrum with broadband emission. The parameters obtained from the SED fitting confirmed the larger dust luminosity for the redder optical continuum.

Three-dimensional Distribution Map of H i Gas and Galaxies around an Enormous Lyα Nebula and Three QSOs at z = 2.3 Revealed by the H i Tomographic Mapping Technique

Abstract We present an IGM H i tomographic map in a survey volume of (cMpc3) centered at MAMMOTH-1 nebula and three neighboring quasars at z = 2.3. The MAMMOTH-1 nebula is an enormous Lyα nebula (ELAN), hosted by a type-II quasar dubbed MAMMOTH1-QSO, that extends over cMpc with no clear physical origin. Here we investigate the H i-gas distribution around MAMMOTH1-QSO with the ELAN and three neighboring type-I quasars, making the IGM H i tomographic map with a spatial resolution of 2.6 h −1 cMpc. Our H i tomographic map is reconstructed with H i Lyα forest absorption of bright background objects at z = 2.4–2.9: one eBOSS quasar and 16 Keck/LRIS galaxy spectra. We estimate the radial profile of H i overdensity for MAMMOTH1-QSO, and find that MAMMOTH1-QSO resides in a volume with fairly weak H i absorption. This suggests that MAMMOTH1-QSO may have a proximity zone where quasar illuminates and photoionizes the surrounding H i gas and suppresses H i absorption, and that the ELAN is probably a photoionized cloud embedded in the cosmic web. The H i radial profile of MAMMOTH1-QSO is very similar to those of three neighboring type-I quasars at z = 2.3, which is compatible with the AGN unification model. We compare the distributions of the H i absorption and star-forming galaxies in our survey volume, and identify a spatial offset between density peaks of star-forming galaxies and H i gas. This segregation may suggest anisotropic UV background radiation created by star-forming galaxy density fluctuations.

Peculiar emission line spectra of core extremely red BOSS quasars at z∼2–3: orientation and/or evolution?

Context. Core extremely red quasars (core ERQ) have been proposed to represent an intermediate evolutionary phase in which a heavily obscured quasar blows out the circumnuclear interstellar medium with very energetic outflows before it becomes an optical quasar. Aims. We investigate whether the properties of core ERQ fit the AGN orientation-based unification scenario. Methods. We revised the general UV and optical emission line properties of core ERQ in the context of the orientation-based scenario. We used diagnostic diagrams based on UV emission line ratios and UV-to-optical line kinematic information to compare the physical and kinematic gas properties of core ERQ with those of other luminous narrow- and broad-line AGN. In particular, we provide a revised comparison of the [OIII] kinematics in 21 core ERQ (20 from Perrotta et al. 2019, MNRAS, 488, 4126 and SDSS J171420.38+414815.7, based on GTC EMIR near-infrared spectroscopy) with other samples of quasars with matching luminosity with the aim of evaluating whether core ERQ host the most extreme [OIII] outflows. Results. The UV line ratios suggest that the physical properties (e.g., density and metallicity) of the ionised gas in core ERQ are similar to those observed in the broad-line region of blue nitrogen-loud quasars. The [OIII] outflow velocities of core ERQ are on average consistent with those of very luminous blue type 1 quasars, although extreme outflows are much more frequent in core ERQ. These similarities can be explained in the context of the AGN unification model under the assumption that core ERQ are viewed with an intermediate orientation between type 2 (edge-on) and type 1 (face-on) quasars. Conclusions. We propose that core ERQ are very luminous but otherwise normal quasars viewed at an intermediate orientation. This orientation allows a direct view of the outer part of the large broad-line region from which core ERQ UV line emission originates; the extreme [OIII] outflow velocities are instead a consequence of the very high luminosity of core ERQ.

Discovery of 21 New Changing-look AGNs in the Northern Sky

Abstract The rare case of changing-look (CL) active galactic nuclei (AGNs), with the appearance or disappearance of broad Balmer emission lines within a few years, challenges our understanding of the AGN unified model. We present a sample of 21 new CL AGNs at 0.08 < z < 0.58 , which doubles the number of such objects known to date. These new CL AGNs were discovered in various ways, from (1) repeat spectra in the SDSS, (2) repeat spectra in the Large Sky Area Multi-Object Fiber Spectroscopic Telescope (LAMOST) and SDSS, and (3) photometric variability and new spectroscopic observations. We use the photometric data from surveys, including the SDSS imaging survey, the Pan-STARRS1, the DESI Legacy imaging survey, the Wide-field Infrared Survey Explorer (WISE), the Catalina Real-time Transient Survey, and the Palomar Transient Factory. The estimated upper limits of the transition timescale of the CL AGNs in this sample spans from 0.9 to 13 years in the rest-frame. The continuum flux in the optical and mid-infrared becomes brighter when the CL AGNs turn on, or vice versa. Variations of more than 0.2 mag in the W1 band were detected in 15 CL AGNs during the transition. The optical and mid-infrared variability is not consistent with the scenario of variable obscuration in 10 CL AGNs at more than the 3σ confidence level. We confirm a bluer-when-brighter trend in the optical. However, the mid-infrared WISE colors W1–W2 become redder when the objects become brighter in the W1 band, possibly due to a stronger hot dust contribution in the W2 band when the AGN activity becomes stronger. The physical mechanism of type transition is important for understanding the evolution of AGNs.

Galaxy and Mass Assembly (GAMA): active galactic nuclei in pairs of galaxies

There exist conflicting observations on whether or not the environment of broad and narrow line AGN differ and this consequently questions the validity of the AGN unification model. The high spectroscopic completeness of the GAMA survey makes it ideal for a comprehensive analysis of the close environment of galaxies. To exploit this, and conduct a comparative analysis of the environment of broad and narrow line AGN within GAMA, we use a double-Gaussian emission line fitting method to model the more complex line profiles associated with broad line AGN. We select 209 type 1 (i.e., unobscured), 464 type 1.5-1.9 (partially obscured), and 281 type 2 (obscured) AGN within the GAMA II database. Comparing the fractions of these with neighbouring galaxies out to a pair separation of $350\,\text{kpc }h^{-1}$ and $\Delta z < 0.012$ shows no difference between AGN of different type, except at separations less than $20\,\text{kpc }h^{-1}$ where our observations suggest an excess of type 2 AGN in close pairs. We analyse the properties of the galaxies neighbouring our AGN and find no significant differences in colour or the star formation activity of these galaxies. Further to this we find that $\Sigma_5$ is also consistent between broad and narrow line AGN. We conclude that the observations presented here are consistent with AGN unification.

Seyfert 1 composite spectrum using SDSS Legacy survey data

We present a rest-frame composite spectrum for Seyfert 1 galaxies using spectra obtained from the DR12 release of the Sloan Digital Sky Survey (SDSS). The spectrum is constructed by combining data from a total of 10,112 galaxies, spanning a redshift range of 0 to 0.793. We produce an electronic table of the median and geometric mean composite Seyfert 1 spectrum. We measure the spectral index of the composite spectrum, and compare it with that of the composite quasar spectrum. We also measure the flux and width of the strong emission lines present in the composite spectrum. We compare the entire spectrum with the quasar spectrum in the context of the AGN unification model. The two composite spectra match extremely well in the blue part of the spectrum, while there is an offset in flux in the red portion of the spectrum.

A DEEPHERSCHEL/PACS OBSERVATION OF CO(40-39) IN NGC 1068: A SEARCH FOR THE MOLECULAR TORUS

Emission from high-J CO lines in galaxies has long been proposed as a tracer of X-ray dominated regions (XDRs) produced by AGN. Of particular interest is the question of whether the obscuring torus, which is required by AGN unification models, can be observed via high-J CO cooling lines. Here we report on the analysis of a deep Herschel-PACS observation of an extremely high J CO transition (40-39) in the Seyfert 2 galaxy NGC 1068. The line was not detected, with a derived 3$\sigma$ upper limit of $2 \times 10^{-17}\,\text{W}\,\text{m}^{-2}$. We apply an XDR model in order to investigate whether the upper limit constrains the properties of a molecular torus in NGC 1068. The XDR model predicts the CO Spectral Line Energy Distributions for various gas densities and illuminating X-ray fluxes. In our model, the CO(40-39) upper limit is matched by gas with densities $\sim 10^{6}-10^{7}\,\text{cm}^{-3}$, located at $1.6-5\,\text{pc}$ from the AGN, with column densities of at least $10^{25}\,\text{cm}^{-2}$. At such high column densities, however, dust absorbs most of the CO(40-39) line emission at $\lambda = 65.69\, \mu$m. Therefore, even if NGC 1068 has a molecular torus which radiates in the CO(40-39) line, the dust can attenuate the line emission to below the PACS detection limit. The upper limit is thus consistent with the existence of a molecular torus in NGC 1068. In general, we expect that the CO(40-39) is observable in only a few AGN nuclei (if at all), because of the required high gas column density, and absorption by dust.

Extremely red quasars from SDSS, BOSS andWISE: classification of optical spectra

Quasars with extremely red infrared-to-optical colours are an interesting population that can test ideas about quasar evolution as well as orientation, obscuration and geometric effects in the so-called AGN unified model. To identify such a population we match the quasar catalogues of the Sloan Digital Sky Survey (SDSS), the Baryon Oscillation Spectroscopic Survey (BOSS) to the Wide-Field Infrared Survey Explorer (WISE) to identify quasars with extremely high infrared-to-optical ratios. We identify 65 objects with r(AB)-W4(Vega)>14 mag (i.e., F_nu(22um)/F_nu(r) > ~1000). This sample spans a redshift range of 0.28<z<4.36 and has a bimodal distribution, with peaks at z~0.8 and z~2.5. It includes three z>2.6 objects that are detected in the W4-band but not W1 or W2 (i.e., W1W2-dropouts). The SDSS/BOSS spectra show that the majority of the objects are reddened Type 1 quasars, Type 2 quasars (both at low and high redshift) or objects with deep low-ionization broad absorption lines (BALs) that suppress the observed r-band flux. In addition, we identify a class of Type 1 permitted broad-emission line objects at z~2-3 which are characterized by emission line rest-frame equivalent widths (REWs) of >~150Ang , much larger than those of typical quasars. In particular, 55% (45%) of the non-BAL Type 1s with measurable CIV in our sample have REW(CIV) > 100 (150)Ang, compared to only 5.8% (1.3%) for non-BAL quasars in BOSS. These objects often also have unusual line ratios, such as very high NV/Ly-alpha ratios. These large REWs might be caused by suppressed continuum emission analogous to Type 2 quasars; however, there is no obvious mechanism in standard Unified Models to suppress the continuum without also obscuring the broad emission lines.

Dust in the torus of the AGN unified model

In these proceedings I briefly summarize our current understanding of the dusty torus of the AGN unified model: its structure, composition, and ideas about its origin and evolution. The paper is based on an invited review talk at the 7th Meeting on Cosmic Dust, which covered dust-related topics in areas from comets to debris disks and high-redshift galaxies. This is therefore intended to be an accessible, introductory overview of the torus with some emphasis on the solid-state spectral features observed.

THE ANGULAR CLUSTERING OFWISE-SELECTED ACTIVE GALACTIC NUCLEI: DIFFERENT HALOS FOR OBSCURED AND UNOBSCURED ACTIVE GALACTIC NUCLEI

We calculate the angular correlation function for a sample of 170,000 AGN extracted from the Wide-field Infrared Survey Explorer (WISE) catalog, selected to have red mid-IR colors (W1 - W2 > 0.8) and 4.6 micron flux densities brighter than 0.14 mJy). The sample is expected to be >90% reliable at identifying AGN, and to have a mean redshift of z=1.1. In total, the angular clustering of WISE-AGN is roughly similar to that of optical AGN. We cross-match these objects with the photometric SDSS catalog and distinguish obscured sources with (r - W2) > 6 from bluer, unobscured AGN. Obscured sources present a higher clustering signal than unobscured sources. Since the host galaxy morphologies of obscured AGN are not typical red sequence elliptical galaxies and show disks in many cases, it is unlikely that the increased clustering strength of the obscured population is driven by a host galaxy segregation bias. By using relatively complete redshift distributions from the COSMOS survey, we find obscured sources at mean redshift z=0.9 have a bias of b = 2.9 \pm 0.6 and are hosted in dark matter halos with a typical mass of log(M/M_odot)~13.5. In contrast, unobscured AGN at z~1.1 have a bias of b = 1.6 \pm 0.6 and inhabit halos of log(M/M_odot)~12.4. These findings suggest that obscured AGN inhabit denser environments than unobscured AGN, and are difficult to reconcile with the simplest AGN unification models, where obscuration is driven solely by orientation.

VARIABILITY-BASED ACTIVE GALACTIC NUCLEUS SELECTION USING IMAGE SUBTRACTION IN THE SDSS AND LSST ERA

With upcoming all sky surveys such as LSST poised to generate a deep digital movie of the optical sky, variability-based AGN selection will enable the construction of highly-complete catalogs with minimum contamination. In this study, we generate $g$-band difference images and construct light curves for QSO/AGN candidates listed in SDSS Stripe 82 public catalogs compiled from different methods, including spectroscopy, optical colors, variability, and X-ray detection. Image differencing excels at identifying variable sources embedded in complex or blended emission regions such as Type II AGNs and other low-luminosity AGNs that may be omitted from traditional photometric or spectroscopic catalogs. To separate QSOs/AGNs from other sources using our difference image light curves, we explore several light curve statistics and parameterize optical variability by the characteristic damping timescale ($\tau$) and variability amplitude. By virtue of distinguishable variability parameters of AGNs, we are able to select them with high completeness of 93.4% and efficiency (i.e., purity) of 71.3%. Based on optical variability, we also select highly variable blazar candidates, whose infrared colors are consistent with known blazars. One third of them are also radio detected. With the X-ray selected AGN candidates, we probe the optical variability of X-ray detected optically-extended sources using their difference image light curves for the first time. A combination of optical variability and X-ray detection enables us to select various types of host-dominated AGNs. Contrary to the AGN unification model prediction, two Type II AGN candidates (out of 6) show detectable variability on long-term timescales like typical Type I AGNs. This study will provide a baseline for future optical variability studies of extended sources.

A mid-IR comparative analysis of the Seyfert galaxies NGC 7213 and NGC 1386

New Gemini mid-infrared spectroscopic observations together with Spitzer Space telescope archival data, are used to study the properties of the dusty torus and circumnuclear star formation in the active galaxies NGC 7213 and NGC 1386. Our main conclusions can be summarised as follows. Polycyclic aromatic hydrocarbon (PAH) emission is absent in the T-ReCS nuclear spectra but is ubiquitous in the data from Spitzer at distances above 100 pc. Star formation rates surface densities are estimated from the 12.8 $\mu m$ [Ne{\sc ii}] line strengths leading to values close to 0.1M$_\odot\,\,{\rm yr}^{-1}\,\,{\rm kpc}^{-2}$. Analogous estimates based on photometric fluxes of IRAC's 8 $\mu m$ images are higher by a factor of almost 15, which could be linked to excitation of PAH molecules by older stellar populations. T-ReCS high spatial resolution data reveal silicate absorption at $\lambda$ 9.7 $\mu m$ in the central tens of parsecs of the Seyfert 2 NGC 1386, and silicate emission in the Seyfert 1 galaxy NGC 7213. In the case of NGC 1386 this feature is confined to the inner 20 pc, implying that the silicate might be linked to the putative dusty torus. Finally, by fitting CLUMPY models to the T-ReCS nuclear spectra we estimate the torus physical properties for both galaxies, finding line of sight inclinations consistent with the AGN unified model.

AGN Absorption Linked to Host Galaxies

AbstractMultiwavelength identification of AGN is crucial not only to obtain a more complete census, but also to learn about the physical state of the nuclear activity (obscuration, efficiency, etc.). A panchromatic strategy plays an especially important role when the host galaxies are star-forming. Selecting far-Infrared galaxies at 0.3&lt;z&lt;1, and using AGN tracers in the X-ray, optical spectra, mid-infrared, and radio regimes, we found a twice higher AGN fraction than previous studies, thanks to the combined AGN identification methods and in particular the recent Mass-Excitation (MEx) diagnostic diagram. We furthermore find an intriguing relation between AGN X-ray absorption and the specific star formation rate (sSFR) of the host galaxies, indicating a physical link between X-ray absorption and either the gas fraction or the gas geometry in the hosts. These findings have implications for our current understanding of both the AGN unification model and the nature of the black hole-galaxy connection.

The origin of the IR emission of low-luminosity AGN

AbstractTo test recent suggestions that the infrared emission of low-luminosity AGN arises in a truncated thin accretion disk, we compare recent, high-resolution IR data with published SED model fits that include emission from the truncated disk. We also fit the data with clumpy torus and optically thin dust shell models. These comparisons suggest that dust can better account for the IR emission of the objects in question than can the truncated disk. That optically thin models give a good fit to the data may support a scenario in which the torus of the AGN unified model does not persist in low accretion rate AGN.

GALAXY CLUSTERS AROUND RADIO-LOUD ACTIVE GALACTIC NUCLEI AT 1.3 &lt;z&lt; 3.2 AS SEEN BYSPITZER

We report the first results from the Clusters Around Radio-Loud AGN (CARLA) program, a Cycle 7 and 8 Spitzer Space Telescope snapshot program to investigate the environments of a large sample of obscured and unobscured luminous radio-loud AGN at 1.2 < z < 3.2. These data, obtained for 387 fields, reach 3.6 and 4.5 um depths of [3.6] (AB) = 22.6 and [4.5] (AB) = 22.9 at the 95% completeness level, which is two to three times fainter than L* in this redshift range. By using the color cut [3.6]-[4.5] > -0.1 (AB), which efficiently selects high-redshift (z > 1.3) galaxies of all types, we identify galaxy cluster member candidates in the fields of the radio-loud AGN. The local density of these IRAC-selected sources is compared to the density of similarly selected sources in blank fields. We find that 92% of the radio-loud AGN reside in environments richer than average. The majority (55%) of the radio-loud AGN fields are found to be overdense at a > 2 {\sigma} level; 10% are overdense at a > 5 {\sigma} level. A clear rise in surface density of IRAC-selected sources towards the position of the radio-loud AGN strongly supports an association of the majority of the IRAC-selected sources with the radio-loud AGN. Our results provide solid statistical evidence that radio-loud AGN are likely beacons for finding high-redshift galaxy (proto-)clusters. We investigate how environment depends on AGN type (unobscured radio-loud quasars vs. obscured radio galaxies), radio luminosity and redshift, finding no correlation with either AGN type or radio luminosity. We find a decrease in density with redshift, consistent with galaxy evolution for this uniform, flux-limited survey. These results are consistent with expectations from the orientation-driven AGN unification model, at least for the high radio luminosity regimes considered in this sample.