Showing 1–50 of 248 results for author: Czerny, B

Chemical abundances along the quasar main sequence

Authors: A. Floris, P. Marziani, S. Panda, M. Sniegowska, M. D'Onofrio, A. Deconto-Machado, A. Del Olmo, B. Czerny

Abstract: The 4D eigenvector 1 (E1) sequence has emerged as a powerful tool for organizing the observational and physical characteristics of type-1 active galactic nuclei (AGNs). In this study, we present a comprehensive analysis of the metallicity of the broad line region gas, incorporating both new data and previously published findings, to assess the presence of any trend along the sequence. We perform a… ▽ More The 4D eigenvector 1 (E1) sequence has emerged as a powerful tool for organizing the observational and physical characteristics of type-1 active galactic nuclei (AGNs). In this study, we present a comprehensive analysis of the metallicity of the broad line region gas, incorporating both new data and previously published findings, to assess the presence of any trend along the sequence. We perform a multi-component analysis on the strongest UV and optical emission lines, compute $\sim 10$ diagnostic ratios, and compare them with the prediction of CLOUDY photoionization simulations, identifying a photoionization solution closest to the data. Our investigation reveals a consistent pattern along the optical plane of the E1. We observe a systematic progression in metallicity, ranging from sub-solar values to metallicity levels several times higher than solar values. These findings underscore the role of metallicity as a fundamental correlate of the 4DE1/main sequence. Extreme values of metallicity, at least several tens solar, are confirmed in low-$z$ AGNs radiating at a high Eddington ratio, although the origin of the extreme enrichment remains open to debate. △ Less

Submitted 7 May, 2024; originally announced May 2024.

Unveiling Energy Pathways in AGN Accretion Flows with the Warm Corona Model for the Soft Excess

Authors: D. R. Ballantyne, V. Sudhakar, D. Fairfax, S. Bianchi, B. Czerny, A. De Rosa, B. De Marco, R. Middei, B. Palit, P. -O. Petrucci, A. Rozanska, F. Ursini

Abstract: The soft excess in active galactic nuclei (AGNs) may arise through a combination of relativistic reflection and the effects of a warm corona at the surface of the accretion disc. Detailed examination of the soft excess can therefore constrain models of the transport and dissipation of accretion energy. Here, we analyze 34 XMM-Newton observations from 14 Type I AGNs with the reXcor spectral model w… ▽ More The soft excess in active galactic nuclei (AGNs) may arise through a combination of relativistic reflection and the effects of a warm corona at the surface of the accretion disc. Detailed examination of the soft excess can therefore constrain models of the transport and dissipation of accretion energy. Here, we analyze 34 XMM-Newton observations from 14 Type I AGNs with the reXcor spectral model which self-consistently combines emission from a warm corona with relativistic reflection assuming a lamppost corona. The model divides accretion energy between the disc, the warm corona, and the lamppost. The XMM-Newton observations span a factor of 188 in Eddington ratio ($λ_{\mathrm{obs}}$) and 350 in black hole mass, and we find that a warm corona is a significant contributor to the soft excess for 13 of the 14 AGNs with a mean warm corona heating fraction of $0.51$. The reXcor fits reveal that the fraction of accretion energy dissipated in the lamppost is anti-correlated with $λ_{\mathrm{obs}}$. In contrast, the relationship between $λ_{\mathrm{obs}}$ and both the optical depth and heating fraction of the warm corona appears to transition from an anti-correlation to a correlation at $λ_{\mathrm{obs,t}} \approx 0.15$. Therefore, at least one other physical process in addition to the accretion rate is needed to explain the evolution of the warm corona. Overall, we find that a warm corona appears to be a crucial depository of accretion energy in AGNs across a broad range of $λ_{\mathrm{obs}}$ and black hole mass. △ Less

Submitted 3 April, 2024; originally announced April 2024.

Comments: 22 pages, 44 figures, accepted by MNRAS

Magnetically arrested disks in FR I radio galaxies

Authors: Han He, Bei You, Ning Jiang, Xinwu Cao, Jingfu Hu, Zhenfeng Sheng, Su Yao, Bozena Czerny

Abstract: A sample of 17 FR I radio galaxies constructed from the 3CR catalog, which is characterized by edge-darkened radio structures, is studied. The optical core luminosities derived from Hubble Space Telescope observation are used to estimate the Eddington ratios which are found to be below $10^{-3.4}$ for this sample. This is supported by the Baldwin-Phillips-Terlevich optical diagnostic diagrams deri… ▽ More A sample of 17 FR I radio galaxies constructed from the 3CR catalog, which is characterized by edge-darkened radio structures, is studied. The optical core luminosities derived from Hubble Space Telescope observation are used to estimate the Eddington ratios which are found to be below $10^{-3.4}$ for this sample. This is supported by the Baldwin-Phillips-Terlevich optical diagnostic diagrams derived with the spectroscopic observation of Telescopio Nazionale Galileo, suggesting that these sources are of low ionization nuclear Emission-line Regions (LINERs). It implies that the accretion in these FR I sources can be modeled as advection-dominated accretion flows (ADAFs). Given the low accretion rate, the predicted jet power with a fast-spinning black hole (BH) $a=0.95$ in the Blandford-Znajek mechanics is lower than the estimated one for almost all the sources in our sample. Such powerful jets indicate the presence of magnetically arrested disks (MAD) in the inner region of the ADAF, in the sense that the magnetic fields in the inner accretion zone are strong. Moreover, we show that, even in the MAD scenario, the BH spins in the sample are most likely moderate and/or fast with $a\gtrsim0.5$. △ Less

Submitted 22 March, 2024; originally announced March 2024.

Comments: 10 pages, 10 figures, 3 tables, Accepted for publication in MNRAS

Reevaluating LSST's Capability for Time Delay Measurements in Quasar Accretion Discs

Authors: F. Pozo Nuñez, B. Czerny, S. Panda, A. Kovacevic, W. Brandt, K. Horne

Abstract: The Legacy Survey of Space and Time (LSST) at the Vera C. Rubin Observatory is poised to observe thousands of quasars using the Deep Drilling Fields (DDF) across six broadband filters over a decade. Understanding quasar accretion disc (AD) time delays is pivotal for probing the physics of these distant objects. Pozo Nuñez et al. (2023) has recently demonstrated the feasibility of recovering AD tim… ▽ More The Legacy Survey of Space and Time (LSST) at the Vera C. Rubin Observatory is poised to observe thousands of quasars using the Deep Drilling Fields (DDF) across six broadband filters over a decade. Understanding quasar accretion disc (AD) time delays is pivotal for probing the physics of these distant objects. Pozo Nuñez et al. (2023) has recently demonstrated the feasibility of recovering AD time delays with accuracies ranging from 5\% to 20\%, depending on the quasar's redshift and time sampling intervals. Here we reassess the potential for measuring AD time delays under the current DDF observing cadence, which is placeholder until a final cadence is decided. We find that contrary to prior expectations, achieving reliable AD time delay measurements for quasars is significantly more challenging, if not unfeasible, due to the limitations imposed by the current observational strategies. △ Less

Submitted 14 February, 2024; originally announced February 2024.

Comments: Accepted as a research note in RNAAS; 3 pages, 1 figure

Journal ref: 2024 Res. Notes AAS 8 47

Lense-Thirring Precession after a Supermassive Black Hole Disrupts a Star

Authors: Dheeraj R. Pasham, Michal Zajacek, C. J. Nixon, Eric R. Coughlin, Marzena Sniegowska, Agnieszka Janiuk, Bozena Czerny, Thomas Wevers, Muryel Guolo, Yukta Ajay, Michael Loewenstein

Abstract: An accretion disk formed around a supermassive black hole (SMBH) after it disrupts a star is expected to be initially misaligned with respect to the black hole's equatorial plane. This misalignment induces relativistic torques (the Lense-Thirring effect) on the disk, causing the disk to precess at early times, while at late times the disk aligns with the black hole and precession terminates. Here,… ▽ More An accretion disk formed around a supermassive black hole (SMBH) after it disrupts a star is expected to be initially misaligned with respect to the black hole's equatorial plane. This misalignment induces relativistic torques (the Lense-Thirring effect) on the disk, causing the disk to precess at early times, while at late times the disk aligns with the black hole and precession terminates. Here, using high-cadence X-ray monitoring observations of a TDE, we report the discovery of strong, quasi-periodic X-ray flux and temperature modulations from a TDE. These X-ray modulations are separated by 17.0$^{+1.2}_{-2.4}$ days and persist for roughly 130 days during the early phase of the TDE. Lense-Thirring precession of the accretion flow can produce this X-ray variability, but other physical mechanisms, such as the radiation-pressure instability, cannot be ruled out. Assuming typical TDE parameters, i.e., a solar-like star with the resulting disk extending at-most to so-called circularization radius, and that the disk precesses as a rigid body, we constrain the disrupting black hole's dimensionless spin parameter to be 0.05<|a|<0.5. △ Less

Submitted 14 February, 2024; originally announced February 2024.

Comments: Accepted for publication in Nature

New theoretical Fe II templates for bright quasars

Authors: Ashwani Pandey, Mary Loli Martínez-Aldama, Bożena Czerny, Swayamtrupta Panda, Michal Zajaček

Abstract: We present a set of new theoretical Fe II templates for bright quasars covering a wavelength range of 1000-10000 Å, based on the recent atomic database available in the C23.00 version of the photoionization code CLOUDY. We compute a grid of models for a range of incident photon flux and gas density and for multiple microturbulence velocities. We analyze the ratios of Fe II emission over a variety… ▽ More We present a set of new theoretical Fe II templates for bright quasars covering a wavelength range of 1000-10000 Å, based on the recent atomic database available in the C23.00 version of the photoionization code CLOUDY. We compute a grid of models for a range of incident photon flux and gas density and for multiple microturbulence velocities. We analyze the ratios of Fe II emission over a variety of wavebands and compare them with observations. Our key results are: (1) Despite the use of the newest atomic data we still confirm the long-standing problem that the predicted Fe II UV/optical ratio is significantly larger than that observed in the AGN spectra. (2) The ratio is not significantly affected by the variations in the microturbulence and the metallicity. (3) The turbulence can create an additional apparent velocity shift of up to 1000 km/s in the spectra. (4) There is no single Fe II template that can fit the observational data covering UV to optical wavelength range. We shortly discuss the most likely effects responsible for the Fe II UV/optical mismatch problem: the assumption of the constant column density and the assumption of the isotropic emission implying equal contribution of the bright irradiated faces and the dark shielded faces of the clouds. △ Less

Submitted 31 January, 2024; originally announced January 2024.

Comments: 21 pages, 16 figures, 4 tables, submited to the Astrophysical Journal Supplement (ApJS) series

Covering factor of the dust-driven broad-line region clouds

Authors: Mohammad Hassan Naddaf, Bozena Czerny

Abstract: The origin of the broad-line region (BLR) clouds in active galactic nuclei is still under discussion. We develop a scenario in which the clouds in the outer, less ionized part of the BLR are launched by the radiation pressure acting on dust. Most of the outflow forms a failed wind, so we refer to it as failed radiatively accelerated dusty outflow (FRADO), but, for a certain parameter range, actual… ▽ More The origin of the broad-line region (BLR) clouds in active galactic nuclei is still under discussion. We develop a scenario in which the clouds in the outer, less ionized part of the BLR are launched by the radiation pressure acting on dust. Most of the outflow forms a failed wind, so we refer to it as failed radiatively accelerated dusty outflow (FRADO), but, for a certain parameter range, actual outflow also takes place. We aim to test the model predictions. In this paper, we present the calculation of the angular distribution of clouds and the net covering factor as this affects the fraction of radiation that can be intercepted and reprocessed in the form of the H-beta or MgII emission line. The results reveal that the covering factor is intricately linked to the mass, accretion rate, and metallicity of the clouds. Notably, as these parameters increase, so does the covering factor, shedding light on the dynamic interplay between the central engine and the surrounding material in AGNs. △ Less

