Search

The local (z=0.0315) AGN Mrk 817, was monitored over more than 500 days with space-borne and ground-based instruments as part of a large international campaign AGN STORM 2. Here, we present a comprehensive analysis of the broad-band continuum variations using detailed modeling of the broad line region (BLR), several types of disk winds classified by their optical depth, and new numerical simulations. We find that diffuse continuum (DC) emission, with additional contributions from strong and broad emission lines, can explain the continuum lags observed in this source during high and low luminosity phases. Disk illumination by the variable X-ray corona contributes only a small fraction of the observed continuum lags. Our BLR models assume radiation pressure-confined clouds distributed over a distance of 2-122 light days. We present calculated mean-emissivity radii of many emission lines, and DC emission, and suggest a simple, transfer-function-dependent method that ties them to cross-correlation lag determinations. We do not find clear indications for large optical depth winds but identify the signature of lower column density winds. In particular, we associate the shortest observed continuum lags with a combination of tau(1 Ryd) approx. 2 wind and a partly shielded BLR. Even smaller optical depth winds may be associated with X-ray absorption features and with noticeable variations in the width and lags of several high ionization lines like HeII and CIV. Finally, we demonstrate the effect of torus dust emission on the observed lags in the i and z bands., Comment: 23 pages, 18 figures. Corrected typographical error in the title of the paper as it appeared in the Metadata

X-ray reverberation mapping is a powerful technique for probing the innermost accretion disk, whereas continuum reverberation mapping in the UV, optical, and infrared (UVOIR) reveals reprocessing by the rest of the accretion disk and broad-line region (BLR). We present the time lags of Mrk 817 as a function of temporal frequency measured from 14 months of high-cadence monitoring from Swift and ground-based telescopes, in addition to an XMM-Newton observation, as part of the AGN STORM 2 campaign. The XMM-Newton lags reveal the first detection of a soft lag in this source, consistent with reverberation from the innermost accretion flow. These results mark the first simultaneous measurement of X-ray reverberation and UVOIR disk reprocessing lags$\unicode{x2013}$effectively allowing us to map the entire accretion disk surrounding the black hole. Similar to previous continuum reverberation mapping campaigns, the UVOIR time lags arising at low temporal frequencies are longer than those expected from standard disk reprocessing by a factor of 2-3. The lags agree with the anticipated disk reverberation lags when isolating short-timescale variability, namely timescales shorter than the H$\beta$ lag. Modeling the lags requires additional reprocessing constrained at a radius consistent with the BLR size scale inferred from contemporaneous H$\beta$-lag measurements. When we divide the campaign light curves, the UVOIR lags show substantial variations, with longer lags measured when obscuration from an ionized outflow is greatest. We suggest that, when the obscurer is strongest, reprocessing by the BLR elongates the lags most significantly. As the wind weakens, the lags are dominated by shorter accretion disk lags., Comment: 28 pages, 11 figures, 3 tables. Accepted for publication in the Astrophysical Journal

We observed the Seyfert 1 galaxy Mrk817 during an intensive multi-wavelength reverberation mapping campaign for 16 months. Here, we examine the behavior of narrow UV absorption lines seen in HST/COS spectra, both during the campaign and in other epochs extending over 14 years. We conclude that while the narrow absorption outflow system (at -3750 km/s with FWHM=177 km/s) responds to the variations of the UV continuum as modified by the X-ray obscurer, its total column density (logNH =19.5 cm-2) did not change across all epochs. The adjusted ionization parameter (scaled with respect to the variations in the Hydrogen ionizing continuum flux) is log UH =-1.0. The outflow is located at a distance smaller than 38 parsecs from the central source, which implies a hydrogen density of nH > 3000 cm-3. The absorption outflow system only covers the continuum emission source and not the broad emission line region, which suggests that its transverse size is small (< 1e16 cm), with potential cloud geometries ranging from spherical to elongated along the line of sight., Comment: 18 pages, 12 Figures, Accepted for publication in ApJ

We present the results of the XMM-Newton and NuSTAR observations taken as part of the ongoing, intensive multi-wavelength monitoring program of the Seyfert 1 galaxy Mrk 817 by the AGN Space Telescope and Optical Reverberation Mapping 2 (AGN STORM 2) Project. The campaign revealed an unexpected and transient obscuring outflow, never before seen in this source. Of our four XMM-Newton/NuSTAR epochs, one fortuitously taken during a bright X-ray state has strong narrow absorption lines in the high-resolution grating spectra. From these absorption features, we determine that the obscurer is in fact a multi-phase ionized wind with an outflow velocity of $\sim$5200 km s$^{-1}$, and for the first time find evidence for a lower ionization component with the same velocity observed in absorption features in the contemporaneous HST spectra. This indicates that the UV absorption troughs may be due to dense clumps embedded in diffuse, higher ionization gas responsible for the X-ray absorption lines of the same velocity. We observe variability in the shape of the absorption lines on timescales of hours, placing the variable component at roughly 1000 $R_g$ if attributed to transverse motion along the line of sight. This estimate aligns with independent UV measurements of the distance to the obscurer suggesting an accretion disk wind at the inner broad line region. We estimate that it takes roughly 200 days for the outflow to travel from the disk to our line of sight, consistent with the timescale of the outflow's column density variations throughout the campaign.

To investigate the short-term variability and determine the size of the optical continuum emitting size of intermediate-mass black holes (IMBHs), we carried out high-cadence, multi-band photometric monitoring of a Seyfert 1 galaxy J0249-0815 across two nights, together with a one-night single-band preliminary test. The presence of the broad Ha component in our target was confirmed by recent Palomar/P200 spectroscopic observations, 23 years after Sloan Digital Sky Survey, ruling out the supernovae origin of the broad Ha line. The photometric experiment was primarily conducted utilizing four-channel imagers MuSCAT 3 & 4 mounted on 2-meter telescopes within the Las Cumbres Observatory Global Telescope Network. Despite the expectation of variability, we observed no significant variation (<1.4%) on timescales of 6-10 hours. This non-detection is likely due to substantial host galaxy light diluting the subtle AGN variability. Dual-band preliminary tests and tailored simulations may enhance the possibility of detecting variability and lag in future IMBH reverberation campaigns., Comment: 14 pages, 6 figures, submitted to ApJ, comments welcome

