Radio Galaxies at TeV Energies

: Unlike blazars, radio galaxies have jets that are misaligned relative to our line-of-sight. This misaligned geometry provides us with a unique view of both the jet and super massive black hole. To date, four radio galaxies have been detected at TeV energies with an additional two active galactic nuclei shown to exhibit both radio galaxy and BL Lac-type properties. TeV observations of radio galaxies have revealed these objects to be fascinating, displaying ultra-fast variability and often relatively hard spectral energy distributions. This work aims to provide a review of the current state of radio galaxy observations within the context of very-high-energy γ -ray astronomy, while also highlighting that radio galaxies are excellent targets for multi-wavelength observations. A number of motivations for the continued study of radio galaxies are provided, and these are discussed with a focus on the key observational results, including implications for future observations with next-generation instruments soon to be operational.


Introduction
Active Galactic Nuclei (AGNs) can be defined as the compact regions at the centre of a galaxy in the process of accreting gas through a disk onto the supermassive black holes (SMBHs) residing at their cores. Detected across a broad band of electromagnetic frequencies/wavelengths, AGNs are typically characterised by their high non-stellar luminosity, and these objects can form jets of plasma extending outward from their SMBHs. Strong radio emission from a small fraction of AGNs indicates the presence of a relativistic jet [1], but it is still unknown as to why only some fraction of AGNs produce jets, while others do not. It is believed that a few AGN properties-mass, spin, accretion rate and mode, and spin orientation-along with the environment, determine whether an AGN is jetted or not, but a clear consensus is still to be reached [2][3][4][5][6][7].
Radio galaxies are a subclass of AGNs that are powerful emitters of radio waves, delivering approximately a million times more power (∼10 38 watts) than our own Milky Way Galaxy (∼10 32 watts). Typically, radio galaxies are characterised as having a bright compact core that is coincident with the visible nucleus of the galaxy. However, the radio emission extends beyond this core along jets that propagate outward from the central core to distances of pc or even kpc scales, as seen in NGC 6251 for example [8]. At the termination site of these jets, large plumes are often visible far beyond the central confines of the core. For a comprehensive discussion of relativistic jets from AGNs, which is beyond the scope of this work, please see the excellent review by [9] and all the references therein.
A large number of radio-loud AGN sub-classes are defined based on observational properties seen in specific wavebands, and it is suggested that radio galaxies form the parent population of blazars (AGNs with a relativistic jet pointing toward us). The radio galaxies with an edge-darkened morphology known as FR-Is (Fanaroff and Riley classification, [10]) correspond to BL Lacs observed at large jet viewing angles, while edgebrightened lobe morphology FR-II radio galaxies are associated with flat-spectrum radio quasars (FSRQs) [11]. The distinction between FR-I and FR-II radio galaxies is largely their radio power [12], where FR-Is are weaker, but more common and, hence, more likely to be and for the purpose of this review, these are referred to as radio galaxies (in line with the convention of the VHE community), even though these AGNs may be close to the boundary with blazars. For an earlier review of radio galaxies addressing TeV challenges, please see [35]. Section 2 provides some context by introducing some key characteristics of radio galaxies in terms of blazar studies and long-term γ-ray observations at GeV energies, as well as highlighting some motivations for the continued study of radio galaxies at VHE. Section 3 focuses on the results and findings from observations of radio galaxies detected at VHE, while finally, Section 4 discusses the implications of these key results, as well as addressing some of the future opportunities relating to the observation of radio galaxies with the next-generation of radio and TeV γ-ray observatories.

Radio Galaxies in Context 2.1. Blazar Studies
Understanding how jets are formed, collimated, and accelerated to relativistic speeds is a key goal of blazar studies. The jet emission from blazars and radio galaxies is characterized by a double-peaked, non-thermal spectral energy distribution (SED) [36]. The lowerfrequency peak, which in blazars has a peak frequency ranging from 10 12 to 10 18 Hz, is well described as synchrotron emission. The higher-frequency peak (10 18 Hz-10 27 Hz) is generally attributed to inverse-Compton (IC) emission arising from synchrotron self-Compton (SSC) [37][38][39] or external Compton (EC) processes 1 , or a combination of the two. Typically, the sources with higher-frequency synchrotron peaks have IC peaks in the ∼10-100 GeV energy range. Thus, high-synchrotron-peak (HSP) blazars (ν peak > 10 15 Hz) are the brightest and best-studied at VHE (currently, 53 of the AGNs detected at VHE), although they are the least powerful. Only nine blazars in the current TeV catalogue are low-synchrotron-peak (LSP) objects (ν peak 10 14 Hz), including both FSRQs and BL Lacs.
Unlike blazars, radio galaxy jets are thought to be oriented at larger angles to the line-ofsight ( 10 • ), implying they are are much less Doppler-boosted than blazars, thus exhibiting SEDs similar to LSP blazars [40]. There does appear to be an observed relationship between jet and accretion power [41]. It is also postulated that, in addition to differentiating source classes by accretion modes, a range of jet viewing angles is needed to explain a blazar sequence; low-power jets with inefficient accretion have structured jets [42], whereby the synchrotron peak frequency decreases as the "misalignment" increases and observations probe the lower-speed sheath of the jet. In this scenario, radio galaxies should then cluster at the low end of the synchrotron peak frequency distribution [40].
Keenan et al. [43] extended this idea, with a resulting catalogue of over 2000 radio-loud AGNs with well-sampled synchrotron SEDs (including over 1000 with jet power estimates) showing the distribution of sources on the synchrotron peak-luminosity (L) peak plane is driven by a dichotomy in jet strength where weak jet sources associated with inefficient accretion dominate at log(ν peak ) > 15 and strong jet sources with efficient accretion are primarily found with log(L peak ) > 45 and log(ν peak ) < 15. This shows that the blazar sequence is at least in part due to a selection effect.

