Search

We present a comprehensive study of Ultraluminous Infrared Galaxies (ULIRGs), leveraging data from the IRAS Faint Source Catalogue (FSC) and the spectroscopic catalog in the Sloan Digital Sky Survey (SDSS) DR16. Our meticulous cross-matching technique significantly enhances the reliability of ULIRG identification, resulting in the identification of 283 reliable ULIRGs, including 102 new detections, while discarding 120 previously reported false sources. Covering a redshift range of $z = 0.018 - 0.996$, with a median redshift of $\bar{z} = 0.259$, our uniform sample reveals apparent interaction features in approximately 40\% of ULIRGs, increasing to 92\% for those with $z < 0.1$. Through optical spectra analysis, it is indicated that over 58\% of ULIRGs host an AGN, which is twice as high as the detections based solely on infrared colors. Moreover, a pronounced excess of radio emissions associated with AGN activity results in a steeper radio-far-infrared correlation. Notably, Type I ULIRGs exhibit properties similar to those of narrow-line Seyfert 1 galaxies (NLS1s), with an elevated incidence rate of \ion{Mg}{2} BALs (16.7\%), surpassing that of typical optically selected quasars by over tenfold, consistent with current evolutionary models. We anticipate that forthcoming telescopes such as the China Space Station Telescope (CSST) and Leighton Chajnantor Telescope (LCT) will provide deeper insights into ULIRG morphology, dust distribution, molecular gas, and AGN activity., Comment: Accepted for publication in ApJS, 33 pages, 13 figures, and 1 table

The particle acceleration of blazar jets is crucial to high-energy astrophysics, yet the acceleration mechanism division in blazar subclasses and the underlying nature of these mechanisms remain elusive. In this work, we utilized the synchrotron spectral information (synchrotron peak frequency, $\log \nu_{\rm sy}$, and corresponding curvature, $b_{\rm sy}$) of 2705 blazars from the literature and studied the subject of particle acceleration in blazar jets by analysing the correlation between $\log \nu_{\rm sy}$ and $1/b_{\rm sy}$. Our results suggested that the entire sample follows an energy-dependent probability acceleration (EDPA). Specifically, the low inverse Compton peak sources (LCPs) follow the mechanism that fluctuations of fractional gain acceleration (FFGA), while the high inverse Compton peak sources (HCPs) follow an acceleration mechanism of EDPA. Our results indicated that the separation between LCPs and HCPs results from the electron peak Lorentz factor ($\gamma_{\rm p}$), and the differentiation should originate from different acceleration mechanisms. Moreover, our study revealed a transition in the acceleration mechanism from FFGA to EDPA around $\log \nu_{\rm sy} \sim 15$ through a detailed analysis of binned-$\log \nu_{\rm sy}$. The mechanism of FFGA dominates the particle acceleration in LCP jets because of stronger jets and the EDPA dominates the particle energy gain in the HCPs due to a more efficient acceleration process., Comment: Accepted to ApJ

In the forthcoming era of big astronomical data, it is a burden to find out target sources from ground-based and space-based telescopes. Although Machine Learning (ML) methods have been extensively utilized to address this issue, the incorporation of in-depth data analysis can significantly enhance the efficiency of identifying target sources when dealing with massive volumes of astronomical data. In this work, we focused on the task of finding AGN candidates and identifying BL Lac/FSRQ candidates from the 4FGL DR3 uncertain sources. We studied the correlations among the attributes of the 4FGL DR3 catalogue and proposed a novel method, named FDIDWT, to transform the original data. The transformed dataset is characterized as low-dimensional and feature-highlighted, with the estimation of correlation features by Fractal Dimension (FD) theory and the multi-resolution analysis by Inverse Discrete Wavelet Transform (IDWT). Combining the FDIDWT method with an improved lightweight MatchboxConv1D model, we accomplished two missions: (1) to distinguish the Active Galactic Nuclei (AGNs) from others (Non-AGNs) in the 4FGL DR3 uncertain sources with an accuracy of 96.65%, namely, Mission A; (2) to classify blazar candidates of uncertain type (BCUs) into BL Lacertae objects (BL Lacs) or Flat Spectrum Radio Quasars (FSRQs) with an accuracy of 92.03%, namely, Mission B. There are 1354 AGN candidates in Mission A, 482 BL Lacs candidates and 128 FSRQ candidates in Mission B were found. The results show a high consistency of greater than 98% with the results in previous works. In addition, our method has the advantage of finding less variable and relatively faint sources than ordinary methods.

