Researcher profile

Marko Stalevski

Marko Stalevski contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
11works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

11 published item(s)

preprint2022arXiv

A multiwavelength-motivated X-ray model for the Circinus Galaxy

Reprocessed X-ray emission in Active Galactic Nuclei (AGN) can provide fundamental information about the circumnuclear environments of supermassive black holes. Recent mid-infrared studies have shown evidence of an extended dusty structure perpendicular to the torus plane. In this work, we build a self-consistent X-ray model for the Circinus Galaxy including the different physical components observed at different wavelengths and needed to reproduce both the morphological and spectral properties of this object in the mid-infrared. The model consists of four components: the accretion disk, the broad line region (BLR), a flared disk in the equatorial plane and a hollow cone in the polar direction. Our final model reproduces well the 3--70 keV Chandra and NuSTAR spectra of Circinus, including the complex Fe K$α$ zone and the spectral curvature, although several additional Gaussian lines, associated to either ionized iron or to broadened Fe K$α$/K$β$ lines, are needed. We find that the flared disk is Compton thick ($ N_{\rm H,d}= \rm 1.01^{+0.03}_{-0.24}\times 10^{25}\: cm^{-2}$) and geometrically thick ($CF=0.55^{+0.01}_{-0.05}$), and that the hollow cone has a Compton-thin column density ($ N_{\rm H,c}= \rm 2.18^{+0.47}_{-0.43}\times 10^{23}\: cm^{-2}$), which is consistent with the values inferred by mid-infrared studies. Including also the BLR, the effective line of sight column density is $ N_{\rm H}= \rm 1.47^{+0.03}_{-0.24}\times 10^{25}\: cm^{-2}$. This approach to X-ray modelling, i.e. including all the different reprocessing structures, will be very important to fully exploit data from future X-ray missions.

preprint2022arXiv

BASS XXVIII: Near-infrared Data Release 2, High-Ionization and Broad Lines in Active Galactic Nuclei

We present the BAT AGN Spectroscopic Survey (BASS) Near-infrared Data Release 2 (DR2), a study of 168 nearby ($\bar z$ = 0.04, $z$ < 0.6) active galactic nuclei (AGN) from the all-sky Swift Burst Array Telescope X-ray survey observed with Very Large Telescope (VLT)/X-shooter in the near-infrared (NIR; 0.8 - 2.4 $μ$m). We find that 49/109 (45%) Seyfert 2 and 35/58 (60%) Seyfert 1 galaxies observed with VLT/X-shooter show at least one NIR high-ionization coronal line (CL, ionization potential $χ$ > 100 eV). Comparing the emission of the [Si vi] $λ$1.9640 CL with the X-ray emission for the DR2 AGN, we find a significantly tighter correlation, with a lower scatter (0.37 dex) than for the optical [O iii] $λ$5007 line (0.71 dex). We do not find any correlation between CL emission and the X-ray photon index $Γ$. We find a clear trend of line blueshifts with increasing ionization potential in several CLs, such as [Si vi] $λ$1.9640, [Si x] $λ$1.4300, [S viii] $λ$0.9915, and [S ix] $λ$1.2520, indicating the radial structure of the CL region. Finally, we find a strong underestimation bias in black hole mass measurements of Sy 1.9 using broad H$α$ due to the presence of significant dust obscuration. In contrast, the broad Pa$α$ and Pa$β$ emission lines are in agreement with the $M$-$σ$ relation. Based on the combined DR1 and DR2 X-shooter sample, the NIR BASS sample now comprises 266 AGN with rest-frame NIR spectroscopic observations, the largest set assembled to date.

preprint2022arXiv

Fitting AGN/galaxy X-ray-to-radio SEDs with CIGALE and improvement of the code

Modern and future surveys effectively provide a panchromatic view for large numbers of extragalactic objects. Consistently modeling these multiwavelength survey data is a critical but challenging task for extragalactic studies. The Code Investigating GALaxy Emission (CIGALE) is an efficient PYTHON code for spectral energy distribution (SED) fitting of galaxies and active galactic nuclei (AGNs). Recently, a major extension of CIGALE (named X-CIGALE) has been developed to account for AGN/galaxy X-ray emission and improve AGN modeling at UV-to-IR wavelengths. Here, we apply X-CIGALE to different samples, including COSMOS spectroscopic type 2 AGNs, CDF-S X-ray detected normal galaxies, SDSS quasars, and COSMOS radio objects. From these tests, we identify several weaknesses of X-CIGALE and improve the code accordingly. These improvements are mainly related to AGN intrinsic X-ray anisotropy, X-ray binary emission, AGN accretion-disk SED shape, and AGN radio emission. These updates improve the fit quality and allow new interpretation of the results, based on which we discuss physical implications. For example, we find that AGN intrinsic X-ray anisotropy is moderate, and can be modeled as $L_X(θ) \propto 1+\cos θ$, where $θ$ is the viewing angle measured from the AGN axis. We merge the new code into the major branch of CIGALE, and publicly release this new version as CIGALE v2022.0 on https://cigale.lam.fr