Radio Galaxies at GeV Energies
Over the past decade, Fermi-LAT has provided a high-quality view of the γ-ray sky at GeV energies, detecting a number of radio galaxies during the course of its mission. These observations have acted as a pathfinder for TeV observations, for using the current generation of ground-based IACTs, as well as for predicting future observations with CTA, e.g., [44]. In the recently released incremental version of the fourth Fermi-LAT point-source catalogue, which covers the first 12 years of the Fermi-LAT mission (4FGL-DR3) [45], there are a total of 43 non-blazar AGN classified as radio galaxies. The 4FGL-DR3 details for these radio galaxies are highlighted in Tables 1 and 2. Table 1 shows the radio galaxies that have Fermi-LAT-recorded SEDs best fitted with a power-law model, whereas Table 2 shows the radio galaxies that have Fermi-LAT-recorded SEDs best fitted with a log-parabola model. Table 1. The Fermi-LAT 4FGL-DR3 radio galaxies, sorted by right ascension, that exhibit SEDs best fitted with a power-law model. Although Cen A is listed in this table, its SED is actually best fitted with a broken power-law model, as discussed in this review and elsewhere, e.g., [46,47]. The columns include the 4FGL catalogue name, an alternative object name, the type of model best fitted to the SED where PL is for power-law, the spectral index value (Γ ), the photon flux between 1-100 GeV [photons cm −2 s −1 ], the 4FGL-DR3 variability index where an index value >18.48 indicates a <1% chance of being a steady source, and a TeVCat flag indicating whether the object is detected at TeV energies. The photon indices of the 4FGL-DR3 radio galaxies best fitted with a power-law model are found to be in the range between 1.4 and 2.9. In addition, the photon fluxes of all the radio galaxies at GeV energies vary between ∼10 −11 photons cm −2 s −1 and ∼10 −8 photons cm −2 s −1 . As can be seen in Tables 1 and 2, the majority of Fermi-LATdetected radio galaxies have an SED adequately described by a simple power-law model, at least over a long time period such as the 12 years covered by the 4FGL-DR3 catalogue. As defined in the original release of the 4FGL eight-year point source catalogue [48], a source is considered to be variable over yearly binned timescales if its Variability Index (shown in Tables 1 and 2) is larger than 18.48. Based on this measure, ∼72% of the Fermi-LAT radio galaxies are not considered variable and ∼28% are, with NGC 1275 and 3C 120 displaying the strongest variability of all the Fermi-LAT-detected radio galaxies. However, given that many of these radio galaxies are faint (i.e., weak) γ-ray sources, it cannot be concluded for certain that these objects are not variable, but rather that Fermi-LAT cannot measure their variability. Table 2. The Fermi-LAT 4FGL-DR3 radio galaxies, sorted by right ascension, that exhibit SEDs best fitted with a log-parabola model. The columns include the 4FGL catalogue name, an alternative object name, the type of model best fitted to the SED where LP is for log-parabola, the log-parabola spectral index value (α), the log-parabola spectral curvature value (β), the photon flux between 1-100 GeV [photons cm −2 s −1 ], the 4FGL-DR3 variability index where an index value >18.48 indicates a <1% chance of being a steady source, and a TeVCat flag indicating whether the object is detected at TeV energies.

Motivations for Studying Radio Galaxies at TeV Energies
Within the context of TeV γ-ray astronomy, radio galaxies are typically thought of as "misaligned" AGNs with much reduced Doppler boosting factors. The detection of these objects at TeV energies was therefore somewhat of a surprise and opened up many questions about the γ-ray emission processes, γ-ray emission sites, and our general understanding about the nature of these extragalactic objects. Observations of radio galaxies are not overwhelmed with very strong Doppler-boosted emission (as is the case with blazars), and instead, the misaligned geometry provides us with a unique view of the AGN's jet and SMBH. Despite the small sample size, the detection of γ-rays from radio galaxies provides a number of motivations for future observations and studying these objects further. The following is a small selection of the outstanding questions relating to radio galaxies at TeV energies: • Do all TeV-detected radio galaxies exist in dense clumpy environments, for example? What role does the environment play in decelerating radio galaxy jets, for example, and does this interaction lead to γ-ray emission and/or γ-ray orphan flares? (See [49] and the references therein.) • Do TeV-detected radio galaxies exhibit any common kinematic characteristics? How do radio galaxies accelerate particles in their jets, and can regions with different apparent velocity profiles (over all jet scales) be correlated with γ-ray emission? What mechanisms give rise to features in radio galaxy jets, e.g., internal shocks, and are these responsible for gamma-ray emission? • Is it right to classify the TeV-detected radio galaxies as misaligned blazars, or are these objects simply blazar-like, with much smaller jet inclination angles to the line-of-sight than currently derived or assumed? • In TeV-detected radio galaxies, where does gamma-ray emission occur? Is it close to the black hole in the inner-jet region or at kpc scales further down the jet? Are there multiple γ-ray emission sites and components, for example, both jet and large-scale emission from the up-scatter of external photons? • Are TeV-detected radio galaxies best described by a singular leptonic, hadronic, or mixed emission model? Which process dominates? • Are TeV-detected radio galaxies cosmic ray sources? If γ-ray emission from radio galaxies can be correlated with neutrino observations, can the population of as-yet-undetected radio galaxies in our universe make a substantial contribution towards the cosmic ray background?
Section 4 attempts to address some of the themes covered in these questions, but in order to resolve any of them, it is important to state that TeV observations alone cannot fully answer all of these questions. Instead, multi-wavelength observations from radio to optical to γ-rays are all needed, and we live in an exciting time with many new large and improved observatories on the horizon, all capable of making major contributions toward a better understanding of AGN science. Taking into account unification considerations [50], there is a motivation to explore these objects deeper and with higher precision. For example, to construct better theoretical emission models for radio galaxies, there is a need for wellsampled broadband spectra of these objects in different states. In order to do this at the highest energies, deeper observations are required. In addition, TeV observations provide us with the best opportunity to probe the fastest variability at the highest energies, as this will enable us to gain a better understanding of the nature of these objects and tell us more about the location of gamma-ray emission sites, for example whether the observed gamma-ray emission originates close to the core or further down the jet. Furthermore, it is suggested that given radio galaxies are now detected at both GeV and TeV energies, their large population likely provides a significant contribution toward the cosmic rays that permeate our universe. Therefore, multi-messenger observations of radio galaxies may provide us with insights into their overall contribution toward the cosmic ray background.

Results
As already mentioned above, four radio galaxies have been detected with groundbased IACTs. The discovery of these objects with the current generation of TeV instruments is exciting as this paves the way for some promising observations with future IACTs to provide new insights into the mechanisms driving gamma-ray emission. In this section, we highlight the key findings from the observations of the four TeV-detected radio galaxies including the additional two non-blazar AGN that display both radio galaxy and blazar characteristics. Table 3 highlights the key characteristics of these objects. Table 3. Highlighted here are the key characteristics of the 6 non-blazar AGNs detected at TeV energies. The table is sorted by redshift in ascending order. The column details are as follows: the common name given to the AGN (Name), an alternative name for the AGN (Alt. Name), redshift (z), the jet inclination angle to the line-of-sight (θ i (degrees)), jet power (P jet 1.4 GHz (W/Hz)), radio luminosity at 5 GHz (L 5 GHz (ergs/s)), synchrotron peak frequency (ν peak (Hz)), VHE discovery date (dd VHE (year-month)), and finally, a column containing references for each object (Refs.).

Name
Alt. Name z θ i log P