We test the potential of Bayesian synthesis of upcoming multi-instrument data to extract orbital parameters and individual light curves of close binary supermassive black holes (CB-SMBH) with subparsec separations. Next generation (ng) interferometers, will make possible the observation of astrometric wobbles in CB-SMBH. Combining them with periodic variable time-domain data from surveys like the Vera C. Rubin Legacy Survey of Space and Time (LSST), allows for a more information on CB-SMBH candidates compared to standalone observational methods. Our method reliably determines binary parameters and component fluxes from binary total flux across long-term, intermediate and short-term binary dynamics and observational configurations, assuming ten annual observations, even in short period "q-accrete" objects. Expected CB-SMBH astrometric wobbles constructed from binary dynamical parameters, might serve in refining observational strategies for CB-SMBH. Combination of inferred mass ratio, light curves of binary components, and observed photocenter wobbles can be a proxy for the activity states of CB-SMBH components., Comment: The article was accepted for publication in the ApJ

We fit the UV/optical lightcurves of the Seyfert 1 galaxy Mrk 817 to produce maps of the accretion disk temperature fluctuations $\delta T$ resolved in time and radius. The $\delta T$ maps are dominated by coherent radial structures that move slowly ($v \ll c$) inwards and outwards, which conflicts with the idea that disk variability is driven only by reverberation. Instead, these slow-moving temperature fluctuations are likely due to variability intrinsic to the disk. We test how modifying the input lightcurves by smoothing and subtracting them changes the resulting $\delta T$ maps and find that most of the temperature fluctuations exist over relatively long timescales ($\sim$100s of days). We show how detrending AGN lightcurves can be used to separate the flux variations driven by the slow-moving temperature fluctuations from those driven by reverberation. We also simulate contamination of the continuum emission from the disk by continuum emission from the broad line region (BLR), which is expected to have spectral features localized in wavelength, such as the Balmer break contaminating the $U$ band. We find that a disk with a smooth temperature profile cannot produce a signal localized in wavelength and that any BLR contamination should appear as residuals in our model lightcurves. Given the observed residuals, we estimate that only $\sim$20% of the variable flux in the $U$ and $u$ lightcurves can be due to BLR contamination. Finally, we discus how these maps not only describe the data, but can make predictions about other aspects of AGN variability., Comment: 23 pages, 18 figures, submitting to ApJ, comments welcome

We explore the influence of optical Fe II quasi-continuum on the measured spectral parameters in the 4150-5500 A range for the spectra of Type 1 active galactic nuclei (AGNs). We assume that the broad line region is composed of two sub-regions: the very broad line region (VBLR) and the intermediate line region (ILR). We constructed a large set of synthetic AGN spectra by taking different portions of the VBLR and ILR contributions, where initially the VBLR and ILR model spectra were constructed on the basis of prototypes of two observed spectra with dominant VBLR (i.e. ILR) emission. To investigate the influence of the optical Fe II quasi-continuum on the AGN measured spectral parameters, we fit the power-law continuum and emission lines in a set of model spectra, as commonly done for observed AGN spectra. We then compared the spectral parameters obtained after the fitting procedure with those of the model. We find that the optical Fe II quasi-continuum can be very strong in the case of spectra with strong and very broad Fe II lines and it is difficult to fully separate it from the power-law continuum. This gives the effect of a slightly underestimated H$\beta$ width and underestimated fluxes of the H$\beta$ and Fe II lines, while the continuum flux is then slightly overestimated. The most affected spectral parameters are the line equivalent widths (EWs), especially EW Fe II, which may be strongly underestimated. We discuss the possible underlying physics in the quasar main sequence, as implied by the results of our spectral modelling. We find that the set of AGN model spectra assuming different ILR and VBLR contributions can aptly reproduce the quasar main sequence, that is, the full width at half maximum (FWHM) H$\beta$ versus Fe II/H$\beta$ anti-correlation, where both parameters in this anti-correlation are strongly dependent on the ILR and VBLR contribution rate., Comment: 12 pages, 15 figures, accepted in Astronomy & Astrophysics

An intensive reverberation mapping campaign on the Seyfert 1 galaxy Mrk817 using the Cosmic Origins Spectrograph (COS) on the Hubble Space Telescope (HST) revealed significant variations in the response of the broad UV emission lines to fluctuations in the continuum emission. The response of the prominent UV emission lines changes over a $\sim$60-day duration, resulting in distinctly different time lags in the various segments of the light curve over the 14 months observing campaign. One-dimensional echo-mapping models fit these variations if a slowly varying background is included for each emission line. These variations are more evident in the CIV light curve, which is the line least affected by intrinsic absorption in Mrk817 and least blended with neighboring emission lines. We identify five temporal windows with distinct emission line response, and measure their corresponding time delays, which range from 2 to 13 days. These temporal windows are plausibly linked to changes in the UV and X-ray obscuration occurring during these same intervals. The shortest time lags occur during periods with diminishing obscuration, whereas the longest lags occur during periods with rising obscuration. We propose that the obscuring outflow shields the ultraviolet broad lines from the ionizing continuum. The resulting change in the spectral energy distribution of the ionizing continuum, as seen by clouds at a range of distances from the nucleus, is responsible for the changes in the line response., Comment: 20 pages, 8 figures, Accepted for publication in ApJ

Development of the Rubin Observatory Legacy Survey of Space and Time (LSST) includes a series of Data Challenges (DC) arranged by various LSST Scientific Collaborations (SC) that are taking place during the projects preoperational phase. The AGN Science Collaboration Data Challenge (AGNSCDC) is a partial prototype of the expected LSST AGN data, aimed at validating machine learning approaches for AGN selection and characterization in large surveys like LSST. The AGNSC-DC took part in 2021 focusing on accuracy, robustness, and scalability. The training and the blinded datasets were constructed to mimic the future LSST release catalogs using the data from the Sloan Digital Sky Survey Stripe 82 region and the XMM-Newton Large Scale Structure Survey region. Data features were divided into astrometry, photometry, color, morphology, redshift and class label with the addition of variability features and images. We present the results of four DC submitted solutions using both classical and machine learning methods. We systematically test the performance of supervised (support vector machine, random forest, extreme gradient boosting, artificial neural network, convolutional neural network) and unsupervised (deep embedding clustering) models when applied to the problem of classifying/clustering sources as stars, galaxies or AGNs. We obtained classification accuracy 97.5% for supervised and clustering accuracy 96.0% for unsupervised models and 95.0% with a classic approach for a blinded dataset. We find that variability features significantly improve the accuracy of the trained models and correlation analysis among different bands enables a fast and inexpensive first order selection of quasar candidates, Comment: Accepted by ApJ. 21 pages, 14 figures, 5 tables