Submitted 9 January, 2024; v1 submitted 18 December, 2023; originally announced December 2023.

Comments: 12 pages, 3 figures, 1 table, Accepted for publication in the Journal of Universe (MDPI)

Accretion processes onto black holes: theoretical problems, observational constraints

Authors: Bozena Czerny, Marzena Sniegowska, Agnieszka Janiuk, Bei You

Abstract: We shortly summarize the standard current knowledge on the structure of the accretion flow onto black holes in galactic binary systems and in active galactic nuclei. We stress the similarities and differences between the two types of systems, and we highlight the complementarity of the data caused by these differences. We highlight some new developments and list the unsolved problems. We shortly summarize the standard current knowledge on the structure of the accretion flow onto black holes in galactic binary systems and in active galactic nuclei. We stress the similarities and differences between the two types of systems, and we highlight the complementarity of the data caused by these differences. We highlight some new developments and list the unsolved problems. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: submitted to the Proceedings of the Polish Astronomical Society, comments welcome within the allowed 6 page limit

Spectral variability studies in Active Galactic Nuclei: Exploring continuum and emission line regions in the age of LSST and JWST

Authors: Swayamtrupta Panda, Paola Marziani, Bozena Czerny, Alberto Rodriguez-Ardila, Francisco Pozo Nunez

Abstract: The investigation of emission line regions within active galaxies (AGNs) has a rich and extensive history, now extending to the use of AGNs and quasars as "standardizable" cosmological indicators, shedding light on the evolution of our universe. As we enter the era of advanced observatories, such as the successful launch of JWST and the forthcoming Vera C. Rubin Observatory's Legacy Survey of Spac… ▽ More The investigation of emission line regions within active galaxies (AGNs) has a rich and extensive history, now extending to the use of AGNs and quasars as "standardizable" cosmological indicators, shedding light on the evolution of our universe. As we enter the era of advanced observatories, such as the successful launch of JWST and the forthcoming Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST), the landscape of AGN exploration across cosmic epochs is poised for exciting advancements. In this work, we delve into recent developments in AGN variability research, anticipating the substantial influx of data facilitated by LSST. The article highlights recent strides made by the AGN Polish Consortium in their contributions to LSST. The piece emphasizes the role of quasars in cosmology, dissecting the intricacies of their calibration as standard candles. The primary focus centers on the relationship between the broad-line region size and luminosity, showcasing recent breakthroughs that enhance our comprehension of this correlation. These breakthroughs encompass a range of perspectives, including spectroscopic analyses, photoionization modeling, and collaborative investigations with other cosmological tools. The study further touches on select studies, underlining how the synergy of theoretical insights and advancements in observational capabilities has yielded deeper insights into these captivating cosmic entities. △ Less

Submitted 22 November, 2023; originally announced November 2023.

Comments: 34 pages, 5 figures, accepted for publication as a review in MDPI/Universe

Journal ref: Universe 2023, 9(12), 492

Investigating the origin of optical flares from the TeV blazar S4 0954+65

Authors: Ashwani Pandey, Rumen Bachev, Bozena Czerny, Paul J. Wiita, Alok C. Gupta, Anton Strigachev, Adam Popowicz

Abstract: Aims. We aim to investigate the extreme variability properties of the TeV blazar S4 0954+65 using optical photometric and polarisation observations carried out between 2017 and 2023 using three ground-based telescopes. Methods. We examined an extensive dataset comprised of 138 intraday (observing duration shorter than a day) light curves (LCs) of S4 0954+65 for flux, spectral, and polarisation v… ▽ More Aims. We aim to investigate the extreme variability properties of the TeV blazar S4 0954+65 using optical photometric and polarisation observations carried out between 2017 and 2023 using three ground-based telescopes. Methods. We examined an extensive dataset comprised of 138 intraday (observing duration shorter than a day) light curves (LCs) of S4 0954+65 for flux, spectral, and polarisation variations on diverse timescales. For the variable LCs, we computed the minimum variability timescales. We investigated flux-flux correlations and colour variations to look for spectral variations on long (several weeks to years) timescales. Additionally, we looked for connections between optical R-band flux and polarisation degree. Results. We found significant variations in 59 out of 138 intraday LCs. We detected a maximum change of 0.58$\pm$0.11 in V-band magnitude within $\sim$2.64 h and a corresponding minimum variability timescale of 18.21$\pm$4.87 mins on 2017 March 25. During the course of our observing campaign, the source brightness changed by $\sim$4 magnitudes in V and R bands; however, we did not find any strong spectral variations. The slope of the relative spectral energy distribution was 1.37$\pm$0.04. The degree of polarisation varied from $\sim$ 3% to 39% during our monitoring. We observed a change of $\sim$120 degrees in polarisation angle (PA) within $\sim$3 h on 2022 April 13. No clear correlation was found between optical flux and the degree of polarisation. Conclusions. The results of our optical flux, colour, and polarisation study provide hints that turbulence in the relativistic jet could be responsible for the intraday optical variations in the blazar S4 0954+65. However, the long-term flux variations may be caused by changes in the Doppler factor. △ Less

Submitted 8 October, 2023; originally announced October 2023.

Comments: 9 pages, 10 figures, 4 tables, 4 appendix, Astronomy & Astrophysics journal (in press)

Journal ref: A&A 679, A28 (2023)

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

Broad-line region in active galactic nuclei: Dusty or dustless?

Authors: Ashwani Pandey, Bozena Czerny, Swayamtrupta Panda, Raj Prince, Vikram Kumar Jaiswal, Mary Loli Martinez-Aldama, Michal Zajacek, Marzena Sniegowska

Abstract: Context. Dust in active galactic nuclei is clearly present right outside the broad-line region (BLR) in the form of a dusty molecular torus. However, some models of the BLR predict that dust may also exist within the BLR. Aims. We study the reprocessing of radiation by the BLR with the aim of observing how the presence of dust affects the reprocessed continuum and the line properties. Methods. We… ▽ More Context. Dust in active galactic nuclei is clearly present right outside the broad-line region (BLR) in the form of a dusty molecular torus. However, some models of the BLR predict that dust may also exist within the BLR. Aims. We study the reprocessing of radiation by the BLR with the aim of observing how the presence of dust affects the reprocessed continuum and the line properties. Methods. We calculated a range of models using the CLOUDY photoionisation code for dusty and dustless plasma. We paid particular attention to the well-studied object NGC 5548, and we compared the line equivalent width predictions with the data from observations for this object. Results. We obtained a rough agreement between the expected equivalent widths of the H$β$ and Mg II lines and the observed values for NGC 5548 for the line distances implied by the time-delay measurement (for H$β)$ and the radius-luminosity relation (for Mg II) when the medium is dusty. We found the incident radiation to be consistent with the radiation seen by the observer, so no shielding between the inner disc and the BLR is required. High ionisation lines such as He II, however, clearly form in the inner dustless region. When the additional absorber is present, the H$β$ emitting region moves closer to the dustless part of the accretion disc surface. △ Less

Submitted 25 October, 2023; v1 submitted 8 October, 2023; originally announced October 2023.

Comments: 10 pages, 7 figures, 1 table, accepted for publication in the Astronomy and Astrophysics journal

Journal ref: A&A 680, A102 (2023)

UV Fe II emission model of HE 0413-4031 and its relation to broad-line time delays

Authors: Michal Zajaček, Swayamtrupta Panda, Ashwani Pandey, Raj Prince, Alberto Rodríguez-Ardila, Vikram Jaiswal, Bożena Czerny, Krzysztof Hryniewicz, Maciej Urbanowicz, Piotr Trzcionkowski, Marzena Śniegowska, Zuzanna Fałkowska, Mary Loli Martínez-Aldama, Norbert Werner

Abstract: Fe II emission is a well-known contributor to the UV spectra of active galactic nuclei and the modeling of this part may affect the results obtained for the MgII$\lambda2800$ emission, which is one of the lines used for black hole mass measurements and cosmological applications. We use the 11-year monitoring of the selected quasar HE 0413-4031 with the South African Large Telescope (SALT) and we s… ▽ More Fe II emission is a well-known contributor to the UV spectra of active galactic nuclei and the modeling of this part may affect the results obtained for the MgII$\lambda2800$ emission, which is one of the lines used for black hole mass measurements and cosmological applications. We use the 11-year monitoring of the selected quasar HE 0413-4031 with the South African Large Telescope (SALT) and we supplement this monitoring with the near-IR spectrum taken with the SOAR telescope. A new redshift determination ($z=1.39117 \pm 0.00017$) using [OIII]$λλ4959,5007$ gave a very different value than the previous determination based only on the UV FeII pseudocontinuum ($z=1.3764$). It favors a different decomposition of the spectrum into Mg II and UV Fe II emissions. The line characteristics and the time delay of the Mg II emission ($224^{+21}_{-23}$ days) are not significantly affected. However, in comparison with the previous analysis, the rest-frame UV FeII time delay ($251^{+9}_{-7}$ days) is consistent with the inferred UV FeII line FWHM of $4200\,{\rm km/s}$ that is only slightly smaller than the MgII line FWHM. Hence the FeII-emitting material is more distant than the MgII-emitting gas in HE 0413-4031 by $\sim 0.023$ pc (4700 AU). The inferred velocity shift of both Mg II and UV Fe II lines with respect to the systemic redshift is now rather low, below 300 km s$^{-1}$. In addition, we construct an updated MgII radius-luminosity ($R-L$) relation from 194 sources, which is more than double the previous sample. The MgII $R-L$ relation is flatter than the UV FeII, optical FeII, and H$β$ $R-L$ relations. While the new decomposition of the spectrum is satisfactory, we see a need to create better Fe II templates using the newest version of the code CLOUDY. △ Less

Submitted 27 December, 2023; v1 submitted 5 October, 2023; originally announced October 2023.

Comments: 19 pages, 12+8 Figures, 5+1 Tables (Main Text+Appendix); accepted by A&A; minor corrections applied

Journal ref: Astronomy & Astrophysics, Volume 683, id.A140, 18 pp., 2024

Effects of heterogeneous data sets and time-lag measurement techniques on cosmological parameter constraints from MgII and CIV reverberation-mapped quasar data

Authors: Shulei Cao, Michal Zajaček, Bożena Czerny, Swayamtrupta Panda, Bharat Ratra

Abstract: Previously, we demonstrated that MgII and CIV reverberation-mapped quasars (RM QSOs) are standardizable and that the cosmological parameters inferred using the broad-line region radius-luminosity (R-L) relation are consistent with those determined from better-established cosmological probes. With more data expected from ongoing and future spectroscopic and photometric surveys, it is imperative to… ▽ More Previously, we demonstrated that MgII and CIV reverberation-mapped quasars (RM QSOs) are standardizable and that the cosmological parameters inferred using the broad-line region radius-luminosity (R-L) relation are consistent with those determined from better-established cosmological probes. With more data expected from ongoing and future spectroscopic and photometric surveys, it is imperative to examine how new QSO data sets of varied quality, with their own specific luminosity and time-delay distributions, can be best used to determine more restrictive cosmological parameter constraints. In this study, we test the effect of adding 25 OzDES MgII RM QSOs as well as 25 lower-quality SDSS RM CIV QSOs, which increases the previous sample of RM QSOs by $\sim 36\%$. Although cosmological parameter constraints become tighter for some cosmological models after adding these new QSOs, the new combined data sets have increased differences between R-L parameter values obtained in different cosmological models and thus a lower standardizability for the larger MgII + CIV compilation. Different time-delay methodologies, particularly the ICCF and CREAM methods used for inferring time delays of SDSS RM QSOs, slightly affect cosmological and R-L relation parameter values, however, the effect is negligible for (smaller) compilations of robust time-delay detections. Our analysis indicates that increasing the sample size is not sufficient for tightening cosmological constraints and a quality cut is necessary to obtain a standardizable RM QSO sample. △ Less

Submitted 8 February, 2024; v1 submitted 28 September, 2023; originally announced September 2023.