NGC 1275
NGC 1275 (also commonly known as 3C 84) is the brightest cluster galaxy at the centre of the Perseus cluster of galaxies. It is located at ∼75 Mpc and was detected by both MAGIC [61] and VERITAS [62]. Radio observations of NGC 1275 show a sub-parsec-scale radio jet oriented between ∼30 • and 60 • with possible bending in the jet. Figure 1 shows NGC 1275 at 43 GHz (left panel) and 22 GHz (right panel) radio frequencies. Here, we see the sub-parsec-scale jet appears to have a complex structure and is dominated by two features 2 : the core (A0) and a knot (C2) [63]. There does not appear to be evidence of superluminal apparent speeds in the parsec-scale jet; indeed; C2 moves at subluminal speeds of ∼0.2c relative to the core. The left panel of Figure 1 also highlights a number of other knots (C1...C8), but the one with the highest recorded proper motion is C6 at ∼0.36 c. Knot C3 appears to move upstream, while knot C4 appears to be at a significantly different, yet stable position angle. Recent observations suggest the component C3 is a remnant of previous emission [49,64].
The knot C2 emerged following a radio outburst in 2005 located in the central 1pc region [65]. After this outburst, limb-brightening of the jet was discovered, indicating that the transverse structure of the jet had changed. This morphological change appears to have coincided with an increased γ-ray flux level at GeV energies [66,67]. Furthermore, the radio emission from knot C2 was significantly polarized, and there appeared to be a location change from southwest to southeast. It is suggested that this movement and enhanced polarization indicated that the jet jittered in an inhomogeneous ambient medium [68]. A ∼0.4 mas directional flip of the knot C2 was confirmed with radio observations taken in 2015 and 2016, providing the first evidence of jet interaction with the environment [49]. It is believed the jet of NGC 1275 interacted/collided with a compact dense cloud, which led to a year-long frustration of the jet before breaking out in 2018. Subsequently, the knot C2 continued to move southward until eventually becoming distorted and returning NGC 1275 to a darkened lobe morphology state [64].
The high-resolution radio observations shown in the right panel of Figure 1 provide evidence of a limb-brightened jet. The large intensity ratio between the bright outer layer and the dim central region can be explained by transverse velocity structures and supports γ-ray emission scenarios via spine-sheath-structured jet models [42,69]. Sub-arcsecond optical images of NGC 1275 reveal a possible recent merger event [70], which is interesting because galaxy mergers are postulated to be potential sources of cosmic rays, neutrinos, and γ-rays [71].
NGC 1275 exhibits significant flux variability with flux doubling times of around 10 h observed, suggesting γ-ray emission originates from a compact region. In its quiescent state, the spectral shape of NGC 1275, determined from MAGIC observations, appears to follow a log-parabola (or power-law with sub-exponential cut-off) with a somewhat steep spectral index of Γ∼−4.1 and no excess γ-ray photons above ∼650 GeV detected. In January 2017, significant increases in the flux levels seen from NGC 1275 were reported by MAGIC (∼1.5-times the flux from the Crab Nebula) [72] and VERITAS (∼65% the flux from the Crab Nebula) [73]. Significant variability was observed in the day-to-day light curve of NGC 1275 with a fastest flux-doubling timescale of (611 ± 101) min. seen [74]. Such a flux-doubling timescale implies that the emission region responsible for VHE γ-rays is larger than the event horizon light crossing time, but it still represents a compact region of the jet as opposed to a large-scale emission zone.
From these MAGIC observations of NGC 1275 during a flaring state, the MAGIC Collaboration et al. [74] attempted to model the observed SED assuming a classical SSC leptonic scenario. However, this was not possible due to the very-fast-variability scale constraining the size of emission zone at one to two orders of magnitude below that typically used in the SSC scenario adopted by the authors. In attempting to model the γ-ray emission from NGC 1275 by assuming an SSC leptonic scenario, and as hinted elsewhere, a higher Doppler factor would be needed to circumvent absorption, thus necessitating a jet inclination angle (to the line-of-sight) that is at odds with the jet inclination angle determined from radio observations. The parsec-scale jet shows a complex structure and is dominated by two features: the core (A0) and a knot (C2). Knot C6 has the highest recorded proper motion of ∼0.36 c. The right panel taken from [67] shows the 22 GHz space-VLBI image of NGC 1275's jet from the core to ∼1 pc. These high-resolution radio observations provide evidence of a limb-brightened jet. The large intensity ratio between the bright outer layer and the dim central region can be explained by transverse velocity structures and supports γ-ray emission scenarios via spine-sheath-structured jet models.

Cen A
Cen A (NGC 5128) is the nearest (d∼3.8 Mpc [75]) and one of the most well-studied extragalactic objects over a wide frequency domain. A misaligned FR-I radio galaxy, Cen A is often considered the prototypical Fanaroff-Riley FR-I classification. Cen A exhibits a complex and extended radio morphology with two giant lobes, as shown by the radio contours (black lines) in Figure 2. The latest sub-parsec-scale VLBI observations of Cen A obtained with the Event Horizon Telescope (EHT) at 228 GHz reveal a highly collimated, asymmetrically edge-brightened jet along with a faint counter-jet [76]. This strong limb-brightening seen in the sub-parsec-scale jet of Cen A is similar to that seen in other TeV-detected radio galaxies M 87 and NGC 1275, which suggests the dominant emission originates from the jet sheath and could potentially indicate that some component of γ-ray emission from Cen A takes place via spine-sheath-structured jet models.
Faint TeV γ-ray emission from Cen A was detected by H.E.S.S. after approximately ∼120 h of observations. The TeV emission was detected with a statistical significance of 5σ above the background coincident within a region including both the radio core and the inner kiloparsec jet. The integrated flux above ∼250 GeV was estimated to be ∼0.8% that of the flux from the Crab Nebula. A power-law shape was best fitted to the spectral energy distribution with a photon index of 2.7 ± 0.5 stat ± 0.2 sys , and no significant flux variability was detected [77].
Using the first four years of Fermi-LAT data, Sahakyan et al. [46] first reported evidence (albeit not statistically significant at < 3σ) for an observed hardening of Cen A's energy spectrum. They reported a break in the spectrum at E b ∼4 GeV with a spectral index of Γ E<4 GeV = 2.74 ± 0.03 below and Γ E>4 GeV = 2.09 ± 0.20 above the break. Sahakyan et al. [46] suggested that the best-fitted broken power-law model could be caused by the contribution of an additional high-energy component beyond common synchrotron self-Compton jet emission. With eight years worth of Fermi-LAT data, Brown et al. [47] confirmed the existence of a statistically significant spectral hardening in Cen A's spectra, a likelihood-ratio test with TS = 28.6 ≈ 5.3σ demonstrating that Cen A's spectra are best described by a broken power-law model compared to a standard power-law model. The spectral index values of their study agreed (within the reported errors) with those reported by Sahakyan et al. [46] both above and below the break energy. In addition, Brown et al. [47] found no evidence for significant flux variability from Cen A in either the low-or highenergy spectral components when considering the full eight-year dataset. However, they did note that their results showed some flux variability (at the 3.3σ level) in the low-energy spectral component during the first 3.5 years of Fermi-LAT data. The lack of variability, in particular for the high-energy spectral component, rules out the possibility that the spectral feature of Cen A is a result of jet-induced leptonic processes. Instead, Brown et al. [47] suggested that the feature seen in the combined GeV-TeV SED is consistent with a spike in the dark matter halo profile, but also postulated that a population of unresolved millisecond pulsars or some other population of energetic particles could be responsible for the spectral feature.
Single-zone SSC models do not adequately describe Cen A's observed broadband spectral energy distribution. As pointed out in [78], a variety of different scenarios exist, and none of them mutually exclusive, for the physical origin of the observed spectral energy break. To avoid going into great detail for each scenario here, the appropriate references are included again: • Magnetospheric pulsar-like scenarios [79,80]; • Inner-jet models involving a number of leptonic SSC-emitting components that propagate at different angles relative to the line-of-sight [39], inverse-Compton emission from structured jets as discussed elsewhere in this article [42], photo-meson decay processes [81][82][83][84], lepto-hadronic emission models [85,86], and γ-ray-induced pair cascades from different regions [87][88][89]; • Extended emission scenarios such as the interaction of energetic protons with ambient matter (proton-proton interactions) at kiloparsec scales [46], millisecond pulsar contributions as already highlighted above [47], the inverse-Compton upscattering of photons on kiloparsec scales [90], or host galaxy starlight [91]; • The self-annihilation of dark matter, as already highlighted above [47].
With the current generation of IACTs, the jets of Cen A cannot be resolved. Therefore, distinguishing between these scenarios is currently difficult, and the observed γ-ray emission has tended to be associated with a region close to the black hole. Recent observations of Cen A by H.E.S.S. have focused on looking for extended emission and by comparing different models using a test statistic, and they found a strong preference for an elliptical Gaussian model (TS = 19.4) compared to a radially symmetric Gaussian model (TS = 6.1) [92]. This result, shown in Figure 2, implies that a major part of the VHE γ-ray emission originates on large scales at great distances from the SMBH.

