DA 362: A Gamma-ray Emitting Compact Symmetric Object

Abstract

The $\gamma$-ray detection from an astrophysical object indicates the presence of an extreme environment where high-energy radiation is produced. With the continuous monitoring of the $\gamma$-ray sky by the Fermi Large Area Telescope (LAT), leading to deeper sensitivity, the high-energy $\gamma$-ray emission has now been detected from a diverse class of jetted active galactic nuclei (AGN). Here, we present the results of a multiwavelength study of the radio source DA 362, which was reported to be a blazar candidate of uncertain type. However, it was recently identified as a bona fide compact symmetric object (CSO) based on its sub-kpc, bi-polar radio morphology, and lack of radio variability. This makes DA 362 the only fourth $\gamma$-ray emitting object of this enigmatic class of radio-loud AGN. Using five very long baseline interferometry observations covering 1996-2018, we found the jet separation velocity to be subluminal ($v_{\rm app}\sim 0.2c$), thus supporting its CSO nature. Its Fermi-LAT observations revealed a $\gamma$-ray flaring activity, a phenomenon never detected from the other three $\gamma$-ray detected CSOs. This object is bright in the near-infrared band but extremely faint in the optical-ultraviolet filters, hinting at possible obscuration. The Swift X-Ray Telescope observation of DA 362 reveals an extremely hard X-ray spectrum, though a strong claim cannot be made due to large uncertainties. We conclude that deeper observations are needed to probe the broadband properties of this enigmatic object and to understand the origin of high-energy $\gamma$-ray emission.

1 Introduction

Compact Symmetric Objects (CSOs) are a special class of active galactic nuclei (AGN) hosting sub-kiloparsec-scale jets and exhibiting symmetric radio morphologies (e.g., Orienti, 2016). They are thought to host misaligned jets, and hence, the observed broadband emission is not expected to be beamed (cf. Wilkinson et al., 1994; Readhead et al., 1996a). These enigmatic objects are likely to be in the early stage of their evolution with kinematic ages smaller than a few thousand years (e.g., Readhead et al., 1996b). The small sizes CSOs could be due to their young age, dense galactic environment inhibiting the jet propagation, and/or recurrent/transient episodes of the nuclear jet activity (see, O’Dea & Saikia, 2021, for a review). From the analysis of their relative numbers, and redshift and linear size distributions, Kiehlmann et al. (2024a) presented strong evidence that CSOs do not evolve into larger-scale radio sources and that they should be considered as a distinct ‘short-lived’ jetted AGN population rather than ‘young’ AGN (see also, Readhead et al., 2024).

The number of $\gamma$-ray emitting CSOs remained tiny compared to more common Fanaroff-Riley Type I and II misaligned jetted AGN and blazars. Only three CSOs are reported to be detected with the Fermi Large Area Telescope (LAT): TXS 0128+554 ($z=0.036$; associated with the $\gamma$-ray source 4FGL J0131.2+5547; Lister et al., 2020), NGC 3894 ($z=0.012$; counterpart of 4FGL J1149.0+5924; Principe et al., 2020), and NGC 6328 ($z=0.015$; associated with 4FGL J1724.2$-$6501; Migliori et al., 2016). All of them are located in the nearby Universe ($z<0.05$), and their proximity could be the primary reason for their Fermi-LAT detection if the $\gamma$-ray emission is unbeamed (e.g., Principe et al., 2021). The $\gamma$ rays can be produced by the hadronic mechanisms or due to the interaction of the relativistic electrons present in the lobes with the low-energy optical-UV photons originating from the accretion disk (cf. Stawarz et al., 2008; Kino & Asano, 2011). The radio lobes typically expand with sub-relativistic velocities; therefore, the $\gamma$-ray radiation is not expected to show significant flux variations, especially on short timescales ($\sim$weeks-to-months). Indeed, none of the three $\gamma$-ray detected CSOs have exhibited significant flux variability as of now (Abdollahi et al., 2023). However, a definite conclusion cannot be drawn given the small number of known $\gamma$-ray emitting CSOs. Increasing the sample size of these peculiar sources is also crucial to investigate the radiative processes powering their jets and their interaction with the surrounding environment, compare them with other non-$\gamma$-ray detected CSOs, and understand their evolution.

