Source author record

Luka C. Popovic

Luka C. Popovic appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

19 published item(s)

preprint2024arXiv

AGN STORM 2: V. Anomalous Behavior of the CIV Light Curve in Mrk 817

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.

preprint2022arXiv

Detection of eccentric close-binary supermassive black holes with incomplete interferometric data

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.

preprint2021arXiv

Conditional Neural Process for non-parametric modeling of AGN light curve

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).

preprint2021arXiv

Tracing the outflow kinematics in Type 2 Active Galactic Nuclei

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.

preprint2020arXiv

Confirmed short periodic variability of subparsec supermassive binary black hole candidate Mrk 231

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.

preprint2020arXiv

Two-dimensional correlation analysis of periodicity in active galactic nuclei time series

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.

preprint2016arXiv

Line shifts and sub-pc super-massive binary black holes

Here we discuss the possibility of super-massive binary black hole (SMBBH) detection, using the shifts of the broad lines emitted from a binary system. We perform a number of simulations of shapes and shifts of $H_β$ lines emitted from SMBBHs, taking into account the emission from two different regions located around both black holes, and kinematical effects which should be present in a SMBBH. In the model we connect the parameters of the lines with the mass of black holes and find that the peak shift depends, not only on kinematical effects of system rotation and black hole mass ratio, but it is also a function of the individual masses of the two black holes (BHs).

preprint2016arXiv

The UV and optical Fe II emission lines in type 1 AGNs

We investigate the spectral properties of the UV ($λλ$2650-3050 Å) and optical ($λλ$4000-5500 Å) Fe II emission features in a sample of 293 type 1 active galactic nuclei (AGNs) from Sloan Digital Sky Survey (SDSS) database. We explore different correlations between their emission line properties, as well as the correlations with the other emission lines from the spectral range. We find several interesting correlations and we can outline the most interesting results as follows. (i) There is a kinematical connection between the UV and optical Fe II lines, indicating that the UV and optical Fe II lines originate from the outer part of the broad line region, so-called intermediate line region; (ii) The unexplained anticorrelations of the optical Fe II (EW Fe II$_{opt}$) versus EW [O III] 5007 Å and EW Fe II$_{opt}$ versus FWHM Hbeta have not been detected for the UV Fe II lines; (iii) The significant averaged redshift in the UV Fe II lines, which is not present in optical Fe II, indicates an inflow in the UV Fe II emitting clouds, and probably their asymmetric distribution. (iv) Also, we confirm the anticorrelation between the intensity ratio of the optical and UV Fe II lines and FWHM of Hbeta, and we find the anticorrelations of this ratio with the widths of Mg II 2800 Å, optical Fe II and UV Fe II. This indicates a very important role for the column density and microturbulence in the emitting gas. We discuss the starburst activity in high-density regions of young AGNs as a possible explanation of the detected optical Fe II correlations and intensity line ratios of the UV and optical Fe II lines.

preprint2016arXiv

Virilization of the Broad Line Region in Active Galactic Nuclei - connection between shifts and widths of broad emission lines

We investigate the virilization of the emission lines Hbeta and Mg II in the sample of 287 Type 1 Active Galactic Nuclei taken from the Sloan Digital Sky Survey database. We explore the connections between the intrinsic line shifts and full widths at different levels of maximal intensity. We found that: (i) Hbeta seems to be a good virial estimator of black hole masses, and an intrinsic redshift of Hbeta is dominantly caused by the gravitational effect, (ii) there is an anti-correlation between the redshift and width of the wings of the Mg II line, (iii) the broad Mg II line can be used as virial estimator only at 50% of the maximal intensity, while the widths and intrinsic shifts of the line wings can not be used for this purpose.

preprint2015arXiv

On the time delay evolution of five Active Galactic Nuclei

Here we investigate light curves of the continuum and emission lines of five type 1 active galactic nuclei (AGN) from our monitoring campaign, to test time-evolution of their time delays.Using both modeled and observed AGN light curves we apply Gaussian-kernel based estimator to capture variation of local patterns of their time evolving delays. The largest variations of time delays of all objects occur in the period when continuum or emission lines luminosity is the highest. However, Gaussian kernel based method shows instability in the case of NGC 5548, 3C 390.3, E1821+643 and NGC 4051 possible due to numerical discrepancies between Damped Random Walk (DRW) time scale of light curves and sliding time windows of the method. The temporal variations of time lags of Arp 102B can correspond to the real nature of the time lag evolution.

preprint2015arXiv

Polarization in lines - a new method for measuring black hole masses in active galaxies