M 87
M 87 (3C 274, NGC 4486) is a radio galaxy located in the Virgo cluster at a distance of d ∼ 16.7 Mpc [93], making this the second closest radio galaxy after Cen A, but the closest kiloparsec-scale radio-optical jet. The jet of M 87 is misaligned with an angle i ∼ 15 • -25 • to the line-of-sight. M 87 is classified as a low-excitation weak-power FR-I-type radio galaxy displaying a radiatively inefficient accretion mode.
HST observations of M 87 provide evidence for superluminal speeds up to ∼6c out to ∼0.5 kpc [94], and radio observations of the feature known as HST-1 (located at a deprojected distance > 100 pc) show superluminal speeds of ∼4c [95][96][97]. However, Kovalev et al. [98] found no evidence for motions faster than 0.07c within the inner ∼1.6 pc of the jet, and using radio images from the European VLBI Network, Asada et al. [99] suggested that the jet of M 87 is gradually accelerated from these subluminal speeds seen at the milliarcsecond scales of the inner-jet out to the superluminal speeds seen at HST-1. Mertens et al. [53] raised concerns with this scenario including the implied need for a jet-to-counter-jet intensity ratio of 1 when the jet of M 87 is essentially one-sided. Furthermore, the observed limbbrightened morphology of the jet strengthens the argument that the subluminal apparent speeds reported for the inner-jet point to a pattern speed or a wind from an outer sheath in the jet.
Using the 10 m Whipple Cherenkov telescope, in the years 2000/2001, early efforts were made by the VERITAS collaboration to detect VHE γ-ray emission from M 87 [100]. The question of whether M 87 emitted TeV γ-rays was again raised by Aharonian et al. [101] after finding excess photons (with a statistical significance of 4σ) at energies E > 730 GeV coincident with the position of the radio galaxy using data taken with HEGRA [102]. VHE γ-ray emission from M 87 was detected by H.E.S.S. with a statistical significance of 13σ [103] and thereby establishing M 87 as the first non-blazar extragalactic TeV source. In addition, the observations providing this first detection of M 87 at TeV energies also provided evidence for statistically significant spectral variability on the timescale of days. This early detection strongly constrained the size of the γ-ray-emitting region and challenged most scenarios of VHE γ-ray production theorised at that time for extragalactic sources. The VERITAS collaboration also reported the detection of VHE γ-rays from M 87 showing a correlation between the VHE γ-ray flux and M 87's X-ray emission detected by the ASM RXTE observatory [104]. Observations by MAGIC also confirmed the detection of M 87 as a TeV source [105]. Again, these early observations seemed to suggest that the γ-ray emission was most likely associated with the core of the radio galaxy as opposed to a knot in the jet, for example. However, a detailed multi-wavelength study of M 87 by Abramowski et al. [106] showed two of three detected flares exhibited rapid variability, suggesting γ-ray emission close to the core and the remaining flare coincident with the brightening of an HST-1 knot in the jet. This study used both the historical VHE observations, as well as new VHE data obtained during several coordinated monitoring campaigns between 2008 and 2010. Figure 3 shows a composite image of M 87 viewed across a broadband spectrum from radio to γ-ray frequencies. The image demonstrates the huge scales between images of the black hole at the heart of M 87 versus images encompassing the whole jet and extent of the radio galaxy itself. Figure 3 should also provide the reader with some insight into the challenges faced by astronomers when trying to understand and model the mechanisms driving the non-thermal emission from radio at one end of the electromagnetic spectrum to γ-rays at the other. Current instruments, at γ-ray energies, are not able to spatially resolve emission from the core and the jet.