Recently, Kiehlmann et al. (2024b) presented a comprehensive catalog of 79 CSOs searched from the literature and by analyzing their multi-frequency radio observations. They adopted the following four criteria: (i) projected jet length $<$1 kpc, (ii) detection of bi-polar radio emission, (iii) non-variable nature, and (iv) no superluminal motion detection in excess of $v_{\rm app}=2.5c$. In order to identify the potential $\gamma$-ray emitting CSOs, we cross-matched the CSO catalog with the fourth data release of the fourth $\gamma$-ray source catalog of Fermi-LAT detected objects (4FGL-DR4; Abdollahi et al., 2020; Ballet et al., 2023). Using a search radius of 5^′′, this exercise led to the identification of four $\gamma$-ray sources. Among them, three are already reported as $\gamma$-ray emitting CSOs in earlier works (Migliori et al., 2016; Principe et al., 2020; Lister et al., 2020). The fourth object, DA 362 (B2 1413+34), which is associated with the $\gamma$-ray source 4FGL J1416.0+3443, is classified as a blazar candidate of uncertain type in the 4FGL-DR4 catalog. However, as shown by Kiehlmann et al. (2024b), this object is a bona fide CSO, thereby making it only the fourth $\gamma$-ray detected object of this class. Moreover, Baldini et al. (2021) reported the detection of a transient $\gamma$-ray elevated activity from this object. In this article, we present the results of an investigation of its multi-wavelength properties utilizing $>$15 years of the Fermi-LAT data and other low-frequency observations and compare them with other $\gamma$-ray detected CSOs. Section 2 elaborates on the data reduction steps. The results are presented in Section 3 and discussed in Section 4. We summarize our findings in Section 5. Throughout, we adopt the convention $S_{\nu}$ $\sim$ $\nu^{\alpha}$ for spectral index $\alpha$, and use the cosmological parameters $\Omega_{m}$=0.27, $\Omega_{\bigwedge}$=0.73 and H_o=70 km s^-1 Mpc^-1.

2 Data Reduction and Analysis

2.1 Fermi-LAT

We followed the standard data reduction procedure to analyze the Fermi-LAT Pass 8 data of DA 362 covering 2008 August 4 to 2024 April 5 (MJD 54683$-$60405). We considered SOURCE class events in the energy range 0.1$-$300 GeV and lying within a region of interest of radius 15^∘ centered at the target AGN. To avoid contamination from Earth’s albedo, a zenith angle cut of $z_{\rm max}<90^{\circ}$ was also applied. All 4FGL-DR4 cataloged sources located within 20^∘ of DA 362 position were considered to model the $\gamma$-ray sky in the likelihood fitting. To take into account the diffuse background emission, we adopted the Galactic and isotropic background templates provided by the Fermi Science Support Center. The $\gamma$-ray spectral parameters of all sources were first optimized, and then the final likelihood fitting was performed by varying the parameters of all sources with detection significance $>5\sigma$ (test statistic, TS$>$25; Mattox et al., 1996).

We also generated the $\gamma$-ray spectrum (in six energy bins covering 0.1$-$300 GeV) and monthly-binned light curve of DA 362 with the same settings described above. In the time/energy bins of non-detections (TS$<$9), we computed flux upper limits at 95% confidence level.

2.2 Swift

The Neil Gehrels Swift satellite observed DA 362 on 2023 March 31 and 2023 April 5 (target id: 15938, PI: Paliya). The Swift X-Ray Telescope (XRT) data were analyzed using the online Swift-XRT data products generator (Evans et al., 2009). Since there were no previous X-ray observations of the source, we first checked its X-ray detection. An X-ray source was found at the right ascension (RA) of 14h 16m 04.03s and declination (Dec) of +34^∘ 44^′ 34.8^′′ with 90% confidence radius of 5.7 arcseconds. The angular separation between the optimized X-ray and radio positions was estimated to be 2.5 arcseconds, thus confirming that the X-ray source is spatially consistent with DA 362. Given the low-exposure individual pointings, we added both observations to generate a combined X-ray spectrum with a net exposure of 4.2 ksec in which a net 22 counts were detected from the source. The obtained spectrum was grouped using the task grppha to have at least one count per bin, and the C-statistic was employed for the spectral fitting in XSPEC (Arnaud, 1996). Uncertainties were estimated at 90% confidence level. For the plotting purpose, the fitted X-ray spectrum was re-binned to have at least 3$\sigma$ detection in each bin or grouped in sets of 3 bins.