Comments: 25 pages, 14 figures, 6 tables; accepted by the MNRAS Main Journal

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 528, Issue 4, March 2024, Pages 6444-6469

Multiwavelength study of extreme variability in LEDA 1154204: A changing-look event in a type 1.9 Seyfert

Authors: T. Saha, A. Markowitz, D. Homan, M. Krumpe, S. Haemmerich, B. Czerny, M. Graham, S. Frederick, M. Gromadzki, S. Gezari, H. Winkler, D. A. H. Buckley, J. Brink, M. H. Naddaf, A. Rau, J. Wilms, A. Gokus, Z. Liu, I. Grotova

Abstract: Context. Multiwavelength studies of transients in actively accreting supermassive black holes have revealed that large-amplitude variability is frequently linked to significant changes in the optical spectra -- a phenomenon referred to as changing-look AGN (CLAGN). Aim. In 2020, the Zwicky Transient Facility detected a transient flaring event in the type-1.9 AGN 6dFGS~gJ042838.8-000040, wherein… ▽ More Context. Multiwavelength studies of transients in actively accreting supermassive black holes have revealed that large-amplitude variability is frequently linked to significant changes in the optical spectra -- a phenomenon referred to as changing-look AGN (CLAGN). Aim. In 2020, the Zwicky Transient Facility detected a transient flaring event in the type-1.9 AGN 6dFGS~gJ042838.8-000040, wherein a sharp increase in magnitude of $\sim$0.55 and $\sim$0.3 in the $g$- and $r$-bands, respectively, occurred over $\sim$40 days. Spectrum Roentgen Gamma (SRG)/eROSITA also observed the object in X-rays as part of its all-sky survey, but only after the flare had started decaying. Methods. We performed a three-year, multiwavelength follow-up campaign of the source to track its spectral and temporal characteristics. This campaign included multiple ground-based facilities for optical spectroscopic monitoring and space-based observatories including \textit{XMM-Newton} and \textit{Swift} for X-ray and UV observations. Results. An optical spectrum taken immediately after the peak revealed a changing-look event wherein the source had transitioned from type 1.9 to 1, with the appearance of a double-peaked broad H$β$ line and a blue continuum, both absent in an archival spectrum from 2005. The X-ray emission exhibits dramatic flux variation: a factor of $\sim$17, but with no spectral evolution, as the power-law photon index remained $\sim$1.9. There is no evidence of a soft X-ray excess. Overall the object exhibits no apparent signatures of a tidal disruption event. Conclusions. The transient event was likely triggered by a disk instability in a pre-existing accretion flow, culminating in the observed multi-wavelength variability and CLAGN event. △ Less

Submitted 16 September, 2023; originally announced September 2023.

Comments: 34 pages, 24 figures, Submitted to Astronomy & Astrophysics

Observations of a black hole X-ray binary indicate formation of a magnetically arrested disk

Authors: Bei You, Xinwu Cao, Zhen Yan, Jean-Marie Hameury, Bozena Czerny, Yue Wu, Tianyu Xia, Marek Sikora, Shuang-Nan Zhang, Pu Du, Piotr T. Zycki

Abstract: Accretion of material onto a black hole drags any magnetic fields present inwards, increasing their strength. Theory predicts that sufficiently strong magnetic fields can halt the accretion flow, producing a magnetically arrested disk (MAD). We analyze archival multi-wavelength observations of an outburst from the black hole X-ray binary MAXI J1820+070 in 2018. The radio and optical fluxes are del… ▽ More Accretion of material onto a black hole drags any magnetic fields present inwards, increasing their strength. Theory predicts that sufficiently strong magnetic fields can halt the accretion flow, producing a magnetically arrested disk (MAD). We analyze archival multi-wavelength observations of an outburst from the black hole X-ray binary MAXI J1820+070 in 2018. The radio and optical fluxes are delayed by about 8 and 17 days respectively, compared to the X-ray flux. We interpret this as evidence for the formation of a MAD. In this scenario, the magnetic field is amplified by an expanding corona, forming a MAD around the time of the radio peak. The optical delay is then due to thermal viscous instability in the outer disk. △ Less

Submitted 11 October, 2023; v1 submitted 31 August, 2023; originally announced September 2023.

Comments: The author's version of the article which will appear in Science on 31 August 2023, 49 pages including the extended data. The online publication version can be found at https://doi.org/10.1126/science.abo4504

Journal ref: Science 381, 961-964 (2023)

Science with a small two-band UV-photometry mission III: Active Galactic Nuclei and nuclear transients

Authors: M. Zajaček, B. Czerny, V. K. Jaiswal, M. Štolc, V. Karas, A. Pandey, D. R. Pasham, M. Śniegowska, V. Witzany, P. Suková, F. Münz, N. Werner, J. Řípa, J. Merc, M. Labaj, P. Kurfürst, J. Krtička

Abstract: In this review (the third in the series focused on a small two-band UV-photometry mission), we assess possibilities for a small UV two-band photometry mission in studying accreting supermassive black holes (SMBHs; mass range $\sim 10^6$-$10^{10}\,M_{\odot}$). We focus on the following observational concepts: (i) dedicated monitoring of selected type-I Active Galactic Nuclei (AGN) in order to measu… ▽ More In this review (the third in the series focused on a small two-band UV-photometry mission), we assess possibilities for a small UV two-band photometry mission in studying accreting supermassive black holes (SMBHs; mass range $\sim 10^6$-$10^{10}\,M_{\odot}$). We focus on the following observational concepts: (i) dedicated monitoring of selected type-I Active Galactic Nuclei (AGN) in order to measure the time delay between the far-UV, the near-UV, and other wavebands (X-ray and optical), (ii) nuclear transients including (partial) tidal disruption events and repetitive nuclear transients, and (iii) the study of peculiar sources, such as changing-look AGN, hollows and gaps in accretion disks, low-luminosity AGN, and candidates for Intermediate-Mass Black Holes (IMBHs; mass range $\sim 10^2$-$10^5\,M_{\odot}$) in galactic nuclei. For tidal disruption events (TDEs), high-cadence UV monitoring is crucial for distinguishing among different scenarios for the origin of the UV emission. The small two-band UV space telescope will also provide the information about the near- and far-UV continuum variability for rare transients, such as repetitive partial TDEs and jetted TDEs. We also discuss the possibilities to study and analyze sources with non-standard accretion flows, such as AGN with gappy disks, low-luminosity active galactic nuclei with intermittent accretion, and SMBH binaries potentially involving intermediate-mass black holes. △ Less

Submitted 12 March, 2024; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: 58 pages, 23 figures, 8 tables; accepted by Space Science Reviews

BAL effect in quasars due to source orientation

Authors: M. Sniegowska, M. H. Naddaf, M. L. Martinez-Aldama, P. Marziani, S. Panda, B. Czerny

Abstract: We investigated a scenario where the presence of a broad absorption line (BAL) feature in quasars (QSOs) is contingent upon the line of sight being situated within an outflow cone emanating from the source. We examined the mechanism of dust-driven winds based on the failed radiatively accelerated dusty outflow (FRADO) model proposed by Czerny & Hryniewicz, letting it be responsible for the formati… ▽ More We investigated a scenario where the presence of a broad absorption line (BAL) feature in quasars (QSOs) is contingent upon the line of sight being situated within an outflow cone emanating from the source. We examined the mechanism of dust-driven winds based on the failed radiatively accelerated dusty outflow (FRADO) model proposed by Czerny & Hryniewicz, letting it be responsible for the formation of massive outflow. We calculated the probability of observing the BAL effect from the geometry of outflow which is a function of global parameters of black hole mass (M$_{BH}$), Eddington ratio ($α_{Edd}$), and metallicity (Z). We then compared the results with prevalence of BAL QSOs in a sample of observational data from SDSS. The consistency of our model with the data supports the interpretation of the BAL phenomenon as a result of source orientation, rather than a transitory stage in AGN evolution △ Less

Submitted 5 June, 2023; originally announced June 2023.

Comments: 5 pages, 3 plots, Submitted to the Proceedings of IAU Symposium 378: Black hole winds at all scales

Effect of extinction on quasar luminosity distances determined from UV and X-ray flux measurements

Authors: Michal Zajaček, Bożena Czerny, Narayan Khadka, Mary Loli Martínez-Aldama, Raj Prince, Swayamtrupta Panda, Bharat Ratra

Abstract: In Khadka et al. (2023), a sample of X-ray-detected reverberation-mapped quasars was presented and applied for the comparison of cosmological constraints inferred using two well-established relations in AGN -- the X-ray/UV luminosity ($L_{X}-L_{UV}$) relation and the broad-line region radius-luminosity ($R-L$) relation. $L_{X}-L_{UV}$ and $R-L$ luminosity distances to the same quasars exhibit a di… ▽ More In Khadka et al. (2023), a sample of X-ray-detected reverberation-mapped quasars was presented and applied for the comparison of cosmological constraints inferred using two well-established relations in AGN -- the X-ray/UV luminosity ($L_{X}-L_{UV}$) relation and the broad-line region radius-luminosity ($R-L$) relation. $L_{X}-L_{UV}$ and $R-L$ luminosity distances to the same quasars exhibit a distribution of their differences that is generally asymmetric and positively shifted for the six cosmological models we consider. We demonstrate that this behaviour can be interpreted qualitatively to arise as a result of the dust extinction of UV/X-ray quasar emission. We show that the extinction always contributes to the non-zero difference between $L_{X}-L_{UV}$-based and $R-L$-based luminosity distances and we derive a linear relationship between the X-ray/UV colour index $E_{X-UV}$ and the luminosity-distance difference, which also depends on the value of the $L_{X}-L_{UV}$ relation slope. Taking into account the median and the peak values of the luminosity-distance difference distributions, the average X-ray/UV colour index falls in the range of $\overline{E}_{X-UV}=0.03-0.28$ mag for the current sample of 58 sources. This amount of extinction is typical for the majority of quasars and it can be attributed to the circumnuclear and interstellar media of host galaxies. After applying the standard hard X-ray and far-UV extinction cuts, heavily extincted sources are removed but overall the shift towards positive values persists. The effect of extinction on luminosity distances is more pronounced for the $L_{X}-L_{UV}$ relation since the extinction of UV and X-ray emissions both contribute. △ Less

Submitted 2 February, 2024; v1 submitted 14 May, 2023; originally announced May 2023.

Comments: 22 pages, 11 figures, 4 tables; accepted by AAS journals

Journal ref: The Astrophysical Journal, 961:229 (18pp), 2024 February 1

Wavelength-resolved reverberation mapping of intermediate redshift quasars HE 0413-4031 and HE 0435-4312: Dissecting Mg II, optical Fe II, and UV Fe II emission regions

Authors: Raj Prince, Michal Zajaček, S. Panda, K. Hryniewicz, V. K. Jaiswal, Bożena Czerny, P. Trzcionkowski, M. Bronikowski, M. Rałowski, C. S. Figaredo, M. L. Martinez-Aldama, M. Śniegowska, J. Średzińska, M. Bilicki, M-H Naddaf, A. Pandey, M. Haas, M. J. Sarna, G. Pietrzyński, V. Karas, A. Olejak, R. Przyłuski, R. R. Sefako, A. Genade, H. L. Worters , et al. (2 additional authors not shown)