3C 264
The fourth radio galaxy detected at VHE energies is 3C 264. This was detected by VERITAS [20] between February 2017 and May 2019. At a distance of ∼93 Mpc, 3C 264 is the most distant radio galaxy detected at TeV energies, but is also one of the brightest and closest radio-optical jets known [108]. The excess γ-rays coincident with the location of 3C 264 were detected with a statistical significance of 7.8σ above the background. The detection of 3C 264 came from deep observations comprising some 57 h of data and is believed to have occurred during an elevated flux state during the 2018 observations [109]. Like M 87, 3C 264 exhibits superluminal motions on kiloparsec scales. When comparing the jet structure, 3C 264 closely resembles M 87. Both galaxies have one-sided FR-I-type jets with a kinetic luminosity of ∼10 43.8 erg s −1 [40]. In addition, their jets display quite similar morphological traits, for example they both have multiple knots. Indeed, a one-to-one correspondence between radio and optical features was found in the jet of 3C 264 [110]. Like the other TeV-detected radio galaxies, 3C 264 exhibits a limb-brightened jet [111] indicating the possibility of γ-ray emission via spine-sheath-structured jet models. Using HST observations, Meyer et al. [112] found apparent speeds as high as β app = (7.0 ± 0.8)c, exceeding the apparent speeds measured in the jet of M 87 highlighted in Section 3.3 above. Thus, 3C 264 could be considered as a more distant analogue of M 87, based on both its beamed and unbeamed radio emission, as well as its kinematic profile. 3C 264 is most likely observed at a smaller angle to the line-of-sight however.
The VHE flux observed by VERITAS during 2017-2019 shows slight variability, which suggests the detection did not result from a strong flare, but instead from a modestly enhanced flux state over a period of time. At GeV energies, 3C 264 shows some variability on six-monthly timescales, but compared to other Fermi-LAT-detected radio galaxies such as NGC 1275 or even 3C 120, it does not exhibit strong variability on these timescales [113]. At VHE energies, Archer et al. [109] reported 3C 264 showing variability on annual and monthly timescales. When comparing HST images of the jet in 3C 264 between 1994 and 2014, Meyer et al. [112] found two superluminal knots to be on a collision course. Indeed, the authors reported that these two knots began interacting during 2014 with the collision still ongoing today. Such an ongoing collision between two knots could be responsible for the observed VHE γ-ray emission. However, Archer et al. [109] concluded that the unresolved core was the most likely source of the enhanced VHE flux, particularly as there is no evidence for any notable changes in any of the contemporaneous high-resolution Chandra or HST images. The authors also stressed that no large changes in the core flux were seen at lower frequencies. Figure 4 shows the broadband SED of 3C 264 compared to that of M 87. Archer et al. [109] were able to model the SED of 3C 264 successfully using a single-zone SSC model, although they did point out that a better fit could be achieved using a multi-zone model. However, from simple visual inspection, it is immediately noticeable that 3C 264 has a relatively highfrequency synchrotron peak near the X-ray frequency band (∼10 17 Hz), which appears to be somewhat unusual for a radio galaxy, as discussed in Section 2. In another study of 3C 264, the broadband SED and jet power were found to be reasonably well described by simple one-zone SSC leptonic models, but only when the high-energy component was associated with a second emitting region located at the end of the acceleration region (∼11 pc) in the spine or sheath of the jet. Furthermore, when the SED was modelled assuming synchrotron emission from the core, the emission at higher energies was significantly underpredicted [111]. , which for a radio galaxy appears consistent with that expected according to the blazar sequence, the SED for 3C 264 shows a relatively high-frequency synchrotron peak near the X-ray frequency band (∼10 17 Hz).

PKS 0625-35 and IC 310
The two remaining non-blazar TeV-detected AGNs, PKS 0625-35 and IC 310, display both radio galaxy and blazar characteristics. TeV γ-rays (above an energy threshold of ∼250 GeV) from PKS 0625-35 were first detected in 2012 by H.E.S.S. with a statistical significance of ∼6σ above the background [114]. Located at a distance d ≈ 235 Mpc (z = 0.055), PKS 0625-35 is much further away relative to other TeV-detected radio galaxies. Radio observations revealed a one-sided jet that extends to a distance of ∼143pc and exhibits superluminal apparent velocities of ∼3c [115][116][117]. Evidence for blazar characteristics was found in optical observations where the presence and intensity of the [OIII] emission line was consistent with those of other BL Lac objects [118]. In further work by the H. E. S. S. Collaboration et al. [119], it was suggested that PKS 0625-35 displays other unusual characteristics for a radio galaxy such as the synchrotron peak lying between ultraviolet and X-ray frequencies, which is typical for BL Lac objects. Furthermore, X-ray variability studies constrain the maximum possible misalignment to θ ≈ 15 • , and the authors therefore argued that the object be classified as a blazar or moderately misaligned blazar.
IC 310 is located in the Perseus cluster along with NGC 1275 and is at a distance of d ≈ 80 Mpc (z = 0.019) from Earth. MAGIC made a serendipitous discovery of this headtail radio galaxy during 2009/2010 while observing NGC 1275 located in the same field of view (∼0.6 • away). An initial detection corresponding to a statistical significance of 7.6σ above the background was obtained from 20.6 h of stereoscopic observations [120,121]. The observed energy spectrum as seen by MAGIC was relatively hard with a spectral index Γ = −2.00 ± 0.14. MAGIC estimated the average flux above 300 GeV to be (3.1 ± 0.5) × 10 −12 cm −2 s −1 , which corresponds to ∼2.5% of the Crab Nebula flux [121]. An extensive multi-wavelength campaign observing IC 310 was conducted between 2012 and 2013. After data quality cuts, MAGIC acquired an additional ∼39 h on IC 310 including ∼3.7 h when MAGIC detected a bright flare. Thus, MAGIC was able to separate out this flare data to estimate energy spectra for two distinct states; during the low state, the energy spectrum is best described by a standard powerlaw model with normalisation at 1 TeV of N 0 = (3.12 ± 0.91) × 10 −13 TeV −1 cm −2 s −1 and spectral index Γ ≈ (2.36 ± 0.30), and during the bright flare, the energy spectrum reached a peak normalisation at 1 TeV of N 0 = (44.2 ± 2.1) × 10 −12 TeV −1 cm −2 s −1 with spectral index Γ ≈ (1.51 ± 0.06) [56].
Therefore, at TeV energies, the spectral index appears to display harder-when-brighter behaviour. Contemporaneous soft X-ray data taken during this multi-wavelength campaign appear to show a similar harder-when-brighter trend, and the authors concluded that this relationship is consistent with an SSC emission scenario, successfully fitting a single-zone SSC model to the broadband SED during the non-flaring state. Similarly, at GeV energies, IC 310 appears highly variable in both flux and spectral state [122].
As mentioned above, the classification of IC 310 is still the subject of some debate. Originally classified as a narrow-angle tail (NAT) radio galaxy, early γ-ray detections of IC 310 categorized this object as an AGN of unknown type before being assigned as an HBL in the Fermi-LAT 2nd AGN catalogue. However, the discovery of a one-sided core-jet structure with blazar-like, beamed radio emission oriented along the kiloparsec tail was not consistent with the NAT classification. Instead, IC 310 appears to be a low-luminosity FR-I radio galaxy that is moderately misaligned i ∼ 10 • -20 • still revealing its BL Lac-type central engine [123]. Moreover, Figure 5 shows new low-frequency radio observations that provide a resolved view of the bent jets and have led to a new estimation of the misalignment to be i ∼ 8 • . This evidence means that the γ-rays most likely originate from a small region close to the SMBH and are beamed toward us along the line-of-sight, within a region of the jet before it is bent into the narrow tail-like structure [55]. Furthermore, rapid minute-scale variability seen at TeV energies is consistent with such an interpretation. Figure 5. Taken from [55], these high-resolution VLA images of IC 310 show a zoomed-in region at the base of the jets. The left panel image from 320-470 MHz radio observations shows two distinct jets that are collimated, narrow, and oriented in the same general direction. The right panel image taken at higher radio frequencies of 4-8 GHz shows the same jet features, but closer to the nucleus, these are less distinguishable compared to the low-frequency observations. γ-ray emission most likely originates close to the SMBH in a region oriented close to our line-of-sight, where Doppler boosting takes place and before the jets are bent.