preprint2022arXiv

The dusty heart of Circinus: I. Imaging the circumnuclear dust in N-band

Active galactic nuclei play a key role in the evolution of galaxies, but their inner workings and physical connection to the host are poorly understood due to a lack of angular resolution. Infrared interferometry makes it possible to resolve the circumnuclear dust in the nearby Seyfert 2 galaxy, Circinus. Previous observations have revealed complex structures and polar dust emission but interpretation was limited to simple models. MATISSE makes it possible to image these structures for the first time. We observed the Circinus Galaxy with VLTI/MATISSE, producing 150 correlated flux spectra and 100 closure phase spectra. We reconstructed images in the N-band at ~10 mas resolution. We fit blackbody functions with dust extinction to several aperture-extracted fluxes from the images to produce a temperature distribution of central dusty structures. We find significant substructure in the circumnuclear dust: central unresolved flux of ~0.5 Jy, a thin disk 1.9 pc in diameter oriented along ~45 deg,and a ~4x1.5 pc polar emission extending orthogonal to the disk. The polar emission exhibits patchiness, which we attribute to clumpy dust. Flux enhancements to the east and west of the disk are seen for the first time. We distinguish the temperature profiles of the disk and of the polar emission: the disk shows a steep temperature gradient indicative of denser material; the polar profile is flatter, indicating clumpiness and/or lower dust density. The unresolved flux is fitted with a high temperature, ~370 K. The polar dust remains warm (~200 K) out to 1.5 pc from the disk. The recovered morphology and temperature distribution resembles modeling of accretion disks with radiation-driven winds at large scales, but we placed new constraints on the subparsec dust. The subparsec features imaged here place new constraints on the physical modeling of circumnuclear dust in active galaxies.

preprint2021arXiv

BAT AGN Spectroscopic Survey-XXIII. A New Mid-Infrared Diagnostic for Absorption in Active Galactic Nuclei

In this study, we use the SWIFT/BAT AGN sample, which has received extensive multiwavelength follow-up analysis as a result of the BAT AGN Spectroscopic Survey (BASS), to develop a diagnostic for nuclear obscuration by examining the relationship between the line-of-sight column densities ($N_{\rm{H}}$), the 2-10 keV-to-$12\,\rm{μm}$ luminosity ratio, and WISE mid-infrared colors. We demonstrate that heavily obscured AGNs tend to exhibit both preferentially &#39;&#39;redder&#39;&#39; mid-infrared colors and lower values of $L_{\rm{X,\,Obs.}}$/$L_{12\,\rm{μm}}$ than less obscured AGNs, and we derive expressions relating $N_{\rm{H}}$ to the $L_{\rm{X,\,Obs.}}$/$L_{12\,\rm{μm}}$ and $L_{22\,\rm{μm}}$/$L_{4.6\,\rm{μm}}$ luminosity ratios as well as develop diagnostic criteria using these ratios. Our diagnostic regions yield samples that are $\gtrsim80$% complete and $\gtrsim60$% pure for AGNs with log($N_{\rm{H}})\geq24$, as well as $\gtrsim85$% pure for AGNs with $\rm{log}(N_{\rm{H}})\gtrsim23.5$. We find that these diagnostics cannot be used to differentiate between optically star forming galaxies and active galaxies. Further, mid-IR contributions from host galaxies that dominate the observed $12~\rm{μm}$ emission can lead to larger apparent X-ray deficits and redder mid-IR colors than the AGNs would intrinsically exhibit, though this effect helps to better separate less obscured and more obscured AGNs. Finally, we test our diagnostics on two catalogs of AGNs and infrared galaxies, including the XMM-Newton XXL-N field, and we identify several known Compton-thick AGNs as well as a handful of candidate heavily obscured AGNs based upon our proposed obscuration diagnostics.

preprint2014arXiv

The narrow Fe K$α$ line and the molecular torus in active galactic nuclei - an IR/X-ray view