The AGN STORM 2 campaign is a large, multiwavelength reverberation mapping project designed to trace out the structure of Mrk 817 from the inner accretion disk to the broad emission line region and out to the dusty torus. As part of this campaign, Swift performed daily monitoring of Mrk 817 for approximately 15 months, obtaining observations in X-rays and six UV/optical filters. The X-ray monitoring shows that Mrk 817 was in a significantly fainter state than in previous observations, with only a brief flare where it reached prior flux levels. The X-ray spectrum is heavily obscured. The UV/optical light curves show significant variability throughout the campaign and are well correlated with one another, but uncorrelated with the X-rays. Combining the Swift UV/optical light curves with Hubble UV continuum light curves, we measure interband continuum lags, $\tau(\lambda)$, that increase with increasing wavelength roughly following $\tau(\lambda) \propto \lambda^{4/3}$, the dependence expected for a geometrically thin, optically thick, centrally illuminated disk. Modeling of the light curves reveals a period at the beginning of the campaign where the response of the continuum is suppressed compared to later in the light curve - the light curves are not simple shifted and scaled versions of each other. The interval of suppressed response corresponds to a period of high UV line and X-ray absorption, and reduced emission line variability amplitudes. We suggest that this indicates a significant contribution to the continuum from the broad line region gas that sees an absorbed ionizing continuum., Comment: 20 pages, 13 figures, 3 tables, accepted for publication in ApJ

In this study, a refined approach for multicomponent fitting of active galactic nuclei (AGN) spectra is presented utilizing the newly developed Python code $fantasy$ (fully automated python tool for AGN spectra analysis). AGN spectra are modeled by simultaneously considering the underlying broken power-law continuum, predefined emission line lists, and an Fe II model, which is here extended to cover the wavelength range 3700 - 11000 A. The Fe II model, founded solely on atomic data, effectively describes the extensive emission of the complex iron ion in the vicinity of the H$\gamma$ and H$\beta$ lines, as well as near the H$\alpha$ line, which was previously rarely studied. The proposed spectral fitting approach is tested on a sample of high-quality AGN spectra from the Sloan Digital Sky Survey (SDSS) Data Release 17. The results indicate that when Fe II emission is present near H$\beta$, it is also detected redward from H$\alpha$, potentially contaminating the broad H$\alpha$ line wings and thus affecting the measurements of its flux and width. The production of Fe II emission is found to be strongly correlated with Eddington luminosity and appears to be controlled by the similar mechanism as the hydrogen Balmer lines. The study highlights the benefits of fitting AGN type 1 spectra with the $fantasy$ code, pointing that it may be used as a robust tool for analyzing a large number of AGN spectra in the coming spectral surveys., Comment: Accepted in ApJS

Deep learning techniques are required for the analysis of synoptic (multi-band and multi-epoch) light curves in massive data of quasars, as expected from the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). In this follow-up study, we introduced an upgraded version of a conditional neural process (CNP) embedded in a multistep approach for analysis of large data of quasars in the LSST Active Galactic Nuclei Scientific Collaboration data challenge database. We present a case study of a stratified set of the u-band light curves for 283 quasars with very low variability $\sim 0.03$. In this sample, CNP average mean square error is found to be $\sim 5\% $($\sim 0.5$ mag). Interestingly, beside similar level of variability there are indications that individual light curves show flare like features. According to preliminary structure function analysis, these occurrences may be associated to microlensing events with larger time scales $5-10$ years., Comment: accepted for publication by the MDPI Universe Special Issue A Multimessenger View of Supermassive Black Holes and the Quasar Main Sequence, Eds. Ascension Del Olmo and Paola Marziani

We analyze the broad H$\beta$ line profile variability of the "changing look" active galactic nucleus (CL-AGN) NGC 3516 over a long period of 25 years. The observed change in the broad line profile may indicate a change in the geometry of the broad line region (BLR). Using spectral line profiles, we aim to explore changes in the kinematics and dimensions of the BLR in NGC 3516. We consider two possible scenarios, i.e. changes in the broad-line emission are caused by a decrease of ionization continuum emission or by the BLR obscuration by outer dusty regions. With this investigation we aim to clarify the CL mechanism of this AGN. We analyze the spectral band around the H$\beta$ line as well as the broad H$\beta$ line parameters, and how they change in time. We model the broad-line profiles assuming that there is an emission from the accretion disc superposed with an emission from a surrounding region that is outside the disc. We find that in the Type 1 activity phase, the BLR is very complex. There is a clear disc-like BLR that contributes to the broad line wings and an additional intermediate line region (ILR) that contributes to the line core. In the high activity phase, the ILR emission is close to the center of the line (in some cases slightly shifted to the red), whereas in the low activity phase (i.e., Type 2 phase), the ILR component has a significant shift to the blue, indicating an outflow. We propose that the changing look mechanism in NGC 3516 is rather connected with the intrinsic effects than with an outer obscuring region. It may still be possible that the dust has an important role in the low activity phase when it is coming inside of the BLR, making a dusty BLR. In this way, it causes a decrease in the ionization and recombination rates., Comment: 16 pages, accepted in A&A

We report the results of a search for long-period ($100

The AGN STORM 2 collaboration targeted the Seyfert 1 galaxy Mrk 817 for a year-long multiwavelength, coordinated reverberation mapping campaign including HST, Swift, XMM-Newton, NICER, and ground-based observatories. Early observations with NICER and XMM revealed an X-ray state ten times fainter than historical observations, consistent with the presence of a new dust-free, ionized obscurer. The following analysis of NICER spectra attributes variability in the observed X-ray flux to changes in both the column density of the obscurer by at least one order of magnitude ($N_\mathrm{H}$ ranges from $2.85\substack{+0.48\\ -0.33} \times 10^{22}\text{ cm}^{-2}$ to $25.6\substack{+3.0\\ -3.5} \times 10^{22} \text{ cm}^{-2}$) and the intrinsic continuum brightness (the unobscured flux ranges from $10^{-11.8}$ to $10^{-10.5}$ erg s$^{-1}$ cm$^{-2}$ ). While the X-ray flux generally remains in a faint state, there is one large flare during which Mrk 817 returns to its historical mean flux. The obscuring gas is still present at lower column density during the flare but it also becomes highly ionized, increasing its transparency. Correlation between the column density of the X-ray obscurer and the strength of UV broad absorption lines suggests that the X-ray and UV continua are both affected by the same obscuration, consistent with a clumpy disk wind launched from the inner broad line region., Comment: 19 pages, 9 figures, Accepted to ApJ

We present reverberation mapping measurements for the prominent ultraviolet broad emission lines of the active galactic nucleus Mrk817 using 165 spectra obtained with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Our ultraviolet observations are accompanied by X-ray, optical, and near-infrared observations as part of the AGN Space Telescope and Optical Reverberation Mapping Program 2 (AGN STORM 2). Using the cross-correlation lag analysis method, we find significant correlated variations in the continuum and emission-line light curves. We measure rest-frame delayed responses between the far-ultraviolet continuum at 1180 A and Ly$\alpha$ $\lambda1215$ A ($10.4_{-1.4}^{+1.6}$ days), N V $\lambda1240$ A ($15.5_{-4.8}^{+1.0}$days), SiIV + OIV] $\lambda1397$ A ($8.2_{-1.4}^{+1.4}$ days), CIV $\lambda1549$ A ($11.8_{-2.8}^{+3.0}$ days), and HeII $\lambda1640$ A ($9.0_{-1.9}^{+4.5}$ days) using segments of the emission-line profile that are unaffected by absorption and blending, which results in sampling different velocity ranges for each line. However, we find that the emission-line responses to continuum variations are more complex than a simple smoothed, shifted, and scaled version of the continuum light curve. We also measure velocity-resolved lags for the Ly$\alpha$, and CIV emission lines. The lag profile in the blue wing of Ly$\alpha$ is consistent with virial motion, with longer lags dominating at lower velocities, and shorter lags at higher velocities. The CIV lag profile shows the signature of a thick rotating disk, with the shortest lags in the wings, local peaks at $\pm$ 1500 $\rm km\,s^{-1}$, and a local minimum at line center. The other emission lines are dominated by broad absorption lines and blending with adjacent emission lines. These require detailed models, and will be presented in future work., Comment: Submitted to ApJ. 25 pages, 8 figures, and 6 tables

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., Comment: Astronomy & Astrophysics in press