Radio Galaxies with Water Cherenkov Detectors
In addition to next-generation IACTs such as CTA, other techniques for exploring the VHE space, such as using water Cherenkov detectors (WCDs), were used in the past such as the VHE surveys conducted by Milagro [124,125] and ARGO [126]. At present, the High Altitude Water Cherenkov Observatory (HAWC) [127] located in central Mexico and a hybrid array called the Large High Altitude Air Shower Observatory (LHAASO) [128] in Sichuan province of China are operational, and another future observatory using WCDs called the Southern Wide-field Gamma-ray Observatory (SWGO) [129] (formerly the Southern Gamma-ray Survey Observatory) is also proposed. Observatories employing WCDs typically have a larger angular resolution compared to IACTs. However, compared to IACTs, their sensitivity for larger extended sources is better, and they are able to operate both day and night under all weather conditions with a large instantaneous field of view > 2 sr. HAWC's angular resolution (defined as the 68% photon containment radius) is both energy-and zenith-angle-dependent and ranges from 0.1 • -1 • [127]. Furthermore, WCDs such as HAWC (and the proposed SWGO) are more sensitive to sources that transit directly overhead and also to sources with a hard energy spectrum. Thus, HAWC's energy threshold for detection is source-dependent. Naturally, this imposes a selection bias on the observable radio galaxies, but may be extremely useful in detecting any steady VHE γ-ray emission from radio galaxies, such as Cen A, that display spectral hardening in the GeV/TeV band and also display extended emission. Unfortunately, Cen A is not visible to HAWC, so undoubtedly will be a priority target for SWGO.
The relative closeness of radio galaxies means smaller γ-ray photon attenuation on the extragalactic background light (EBL), making these potentially great targets for studying the far-infrared band of the EBL above ∼10 TeV. In the latest results from the HAWC survey of active galaxies, six radio galaxies with favourable transits for HAWC were considered: NGC 1218, IC 310, NGC 1275, 4C +39.12, 3C 264, and M 87. At GeV energies, all of these radio galaxies have hard spectra; in the case of NGC 1275, this is true for active states. From live data spanning a period of approximately 4.5 years, HAWC sees only marginal evidence of γ-ray emission from M 87 with a statistical significance of ∼3.6σ at E > 0.5 TeV [130]. Similarly, M 87, NGC 1275, and IC 310 are all nearby sources within the LHAASO field of view [131], and since the observatory was only completed in 2021, observational findings for these sources will undoubtedly be published in the near future. Finally, even if WCDs require a longer integration time to reach the accepted statistical threshold to claim a detection, these types of detectors have a part to play in improving our understanding of radio galaxies, particularly if SWGO is able to bring advances that greatly improve the achievable angular resolution of current WCDs.

Discussion
The above sections have provided some context for radio galaxies within γ-ray astronomy by highlighting some key characteristics of these objects relating to wider blazar studies, as well as highlighting key findings from multi-wavelength observations of radio galaxies detected at VHE. This section will address some of the themes highlighted in Section 2.3 by discussing the principal findings from these observations and what implications if any these have for future observations of these fascinating extragalactic objects.

Radio Galaxy Environment
The first and probably most obvious finding from the TeV observations is that all four radio galaxies detected at VHE energies have low-power FR-I-type jets. This may seem counter-intuitive for accelerating particles to VHE, but equally, this may tell us that FR-I-type jets live in denser environments and that it is the interaction of the jet with these dense environments, rather than the intrinsic power of the jet, that provides the conditions for accelerating particles to VHE. For example, it is suggested that changes in the direction and flux of NGC 1275's bright C2 component (as shown in Figure 1) during August-September of 2015 are consistent with hydrodynamical simulations of jets in a clumpy ambient medium [49]. Furthermore, observations of 3C 264 show that the radio galaxy is located within a large region exhibiting low-surface-brightness emission, which could indicate that the radio galaxy resides in a dense environment [111]. Therefore, if radio monitoring reveals potential jet interactions/collisions with dense clumps in the environment (as in the case of NGC 1275 [49,64]), this could indicate an increased likelihood for enhanced γ-ray emission via internal shock models, leading to flares or possibly even orphan flares at TeV energies.
At the time of writing, no high-power FR-II-type radio galaxies are listed on TeV-Cat [22], but it is reasonable to assume that the current generation of ground-based IACTs have attempted to detect these objects without success. Perhaps this is not completely surprising given that FR-IIs are systematically farther away [132,133] and their VHE γ-ray photons are more likely to be absorbed by the EBL, coupled with the fact that the current generation of IACTs simply does not possess the necessary sensitivity to detect γ-rays produced at their jet termination sites. If CTA is able to detect FR-II radio galaxies, this will expand our understanding of these objects further and provide better insight into the emission mechanisms at play, the emission locations, and perhaps, even their environments.

Kinematic Studies and Classification Considerations
Another obvious question one might ask is whether the TeV-detected radio galaxies are indeed radio galaxies in the "true sense", particularly as most of these objects exhibit properties such as a one-sided jet, fast variability, and superluminal motion, all characteristics of blazars, but less extreme in these cases. This of course is an argument of classification, and as already highlighted above, classification is often a product of the observed wavelength/frequency. Given that TeV emission from these extragalactic AGNs requires a relatively high Doppler factor, any blazars with jets not aligned closely to the line-of-sight are designated misaligned blazars. Thus, in the context of VHE γ-ray observations, this includes both "true" radio galaxies, as well as quasars, but also less extreme blazars falling near the radio galaxy boundary. There is no absolute angle used to precisely define the exact boundary of when a blazar becomes a radio galaxy; instead, if the jet is inclined with respect to the line-of-sight at angles θ 10 • , then it is generally accepted that these are not BL Lacs, nor FSRQs in the strict sense. This corresponds to the critical angle used to distinguish the division between BL Lacs and FR-I-type radio galaxies [134].
In the case of IC 310, perhaps our understanding of the true nature of these objects is not complete. Perhaps γ-rays are emitted in regions that are closely aligned to our line-of-sight where large Doppler factors are attainable. What this result shows us is that high-resolution parsec-and kiloparsec-scale radio images are needed (in particular for radio galaxies that display rapid VHE variability) in order to better determine line-of-sight angles. Indeed, blazar multi-wavelength monitoring programs such as the MOJAVE survey program [135], the BEAM-ME program (successor to the VLBA-BU-BLAZAR program) [136], and the TANAMI program [59] are crucial as we enter a new era of observations with new instruments coming online soon, such as the James Webb Space Telescope (JWST) [137], the Square Kilometer Array (SKA) [138], and CTA [26], or planned, such as the Next-Generation Very Large Array (ngVLA) [139], as well as successors of Chandra [140] and Fermi-LAT.
Programs such as those highlighted above include many of the Fermi-LAT-detected radio galaxies and allow regular monitoring of radio amplitudes that enable correlation studies between radio and γ-ray frequencies. In addition, these monitoring programs can provide polarisation data that underline the magnetic field structure of the radio-emitting regions. All of this is crucial information that is needed when attempting to associate γ-ray emission with a particular region or emission process for example. Additionally, given the relatively poor angular resolution of γ-ray instruments, both from space and ground, very-long-baseline interferometry (VLBI) monitoring programs play a vital role in charting the evolution of parsec-scale radio emission, which again enables the study of any relationship between radio and γ-ray emission.
A further point for consideration is that we might only detect these four radio galaxies at TeV energies due to their proximity, i.e., even though they are misaligned, the combination of their proximity coupled with strong-enough Doppler boosting means they are detectable at TeV energies. However, it should be noted that at least two other Fermi-LAT-detected radio galaxies reside closer than 3C 264, and these are not detected at TeV energies. Furthermore, even though PKS 0625-35 displays blazar characteristics, it appears misaligned nonetheless, and there are many radio galaxies undetected at TeV energies closer than this object.
One last observation is that all four of the TeV-detected radio galaxies exhibit limbbrightened jets. Perhaps this is an important feature that could indicate an increased likelihood for gamma-ray emission via spine-sheath-structured jet models, where a smallenough beaming angle is maintained to provide sufficient Doppler boosting.