Abstract: We present the wavelength-resolved reverberation mapping (RM) of combined MgII and UV FeII broad-line emissions for two intermediate redshifts (z$\sim$1), luminous quasars - HE 0413-4031 and HE 0435-4312, monitored by the SALT and 1-m class telescopes between 2012-2022. Through this technique, we aim to disentangle the Mg II and FeII emission regions and to build a radius-luminosity relation for U… ▽ More We present the wavelength-resolved reverberation mapping (RM) of combined MgII and UV FeII broad-line emissions for two intermediate redshifts (z$\sim$1), luminous quasars - HE 0413-4031 and HE 0435-4312, monitored by the SALT and 1-m class telescopes between 2012-2022. Through this technique, we aim to disentangle the Mg II and FeII emission regions and to build a radius-luminosity relation for UV FeII emission, which has so far remained unconstrained. Several methodologies have been applied to constrain the time delays for total MgII and FeII emissions. In addition, this technique is performed to quantify the inflow or outflow of broad-line region gas around the supermassive black hole and to disentangle the emission/emitting regions from lines produced in proximity to each other. The mean total FeII time delay is nearly equal to the mean total Mg II time delay for HE 0435-4312 suggesting a co-spatiality of their emissions. However, in HE 0413-4031, the mean FeII time delay is found to be longer than the mean MgII time delay, suggesting that FeII is produced at longer distances from the black hole. The UV Fe II R-L relation is updated with these two quasars and compared with the optical FeII relation, which suggests that the optical FeII region is located further than the UV FeII by a factor of 1.7-1.9 i.e. $R_{\rm FeII-opt}\sim(1.7-1.9)R_{\rm FeII-UV}$. We detected a weak pattern in the time delay vs. wavelength relation, suggesting that the MgII broad-line originates a bit closer to the SMBH than the UV FeII, however, the difference is not very significant. Comparison of MgII, UV, and optical FeII R-L relations suggests that the difference may be larger for lower-luminosity sources, possibly with the MgII emission originating further from the SMBH. In the future, more RM data will be acquired to put better constraints on these trends, in particular the UV FeII R-L relation. △ Less

Submitted 15 July, 2023; v1 submitted 26 April, 2023; originally announced April 2023.

Comments: 24 pages, 16 figures, 5 tables, Accepted for publication in A&A, in Press

Journal ref: A&A 678, A189 (2023)

Spectral Energy Distribution profiles from AGN accretion disc in multi-gap setup

Authors: Marcel Štolc, Michal Zajaček, Bożena Czerny, Vladimír Karas

Abstract: Spectral Energy Distribution (SED) of the broad-band continuum emission from black-hole accretion discs can serve as a tool to measure parameters of the central body and constrain the geometry of the inner accretion flow. We focus on the case of an active galactic nucleus (AGN), with an accretion disc dominating the UV/optical bands. We parameterize the changes in the thermal and power-law compone… ▽ More Spectral Energy Distribution (SED) of the broad-band continuum emission from black-hole accretion discs can serve as a tool to measure parameters of the central body and constrain the geometry of the inner accretion flow. We focus on the case of an active galactic nucleus (AGN), with an accretion disc dominating the UV/optical bands. We parameterize the changes in the thermal and power-law components, which can reveal the diminution of the emissivity. To this end we explore the effects of gaps in the accretion disc and the emerging SED that can be caused by the presence of either (i) the inner, optically thin, radiatively inefficient hot flow; (ii) a secondary black hole embedded within the accretion disc; or (iii) a combination of both components. We suggest that the resulting changes in the SED of the underlying continuum can help us to understand some departures from the standard-disc scenario. We estimate that the data required for such a project must be sampled in detail over the far-UV to soft X-ray bands during the interval of about a month corresponding to the characteristic variability timescale of an AGN. Detecting a gap at intermediate radii of a few 100 gravitational radii would require quality photometry with uncertainties up to $\sim$ 1%. The presence of the central cavity in the standard disc can be recovered in UV photometric data with an accuracy of 5% and better. We show the effect of the intrinsic reddening of the source and demonstrate when it can be disentangled. △ Less

Submitted 16 June, 2023; v1 submitted 6 April, 2023; originally announced April 2023.

Comments: 17 pages, 10 figures and 5 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Expectations for time-delay measurements in active galactic nuclei with the Vera Rubin Observatory

Authors: Bozena Czerny, Swayamtrupta Panda, Raj Prince, Vikram Kumar Jaiswal, Michal Zajacek, Mary Loli Martinez Aldama, Szymon Kozlowski, Andjelka B. Kovacevic, Dragana Ilic, Luka C. Popovic, Francisco Pozo Nunez, Sebastian F. Hoenig, William N. Brandt

Abstract: The Vera Rubin Observatory will provide an unprecedented set of time-dependent observations of the sky. The planned Legacy Survey of Space and Time (LSST) operating for 10 years will provide dense lightcurves for thousands of active galactic nuclei (AGN) in Deep Drilling Fields (DDFs) and less dense lightcurves for millions of AGN. We model the prospects for measuring time delays for emission line… ▽ More The Vera Rubin Observatory will provide an unprecedented set of time-dependent observations of the sky. The planned Legacy Survey of Space and Time (LSST) operating for 10 years will provide dense lightcurves for thousands of active galactic nuclei (AGN) in Deep Drilling Fields (DDFs) and less dense lightcurves for millions of AGN. We model the prospects for measuring time delays for emission lines with respect to the continuum, using these data. We model the artificial lightcurves using Timmer-Koenig algorithm, we use the exemplary cadence to sample them, we supplement lightcurves with the expected contamination by the strong emission lines (Hbeta, Mg II and CIV as well as with Fe II pseudo-continuum and the starlight). We choose the suitable photometric bands appropriate for the redshift and compare the assumed line time delay with the recovered time delay for 100 statistical realizations of the light curves. We show that time delays for emission lines can be well measured from the Main Survey for the bright tail of the quasar distribution (about 15% of all sources) with the accuracy within 1 sigma error, for DDFs results for fainter quasars are also reliable when all 10 years of data are used. There are also some prospects to measure the time delays for the faintest quasars at the smallest redshifts from the first two years of data, and eventually even from the first season. The entire quasar population will allow obtaining results of apparently high accuracy but in our simulations, we see a systematic offset between the assumed and recovered time delay depending on the redshift and source luminosity which will not disappear even in the case of large statistics. Such a problem might affect the slope of the radius-luminosity relation and cosmological applications of quasars if simulations correcting for such effects are not performed. △ Less

Submitted 15 June, 2023; v1 submitted 21 January, 2023; originally announced January 2023.

Comments: Astronomy & Astrophysics in press

Journal ref: A&A 675, A163 (2023)

Dissecting the broadband emission from γ-ray blazar PKS 0735+178 in search of neutrinos

Authors: Raj Prince, Saikat Das, Nayantara Gupta, Pratik Majumdar, Bożena Czerny

Abstract: The origin of the diffuse flux of TeV-PeV astrophysical neutrinos is still unknown. The $γ$-ray blazar PKS 0735+178, located outside the 90\% localization region at 2.2 deg from the best-fit IC-211208A event, was found to be flaring across all wavebands. In addition to leptonic synchrotron (SYN) and synchrotron self-Compton (SSC) emission, we invoke photohadronic ($pγ$) interactions inside the jet… ▽ More The origin of the diffuse flux of TeV-PeV astrophysical neutrinos is still unknown. The $γ$-ray blazar PKS 0735+178, located outside the 90\% localization region at 2.2 deg from the best-fit IC-211208A event, was found to be flaring across all wavebands. In addition to leptonic synchrotron (SYN) and synchrotron self-Compton (SSC) emission, we invoke photohadronic ($pγ$) interactions inside the jet to model the spectral energy distribution (SED) and neutrino emission. We analyze the 100 days $γ$-ray and X-ray data and 10 days around the neutrino event is chosen to generate the broadband SED. The temporal light curve indicates that the source was in a high state in optical, UV, $γ$-ray, and X-ray frequencies during the neutrino detection epoch. In the one-zone lepto-hadronic model, the SSC photons do not provide enough seed photons for $pγ$ interactions to explain the neutrino event. However, including an external photon field yields a neutrino event rate of 0.12 in 100 days, for the IceCube detector, using physically motivated values of the magnetic field, an external photon field peaking at optical wavelength, and other jet parameters. The radiation from secondary electrons at X-ray energies severely constrains the neutrino flux to a lower value than found in previous studies. Moreover, the flux of high-energy $γ$-rays at GeV energies from the decay of neutral pions is subdominant at the high-energy peak of the SED, suggesting a higher correlation of neutrinos flux with X-ray flux is plausible. △ Less

Submitted 6 December, 2023; v1 submitted 16 January, 2023; originally announced January 2023.

Comments: 9 pages, 5 figures, 1 table, accepted in mnras

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 527, Issue 3, January 2024, Pages 8746-8754

Quasar UV/X-ray relation luminosity distances are shorter than reverberation-measured radius-luminosity relation luminosity distances

Authors: Narayan Khadka, Michal Zajaček, Raj Prince, Swayamtrupta Panda, Bożena Czerny, Mary Loli Martínez-Aldama, Vikram Kumar Jaiswal, Bharat Ratra

Abstract: We use measurements of 59/58 quasars (QSOs), over a redshift range $0.0041\leq z \leq 1.686$, to do a comparative study of the radius--luminosity ($R-L$) and X-ray$-$UV luminosity ($L_X-L_{UV}$) relations and the implication of these relations for cosmological parameter estimation. By simultaneously determining $R-L$ or $L_X-L_{UV}$ relation parameters and cosmological parameters in six different… ▽ More We use measurements of 59/58 quasars (QSOs), over a redshift range $0.0041\leq z \leq 1.686$, to do a comparative study of the radius--luminosity ($R-L$) and X-ray$-$UV luminosity ($L_X-L_{UV}$) relations and the implication of these relations for cosmological parameter estimation. By simultaneously determining $R-L$ or $L_X-L_{UV}$ relation parameters and cosmological parameters in six different cosmological models, we find that both $R-L$ and $L_X-L_{UV}$ relations are standardizable but provide only weak cosmological parameter constraints, with $L_X-L_{UV}$ relation data favoring larger current non-relativistic matter density parameter $Ω_{m0}$ values than $R-L$ relation data and most other available data. We derive $L_X-L_{UV}$ and $R-L$ luminosity distances for each of the sources in the six cosmological models and find that $L_X-L_{UV}$ relation luminosity distances are shorter than $R-L$ relation luminosity distances as well as standard flat $Λ$CDM model luminosity distances. This explains why $L_X-L_{UV}$ relation QSO data favor larger $Ω_{m0}$ values than do $R-L$ relation QSO data or most other cosmological measurements. While our sample size is small and only spans a small $z$ range, these results indicate that more work is needed to determine whether the $L_X-L_{UV}$ relation can be used as a cosmological probe. △ Less

Submitted 4 April, 2023; v1 submitted 20 December, 2022; originally announced December 2022.

Comments: 20 pages, 10 figures, MNRAS in press

Journal ref: Mon.Not.Roy.Astron.Soc. 522 (2023) 1, 1247

Modeling photometric reverberation mapping data for the next generation of big data surveys. Quasar accretion disks sizes with the LSST

Authors: F. Pozo Nuñez, C. Bruckmann, S. Desamutara, B. Czerny, S. Panda, A. P. Lobban, G. Pietrzyński, K. L. Polsterer

Abstract: Photometric reverberation mapping can detect the radial extent of the accretion disc (AD) in Active Galactic Nuclei by measuring the time delays between light curves observed in different continuum bands. Quantifying the constraints on the efficiency and accuracy of the delay measurements is important for recovering the AD size-luminosity relation, and potentially using quasars as standard candles… ▽ More Photometric reverberation mapping can detect the radial extent of the accretion disc (AD) in Active Galactic Nuclei by measuring the time delays between light curves observed in different continuum bands. Quantifying the constraints on the efficiency and accuracy of the delay measurements is important for recovering the AD size-luminosity relation, and potentially using quasars as standard candles. We have explored the possibility of determining the AD size of quasars using next-generation Big Data surveys. We focus on the Legacy Survey of Space and Time (LSST) at the Vera C. Rubin Observatory, which will observe several thousand quasars with the Deep Drilling Fields and up to 10 million quasars for the main survey in six broadband filter during its 10-year operational lifetime. We have developed extensive simulations that take into account the characteristics of the LSST survey and the intrinsic properties of the quasars. The simulations are used to characterise the light curves from which AD sizes are determined using various algorithms. We find that the time delays can be recovered with an accuracy of 5 and 15% for light curves with a time sampling of 2 and 5 days, respectively. The results depend strongly on the redshift of the source and the relative contribution of the emission lines to the bandpasses. Assuming an optically thick and geometrically thin AD, the recovered time-delay spectrum is consistent with black hole masses derived with 30% uncertainty. △ Less