Variability is a prominent observational feature of blazars. The high-energy radiation mechanism of jets has always been important but still unclear. In this work, we performed a detailed analysis using Fermi-LAT data across 15 years and obtained GeV light curve information for 78 TeV blazars detected by Fermi. We provided annual GeV fluxes and corresponding spectral indices for the 78 TeV blazars and thorough monthly GeV fluxes for a subsample of 41 bright blazars. Our results suggest a strong correlation between the ${\gamma}$-ray photon index and $\log L_{\rm \gamma}$ for the flat spectrum radio quasars (FSRQs) and high-energy-peaked BL Lacs (HBLs). 14 sources in our sample show significant GeV outbursts/flares above the relatively stable, low-flux light curve, with 6 of them showing a clear sharp peak profile in their 5-day binned light curves. We quantified the variability utilizing the fractional variability parameter $F_{\mathrm{var}}$, and found that the flux of the FSRQs showed significantly stronger variability than that of the BL Lacs. The 41 bright blazars in this work are best fit by a log-normal flux distribution. We checked the spectral behavior and found 11 out of the 14 sources show a 'bluer-when-brighter (BWB)' trend, suggesting this spectral behavior for these TeV blazars at the GeV band arises from the mechanism that the synchrotron-self Compton (SSC) process dominates the GeV emission. Our research offers a systematic analysis of the GeV variability properties of TeV blazars and serves as a helpful resource for further associated blazar studies., Comment: 22 pages, 9 figures, accepted for publication in ApJS

The studies and constraints on the emission region are crucial to the blazar radiation mechanism. Yet the previous works mainly focus on individual sources. In this work, we make use of the largest and the latest spectral energy distribution (SED) fitting results in the literature to statistically study the blazar emission region property in the framework of leptonic one-zone. Our results reveal (1) FSRQs show lower electron energy ($\gamma_{\rm p} \lesssim 1.6 \times 10^{3}$) than BL Lacs and tend to have a stronger magnetic field ($B$) and smaller electron-to-magnetic energy ratio ($U_{\rm e}/U_{\rm B}$) than BL Lacs; (2) we find the electro-magnetic equipartition would rather happen in the jets of BL Lacs than happen in the jets of FSRQs; (3) there are 682 blazars with a magnetic field weaker critical value of generating the Kelvin-Helmholtz instability, thus one-third of the blazars in our sample are able to produce this instability; (4) the distance ($d_{\rm em}$) between the emission region and the central black hole (BH) is in the scale of $\sim$0.1 pc, the location of the emission region may be evenly distributed inside and outside the broad line region (BLR)., Comment: 13 pages, 5 figures, 2 tables. Accepted to ApJS