The Swift UltraViolet and Optical Telescope (UVOT) data were reduced following the recommended guidelines. DA 362 was observed in two UV filters, namely UVW1 and UVM2. We combined the individual frames using the task uvotimsum. For the photometry, we applied the command uvotsource. The source was undetected, with UVM2 and UVW1 magnitudes fainter than 19.21 and 21.01 (Vega system), respectively.

2.3 Radio Observations

The five epochs of Very Large Baseline Interferometry (VLBI) observations of DA 362 covering 22 years (1996 to 2018) are provided on the Astrogeo website (Petrov, 2021). Four observations were in the X and S bands, and one epoch of data was taken in the C band. Additionally, the source was also detected in the Low-Frequency Array Two-metre Sky Survey, Very Large Array Sky Survey, and Rapid ASKAP Continuum Survey (Lacy et al., 2020; Shimwell et al., 2022; Duchesne et al., 2024).

2.4 Other Observations

We collected archival flux measurements from Space Science Data Center. These datasets were published in several radio-to-optical/UV catalogs (Neugebauer et al., 1984; Wright & Otrupcek, 1990; Wright et al., 1994; Epchtein et al., 1994; Gregory et al., 1996; Condon et al., 1998; Mauch et al., 2003; Skrutskie et al., 2006; Healey et al., 2007; Wright et al., 2010; Yamamura et al., 2010; Yershov, 2014; Helfand et al., 2015; Ahumada et al., 2020). Since DA 362 was undetected in the Swift-UVOT data analysis, we checked whether it was included in the Panoramic Survey Telescope and Rapid Response System catalog (PanSTARRS; Chambers et al., 2016). A faint optical source positionally consistent with DA 362 was identified.

Figure 1: The test statistic map of the $\gamma$-ray sky centered at the optimized $\gamma$-ray position of DA 362. The black circle shows the 95% uncertainty region for the $\gamma$-ray position. The radio position of DA 362 is highlighted with the ‘+’ mark.

Table 1: The $\gamma$- and X-ray spectral parameters of DA 362 and other $\gamma$-ray emitting CSOs as derived from the Fermi-LAT and Swift-XRT data analyses, respectively.

		Gamma-ray Results
4FGL Name	Counterpart	$F_{\rm 0.1-300~{}GeV}$	$\Gamma_{\rm 0.1-300~{}GeV}$	TS
		10-8 $\mathrm{ph}\,\mathrm{s}^{-1}\,\mathrm{cm}^{-2}$
J1416.0+3443	DA 362	1.26$\pm$0.07	2.71$\pm$0.02	166
J0131.2+5547	TXS 0128+554	0.56$\pm$0.02	2.05$\pm$0.02	211
J1149.0+5924	NGC 3894	0.32$\pm$0.09	2.16$\pm$0.10	133
J1724.2$-$6501	NGC 6328	0.37$\pm$0.12	2.28$\pm$0.12	48
		X-ray Results
	Exposure	$F_{\rm 0.3-10~{}keV}$	$\Gamma_{\rm 0.3-10~{}keV}$	C-stat/dof
	ksec	10-13 $\mathrm{erg}\,\mathrm{s}^{-1}\,\mathrm{cm}^{-2}$
J1416.0+3443	4.2	7.25${}^{+5.38}_{-3.17}$	0.79${}^{+0.52}_{-0.46}$	22.17/22
J0131.2+5547	19.2	4.08${}^{+0.81}_{-0.69}$	3.11${}^{+0.37}_{-0.36}$	107.94/103
J1149.0+5924	5.5	1.74${}^{+1.41}_{-1.74}$	1.57${}^{+2.20}_{-3.20}$	12.71/11
J1724.2$-$6501	35.6	7.44${}^{+1.34}_{-0.71}$	1.66${}^{+0.23}_{-0.22}$	260.59/263

3 Results