Submitted 24 January, 2023; v1 submitted 18 December, 2022; originally announced December 2022.

Comments: 19 pages, 17 figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society

Dust-driven wind as a model of Broad Absorption Line quasars

Authors: Mohammad Hassan Naddaf, Mary Loli Martinez-Aldama, Paola Marziani, Swayamtrupta Panda, Marzena Sniegowska, Bozena Czerny

Abstract: We test a scenario claiming that the broad absorption line (BAL) phenomenon in quasars (QSOs) is not a temporary stage of their life. In this scenario, we see the BAL effect only if the line of sight is within a spatially limited and collimated massive outflow cone covering only a fraction of sky from the point of view of the nucleus. The aim is to understand the theoretical mechanism behind the m… ▽ More We test a scenario claiming that the broad absorption line (BAL) phenomenon in quasars (QSOs) is not a temporary stage of their life. In this scenario, we see the BAL effect only if the line of sight is within a spatially limited and collimated massive outflow cone covering only a fraction of sky from the point of view of the nucleus. The aim is to understand the theoretical mechanism behind the massive outflow in BAL QSOs which is important for modelling the impact of quasars onto the star formation rate in the host galaxy, and, subsequently, onto the galaxy evolution. We apply the specific theoretical model of dust-driven wind. The model has considerable predictive power. The 2.5D version of FRADO model of Czerny & Hryniewicz gives rise to the formation of fast funnel-shaped outflow from the disk for a certain range of black hole masses, Eddington ratios and metallicities. We now interpret BAL QSO as sources viewed along the outflowing stream and we calculate the probabilities of seeing the BAL phenomenon as functions of these global parameters, and we compare these probabilities to those seen in the observational data. We include considerations on the presence/absence of obscuring torus. Comparing our theoretical results with observational data for a sample consisting of two sub-populations of BAL and non-BAL QSOs we found that both in the model and in the data the BAL phenomenon mostly happens for sources with black hole masses larger than 10^{8} solar mass, the effect get stronger with accretion rate, and also high metallicities are more likely in QSOs showing BAL features if the presence of torus is taken into account. The consistency of the model with the data provides support for the interpretation of the BAL phenomenon as the result of the orientation of the source. It also supports the underlying theoretical model although more consistency checks should be done in the future. △ Less

Submitted 21 June, 2023; v1 submitted 15 December, 2022; originally announced December 2022.

Comments: 15 pages, 13 figures, 3 tables, Accepted for publication in the Journal of Astronomy and Astrophysics (A&A)

Journal ref: A&A 675, A43 (2023)

Dusty plasma in active galactic nuclei

Authors: Bożena Czerny, Michal Zajaček, Mohammad-Hassan Naddaf, Marzena Sniegowska, Swayamtrupta Panda, Agata Różanska, Tek P. Adhikari, Ashwani Pandey, Vikram Kumar Jaiswal, Vladimír Karas, Abhijeet Borkar, Mary Loli Martínez-Aldama, Raj Prince

Abstract: Since many years we know that dust in the form of the dusty-molecular torus is responsible for the obscuration in active galactic nuclei (AGN) at large viewing angles and thus for the classification of AGN. Recently, we gained some observational and theoretical insight into geometry of the region and the role of the dust in the dynamics of the outflow and failed winds. We will briefly touch on all… ▽ More Since many years we know that dust in the form of the dusty-molecular torus is responsible for the obscuration in active galactic nuclei (AGN) at large viewing angles and thus for the classification of AGN. Recently, we gained some observational and theoretical insight into geometry of the region and the role of the dust in the dynamics of the outflow and failed winds. We will briefly touch on all these aspects, including our dust-based model (FRADO - Failed Radiatively Accelerated Dusty Outflow) of the formation of the Balmer lines in AGN. △ Less

Submitted 7 March, 2023; v1 submitted 20 November, 2022; originally announced November 2022.

Comments: Partially based on talk at 31st Symposium on the Physics of Ionized Gases in Belgrade, accepted for European Physical Journal D Topical Issue: "Physics of Ionized Gases and Spectroscopy of Isolated Complex Systems: Fundamentals and Applications"

The Wind Dynamics of Super-Eddington Sources in FRADO

Authors: Mohammad-Hassan Naddaf, Bożena Czerny, Michal Zajaček

Abstract: We perform non-hydrodynamical 2.5D simulations to study the dynamics of material above accretion disk based on the disk radiation pressure acting on dust. We assume a super-accreting underlying disk with the accretion rate of 10 times the Eddington rate with central black hole mass ranging from $10^7$ up to $10^9 M_{\odot}$. Such high accretion rates are characteristic for extreme sources. We show… ▽ More We perform non-hydrodynamical 2.5D simulations to study the dynamics of material above accretion disk based on the disk radiation pressure acting on dust. We assume a super-accreting underlying disk with the accretion rate of 10 times the Eddington rate with central black hole mass ranging from $10^7$ up to $10^9 M_{\odot}$. Such high accretion rates are characteristic for extreme sources. We show that for high accretors radiatively dust-driving mechanism based on FRADO model always leads to a massive outflow from the disk surface, and the failed wind develops only at larger radii. The outflow rate strongly depends on the black hole mass, and in optically-thick energy-driven solution can exceed the accretion rate for masses larger than $10^ 8 M_{\odot}$ but momentum-driven outflow does not exceed the accretion rate even for super-Eddington accretion, therefore not violating the adopted stationarity of the disk. However, even in this case the outflow from the disk implies a strong mechanical feedback. △ Less

Submitted 19 September, 2022; originally announced September 2022.

Comments: 13 pages, 3 figures, Accepted for publication in the journal of Dynamics (MDPI)

Journal ref: Dynamics 2022, 2(3), 295-305

Accretion disks, quasars and cosmology: meandering towards understanding

Authors: Bozena Czerny, Shulei Cao, Vikram Kumar Jaiswal, Vladimír Karas, Narayan Khadka, Mary Loli Martínez-Aldama, Mohammad Hassan Naddaf, Swayamtrupta Panda, Francisco Pozo Nuñez, Raj Prince, Bharat Ratra, Marzena Sniegowska, Zhefu Yu, Michal Zajaček

Abstract: As Setti & Woltjer noted back in 1973, quasars could be used to construct the Hubble diagram but the actual application was not that straightforward. It took years to implement the idea successfully. Most of the ways to use quasars for cosmology now require an advanced understanding of their structure, step by step. We briefly review this progress, with unavoidable personal bias, and concentrate o… ▽ More As Setti & Woltjer noted back in 1973, quasars could be used to construct the Hubble diagram but the actual application was not that straightforward. It took years to implement the idea successfully. Most of the ways to use quasars for cosmology now require an advanced understanding of their structure, step by step. We briefly review this progress, with unavoidable personal bias, and concentrate on bright unobscured sources. We will mention the problem of the gas flow character close to the innermost stable circular orbit close to the black hole, discussed 50 years ago, which later led to the development of the slim disk scenario, but was recently revived in the context of Magnetically Arrested Disks (MAD) and Standard and Normal Evolution (SANE) disk models. We also discuss the hot/warm corona issue, which is still under discussion and complicates the analysis of the X-ray reflection. We present the scenario of the formation of the low ionization part of the Broad Line Region as a mostly failed wind powered by radiation pressure acting on dust (FRADO - Failed Radiatively Driven Dusty Outflow model). Next, we discuss the cosmological constraints that are currently achievable with quasars, mostly concentrating on light echo methods (continuum time delays and spectral line time delays with respect to the continuum) which are (or should be) incorporating the progress mentioned above. Finally, we briefly mention future prospects in this direction. △ Less

Submitted 8 February, 2023; v1 submitted 14 September, 2022; originally announced September 2022.

Comments: Astrophysics and Space Science (in press)

Black-hole activity feedback across vast scales

Authors: Michal Zajaček, Bożena Czerny, Rainer Schödel, Norbert Werner, Vladimír Karas

Abstract: Both observational and theoretical studies of black-hole activity or active galactic nucleus (AGN) feedback have been ongoing since the first indication of supermassive black holes powering quasar activity in the 1960s. Although several crucial astrophysical questions have been answered in the following decades, a number of open problems remain, in particular how AGN feedback operates over nearly… ▽ More Both observational and theoretical studies of black-hole activity or active galactic nucleus (AGN) feedback have been ongoing since the first indication of supermassive black holes powering quasar activity in the 1960s. Although several crucial astrophysical questions have been answered in the following decades, a number of open problems remain, in particular how AGN feedback operates over nearly eight orders of magnitude - from scales of $\sim 10^{-3}\,{\rm pc}$ to the galaxy-cluster scales of a few hundred kiloparsecs. At the beginning of June 2022, about 50 junior as well as senior researchers met in Brno for the post-lockdown edition of the Cologne-Prague-Brno meeting to try to connect the dots. △ Less

Submitted 12 September, 2022; originally announced September 2022.

Comments: in print in Nature Astronomy as a meeting report; 5 pages, 1 figure

The LSST era of supermassive black holes accretion-disk reverberation mapping

Authors: Andjelka B. Kovacevic, Viktor Radovic, Dragana Ilic, Luka C. Popovic, Roberto J. Assef, Paula Sanchez-Saez, Robert Nikutta, Claudia M. Raiteri, Ilsang Yoon, Yasaman Homayouni, Yan-Rong Li, Neven Caplar, Bozena Czerny, Swayamtrupta Panda, Claudio Ricci, Isidora Jankov, Hermine Landt, Christian Wolf, Jelena Kovacevic-Dojcinovic, Masa Lakicevic, Djorđe Savic, Oliver Vince, Sasa Simic, Iva Cvorovic-Hajdinjak, Sladjana Marceta-Mandic

Abstract: The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will detect an unprecedentedly large sample of actively accreting supermassive black holes with typical accretion disk (AD) sizes of a few light days. This brings us to face challenges in the reverberation mapping (RM) measurement of AD sizes in active galactic nuclei (AGNs) using interband continuum delays. We examine the effect… ▽ More The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) will detect an unprecedentedly large sample of actively accreting supermassive black holes with typical accretion disk (AD) sizes of a few light days. This brings us to face challenges in the reverberation mapping (RM) measurement of AD sizes in active galactic nuclei (AGNs) using interband continuum delays. We examine the effect of LSST cadence strategies on AD RM using our metric AGNTimeLagMetric. It accounts for redshift, cadence, the magnitude limit, and magnitude corrections for dust extinction. Running our metric on different LSST cadence strategies, we produce an atlas of the performance estimations for LSST photometric RM measurements. We provide an upper limit on the estimated number of quasars for which the AD time lag can be computed within 0<z<7 using the features of our metric. We forecast that the total counts of such objects will increase as the mean sampling rate of the survey decreases. The AD time lag measurements are expected for >1000 sources in each Deep Drilling field (DDF, 10 sq. deg) in any filter, with the redshift distribution of these sources peaking at z~1. We find the LSST observation strategies with a good cadence (~ 5 days) and a long cumulative season (~9 yr), as proposed for LSST DDF, are favored for the AD size measurement. We create synthetic LSST light curves for the most suitable DDF cadences and determine RM time lags to demonstrate the impact of the best cadences based on the proposed metric. △ Less

Submitted 12 August, 2022; originally announced August 2022.

Comments: accepted for publication in the Astrophysical Journal Focus Issue on Rubin LSST cadence and survey strategy

Monitoring AGNs with H$β$ Asymmetry. III. Long-term Reverberation Mapping Results of 15 Palomar-Green Quasars