Earlier, we suggested the "reload" concept of the polarimetric reverberation mapping of active galactic nuclei (AGN), proposed for the first time more than 10 years ago. We have successfully tested this approach of reverberation mapping of the broad emission line on the galaxy Mrk 6. It was shown that such an idea allows one to look at the AGN central parsec structure literally in a new light. However, the method originally assumed the use of spectropolarimetric observations, expensive in terms of telescope time, and implemented on rare large telescopes. Currently, we propose an adaptation of the polarimetric reverberation mapping of broad lines in medium-band filters following the idea of the photometric reverberation mapping, when filters are selected so that their bandwidth is oriented to the broad line and the surrounding continuum near. In this paper, we present the progress status of such monitoring conducted jointly at the Special astrophysical observatory and Asiago Cima Ekar observatory (OAPd/INAF) with support from Rozhen National Astronomical Observatory (NAO), some first results for the most frequently observed AGNs Mrk 335, Mrk 509, and Mrk 817, and the discussion of the future perspectives of the campaign., Comment: Published in Universe

Extraction of information in the form of oscillations from noisy data of natural phenomena such as sounds, earthquakes, ionospheric and brain activity, and various emissions from cosmic objects is extremely difficult. As a method for finding periodicity in such challenging data sets, the 2D Hybrid approach, which employs wavelets, is presented. Our technique produces a wavelet transform correlation intensity contour map for two (or one) time series on a period plane defined by two independent period axes. Notably, by spreading peaks across the second dimension, our method improves apparent resolution of detected oscillations in the period plane and identifies the direction of signal changes using correlation coefficients. We demonstrate the performance of the 2D Hybrid technique on a very low frequency (VLF) signal emitted in Italy and recorded in Serbia in time vicinity of the occurrence of an earthquake on November 3, 2010, near Kraljevo, Serbia. We identified a distinct signal in the range 120-130 s that appears only in association with the considered earthquake. Other wavelets, such as Superlets, which may detect fast transient oscillations, will be employed in the future analysis., Comment: published in Mathematics MDPI

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 01000 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., Comment: accepted for publication in the Astrophysical Journal Focus Issue on Rubin LSST cadence and survey strategy

The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) dataset will dramatically alter our understanding of the Universe, from the origins of the Solar System to the nature of dark matter and dark energy. Much of this research will depend on the existence of robust, tested, and scalable algorithms, software, and services. Identifying and developing such tools ahead of time has the potential to significantly accelerate the delivery of early science from LSST. Developing these collaboratively, and making them broadly available, can enable more inclusive and equitable collaboration on LSST science. To facilitate such opportunities, a community workshop entitled "From Data to Software to Science with the Rubin Observatory LSST" was organized by the LSST Interdisciplinary Network for Collaboration and Computing (LINCC) and partners, and held at the Flatiron Institute in New York, March 28-30th 2022. The workshop included over 50 in-person attendees invited from over 300 applications. It identified seven key software areas of need: (i) scalable cross-matching and distributed joining of catalogs, (ii) robust photometric redshift determination, (iii) software for determination of selection functions, (iv) frameworks for scalable time-series analyses, (v) services for image access and reprocessing at scale, (vi) object image access (cutouts) and analysis at scale, and (vii) scalable job execution systems. This white paper summarizes the discussions of this workshop. It considers the motivating science use cases, identified cross-cutting algorithms, software, and services, their high-level technical specifications, and the principles of inclusive collaborations needed to develop them. We provide it as a useful roadmap of needs, as well as to spur action and collaboration between groups and individuals looking to develop reusable software for early LSST science., Comment: White paper from "From Data to Software to Science with the Rubin Observatory LSST" workshop

Recent studies have proposed that GRAVITY+ instrument is able to trace the circular orbit of the subparsec close-binary supermassive black holes (CB-SMBHs) by measuring the photocentre variation of the hot dust emission. However, the CB-SMBHs orbit may become highly eccentric throughout the evolution of these objects, and the orbital period may be far longer than the observational time baseline. We investigate the problem of detecting the CB-SMBH with hot dust emission and high eccentricity (eCBSMBH, e=0.5) when the observed time baselines of their astrometric data and radial velocities are considerably shorter than the orbital period. The parameter space of the Keplerian model of the eCBSMBH is large for exploratory purposes. We therefore applied the Bayesian method to fit orbital elements of the eCBSMBH to combined radial velocity and astrometric data covering a small fraction of the orbital period. We estimate that a number of potential eCBSMBH systems within reach of GRAVITY + will be similar to the number of the planned circular targets. We show that using observational time baselines that cover ~ 10% of the orbit increases the possibility of determining the period, eccentricity, and total mass of an eCBSMBH. When the observational time baseline becomes too short (~ 5%), the quality of the retrieved eCBSMBH parameters degrades. We also illustrate how interferometry may be used to estimate the photo-centre at the eCBSMBH emission line, which could be relevant for GRAVITY+ successors. Even if the astrometric signal for eCBSMBH systems is reduced by a factor of sqrt{1-e^{2}} compared to circular ones, we find that the hot dust emission of eCBSMBHs can be traced by GRAVITY+ at the elementary level., Comment: accepted for publication in the journal Astronomy & Astrophysics