Only three CSOs were identified in the $\gamma$-ray band prior to this work, making DA 362 the only fourth $\gamma$-ray emitting object of this class. It is included in the 4FGL-DR4 catalog but was missing from previous Fermi-LAT catalogs. We carried out a dedicated data reduction covering the first $\sim$15.75 years of the Fermi-LAT operation. We also optimized the $\gamma$-ray source position of 4FGL J1416.0+3443 and estimated it to be RA = 14h 15m 55s and Dec = 34^∘ 41^′ 18^′′. The 95% uncertainty in the measured position is 0^∘.07. We show the TS map of the $\gamma$-ray region in Figure 1 where the radio and optimized $\gamma$-ray positions are also overplotted. Within the 95% uncertainty region, both positions are consistent, thus confirming the association of DA 362 with the $\gamma$-ray source 4FGL J1416.0+3443. The spectral parameters obtained from the power law fit are reported in Table 1. For comparison, we also provide the spectral parameters of the other three $\gamma$-ray detected CSOs obtained following the same methodology outlined in Section 2.1. The computed $\gamma$-ray spectral parameters are on average consistent with that published in previous works (cf. Principe et al., 2021). The minor differences, if any, could be mainly due to the different time periods for the Fermi-LAT data reduction done in our and previous works.

Furthermore, the monthly binned $\gamma$-ray light curve of DA 362 is shown in Figure 2. Though sporadically detected, the source remained mostly in quiescence during the first $\sim$12 years of the Fermi-LAT operation, thus explaining its absence in 4FGL-DR3 and earlier $\gamma$-ray catalogs. Interestingly, a flaring activity was identified during MJD 59075-59287. The brightest $\gamma$-ray flux was found to be ($1.3\pm 0.2)\times 10^{-7}$ $\mathrm{ph}\,\mathrm{s}^{-1}\,\mathrm{cm}^{-2}$ which is ten times larger than its mission averaged flux value. The detection significance of the flare peak was found to be $\sim$8$\sigma$ (TS=61).

We fitted the Swift-XRT spectrum of DA 362 with a power law model modified for the Galactic absorption fixed to the neutral hydrogen column density of $N_{\rm H}=1.32\times 10^{20}$ cm^-2 (Kalberla et al., 2005). In the energy range of 0.3$-$10 keV, the estimated photon index and energy flux are 0.79${}^{+0.52}_{-0.46}$ and 7.25${}^{+5.38}_{-3.17}\times 10^{-13}$ $\mathrm{erg}\,\mathrm{s}^{-1}\,\mathrm{cm}^{-2}$, respectively. Given the fact that CSOs are typically misaligned jetted AGN, the observation of such a flat X-ray spectrum indicates a strong absorption of soft X-ray emission usually observed in Compton thick sources (e.g., Georgantopoulos et al., 2007; Marchesi et al., 2018). Therefore, we attempted adding a redshifted absorption component (ZTBABS in XSPEC) during the fit; however, the parameters could not be constrained due to fitting failure. To get an idea about the intrinsic absorption, we froze the photon index to 1.8 and 2, i.e., to a typical value estimated for AGN, during the XSPEC fit. We obtained the fitted $N_{\rm H}$ value to be $1.05^{+1.02}_{-0.63}\times 10^{22}$ cm^-2 and $1.22^{+1.09}_{-0.69}\times 10^{22}$ cm^-2, respectively. These results provide supporting evidence about the possible obscured nature of DA 362. Deeper spectral observations will be needed to characterize the possible X-ray obscured nature of this $\gamma$-ray detected CSO.

DA 362 is extremely faint in the optical-UV band. It remained undetected in the Swift-UVOT observations. PanSTARRS catalog reports a faint optical source with $g$ and $r$-band magnitudes being 21.59$\pm$0.14 and 21.32$\pm$0.18, respectively, positionally consistent with the radio source. However, it remained undetected in the $i$ and $z$ filters. White (1992) reported its redshift to be $z=0.26$, though its origin, whether spectroscopic or photometric, is unclear. The source lies in the Sloan Digital Sky Survey coverage area; however, it lacks spectroscopic measurement, likely due to its faintness. Moreover, Stanghellini et al. (1993) provided a limit of 24 and 23 magnitudes in the $r$ and $i$ bands, respectively. In contrast, it is well detected in the mid-infrared (MIR) band covered by the Wide-field Infrared Survey Explorer (WISE).