Emission Sites and Variability
The detection of synchrotron X-ray emission from the jets of Cen A and M 87 implies the presence of electrons at and above TeV energies, and it is possible that the TeV γ-ray emission is produced (or partly produced) by inverse-Compton scattering of photons by the high-energy electrons producing the synchrotron emission on kiloparsec scales [90]. Indeed, the detection of extended gamma-ray emission in Cen A with Fermi-LAT [141] coupled with a lack of evidence for strong variability means that an extended leptonic origin for its emission cannot be ruled out. As highlighted in Section 3.2, the latest H.E.S.S. observations of Cen A suggest that a major part of the VHE emission arises on large scales [92]. Moreover, the H.E.S.S. findings are consistent with (and provide the first evidence of) VHE γ-ray emission resulting from the inverse-Compton scattering of low-energy photons by ultrarelativistic electrons in the jet. This inverse-Compton emission potentially contributes toward the spectral hardening (discussed above) observed in Cen A and also suggests that radio galaxy jets are potentially a steady source of VHE gamma-rays.
However, in the case of M 87, the picture is not so clear. As highlighted above, two out of three TeV-detected flares exhibit rapid variability with timescales on the order of approximately half a day, which suggests that γ-ray emission originates from a region close to the core. Yet, the remaining TeV-detected flare was correlated (using multi-wavelength data) with the HST-1 knot located downstream in the jet of M 87 [106].
Fermi-LAT observations of M 87 found flux variability on six-monthly timescales at energies below E < 10 GeV [142]. Similarly, in a study of 26 radio galaxies at GeV energies using Fermi-LAT data [113], only 7 were seen to exhibit any statistically significant flux variability on six-monthly timescales. Of these seven, two (3C 264 and PKS 0625+35) were detected at TeV energies. On the other hand, NGC 1275 has displayed rapid variability at GeV energies during flares with flux-doubling timescales of 1.21 ± 0.22 h and 8.03 ± 1.24 h reported [143]. Such flux-doubling timescales are more usual for blazars rather than radio galaxies, but as already alluded to before, radio galaxies are weak γ-ray sources and the sensitivity of the current space-and ground-based γ-ray observatories limits us to seeing variability in the brightest radio galaxies.

Emission Mechanisms and Modelling
When modelling the broadband spectral energy distributions of blazars and indeed radio galaxies, often a single-zone synchrotron self-Compton (SSC) jet model is used to provide a reasonable description of the emission mechanisms. However, for radio galaxies such as Cen A with very modest Doppler boosting, it is clear that the single-zone SSC models do not adequately describe the VHE γ-ray emission [39,144,145]. Yet, for other radio galaxies such as 3C 264, its broadband emission spectrum can be described by a single-zone SSC model, as shown in Figure 4, but is probably better fit using multi-zone models [109].
When modelling the broadband SEDs of radio galaxies, the jet misalignment value used (or assumed) has a profound impact on the model used to describe the emission. Assuming TeV-detected radio galaxies displaying rapid variability are truly misaligned, then exotic scenarios such as "mini-jets in jets" [146] or "black-hole lightning" [147] are necessary in order to produce VHE γ-rays. In-depth discussion of all the models applicable to radio galaxies along with the intricacies of applying those models to SEDs is beyond the scope of this work. For a comprehensive review of blazar modelling in the context of TeV extragalactic astronomy please, see [16], and for an earlier review of radio galaxies at VHE energies that includes discussion of the models relevant to radio galaxies, please see [148]. Instead, what is important to take away is that radio galaxies are excellent targets for multi-wavelength observations, which provide us with an opportunity for correlation studies across wavebands in order to single out a preferred model (or set of models) that explain γ-ray emission from these objects.

Radio Galaxies as Cosmic Ray Sources
It is well established that cosmic ray particles permeating our galaxy cannot be traced back to their origins (except perhaps at the highest energies) as they are deflected by magnetic fields. We know that the origin of cosmic rays above a few EeV must be extragalactic, but as of yet, no specific type of source is positively identified as the primary cosmic ray contributor [149,150]. However, there is growing evidence in the Pierre Auger Observatory data (independent of the inconsistency arguments highlighted regarding the connection between ultra-high-energy cosmic rays (UHECRs) and AGNs) that Cen A is a potential source of UHECRs and neutrinos [150][151][152][153][154]. Indeed, the work by Koers and Tinyakov [155] highlights the possibility of Cen A being an exceptionally efficient neutrino source. This is important because if radio galaxies like Cen A are shown to be neutrino sources and can be correlated to TeV γ-ray emission in these objects, this provides very strong evidence that hadronic or lepto-hadronic emission processes are responsible.
Numerous authors have considered radio galaxies as candidates for neutrino events [156][157][158][159], and even FR-0-type radio galaxies were considered as well [160]. FR-0type radio galaxies are a relatively new class of radio galaxy that display a different radio morphology to the observed FR-I-and FR-II-type radio galaxies. They lack the large-scale extended radio structures of the latter two classes, instead displaying jets that are both weak and short, extending only a few tens of kiloparsecs [161]. FR-0s share similar properties to those of FR-Is except the ratio of their core to total emission is ∼30-times higher than that of FR-I-type radio galaxies [162]. Despite the abundance of FR-0s in the local universe (z < 0.05), these are difficult objects to study because of their weak and compact nature [163]. No TeV emission is currently detected from an FR-0-type radio galaxy, but there are many open questions relating to these radio galaxies, not least the fact that at least one, 4C +39.12, was detected with Fermi-LAT at GeV energies displaying a relatively hard spectrum similar to that of 3C 264.
There is also evidence that a non-negligible positional association exists between AGNs with strong parsec-scale cores and neutrino events. This implies that AGNs with strongly Doppler-boosted jets, as opposed to radio galaxies, are an important source of neutrinos [164]. In spite of this, radio galaxies are important multi-messenger targets, in particular the closest radio galaxies Cen A and M 87.
As already discussed above, single-zone SSC models can be inadequate to explain the observed γ-ray emission from radio galaxies. The rising TeV components seen in the broadband SEDs of Cen A (and IC 310 during a flaring state) have led to suggestions of either a hadronic origin or leptonic origin with multiple electron distributions [47,165]. Using Fermi-LAT data, Rulten et al. [113] analysed 26 AGNs classified as radio galaxies in the Fermi-LAT 4FGL catalogue [48]. The goal of this work was to search for Cen A-like spectral features; no other radio galaxies exhibited a Cen A-like break in their spectral energy distributions. Whether other radio galaxies possess such a break is still an open question because, given Fermi-LAT's limited sensitivity at energies above 100 GeV coupled with the poor sensitivity of the current generation of IACTs for detecting weak γ-ray sources, the possibility still exists that other radio galaxies do indeed exhibit Cen A-like spectral breaks, but we just happen to see Cen A's because of the AGN's proximity to us. Discovery of such hardening amongst a greater number of radio galaxies would strongly implicate radio galaxies as a major source of cosmic rays within the universe.