Authors: Dong-Wei Bao, Michael S. Brotherton, Pu Du, Jacob N. McLane, T. E. Zastrocky, Kianna A. Olson, Feng-Na Fang, Shuo Zhai, Zheng-Peng Huang, Kai Wang, Bi-Xuan Zhao, Sha-Sha Li, Sen Yang, Yong-Jie Chen, Jun-Rong Liu, Zhu-Heng Yao, Yue-Chang Peng, Wei-Jian Guo, Yu-Yang Songsheng, Yan-Rong Li, Bo-Wei Jiang, David H. Kasper, William T. Chick, My L. Nguyen, Jaya Maithil , et al. (20 additional authors not shown)

Abstract: In this third paper of the series reporting on the reverberation mapping (RM) campaign of active galactic nuclei with asymmetric H$β$ emission-line profiles, we present results for 15 Palomar-Green (PG) quasars using spectra obtained between the end of 2016 to May 2021. This campaign combines long time spans with relatively high cadence. For 8 objects, both the time lags obtained from the entire l… ▽ More In this third paper of the series reporting on the reverberation mapping (RM) campaign of active galactic nuclei with asymmetric H$β$ emission-line profiles, we present results for 15 Palomar-Green (PG) quasars using spectra obtained between the end of 2016 to May 2021. This campaign combines long time spans with relatively high cadence. For 8 objects, both the time lags obtained from the entire light curves and the measurements from individual observing seasons are provided. Reverberation mapping of 9 of our targets has been attempted for the first time, while the results for 6 others can be compared with previous campaigns. We measure the H$β$ time lags over periods of years and estimate their black hole masses. The long duration of the campaign enables us to investigate their broad line region (BLR) geometry and kinematics for different years by using velocity-resolved lags, which demonstrate signatures of diverse BLR geometry and kinematics. The BLR geometry and kinematics of individual objects are discussed. In this sample, the BLR kinematics of Keplerian/virialized motion and inflow is more common than outflow. △ Less

Submitted 1 July, 2022; originally announced July 2022.

Comments: 35 pages, 18 figures, 8 tables, accepted for publication in The Astrophysical Journal Supplement

Multiwavelength temporal and spectral study of TeV blazar 1ES 1727+502 during 2014 to 2021

Authors: Raj Prince, Rukaiya Khatoon, Pratik Majumdar, Bożena Czerny, Nayantara Gupta

Abstract: One of the most important questions in blazar physics is the origin of broadband emission and fast-flux variation. In this work, we studied the broadband temporal and spectral properties of a TeV blazar 1ES 1727+502 and explore the one-zone synchrotron-self Compton (SSC) model to fit the broadband spectral energy distribution (SED). We collected the long-term (2014-2021) multiband data which inclu… ▽ More One of the most important questions in blazar physics is the origin of broadband emission and fast-flux variation. In this work, we studied the broadband temporal and spectral properties of a TeV blazar 1ES 1727+502 and explore the one-zone synchrotron-self Compton (SSC) model to fit the broadband spectral energy distribution (SED). We collected the long-term (2014-2021) multiband data which includes both the low and high flux states of the source. The entire light curve is divided into three segments of different flux states and the best-fit parameters obtained by broadband SED modeling corresponding to three flux states were then compared. The TeV blazar 1ES 1727+502 has been observed to show the brightest flaring episode in X-ray followed by optical-UV and gamma-ray. The fractional variability estimated during various segments behaves differently in multiple wavebands, suggesting a complex nature of emission in this source. This source has shown a range of variability time from days scale to month scale during this long period of observations between 2014-2021. A "harder-when-brighter" trend is not prominent in X-ray but seen in optical-UV and an opposite trend is observed in gamma-ray. The complex nature of correlation among various bands is observed. The SED modeling suggests that the one-zone SSC emission model can reproduce the broadband spectrum in the energy range from optical-UV to very high energy gamma-ray. △ Less

Submitted 24 June, 2022; originally announced June 2022.

Comments: 14 pages, 8 figure, 2 table. Accepted for publication in MNRAS

Modelling time delays from two reprocessors in active galactic nuclei

Authors: Vikram Kumar Jaiswal, Raj Prince, Swayamtrupta Panda, Bożena Czerny

Abstract: Context. Continuum time delays from accretion disks in active galactic nuclei (AGN) has been proposed long time ago as a tool for measuring distances to the monitored sources. However, the method faces serious problems as a number of effects must be taken into account, including the contribution from the Broad Line Region (BLR). Aims. In this paper, we model the expected time delays when both the… ▽ More Context. Continuum time delays from accretion disks in active galactic nuclei (AGN) has been proposed long time ago as a tool for measuring distances to the monitored sources. However, the method faces serious problems as a number of effects must be taken into account, including the contribution from the Broad Line Region (BLR). Aims. In this paper, we model the expected time delays when both the disk reprocessing of the incident X-ray flux and further reprocessing by the BLR are included, with the aim to see if the two effects can be disentangled. Methods. We use simple response function for the accretion disk, without relativistic effects, and we use a parametric description to account for the BLR contribution. We include only scattering of the disk emission by the BLR inter-cloud medium. We also use artificial lightcurves with 1-day sampling to check if the effects are likely to be seen in real data. Results. We show that the effect of the BLR scattering on the predicted time delay is very similar to the effect of rising height of the X-ray source, without any BLR contribution. This brings additional degeneracy if we want in the future to recover the parameters of the system from the observed time delays in a specific object. Both effects, however, modify the slope of the delay versus wavelength curve when plotted in log space which opens a way to obtaining bare disk time delay needed for cosmology. In addition, when the disk irradiation is strong, the modification of the predicted delay by the BLR scattering and by X-ray source height become considerably different. However, the required data quality must be high since our artificial curves with 1-day sampling give rather noisy results. △ Less

Submitted 12 January, 2023; v1 submitted 24 June, 2022; originally announced June 2022.

Comments: 14 pages, 16 figures, and 1 table. Published in A&A. Comments are welcome

Standardizing reverberation-measured C IV time-lag quasars, and using them with standardized Mg II quasars to constrain cosmological parameters

Authors: Shulei Cao, Michal Zajaček, Swayamtrupta Panda, Mary Loli Martínez-Aldama, Bożena Czerny, Bharat Ratra

Abstract: We use 38 C IV quasar (QSO) reverberation-mapped (RM) observations, which span eight orders of magnitude in luminosity and the redshift range $0.001064 \leq z \leq 3.368$, to simultaneously constrain cosmological-model and QSO radius-luminosity ($R-L$) relation parameters in six cosmological models, using an improved technique that more correctly accounts for the asymmetric errors bars of the time… ▽ More We use 38 C IV quasar (QSO) reverberation-mapped (RM) observations, which span eight orders of magnitude in luminosity and the redshift range $0.001064 \leq z \leq 3.368$, to simultaneously constrain cosmological-model and QSO radius-luminosity ($R-L$) relation parameters in six cosmological models, using an improved technique that more correctly accounts for the asymmetric errors bars of the time-lag measurements. We find that $R-L$ relation parameters are independent of the cosmological models used in the analysis and so the $R-L$ relation can be used to standardize the C IV QSOs. The C IV QSO cosmological constraints are consistent with those from Mg II QSOs, allowing us to derive joint C IV + Mg II QSO cosmological constraints which are consistent with currently accelerated cosmological expansion, as well as consistent with cosmological constraints derived using better-established baryon acoustic oscillation (BAO) and Hubble parameter [$H(z)$] measurements. When jointly analyzed with $H(z)$ + BAO data, current C IV + Mg II QSO data mildly tighten current $H(z)$ + BAO data cosmological constraints. △ Less

Submitted 16 August, 2022; v1 submitted 31 May, 2022; originally announced May 2022.

Comments: 19 pages, 11 figures, MNRAS accepted version

Journal ref: MNRAS 516, 1721-1740 (2022)

Modified models of radiation pressure instability in application to 10, 10$^5$, and 10$^7$ $M_{\odot}$ accreting black holes

Authors: Marzena Śniegowska, Mikołaj Grzędzielski, Bożena Czerny, Agnieszka Janiuk

Abstract: Some of the accreting black holes exhibit much stronger variability patterns than the usual stochastic variations. Radiation pressure instability is one of the proposed mechanisms which could account for this effect. We aim to model luminosity changes for objects with black hole mass of 10, 10$^5$, and 10$^7$ solar masses, using the time-dependent evolution of an accretion disk unstable due to the… ▽ More Some of the accreting black holes exhibit much stronger variability patterns than the usual stochastic variations. Radiation pressure instability is one of the proposed mechanisms which could account for this effect. We aim to model luminosity changes for objects with black hole mass of 10, 10$^5$, and 10$^7$ solar masses, using the time-dependent evolution of an accretion disk unstable due to the dominant radiation pressure. We use a 1-dimensional, vertically integrated time-dependent numerical scheme which models simultaneous evolution of the disk and corona, coupled by the vertical mass exchange. We also discuss the possibility of presence of an inner optically thin flow, namely the Advection-Dominated Accretion Flow (ADAF). We found that the outburst character strongly depends on the magnetic field and the outer radius of the disk if this radius is smaller (due to TDE phenomenon) than the size of the instability zone in a stationary disk with infinite radius. For microquasars, the dependence on the magnetic field is monotonic, and the period decreases with the field strength. For larger black hole masses, the dependence is non-monotonic, and initial rise of the period is later replaced with the relatively rapid decrease as the magnetic field continues to rise. Still stronger magnetic field stabilizes the disk. Our computations confirm that the radiation pressure instability model can account for heartbeat states in microquasars. Rapid variability detected in IMBH in the form of Quasi-Periodic Ejection can be consistent with the model but only if combined with TDE phenomenon. Yearly repeating variability in Changing Look AGN also requires, in our model, small outer radius either due to the recent TDE or due to the presence of the gap in the disk related to the presence of a secondary black hole. △ Less

Submitted 21 April, 2022; originally announced April 2022.

Comments: 23 pages, 18 figures, Submitted for publication on Astronomy and Astrophysics, comments very welcome

Journal ref: A&A 672, A19 (2023)

Non-thermal emission from fall-back clouds in the Broad-Line Region of Active Galactic Nuclei

Authors: Ana Laura Müller, Mohammad-Hassan Naddaf, Michal Zajaček, Bożena Czerny, Anabella Araudo, Vladimír Karas

Abstract: The spectra of active galactic nuclei exhibit broad-emission lines that presumably originate in the Broad-Line Region (BLR) with gaseous-dusty clouds in a predominantly Keplerian motion around the central black hole. Signatures of both inflow and outflow motion are frequently seen. The dynamical character of BLR is consistent with the scenario that has been branded as the Failed Radiatively Accele… ▽ More The spectra of active galactic nuclei exhibit broad-emission lines that presumably originate in the Broad-Line Region (BLR) with gaseous-dusty clouds in a predominantly Keplerian motion around the central black hole. Signatures of both inflow and outflow motion are frequently seen. The dynamical character of BLR is consistent with the scenario that has been branded as the Failed Radiatively Accelerated Dusty Outflow (FRADO; Czerny & Hryniewicz 2011). In this scheme, frequent high-velocity impacts of BLR clouds falling back onto the underlying accretion disk are predicted. The impact velocities depend mainly on the black-hole mass, accretion rate, and metallicity and they range from a few km s$^{-1}$ up to thousands of km s$^{-1}$. Formation of strong shocks due to the collisions can give rise to the production of relativistic particles and associated radiation signatures. In this work, the non-thermal radiation generated in this process is investigated, and the spectral energy distributions for different parameter sets are presented. We find that the non-thermal processes caused by the impacts of clouds can lead to emission in the X-ray and the gamma-ray bands, playing the cloud density and metallicity a key role. △ Less

Submitted 11 April, 2022; originally announced April 2022.

Comments: 21 pages, 14 figures, and 2 tables. Accepted for publication in ApJ; extended Appendix

Cosmology Intertwined: A Review of the Particle Physics, Astrophysics, and Cosmology Associated with the Cosmological Tensions and Anomalies