In spectra of the Active Galactic Nuclei (AGNs), the [N II] 6548, 6583 A lines are commonly fitted using the fixed intensity ratio of these two lines (R[N II]=I$_{6583}$/I$_{6548}$). However, the used values for fixed intensity ratio are slightly different through literature. There are several theoretical calculations of the transition probabilities which can be used for the line ratio estimation, but there are no experimental measurements of this ratio, since the [N II] lines are extremely weak in laboratory plasma. Therefore, the intensity ratio of [N II] lines can be measured only in the spectra of astrophysical objects. However, precise and systematic measurements have not be done so far, because of difficulties in measurement of the [N II] ratio in various spectra (overlapping with H$\alpha$, weak intensity of [N II], influence of the continuum noise and outflow contribution, etc.). Here we present the measurements of the flux ratio of the [N II]$\lambda\lambda$ 6548, 6583 A emission lines for a sample of 250 Type 2 AGNs spectra taken form Sloan Digital Sky Survey (SDSS) data base. The spectra are chosen to have high signal-to-noise ratio and to [N II] and H$\alpha$ lines do not overlap. The obtained mean flux ratio from measurements is 3.049 $\pm$ 0.021. Our result is in agreement with theoretical result obtained by taking into account the relativistic corrections to the magnetic dipole operator., Comment: Accepted in Advances in Space Research, 8 pages, 4 figures, 1 table

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$\alpha$ 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., Comment: 25 pages, 21 figures, 11 tables, to be published in Astronomy & Astrophysics, version after proofs

We have used the sample of 577 active galactic nuclei Type 1.8-2 spectra (z < 0.25), taken from Sloan Digital Sky Survey, to trace the influence of the outflow kinematics to the profiles of different emission lines (Hbeta, [O III], Halpha, [N II], [S II]). All considered lines were fitted with two Gaussian components: one which fits the core of the line, and another which fits the wings. We gave the procedure for decomposition of Halpha+[N II] wavelength band, for the spectra where these lines overlap. The influence of the gravitational/non-gravitational kinematics to the line components is investigated by comparing the dispersions of the line components with stellar velocity dispersion. We found that wing components of all considered emission lines have pure non-gravitational kinematics, the core components are consistent with gravitational kinematics for the Halpha, [N II] and [S II] lines, while in the [O III] there is evidence for contribution from non-gravitational kinematics. We adopted the wing components as proxy of the outflow contribution and we investigated the outflow kinematics by analysing the correlations between widths and between shifts of the wing components of different lines. We found the strong correlations between shifts and between wing component widths of all considered lines, with exception of the Hbeta wing component width. These correlations indicate that outflow dynamics systemically affects all emission lines in spectrum. However, it reflects with different strength in their profiles, which is observed as different widths of the wing components. The strongest outflow signature is observed in the [O III] lines, which have the broadest wing components, weaker in Halpha and [N II], and the weakest in [S II]. These results imply that considered lines arise in different parts of an outflowing region., Comment: Accepted for publication in Astronomy and Astrophysics, 18 pages, 15 figures, 2 tables

The upcoming photometric surveys, such as the Rubin Observatory's Legacy Survey of Space and Time (LSST) will monitor unprecedented number of active galactic nuclei (AGN) in a decade long campaign. Motivated by the science goals of LSST, which includes the harnessing of broadband light curves of AGN for photometric reverberation mapping (PhotoRM), we implement the existing formalism to estimate the lagged response of the emission line flux to the continuum variability using only mutli-band photometric light curves. We test the PhotoRM method on a set of 19 artificial light curves simulated using a stochastic model based on the Damped Random Walk process. These light curves are sampled using different observing strategies, including the two proposed by the LSST, in order to compare the accuracy of time-lag retrieval based on different observing cadences. Additionally, we apply the same procedure for time-lag retrieval to the observed photometric light curves of NGC 4395, and compare our results to the existing literature., Comment: The manuscript was accepted for publication in the Astronomische Nachrichten Special Issue

We investigate the broad line profiles emitted from a system supermassive binary black hole (SMBBH) having elliptical orbits and low mass ratio of $m_2/m_1\sim 0.1$. Our model assumes a super Eddington accretion flow in the case of a smaller component, whereas the massive component has very small or negligible accretion, therefore supposing that no broad line region (BLR) is attached to it. Thus, the proposed SMBBH system contains one moving BLR, associated with the less massive component and one circum-binary BLR. We study the effect of different total mass of the system (ranging from 10$^6$ to 10$^8$ Solar masses) to the $\mathrm{H\beta}$ line profiles and to the continuum and line light curves. The resulted broad line profiles are asymmetric and shifted, and are varying during the orbital period. The asymmetry in the broad line profiles is discussed in terms of expected differences between the proposed model of the SMBBH with one active component and the scenario of a recoiling black hole. We discuss the periodicity detected in the line and continuum light curves, as well as in the variations of the line asymmetry and shift., Comment: 10 pages, 8 figures, 1 table

The consequences of complex disturbed environments in the vicinity of a supermassive black hole are not well represented by standard statistical models of optical variability in active galactic nuclei (AGN). Thus, developing new methodologies for investigating and modeling AGN light curves is crucial. Conditional Neural Processes (CNPs) are nonlinear function models that forecast stochastic time-series based on a finite amount of known data without the use of any additional parameters or prior knowledge (kernels). We provide a CNP algorithm that is specifically designed for simulating AGN light curves. It was trained using data from the All-Sky Automated Survey for Supernovae, which included 153 AGN. We present CNP modeling performance for a subsample of five AGNs with distinctive difficult-to-model properties. The performance of CNP in predicting temporal flux fluctuation was assessed using a minimizing loss function, and the results demonstrated the algorithm's usefulness. Our preliminary parallelization experiments show that CNP can efficiently handle large amounts of data. These results imply that CNP can be more effective than standard tools in modeling large volumes of AGN data (as anticipated from time-domain surveys such as the Vera C. Rubin Observatory's Legacy Survey of Space and Time)., Comment: The manuscript was accepted for publication in the Astronomische Nachrichten Special Issue

The differences between Narrow Line Seyfert 1 galaxies (NLS1s) and Broad Line AGNs (BLAGNs) are not completely understood; it is thought that they may have different inclinations and/or physical characteristics. The FWHM(Hb)-luminosities correlations are found for NLS1s and their origin is the matter of debate. Here we investigated the spectroscopic parameters and their correlations considering a dusty, cone model of AGN. We apply a simple conical dust distribution (spreading out of broad line region, BLR), assuming that the observed surface of the model is in a good correlation with MIR emission. The dusty cone model in combination with a BLR provides the possibility to estimate luminosity dependence on the cone inclination. The FWHM(Hb)-luminosities correlations obtained from model in comparison with observational data show similarities which may indicate the influence of AGN inclination and structure to this correlation. An alternative explanation for FWHM(Hb)-luminosities correlations is the selection effect by the black hole mass. These FWHM(Hb)-luminosities correlations may be related to the starburst in AGNs, as well. The distinction between spectral properties of the NLS1s and BLAGNs could be caused by multiple effects: beside physical differencies between NLS1s and BLAGNs (NLS1s have lighter black hole mass than BLAGNs), inclination of the conical AGN geometry may have important role as well, where NLS1s may be seen in lower inclination angles., Comment: Accepted to MNRAS