The radio structure of DA 362 has been known to consist of a prominent jet towards the north-east, a weaker counter lobe on the western side, and a prominent central component (cf. Dallacasa et al., 2013, and references therein). We utilized the five epochs of VLBI images covering 1996-2018 available in Astrogeo (Petrov, 2021) to determine the jet velocity. In the S-band and C-band observations, we determined the separation velocity between the central component C and the jet component E, which are well-imaged in all the maps at these frequencies (Figure 4). In the higher-frequency X-band, where the source has been observed in four epochs, the central component is resolved into a double, labeled C1 and C2 in Figure 4. It is unclear which of these two components may be hosting the nucleus of the galaxy. The apparent transverse velocity $\beta_{\rm app}=\mu d_{\theta}(1+z)/c$, where $\mu$ is the observed proper motion, $d_{\theta}$ the angular size distance, $z$ is the redshift and $c$ the speed of light. The least squares fit to the data show $\beta_{\rm app}$ between C and E, and C1 and C2 to be 0.21$\pm$0.30 and 0.24$\pm$0.16, respectively, consistent with a CSO and not a blazar. The separation between C and E is 102 pc, indicating a kinematic age of $\sim$1600 yr, assuming this velocity to be a constant. Extended emission beyond component E is seen, for example, in the lower-frequency L-band image (Dallacasa et al., 1995), where they quote a total angular extent of 40 mas (161 pc) suggesting a kinematic age of $\sim$2500 yr for their outermost component. For DA 362, Kiehlmann et al. (2024b) reported the upper limit to the projected linear size to be 693 pc, considering the second lowest contours. On the other hand, we have reported the angular separation of the eastern component from the core by measuring it from the peak-to-peak positions of the core and the eastern hotspot.

4 Discussion

The detection of small, parsec-scale, bi-polar radio emission and a subluminal motion, as revealed by the VLBI datasets, provide unambiguous confirmation that DA 362 is a bona fide CSO, thereby making it only the fourth $\gamma$-ray detected object of this class. On comparing its $\gamma$-ray spectral properties with the other three $\gamma$-ray detected CSOs, we found that DA 362 is the brightest among them and exhibits a spectrum steeper than the other sources. In the $\gamma$-ray luminosity versus photon index plane, DA 362 appears to lie in a region of high luminosity and soft spectrum (Figure 5, top left panel). The $\gamma$-ray luminosity of the source was calculated assuming $z=0.26$ (White, 1992). Given that the origin of the redshift is uncertain, a $\gamma$-ray luminosity more than an order of magnitude larger than other CSOs should be treated with caution. We also compared the radio and $\gamma$-ray luminosities of CSOs with other jetted AGNs in the top middle panel of Figure 5 adopting the 8 GHz VLBI flux densities reported in the Radio Fundamental Catalog (Petrov & Kovalev, 2024, under review). In this diagram, DA 362 appears to follow the observed correlation between the radio and $\gamma$-ray luminosities that have larger values compared to other $\gamma$-ray detected CSOs.

Figure 2: The monthly binned $\gamma$-ray light curve of DA 362. The 95% flux upper limits are shown with the downward arrows. The horizontal dashed line refers to the mission averaged $\gamma$-ray flux of the source.

Figure 3: The residuals of the power law model fitting with Galactic absorption.

Figure 4: The VLBA images of DA 362 at S-band (2.29 GHz, upper left) and X-band (8.65 GHz, lower left) from observations on 2000 May 22 (from astrogeo.org uploaded by Alexandr Pushkarev). The contour levels are 3$\sigma\times-1,1,2,4,...$ where $\sigma$ is 1.1 and 0.7 mJy/beam for the S- and X-band images respectively. The linear least squares fit to the separations between components C and E (upper right) and C1 and C2 (lower right) are shown for the available data from astrogeo.org.

             