Radio Galaxies with Future Instruments
It is imperative that TeV monitoring of radio galaxies continues, both with the current generation of IACTs, as well as with CTA when that instrument comes online. Deeper observations are necessary in order to better sample the energy spectra of these objects and, importantly, detect them in different states. Ideally, multi-wavelength campaigns should be run in order to capture data simultaneously, particularly important for correlating γ-ray emission with knots downstream in the jet. For example, a quasi-simultaneous multi-wavelength campaign monitoring M 87 produced the most comprehensive multi-wavelength spectrum of the radio galaxy. Combining data from radio, optical, ultraviolet, X-ray, and gamma-ray wavebands, a near-complete broadband picture was produced during a historically low state for the AGN. In this work, the authors concluded that a structured jet is necessary to explain the observed spectrum of M 87 [166]. From modelling of the spectrum, the authors were able to exclude simultaneous γ-ray emission from inverse-Compton emission originating in the same region that produces the mm-band emission seen by the EHT. Moreover, the authors also concluded that γ-rays originating from the inner-jets (inward of HST-1) are only produced in the presence of strong particle-dominated regions [166]. However, they were not yet able to exclude direct synchrotron emission from accelerated protons (and secondaries). Such campaigns allow better modelling of radio galaxy broadband SEDs and provide the necessary information to favour a particular emission scenario for different AGN states.
Furthermore, the unequalled performance of CTA (relative to the current generation of IACTs) will allow γ-ray astronomers to investigate VHE emission from Cen A, for example with much better spatial and temporal resolution and sensitivity than ever before. Undoubtedly, this will help to better constrain the ratio of γ-ray emission between the core and extended regions of this particular radio galaxy.
New observatories such as SKA, JWST, and CTA will bring a new era of astronomy that will provide coverage of radio galaxies across the electromagnetic spectrum from radio to γ-ray frequencies, advancing our understanding of AGN physics. For example, at γ-ray energies, we know radio galaxies and blazars in general exhibit variability. The current generation of radio telescopes have limited survey speed, restricting the number of AGN with detailed variability studies to approximately a thousand of the brightest objects. Phase 1 of SKA will provide variability information at a sub-day-scale cadence for hundreds of thousands of AGN [167]. Such coverage will enable vastly improved multi-wavelength variability studies of radio galaxies between CTA and SKA, and perhaps facilitating the first insights into a radio-VHE connection, which has so far proven elusive. Furthermore, SKA is likely to play a crucial role in counterpart identification and classification of the many unidentified Fermi-LAT extragalactic γ-ray sources [168]. One of the most exciting discovery spaces in this new era regarding radio galaxies will perhaps come from the very-longbaseline interferometry (VLBI) capability of SKA complementing ngVLA advances, which is expected to provide ten-times the continuum imaging sensitivity of e-MERLIN [169] at 5 GHz in 1/10th the integration time, as well as achieving baselines of 1000 km [139]. Being able to conduct astrometry at the microarcsecond level out to tens of kiloparsecs will enable high-precision study of the proper motions of features inside radio galaxy jets [170]. Furthermore, with full linear and circular polarization capabilities, ngVLA will be able to study the magnetic field structure of radio galaxy lobes and hotspots with high precision, providing a better understanding of both the particle acceleration and jet deceleration and how this relates to γ-ray emission [139]. For example, do γ-ray flares occur after the radio emission peak and thus potentially indicating a jet "loading up" before emitting γ-rays, such as the creation of an emerging radio knot that interacts with an internal shock downstream in the jet [112]? Likewise, if no correlations are seen between γ-ray flares and knot brightening and the radio emission lags any γ-ray emission, this could be evidence that the primary γ-ray emission zone resides very close to the SMBH.

Conclusions
γ-ray observations of radio galaxies over the past decade or so have shown these objects to be fascinating, exhibiting ultra-fast variability or displaying relatively hard spectral energy distributions. These observations have also highlighted the need for instruments with improved sensitivity and resolution in order to better classify objects or to better resolve their morphological and timing characteristics, for example.
Although our understanding of radio galaxies as γ-ray emitters has advanced in the past decade or so, there still remains many outstanding questions, and a number of these were highlighted in Section 2.3. These relate to the radio galaxy environment, kinematic studies, classification considerations, γ-ray emission sites, γ-ray variability, γ-ray emission mechanisms, γ-ray emission modelling, and radio galaxies as cosmic ray sources. From the discussion in Section 4, it is clear that many of these outstanding questions will need to be addressed with multi-wavelength studies, employing a number of different instruments.
Looking forward, the near-term horizon appears very promising with many new and improved instruments arriving, such as JWST, SKA, and CTA, to name a few. They will provide astronomers with many opportunities for future observations of radio galaxies. What is more, the unprecedented coverage and resolution of these instruments will help to advance our understanding of AGN science, for example by identifying γ-ray emission sites with better precision, by constructing high-quality apparent velocity profiles over all jet scales or through deeper observations with well-sampled broadband spectra to enable better modelling of γ-ray emission.

Funding:
The author acknowledges the financial support of the UK Science and Technology Facilities Council Grant ST/V006347/1.

Data Availability Statement: Not applicable
Acknowledgments: The author would like to thank Professor Lucy Fortson (University of Minnesota) and Professor Paula Chadwick (University of Durham) for taking the time to read the article and provide valuable feedback and discussion, which has helped to improve the overall quality of the review. Finally, this research made use of the NASA/IPAC Extragalactic Database (NED), which is funded by the National Aeronautics and Space Administration and operated by the California Institute of Technology.

Conflicts of Interest:
The author declares no conflict of interest.

Abbreviations
The following abbreviations are used in this manuscript:

AGNs
Active In external Compton processes, the seed photons originate in regions external to the jet such as the accretion disk, the broad-line region, the AGN's dust torus, or molecular clouds, for example. 2 In the literature, these two features are also referred to as C1 and C3, respectively, and under this alternative naming scheme, C4 (shown in Figure 1) is referred to as C2.