Spectropolarimetric efforts in the last few years have provided an efficient method that is based on the profiles of the polarization plane position angle of broad emission lines in active galactic nuclei (AGNs). Here we present black hole measurements of SBS 1419+538 using spectropolarimetric observations in the Mg II spectral band. The observations are performed by 6m telescope of SAO RAS using SCORPIO-2. We found a good agreement for the estimated supermassive black hole (SMBH) mass for this object using spectropolarimetry when compared with the mass obtained using other methods., Comment: Accepted by ApJL on the 24th of September 2021, 7 pages, 3 figures

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$\beta$ and H$\alpha$ 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., Comment: 22 pages, 3 tables and 14 figures. Accepted for publication in MNRAS

Motivated by upcoming photometric and spectroscopic surveys (Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), Manuakea Spectroscopic Explorer), we design the statistical proxies to measure the cadence effects on active galactic nuclei (AGN) variability-observables (time-lags, periodicity, and structure-function (SF)). We constructed a multiple-regression model to statistically identify the cadence-formal error pattern knowing AGN time-lags and periodicity from different surveys. We defined the simple metric for the SF's properties, accounting for the 'observed' SF's deviation relative to those obtained from the homogenously-sampled light curves. We tested the regression models on different observing strategies: the optical dataset of long light-curves of eight AGN with peculiarities and the artificial datasets based on several idealized and LSST-like cadences. The SFs metric is assessed on synthetic datasets. The regression models (for both data types) predict similar cadences for time-lags and oscillation detection, whereas for light curves with low variability ($\sim 10\%$), cadences for oscillation detection differ. For higher variability ($\sim20\%$), predicted cadences are larger than for $F_{var}\sim 10\%$. The predicted cadences are decreasing with redshift. SFs with dense and homogenous cadences are more likely to behave similarly. SFs with oscillatory signals are sensitive to the cadences, possibly impacting LSST-like operation strategy. The proposed proxies can help to select spectroscopic and photometric-surveys cadence strategies, and they will be tested further in larger samples of objects., Comment: accepted for publication in the MNRAS

We have developed two metrics related to AGN variability observables (time-lags, periodicity, and Structure Function (SF)) to evaluate LSST OpSim FBS 1.5, 1.6, 1.7 performance in AGN time-domain analysis. For this purpose, we generate an ensemble of AGN light curves based on AGN empirical relations and LSST OpSim cadences. Although our metrics show that denser LSST cadences produce more reliable time-lag, periodicity, and SF measurements, the discrepancies in the performance between different LSST OpSim cadences are not drastic based on Kullback-Leibler divergence. This is complementary to Yu and Richards results on DCR and SF metrics, extending them to include the point of view of AGN variability., Comment: Cadence Notes AGN and TVS SC

Here we consider the perspective to detect sub-pc super-massive binary black-hole (SMBBH) systems using long-term photometric and spectroscopic monitoring campaigns of active galactic nuclei. This work explores the nature of long-term spectral variability caused by the dynamical effects of SMBBH systems. We describe in great detail a model of SMBBH system which considers that both black holes have their accretion disc and additional line emitting region(s). We simulate the H$\beta$ spectral band (continuum+broad H$\beta$ line) for different mass ratios of components and different total masses of the SMBBH systems ($10^6-10^8\mathrm{M\odot}$). We analyze the set of continuum and broad line light curves for several full orbits of SMBBHs with different parameters, to test the possibility to extract the periodicity of the system. We consider different levels of the signal-to-noise ratio, which is added to the simulated spectra. Our analysis showed that the continuum and broad line profiles emitted from an SMBBH system are strongly dependent, not only on the mass ratio of the components but also on the total mass of the system. We found that the mean broad line profile and its rms could indicate the presence of an SMBBH. However, some effects caused by the dynamics of a binary system could be hidden due to a low signal-to-noise ratio. Finally, we can conclude that the long-term AGN monitoring campaigns could be beneficial for the detection of SMBBH candidates., Comment: 20 pages, 22 figures

We present the first results from the ongoing, intensive, multi-wavelength monitoring program of the luminous Seyfert 1 galaxy Mrk 817. While this AGN was, in part, selected for its historically unobscured nature, we discovered that the X-ray spectrum is highly absorbed, and there are new blueshifted, broad and narrow UV absorption lines, which suggest that a dust-free, ionized obscurer located at the inner broad line region partially covers the central source. Despite the obscuration, we measure UV and optical continuum reverberation lags consistent with a centrally illuminated Shakura-Sunyaev thin accretion disk, and measure reverberation lags associated with the optical broad line region, as expected. However, in the first 55 days of the campaign, when the obscuration was becoming most extreme, we observe a de-coupling of the UV continuum and the UV broad emission line variability. The correlation recovers in the next 42 days of the campaign, as Mrk 817 enters a less obscured state. The short CIV and Ly alpha lags suggest that the accretion disk extends beyond the UV broad line region., Comment: 28 pages, 14 figures, submitted to ApJ. Comments welcome

On 2019/07/30.86853 UT, IceCube detected a high-energy astrophysical neutrino candidate. The Flat Spectrum Radio Quasar PKS 1502+106 is located within the 50 percent uncertainty region of the event. Our analysis of 15 GHz Very Long Baseline Array (VLBA) and astrometric 8 GHz VLBA data, in a time span prior and after the IceCube event, reveals evidence for a radio ring structure which develops with time. Several arc-structures evolve perpendicular to the jet ridge line. We find evidence for precession of a curved jet based on kinematic modelling and a periodicity analysis. An outflowing broad line region (BLR) based on the C IV line emission (Sloan Digital Sky Survey, SDSS) is found. We attribute the atypical ring to an interaction of the precessing jet with the outflowing material. We discuss our findings in the context of a spine-sheath scenario where the ring reveals the sheath and its interaction with the surroundings (narrow line region, NLR, clouds). We find that the radio emission is correlated with the $\gamma$-ray emission, with radio lagging the $\gamma$-rays. Based on the $\gamma$-ray variability timescale, we constrain the $\gamma$-ray emission zone to the BLR (30-200 $r_{\rm g}$) and within the jet launching region. We discuss that the outflowing BLR provides the external radiation field for $\gamma$-ray production via external Compton scattering. The neutrino is most likely produced by proton-proton interaction in the blazar zone (beyond the BLR), enabled by episodic encounters of the jet with dense clouds, i.e. some molecular cloud in the NLR., Comment: 35 pages, 33 figures, 3 tables; accepted by the MNRAS Main Journal