Figure 5: The plots of the $\gamma$-ray luminosity vs. spectral index, and that of radio vs. $\gamma$-ray luminosities are shown in the top left and middle panels. The top right panel shows the WISE color-color diagram. Blue and pink circles denote BL Lac objects and flat spectrum radio quasars, respectively. The yellow rectangle, green circles, and brown circles refer to steep spectrum radio quasars, radio galaxies and compact steep spectrum sources, respectively. These $\gamma$-ray sources were selected from the 4FGL-DR4 catalog (Abdollahi et al., 2020; Ballet et al., 2023). We also show $\gamma$-ray detected CSOs, as labeled. The broadband SEDs of $\gamma$-ray detected CSOs are plotted in the bottom panels.

The $\gamma$-ray emission detected from CSOs is not expected to exhibit significant variability if it originated from the radio lobes. Indeed, none of the three $\gamma$-ray emitting CSOs have displayed any flux variability so far. However, a $\gamma$-ray flaring activity of DA 362 was identified with the Fermi-LAT (Figure 2; see also Baldini et al., 2021). This observation suggests the $\gamma$-ray emission to be likely produced in the inner regions of the jet/core and not by the radio lobes. Similar results were found for another $\gamma$-ray emitting CSO TXS 0128+554 (Lister et al., 2020). However, given the low photon statistics, it is not possible to quantify the flux variability timescales to derive meaningful constraints on the location of the emission region.

We considered the Swift-XRT data of other three $\gamma$-ray emitting CSOs for a similar comparison and adopted the same data reduction and spectral fitting methodology described in Section 2.2. In particular, TXS 0128+554 was serendipitously detected when Swift-XRT was observing GRB 190203a (target id: 88751). The Swift-XRT was pointed to NGC 3894 once (target id: 89108) and it was also serendipitously observed during Swift Gravitational Wave Galaxy Survey (target ids: 3107061, 3105839). On the other hand, NGC 6328 was observed 13 times (target id: 31815, 89109). We added individual observations to generate a combined X-ray spectrum for every CSO. The obtained net counts were 124, 12, and 413, for TXS 0128+554, NGC 3894, and NGC 6328, respectively.

We fitted a model including a power law component and Galactic and intrinsic absorption. For TXS 128+554 and NGC 3894, the intrinsic neutral hydrogen column density could not be constrained, hence frozen to the values reported in recent works (Lister et al., 2020; Balasubramaniam et al., 2021). For NGC 6328, on the other hand, the data quality was good enough to determine the intrinsic column density. It was estimated to be $1.39^{+0.80}_{-0.70}\times 10^{21}$ cm^-2 which is consistent with that obtained by Bronzini et al. (2024). The obtained parameters are provided in Table 1 and we show the residuals of the fit in Figure 3. The unfolded spectra, in $\nu F_{\nu}$ versus $\nu$ format, were extracted by rebinning them to have at least 5$\sigma$ detection in each bin or grouped in sets of 5 bins. The Swift-XRT spectral fitting results for NGC 6328 were published by Matake & Fukazawa (2023) and the spectral parameters estimated in this work are fully consistent with their published values. For all other three sources, the Swift-XRT data fitting results are reported in this work for the first time. We note that NGC 3894 and NGC 6328 have also been observed with XMM-Newton, Chandra, and/or NuSTAR, and more complex spectral modeling have been performed on them (e.g., Balasubramaniam et al., 2021; Sobolewska et al., 2022; Bronzini et al., 2024). High-quality X-ray observations of DA 362, e.g., with XMM-Newton, will be needed to test physically motivated spectral models as done for other $\gamma$-ray detected CSOs.

On comparing the Swift-XRT spectral parameters of DA 362 with those estimated for other $\gamma$-ray detected CSOs, we found that it has absorption-corrected 0.3$-$10 keV flux brightness similar to other sources, though the uncertainties in the flux value are large, possibly due to low exposure. While other $\gamma$-ray emitting CSOs have X-ray photon indices similar to that typically observed from $\gamma$-ray emitting radio galaxies (cf. Matake & Fukazawa, 2023), DA 362 exhibits an extremely hard X-ray spectrum. Usually, such a flat X-ray spectrum is observed from Compton-thick AGN due to severe absorption of the soft X-ray photons. Interestingly, the WISE observations indicate it to be bright in the MIR band, while the source is extremely faint in the optical-UV band, thus making it a very red object. Combining the MIR-to-UV photometric results with the X-ray spectral parameters, a picture emerges hinting at strong dust obscuration supporting the possible obscured nature of DA 362. Deeper X-ray observations will be crucial to test this scenario.

