Astrophysics
Comparison of AGN Properties in Luminous Radio Galaxies to their Radio-Quiet (RQ) Counterparts
My PhD thesis focuses on radiative properties of the accretion flows in radio-loud (RL) and RQ AGN and their relation to jet production efficiency. RL AGN have on average larger back hole masses (MBH) and lower Eddington ratios λE ≡ LBol/LEdd) than RQ AGN, so to avoid biases our samples are selected from the similar ranges of these parameters. Results of our studies have been published in several papers.
In Gupta et al. (2016) we studied the dependence of jet production efficiency on the outer boundary conditions of the accretion discs for RL and RQ quasars. The RL and RQ quasars were selected from the DR7 quasar catalogue, where we chose FRII quasars to represent the RL sample. We studied the variation of the covering factor (defined to be the fraction of the optical-UV radiation obscured by the circumnuclear dusty tori), with z, MBH and λE between the RL and RQ samples. We found that those tori have on average similar geometrical thickness for both RL and RQ quasars and hence jet activity is not conditioned by actual properties of accretion flows at distances larger than the dust sublimation radius. We also found a lack of statistically significant dependence of the covering factors on the Eddington ratio down to values λE < 0.03 which excludes the possibility that dusty obscurers could be dominated by the winds being powered by radiation pressure exerted on dust.
In Gupta et al. (2018) we studied the differences in X-ray properties of RL and RQ AGN. Our studies were based on the X-ray data from the Swift/BAT AGN catalogue. We compared and contrasted the hard X-ray spectral slope, X-ray loudness (defined as the ratio of the hard X-ray luminosities in the band 14-195 keV to the mid-IR luminosities at λ = 12μm), column densities, cut off energy and reflection parameter between the RL and RQ populations. We found that RL AGN are on average 2 times X-ray-louder than RQ AGN, but found their spectral slopes, reflection parameter and high energy cutoff to be very similar. These results suggest that both in RL AGN and RQ AGN samples, the X-ray emitting region can be associated with the hot, geometrically thick, central portion of the accretion flow, where production of hard X-rays is likely to be dominated by Comptonization of the optical/UV radiation of the 'cold' accretion disc by electrons heated by the magnetic reconnection. The more efficient X-ray production in RL AGN than in RQ AGN can be associated with having a larger magnetization of the innermost portions of the accretion flows and larger black-hole spins.
We built upon this work in Gupta et al. (2020) where we expanded the sample by calculating M for the objects that lacked it in literature and we collected multispectral data for the objects in our sample. We studied the X-ray loudness in the Type-1 and Type-2 samples and found that the X-ray production is isotropic in both RL and RQ AGN. The found isotropy implies that the X-ray coronas in both RQ and RL AGN cannot be too compact, otherwise their isotropic emission would be strongly affected by gravitational lensing. In addition, in our UV/mid-IR and Spectral Energy Distribution analysis we found that the integrated mid-IR luminosity reaches or even exceeds UV luminosities of Type 1 AGN, strongly implying that a significant fraction of UV radiation undergoes extinction and is reprocessed into IR radiation by the dust located within the ionization cone.
Results of all our studies seem to be consistent with the MAD (magnetically-arrested-disc, Narayan, Igumenshchev, & Abramowicz 2003) scenario, according to which powerful jets of RL AGN are powered by the rotating black-holes immersed in magnetic fields confined on them by the ram pressure of the accretion flow.
The distribution of the radio loudness, R = F1.4/FνW3 , after removing AGN with MBH < 10^8.5 M_sun is shown. Besides showing objects with radio detections and categorized as: R < 1 as RQ in blue; 1 < R < 10 as RI in yellow; and R > 10 as RL in red, we also present the radio undetected sources classified as RQU in hatched blue and RIU in hatched yellow.
Medians of the dominant SED components of the RL (red) and RQ (blue) Type 1 objects of the sample. All luminosities have been normalized w.r.t. to the Eddington luminosity. In case of FUV and NUV data the objects without detection have been presented with the flux limit on detection. We also present here SED template of the giant elliptical which has been scaled to K-band luminosity. Centre frequencies of each band are labelled.
Studying X-ray variability using XMM-Newton data
Using telescope data to study X-ray variability is one of my core interests. To wit, XMM-Newton has a high sensitivity, and the EPIC instruments possess good temporal resolution that enable the recording of rapid X-ray transients. However, the current pipeline is designed to extract lightcurves and search for variability only for sources with at least 100 EPIC counts. Since almost 2/3 of the sources have EPIC counts below this threshold, their variability is not calculated. Further, very faint sources showing very short variability may be difficult to detect. To address this, we developed EXODUS, an algorithm that identifies variability in the whole EPIC field of view irrespective of source detection and counts. It accomplishes this by binning the observations into short-time windows and comparing pixel counts per window to the median pixel counts to detect variable sources within the observation, making EXODUS computationally inexpensive.
By applying EXODUS to all observations composing the XMM catalogue, we were able to detect a plethora of new variable sources, including compact objects in nearby galaxies and other extragalactic sources that were not identified using traditional source detection methods. Our study revealed a large number of faint variable sources, including at least 14 new sources, and we showed that running EXODUS to search for short-term variability helps to detect low-mass AGN. EXODUS can be applied to future XMM-Newton observations as well as adapted to other X-ray observatories with similar detection techniques. While working on EXODUS, I developed keen familiarity with XMM-Newton, its instruments, its data processing pipeline, and how to approach its unique challenges and opportunities. I have a continued interest in obtaining and using data collected from XMM-Newton and other space telescopes.
Variability computation with the different stages of the algorithm.
PN, MOS1 and MOS2 variability of observation 0831790701. The one variable source above the variability threshold is detected at the position of the white circle. As we can observe the same source has been identified using all three EPIC cameras.
X-ray and Gamma-ray variability studies of AGN
In Nalewajko et al. (2019) we studied the long-term gamma-ray and optical spectroscopic variations in blazar 3C 454.3.
We estimated the line flux and found that despite large variations in the gamma-ray and optical continua, the line flux changes only moderately, which could be due to a change from a very high state to a deep low state. Therefore, the line flux is anti-correlated with the continuum or the increase is delayed by 600 days.
During my previous predoctoral fellowship at CfA I reduced the Chandra observational data for SWIFT J1708.6+2155, an RQ object from Gupta et al. (2020), and performed broadband fitting of the Chandra + Swift/BAT and Swift/XRT + Swift/BAT data from two epochs to determine the X-ray flux, spectral properties, and observed variability between the two observations. Based on the results, we concluded that it is unlikely that the variability comes from an obscurer, and is rather intrinsic in nature [Gupta et al. in prep].
I am also pursuing an additional project at CfA where we are studying the X-ray properties of J1548+3335, a high redshift quasar (z > 4.5). This particular source was chosen due to its steep spectra (Gamma > 2) and demonstrated variability.
Lightcurvesofblazar 3C454.3. γ-ray photon flux for photon energies above 100 MeV obtained from the FermiLAT data in 10-day bins.
AGN and Host Galaxy relations
Rusinek et al. (2020) again uses the sample from Gupta et al. (2020) in converting radio loudness to jet production efficiency, and found that a deficiency of jets with intermediate jet production efficiency implies the existence of threshold conditions
for the production of powerful jets. By studying the optical morphology, we postulated that RL AGNs hosted by spiral galaxies and RQ AGNs hosted by either spiral or elliptical galaxies, all of which are extremely rare cases, seem to favor the scenario, according to which the MAD is built up by the Cosmic Battery and can be accomplished within the AGN lifetime only in AGNs with accretion discs rotating in the opposite direction to the BH.
Additionally in Shastri et al. (2015) we studied the morphology, kinematics, and radio structure of about 120 nearby targets, and investigated the connections between their nuclear properties and the extended emission-line regions and star-formation regions to place constraints on models of the interplay between AGN outflows and their host galaxies.
During my MSc, I studied the clustering properties of galaxies (using the 2-point and projected correlation functions) and the correlation of the galaxy type with mass, colour, and UV properties of the galaxies. The objects were selected from the Galaxy and Mass Assembly Survey (GAMA) survey. We observed that the AGN in our sample have a stronger clustering property and are more massive and redder as compared to the composite and star-forming galaxies.
Distribution of radio-detected sources in the logR vs. projected size plane. Morphological groups are presented as follows : complex correspond to black crosses, triples and doubles are presented together as purple diamonds, knotty objects are marked as green stars, and the group of compact resolved sources is shown as orange circles while compact unresolved ones correspond to gray triangles. The left panel presents all radio-detected objects (257AGNs). Compact sources on the right panel have been limited to those with FIRST data (122 AGNs total)
Multiwavelength studies of Blazar
In Nalewajko & Gupta (2017) we studied blazars at multiple frequencies. Our results indicate a clean separation between FSRQs and BL Lacs that can be seen in the parameter space of MIR photon index and q = L_IC/L_syn; MIR photon index is correlated with the gamma-ray photon index, MIR luminosity, and q; and that the observed distribution of MIR luminosity versus q can be modeled as a sequence of lepto-magnetic jet powers with the preference for sub equipartition magnetic fields.
Additionally in Shastri et al. (2015) we studied the systematics of relativistic jets from blazing AGN from the third Fermi LAT Catalog (3FGL) using the multi-wavelength data in order to discern the impact of synchrotron emitting bulk relativistic jets. We observed correlations between the blazar type and their properties such as luminosity and apparent velocity of the jets.
Distribution of Fermi+WISE blazars in the parameter space of Compton dominance qFW = Lγ/Lsyn and synchrotron luminosity Lsyn, where Lγ=LSSC+LERC (FSRQs–red, BLLacs–blue). Thick solid lines show our reference model of the blazar sequence for Γj=15, PB/Pe=0.2, Γjθj=0.3, Epsilon em=0.5, and κ=1.6. Top left panel: dependence of the model on the jet magnetisation PB/Pe=0.05(dashed), 1(thin solid). Top right panel: dependence of the model on the jet Lorentz factor Γj=7.5 (dashed), 30 (thin solid). Bottom left panel: dependence of the model on the jet collimation factor Γjθj=0.15 (dashed), 0.6 (thin solid). Bottom right panel: dependence of the model on the power-law index κ in the scaling of external radiation luminosity Lext∝ Pκ eB . The black stars along each track indicate the lepto-magnetic jet power values log10 PeB= 42,43,44,45.
Refereed Publications
E. Quintin, N. Webb, I. Georgantopoulos, M. Gupta, E. Kammoun, L. Michel, A. Schwope, H. Tranin, & I. Traulsen: STONKS: Quasi-real time XMM-Newton transient detection system; Under review Astronomy & Astrophysics (2024).
N. Webb, D. Barret, O. Godet, M. Gupta, D. Lin, E. Quintin, & H. Tranin: Tidal disruption events and quasi‑periodic eruptions; Astronomische Nachrichten,344, 4 (2023).
N. Webb, .., M. Gupta, .., et. al (35 authors): XMM2ATHENA, the H2020 project to improve XMM-Newton analysis software and prepare for Athena; Astronomische Nachrichten,344, 7 (2023).
M. Gupta , M. Sikora, & K. Rusinek: Comparison of SEDs of very massive radio-loud and radio-quiet AGN; MNRAS, 492, 315 (2020) .
K. Nalewajko, A. C. Gupta, M. Liao, K. Hryniewicz, M. Gupta , & M. Gu: Long-term optical spectroscopic variations in blazar 3C 454.3; A&A, 631, A4 (2019) .
M. Gupta , M. Sikora, K. Rusinek, & GM Madejski: Comparison of hard X-rays spectra of luminous radio galaxies and their radio-quiet counterparts; MNRAS, 480, 2861 (2018).
K. Nalewajko, & M. Gupta : The sequence of Compton dominance in blazars based on data from WISE and Fermi-LAT; A&A, 606, A44 (2017) .
M. Gupta , M. Sikora, & K. Nalewajko: Covering factors of the dusty obscurers in radio-loud and radio-quiet quasars; MNRAS, 461, 234 (2016).
Conference Proceedings
M. Gupta, M. Sikora, & K. Rusinek: Comparing radio-loud Swift/BAT AGN with their radio-quiet counterparts; Bulletin of the American Astronomical Society, Vol. 53, No.1 (2021).
M. Gupta , & M. Sikora: Dependence of X-ray loudness of AGN on Eddington ratio and black hole mass; XXXVIII PAN Meeting, 7, 330 (2018) .
P. Shastri, … ., M. Gupta , et. al. (20 authors): Probing the Physics of Seyfert Galaxies Using their Emission-line Regions; Refereed Proceedings , American Institute of Physics, AIP Conference Proceedings 1697, 110002 (2015) .
P. Shastri, M. Gupta , E. Jimenez-Gomez, & G. Madejski: The Systematics of Relativistically Beamed Jets from Active Galaxies and the Blazar Divide; APS April Meeting (2015) .
P. Shastri, M. Gupta , E. Jimenez-Gomez & G. Madejski: Doppler-beamed AGN Jets and the Blazar Divide: Insights from Multi-wavelength Systematics; 29th IAU General Assembly, Honolulu, USA (2015).
P. Shastri, … ., M. Gupta , et. al. (20 authors): Probing the Interplay between AGN Outflows and their Host Galaxies: – Optical Integral Field Unit and Radio Imaging; 29th IAU General Assembly, Honolulu, USA (2015) .
P. Shastri, … ., M. Gupta , et. al. (15 authors): Probing the Physics of Seyfert Galaxies using a Spectroscopic and Multiwavelength data, Extragalactic Relativistic Jets: Cause and Effect; India (2015).
P. Shastri, E. Jimenez-Gomez, M. Gupta , & G. Madejski: Systematics of jets using multi-wavelength data; Extragalactic jets from every angle, Proceedings of the International Astronomical Union, IAU Symposium, Volume 313, pp. 164-168 (2015) .