The Very Large Telescope Interferometer (VLTI), and the Extremely Large Telescope (ELT) will be a robust astrophysics suite offering the opportunity of probing the structure and dynamics of CB-SMBH at high spectral and angular resolution. Here, we explore and illustrate the application of differential interferometry on unresolved the CB-SMBH systems in elliptical orbital configurations and a single SMBH with clouds in elliptical orbital motion. Photocenter displacements between each SMBH and regions in their disc-like broad line regions (BLR) appear as small interferometric differential phase variability. To investigate the application of interferometric phases for the detection of CB-SMBH systems, we simulate a series of differential interferometry signatures, based on our model comprising ensembles of clouds surrounding each of supermassive black hole in a CB-SMBH. Setting model to the parameters of a single SMBH with elliptical cloud motion, we also calculated a series of differential interferometry observables for this case. We found various deviations from the canonical S-shaped of CB-SMBH phases profile for elliptically configured CB-SMBH systems. The amplitude and specific shape of the interferometry observables depend on orbital configurations of the CB-SMBH system. We get distinctive results when considering antialigned angular momenta of cloud orbits regarding total CB-SMBH angular momentum. Some simulated spectral lines from our model closely resemble observations of Pa \alpha line got from near-infrared AGN surveys. We found differences between differential phases zoo of a single SMBH and CB-SMBH systems. The differential phases zoo for a single SMBH comprises deformed S shape. We also showed how their differential phase shape, amplitude, and slope evolve with various sets of cloud orbital parameters and observer position., Comment: accepted by A&A Journal

For type-1 active galactic nuclei (AGNs) for which the equatorial scattering is the dominant broad line polarization mechanism, it is possible to measure the supermassive black hole mass by tracing the Keplerian motion across the polarization plane positionangle. So far this method has been used for 30 objects but only for H${\alpha}$ emission line. We explore the possibilities this method for determining SMBH masses using polarization in broad emission lines by applying it for the first time to Mg II${\lambda}$2798AA spectral line. We use 3D Monte Carlo radiative transfer code stokes for simultaneous modeling of equatorial scattering of H${\alpha}$, H${\beta}$ and Mg II lines. We included vertical inflows and outflows in the MgII broad line region (BLR). We find that polarization states of H${\alpha}$ and H${\beta}$ lines are almost identical and SMBH mass estimates differ by 7%. For Mg II line, we find that polarization angle ${\varphi}$ exhibits an additional "plateau" with a constant ${\varphi}$, which deviates than the profiles expected for pure Keplerian motion. SMBH mass estimates using Mg II line are higher by up to 35% than those obtained from H${\alpha}$ and H${\beta}$ lines. Our model shows that for vertical inflows and outflows in the BLR that is higher or comparable to Keplerian velocity, this method can be applied as a first approximation for obtaining SMBH mass., Comment: 8 pages, 4 figures, 3 tables, accepted for publishing in MNRAS

The active galactic nuclei (AGN) are among the most powerful sources with an inherent, pronounced and random variation of brightness. The randomness of their time series is so subtle as to blur the border between aperiodic fluctuations and noisy oscillations. This poses challenges to analysing of such time series because neither visual inspection nor pre-exisitng methods can identify well oscillatory signals in them. Thus, there is a need for an objective method for periodicity detection. Here we review our a new data analysis method that combines a two-dimensional correlation (2D) of time series with the powerful methods of Gaussian processes. To demonstrate the utility of this technique, we apply it to two example problems which were not exploited enough: damped rednoised artificial time series mimicking AGN time series and newly published observed time series of changing look AGN (CL AGN) NGC 3516. The method successfully detected periodicities in both types of time series. Identified periodicity of $\sim 4$ yr in NGC 3516 allows us to speculate that if the thermal instability formed in its accretion disc (AD) on a time scale resembling detected periodicity then AD radius could be $\sim 0.0024$ pc., Comment: accepted by Open astronomy

Here we confirm the short periodic variability of a subparsec supermassive binary black hole (SMBBH) candidate Mrk 231 in the extended optical photometric data set collected by the Catalina Real-Time Transient Survey (CRTS) and All-Sky Automated Survey for Supernovae (ASAS-SN). Using the Lomb-Scargle periodogram and 2DHybrid method, we detected the significant periodicity of ~ 1.1 yr beyond a damped random walk model in the CRTS+ASAS-SN optical data set. Mrk 231 has been previously proposed as a SMBBH candidate with a highly unequal mass ratio (q~ 0.03), very tight mutual separation of ~590 AU, and an orbital period of ~1.2 yr. Hence, our result further supports, even though not prove, the intriguing hypothesis that SMBBHs with low mass ratios may be more common than close-equal mass SMBBHs. This result, however, was obtained from the contribution of CRTS data with limited sampling cadence and photometric accuracy, and further monitoring of Mrk 231 is crucial to confirm the periodicity., Comment: Accepted for publication in MNRAS

We present the variability study of the lowest-luminosity Seyfert 1 galaxy NGC 4395 based on the photometric monitoring campaigns in 2017 and 2018. Using 22 ground-based and space telescopes, we monitored NGC 4395 with a $\sim$5 minute cadence during a period of 10 days and obtained light curves in the UV, V, J, H, and K/Ks bands as well as the H$\alpha$ narrow-band. The RMS variability is $\sim$0.13 mag on \emph{Swift}-UVM2 and V filter light curves, decreasing down to $\sim$0.01 mag on K filter. After correcting for continuum contribution to the H$\alpha$ narrow-band, we measured the time lag of the H$\alpha$ emission line with respect to the V-band continuum as ${55}^{+27}_{-31}$ to ${122}^{+33}_{-67}$ min. in 2017 and ${49}^{+15}_{-14}$ to ${83}^{+13}_{-14}$ min. in 2018, depending on the assumption on the continuum variability amplitude in the H$\alpha$ narrow-band. We obtained no reliable measurements for the continuum-to-continuum lag between UV and V bands and among near-IR bands, due to the large flux uncertainty of UV observations and the limited time baseline. We determined the AGN monochromatic luminosity at 5100\AA\ $\lambda L_\lambda = \left(5.75\pm0.40\right)\times 10^{39}\,\mathrm{erg\,s^{-1}}$, after subtracting the contribution of the nuclear star cluster. While the optical luminosity of NGC 4395 is two orders of magnitude lower than that of other reverberation-mapped AGNs, NGC 4395 follows the size-luminosity relation, albeit with an offset of 0.48 dex ($\geq$2.5$\sigma$) from the previous best-fit relation of Bentz et al. (2013)., Comment: Accepted for publication in ApJ (Feb. 23rd, 2020). 18 pages, 10 figures