Measuring of the masses of galactic supermassive black holes (SMBHs) is an important task, since they correlate with the host galaxy properties and play an important role in evolution of galaxies. Here we present a new method for measuring of SMBH masses using the polarization of the broad lines emitted from active galactic nuclei (AGNs). We performed spectropolarometric observations of 9 AGNs and find that this method gives measured masses which are in a good agreement with reverberation measurements. An advantage of this method is that it can be used to measure the masses of SMBHs in a consistent way at different cosmological epochs.

preprint2013arXiv

A model for the Balmer pseudocontinuum in spectra of type 1 AGNs

Here we present a new method for subtracting the Balmer pseudocontinuum in the UV part of type 1 AGN spectra. We calculate the intensity of the Balmer pseudocontinuum using the prominent Balmer lines in AGN spectra. We apply the model on a sample of 293 type 1 AGNs from SDSS database, and found that our model of Balmer pseudocontinuum + power law continuum very well fits the majority of the AGN spectra from the sample, while in 15% of AGNs, the model fits reasonable the UV continuum, but a discrepancy between the observed and fitted spectra is noted. Some of the possible reasons for the discrepancy may be a different value for the optical depth in these spectra than used in our model or the influence of the intrinsic reddening.

preprint2013arXiv

Spectral properties of a sample of type 1 AGNs: influence of star formation

To find the spectral properties of AGNs in optical spectral band (around the Hb line) we constructed an Fe II template, that is covering Fe II emission from 4100 A to 5600 A. Using the new Fe II template we explored spectral properties of 302 type 1 AGNs. The most interesting results we found is the correlation of the Hb Full Width at Half Maximum (FWHM) and luminosity for a subsample of type 1 AGNs where the ratio of narrow lines indicates a significant starburst contribution.

preprint2011arXiv

Optical emission line properties of a sample of the broad-line AGNs: the Baldwin effect and eigenvector 1

We divide a sample of 302 type-1 AGNs into two subsamples based on the narrow line [OIII]/Hbeta_{NLR} ratio, since we expect that there will be a stronger starburst (HII region) contribution to the narrow line emission for R=log([OIII]/Hbeta_{NLR})<0.5. For the two samples we {find significant differences in correlations between} spectral properties of objects with $R<0.5$ and R>0.5. {We find similar differences when we divided the sample based on the FWHM ratios of [OIII] and broad Hbeta lines (R_1=log(FWHM[OIII]/FWHM Hbeta_broad)^>_<-0.8), i.e. similar correlations between R>0.5 and R_1<-0.8 subsamples from one side and R<0.5 and R_1>-0.8 subsamples from the other side.} The most interesting difference is in the correlation between the broad Hbeta FWHM and luminosity in the R<0.5 (R_1>-0.8) sample that indicates a connection between the BLR kinematics and photoionization source. We discuss possible effects which can cause such differences in spectral properties of two subsamples.

preprint2011arXiv

Photo-centric variability of quasars caused by variations in their inner structure: Consequences on Gaia measurements

We study the photocenter position variability due to variations in the quasar inner structure. We consider variability in the accretion disk emissivity and torus structure variability due to different illumination by the central source. We discuss possible detection of these effects by Gaia. Observations of the photocenter variability in two AGNs, SDSS J121855+020002 and SDSS J162011+1724327 have been reported and discussed. With investigation of the variations in the quasar inner structure we explore how much this effect can affect the position determination and whether it can be (or not) detected with Gaia mission. We used (a) a model of a relativistic disk, including the perturbation that can increase brightness of a part of the disk, and consequently offset the photocenter position, and (b) a model of a dusty torus which absorbs and re-emits the incoming radiation from accretion disk. We estimated the value of the photocenter offset due to these two effects. We found that perturbations in the inner structure can significantly offset the photocenter. It depends on the characteristics of perturbation and accretion disk and structure of the torus. In the case of two considered QSOs the observed photocenter offsets cannot be explained by variations in the accretion disk and other effects should be considered. We discussed possibility of exploding stars very close to the AGN source, and also possibility that there are two variable sources in the center of these two AGNs that may indicate a binary super-massive black hole system on a kpc (pc) scale. The Gaia mission seems to be very perspective, not only for astrometry, but also for exploring the inner structure of AGNs. We conclude that variations in the quasar inner structure can affect the observed photocenter (up to several mas). There is a chance to observe such effect in the case of bright and low-redshifted QSOs.

preprint2011arXiv

Super-massive binary black holes and emission lines in active galactic nuclei

