Showing 1–50 of 133 results for author: Stalin, C

Spectro-polarimetric view of the gamma-ray emitting NLS1 1H0323+342

Authors: Jincen Jose, Suvendu Rakshit, Swayamtrupta Panda, Jong-Hak Woo, C. S. Stalin, Neha Sharma, Shivangi Pandey

Abstract: The gamma-ray emitting narrow-line Seyfert 1 galaxies are a unique class of objects that launch powerful jets from relatively lower-mass black hole systems compared to the Blazars. However, the black hole masses estimated from the total flux spectrum suffer from the projection effect, making the mass measurement highly uncertain. The polarized spectrum provides a unique view of the central engine… ▽ More The gamma-ray emitting narrow-line Seyfert 1 galaxies are a unique class of objects that launch powerful jets from relatively lower-mass black hole systems compared to the Blazars. However, the black hole masses estimated from the total flux spectrum suffer from the projection effect, making the mass measurement highly uncertain. The polarized spectrum provides a unique view of the central engine through scattered light. We performed spectro-polarimetric observations of the gamma-ray emitting narrow-line Seyfert 1 galaxy 1H0323+342 using SPOL/MMT. The degree of polarization and polarization angle is 0.122 $\pm$ 0.040 % and 142 $\pm$ 9 degrees, while the H$α$ line is polarized at 0.265 $\pm$ 0.280 %. We decomposed the total flux spectrum and estimated broad H$α$ FWHM of 1015 km/s. The polarized flux spectrum shows a broadening similar to the total flux spectrum, with a broadening ratio of 1.22. The Monte Carlo radiative transfer code `STOKES' applied to the data provides the best fit for a small viewing angle of 9-24 degrees and a small optical depth ratio between the polar and the equatorial scatters. A thick BLR with significant scale height can explain a similar broadening of the polarized spectrum compared to the total flux spectrum with a small viewing angle. △ Less

Submitted 11 July, 2024; v1 submitted 9 July, 2024; originally announced July 2024.

Comments: Accepted for publication in MNRAS (Jul. 08th, 2024). 12 pages, 6 figures

UVIT Survey of AGN host Galaxies -- I: Star Formation Scenarios

Authors: Payel Nandi, C. S. Stalin, Poulomi Dam, D. J. Saikia

Abstract: Circumnuclear star formation (SF) is generally seen in galaxies hosting active galactic nuclei (AGN); however, the connection between the AGN activity and SF in them is less well understood. To explore this connection on scales of a few tens of parsec to a few tens of kiloparsec and larger, we carried out an investigation of SF in seven Seyfert type AGN and one LINER galaxy, using observations wit… ▽ More Circumnuclear star formation (SF) is generally seen in galaxies hosting active galactic nuclei (AGN); however, the connection between the AGN activity and SF in them is less well understood. To explore this connection on scales of a few tens of parsec to a few tens of kiloparsec and larger, we carried out an investigation of SF in seven Seyfert type AGN and one LINER galaxy, using observations with the Ultra-Violet Imaging Telescope (UVIT) on board {\it AstroSat} in the near ultra-Violet (NUV; 2000-3000 Å) and far ultra-Violet (FUV; 2000-3000 Å) bands. A total of 1742 star-forming regions were identified, having size scales of 0.010 to 63.642 kpc$^2$. Considering all the galaxies, we found a positive correlation between their total surface density of SF ($Σ_{SFR}$) and extinction. For five galaxies, namely NGC 1365, NGC 4051, NGC 4321, NGC 5033 and NGC 6814, we found a gradual decrease of both extinction and $Σ_{SFR}$ from the centre to the outer regions. Four sources are found to lie in the main sequence (MS) of star-forming galaxies, and the other four are away from MS. We found the ratio of the star formation rate (SFR) in the nuclear region to the total SFR to be positively correlated with the Eddington ratio. This points to the influence of AGN on enhancing the SF characteristics of the hosts. However, the impact is dominant only in the central nuclear region with no significant effect on the larger scales probed in this work. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: Accepted for the publication of The Astrophysical Journal; Comments are welcome

A comparative analysis of the active galactic nucleus and star formation characteristics of broad and narrow line Seyfert 1 galaxies

Authors: Kshama S. Kurian, C. S. Stalin, S. Rakshit, G. Mountrichas, D. Wylezalek, R. Sagar, M. Kissler-Patig

Abstract: We report here our comparative analysis of active galactic nucleus (AGN) and star formation (SF) characteristics of a sample of narrow-line Seyfert 1 (NLS1) and broad-line Seyfert 1 (BLS1) galaxies. Our sample consisted of 373 BLS1 and 240 NLS1 galaxies spanning the redshift 0.02 to 0.8. The broad-band spectral energy distribution, constructed using data from the ultra-violet to the far-infrared,… ▽ More We report here our comparative analysis of active galactic nucleus (AGN) and star formation (SF) characteristics of a sample of narrow-line Seyfert 1 (NLS1) and broad-line Seyfert 1 (BLS1) galaxies. Our sample consisted of 373 BLS1 and 240 NLS1 galaxies spanning the redshift 0.02 to 0.8. The broad-band spectral energy distribution, constructed using data from the ultra-violet to the far-infrared, was modelled using CIGALE to derive the basic properties of our sample. We searched for differences in stellar mass, star formation rate (SFR), and AGN luminosity in the two populations. We also estimated new radiation-pressure corrected black hole masses for our sample of BLS1 and NLS1 galaxies. While the virial black hole mass (MBH) of BLS1 galaxies is similar to their radiation-pressure corrected MBH values, the virial MBH values of NLS1 galaxies are underestimated. We found that NLS1 galaxies have a lower MBH of log (MBH) = 7.45 +/- 0.27 and a higher log (Eddington ratio) of -0.72 +/- 0.22 than BLS1 galaxies, which have log (MBH) and log (Eddington ratio) of 8.04 +/- 0.26 and -1.08 +/- 0.24, respectively. The distributions of stellar mass, SFR, and specific star formation (sSFR = SFR/stellar mass) for the two populations are indistinguishable. This analysis is based on an independent approach and contradicts reports in the literature that NLS1 galaxies have a higher SF than BLS1 galaxies. While we found that AGN luminosity increases with stellar mass, luminosity of SF flattens at high stellar mass for both BLS1 and NLS1 galaxies. The reason may be that SF is suppressed by AGN feedback at stellar mass higher than 10^11 solar mass or that the AGN fuelling mechanism is decoupled from SF. Separating the sample into radio-detected and radio-undetected subsamples, we found no difference in their SF properties suggesting that the effect of AGN jets on SF is negligible. △ Less

Submitted 17 June, 2024; originally announced June 2024.

Comments: 15 pages, accepted for publication in A&A

Revisiting the dust torus size-luminosity relation based on a uniform reverberation mapping analysis

Authors: Amit Kumar Mandal, Jong-Hak Woo, Shu Wang, Suvendu Rakshit, Hojin Cho, Donghoon Son, C. S. Stalin

Abstract: We investigate the torus size -- luminosity relation of Type 1 AGNs based on the reverberation-mapping analysis using the light curves of the optical continuum and the IR continuum obtained with the W1 and W2-bands of the Wide-field Infrared Survey Explorer (WISE) survey. The final sample consists of 446 and 416 AGNs, respectively, for W1 and W2-band light curves, covering a large dynamic range of… ▽ More We investigate the torus size -- luminosity relation of Type 1 AGNs based on the reverberation-mapping analysis using the light curves of the optical continuum and the IR continuum obtained with the W1 and W2-bands of the Wide-field Infrared Survey Explorer (WISE) survey. The final sample consists of 446 and 416 AGNs, respectively, for W1 and W2-band light curves, covering a large dynamic range of bolometric luminosity from $10^{43.4}$ to $10^{47.6}$ $erg \, s^{-1}$, which show reliable lag measurements based on our quality assessment analysis. After correcting for the accretion disk contamination in the observed IR flux, we constrain the torus size ($R_{dust}$) and AGN bolometric luminosity ($L_{bol}$) relationship with the best-fit slope of 0.39 (0.33) for the W1- (W2-) band, which is shallower than expected from the dust radiation equilibrium model. By combining the previous K-band lag measurements, we find that the measured torus size depends on the observed wavelength of the dust radiation, as $R_{dust,K}:R_{dust,W1}:R_{dust,W2}$ = 1.0:1.5:1.8 ($R_{dust} \, \propto \, λ^{0.80}$) at $L_{bol}$ = $10^{46} \, erg \, s^{-1}$, confirming a stratified structure of the torus, where wavelength-dependent emissions originate from distinct regions of the torus. By investigating the deviation from the best-fit torus size -- luminosity relation, we find a moderate correlation between the offset from the $R_{dust}$--$L_{bol}$ relation and Eddington ratio. This suggests a possible influence of the Eddington ratio on the observed flattening of the $R_{dust}$--$L_{bol}$ relationship. △ Less

Submitted 19 April, 2024; v1 submitted 4 March, 2024; originally announced March 2024.

Comments: 35 pages, 21 figures, 6 tables. Corrected for typos. Accepted for publication in ApJ

Multi-wavelength spectral modelling of the candidate neutrino blazar PKS 0735+178

Authors: Athira M Bharathan, C. S. Stalin, S. Sahayanathan, Subir Bhattacharyya, Blesson Mathew

Abstract: The BL Lac object PKS 0735+178 was in its historic $γ$-ray brightness state during December 2021. This period also coincides with the detection of a neutrino event IC211208A, which was localized close to the vicinity of PKS 0735+178. We carried out detailed $γ$-ray timing and spectral analysis of the source in three epochs (a) quiescent state ($E_{1}$), (b) moderate activity state ($E_{2}$) and (c… ▽ More The BL Lac object PKS 0735+178 was in its historic $γ$-ray brightness state during December 2021. This period also coincides with the detection of a neutrino event IC211208A, which was localized close to the vicinity of PKS 0735+178. We carried out detailed $γ$-ray timing and spectral analysis of the source in three epochs (a) quiescent state ($E_{1}$), (b) moderate activity state ($E_{2}$) and (c) high activity state ($E_{3}$) coincident with the epoch of neutrino detection. During the epoch of neutrino detection ($E_{3}$), we found the largest variability amplitude of 95%. The $γ$-ray spectra corresponding to these three epochs are well fit by the power law model and the source is found to show spectral variations with a softer when brighter trend. In the epoch $E_{3}$, we found the shortest flux doubling/halving time of 5.75 hrs. Even though the spectral energy distribution in the moderate activity state and in the high activity state could be modeled by the one-zone leptonic emission model, the spectral energy distribution in the quiescent state required an additional component of radiation over and above the leptonic component. Here we show that a photo-meson process was needed to explain the excess $γ$-ray emission in the hundreds of GeV which could not be accounted for by the synchrotron self-Compton process. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: Accepted by MNRAS

Investigation of the correlation between optical and $γ$-ray flux variation in the blazar Ton 599

Authors: Bhoomika Rajput, Amit Kumar Mandal, Ashwani Pandey, C. S. Stalin, Walter Max-Moerbeck, Blesson Mathew

Abstract: The correlation between optical and $γ$-ray flux variations in blazars reveals a complex behaviour. In this study, we present our analysis of the connection between changes in optical and $γ$-ray emissions in the blazar Ton 599 over a span of approximately 15 years, from August 2008 to March 2023. Ton 599 reached its highest flux state across the entire electromagnetic spectrum during the second w… ▽ More The correlation between optical and $γ$-ray flux variations in blazars reveals a complex behaviour. In this study, we present our analysis of the connection between changes in optical and $γ$-ray emissions in the blazar Ton 599 over a span of approximately 15 years, from August 2008 to March 2023. Ton 599 reached its highest flux state across the entire electromagnetic spectrum during the second week of January 2023. To investigate the connection between changes in optical and $γ$-ray flux, we have designated five specific time periods, labeled as epochs A, B, C, D, and E. During periods B, C, D, and E, the source exhibited optical flares, while it was in its quiescent state during period A. The $γ$-ray counterparts to these optical flares are present during periods B, C, and E, however during period D, the $γ$-ray counterpart is either weak or absent. We conducted a broadband spectral energy distribution (SED) fitting by employing a one-zone leptonic emission model for these epochs. The SED analysis unveiled that the optical-UV emission primarily emanated from the accretion disk in quiescent period A, whereas synchrotron radiation from the jet dominated during periods B, C, D, and E. Diverse correlated patterns in the variations of optical and $γ$-ray emissions, like correlated optical and $γ$-ray flares, could be accounted for by changes in factors such as the magnetic field, bulk Lorentz factor, and electron density. On the other hand, an orphan optical flare could result from increased magnetic field and bulk Lorentz factor. △ Less

Submitted 22 December, 2023; originally announced December 2023.

Comments: 16 pages, 8 figures, Accepted for publishing in the MNRAS

Narrow-line Seyfert 1 galaxies in Sloan Digital Sky Survey: a new optical spectroscopic catalogue

Authors: Vaidehi S. Paliya, C. S. Stalin, Alberto Dominguez, D. J. Saikia

Abstract: Narrow-line Seyfert 1 (NLSy1) galaxies are an enigmatic class of active galactic nuclei (AGN) that exhibit peculiar multiwavelength properties across the electromagnetic spectrum. For example, these sources have allowed us to explore the innermost regions of the central engine of AGN using X-ray observations and have also provided clues about the origin of relativistic jets considering radio and g… ▽ More Narrow-line Seyfert 1 (NLSy1) galaxies are an enigmatic class of active galactic nuclei (AGN) that exhibit peculiar multiwavelength properties across the electromagnetic spectrum. For example, these sources have allowed us to explore the innermost regions of the central engine of AGN using X-ray observations and have also provided clues about the origin of relativistic jets considering radio and gamma-ray bands. Keeping in mind the ongoing and upcoming wide-field, multi-frequency sky surveys, we present a new catalogue of NLSy1 galaxies. This was done by carrying out a detailed decomposition of >2 million optical spectra of quasars and galaxies from the Sloan Digital Sky Survey Data Release 17 (SDSS-DR17) using the publicly available software "Bayesian AGN Decomposition Analysis for SDSS Spectra". The catalogue contains 22656 NLSy1 galaxies which is more than twice the size of the previously identified NLSy1s based on SDSS-DR12. As a corollary, we also release a new catalogue of 52273 broad-line Seyfert 1 (BLSy1) galaxies. The estimated optical spectral parameters and derived quantities confirm the previously known finding of NLSy1 galaxies being AGN powered by highly accreting, low-mass black holes. We conclude that this enlarged sample of NLSy1 and BLSy1 galaxies will enable us to explore the low-luminosity end of the AGN population by effectively utilizing the sensitive, high-quality observations delivered by ongoing/upcoming wide-field sky surveys. The catalogue has been made public at https://www.ucm.es/blazars/seyfert △ Less

Submitted 23 November, 2023; originally announced November 2023.

Comments: MNRAS, in press

Detection of X-ray Polarization in the High Synchrotron Peaked Blazar 1ES 1959+650

Authors: Athira M. Bharathan, C. S. Stalin, Rwitika Chatterjee, S. Sahayanathan, Indrani Pal, Blesson Mathew, Vivek K. Agrawal

Abstract: We report the first time measurement of X-ray polarization in the 2$-$8 keV band from the high synchrotron peaked (HSP) blazar 1ES 1959+650. The low energy hump in the broadband spectral energy distribution (SED) of blazars in the log$ν$ versus $νF_ν$ plane is believed to be due to the synchrotron emission process from relativistic particles in their jets. In HSP blazars, the observed X-ray emissi… ▽ More We report the first time measurement of X-ray polarization in the 2$-$8 keV band from the high synchrotron peaked (HSP) blazar 1ES 1959+650. The low energy hump in the broadband spectral energy distribution (SED) of blazars in the log$ν$ versus $νF_ν$ plane is believed to be due to the synchrotron emission process from relativistic particles in their jets. In HSP blazars, the observed X-ray emission lies in the high energy tail of the synchrotron part of the SED, and X-ray polarization measurements can constrain the processes by which particles are accelerated in their jets. We present here our results on the X-ray polarization characteristics of 1ES 1959+650 using observations carried out by the {\it Imaging X-ray Polarimetry Explorer (IXPE)}. {\it IXPE} observed this source for four epochs between May 2022 and August 2023, with observation time ranging from about 50 to 300 ksec. Of the four epochs of observations, we detected polarization on two epochs, significantly larger than the minimum detectable polarization values. From model-independent analysis during the observations on 28 October 2022, in the 2$-$8 keV band, we found the degree of polarization of $Π_X$ = 9.0 $\pm$ 1.6 \% and an electric vector position angle of $ψ_X$ = 53 $\pm$ 5 deg. Similarly, from the observations on 14 August 2023, we found $Π_X$ and $ψ_X$ values as 12.5 $\pm$ 0.7 \% and 20 $\pm$ 2 deg respectively in the 2$-$8 keV band. These values obtained from model-independent analysis are also in agreement with the values obtained from spectro-polarimetric analysis of the I, Q, and U spectra in the 2$-$8 keV range. The measured X-ray polarization is larger than that known in the optical that ranges between 2.5$-$9\% during the period 2008 to 2018. This result can be inferred as the shock acceleration of the particles in the jet of 1ES 1959+650. △ Less

Submitted 3 November, 2023; originally announced November 2023.

Comments: Submitted to PRD

On the properties of X-ray corona in Seyfert 1 galaxies

Authors: Indrani Pal, Anju A., H. Sreehari, Gitika Rameshan, C. S. Stalin, Claudio Ricci, Stefano Marchesi

Abstract: We carried out a uniform and systematic analysis of a sample of 112 nearby bright Seyfert 1 type AGN, the observations of which were carried out by the {\it Nuclear Spectroscopic Telescope Array (NuSTAR)} between August 2013 and May 2022. The main goal of this analysis is to investigate the nature of the X-ray corona in Seyfert 1 galaxies. From the physical model that fits the {\it NuSTAR} spectra… ▽ More We carried out a uniform and systematic analysis of a sample of 112 nearby bright Seyfert 1 type AGN, the observations of which were carried out by the {\it Nuclear Spectroscopic Telescope Array (NuSTAR)} between August 2013 and May 2022. The main goal of this analysis is to investigate the nature of the X-ray corona in Seyfert 1 galaxies. From the physical model that fits the {\it NuSTAR} spectra, we could constrain the high energy cut-off ($\rm{E_{cut}}$) for 73 sources in our sample. For those 73 sources, we fitted the Comptonization model to estimate the temperature ($\rm{kT_{e}}$) of their corona. $\rm{kT_{e}}$ could be constrained in 42 sources. We investigated for possible correlations between various properties of the corona obtained from physical model fits to the observed spectra and between various coronal parameters and physical properties of the sources such as Eddington ratio and black hole mass. We found (a) a strong correlation between $\rm{E_{cut}}$ and the photon index and (b) a significant negative correlation between $\rm{kT_{e}}$ and the optical depth. △ Less

Submitted 27 October, 2023; originally announced October 2023.

Comments: 33 pages, 14 figures, Submitted to ApJ, comments are welcome

Origin of the broadband emission from the transition blazar B2 1308+326

Authors: Ashwani Pandey, Pankaj Kushwaha, Paul J. Wiita, Raj Prince, Bozena Czerny, C. S. Stalin

Abstract: Transition blazars exhibit a shift from one subclass to the next during different flux states. It is therefore crucial to study them to understand the underlying physics of blazars. We probe the origin of the multi-wavelength emission from the transition blazar B2 1308+326 using 14-year-long gamma-ray light curve from Fermi and the quasi-simultaneous data from Swift. We used the Bayesian block alg… ▽ More Transition blazars exhibit a shift from one subclass to the next during different flux states. It is therefore crucial to study them to understand the underlying physics of blazars. We probe the origin of the multi-wavelength emission from the transition blazar B2 1308+326 using 14-year-long gamma-ray light curve from Fermi and the quasi-simultaneous data from Swift. We used the Bayesian block algorithm to identify epochs of flaring and quiescent flux states and modelled the broadband SEDs for these epochs. We employed the one-zone leptonic model in which the synchrotron emission causes the low-energy part of the SED and the high-energy part is produced by the IC emission of external seed photons. We also investigated its multi-band variability properties and gamma-ray flux distribution, and the correlation between optical and gamma-ray emissions. We observed a historically bright flare from B2 1308+326 across the optical to gamma-ray bands in June and July 2022. The highest daily averaged gamma-ray flux was (14.24$\pm$2.36) $\times$ 10$^{-7}$ ph cm$^{-2}$ s$^{-1}$ and was detected on 1 July 2022. The gamma-ray flux distribution was found to be log-normal. The optical and gamma-ray emissions are well correlated with zero time lag. The synchrotron peak frequency changes from $\sim 8 \times$ 10$^{12}$ Hz (in the quiescent state) to $\sim 6 \times$ 10$^{14}$ Hz (in the flaring state), together with a decrease in the Compton dominance providing a hint that the source transitions from a LSP to an ISP. The SEDs for these two states are well-fitted by one-zone leptonic models. The parameters in the model fits are essentially consistent between both SEDs, except for the Doppler-beaming factor, which changes from $\sim$15.6 to $\sim$27 during the transition. An increase in the Doppler factor might cause both the flare and the transition of B2 1308+326 from an LSP to an ISP blazar. △ Less

Submitted 16 November, 2023; v1 submitted 8 October, 2023; originally announced October 2023.

Comments: 10 pages, 7 figures, 2 tables, accepted for publication in the Astronomy & Astrophysics journal

Star formation characteristics of galaxies hosting AGN

Authors: Payel Nandi, C. S. Stalin

Abstract: We present an analysis of the ultra-violet (UV) observations of two Seyfert type active galactic nuclei (AGN), namely NGC$~$4051 and NGC$~$4151. These observations aimed at studying the star formation in the hosts of these AGN were carried out with the ultra-violet imaging telescope on board AstroSat in far-UV. A total of 193 and 328 star-forming regions (SF) were identified using SExtractor in NG… ▽ More We present an analysis of the ultra-violet (UV) observations of two Seyfert type active galactic nuclei (AGN), namely NGC$~$4051 and NGC$~$4151. These observations aimed at studying the star formation in the hosts of these AGN were carried out with the ultra-violet imaging telescope on board AstroSat in far-UV. A total of 193 and 328 star-forming regions (SF) were identified using SExtractor in NGC$~$4051 and NGC$~$4151, respectively. Using aperture photometry of the identified SF regions, we estimated the star formation rates (SFRs). We found NGC$~$4051 to have the lowest SFR with a median value of 3.16 $\times$ 10$^{-5}$ M$_{\odot}$ yr$^{-1}$ while for NGC$~$4151, we found a median SFR of 0.012 M$_{\odot}$ yr$^{-1}$. △ Less

Submitted 17 August, 2023; originally announced August 2023.

Comments: Comments are welcome. Accepted for the 3rd BINA Proceedings

Evidence for low power radio jet-ISM interaction at 10 parsec in the dwarf AGN host NGC 4395

Authors: Payel Nandi, C. S. Stalin, D. J. Saikia, Rogemar A. Riffel, Arijit Manna, Sabyasachi Pal, O. L. Dors, Dominika Wylezalek, Vaidehi S. Paliya, P. Saikia, Pratik Dabhade, Markus-Kissler Patig, Ram Sagar

Abstract: Black hole driven outflows in galaxies hosting active galactic nuclei (AGN) may interact with their interstellar medium (ISM) affecting star formation. Such feedback processes, reminiscent of those seen in massive galaxies, have been reported recently in some dwarf galaxies. However, such studies have usually been on kiloparsec and larger scales and our knowledge on the smallest spatial scales to… ▽ More Black hole driven outflows in galaxies hosting active galactic nuclei (AGN) may interact with their interstellar medium (ISM) affecting star formation. Such feedback processes, reminiscent of those seen in massive galaxies, have been reported recently in some dwarf galaxies. However, such studies have usually been on kiloparsec and larger scales and our knowledge on the smallest spatial scales to which these feedback processes can operate is unclear. Here we demonstrate radio jet$-$ISM interaction on the scale of an asymmetric triple radio structure of $\sim$ 10 parsec size in NGC 4395. This triple radio structure is seen in the 15 GHz continuum image and the two asymmetric jet-like structures are situated on either side of the radio core that coincides with the optical {\it Gaia} position. The high resolution radio image and the extended [OIII]$λ$5007 emission, indicative of an outflow, are spatially coincident and are consistent with the interpretation of a low power radio jet interacting with the ISM. Modelling of the spectral lines using {\tt MAPPINGS}, and estimation of temperature using optical integral field spectroscopic data suggest shock ionization of the gas. The continuum emission at 237 GHz, though weak, was found to spatially coincide with the AGN. However, the CO(2$-$1) line emission was found to be displaced by around 20 parsec northward of the AGN core. The spatial coincidence of molecular H$_2$$λ$2.4085 along the jet direction, the morphology of ionised [OIII]$λ$5007 and displacement of the CO(2$-$1) emission argues for conditions less favourable for star formation in the central $\sim$ 10 parsec region. △ Less

Submitted 14 November, 2023; v1 submitted 25 July, 2023; originally announced July 2023.

Comments: 20 pages, 18 figures, accepted for the publication of The Astrophysical Journal

X-ray polarization observations of IC 4329A with IXPE: Constraining the geometry of the X-ray corona

Authors: Indrani Pal, C. S. Stalin, Rwitika Chatterjee, Vivek K. Agrawal

Abstract: X-ray polarimetry is a powerful tool to probe the geometry of the hot X-ray corona in active galactic nuclei (AGN). Here, we present our results on the characterisation of the X-ray polarization of the radio-quiet Seyfert-type AGN IC 4329A at a redshift of $z$ = 0.016. This is based on observations carried out by the {\it Imaging X-ray Polarimeter (IXPE)}. {\it IXPE} observed IC 4329A on January 5… ▽ More X-ray polarimetry is a powerful tool to probe the geometry of the hot X-ray corona in active galactic nuclei (AGN). Here, we present our results on the characterisation of the X-ray polarization of the radio-quiet Seyfert-type AGN IC 4329A at a redshift of $z$ = 0.016. This is based on observations carried out by the {\it Imaging X-ray Polarimeter (IXPE)}. {\it IXPE} observed IC 4329A on January 5, 2023, for a total observing time of 458 ks. From the model-independent analysis, we found a polarization degree ($Π_{X}$) of 3.7$\pm$1.5$\%$ and a polarization position angle ($Ψ_{X}$) of 61$^{\circ}$$\pm$12$^{\circ}$ in the 2$-$8 keV energy range (at 1$σ$ confidence level). This is also in agreement with the values of $Π_{X}$ and $Ψ_{X}$ of 4.7$\pm$2.2$\%$ and 71$^{\circ}$ $\pm$14$^{\circ}$ respectively obtained from spectro-polarimetric analysis of the I, Q and U Stokes spectra in the 2$-$8 keV energy band (at the 90$\%$ confidence). The value of $Π_X$ in the 2-8 keV band obtained from the model-independent analysis is lower than the minimum detectable polarization (MDP) value of 4.5$\%$. However, $Π_X$ obtained from spectro-polarimetric analysis in the 2-8 keV band is larger than the MDP value. In the 3-5 keV band, we found $Π_X$ of 6.5 $\pm$ 1.8, which is larger than the MDP value of 5.5$\%$. The observed moderate value of $Π_{X}$ obtained from the analysis of the {\it IXPE} data in the 3$-$5 keV band argues against a spherical lamp$-$post geometry for the X-ray corona in IC 4329A; however, considering simulations, the observed polarization measurements tend to favour a conical shape geometry for the corona. This is the first time measurement of X-ray polarization in IC 4329A. Measurements of the X-ray polarization in many such radio-quiet AGN will help in constraining the geometry of the X-ray corona in AGN. △ Less

Submitted 16 May, 2023; originally announced May 2023.

Comments: 7 pages, 4 figures, Submitted to the Journal of Astronomy and Astrophysics

Star formation in the dwarf Seyfert galaxy NGC 4395: Evidence for both AGN and SNe feedback?

Authors: Payel Nandi, C. S. Stalin, D. J. Saikia, S. Muneer, George Mountrichas, Dominika Wylezalek, R. Sagar, Markus Kissler-Patig

Abstract: We present a detailed multi-wavelength study of star formation in the dwarf galaxy NGC 4395 which hosts an active galactic nucleus (AGN). From our observations with the Ultra-Violet Imaging Telescope, we have compiled a catalogue of 284 star forming (SF) regions, out of which we could detect 120 SF regions in H$α$ observations. Across the entire galaxy, we found the extinction corrected star forma… ▽ More We present a detailed multi-wavelength study of star formation in the dwarf galaxy NGC 4395 which hosts an active galactic nucleus (AGN). From our observations with the Ultra-Violet Imaging Telescope, we have compiled a catalogue of 284 star forming (SF) regions, out of which we could detect 120 SF regions in H$α$ observations. Across the entire galaxy, we found the extinction corrected star formation rate (SFR) in the far ultra-violet (FUV) to range from 2.0 $\times$ 10$^{-5}$ M$_\odot$yr$^{-1}$ to 1.5 $\times$ 10$^{-2}$ M$_\odot$yr$^{-1}$ with a median of 3.0 $\times$ 10$^{-4}$ M$_\odot$yr$^{-1}$ and age to lie in the range of $\sim$ 1 to 98 Myr with a median of 14 Myr. In H$α$ we found the SFR to range from 7.2 $\times$ 10$^{-6}$ M$_\odot$yr$^{-1}$ to 2.7 $\times$ 10$^{-2}$ M$_\odot$yr$^{-1}$ with a median of 1.7 $\times$ 10$^{-4}$ M$_\odot$yr$^{-1}$ and age to lie between 3 to 6 Myr with a median of 5 Myr. The stellar ages derived from H$α$ show a gradual decline with galactocentric distance. We found three SF regions close to the center of NGC~4395 with high SFR both from H$α$ and UV which could be attributed to feedback effects from the AGN. We also found six other SF regions in one of the spiral arms having higher SFR. These are very close to supernovae remnants which could have enhanced the SFR locally. We obtained a specific SFR (SFR per unit mass) for the whole galaxy 4.64 $\times$ 10$^{-10}$ yr$^{-1}$. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: 10 pages, 11 figures, Accepted for publication in ApJ

Multiwavelength study of radio galaxy Pictor A: detection of western hotspot in far-UV and possible origin of high energy emissions

Authors: Sanna Gulati, Debbijoy Bhattacharya, M. C. Ramadevi, C. S. Stalin, P. Sreekumar

Abstract: A comprehensive study of the nucleus and western hotspot of Pictor A is carried out using AstroSat observations, 13 years of Fermi, and archival Swift observations, along with other published data. We report the first detection of the western hotspot of Pictor A in the far-UV band using observations from AstroSat-UVIT. The broad-band SED of the western hotspot is explained by a multizone emission… ▽ More A comprehensive study of the nucleus and western hotspot of Pictor A is carried out using AstroSat observations, 13 years of Fermi, and archival Swift observations, along with other published data. We report the first detection of the western hotspot of Pictor A in the far-UV band using observations from AstroSat-UVIT. The broad-band SED of the western hotspot is explained by a multizone emission scenario where X-ray emission is caused by synchrotron emission process in the substructures embedded in the diffuse region, while the emission in radio to optical is caused by synchrotron emission process in the diffuse region. We do not notice any excess in the IR band and an additional zone (beyond 2-zone) is not required to account for the X-ray emission. Our broad-band spectro-temporal study and associated modelling of the core and hotspot of Pictor A suggests that (a) gamma-rays originate in the nuclear jet and not from the hotspot (b) X-ray emission from the core of Pictor A has nuclear jet-origin instead of previously reported disk-origin. △ Less

Submitted 8 March, 2023; originally announced March 2023.

Comments: 11 pages, 4 figures, Accepted for publication in MNRAS

An empirical relation to estimate host galaxy stellar light from AGN spectra

Authors: Priyanka Jalan, Suvendu Rakshit, Jong-Hak Woo, Jari Kotilainen, C. S. Stalin

Abstract: Measurement of black hole mass for low-$z$ ($z\leq 0.8$) Active Galactic Nuclei (AGNs) is difficult due to the strong contribution from host galaxy stellar light necessitating detailed spectral decomposition to estimate the AGN luminosity. Here, we present an empirical relation to estimate host galaxy stellar luminosity from the optical spectra of AGNs at $z\leq 0.8$. The spectral data were select… ▽ More Measurement of black hole mass for low-$z$ ($z\leq 0.8$) Active Galactic Nuclei (AGNs) is difficult due to the strong contribution from host galaxy stellar light necessitating detailed spectral decomposition to estimate the AGN luminosity. Here, we present an empirical relation to estimate host galaxy stellar luminosity from the optical spectra of AGNs at $z\leq 0.8$. The spectral data were selected from the fourteenth data release of the Sloan Digital Sky Survey (SDSS-DR14) quasar catalog having a signal-to-noise ratio at 5100 Å(SNR$_{5100}$) $>$10 containing 11415 quasars. The median total luminosity (log ($L_\text{total}$/[erg s$^{-1}$])), stellar luminosity (log ($L_\text{star}$/[erg s$^{-1}$])), and AGN continuum luminosity (log ($L_\text{cont}$/[erg s$^{-1}$])) in our sample are 44.52, 44.06, and 44.30, respectively. We fit the AGN power-law continuum, host galaxy, and iron blend contribution, simultaneously over the entire available spectrum. We found the host galaxy fraction to anti-correlate with continuum luminosity and can be well-represented by a polynomial function, which can be used to correct the stellar light contribution from AGN spectra. We also found anti-correlation between host galaxy fraction and iron strength, Eddington ratio, and redshift. The empirical relation gives comparable results of host-fraction with the image decomposition method. △ Less

Submitted 28 March, 2023; v1 submitted 3 February, 2023; originally announced February 2023.

Comments: 6 pages, 6 figures, published in MNRASL

Journal ref: MNRAS; 521; L11-L16; (2023)

UVIT observations of the Small Magellanic Cloud: Point source catalogue

Authors: Ashish Devaraj, Prajwel Joseph, C. S. Stalin, Shyam N. Tandon, Swarna K. Ghosh

Abstract: Three fields in the outskirts of the Small Magellanic Cloud were observed by the Ultra-Violet Imaging Telescope (UVIT) on board AstroSat, during 31 December 2017 and 01 January 2018. The observations were carried out on a total of seven filters, three in the far ultra-violet (FUV; 1300$-$1800 Å) band and four in the near ultra-violet (NUV; 2000$-$3000 Å) band. We carried out photometry of these ob… ▽ More Three fields in the outskirts of the Small Magellanic Cloud were observed by the Ultra-Violet Imaging Telescope (UVIT) on board AstroSat, during 31 December 2017 and 01 January 2018. The observations were carried out on a total of seven filters, three in the far ultra-violet (FUV; 1300$-$1800 Å) band and four in the near ultra-violet (NUV; 2000$-$3000 Å) band. We carried out photometry of these observations that have a spatial resolution better than 1.5$^{\prime\prime}$. We present here the first results of this work, which is a matched catalogue of 11,241 sources detected in three FUV and four NUV wavelengths. We make the catalogue available online, which would be of use to the astronomical community to address a wide variety of astrophysical problems. We provide an expression to estimate the total count rate in the full point spread function of UVIT that also incorporate the effect of saturation. △ Less

Submitted 2 February, 2023; originally announced February 2023.

TXS 1433+205: The most distant gamma-ray emitting FR II radio galaxy

Authors: Vaidehi S. Paliya, D. J. Saikia, C. S. Stalin

Abstract: The orientation of the jet axis to the line of sight of the observer plays a major role in explaining the phenomena observed from blazars and radio galaxies. In the gamma-ray band, only a handful of radio galaxies have been identified, all being located in the nearby Universe (z<0.5). Here we report the identification of 4FGL J1435.5+2021, associated with TXS 1433+205, as a Fanaroff-Riley type II… ▽ More The orientation of the jet axis to the line of sight of the observer plays a major role in explaining the phenomena observed from blazars and radio galaxies. In the gamma-ray band, only a handful of radio galaxies have been identified, all being located in the nearby Universe (z<0.5). Here we report the identification of 4FGL J1435.5+2021, associated with TXS 1433+205, as a Fanaroff-Riley type II (FR II) radio galaxy at a considerably higher redshift of z=0.748, thereby making it the most distant gamma-ray detected radio galaxy known as of now. The Very Large Array Sky Survey data at 3 GHz resolves the source morphology into a bright core, a jet and two hotspots, with a total end-to-end projected length between lobe extremities of ~170 kpc. The optical and radio properties of this enigmatic object suggest it to be a high-excitation FR II radio galaxy. The multi-wavelength behaviour of TXS 1433+205 is found to be similar to other gamma-ray detected FR II sources but is at the high luminosity end. We suggest that the ongoing and upcoming high-resolution radio surveys will lead to the identification of many more high-redshift radio galaxies in the gamma-ray sky, thus allowing comprehensive studies of misaligned relativistic jets. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Comments: MNRAS Letters, in press

Search for coronal temperature variation in Seyfert galaxies

Authors: Indrani Pal, C. S. Stalin

Abstract: While the temperature of the X-ray corona ($\rm{kT_e}$) in active galactic nuclei (AGN) are known for many sources, its variation, if any, is limited to a handful of objects. This is in part due to the requirement of good signal-to-noise X-ray spectra covering a wide range of energies. We present here results on the X-ray spectral analysis of 18 Seyferts, having more than one epoch of observations… ▽ More While the temperature of the X-ray corona ($\rm{kT_e}$) in active galactic nuclei (AGN) are known for many sources, its variation, if any, is limited to a handful of objects. This is in part due to the requirement of good signal-to-noise X-ray spectra covering a wide range of energies. We present here results on the X-ray spectral analysis of 18 Seyferts, having more than one epoch of observations to look for variation in $\rm{kT_e}$. The data for a total of 52 epochs on these 18 AGN were taken from observations carried out by NuSTAR in the 3$-$79 keV energy band. From phenomenological and physical model fits to the multi-epoch data on these 18 sources from {\it NuSTAR}, we could constrain the cut-off energy ($E_{cut}$) in a large fraction of the sources. Also, from Comptonized model fits, we could obtain $\rm{kT_e}$ for our sample. Of the 18 sources, at the 90 per cent confidence level, evidence for variation in $\rm{kT_e}$ was found for only one source, namely MCG+08-11-011. For this source, between two epochs, separated by about five years, we found $\rm{kT_e}$ to decrease from 57$^{+29}_{-16}$ keV to 30$^{+11}_{-7}$ keV. During the same period, the flux decreased from (12.60 to 14.02) $\times$ 10$^{-11}$ erg cm$^{-2}$ s$^{-1}$ and the optical depth increased from 1.68 to 2.73. We thus found a positive correlation between flux and coronal temperature with a reduction of about 40 per cent in optical depth. Our observations tend to favour the vertically outflowing corona scenario for the observed variation in $\rm{kT_e}$ in MCG+08-11-011. △ Less

Submitted 2 November, 2022; originally announced November 2022.

Comments: 17 pages, 7 figuers, accepted for publication in MNRAS-MJ

RAD@home citizen science discovery of an AGN spewing a large unipolar radio bubble onto its merging companion galaxy

Authors: Ananda Hota, Pratik Dabhade, Sravani Vaddi, Chiranjib Konar, Sabyasachi Pal, Mamta Gulati, C S. Stalin, Ck Avinash, Avinash Kumar, Megha Rajoria, Arundhati Purohit

Abstract: AGN feedback during galaxy merger has been the most favoured model to explain black hole-galaxy co-evolution. However, how the AGN-driven jet/wind/radiation is coupled with the gas of the merging galaxies, which leads to positive feedback, momentarily enhanced star formation, and subsequently negative feedback, a decline in star formation, is poorly understood. Only a few cases are known where the… ▽ More AGN feedback during galaxy merger has been the most favoured model to explain black hole-galaxy co-evolution. However, how the AGN-driven jet/wind/radiation is coupled with the gas of the merging galaxies, which leads to positive feedback, momentarily enhanced star formation, and subsequently negative feedback, a decline in star formation, is poorly understood. Only a few cases are known where the jet and companion galaxy interaction leads to minor off-axis distortions in the jets and enhanced star formation in the gas-rich minor companions. Here, we briefly report one extraordinary case, RAD12, discovered by RAD@home citizen science collaboratory, where for the first time a radio jet-driven bubble ~137 kpc is showing a symmetric reflection after hitting the incoming galaxy which is not a gas-rich minor but a gas-poor early-type galaxy in a major merger. Surprisingly, neither positive feedback nor any radio lobe on the counter jet side, if any, is detected. It is puzzling if RAD12 is a genuine one-sided jet or a case of radio lobe trapped, compressed and re-accelerated by shocks during the merger. This is the first imaging study of RAD12 presenting follow-up with the GMRT, archival MeerKAT radio data and CFHT optical data. △ Less

Submitted 12 October, 2022; originally announced October 2022.

Comments: Published in the MNRAS Letters. Contains 5 pages and 2 figures. Comments and Collaborations for follow up of this #RADatHomeIndia #CitizenScience discovery, most welcome

Journal ref: MNRAS 2022

Detection of minute-timescale gamma-ray variability in BL Lacertae by Fermi-LAT

Authors: Ashwani Pandey, C. S. Stalin

Abstract: BL Lacertae, the prototype of the BL Lacertae (BL Lac) category of blazars, underwent a giant $γ-$ray flare in April 2021. The Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope (hereafter Fermi-LAT) observed a peak $γ-$ray (0.1$-$500 GeV) flux of $\sim$2 $\times$ 10$^{-5}$ photons cm$^{-2}$ s$^{-1}$ within a single orbit on 2021 April 27, which is historically the brightest… ▽ More BL Lacertae, the prototype of the BL Lacertae (BL Lac) category of blazars, underwent a giant $γ-$ray flare in April 2021. The Large Area Telescope (LAT) onboard the Fermi Gamma-ray Space Telescope (hereafter Fermi-LAT) observed a peak $γ-$ray (0.1$-$500 GeV) flux of $\sim$2 $\times$ 10$^{-5}$ photons cm$^{-2}$ s$^{-1}$ within a single orbit on 2021 April 27, which is historically the brightest $γ-$ray flux ever detected from the source. Here, we report, for the first time, the detection of significant minute-timescale GeV $γ-$ray flux variability in the BL Lac subclass of blazars by the Fermi-LAT. We resolved the source variability down to 2-min binned timescales with a flux halving time of $\sim$1 minute, which is the shortest GeV variability timescale ever observed from blazars. The detected variability timescale is much shorter than the light-crossing time ($\sim 14$ minutes) across the central black hole of BL Lac indicating a very compact $γ-$ray emission site within the outflowing jet. Such a compact emitting region requires the bulk Lorentz factor of the jet to be larger than 16 so that the jet power is not super Eddington. We found a minimum Doppler factor $δ_{min}$ of 15 using the $δ$ function approximation for the $γγ$ opacity constraint. For a conical jet geometry, considering $Γ= δ_{min}$, the observed short variability timescale suggests the very compact emission region to lie at a distance of about 8.62 $\times$ 10$^{14}$ cm from the central engine of BL Lac. △ Less

Submitted 3 October, 2022; originally announced October 2022.

Comments: 6 pages, 5 figures, 2 tables, accepted for publication in the Astronomy & Astrophysics (A&A) main journal

Journal ref: A&A 668, A152 (2022)

X-ray spectral and timing analysis of the Compton Thick Seyfert 2 galaxy NGC 1068

Authors: Indrani Pal, C. S. Stalin, M. L. Parker, Vivek. K. Agrawal, S. Marchesi

Abstract: We present the timing and spectral analysis of the Compton Thick Seyfert 2 active galactic nuclei NGC 1068 observed using {\it NuSTAR} and {\it XMM-Newton}. In this work for the first time we calculated the coronal temperature ($\rm{kT_{e}}$) of the source and checked for its variation between the epochs if any. The data analysed in this work comprised of (a) eight epochs of observations with {\it… ▽ More We present the timing and spectral analysis of the Compton Thick Seyfert 2 active galactic nuclei NGC 1068 observed using {\it NuSTAR} and {\it XMM-Newton}. In this work for the first time we calculated the coronal temperature ($\rm{kT_{e}}$) of the source and checked for its variation between the epochs if any. The data analysed in this work comprised of (a) eight epochs of observations with {\it NuSTAR} carried out during the period December 2012 to November 2017, and, (b) six epochs of observations with {\it XMM-Newton} carried out during July 2000 to February 2015. From timing analysis of the {\it NuSTAR} observations, we found the source not to show any variations in the soft band. However, on examination of the flux at energies beyond 20 keV, during August 2014 and August 2017 the source was brighter by about 20\% and 30\% respectively compared to the mean flux of the three 2012 {\it NuSTAR} observations as in agreement with earlier results in literature. From an analysis of {\it XMM-Newton} data we found no variation in the hard band (2 $-$ 4 keV) between epochs as well as within epochs. In the soft band (0.2 $-$ 2 keV), while the source was found to be not variable within epochs, it was found to be brighter in epoch B relative to epoch A. By fitting physical models we determined $\rm{kT_{e}}$ to range between 8.46$^{+0.39}_{-0.66}$ keV and 9.13$^{+0.63}_{-0.98}$ keV. From our analysis, we conclude that we found no variation of $\rm{kT_{e}}$ in the source. △ Less

Submitted 21 September, 2022; originally announced September 2022.

Comments: 13 pages, 14 figures, accepted for publication in MNRAS

Study of correlation between optical flux and polarization variations in BL Lac objects

Authors: Bhoomika Rajput, Ashwani Pandey, C. S. Stalin, Blesson Mathew

Abstract: Polarized radiation from blazars is one key piece of evidence for synchrotron radiation at low energy, which also shows variations. We present here our results on the correlation analysis between optical flux and polarization degree (PD) variations in a sample of 11 BL Lac objects using $\sim$ 10 years of data from the Steward Observatory. We carried out the analysis on long-term ($\sim$ several m… ▽ More Polarized radiation from blazars is one key piece of evidence for synchrotron radiation at low energy, which also shows variations. We present here our results on the correlation analysis between optical flux and polarization degree (PD) variations in a sample of 11 BL Lac objects using $\sim$ 10 years of data from the Steward Observatory. We carried out the analysis on long-term ($\sim$ several months) as well as on short-term timescales ($\sim$ several days). On long-term timescales, for about 85% of the observing cycles, we found no correlation between optical flux and PD. On short-term timescales, we found a total of 58 epochs with a significant correlation between optical flux and PD, where both positive and negative correlation were observed. In addition, we also found a significant correlation between optical flux and $γ$-ray flux variations on long-term timescales in 11% of the observing cycles. The observed PD variations in our study cannot be explained by changes in the power-law spectral index of the relativistic electrons in the jets. The shock-in-jet scenario is favoured for the correlation between optical flux and PD, whereas the anti-correlation can be explained by the presence of multi-zone emission regions. The varying correlated behaviour can also be explained by the enhanced optical flux caused by the newly developed radio knots in the jets and their magnetic field alignment with the large scale jet magnetic field. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: 20 pages, 6 figures, accepted for publication in MNRAS

Spectroscopic reverberation mapping of Quasar PKS 0736+017: Broad-Line Region and Black-hole Mass

Authors: Shivangi Pandey, Suvendu Rakshit, Jong-Hak Woo, C. S. Stalin

Abstract: To understand the mass distribution and co-evolution of supermassive black holes with their host galaxy, it is crucial to measure the black hole mass of AGN. Reverberation mapping is a unique tool to estimate the black hole masses in AGN. We performed spectroscopic reverberation study using long-term monitoring data with more than 100 spectra of a radio-loud quasar PKS 0736+017 to estimate the siz… ▽ More To understand the mass distribution and co-evolution of supermassive black holes with their host galaxy, it is crucial to measure the black hole mass of AGN. Reverberation mapping is a unique tool to estimate the black hole masses in AGN. We performed spectroscopic reverberation study using long-term monitoring data with more than 100 spectra of a radio-loud quasar PKS 0736+017 to estimate the size of the broad line region (BLR) and black hole mass. The optical spectrum shows strong H$\mathrmβ$ and H$\mathrmγ$ emission lines. We generated the light curves of 5100Å continuum flux ($f_{5100}$), H$β$, and H$γ$. All the light curves are found to be strongly variable with fractional variability of 69$\%$, 21$\%$, 30$\%$ for V-band, H$β$, and H$γ$ light curves, respectively. Along with the thermal contribution, non-thermal emission contributes to the estimated continuum luminosity at 5100Å. Using different methods, e.g., CCF, {\small JAVELIN}, von-neumann, we estimated the size of the BLR, which is found to be 66.4$^{+6.0}_{-4.2}$ light days in the rest frame. The BLR size combined with the line width of H$β$ provides a black hole mass of 7.32$^{+0.89}_{-0.91} \times 10^{7}M_{\odot}$. The source closely follows the BLR size-luminosity relation of AGN. △ Less

Submitted 4 September, 2022; originally announced September 2022.

Comments: published in MNRAS, 13 pages, 13 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, 2022

UVIT view of Centaurus A; a detailed study on positive AGN feedback

Authors: Prajwel Joseph, P. Sreekumar, C. S. Stalin, K. T. Paul, Chayan Mondal, Koshy George, Blesson Mathew

Abstract: $\require{mediawiki-texvc}… ▽ More $\require{mediawiki-texvc}$ Supermassive black holes at the centre of active galactic nuclei (AGN) produce relativistic jets that can affect the star formation characteristics of the AGN hosts. Observations in the ultraviolet (UV) band can provide an excellent view of the effect of AGN jets on star formation. Here, we present a census of star formation properties in the Northern Star-forming Region (NSR) that spans about 20 kpc of the large radio source Centaurus A hosted by the giant elliptical galaxy NGC 5128. In this region, we identified 352 UV sources associated with Cen A using new observations at an angular resolution of $<$1.5 arcseconds observed with the Ultra-Violet Imaging Telescope (UVIT) onboard AstroSat. These observations were carried out in one far-ultraviolet (FUV; $λ_{\text{mean}}$ = 1481 $Å$) and three near-ultraviolet (NUV; with $λ_{\text{mean}}$ of 2196 $Å$, 2447 $Å$, and 2792 $Å$, respectively) bands. The star-forming sources identified in UV tend to lie in the direction of the jet of Cen A, thereby suggesting jet triggering of star formation. Separating the NSR into Outer and Inner regions, we found the stars in the Inner region to have a relatively younger age than the Outer region, suggesting that the two regions may have different star formation histories. We also provide the UVIT source catalogue in the NSR. △ Less

Submitted 22 August, 2022; originally announced August 2022.

Comments: Accepted for publication by MNRAS. Data and Jupyter notebooks are available at https://github.com/prajwel/Centaurus-A_NSR . 15 pages, 19 figures, 6 tables

A multiwavelength study of the flat spectrum radio-quasar NVSS J141922-083830 covering four flaring episodes

Authors: D. A. H. Buckley, R. J. Britto, S. Chandra, V. Krushinsky, M. Böttcher, S. Razzaque, V. Lipunov, C. S. Stalin, E. Gorbovskoy, N. Tiurina, D. Vlasenko, A. Kniazev

Abstract: We present multiwavelength observations and a model for flat spectrum radio quasar NVSS J141922-083830, originally classified as a blazar candidate of unknown type (BCU II object) in the Third Fermi-LAT AGN Catalog (3LAC). Relatively bright flares (>3 magnitudes) were observed on 21 February 2015 (MJD 57074) and 8 September 2018 (MJD 58369) in the optical band with the MASTER Global Robotic Net (M… ▽ More We present multiwavelength observations and a model for flat spectrum radio quasar NVSS J141922-083830, originally classified as a blazar candidate of unknown type (BCU II object) in the Third Fermi-LAT AGN Catalog (3LAC). Relatively bright flares (>3 magnitudes) were observed on 21 February 2015 (MJD 57074) and 8 September 2018 (MJD 58369) in the optical band with the MASTER Global Robotic Net (MASTER-Net) telescopes. Optical spectra obtained with the Southern African Large Telescope (SALT) on 1 March 2015 (MJD 57082), during outburst, and on 30 May 2017 (MJD 57903), during quiescence, showed emission lines at 5325Å and at $\approx$3630Å that we identified as the Mg II 2798Å and C III] 1909Å lines, respectively, and hence derived a redshift z = 0.903. Analysis of Fermi-LAT data was performed in the quiescent regime (5 years of data) and during four prominent flaring states in February-April 2014, October-November 2014, February-March 2015 and September 2018. We present spectral and timing analysis with Fermi-LAT. We report a hardening of the gamma-ray spectrum during the last three flaring periods, with a power-law spectral index $Γ= 2.0$-$2.1$. The maximum gamma-ray flux level was observed on 24 October 2014 (MJD 56954) at $(7.57 \pm 1.83) \times 10^{-7}$ ph~cm$^{-2}$s$^{-1}$. The multi-wavelength spectral energy distribution during the February-March 2015 flare supports the earlier evidence of this blazar to belong to the FSRQ class. The SED can be well represented with a single-zone leptonic model with parameters typical of FSRQs, but also a hadronic origin of the high-energy emission can not be ruled out. △ Less

Submitted 29 July, 2022; originally announced July 2022.

Comments: 15 pages, 11 figures, 8 tables, accepted for publication in MNRAS (29 July 2022)

On the origin of core radio emissions from black hole sources in the realm of relativistic shocked accretion flow

Authors: Santabrata Das, Anuj Nandi, C. S. Stalin, Suvendu Rakshit, Indu Kalpa Dihingia, Swapnil Singh, Ramiz Aktar, Samik Mitra

Abstract: We study the relativistic, inviscid, advective accretion flow around the black holes and investigate a key feature of the accretion flow, namely the shock waves. We observe that the shock-induced accretion solutions are prevalent and such solutions are commonly obtained for a wide range of the flow parameters, such as energy (${\cal E}$) and angular momentum ($λ$), around the black holes of spin v… ▽ More We study the relativistic, inviscid, advective accretion flow around the black holes and investigate a key feature of the accretion flow, namely the shock waves. We observe that the shock-induced accretion solutions are prevalent and such solutions are commonly obtained for a wide range of the flow parameters, such as energy (${\cal E}$) and angular momentum ($λ$), around the black holes of spin value $0\le a_{\rm k} < 1$. When the shock is dissipative in nature, a part of the accretion energy is released through the upper and lower surfaces of the disc at the location of the shock transition. We find that the maximum accretion energies that can be extracted at the dissipative shock ($Δ{\cal E}^{\rm max}$) are $\sim 1\%$ and $\sim 4.4\%$ for Schwarzschild black holes ($a_{\rm k}\rightarrow 0$) and Kerr black holes ($a_{\rm k}\rightarrow 1$), respectively. Using $Δ{\cal E}^{\rm max}$, we compute the loss of kinetic power (equivalently shock luminosity, $L_{\rm shock}$) that is enabled to comply with the energy budget for generating jets/outflows from the jet base ($i.e.$, post-shock flow). We compare $L_{\rm shock}$ with the observed core radio luminosity ($L_R$) of black hole sources for a wide mass range spanning $10$ orders of magnitude with sub-Eddington accretion rate and perceive that the present formalism seems to be potentially viable to account $L_R$ of $16$ Galactic black hole X-ray binaries (BH-XRBs) and $2176$ active galactic nuclei (AGNs). We further aim to address the core radio luminosity of intermediate-mass black hole (IMBH) sources and indicate that the present model formalism perhaps adequate to explain core radio emission of IMBH sources in the sub-Eddington accretion limit. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 15 pages, 4 figures, and 3 tables; accepted for publication in MNRAS

Variable mass accretion and failed wind explain changing look phenomena in NGC 1365

Authors: Santanu Mondal, Tek P. Adhikari, Krzysztof Hryniewicz, C. S. Stalin, Ashwani Pandey

Abstract: Changing look active galactic nuclei (CLAGNs) show complex nature in their X-ray spectral shape and line of sight column density variation. The physical mechanisms responsible for these variations are unclear. Here, we study the spectral properties of a CLAGN, NGC\,1365 using combined {\it XMM-Newton} and {\it NuSTAR} observations to understand the CL behavior. The model fitted mass accretion rate… ▽ More Changing look active galactic nuclei (CLAGNs) show complex nature in their X-ray spectral shape and line of sight column density variation. The physical mechanisms responsible for these variations are unclear. Here, we study the spectral properties of a CLAGN, NGC\,1365 using combined {\it XMM-Newton} and {\it NuSTAR} observations to understand the CL behavior. The model fitted mass accretion rate varied between $0.003\pm 0.001$ and $0.009\pm0.002$ $\dot M_{\rm Edd}$ and the dynamic corona changed from $28\pm 3$ to $10\pm1$ $r_g$. We found that the variable absorption column density correlates with the mass accretion rate and the geometry of the corona. The derived wind velocity was sufficiently low compared to the escape velocity to drive the wind away from the disc for the epochs when column densities were high. This suggests that the high and variable absorption can be due to failed winds from the disc. Our estimated ratio of mass outflow to inflow rate from the inner region of the disc lies between $0.019\pm0.006$ and $0.12\pm0.04$. From spectral fitting of the combined data, we found the mass of the central black hole to be constant $4.38\pm0.34 - 4.51\pm0.29 \times10^{6} M_\odot$, consistent with earlier findings. The confidence contours of $N_H$ with other model parameters show that the model fitted parameters are robust and non-degenerate. Our study construed that the changing accretion rate, which is a fundamental physical quantity and the geometry of the corona driving the CL phenomena in NGC\,1365. The physical picture considered in this work connects both variable continuum and variable absorbing medium scenarios. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 8 pages, 7 figures, 1 table, accepted for publication in A&A, comments welcome

Journal ref: A&A 662, A77 (2022)

Flux and spectral variability of Mrk 421 during its moderate activity state using NuSTAR: Possible accretion disc contribution?

Authors: Santanu Mondal, Priyanka Rani, C. S. Stalin, Sandip K. Chakrabarti, Suvendu Rakshit

Abstract: The X-ray emission in BL Lac objects is believed to be dominated by synchrotron emission from their relativistic jets. However, when the jet emission is not strong, one could expect signatures of X-ray emission from inverse Compton scattering of accretion disc photons by the corona. Moreover, the observed X-ray variability can also originate in the disc, and gets propagated and amplified by the je… ▽ More The X-ray emission in BL Lac objects is believed to be dominated by synchrotron emission from their relativistic jets. However, when the jet emission is not strong, one could expect signatures of X-ray emission from inverse Compton scattering of accretion disc photons by the corona. Moreover, the observed X-ray variability can also originate in the disc, and gets propagated and amplified by the jet. Here, we present results on the BL Lac object Mrk 421 using the NuSTAR data acquired during 2017 when the source was in a moderate X-ray brightness state. For comparison with high jet activity, we also considered one epoch data in April 2013 during a very high X-ray brightness state. Our aim is to explore the possibility of the signature of accretion disc emission in the overall X-ray emission from Mrk 421. The spectral fitting of the data using the two component advective flow model shows gives (a) the size of the dynamic corona at the base of the jet from ~28 to 10 r$_s$, (b) the disc mass accretion rate from 0.021 to 0.051 $\dot M_{\rm Edd}$, (c) the halo mass accretion rate from 0.22 to 0.35 $\dot M_{\rm Edd}$, and (d) the viscosity parameter of the Keplerian accretion disc from 0.18 $-$ 0.25. In the assumed model, the total flux, disc and jet flux correlate with the radio flux observed during these epochs. We conclude that the spectra of all the epochs of Mrk 421 in 2017 are well described by the accretion disc based two component advective flow model. The estimated disc and jet flux relations with radio flux show that accretion disc can contribute to the observed X-ray emission, when X-ray data (that covers a small portion of the broad band spectral energy distribution of Mrk 421) is considered in isolation. However, the present disc based models are disfavoured with respect to the relativistic jet models when considering the X-ray data in conjunction with data at other wavelengths. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 15 pages, 9 figures, 9 tables, accepted for publication in A&A, comments welcome

Journal ref: A&A 663, A178 (2022)

An Automated Pipeline for Ultra-Violet Imaging Telescope (UVIT)

Authors: S. K. Ghosh, S. N. Tandon, S. K. Singh, D. S. Shelat, P. Tahlani, A. K. Singh, T. P. Srinivasan, P. Joseph, A. Devaraj, K. George, R. Mohan, J. Postma, C. S. Stalin

Abstract: We describe a versatile pipeline for processing the data collected by the Ultra-Violet Imaging Telescope (UVIT) on board Indian Multi-wavelength astronomical satellite AstroSat.The UVIT instrument carries out simultaneous astronomical imaging through selected filters / gratings in Far-Ultra-Violet (FUV), Near-Ultra-Violet & visible (VIS) bands of the targeted circular sky field (~ 0.5 deg dia). Th… ▽ More We describe a versatile pipeline for processing the data collected by the Ultra-Violet Imaging Telescope (UVIT) on board Indian Multi-wavelength astronomical satellite AstroSat.The UVIT instrument carries out simultaneous astronomical imaging through selected filters / gratings in Far-Ultra-Violet (FUV), Near-Ultra-Violet & visible (VIS) bands of the targeted circular sky field (~ 0.5 deg dia). This pipeline converts the data (Level-1) emanating from UVIT in their raw primitive format supplemented by inputs from the spacecraft sub-systems into UV sky images (& slitless grating spectra) and associated products readily usable by astronomers (Level-2). The primary products include maps of Intensity (rate of photon arrival), error on Intensity and effective Exposure. The pipeline is open source, extensively user configurable with many selectable parameters and its execution is fully automated. The key ingredients of the pipeline includes - extraction of drift in pointing of the spacecraft, and disturbances in pointing due to internal movements; application of various corrections to measured position in the detector for each photon - e.g. differential pointing with respect to a reference frame for shift and add operation, systematic effects and artifacts in the optics of the telescopes and detectors, exposure tracking on the sky, alignment of sky products from multi-episode exposures to generate a consolidated set and astrometry. Detailed logs of operations and intermediate products for every processing stage are accessible via user selectable options. While large number of selectable parameters are available for the user, a well characterized standard default set is used for executing this pipeline at the Payload Operation Centre (POC) for UVIT and selected products are archived and disseminated by the Indian Space Research Organization (ISRO) through its ISSDC portal. △ Less

Submitted 15 March, 2022; originally announced March 2022.

Comments: Accepted for publication in Journal of Astrophysics & Astronomy, 50 pages, 16 figures

NuSTAR spectral analysis of three Seyfert galaxies: NGC 3227, NGC 5548 and MR 2251$-$178

Authors: Indrani Pal, C. S. Stalin, L. Mallick, Priyanka Rani

Abstract: The observed nuclear X-ray emission in the radio-quiet category of active galactic nuclei (AGN) is believed to be from a compact region, the corona situated in the vicinity of the central supermassive black holes (SMBH). The shape of the X-ray continuum, among other factors, depends on the temperature of the corona ($\rm{kT_{e}}$). The launch of the Nuclear Spectroscopic Telescope Array (NuSTAR) h… ▽ More The observed nuclear X-ray emission in the radio-quiet category of active galactic nuclei (AGN) is believed to be from a compact region, the corona situated in the vicinity of the central supermassive black holes (SMBH). The shape of the X-ray continuum, among other factors, depends on the temperature of the corona ($\rm{kT_{e}}$). The launch of the Nuclear Spectroscopic Telescope Array (NuSTAR) has led to the determination of the high energy cut-off ($\rm{E_{cut}}$; and thereby $\rm{kT_{e}}$) in many AGN. In a handful of sources, multiple observations with NuSTAR have also revealed changes in $\rm{E_{cut}}$. In this work we aimed to investigate the variation in $\rm{kT_{e}}$ in three AGN, namely NGC 3227, NGC 5548 and MR 2251$-$178 using more than one epoch of data on a source from NuSTAR. We carried out spectral analysis of multiple epochs of data acquired using NuSTAR on the three sources including a few new observations not published so far. By fitting Comptonization model to the data we determined the temperature of the corona and also investigated changes in $\rm{kT_{e}}$ if any in these sources. In NGC 3227, we found evidence for variation in $\rm{kT_{e}}$. We found no correlation of $\rm{kT_{e}}$, photon index ($Γ$), reflection fraction ($R$) and optical depth ($τ$) with flux, while, $τ$ is found to anti-correlate with $\rm{kT_{e}}$. These could be due to more than one physical process at work in the source that causes the change in $\rm{kT_{e}}$. Conclusive evidence for the variation in $\rm{kT_{e}}$ is not found in MR 2251$-$178 and NGC 5548. △ Less

Submitted 3 March, 2022; originally announced March 2022.

Comments: 18 pages, 18 figures, Accepted for publication in Astronomy & Astrophysics (A&A)

Journal ref: A&A 662, A78 (2022)

Intra-night optical variability monitoring of $γ$-ray emitting blazars

Authors: K. Subbu Ulaganatha Pandian, A. Natarajan, C. S. Stalin, Ashwani Pandey, S. Muneer, B. Natarajan

Abstract: We present the results obtained from our campaign to characterize the intra-night-optical variability properties of blazars detected by the {\it Fermi} Large Area Telescope. This involves R-band monitoring observations of a sample of 18 blazars, that includes five flat-spectrum radio quasars (FSRQs) and thirteen BL Lac objects (BL Lacs) covering the redshift range z = 0.085$-$1.184. Our observatio… ▽ More We present the results obtained from our campaign to characterize the intra-night-optical variability properties of blazars detected by the {\it Fermi} Large Area Telescope. This involves R-band monitoring observations of a sample of 18 blazars, that includes five flat-spectrum radio quasars (FSRQs) and thirteen BL Lac objects (BL Lacs) covering the redshift range z = 0.085$-$1.184. Our observations, carried out using the 1.3 m J.C. Bhattacharya Telescope cover a total of 40 nights ($\sim$200 hrs) between the period 2016 December and 2020 March. We characterized variability using the power enhanced $F-$test. We found a duty cycle (DC) of the variability of about 11\% for FSRQs and 12\% for BL Lacs. Dividing the sample into different sub-classes based on the position of the synchrotron peak in their broadband spectral energy distribution (SED), we found DC of $\sim$16\%, $\sim$10\% and $\sim$7\% for low-synchrotron peaked (LSP), intermediate synchrotron peaked (ISP) and high synchrotron peaked (HSP) blazars. Such high DC of variability in LSP blazars could be understood in the context of the R-band tracing the falling part (contributed by high energy electrons) of the synchrotron component of the broadband SED. Also, the R-band tracing the rising synchrotron part (produced by low energy electrons) in the case of ISP and HSP blazars, could cause lesser variability in them. Thus, the observed high DC of variability in LSP blazars relative to ISP and HSP blazars is in accordance with the leptonic model of emission from blazar jets. △ Less

Submitted 15 February, 2022; originally announced February 2022.

Comments: 18 pages, 9 figures, 4 tables, accepted for publication in The Journal of Astrophysics and Astronomy

Correlation between optical flux and polarization variations in Flat Spectrum Radio Quasars on diverse timescales

Authors: Ashwani Pandey, Bhoomika Rajput, C. S. Stalin

Abstract: Study of the polarization behaviour in blazars is a powerful tool to discern the role of magnetic field in the variable emission process in their relativistic jets. We present here results of our systematic investigation on the correlation between optical flux and polarization variations for eight flat-spectrum radio quasars on various timescales using data from the Steward Observatory that covers… ▽ More Study of the polarization behaviour in blazars is a powerful tool to discern the role of magnetic field in the variable emission process in their relativistic jets. We present here results of our systematic investigation on the correlation between optical flux and polarization variations for eight flat-spectrum radio quasars on various timescales using data from the Steward Observatory that covers a period of $\sim$10 years. On long time scales ($\sim$several months), from a total of 79 observing cycles, in 34 observing cycles, we found a significant positive correlation between optical flux and optical polarization degree (PD), negative correlation in 3 cycles and no correlation in 42 cycles. On short time scales ($\sim$few days), in 47 out of a total of 55 epochs, we found a positive correlation between optical flux and PD, while, on the remaining 8 epochs, an anti-correlation was detected between the two quantities. Moreover, we noticed a significant positive correlation between optical and $γ-$ray fluxes on 14 epochs and a negative correlation between the two on one epoch. While the observed optical flux changes well fit in the shock-in-jet model, the observed changes in PD are not explainable by changes in the power-law spectral index of the relativistic electrons in the jet. Instead, the observed varied correlations between optical flux and PD could be due to multi-zone emission regions or the enhanced flux coinciding with the emergence of a new emission knot with its magnetic field either aligned or misaligned with the large scale magnetic field. △ Less

Submitted 16 November, 2021; originally announced November 2021.

Comments: 27 pages, 5 figures, 7 tables, accepted for publication in MNRAS main journal

On the Origin of Gamma-ray Flares from Bright Fermi Blazars

Authors: Vaidehi S. Paliya, Markus Boettcher, Mark Gurwell, C. S. Stalin

Abstract: The origin of gamma-ray flares observed from blazars is one of the major mysteries in jet physics. We have attempted to address this problem following a novel spectral energy distribution (SED) fitting technique that explored the flaring patterns identified in the broadband SEDs of two gamma-ray bright blazars 3C 279 (z=0.54) and 3C 454.3 (z=0.86), using near-simultaneous radio-to-gamma-ray observ… ▽ More The origin of gamma-ray flares observed from blazars is one of the major mysteries in jet physics. We have attempted to address this problem following a novel spectral energy distribution (SED) fitting technique that explored the flaring patterns identified in the broadband SEDs of two gamma-ray bright blazars 3C 279 (z=0.54) and 3C 454.3 (z=0.86), using near-simultaneous radio-to-gamma-ray observations. For both sources, the gamma-ray flux strongly correlates with the separation of the SED peaks and the Compton dominance. We propose that spectral hardening of the radiating electron population and/or enhancement of the Doppler factor can naturally explain these observations. In both cases, magnetic reconnection may play a pivotal role in powering the luminous gamma-ray flares. △ Less

Submitted 8 November, 2021; originally announced November 2021.

Comments: ApJS, in press

Changing look AGN Mrk 590: Broad line region and black hole mass from photometric reverberation mapping

Authors: Amit Kumar Mandal, Malte Schramm, Suvendu Rakshit, C. S. Stalin, Bovornpratch Vijarnwannaluk, Wiphu Rujopakarn, Saran Poshyachinda, Vladimir V. Kouprianov, Joshua B. Haislip, Daniel E. Reichart, Ram Sagar, Blesson Mathew

Abstract: We present the results of photometric reverberation mapping observations on the changing look active galactic nucleus Mrk 590 at z = 0.026. The observations were carried out from July to December, 2018 using broad band B, R and narrow band Hα and S II filters. B-band traces the continuum emission from the accretion disk, R-band encompasses both the continuum emission from the accretion disk and th… ▽ More We present the results of photometric reverberation mapping observations on the changing look active galactic nucleus Mrk 590 at z = 0.026. The observations were carried out from July to December, 2018 using broad band B, R and narrow band Hα and S II filters. B-band traces the continuum emission from the accretion disk, R-band encompasses both the continuum emission from the accretion disk and the redshifted Hα line from the broad line region (BLR), S II band contains the redshifted Hα emission and the Hα band traces the continuum emission underneath the S II band. All the light curves showed strong variation with a fractional root-mean-square variation of 0.132 (+/-) 0.001 in B-band and 0.321 (+/-) 0.001 in Hα line. From cross-correlation function analysis, we obtained a delayed response of Hα line emission to the optical B-band continuum emission of 21.44(+1.49/-2.11) days in the rest-frame of the source, corresponding to a linear size of the BLR of 0.018 pc. This is consistent with previous estimates using Hβ. By combining the BLR size with the Hα line full width at half maximum of 6478 (+/-) 240 km/s measured from a single-epoch spectrum obtained with the Subaru telescope, we derived a black hole mass of 1.96 (+0.15 / -0.21) X 10^8 Msun. △ Less

Submitted 5 October, 2021; originally announced October 2021.

Comments: 16 pages, 12 figures, 4 tables, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS), 5 October 2021

High-redshift Narrow-line Seyfert 1 Galaxies: A Candidate Sample

Authors: Suvendu Rakshit, C. S. Stalin, Jari Kotilainen, Jaejin Shin

Abstract: The study of narrow-line Seyfert 1 galaxies (NLS1s) is now mostly limited to low redshift ($z<0.8$) because their definition requires the presence of the H$β$ emission line, which is redshifted out of the spectral coverage of major ground-based spectroscopic surveys at $z>0.8$. We studied the correlation between the properties of H$β$ and Mg II lines of a large sample of SDSS DR14 quasars to find… ▽ More The study of narrow-line Seyfert 1 galaxies (NLS1s) is now mostly limited to low redshift ($z<0.8$) because their definition requires the presence of the H$β$ emission line, which is redshifted out of the spectral coverage of major ground-based spectroscopic surveys at $z>0.8$. We studied the correlation between the properties of H$β$ and Mg II lines of a large sample of SDSS DR14 quasars to find high-$z$ NLS1 candidates. Based on the strong correlation of $\mathrm{FWHM(MgII)=(0.880\pm 0.005) \times FWHM(Hβ)+ (0.438\pm0.018)}$, we present a sample of high-$z$ NLS1 candidates having FWHM of Mg II $<$ 2000 km s$^{-1}$. The high-$z$ sample contains 2684 NLS1s with redshift $z=0.8-2.5$ with a median logarithmic bolometric luminosity of $46.16\pm0.42$ erg s$^{-1}$, logarithmic black hole mass of $8.01\pm0.35 M_{\odot}$, and logarithmic Eddington ratio of $0.02\pm0.27$. The fraction of radio-detected high-$z$ NLS1s is similar to that of the low-$z$ NLS1s and SDSS DR14 quasars at a similar redshift range, and their radio luminosity is found to be strongly correlated with their black hole mass. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: Published in ApJS on March 2021. Catalog is available on Zenodo (https://doi.org/10.5281/zenodo.4405039)

Journal ref: 2021ApJS..253...28R

The Sharpest Ultraviolet view of the star formation in an extreme environment of the nearest Jellyfish Galaxy IC 3418

Authors: Ananda Hota, Ashish Devaraj, Ananta C. Pradhan, C S Stalin, Koshy George, Abhisek Mohapatra, Soo-Chang Rey, Youichi Ohyama, Sravani Vaddi, Renuka Pechetti, Ramya Sethuram, Jessy Jose, Jayashree Roy, Chiranjib Konar

Abstract: We present the far ultraviolet (FUV) imaging of the nearest Jellyfish or Fireball galaxy IC3418/VCC 1217, in the Virgo cluster of galaxies, using Ultraviolet Imaging Telescope (UVIT) onboard the ASTROSAT satellite. The young star formation observed here in the 17 kpc long turbulent wake of IC3418, due to ram pressure stripping of cold gas surrounded by hot intra-cluster medium, is a unique laborat… ▽ More We present the far ultraviolet (FUV) imaging of the nearest Jellyfish or Fireball galaxy IC3418/VCC 1217, in the Virgo cluster of galaxies, using Ultraviolet Imaging Telescope (UVIT) onboard the ASTROSAT satellite. The young star formation observed here in the 17 kpc long turbulent wake of IC3418, due to ram pressure stripping of cold gas surrounded by hot intra-cluster medium, is a unique laboratory that is unavailable in the Milkyway. We have tried to resolve star forming clumps, seen compact to GALEX UV images, using better resolution available with the UVIT and incorporated UV-optical images from Hubble Space Telescope archive. For the first time, we resolve the compact star forming clumps (fireballs) into sub-clumps and subsequently into a possibly dozen isolated stars. We speculate that many of them could be blue supergiant stars which are cousins of SDSS J122952.66+112227.8, the farthest star (~17 Mpc) we had found earlier surrounding one of these compact clumps. We found evidence of star formation rate (4 - 7.4 x 10^-4 M_sun per yr ) in these fireballs, estimated from UVIT flux densities, to be increasing with the distance from the parent galaxy. We propose a new dynamical model in which the stripped gas may be developing vortex street where the vortices grow to compact star forming clumps due to self-gravity. Gravity winning over turbulent force with time or length along the trail can explain the puzzling trend of higher star formation rate and bluer/younger stars observed in fireballs farther away from the parent galaxy. △ Less

Submitted 29 April, 2021; originally announced April 2021.

Comments: 12 pages, 9 figures. Accepted for publication in the Special Issue of the Journal of Astrophysics and Astronomy on ASTROSAT

Correlation between optical and $γ$-ray flux variations in BL Lacs

Authors: Bhoomika Rajput, Zahir Shah, C. S. Stalin, S. Sahayanathan, Suvendu Rakshit

Abstract: We report here results on the analysis of correlated flux variations between the optical and GeV $γ$-ray bands in three bright BL Lac objects, namely AO\, 0235+164, OJ 287 and PKS 2155$-$304. This was based on the analysis of about 10 years of data from the {\it Fermi} Gamma-ray Space Telescope covering the period between 08 August 2008 to 08 August 2018 along with optical data covering the same p… ▽ More We report here results on the analysis of correlated flux variations between the optical and GeV $γ$-ray bands in three bright BL Lac objects, namely AO\, 0235+164, OJ 287 and PKS 2155$-$304. This was based on the analysis of about 10 years of data from the {\it Fermi} Gamma-ray Space Telescope covering the period between 08 August 2008 to 08 August 2018 along with optical data covering the same period. For all the sources, during the flares analysed in this work, the optical and $γ$-ray flux variations are found to be closely correlated. From broad band spectral energy distribution modelling of different epochs in these sources using the one zone leptonic emission model, we found that the optical-UV emission is dominated by synchrotron emission from the jet. The $γ$-ray emission in the low synchrotron peaked sources AO\, 0235+164 and OJ 287 are found to be well fit with external Compton (EC) component, while, the $γ$-ray emission in the high synchrotron peaked source PKS 2155$-$304 is well fit with synchrotron self Compton component. Further we note that the $γ$-ray emission during the high flux state of AO 0235+164 (epochs A and B) requires seed photons from both the dusty torus and broad line region, while the $γ$-ray emission in OJ 287 and during epochs C and D of AO\,0235+164 can be modelled by EC scattering of infra-red photons from the torus. △ Less

Submitted 8 April, 2021; v1 submitted 1 April, 2021; originally announced April 2021.

Comments: 17 pages, 26 figures, Accepted for publication in MNRAS (In the previous edition, a minor correction was made.)

Changing Accretion Geometry of Seyfert 1 Mrk 335 with {\it NuSTAR}: A Comparative Study

Authors: Santanu Mondal, C. S. Stalin

Abstract: We present a detailed spectral study of the narrow-line Seyfert~1 galaxy, Markarian~335, using eight epoch observations made between 2013 and 2020 with the Nuclear Spectroscopic Telescope Array. The source was variable during this period both in spectral flux and flow geometry. We estimated the height of the Compton cloud from the model fitted parameters for the whole observation period. This allo… ▽ More We present a detailed spectral study of the narrow-line Seyfert~1 galaxy, Markarian~335, using eight epoch observations made between 2013 and 2020 with the Nuclear Spectroscopic Telescope Array. The source was variable during this period both in spectral flux and flow geometry. We estimated the height of the Compton cloud from the model fitted parameters for the whole observation period. This allowed us to investigate the underlying physical processes that drive the variability in X-rays. Our model fitted mass varies in a narrow range, between $(2.44\pm0.45-3.04\pm0.56)\times10^{7}M_\odot$, however, given the large error bars, it is consistent with being constant and is in agreement with that known from optical reverberation mapping observations. The disk mass accretion rate reached a maximum of 10\% of the Eddington rate during June 2013. Our study sheds light on mass outflows from the system and also compares different aspects of accretion with X-ray binaries. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Comments: 15 pages, 4 tables, 6 figures

Journal ref: Galaxies 2021, 9(2), 21

TXS 1206+549: a new $γ$-ray detected narrow-line Seyfert 1 galaxy at redshift 1.34?

Authors: Suvendu Rakshit, Malte Schramm, C. S. Stalin, I. Tanaka, Vaidehi S. Paliya, Indrani Pal, Jari Kotilainen, Jaejin Shin

Abstract: Radio and $γ$-ray loud narrow-line Seyfert 1 galaxies (NLS1s) are unique objects to study the formation and evolution of relativistic jets, as they are believed to have high accretion rates and powered by low mass black holes contrary to that known for blazars. However, only about a dozen $γ$-ray detected NLS1s ($γ$-NLS1s) are known to date and all of them are at $z\le1$. Here, we report the ident… ▽ More Radio and $γ$-ray loud narrow-line Seyfert 1 galaxies (NLS1s) are unique objects to study the formation and evolution of relativistic jets, as they are believed to have high accretion rates and powered by low mass black holes contrary to that known for blazars. However, only about a dozen $γ$-ray detected NLS1s ($γ$-NLS1s) are known to date and all of them are at $z\le1$. Here, we report the identification of a new $γ$-ray emitting NLS1 TXS 1206+549 at $z=1.344$. A near-infrared spectrum taken with the Subaru telescope showed H$β$ emission line with FWHM of $1194\pm77$ km s$^{-1}$ and weak [O III] emission line but no optical Fe II lines, due to the limited wavelength coverage and poor signal-to-noise ratio. However, UV Fe II lines are present in the SDSS optical spectrum. The source is very radio-loud, unresolved, and has a flat radio spectrum. The broadband SED of the source has the typical two hump structure shown by blazars and other $γ$-NLS1s. The source exhibits strong variability at all wavelengths such as the optical, infrared, and $γ$-ray bands. All these observed characteristics show that TXS 1206+549 is the most distant $γ$-NLS1 known to date. △ Less

Submitted 31 March, 2021; v1 submitted 27 March, 2021; originally announced March 2021.

Comments: Accepted in MNRAS Letter

Broadband study of OQ 334 during its flaring state

Authors: Raj Prince, Rukaiya Khatoon, C. S. Stalin

Abstract: The blazar OQ 334 displayed a γ-ray flare in 2018, after being in the long quiescent γ-ray state since 2008. Subsequent to the flare, the source was in a higher γ-ray flux state and again flared in 2020. We present here the first spectral and timing analysis of the source at its various flaring states. During the higher γ-ray state, we found four major peaks identified as P1, P2, P3, and P4. From… ▽ More The blazar OQ 334 displayed a γ-ray flare in 2018, after being in the long quiescent γ-ray state since 2008. Subsequent to the flare, the source was in a higher γ-ray flux state and again flared in 2020. We present here the first spectral and timing analysis of the source at its various flaring states. During the higher γ-ray state, we found four major peaks identified as P1, P2, P3, and P4. From timing analysis, we found the rise and decay time of the order of hours with the fastest variability time of 9.01+/-0.78 hr. We found the highest γ-ray photon of 77 GeV during P4, which suggests the location of the γ-ray emitting region at the outer edge of the broad-line region or the inner edge of the torus. The γ-ray spectral analysis of the source indicates that during P4, the γ-ray spectrum clearly deviates from the power-law behavior. From cross-correlation analysis of the γ-ray and radio lightcurves, we found that the two emission regions are separated by about 11 pc. Our broadband spectral energy distribution modeling of the source during quiescent and active phases indicates that more electron and proton power are required to change the source from low flux to high flux state. The Anderson-Darling test and histogram fitting results suggest that the three days binned γ-ray fluxes follow a lognormal distribution. △ Less

Submitted 6 February, 2021; originally announced February 2021.

Comments: 15 pages, 12 figures, 5 tables, Accepted for publication in MNRAS

Multiwavelength monitoring of NGC 1275 over a decade: Evidence of a shift in synchrotron peak frequency and long-term multi-band flux increase

Authors: Sanna Gulati, Debbijoy Bhattacharya, Subir Bhattacharyya, Nilay Bhatt, C. S. Stalin, V. K. Agrawal

Abstract: We carried out a detailed study of the temporal and broadband spectral behaviour of one of the brightest misaligned active galaxies in gamma-rays, NGC 1275 utilising 11 years of Fermi, and available Swift and AstroSat observations. Based on the cumulative flux distribution of the gamma-ray lightcurve, we identified four distinct activity states and noticed an increase in the baseline flux during t… ▽ More We carried out a detailed study of the temporal and broadband spectral behaviour of one of the brightest misaligned active galaxies in gamma-rays, NGC 1275 utilising 11 years of Fermi, and available Swift and AstroSat observations. Based on the cumulative flux distribution of the gamma-ray lightcurve, we identified four distinct activity states and noticed an increase in the baseline flux during the first three states. Similar nature of the increase in the average flux was also noticed in X-ray and UV bands. A large flaring activity in gamma-rays was noticed in the fourth state. The source was observed twice by AstroSat for shorter intervals (~days) during the longer observing periods (~years) state 3 and 4. During AstroSat observing periods, the source gamma-ray flux was higher than the average flux observed during longer duration states. The increase in the average baseline flux from state 1 to state 3 can be explained considering a corresponding increase of jet particle normalisation. The inverse Comptonisation of synchrotron photons explained the average X-ray and gamma-ray emission by jet electrons during the first three longer duration states. However, during the shorter duration AstroSat observing periods, a shift of the synchrotron peak frequency was noticed, and the synchrotron emission of jet electrons well explained the observed X-ray flux. △ Less

Submitted 10 July, 2021; v1 submitted 27 January, 2021; originally announced January 2021.

Comments: Corrected AstroSat-SXT and LAXPC instrument description in Section 2.2.2 and updated references to match the published version. Results unchanged

Journal ref: MNRAS, Volume 503, Issue 1, May 2021, Pages 446-457

Discovery of soft and hard X-ray time lags in low-mass AGNs

Authors: Labani Mallick, Daniel R. Wilkins, William N. Alston, Alex Markowitz, Barbara De Marco, Michael L. Parker, Anne M. Lohfink, C. S. Stalin

Abstract: The scaling relations between the black hole (BH) mass and soft lag properties for both active galactic nuclei (AGNs) and BH X-ray binaries (BHXRBs) suggest the same underlying physical mechanism at work in accreting BH systems spanning a broad range of mass. However, the low-mass end of AGNs has never been explored in detail. In this work, we extend the existing scaling relations to lower-mass AG… ▽ More The scaling relations between the black hole (BH) mass and soft lag properties for both active galactic nuclei (AGNs) and BH X-ray binaries (BHXRBs) suggest the same underlying physical mechanism at work in accreting BH systems spanning a broad range of mass. However, the low-mass end of AGNs has never been explored in detail. In this work, we extend the existing scaling relations to lower-mass AGNs, which serve as anchors between the normal-mass AGNs and BHXRBs. For this purpose, we construct a sample of low-mass AGNs ($M_{\rm BH}<3\times 10^{6} M_{\rm \odot}$) from the XMM-Newton archive and measure frequency-resolved time delays between the soft (0.3-1 keV) and hard (1-4 keV) X-ray emissions. We report that the soft band lags behind the hard band emission at high frequencies $\sim[1.3-2.6]\times 10^{-3}$ Hz, which is interpreted as a sign of reverberation from the inner accretion disc in response to the direct coronal emission. At low frequencies ($\sim[3-8]\times 10^{-4}$ Hz), the hard band lags behind the soft band variations, which we explain in the context of the inward propagation of luminosity fluctuations through the corona. Assuming a lamppost geometry for the corona, we find that the X-ray source of the sample extends at an average height and radius of $\sim 10r_{\rm g}$ and $\sim 6r_{\rm g}$, respectively. Our results confirm that the scaling relations between the BH mass and soft lag amplitude/frequency derived for higher-mass AGNs can safely extrapolate to lower-mass AGNs, and the accretion process is indeed independent of the BH mass. △ Less

Submitted 6 April, 2021; v1 submitted 23 January, 2021; originally announced January 2021.

Comments: 11 pages, 5 figures, 4 tables, Published in MNRAS

Journal ref: MNRAS 503, 3775-3783 (2021)

Curvit: An open-source Python package to generate light curves from UVIT data

Authors: P. Joseph, C. S. Stalin, S. N. Tandon, S. K. Ghosh

Abstract: Curvit is an open-source Python package that facilitates the creation of light curves from the data collected by the Ultra-Violet Imaging Telescope (UVIT) onboard AstroSat, India's first multi-wavelength astronomical satellite. The input to Curvit is the calibrated events list generated by the UVIT-Payload Operation Center (UVIT-POC) and made available to the principal investigators through the In… ▽ More Curvit is an open-source Python package that facilitates the creation of light curves from the data collected by the Ultra-Violet Imaging Telescope (UVIT) onboard AstroSat, India's first multi-wavelength astronomical satellite. The input to Curvit is the calibrated events list generated by the UVIT-Payload Operation Center (UVIT-POC) and made available to the principal investigators through the Indian Space Science Data Center. The features of Curvit include (i) automatically detecting sources and generating light curves for all the detected sources and (ii) custom generation of light curve for any particular source of interest. We present here the capabilities of Curvit and demonstrate its usability on the UVIT observations of the intermediate polar FO Aqr as an example. Curvit is publicly available on GitHub at https://github.com/prajwel/curvit. △ Less

Submitted 16 January, 2021; originally announced January 2021.

Comments: Accepted for publication in Journal of Astrophysics & Astronomy, 11 pages, 5 figures

Estimation of the size and structure of the broad line region using Bayesian approach

Authors: Amit Kumar Mandal, Suvendu Rakshit, C. S. Stalin, R. G. Petrov, Blesson Mathew, Ram Sagar

Abstract: Understanding the geometry and kinematics of the broad line region (BLR) of active galactic nuclei (AGN) is important to estimate black hole masses in AGN and study the accretion process. The technique of reverberation mapping (RM) has provided estimates of BLR size for more than 100 AGN now, however, the structure of the BLR has been studied for only a handful number of objects. Towards this, we… ▽ More Understanding the geometry and kinematics of the broad line region (BLR) of active galactic nuclei (AGN) is important to estimate black hole masses in AGN and study the accretion process. The technique of reverberation mapping (RM) has provided estimates of BLR size for more than 100 AGN now, however, the structure of the BLR has been studied for only a handful number of objects. Towards this, we investigated the geometry of the BLR for a large sample of 57 AGN using archival RM data. We performed systematic modeling of the continuum and emission line light curves using a Markov Chain Monte Carlo method based on Bayesian statistics implemented in PBMAP (Parallel Bayesian code for reverberation-MAPping data) code to constrain BLR geometrical parameters and recover velocity integrated transfer function. We found that the recovered transfer functions have various shapes such as single-peaked, double-peaked and top-hat suggesting that AGN have very different BLR geometries. Our model lags are in general consistent with that estimated using the conventional cross-correlation methods. The BLR sizes obtained from our modeling approach is related to the luminosity with a slope of 0.583 (+/-) 0.026 and 0.471 (+/-) 0.084 based on Hβ and Hα lines, respectively. We found a non-linear response of emission line fluxes to the ionizing optical continuum for 93\% objects. The estimated virial factors for the AGN studied in this work range from 0.79 to 4.94 having a mean at 1.78 (+/-) 1.77 consistent with the values found in the literature. △ Less

Submitted 4 January, 2021; originally announced January 2021.

Comments: 20 pages, 16 figures, 5 tables, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS), 30 December 2020

Performance of the UVIT Level-2 Pipeline

Authors: S. K. Ghosh, S. N. Tandon, P. Joseph, A. Devaraj, D. S. Shelat, C. S. Stalin

Abstract: Performance of the Level-2 pipeline, which translates the UVIT data created by the ISRO's ground segment processing systems (Level-1) into astronomer ready scientific data products, is described. This pipeline has evolved significantly from experiences during the in orbit mission. With time, the detector modules of UVIT developed certain defects which led to occasional corruption of imaging and ti… ▽ More Performance of the Level-2 pipeline, which translates the UVIT data created by the ISRO's ground segment processing systems (Level-1) into astronomer ready scientific data products, is described. This pipeline has evolved significantly from experiences during the in orbit mission. With time, the detector modules of UVIT developed certain defects which led to occasional corruption of imaging and timing data. This article will describe the improvements and mitigation plans incorporated in the pipeline and report its efficacy and quantify the performance. △ Less

Submitted 25 December, 2020; originally announced December 2020.

Comments: Accepted for publication in Journal of Astrophysics & Astronomy, 6 pages, 4 figures

In orbit performance of UVIT over the 5 years

Authors: S. K. Ghosh, P. Joseph, A. Kumar, J. Postma, C. S. Stalin, A. Subramaniam, S. N. Tandon, I. V. Barve, A. Devaraj, K. George, V. Girish, J. B. Hutchings, P. U. Kamath, S. Kathiravan, J. P. Lancelot, D. Leahy, P. K. Mahesh, R. Mohan, S. Nagabhushana, A. K. Pati, N. Kameswara Rao, K. Sankarasubramanian, P. Sreekumar, S. Sriram

Abstract: Over the last 5 years, UVIT has completed observations of more than 500 proposals with ~ 800 unique pointings. In addition, regular planned monitoring observations have been made and from their analysis various key parameters related to in orbit performance of UVIT have been quantified. The sensitivities of the UV channels have remained steady indicating no effect of potential molecular contaminat… ▽ More Over the last 5 years, UVIT has completed observations of more than 500 proposals with ~ 800 unique pointings. In addition, regular planned monitoring observations have been made and from their analysis various key parameters related to in orbit performance of UVIT have been quantified. The sensitivities of the UV channels have remained steady indicating no effect of potential molecular contamination confirming the adequacy of all the protocols implemented for avoiding contamination. The quality of the PSF through the years confirms adequacy of thermal control measures. The early calibrations obtained during the Performance Verification (PV) phase have been further revised for more subtle effects. These include flat fields and detector distortions with greater precision. The operations of UVIT have also evolved through in orbit experience, e.g. tweaking of operational sequencing, protocol for recovery from bright object detection (BOD) shutdowns, parameters for BOD thresholds, etc. Finally, some effects of charged particle hits on electronics led to opimised strategy for regular resetting. The Near-UV channel was lost in one of such operations. All the above in-orbit experiences are presented here. △ Less

Submitted 25 December, 2020; originally announced December 2020.

Comments: Accepted for publication in Journal of Astrophysics & Astronomy, 11 pages, 10 figures, 1 table

Dust Reverberation Mapping of Z229-15

Authors: Amit Kumar Mandal, Suvendu Rakshit, C. S. Stalin, Dominika Wylezalek, Markus Kissler Patig, Ram Sagar, Blesson Mathew, S. Muneer, Indrani Pal

Abstract: We report results of the dust reverberation mapping (DRM) on the Seyfert 1 galaxy Z229-15 at z = 0.0273. Quasi-simultaneous photometric observations for a total of 48 epochs were acquired during the period 2017 July to 2018 December in B, V, J, H and Ks bands. The calculated spectral index (α) between B and V bands for each epoch was used to correct for the accretion disk (AD) component present in… ▽ More We report results of the dust reverberation mapping (DRM) on the Seyfert 1 galaxy Z229-15 at z = 0.0273. Quasi-simultaneous photometric observations for a total of 48 epochs were acquired during the period 2017 July to 2018 December in B, V, J, H and Ks bands. The calculated spectral index (α) between B and V bands for each epoch was used to correct for the accretion disk (AD) component present in the infrared light curves. The observed α ranges between -0.99 and 1.03. Using cross correlation function analysis we found significant time delays between the optical V and the AD corrected J, H and Ks light curves. The lags in the rest frame of the source are 12.52 (+10.00/-9.55) days (between V and J), 15.63 (+5.05/-5.11) days (between V and H) and 20.36 (+5.82/-5.68) days (between V and Ks). Given the large error bars, these lags are consistent with each other. However, considering the lag between V and Ks bands to represent the inner edge of the dust torus, the torus in Z229-15 lies at a distance of 0.017 pc from the central ionizing continuum. This is smaller than that expected from the radius luminosity (R-L) relationship known from DRM. Using a constant α = 0.1 to account for the AD component, as is normally done in DRM, the deduced radius (0.025 pc) lies close to the expected R-L relation. However, usage of constant α in DRM is disfavoured as the α of the ionizing continuum changes with the flux of the source. △ Less

Submitted 9 December, 2020; originally announced December 2020.

Comments: 12 pages, 13 figures, 3 tables, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS), 08 December 2020

The First Gamma-ray Emitting BL Lacertae Object at the Cosmic Dawn

Authors: Vaidehi S. Paliya, A. Domínguez, C. Cabello, N. Cardiel, J. Gallego, B. Siana, M. Ajello, D. Hartmann, A. Gil de Paz, C. S. Stalin

Abstract: One of the major challenges in studying the cosmic evolution of relativistic jets is the identification of the high-redshift ($z>3$) BL Lacertae objects, a class of jetted active galactic nuclei characterized by their quasi-featureless optical spectra. Here we report the identification of the first $γ$-ray emitting BL Lac object, 4FGL~J1219.0+3653 (J1219), beyond $z=3$, i.e., within the first two… ▽ More One of the major challenges in studying the cosmic evolution of relativistic jets is the identification of the high-redshift ($z>3$) BL Lacertae objects, a class of jetted active galactic nuclei characterized by their quasi-featureless optical spectra. Here we report the identification of the first $γ$-ray emitting BL Lac object, 4FGL~J1219.0+3653 (J1219), beyond $z=3$, i.e., within the first two billion years of the age of the Universe. The optical and near-infrared spectra of J1219 taken from 10.4 m Gran Telescopio Canarias exhibit no emission lines down to an equivalent width of $\sim$3.5 A supporting its BL Lac nature. The detection of a strong Lyman-$α$ break at $\sim$5570 A, on the other hand, confirms that J2119 is indeed a high-redshift ($z\sim3.59$) quasar. Based on the prediction of a recent BL Lac evolution model, J1219 is one of the only two such objects expected to be present within the comoving volume at $z=3.5$. Future identifications of more $z>3$ $γ$-ray emitting BL Lac sources, therefore, will be crucial to verify the theories of their cosmic evolution. △ Less

Submitted 24 October, 2020; originally announced October 2020.

Comments: Accepted for publication in The Astrophysical Journal Letters

Correlation between optical and $γ$-ray flux variations in bright flat spectrum radio quasars

Authors: Bhoomika Rajput, C. S. Stalin, S. Sahayanathan

Abstract: Blazars are known to show flux variations over a range of energies from low energy radio to high energy gamma-rays. Cross-correlation analysis of the optical and $γ$-ray light curves in blazars shows that flux variations are generally correlated in both bands, however, there are exceptions. We explored this optical-GeV connection in four flat spectrum radio quasars (FSRQs) by a systematic investig… ▽ More Blazars are known to show flux variations over a range of energies from low energy radio to high energy gamma-rays. Cross-correlation analysis of the optical and $γ$-ray light curves in blazars shows that flux variations are generally correlated in both bands, however, there are exceptions. We explored this optical-GeV connection in four flat spectrum radio quasars (FSRQs) by a systematic investigation of their long term optical and $γ$-ray light curves. On analysis of the four sources, namely 3C 273, 3C 279, PKS 1510$-$089 and CTA 102 we noticed different behaviours between the optical and GeV flux variations. We found instances when (i) the optical and GeV flux variations are closely correlated (ii) there are optical flares without $γ$-ray counterparts and (iii) $γ$-ray flares without optical counterparts. To understand these diverse behaviours, we carried out broad band spectral energy distribution (SED) modelling of the sources at different epochs using a one-zone leptonic emission model. The optical-UV emission is found to be dominated by emission from the accretion disk in the sources PKS 1510$-$089, CTA 102 and 3C 273, while in 3C 279, the synchrotron radiation from the jet dominates the optical-UV emission. Our SED analysis indicates that (i) correlated optical and $γ$-ray flux variations are caused by changes in the bulk Lorentz factor ($Γ$), (ii) $γ$-ray flares without optical counterparts are due to increase in $Γ$ and/or the electron energy density and (iii) an optical flare without $γ$-ray counterpart is due to increase in the magnetic field strength. △ Less

Submitted 4 September, 2020; originally announced September 2020.

Comments: 22 pages, 32 figures, accepted to MNRAS