Authors: Elcio Abdalla, Guillermo Franco Abellán, Amin Aboubrahim, Adriano Agnello, Ozgur Akarsu, Yashar Akrami, George Alestas, Daniel Aloni, Luca Amendola, Luis A. Anchordoqui, Richard I. Anderson, Nikki Arendse, Marika Asgari, Mario Ballardini, Vernon Barger, Spyros Basilakos, Ronaldo C. Batista, Elia S. Battistelli, Richard Battye, Micol Benetti, David Benisty, Asher Berlin, Paolo de Bernardis, Emanuele Berti, Bohdan Bidenko , et al. (178 additional authors not shown)

Abstract: In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of system… ▽ More In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $σ_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,σ$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $Ω_m$, and the amplitude or rate of the growth of structure ($σ_8,fσ_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged] △ Less

Submitted 24 April, 2022; v1 submitted 11 March, 2022; originally announced March 2022.

Comments: Contribution to Snowmass 2021. 224 pages, 27 figures. Accepted for publication in JHEAp

Journal ref: J. High En. Astrophys. 2204, 002 (2022)

Spectropolarimetry and spectral decomposition of high-accreting Narrow Line Seyfert 1 galaxies

Authors: Marzena Śniegowska, Swayamtrupta Panda, Bożena Czerny, Đorge Savić, Mary Loli Martínez-Aldama, Paola Marziani, Jian-Min Wang, Pu Du, Luka Č. Popović, Chandra Shekhar Saraf

Abstract: Narrow Line Seyfert 1 (NLSy1) galaxies have been shown to have high Eddington ratios and relatively small black hole masses. The measurement of the black hole masses is based on the virial relation which is dependent on the distribution of the line-emitting gas and the viewing angle to the source. Spectropolarimetry enables us to probe the geometry of this line-emitting gas and allows us to estima… ▽ More Narrow Line Seyfert 1 (NLSy1) galaxies have been shown to have high Eddington ratios and relatively small black hole masses. The measurement of the black hole masses is based on the virial relation which is dependent on the distribution of the line-emitting gas and the viewing angle to the source. Spectropolarimetry enables us to probe the geometry of this line-emitting gas and allows us to estimate independently the viewing angle of the source by comparing the spectrum viewed under natural light and in the polarized light. We performed spectropolarimetric observations of three NLSy1 - Mrk 1044, SDSS J080101.41+184840.7, and IRAS 04416+1215 using the European Southern Observatory's Very Large Telescope. We use ESO Reflex workflow to perform standard data reduction and extract the natural and polarized spectra. We estimate the Stokes parameters and the viewing angles of the three sources. We model the Stokes parameters and infer the properties of the scattering media - located in the equatorial and polar regions, and simulate the spectra observed both in natural light and in polarized light using the polarization radiative transfer code STOKES. We confirm that all three sources are high Eddington ratio objects. We are successful in recovering the observed H$α$ line profile both in the natural and polarized light using the STOKES modelling. We recover the polarization fractions of the order of 0.2-0.5% for the three sources. Our principal component analysis shows that the sample of the 25 sources including our sources, Fairall 9 from Jiang et al. (2021), and sources from Capetti et al. (2021) are mainly driven by the black hole mass and Eddington ratio. We re-affirm the connection of the strength of the optical FeII emission with the Eddington ratio, but the dependence on the viewing angle is moderate, more like a secondary effect. △ Less

Submitted 15 August, 2023; v1 submitted 28 February, 2022; originally announced February 2022.

Comments: 25 pages, 21 figures, 11 tables, to be published in Astronomy & Astrophysics, version after proofs

Journal ref: A&A 678, A63 (2023)

Wavelength-resolved Reverberation Mapping of quasar CTSC30.10: Dissecting MgII and FeII emission regions

Authors: Raj Prince, Michal Zajaček, Bożena Czerny, P. Trzcionkowski, M. Bronikowski, C. S. Figaredo, S. Panda, M. L. Martinez-Aldama, K. Hryniewicz, V. K. Jaiswal, M. Śniegowska, M. Naddaf, M. Bilicki, M. Haas, M. J. Sarna, V. Karas, A. Olejak, R. Przyłuski, M. Rałowski, A. Udalski, R. R. Sefako, A. Genade, H. L. Worters

Abstract: We present the results of the reverberation monitoring aimed at MgII broad line and FeII pseudocontinuum for the luminous quasar CTS C30.10 (z = 0.90052) with the Southern African Large Telescope covering the years 2012-2021. We aimed at disentangling the MgII and UV FeII variability and the first measurement of UV FeII time delay for a distant quasar. We used several methods for time-delay measur… ▽ More We present the results of the reverberation monitoring aimed at MgII broad line and FeII pseudocontinuum for the luminous quasar CTS C30.10 (z = 0.90052) with the Southern African Large Telescope covering the years 2012-2021. We aimed at disentangling the MgII and UV FeII variability and the first measurement of UV FeII time delay for a distant quasar. We used several methods for time-delay measurements and determined both FeII and MgII time delays as well as performed a wavelength-resolved time delay study for a combination of MgII and FeII in the 2700 - 2900 Årestframe wavelength range. We obtain the time delay for MgII of $275.5^{+12.4}_{-19.5}$ days in the rest frame, while for FeII we have two possible solutions of $270.0^{+13.8}_{-25.3}$ days and $180.3^{+26.6}_{-30.0}$ in the rest frame. Combining this result with the old measurement of FeII UV time delay for NGC 5548 we discuss for the first time the radius-luminosity relation for UV FeII with the slope consistent with $0.5$ within uncertainties. Since FeII time delay has a shorter time-delay component but lines are narrower than MgII, we propose that the line delay measurement is biased towards the BLR part facing the observer, with the bulk of the Fe II emission may arise from the more distant BLR region, one that is shielded from the observer. △ Less

Submitted 26 January, 2022; originally announced January 2022.

Comments: 22 pages, 19 Figures, 6 Tables, Submitted to Astronomy and Astrophysics, Comments are welcome

Journal ref: A&A 667, A42 (2022)

Mass loss rate of accretion disk in FRADO

Authors: Mohammad-Hassan Naddaf, Bozena Czerny

Abstract: We have developed the 2.5D version of the basic physically motivated 1D model of Czerny & Hryniewicz (2011), i.e. Failed Radiatively Accelerated Dusty Outflow (FRADO) model. This model is based on the idea that radiation pressure acting on dust is responsible for the formation of the low ionized part of the Broad Line Region (BLR). Such radiation pressure is strong enough to form a fast outflow fr… ▽ More We have developed the 2.5D version of the basic physically motivated 1D model of Czerny & Hryniewicz (2011), i.e. Failed Radiatively Accelerated Dusty Outflow (FRADO) model. This model is based on the idea that radiation pressure acting on dust is responsible for the formation of the low ionized part of the Broad Line Region (BLR). Such radiation pressure is strong enough to form a fast outflow from the disk surface in the inner part of low ionized BLR. The outflow properties depend on the basic physical parameters, like black hole mass, Eddington ratio and gas metallicity. We here aim at estimating the disk mass loss rate due to this process, and comparing the results with outflows detected in Broad Absorption Line (BAL) quasars. △ Less

Submitted 15 December, 2021; originally announced December 2021.

Comments: 3 pages, 1 figure, submitted to the Proceedings of the 40th Polish Astronomical Society Meeting, 13-17 September 2021, Szczecin, Poland

Do reverberation-measured H$β$ quasars provide a useful test of cosmology?

Authors: Narayan Khadka, Mary Loli Martínez-Aldama, Michal Zajaček, Bożena Czerny, Bharat Ratra

Abstract: We use 118 H$β$ quasar (QSO) observations in the redshift range $0.0023 \leq z \leq 0.89$ to simultaneously constrain cosmological model parameters and QSO 2-parameter radius-luminosity ($R-L$) relation parameters in six different cosmological models. We find that the $R-L$ relation parameters for these QSOs are independent of the assumed cosmology so these QSOs seem to be standardizable through t… ▽ More We use 118 H$β$ quasar (QSO) observations in the redshift range $0.0023 \leq z \leq 0.89$ to simultaneously constrain cosmological model parameters and QSO 2-parameter radius-luminosity ($R-L$) relation parameters in six different cosmological models. We find that the $R-L$ relation parameters for these QSOs are independent of the assumed cosmology so these QSOs seem to be standardizable through the $R-L$ relation (although there is a complication that might render this untrue). Cosmological constraints obtained using these QSOs are weak, more favor currently decelerated cosmological expansion, and typically are in $\sim 2σ$ tension with those obtained from a joint analysis of baryon acoustic oscillation and Hubble parameter measurements. Extending the $R-L$ relation to a 3-parameter one to try to correct for the accretion rate effect does not result in a reduction of the cosmological constraints discrepancy nor does it result in the hoped-for significant reduction of the intrinsic scatter of the $R-L$ relation. △ Less

Submitted 12 May, 2022; v1 submitted 30 November, 2021; originally announced December 2021.

Comments: 20 pages, 7 figures, H$β$ QSO data are given

Journal ref: Mon.Not.Roy.Astron.Soc. 513 (2022) 2, 1985

Radiation pressure on dust explains the Low Ionized Broad Emission Lines in Active Galactic Nuclei

Authors: Mohammad-Hassan Naddaf, Bozena Czerny

Abstract: Broad emission lines are the most characteristic features in the spectra of active galaxies. They mostly show either a single-peaked or double-peaked profiles; and originate from a complex dynamics of the likely discrete clouds moving in a spatially extended region so-called Broad Line Region (BLR). We follow a non-hydrodynamical single-cloud approach to the BLR dynamics based on the dust-driving… ▽ More Broad emission lines are the most characteristic features in the spectra of active galaxies. They mostly show either a single-peaked or double-peaked profiles; and originate from a complex dynamics of the likely discrete clouds moving in a spatially extended region so-called Broad Line Region (BLR). We follow a non-hydrodynamical single-cloud approach to the BLR dynamics based on the dust-driving model of Czerny & Hryniewicz. We previously showed in details that the 2.5D version of the model could provide us with the 3D geometry of the BLR. In this paper, we provide a large grid of results based on which we aim at testing the model with calculation of the spectral line generic profiles. We show that the shape of profiles not only depends on the accretion rate of the source, the black hole mass, and the viewing angle, but also it is most significantly affected by the adopted dust-to-gas mass ratio that regulates the strength of the radiation pressure. We also show that the model can appropriately explain the low ionized broad emission lines of the mean spectrum of quasars, such as MgII and Hbeta. △ Less

Submitted 13 April, 2022; v1 submitted 29 November, 2021; originally announced November 2021.

Comments: 11 pages, 6 figures, Accepted for publication in A&A

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

The Main Sequence of quasars: the taming of the extremes

Authors: P. Marziani, E. Bon, N. Bon, M. D'Onofrio, B. Punsly, M. Sniegowska, B. Czerny, S. Panda, M. L. Martinez Aldama, A. del Olmo, A. Deconto-Machado, C. A. Negrete, D. Dultzin, T. Buendia, K. Garnica

Abstract: The last few years have seen the confirmation of several trends associated with the quasar main sequence. The idea of a main sequence for quasars is relatively recent, and its full potential for the observational classification and contextualization of quasar properties has yet to be fully exploited. The main sequence drivers are discussed in terms of the properties of extreme objects. We briefly… ▽ More The last few years have seen the confirmation of several trends associated with the quasar main sequence. The idea of a main sequence for quasars is relatively recent, and its full potential for the observational classification and contextualization of quasar properties has yet to be fully exploited. The main sequence drivers are discussed in terms of the properties of extreme objects. We briefly summarize developments that constrain the viewing angle of the accretion disk in a particular class of quasars (extreme Population B, radiating at low Eddington ratio), as well as inferences on the chemical composition of the broad line emitting gas, and on the nature of radio emission along the quasar main sequence. △ Less

Submitted 7 November, 2021; originally announced November 2021.

Comments: Paper presented at the XIII SCSLSA, to appear in Astronomische Nachrichten