The broad emission spectral lines emitted from AGNs are our main probe of the geometry and physics of the broad line region (BLR) close to the SMBH. There is a group of AGNs that emits very broad and complex line profiles, showing two displaced peaks, one blueshifted and one redshifted from the systemic velocity defined by the narrow lines, or a single such peak. It has been proposed that such line shapes could indicate a supermassive binary black hole (SMB) system. We discuss here how the presence of an SMB will affect the BLRs of AGNs and what the observational consequences might be. We review previous claims of SMBs based on broad line profiles and find that they may have non-SMB explanations as a consequence of a complex BLR structure. Because of these effects it is very hard to put limits on the number of SMBs from broad line profiles. It is still possible, however, that unusual broad line profiles in combination with other observational effects (line ratios, quasi-periodical oscillations, spectropolarimetry, etc.) could be used for SMBs detection. Some narrow lines (e.g., [O\,III]) in some AGNs show a double-peaked profile. Such profiles can be caused by streams in the Narrow Line Region (NLR), but may also indicate the presence of a kilo-parsec scale mergers. A few objects indicated as double-peaked narrow line emitters are confirmed as kpc-scale margers, but double-peaked narrow line profiles are mostly caused by the complex NLR geometry. We briefly discuss the expected line profile of broad Fe K$α$ that probably originated in the accretion disk(s) around SMBs. Finally we consider rare configurations where a SMB system might be gravitationally lensed by a foreground galaxy, and discuss the expected line profiles in these systems.

preprint2011arXiv

Three-dimensional radiative transfer modeling of AGN dusty tori as a clumpy two-phase medium

We investigate the emission of active galactic nuclei (AGN) dusty tori in the infrared domain. Following theoretical predictions coming from hydrodynamical simulations, we model the dusty torus as a 3D two-phase medium with high-density clumps and low-density medium filling the space between the clumps. Spectral energy distributions (SED) and images of the torus at different wavelengths are obtained using 3D Monte Carlo radiative transfer code SKIRT. Our approach of generating clumpy structure allows us to model tori with single clumps, complex structures of merged clumps or interconnected sponge-like structure. A corresponding set of clumps-only models and models with smooth dust distribution is calculated for comparison. We found that dust distribution, optical depth, clump size and their actual arrangement in the innermost region, all have an impact on the shape of near- and mid-infrared SED. The 10 micron silicate feature can be suppressed for some parameters, but models with smooth dust distribution are also able to produce a wide range of the silicate feature strength. Finally, we find that having the dust distributed in a two-phase medium, might offer a natural solution to the lack of emission in the near-infrared, compared to observed data, which affects clumpy models currently available in the literature.

preprint2010arXiv

Analysis of optical Fe II emission in a sample of AGN spectra

We present a study of optical Fe II emission in 302 AGNs selected from the SDSS. We group the strongest Fe II multiplets into three groups according to the lower term of the transition (b $^4$F, a $^6$S and a $^4$G terms). These correspond approximately to the blue, central, and red part respectively of the "iron shelf" around Hb. We calculate an Fe II template which takes into account transitions into these three terms and an additional group of lines, based on a reconstruction of the spectrum of I Zw 1. This Fe II template gives a more precise fit of the Fe II lines in broad-line AGNs than other templates. We extract Fe II, Ha, Hb, [O III] and [N II] emission parameters and investigate correlations between them. We find that Fe II lines probably originate in an Intermediate Line Region. We notice that the blue, red, and central parts of the iron shelf have different relative intensities in different objects. Their ratios depend on continuum luminosity, FWHM Hb, the velocity shift of Fe II, and the Ha/Hb flux ratio. We examine the dependence of the well-known anti-correlation between the equivalent widths of Fe II and [O III] on continuum luminosity. We find that there is a Baldwin effect for [O III] but an inverse Baldwin effect for the Fe II emission. The [O III]/Fe II ratio thus decreases with L\lambda5100. Since the ratio is a major component of the Boroson and Green eigenvector 1, this implies a connection between the Baldwin effect and eigenvector 1, and could be connected with AGN evolution. We find that spectra are different for Hb FWHMs greater and less than ~3000 km/s, and that there are different correlation coefficients between the parameters.

preprint2003arXiv

Diagnostics of Plasma Properties in Broad Line Region of AGNs

The Boltzmann-plot (BPT) method for laboratory plasma diagnostic was used for a quick estimate of physical conditions in the Broad Line Region (BLR) of 14 Active Galactic Nuclei (AGNs). For the BLR of nine AGNs, where PLTE exist, the estimated electron temperature are in the range T $\sim$ 130000K - 37000K are in good agreement with previous estimates. The estimated electron densities depend on the opacity of the emitting plasma in the BLR, and they are in range from $10^9\rm cm^{-3}$ (for optically thick plasma) to $10^{14}\rm cm^{-3}$ (for optically thin plasma). Although an alternative of PLTE in some AGNs may be very high intrinsic reddening effect, this method may be used for fast insight into physical processes in the BLR prior to applying more sophisticated physical models.