The narrow component of the iron K$α$ is an almost ubiquitous feature in the X-ray spectra of active galactic nuclei (AGN) and is believed to originate in neutral material, possibly located in the molecular torus. This would imply a tight connection between the Fe K$α$ equivalent width (EW) and the physical properties of the torus. In a recent work we have shown that the decrease of the covering factor of the torus with the luminosity, as expected by luminosity-dependent unification models, would be able to explain the decrease of Fe K$α$ EW with the luminosity (i.e., the X-ray Baldwin effect). Recent developments in the study of the mid-IR (MIR) spectrum of AGN allow important parameters of the torus to be deduced, such as its covering factor ($f_{\rm\,obs}$) and equatorial column density ($N_{\rm\,H}^{\rm\,T}$), by applying clumpy torus models. Using XMM-Newton/EPIC observations of a sample of 24 type-I AGN, we investigate the relation between the physical parameters of the torus obtained by recent MIR works and the properties of the Fe K$α$ line. We correct the values of the Fe K$α$ EW by taking the inclination angle, the photon index, the equatorial column density, and half-opening angle of the torus into account using a physical torus model of X-ray reprocessed radiation. We find that the relation between Fe K$α$ EW and $f_{\rm\,obs}$ shows a slope that is consistent with the expected value, albeit with a low statistical significance. A trend that is consistent with the theoretical prediction is also found when comparing the Fe K$α$ EW to $N_{\rm\,H}^{\rm\,T}$. Our work seems to confirm that the bulk of the narrow Fe K$α$ line is produced by the same material responsible for the MIR emission.

preprint2012arXiv

Gravitational microlensing of AGN dusty tori

We investigated the gravitational microlensing of active galactic nucleus dusty tori in the case of lensed quasars in the infrared domain. The dusty torus is modeled as a clumpy two-phase medium. To obtain spectral energy distributions and images of tori at different wavelengths, we used the 3D Monte Carlo radiative transfer code SKIRT. A ray-shooting technique has been used to calculate microlensing magnification maps. We simulated microlensing by the stars in the lens galaxy for different configurations of the lensed system and different values of the torus parameters, in order to estimate (a) amplitudes and timescales of high magnification events, and (b) the influence of geometrical and physical properties of dusty tori on light curves in the infrared domain. We found that, despite their large size, dusty tori could be significantly affected by microlensing in some cases, especially in the near-infrared domain (rest-frame). The very long time-scales of such events, in the range from several decades to hundreds of years, are limiting the practical use of this method to study the properties of dusty tori. However, our results indicate that, when studying flux ratios between the images in different wavebands of lensed quasars, one should not disregard the possibility that the near and mid-infrared flux ratios could be under the influence of microlensing.

preprint2011arXiv

Efficient 3D NLTE dust radiative transfer with SKIRT

We present an updated version of SKIRT, a 3D Monte Carlo radiative transfer code developed to simulate dusty galaxies. The main novel characteristics of the SKIRT code are the use of a stellar foam to generate random positions, an efficient combination of eternal forced scattering and continuous absorption, and a new library approach that links the radiative transfer code to the DustEM dust emission library. This approach enables a fast, accurate and self-consistent calculation of the dust emission of arbitrary mixtures of transiently heated dust grains and polycyclic aromatic hydrocarbons, even for full 3D models containing millions of dust cells. We have demonstrated the accuracy of the SKIRT code through a set of simulations based on the edge-on spiral galaxy UGC 4754. The models we ran were gradually refined from a smooth, 2D, LTE model to a fully 3D model that includes NLTE dust emission and a clumpy structure of the dusty ISM. We find that clumpy models absorb UV and optical radiation less efficiently than smooth models with the same amount of dust, and that the dust in clumpy models is on average both cooler and less luminous. Our simulations demonstrate that, given the appropriate use of optimization techniques, it is possible to efficiently and accurately run Monte Carlo radiative transfer simulations of arbitrary 3D structures of several million dust cells, including a full calculation of the NLTE emission by arbitrary dust mixtures.

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

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

The complex interplay of dust and star light in spiral galaxy discs

Interstellar dust grains efficiently absorb and scatter UV and optical radiation in galaxies, and therefore can significantly affect the apparent structure of spiral galaxies. We discuss the effect of dust attenuation on the observed structural properties of bulges and discs. We also present some first results on modelling the dust content of edge-on spiral galaxies using both optical and Herschel far-infrared data. Both of these results demonstrate the complex interplay of dust and star light in spiral galaxies.