In the WISE color-color diagram of $\gamma$-ray sources, DA 362 appears to lie in a region mainly populated by flat spectrum radio quasars (Figure 5, top right panel). Therefore, its MIR emission could originate via synchrotron emission produced by compact radio lobes expanding with mildly relativistic velocities that may not be beamed. In contrast, other $\gamma$-ray emitting CSOs are located in an area dominated by elliptical galaxies. A possible explanation could be that the MIR emission observed from these objects is primarily dominated by the thermal radiation from the host galaxy (see also Kosmaczewski et al., 2020). DA 362, on the other hand, lies in a region mainly occupied by quasars/Seyferts (Mateos et al., 2012; Stern et al., 2012), which also includes the 3 sources marked as Compton thick in Kosmaczewski et al. (2020), as well as a few other CSOs.

The broadband spectral energy distribution (SED) of DA 362 is shown in the bottom left panel of Figure 5. For comparison, we also plot the SEDs of the other three $\gamma$-ray emitting CSOs that have also been studied in previous works (Lister et al., 2020; Balasubramaniam et al., 2021; Sobolewska et al., 2022; Bronzini et al., 2024). At GHz frequencies, DA 362 has a brightness similar to TXS 0128+554 and NGC 3894, though fainter than NGC 6328. The primary difference can be seen in the MIR-to-UV energy range, where the observed emission is dominated by the host galaxy for other $\gamma$-ray emitting CSOs, which is not the case for DA 362. Deep optical-MIR photometric and spectroscopic observations are needed to characterize its host galaxy properties and accretion activity. We discussed the comparison of the X- and $\gamma$-ray properties earlier in this section, and the observed SEDs at these energies are consistent with the reported findings. We also show the multiwavelength SEDs of $\gamma$-ray emitting CSOs in $\nu L_{\nu}$ versus $\nu$ plane in the bottom right panel of Figure 5. Given the larger redshift of DA 362, it appeared more luminous compared to other $\gamma$-ray detected CSOs at all wavelengths, except at optical frequencies where it is less luminous. However, since the redshift information of DA 362 is tentative, firm conclusions about its nature cannot be drawn.

5 Summary

In this work, we have studied the multiwavelength properties of a CSO DA 362, which was recently found to be a $\gamma$-ray emitter by the Fermi-LAT, thereby making it only the fourth $\gamma$-ray detected object of this class of AGN. We summarize our key findings below:

1. We confirm the association of the $\gamma$-ray source 4FGL J1416.0+3443 with DA 362 by analyzing $\sim$15.75 years of the Fermi-LAT data. The optimized $\gamma$-ray position was consistent with the radio source within the estimated 95% $\gamma$-ray uncertainty region.

2. The monthly binned $\gamma$-ray light curve of DA 362 revealed a flaring activity during MJD 59075-59287. This is the first detection of a $\gamma$-ray flare from a CSO which was also reported by Baldini et al. (2021). This peculiar flaring activity indicates that the $\gamma$-ray emission to originate from the core/jet rather than from the radio lobes.

3. The source exhibits an extremely hard X-ray spectrum (0.3$-$10 keV photon index = 0.79${}^{+0.52}_{-0.46}$) as found with the analysis of the low-exposure Swift-XRT data. However, a strong claim cannot be made due to large uncertainties.

4. DA 362 is bright in MIR but extremely faint in the optical band, thus suggesting possible dust obscuration. By also considering the observed X-ray spectral shape, these results indicate a possible X-ray obscured nature of the source.

5. We used the calibrated VLBI images from the Astrogeo website and estimated the jet separation velocity to be $v_{\rm app}\sim 0.2c$. This detection of a subluminal motion further supports the CSO nature of DA 362.

6. The available observations have provided tantalizing clues about the enigmatic behavior of this $\gamma$-ray emitting CSO. Deeper observations with sensitive observing facilities will be needed to explore the broadband physical properties of DA 362 and probe the origin of $\gamma$-ray emission.