The `blazar sequence' has been proposed for more than 20 years, yet its nature is still unclear. In this work, for the first time, we expand this topic to the TeV band by using a sample of 58 TeV blazars including 48 blazars in the quiescent state and 21 blazars in the flaring state. We investigate the correlation between the TeV luminosity, which has been compensated for attenuation from extragalactic background light, and the synchrotron peak frequency. We note that there is no correlation between TeV luminosity and peak frequency in the quiescent state and a strong anti-correlation in the flaring state for the observed value. However, there is a strong positive correlation in both the quiescent state and the flaring state for the intrinsic value. This indicates that the blazar sequence is shown in the flaring state rather than in the quiescent state for the observed value and the blazar sequence is not present in both two states after removing the beaming effect. In addition, to confirm whether the beaming effect results in the blazar sequence, we compare the \textit{Fermi} $\gamma$-ray luminosity between the quiescent state and the flaring state. We find the \textit{Fermi} $\gamma$-ray luminosity in the flaring state is greater than that in the quiescent state and the Doppler factor in the flaring state is greater. We suggest the blazar sequence in the flaring state may be due to the stronger beaming effect., Comment: 19 pages, 4 figures and 4 tables. Accepted to APJ

3C 279 showed enhanced flux variations in Fermi-LAT {\gamma}-ray observations from January to June 2018. We present a detailed Fermi-LAT analysis to investigate the variability and spectral behaviors of 3C 279 during the {\gamma}-ray flares in 2018. In this work, we analyzed the {\gamma}-ray spectra and found that the spectra in either the flaring or quiescent states do not show any clear breaks (or cutoffs). This indicates that the dissipation region is outside the broad-line region, and the energy dissipation may be due to the inverse Compton process of scattering the dust torus infrared photons, this result is also consistent with that in Tolamatti et al. An external inverse Compton scattering of dusty torus (DT) photons is employed to calculate the broadband spectral energy distribution (SED). This model was further supported by the fact that we found flare decay timescale was consistent with the cooling time of relativistic electrons through DT photons. During the SED modeling, a relatively harder spectrum for the electron energy distribution (EED) is found and suggests these electrons may not be accelerated by the shock that happened in the dissipation region. Besides, the magnetic reconnection is also ruled out due to a low magnetization ratio. Thus, we suggest an injection of higher-energy electrons from outside the blob and raising the flare., Comment: 12 pages, 6 figures, published in the Publications of the Astronomical Society of the Pacific

It is known that the blazar jet emissions are dominated by non-thermal radiation while the accretion disk jets are normally dominated by thermal emission. In this work, our aim is to study the connection between the two types of emission by investigating the correlation between the blazar emission line intensity property, which embodies the nature of accretion disk, and the $\gamma$-ray flux property, which is the representative of jet emission. We compiled a sample of 656 blazars with available emission line equivalent widths ($EW$), the GeV $\gamma$-ray flux, and the SED information from the literature. In this work, we found 55 previous BCUs are now identified as FSRQs, and found 52 Changing-look blazars based on their $EW$ and 45 of them are newly confirmed. These Changing-look blazars have a larger accretion ratio (${\dot M}/{\dot M}_{\rm Edd}$) than BL Lac objects. In addition, we suggest that the lower synchrotron peak blazars (LSPs) could be the source of Changing-look blazars because 90.7\% of the Changing-look blazars in this work are confirmed as LSPs. An anti-correlation between $EW$ and continuum intensity, the so-called global Baldwin effect (BEff) has been confirmed. We suggest the steeper global BEff observed for blazar than for radio-quiet active galactic nuclei (RQ-AGNs) is caused by the inverse Compton scattering of broad-emission-line photons. This interpretation is further supported by the positive correlation between the emission line $EW$ and intrinsic inverse Compton luminosity., Comment: Accepted to ApJ

To investigate the properties of dust in the nuclear region of quasars, we explored the extinction curves of the iron low-ionization broad absorption line (FeLoBAL) quasar SDSS J163004.29+311957.6 and its two analogues. The parameterized extinction curves indicated the Milky Way-like 2175 \AA~ bump features in underlying extinction, which are similar to those seen in the Local Group and a subset of high-redshift star-forming galaxies. Compared to the bump features in the Large Magellanic Clouds (LMC), the detections in this work are much closer to those in the Milky Way (MW). These bump features, as well as those in the high- and low-ionization broad absorption line (BAL) quasars of Zhang et al., are probably the counterpart of the 2175 \AA~ bump features in the quasar environment. This type of dust grain is generally small, easily disrupted by high-energy photons and has difficulty surviving in the radiation field of the active galactic nucleus (AGN). However, due to the presence of absorption-line outflows, the 2175 \AA~ bump feature in quasars, which should be rare, is seen many times in BAL quasars. The shielding effect of outflow clouds allows the MW-like dust grains to be assembled or extends the survival period in the quasar nuclear region. The process, and physical and chemical conditions deserve further observational study and investigation., Comment: 13 pages, 10 figures. Accepted by Astronomy & Astrophysics

The question of radio dichotomy in the active galactic nuclei (AGNs) is still in debate even it has been proposed for more than forty years. In order to solve the old riddle, we collect a sample of AGNs with optical $B$ band and radio 6cm wavelength data to analyze the radio loudness ${\rm log}R$. Our results indicate a separation of ${\rm log}R = \langle 1.37 \pm 0.02 \rangle$ between radio-loud (RL) AGNs and radio-quiet (RQ) AGNs, suggest the existence of an RL/RQ dichotomy. For the first time, we suggest combining radio luminosity and radio loudness as a double-criterion to divide AGNs into RLs and RQs to avoid misclassification problems that may happen in the single-criterion scenario, we propose the double-criterion dividing line ${\rm log}L_{\rm 6cm} = -2.7{\rm log}R +44.3$ by using a machine learning method. In addition, the key point of the RL/RQ dichotomy is the origin of radio emission for the two classes, we suggest the radio emission from RLs and RQs share the same origin, e.g. jets and mini-jets (aborted-jet or outflow), through a correlation study between radio 6cm luminosity and optical $B$ band luminosity., Comment: 17 pages, 5 figures. Published in PASJ

Jet origination is one of the most important questions of AGN, yet it stays obscure. In this work, we made use of information of emission lines, spectral energy distributions (SEDs), \textit{Fermi}-LAT $\gamma$-ray emission, construct a blazar sample that contains 667 sources. We notice that jet power originations are different for BL Lacs and for FSRQs. The correlation between jet power $P_{\rm jet}$ and the normalized disk luminosity $L_{\rm Disk}/L_{\rm Edd}$ shows a slope of -1.77 for BL Lacs and a slope of 1.16 for FSRQs. The results seem to suggest that BL Lac jets are powered by extracting blackhole rotation energy, while FSRQ jets are mostly powered by accretion disks. Meanwhile, we find the accretion ratio $\dot{M} / \dot{M}_{\rm Edd}$ increase with the normalized $\gamma$-ray luminosity. Base on this, we propose a dividing line, ${\rm log} (L_{\rm BLR}/L_{\rm Edd}) = 0.25 \ {\rm log} (L_{\rm \gamma}/L_{\rm Edd}) - 2.23$, to separate FSRQs and BL Lacs in the diagram of $L_{\rm BLR}/L_{\rm Edd}$ against $L_{\rm \gamma}/L_{\rm Edd}$ through using the machine learning method, the method gives an accuracy of 84.5\%. In addition, we propose an empirical formula, $M_{\rm BH}/M_{\rm \odot} \simeq L_{\rm \gamma}^{0.65}/21.46$, to estimate blackhole mass based on a strong correlation between $\gamma$-ray luminosity and blackhole mass. Strong $\gamma$-ray emission is typical in blazars, and the emission is always boosted by a Doppler beaming effect. In this work, we generate a new method to estimate a lower-limit of Doppler factor $\delta$ and give $\delta_{\rm BL Lac} = 7.94$ and $\delta_{\rm FSRQ} = 11.55$., Comment: Accepted to ApJ

In this work, Machine Learning (ML) methods are used to efficiently identify the unassociated sources and the Blazar Candidate of Uncertain types (BCUs) in the Fermi-LAT Third Source Catalog (3FGL). The aims are twofold: 1) to distinguish the Active Galactic Nuclei (AGNs) from others (non-AGNs) in the unassociated sources; 2) to identify BCUs into BL Lacertae objects (BL Lacs) or Flat Spectrum Radio Quasars (FSRQs). Two dimensional reduction methods are presented to decrease computational complexity, where Random Forest (RF), Multilayer Perceptron (MLP) and Generative Adversarial Nets (GAN) are trained as individual models. In order to achieve better performance, the ensemble technique is further explored. It is also demonstrated that grid search method is of help to choose the hyper-parameters of models and decide the final predictor, by which we have identified 748 AGNs out of 1010 unassociated sources, with an accuracy of 97.04%. Within the 573 BCUs, 326 have been identified as BL Lacs and 247 as FSRQs, with an accuracy of 92.13%., Comment: Accepted for publication in Astronomy and Computing

In our previous work in Xiao et al. (2019), we suggested that 6 superluminal sources could be gamma-ray candidates, and in fact 5 of them have been confirmed in the fourth Fermi-LAT source catalogue (4FGL). In this work, based on the 4FGL, we report a sample of 229 Fermi detected superluminal sources (FDSs) including 40 new FDSs and 62 non-Fermi detected superluminal sources (non-FDSs). Thus, we believe that all superluminal sources should have $\gamma$-ray emissions, and superluminal motion could also be a clue to detect $\gamma$-ray emission from active galactic nuclei (AGN). We present a new approach of Doppler factor estimate through the study of the $\gamma$-ray luminosity ($L_{\rm \gamma}$) and of the viewing angle ($\phi$).

Active galactic nuclei (AGNs) have been attracting research attention due to their special observable properties. Specifically, a majority of AGNs are detected by Fermi-LAT missions, but not by Fermi-LAT, which raises the question of whether any differences exist between the two. To answer this issue, we compile a sample of 291 superluminal AGNs (189 FDSs and 102 non-FDSs) from available multi-wavelength radio, optical, and X-ray (or even $\gamma$-ray) data and Doppler factors and proper motion ($\mu$) (or apparent velocity ($\beta_{\rm{app}}$)); calculated the apparent velocity from their proper motion, Lorentz factor ($\Gamma$), viewing angle ($\phi$) and co-moving viewing angle ($\phi_{co}$) for the sources with available Doppler factor ($\delta$); and performed some statistical analyses for both types. Our study indicated that1. In terms of average values, FDSs have higher proper motions ($\mu$), apparent velocities ($\beta_{\rm app}$), Doppler factor ($\delta$), Lorentz factor ($\Gamma$), and smaller viewing angle ($\phi$). Nevertheless, there is no clear difference in co-moving viewing angles ($\phi_{\rm co}$)., Comment: 19 pages, 9 figures, accepted to SCPMA