Modeling changing-look (CL) AGN phenomenon in 1D using accretion disk instabilities

Authors: Marzena Sniegowska, Mikolaj Grzedzielski, Bozena Czerny, Agnieszka Janiuk

Abstract: Apart from regular, low-level stochastic variability, some AGN occasionally show exceptionally large changes in the luminosity, spectral shape, and/or X-ray absorption. The most notable are the changes of the spectral type when the source classified as a Seyfert 1 becomes a Seyfert 2 galaxy or vice versa. Thus a name was coined of 'Changing-Look AGN' (CL AGN). The origin of this phenomenon is stil… ▽ More Apart from regular, low-level stochastic variability, some AGN occasionally show exceptionally large changes in the luminosity, spectral shape, and/or X-ray absorption. The most notable are the changes of the spectral type when the source classified as a Seyfert 1 becomes a Seyfert 2 galaxy or vice versa. Thus a name was coined of 'Changing-Look AGN' (CL AGN). The origin of this phenomenon is still unknown, but for most of the sources, there are strong arguments in favor of the intrinsic changes. Understanding the nature of such rapid changes is a challenge to the models of black hole accretion flows since the timescales of the changes are much shorter than the standard disk viscous timescales. We aim to model the CL AGN phenomenon assuming that the underlying mechanism is the time-dependent evolution of a black hole accretion disk unstable due to the dominant radiation pressure. We use a 1-dimensional, vertically integrated disk model, but we allow for the presence of the hot coronal layer above the disk and the presence of the inner purely hot flow. We focus on the variability timescales and amplitudes, which can be regulated by the action of large-scale magnetic fields, the description of the disk-corona coupling, and the presence of an inner optically thin flow, like Advection-Dominated Accretion Flow (ADAF). We compare model predictions for the accretion disk around black hole mass 10$^7$M$_{\odot}$. △ Less

Submitted 25 October, 2021; originally announced October 2021.

Comments: 7 pages, 5 figures, accepted for publication as part of the Special Issue featuring the "13th SCSLSA Workshop" in Astronomische Nachrichten

Linear Spectropolarimetric Analysis of Fairall 9 with VLT/FORS2

Authors: Bo-Wei Jiang, Paola Marziani, Đorđe Savić, Elena Shablovinskaya, Luka Č. Popović, Victor L. Afanasiev, Bożena Czerny, Jian-Min Wang, Ascensión del Olmo, Mauro D'Onofrio, Marzena Śniegowska, Paola Mazzei, Swayamtrupta Panda

Abstract: The quasar Main Sequence (MS) appears to be an incredibly powerful tool to organize the diversity in large samples of type-1 quasars but the most important physical parameters governing it are still unclear. Here we investigate the origin of the broadening and of a defining feature of Population B sources: a strong redward asymmetry of the Balmer emission lines. We focus on a prototypical source,… ▽ More The quasar Main Sequence (MS) appears to be an incredibly powerful tool to organize the diversity in large samples of type-1 quasars but the most important physical parameters governing it are still unclear. Here we investigate the origin of the broadening and of a defining feature of Population B sources: a strong redward asymmetry of the Balmer emission lines. We focus on a prototypical source, Fairall 9. Spectropolarimetric data of the Fairall 9 broad H$β$ and H$α$ profiles allowed for a view of the geometric and dynamical complexity of the line emitting regions. Measurements (1) provided evidence of rotational motion; (2) were helpful to test the presence of polar and equatorial scatterers, and their association with non-virial motions. △ Less

Submitted 2 August, 2021; originally announced August 2021.

Comments: 22 pages, 3 tables and 14 figures. Accepted for publication in MNRAS

Standardizing reverberation-measured Mg II time-lag quasars, by using the radius-luminosity relation, and constraining cosmological model parameters

Authors: Narayan Khadka, Zhefu Yu, Michal Zajaček, Mary Loli Martinez-Aldama, Bożena Czerny, Bharat Ratra

Abstract: We use 78 reverberation-measured Mg II time-lag quasars (QSOs) in the redshift range $0.0033 \leq z \leq 1.89$ to constrain cosmological parameters in six different cosmological models. The basis of our method is the use of the radius-luminosity or $R-L$ relation to standardize these 78 Mg II QSOs. In each cosmological model we simultaneously determine $R-L$ relation and cosmological model paramet… ▽ More We use 78 reverberation-measured Mg II time-lag quasars (QSOs) in the redshift range $0.0033 \leq z \leq 1.89$ to constrain cosmological parameters in six different cosmological models. The basis of our method is the use of the radius-luminosity or $R-L$ relation to standardize these 78 Mg II QSOs. In each cosmological model we simultaneously determine $R-L$ relation and cosmological model parameters, thus avoiding the circularity problem. We find that the $R-L$ relation parameter values are independent of the cosmological model used in the analysis thus establishing that current Mg II QSOs are standardizable candles. Cosmological constraints obtained using these QSOs are significantly weaker than, but consistent with, those obtained from a joint analysis of baryon acoustic oscillation (BAO) observations and Hubble parameter [$H(z)$] measurements. So, we also analyse these QSOs in conjunction with the BAO + $H(z)$ data and find cosmological constraints consistent with the standard spatially-flat $Λ$CDM model as well as with mild dark energy dynamics and a little spatial curvature. A larger sample of higher-quality reverberation-measured QSOs should have a smaller intrinsic dispersion and so should provide tighter constraints on cosmological parameters. △ Less

Submitted 26 September, 2021; v1 submitted 21 June, 2021; originally announced June 2021.

Comments: 16 pages, 11 figure, Mg II QSO data are given, MNRAS accepted version. arXiv admin note: text overlap with arXiv:2105.12692

Journal ref: Mon.Not.Roy.Astron.Soc. 508 (2021) 4, 4722

Viewing angle observations and effects of evolution with redshift, black hole mass, and Eddington ratio in quasar based cosmology

Authors: Raj Prince, Krzysztof Hryniewicz, Swayamtrupta Panda, Bożena Czerny, Agnieszka Pollo

Abstract: This study is focused on the observational measurement of the viewing angle of individual quasars by modeling the broadband quasar spectrum ranging from the infra-red (IR) to the soft X-ray band. Sources are selected from various published catalogs, and their broadband quasi-simultaneous spectral data points were collected and used to model. We started with a COSMOS sample of type-1 sources which… ▽ More This study is focused on the observational measurement of the viewing angle of individual quasars by modeling the broadband quasar spectrum ranging from the infra-red (IR) to the soft X-ray band. Sources are selected from various published catalogs, and their broadband quasi-simultaneous spectral data points were collected and used to model. We started with a COSMOS sample of type-1 sources which have broadband photometric points. Then, to include more data points, we cross-matched the COSMOS with the SDSS DR14 quasar catalog, and eventually, we find 90 sources that have broadband data ranging from IR to soft X-ray. The broadband spectral energy distribution (SED) modeling is done in Xspec by using the optxagnf and the SKIRTOR models for the X-ray, UV, Optical, and IR regimes for each source. The whole sample is divided into four bins with respect to redshift, black hole (BH) mass, and Eddington ratio with an equal number of sources in each bin. The viewing angle is estimated in each bin, and its evolution with respect to redshift, BH mass, and Eddington ratio is examined. As a result, we did not find any significant evolution of viewing angle with those parameters within the 95$\%$ confidence interval. We conclude that the use of quasars in cosmology to determine the expansion rate of the universe is therefore justified, and biases are not expected. △ Less

Submitted 2 December, 2021; v1 submitted 7 June, 2021; originally announced June 2021.

Comments: 13 pages, 4 figures, 2 tables, Accepted for publication in ApJ

Multi-wavelength Analysis and Modeling of OJ 287 During 2017-2020

Authors: Raj Prince, Aditi Agarwal, Nayantara Gupta, Pratik Majumdar, Bożena Czerny, Sergio A. Cellone, I. Andruchow

Abstract: Context. The blazar OJ 287 has been proposed as a binary black hole system based on its periodic optical outburst. Among blazars with parsec scale jets, the black hole binary systems are very rare and hence this source is very interesting to study. Aims. The BL Lac OJ 287 is an interesting object for multi-wavelength study due to its periodic outbursts. We have analyzed the optical, X-ray, and gam… ▽ More Context. The blazar OJ 287 has been proposed as a binary black hole system based on its periodic optical outburst. Among blazars with parsec scale jets, the black hole binary systems are very rare and hence this source is very interesting to study. Aims. The BL Lac OJ 287 is an interesting object for multi-wavelength study due to its periodic outbursts. We have analyzed the optical, X-ray, and gamma-ray data of OJ 287 for the period of 2017-2020. There are several high states in optical-UV and X-ray frequencies during this period. Based on the observed variability in optical and X-ray frequencies the entire period 2017-2020 is divided into five segments, referred to as A, B, C, D, & E in this paper. A detailed temporal and spectral analysis is performed to understand the nature of its flaring activities. Methods. To understand the temporal variability in this source we have studied the intra-day, and fractional variability for all the various states, and along with that fast variability time was also estimated to understand the nature of variability. Further, the multi-wavelength SED modeling is performed to know more about the physical processes responsible for the simultaneous broadband emission and the fast variability. Results. The Fermi-LAT observations show a moderate flux level of this source in gamma-ray frequency throughout this period, though flux variability has been observed. The source has shown a strong flux variability in X-ray, optical, and UV during early 2017 and mid-2020 when the source was in a very high state. A single zone SSC emission model is considered to model the spectral energy distributions and this helps us to explore the nature of this BL Lac with binary super-massive black holes. △ Less

Submitted 9 May, 2021; originally announced May 2021.

Comments: 14 pages, 5 figures, 6 tables, Accepted for publication in A&A

Journal ref: A&A 654, A38 (2021)

A comprehensive study of 2019-2020 flare of OJ 287 using AstroSat, Swift, and NuSTAR

Authors: Raj Prince, Gayathri Raman, Rukaiya Khatoon, Aditi Agarwal, Varun, Nayantara Gupta, Bożena Czerny, Pratik Majumdar

Abstract: OJ 287 is a well-studied binary black hole system, that occasionally exhibits bright X-ray and optical flares. Here we present a detailed spectral study of its second brightest X-ray flare observed during 2019-2020 using archival Swift and NuSTAR observations along with ToO observations from AstroSat. The entire flaring period is divided into three states, defined as low, intermediate, and high st… ▽ More OJ 287 is a well-studied binary black hole system, that occasionally exhibits bright X-ray and optical flares. Here we present a detailed spectral study of its second brightest X-ray flare observed during 2019-2020 using archival Swift and NuSTAR observations along with ToO observations from AstroSat. The entire flaring period is divided into three states, defined as low, intermediate, and high states. The variation of hardness ratio (HR) with 0.3-10.0 keV integrated flux suggests a \enquote{softer-when-brighter} behavior, as also previously reported based on flux-index variations. Simultaneous high state X-ray spectra obtained using Swift, NuSTAR and AstroSat are very steep with a power-law index $>$2. A significant spectral change is observed in AstroSat-SXT and LAXPC spectrum which is consistent with Swift-XRT and NuSTAR spectrum. Together, optical-UV and X-ray spectrum during the high flux state, suggesting the emergence of a new high BL Lacertae (HBL) component. We have modeled the synchrotron peak with publicly available code named GAMERA for low, intermediate, and high flux state. Our modeling suggests the need of high magnetic field to explain the high state under the leptonic scenarios. Increase in the magnetic field value inside the jet could be linked to the increase in accretion rate as expected in the BH-disk impact scenario. The color-magnitude diagram reveals a \enquote{bluer-when-brighter} spectral energy distribution chromatism during the flaring period. Different chromatism or no chromatism at various occasions suggests a complex origin of optical emission, which is believed to be produced by disc impact or through synchrotron emission in the jet. △ Less

Submitted 7 September, 2021; v1 submitted 9 May, 2021; originally announced May 2021.

Comments: 14 pages, 13 figures, 3 tables, Accepted in MNRAS