The mass measurement of supermassive black holes (SMBHs) is a very complex task. Between several methods for SMBH mass measurements, some of them use the spectral lines, which indicate the motion of the emitting/absorbing material around an SMBH. Mostly there is assumption of virialization of line emitting gas in the region close to SMBHs that is used for their mass measurements. In this paper we will give an overview of methods for the SMBH mass measurements using broad emission spectral lines observed in Type 1 AGNs. First we give the basic idea to use the parameters of broad lines to SMBH mass measurements. After that we give an overview of broad lines from X-ray (Fe k{\alpha}) to the IR (Pashen and Brecket lines) which have been used for SMBH mass estimates. Additionally, we describe and discuss a new method for SMBH mass measurements using the polarization in the broad lines emitted from Type 1 AGNs., Comment: The review was presented at 12th Serbian Confer- ence on Spectral Line Shapes in Astrophysics in the special session "Broad lines in AGNs: The physics of emission gas in the vicinity of super-massive black hole" (In memory of the life and work of dr Alla Ivanovna Shapovalova). Revised version is accepted to be published in Open Astronomy

An intense search for kinematic signatures of reverberation mapping of close binary SMBH (SMBBH) is currently ongoing. Elliptical SMBBH systems (i.e. both orbits and disc-like broad-line regions (BLR) are elliptic) have not been assessed in two-dimensional transfer functions (2DTF) studies. We aim to numerically infer such a 2DTF because the geometry of the unresolved region is imprinted on their optical emission. Through this, we determine their specific kinematical signatures. We simulated the geometry and kinematics of SMBBH whose components are on elliptical orbits. Each SMBH had a disc-like elliptical BLR. The SMBHs were active and orbited each other tightly at a subparsec distance. Here we calculate for the first time 2DTF, as defined in the velocity-time delay plane, for several elliptical configurations of SMBBH orbits and their BLRs. We find that these very complex configurations are clearly resolved in maps. These results are distinct from those obtained from circular and disc-wind geometry. We calculate the expected line variability for all SMBBH configurations. We show that the line shapes are influenced by the orbital phase of the SMBBH. Some line profiles resemble observed profiles, but they can also be much deformed to look like those from the disc-wind model. First, our results imply that using our 2DTF, we can detect and quantify kinematic signatures of elliptical SMBBH. Second, the calculated expected line profiles share some intriguing similarities with observed profiles, but also with some profiles that are synthesised in disc-wind models. To overcome the non-uniqueness of the spectral line shapes as markers of SMBBH, they must be accompanied with 2DTF., Comment: accepted for publication to A&A

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more., Comment: 9 chapters, 301 pages, 100 figures. This version of the DSC is a comprehensive update of the original version, released in 2016, which can be downloaded at arXiv:1606.00043. A detailed summary of the design of MSE is available in the MSE Book 2018, available at arXiv:1810.08695

The ubiquitous variability of quasars across a wide range of wavelengths and timescales encodes critical information about the structure and dynamics of the circumnuclear emitting regions that are too small to be directly resolved, as well as detailed underlying physics of accretion and feedback processes in these active supermassive black holes. We emphasize the importance of studying quasar variability with time-domain spectroscopy, focusing on two science cases: (1) reverberation mapping (RM) to measure the broad-line region sizes and black hole masses in distant quasars; (2) spectroscopic follow-up of extreme variability quasars that dramatically change their continuum and broad-line flux within several years. We highlight the need for dedicated optical-infrared spectroscopic survey facilities in the coming decades to accompany wide-area time-domain imaging surveys, including: (1) the next phase of the Sloan Digital Sky Survey (SDSS-V; ~2020-2025), an all-sky, time-domain multi-object spectroscopic survey with 2.5m-class telescopes; (2) the planned Maunakea Spectroscopic Explorer, a dedicated 10m-class spectroscopic survey telescope with a 1.5 sq. deg field-of-view and multiplex of thousands of fibers in both optical and near-IR (J+H) to begin operations in 2029; (3) the Time-domain Spectroscopic Observatory (TSO), a proposed Probe-class ~1.3m telescope at L2, with imaging and spectroscopy (R=200, 1800) in 4 bands (0.3 - 5 micron) and rapid slew capability to 90% of sky, which will extend the coverage of Hbeta to z=8., Comment: Science White Paper for the US Astro 2020 Decadal Survey

We investigate the structure of the Mg II broad line emission region in a sample of 284 Type 1 active galactic nuclei (AGNs), through comparing the kinematical parameters of the broad Mg II and broad H$\beta$ lines. We found that the Mg II emitting region has more complex kinematics than the H$\beta$ one. It seems that the \ion{Mg}{ii} broad line originates from two subregions: one which contributes to the line core, which is probably virialized, and the other, 'fountain-like' emitting region, with outflows-inflows nearly orthogonal to the disc, which become suppressed with stronger gravitational influence. This subregion mostly contributes to the emission of the Mg II broad line wings. The kinematics of the Mg II core emitting region is similar to that of the H$\beta$ broad line region (seems to be virialized), and therefore the Full Width at Half Maximum (FWHM) of Mg II still can be used for the black hole (BH) mass estimation in the case where the Mg II core component is dominant. However, one should be careful with using the \ion{Mg}{ii} broad line for the BH mass estimation in the case of very large widths (FWHM$>$ 6000 km s$^{-1}$) and/or in the case of strong blue asymmetry., Comment: accepted in MNRAS

The photometric light curve of PG1302-102 shows periodic variability which makes this object one of the most plausible supermassive black hole binary candidate. Interestingly, the most recent study of its updated optical light curve reports a decrease in significance of periodicity which may suggest that the binary model is less favorable. Here, we model the PG 1302-102 light curve, spanning almost 20 years, with a supermassive black hole binary system in which a perturbation in the accretion disk of more massive component is present. Our model reproduces well the observed light curve with a slight perturbation of a sinusoidal feature and predicts that a slightly larger period than previously reported, of about 1899 days, could arise due to a cold spot in the disk of {\bf more massive} component of a close, unequal-mass ($\frac{\mathrm{m}_{1}}{\mathrm{m}_{2}}=0.1$) black hole binary system. The light curve resembles the pattern of sinusoid-like shape within a few years, which could be confirmed by future observations. In addition, using our hybrid method for periodicity detection, we show that the periods in the observed ($1972\pm 254$ days) and modeled ($1873 \pm \ 250$ days) light curves are within one-sigma, which is also consistent with our physical model prediction and with previous findings. Thus, both the periodic nature and its slight fluctuation of the light curve of PG1302-102 is evident from our physical model and confirmed by the hybrid method for periodicity detection., Comment: Accepted to ApJ