Researcher profile

M. Baloković

M. Baloković contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - Emerging
22works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

22 published item(s)

preprint2016arXiv

A growth-rate indicator for Compton-thick active galactic nuclei

Due to their heavily obscured central engines, the growth rate of Compton-thick (CT) active galactic nuclei (AGN) is difficult to measure. A statistically significant correlation between the Eddington ratio, λ$_{Edd}$, and the X-ray power-law index, Γ, observed in unobscured AGN offers an estimate of their growth rate from X-ray spectroscopy (albeit with large scatter). However, since X-rays undergo reprocessing by Compton scattering and photoelectric absorption when the line-of-sight to the central engine is heavily obscured, the recovery of the intrinsic Γ is challenging. Here we study a sample of local, predominantly Compton-thick megamaser AGN, where the black hole mass, and thus Eddington luminosity, are well known. We compile results on X-ray spectral fitting of these sources with sensitive high-energy (E> 10 keV) NuSTAR data, where X-ray torus models which take into account the reprocessing effects have been used to recover the intrinsic Γ values and X-ray luminosities, L$_X$. With a simple bolometric correction to L$_X$ to calculate λ$_{Edd}$, we find a statistically significant correlation between Γ and λ$_{Edd}$ (p = 0.007). A linear fit to the data yields Γ = (0.41$\pm$0.18)log$_{10}$λ$_{Edd}$+(2.38$\pm$ 0.20), which is statistically consistent with results for unobscured AGN. This result implies that torus modeling successfully recovers the intrinsic AGN parameters. Since the megamasers have low-mass black holes (M$_{BH}\approx10^6-10^7$ M$_{sol}$) and are highly inclined, our results extend the Γ-λ$_{Edd}$ relationship to lower masses and argue against strong orientation effects in the corona, in support of AGN unification. Finally this result supports the use of Γ as a growth-rate indicator for accreting black holes, even for Compton-thick AGN.

preprint2016arXiv

A search for spectral hysteresis and energy-dependent time lags from X-ray and TeV gamma-ray observations of Mrk 421

Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three "target-of-opportunity" (ToO) observations of Mrk 421, one of the brightest TeV blazars, triggered by a strong flaring event at TeV energies in 2014. These observations feature long, continuous, and simultaneous exposures with XMM-Newton (covering X-ray and optical/ultraviolet bands) and VERITAS (covering TeV gamma-ray band), along with contemporaneous observations from other gamma-ray facilities (MAGIC and Fermi-LAT) and a number of radio and optical facilities. Although neither rapid flares nor significant X-ray/TeV correlation are detected, these observations reveal subtle changes in the X-ray spectrum of the source over the course of a few days. We search the simultaneous X-ray and TeV data for spectral hysteresis patterns and time delays, which could provide insight into the emission mechanisms and the source properties (e.g. the radius of the emitting region, the strength of the magnetic field, and related timescales). The observed broadband spectra are consistent with a one-zone synchrotron self-Compton model. We find that the power spectral density distribution at $\gtrsim 4\times 10^{-4}$ Hz from the X-ray data can be described by a power-law model with an index value between 1.2 and 1.8, and do not find evidence for a steepening of the power spectral index (often associated with a characteristic length scale) compared to the previously reported values at lower frequencies.

preprint2016arXiv

First NuSTAR observations of the BL Lac - type blazar PKS~2155-304: constraints on the jet content and distribution of radiating particles

We report the first hard X-ray observations with NuSTAR of the BL Lac type blazar PKS 2155-304, augmented with soft X-ray data from XMM-Newton and gamma-ray data from the Fermi Large Area Telescope, obtained in April 2013 when the source was in a very low flux state. A joint NuSTAR and XMM spectrum, covering the energy range 0.5 - 60 keV, is best described by a model consisting of a log-parabola component with curvature beta = 0.3(+0.2,-0.1) and a (local) photon index 3.04 +/- 0.15 at photon energy of 2 keV, and a hard power-law tail with photon index 2.2 +/- 0.4. The hard X-ray tail can be smoothly joined to the quasi-simultaneous gamma-ray spectrum by a synchrotron self-Compton component produced by an electron distribution with index p = 2.2. Assuming that the power-law electron distribution extends down to the minimum electron Lorentz factor gamma_min = 1 and that there is one proton per electron, an unrealistically high total jet power L_p of roughly 10^47 erg/s is inferred. This can be reduced by two orders of magnitude either by considering a significant presence of electron-positron pairs with lepton-to-proton ratio of at least 30, or by introducing an additional, low-energy break in the electron energy distribution at the electron Lorentz factor gamma_br1 of roughly 100. In either case, the jet composition is expected to be strongly matter-dominated.

preprint2016arXiv

NuSTAR observations of water megamaser AGN

Aims. Study the connection between the masing disk and obscuring torus in Seyfert 2 galaxies. Methods. We present a uniform X-ray spectral analysis of the high energy properties of 14 nearby megamaser Active Galactic Nuclei observed by NuSTAR. We use a simple analytical model to localize the maser disk and understand its connection with the torus by combining NuSTAR spectral parameters with available physical quantities from VLBI mapping. Results. Most of the sources analyzed are heavily obscured, showing a column density in excess of $\sim 10^{23}$ cm$^{-2}$. In particular, $79\%$ are Compton-thick ($N_{\rm H} > 1.5 \times 10^{24}$ cm$^{-2}$). Using column densities measured by NuSTAR, with the assumption that the torus is the extension of the maser disk, and further assuming a reasonable density profile, the torus dimensions can be predicted. They are found to be consistent with mid-IR interferometry parsec-scale observations of Circinus and NGC 1068. In this picture, the maser disk is intimately connected to the inner part of the torus. It is probably made of a large number of molecular clouds connecting the torus and the outer part of the accretion disk, giving rise to a thin disk rotating in most cases in Keplerian or sub-Keplerian motion. This toy model explains the established close connection between water megamaser emission and nuclear obscuration as a geometric effect.

preprint2016arXiv

RoboPol: Do optical polarization rotations occur in all blazars?

We present a new set of optical polarization plane rotations in blazars, observed during the third year of operation of RoboPol. The entire set of rotation events discovered during three years of observations is analysed with the aim of determining whether these events are inherent in all blazars. It is found that the frequency of the polarization plane rotations varies widely among blazars. This variation cannot be explained either by a difference in the relativistic boosting or by selection effects caused by a difference in the average fractional polarization. We conclude that the rotations are characteristic of a subset of blazars and that they occur as a consequence of their intrinsic properties.

preprint2016arXiv

RoboPol: First season rotations of optical polarization plane in blazars

We present first results on polarization swings in optical emission of blazars obtained by RoboPol, a monitoring program of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events. A possible connection of polarization swing events with periods of high activity in gamma rays is investigated using the dataset obtained during the first season of operation. It was found that the brightest gamma-ray flares tend to be located closer in time to rotation events, which may be an indication of two separate mechanisms responsible for the rotations. Blazars with detected rotations have significantly larger amplitude and faster variations of polarization angle in optical than blazars without rotations. Our simulations show that the full set of observed rotations is not a likely outcome (probability $\le 1.5 \times 10^{-2}$) of a random walk of the polarization vector simulated by a multicell model. Furthermore, it is highly unlikely ($\sim 5 \times 10^{-5}$) that none of our rotations is physically connected with an increase in gamma-ray activity.

preprint2016arXiv

RoboPol: optical polarization-plane rotations and flaring activity in blazars

We present measurements of rotations of the optical polarization of blazars during the second year of operation of RoboPol, a monitoring programme of an unbiased sample of gamma-ray bright blazars specially designed for effective detection of such events, and we analyse the large set of rotation events discovered in two years of observation. We investigate patterns of variability in the polarization parameters and total flux density during the rotation events and compare them to the behaviour in a non-rotating state. We have searched for possible correlations between average parameters of the polarization-plane rotations and average parameters of polarization, with the following results: (1) there is no statistical association of the rotations with contemporaneous optical flares; (2) the average fractional polarization during the rotations tends to be lower than that in a non-rotating state; (3) the average fractional polarization during rotations is correlated with the rotation rate of the polarization plane in the jet rest frame; (4) it is likely that distributions of amplitudes and durations of the rotations have physical upper bounds, so arbitrarily long rotations are not realised in nature.

preprint2016arXiv

The Pan-STARRS1 distant z>5.6 quasar survey: more than 100 quasars within the first Gyr of the universe

Luminous quasars at z>5.6 can be studied in detail with the current generation of telescopes and provide us with unique information on the first gigayear of the universe. Thus far these studies have been statistically limited by the number of quasars known at these redshifts. Such quasars are rare and therefore wide-field surveys are required to identify them and multiwavelength data are needed to separate them efficiently from their main contaminants, the far more numerous cool dwarfs. In this paper, we update and extend the selection for z~6 quasars presented in Banados et al. (2014) using the Pan-STARRS1 (PS1) survey. We present the PS1 distant quasar sample, which currently consists of 124 quasars in the redshift range 5.6<z<6.7 that satisfy our selection criteria. Seventy-seven of these quasars have been discovered with PS1, and 63 of them are newly identified in this paper. We present composite spectra of the PS1 distant quasar sample. This sample spans a factor of ~20 in luminosity and shows a variety of emission line properties. The number of quasars at z>5.6 presented in this work almost double the quasars previously known at these redshifts, marking a transition phase from studies of individual sources to statistical studies of the high-redshift quasar population, which was impossible with earlier, smaller samples.

preprint2015arXiv

A combined radio and GeV gamma-ray view of the 2012 and 2013 flares of Mrk 421

In 2012 Markarian 421 underwent the largest flare ever observed in this blazar at radio frequencies. In the present study, we start exploring this unique event and compare it to a less extreme event in 2013. We use 15 GHz radio data obtained with the Owens Valley Radio Observatory 40-m telescope, 95 GHz millimeter data from the Combined Array for Research in Millimeter-Wave Astronomy, and GeV gamma-ray data from the Fermi Gamma-ray Space Telescope. The radio light curves during the flaring periods in 2012 and 2013 have very different appearances, both in shape and peak flux density. Assuming that the radio and gamma-ray flares are physically connected, we attempt to model the most prominent sub-flares of the 2012 and 2013 activity periods by using the simplest possible theoretical framework. We first fit a one-zone synchrotron self-Compton (SSC) model to the less extreme 2013 flare and estimate parameters describing the emission region. We then model the major gamma-ray and radio flares of 2012 using the same framework. The 2012 gamma-ray flare shows two distinct spikes of similar amplitude, so we examine scenarios associating the radio flare with each spike in turn. In the first scenario, we cannot explain the sharp radio flare with a simple SSC model, but we can accommodate this by adding plausible time variations to the Doppler beaming factor. In the second scenario, a varying Doppler factor is not needed, but the SSC model parameters require fine tuning. Both alternatives indicate that the sharp radio flare, if physically connected to the preceding gamma-ray flares, can be reproduced only for a very specific choice of parameters.

preprint2015arXiv

Locating the gamma-ray emission site in Fermi/LAT blazars. II Multifrequency correlations

In an attempt to constrain and understand the emission mechanism of gamma rays, we perform a cross-correlation analysis of 15 blazars using light curves in millimetre, optical and gamma rays. We use discrete correlation function and consider only correlations significant at 99 per cent level. A strong correlation was found between 37 and 95 GHz with a near-zero time delay in most of the sources, and ~1 month or longer in the rest. A similar result was obtained between the optical and gamma-ray bands. Of the 15 sources, less than 50 per cent showed a strong correlation between the millimetre and gamma-ray or millimetre and optical bands. The primary reason for the lack of statistically significant correlation is the absence of a major outburst in the millimetre bands of most of the sources during the 2.5 yr time period investigated in our study. This may indicate that only the long-term variations or large flares are correlated between these bands. The variability of the sources at every waveband was also inspected using fractional rms variability. The fractional rms variability displays an increase with frequency reaching its maximum in the gamma rays.

preprint2015arXiv

NuSTAR Observations of the Compton-thick Active Galactic Nucleus and Ultraluminous X-ray Source Candidate in NGC 5643

We present two NuSTAR observations of the local Seyfert 2 active galactic nucleus (AGN) and an ultraluminous X-ray source (ULX) candidate in NGC 5643. Together with archival data from Chandra, XMM-Newton and Swift-BAT, we perform a high-quality broadband spectral analysis of the AGN over two decades in energy ($\sim$0.5-100 keV). Previous X-ray observations suggested that the AGN is obscured by a Compton-thick (CT) column of obscuring gas along our line-of-sight. However, the lack of high-quality $\gtrsim$ 10 keV observations, together with the presence of a nearby X-ray luminous source, NGC 5643 X-1, had left significant uncertainties in the characterization of the nuclear spectrum. NuSTAR now enables the AGN and NGC 5643 X-1 to be separately resolved above 10 keV for the first time and allows a direct measurement of the absorbing column density toward the nucleus. The new data show that the nucleus is indeed obscured by a CT column of $N_{\rm{H}}$ $\gtrsim$ 5 $\times$ 10$^{24}$ cm$^{-2}$. The range of 2-10 keV absorption-corrected luminosity inferred from the best fitting models is $L_{2-10,\rm{int}} =$ (0.8-1.7) $\times$ 10$^{42}$ erg s$^{-1}$, consistent with that predicted from multiwavelength intrinsic luminosity indicators. We also study the NuSTAR data for NGC 5643 X-1, and show that it exhibits evidence for a spectral cut-off at energy, $E$ $\sim$ 10 keV, similar to that seen in other ULXs observed by NuSTAR. Along with the evidence for significant X-ray luminosity variations in the 3-8 keV band from 2003-2014, our results further strengthen the ULX classification of NGC 5643 X-1.

preprint2015arXiv

Rapid Variability of Blazar 3C 279 during Flaring States in 2013-2014 with Joint Fermi-LAT, NuSTAR, Swift, and Ground-Based Multi-wavelength Observations

We report the results of a multi-band observing campaign on the famous blazar 3C 279 conducted during a phase of increased activity from 2013 December to 2014 April, including first observations of it with NuSTAR. The $γ$-ray emission of the source measured by Fermi-LAT showed multiple distinct flares reaching the highest flux level measured in this object since the beginning of the Fermi mission, with $F(E > 100\,{\rm MeV})$ of $10^{-5}$ photons cm$^{-2}$ s$^{-1}$, and with a flux doubling time scale as short as 2 hours. The $γ$-ray spectrum during one of the flares was very hard, with an index of $Γ_γ= 1.7 \pm 0.1$, which is rarely seen in flat spectrum radio quasars. The lack of concurrent optical variability implies a very high Compton dominance parameter $L_γ/L_{\rm syn} > 300$. Two 1-day NuSTAR observations with accompanying Swift pointings were separated by 2 weeks, probing different levels of source activity. While the 0.5$-$70 keV X-ray spectrum obtained during the first pointing, and fitted jointly with Swift-XRT is well-described by a simple power law, the second joint observation showed an unusual spectral structure: the spectrum softens by $ΔΓ_{\rm X} \simeq 0.4$ at $\sim$4 keV. Modeling the broad-band SED during this flare with the standard synchrotron plus inverse Compton model requires: (1) the location of the $γ$-ray emitting region is comparable with the broad line region radius, (2) a very hard electron energy distribution index $p \simeq 1$, (3) total jet power significantly exceeding the accretion disk luminosity $L_{\rm j}/L_{\rm d} \gtrsim 10$, and (4) extremely low jet magnetization with $L_{\rm B}/L_{\rm j} \lesssim 10^{-4}$. We also find that single-zone models that match the observed $γ$-ray and optical spectra cannot satisfactorily explain the production of X-ray emission.

preprint2015arXiv

The hard X-ray emission of the luminous infrared galaxy NGC 6240 as observed by NuSTAR

We present a broad-band (~0.3-70 keV) spectral and temporal analysis of NuSTAR observations of the luminous infrared galaxy NGC 6240, combined with archival Chandra, XMM-Newton and BeppoSAX data. NGC 6240 is a galaxy in a relatively early merger state with two distinct nuclei separated by ~1.&#34;5. Previous Chandra observations have resolved the two nuclei, showing that they are both active and obscured by Compton-thick material. Although they cannot be resolved by NuSTAR, thanks to the unprecedented quality of the NuSTAR data at energies >10 keV, we clearly detect, for the first time, both the primary and the reflection continuum components. The NuSTAR hard X-ray spectrum is dominated by the primary continuum piercing through an absorbing column density which is mildly optically thick to Compton scattering (tau ~ 1.2, N_H ~ 1.5 x 10^(24) cm^-2). We detect moderate hard X-ray (> 10 keV) flux variability up to 20% on short (15-20 ksec) timescales. The amplitude of the variability is maximum at ~30 keV and is likely to originate from the primary continuum of the southern nucleus. Nevertheless, the mean hard X-ray flux on longer timescales (years) is relatively constant. Moreover, the two nuclei remain Compton-thick, although we find evidence of variability of the material along the line of sight with column densities N_H <~ 2 x 10^(23) cm-2 over long (~3-15 years) timescales. The observed X-ray emission in the NuSTAR energy range is fully consistent with the sum of the best-fit models of the spatially resolved Chandra spectra of the two nuclei.

preprint2015arXiv

The NuSTAR view of reflection and absorption in NGC 7582

NGC 7582 is a well-studied X-ray bright Seyfert 2 with moderately heavy ($N_{\text{H}}\sim10^{23}-10^{24}$~cm$^{-2}$), highly variable absorption and strong reflection spectral features. The spectral shape changed around the year 2000, dropping in observed flux and becoming much more highly absorbed. Two scenarios have been put forth to explain this spectral change: 1) the central X-ray source partially ``shut off&#39;&#39; around this time, decreasing in intrinsic luminosity, with a delayed decrease in reflection features due to the light-crossing time of the Compton-thick material or 2) the source became more heavily obscured, with only a portion of the power law continuum leaking through. NuSTAR observed NGC~7582 twice in 2012, two weeks apart, in order to quantify the reflection using high-quality data above 10 keV. We find that the most plausible scenario is that NGC 7582 has recently become more heavily absorbed by a patchy torus with a covering fraction of $\sim\,80-90\%$ and an equatorial column density of $\sim 3 \times10^{24}$ cm$^{-2}$. We find the need for an additional highly variable full-covering absorber with $N_{\text{H}}= 4-6 \times10^{23}$ cm$^{-2}$ in the line of sight, possibly associated with a hidden broad line region.

preprint2014arXiv

Coronal Properties of the Seyfert 1.9 Galaxy MCG -05-23-016 Determined from Hard X-ray Spectroscopy with NuSTAR

Measurements of the high-energy cut-off in the coronal continuum of active galactic nuclei have long been elusive for all but a small number of the brightest examples. We present a direct measurement of the cut-off energy in the nuclear continuum of the nearby Seyfert 1.9 galaxy MCG -05-23-016 with unprecedented precision. The high sensitivity of NuSTAR up to 79 keV allows us to clearly disentangle the spectral curvature of the primary continuum from that of its reflection component. Using a simple phenomenological model for the hard X-ray spectrum, we constrain the cut-off energy to $116_{-5}^{+6}$ keV with 90% confidence. Testing for more complex models and nuisance parameters that could potentially influence the measurement, we find that the cut-off is detected robustly. We further use simple Comptonized plasma models to provide independent constraints for both the kinetic temperature of the electrons in the corona and its optical depth. At the 90% confidence level, we find $kT_e=29\pm2$ keV and $τ_e=1.23\pm0.08$ assuming a slab (disk-like) geometry, and $kT_e=25\pm2$ keV and $τ_e=3.5\pm0.2$ assuming a spherical geometry. Both geometries are found to fit the data equally well and their two principal physical parameters are correlated in both cases. With the optical depth in the optically thick regime, the data are pushing the currently available theoretical models of the Comptonized plasma to the limits of their validity. Since the spectral features and variability arising from the inner accretion disk have been observed previously in MCG -05-23-016, the inferred high optical depth implies that a spherical or disk-like corona cannot be homogeneous.

preprint2014arXiv

NuSTAR J033202-2746.8: direct constraints on the Compton reflection in a heavily obscured quasar at z~2

We report NuSTAR observations of NuSTAR J033202-2746.8, a heavily obscured, radio-loud quasar detected in the Extended Chandra Deep Field-South, the deepest layer of the NuSTAR extragalactic survey (~400 ks, at its deepest). NuSTAR J033202-2746.8 is reliably detected by NuSTAR only at E>8 keV and has a very flat spectral slope in the NuSTAR energy band (Gamma=0.55^{+0.62}_{-0.64}; 3-30 keV). Combining the NuSTAR data with extremely deep observations by Chandra and XMM-Newton (4 Ms and 3 Ms, respectively), we constrain the broad-band X-ray spectrum of NuSTAR J033202-2746.8, indicating that this source is a heavily obscured quasar (N_H=5.6^{+0.9}_{-0.8}x10^23 cm^-2) with luminosity L_{10-40 keV}~6.4x10^44 erg s^-1. Although existing optical and near-infrared (near-IR) data, as well as follow-up spectroscopy with the Keck and VLT telescopes, failed to provide a secure redshift identification for NuSTAR J033202-2746.8, we reliably constrain the redshift z=2.00+/-0.04 from the X-ray spectral features (primarily from the iron K edge). The NuSTAR spectrum shows a significant reflection component (R=0.55^{+0.44}_{-0.37}), which was not constrained by previous analyses of Chandra and XMM-Newton data alone. The measured reflection fraction is higher than the R~0 typically observed in bright radio-loud quasars such as NuSTAR J033202-2746.8, which has L_{1.4 GHz}~10^27 W Hz^-1. Constraining the spectral shape of AGN, including bright quasars, is very important for understanding the AGN population, and can have a strong impact on the modeling of the X-ray background. Our results show the importance of NuSTAR in investigating the broad-band spectral properties of quasars out to high redshift.

preprint2014arXiv

NuSTAR unveils a Compton-thick Type 2 quasar in Mrk 34

We present Nustar 3-40 keV observations of the optically selected Type 2 quasar (QSO2) SDSS J1034+6001 or Mrk 34. The high-quality hard X-ray spectrum and archival XMM-Newton data can be fitted self-consistently with a reflection-dominated continuum and strong Fe Kalpha fluorescence line with equivalent-width >1 keV. Prior X-ray spectral fitting below 10 keV showed the source to be consistent with being obscured by Compton-thin column densities of gas along the line-of-sight, despite evidence for much higher columns from multiwavelength data. NuSTAR now enables a direct measurement of this column, and shows that Nh lies in the Compton-thick (CT) regime. The new data also show a high intrinsic 2-10 keV luminosity of L_{2-10}~10^{44} erg/s, in contrast to previous low-energy X-ray measurements for which L_{2-10}<~10^{43} erg/s (i.e. X-ray selection below 10 keV does not pick up this source as an intrinsically luminous obscured quasar). Both the obscuring column and the intrinsic power are about an order of magnitude (or more) larger than inferred from pre-NuSTAR X-ray spectral fitting. Mrk34 is thus a &#39;gold standard&#39; CT QSO2 and is the nearest non-merging system in this class, in contrast to the other local CT quasar NGC6240 which is currently undergoing a major merger coupled with strong star-formation. For typical X-ray bolometric correction factors, the accretion luminosity of Mrk34 is high enough to potentially power the total infrared luminosity. X-ray spectral fitting also shows that thermal emission related to star-formation is unlikely to drive the observed bright soft component below ~3 keV, favoring photionization instead.

preprint2014arXiv

The NuSTAR View of Nearby Compton-thick AGN: The Cases of NGC 424, NGC 1320 and IC 2560

We present X-ray spectral analyses for three Seyfert 2 active galactic nuclei, NGC 424, NGC 1320, and IC 2560, observed by NuSTAR in the 3-79 keV band. The high quality hard X-ray spectra allow detailed modeling of the Compton reflection component for the first time in these sources. Using quasi-simultaneous NuSTAR and Swift/XRT data, as well as archival XMM-Newton data, we find that all three nuclei are obscured by Compton-thick material with column densities in excess of ~5 x $10^{24}$ cm$^{-2}$, and that their X-ray spectra above 3 keV are dominated by reflection of the intrinsic continuum on Compton-thick material. Due to the very high obscuration, absorbed intrinsic continuum components are not formally required by the data in any of the sources. We constrain the intrinsic photon indices and the column density of the reflecting medium through the shape of the reflection spectra. Using archival multi-wavelength data we recover the intrinsic X-ray luminosities consistent with the broadband spectral energy distributions. Our results are consistent with the reflecting medium being an edge-on clumpy torus with a relatively large global covering factor and overall reflection efficiency of the order of 1%. Given the unambiguous confirmation of the Compton-thick nature of the sources, we investigate whether similar sources are likely to be missed by commonly used selection criteria for Compton-thick AGN, and explore the possibility of finding their high-redshift counterparts.

preprint2014arXiv

The RoboPol optical polarization survey of gamma-ray - loud blazars

We present first results from RoboPol, a novel-design optical polarimeter operating at the Skinakas Observatory in Crete. The data, taken during the May - June 2013 commissioning of the instrument, constitute a single-epoch linear polarization survey of a sample of gamma-ray - loud blazars, defined according to unbiased and objective selection criteria, easily reproducible in simulations, as well as a comparison sample of, otherwise similar, gamma-ray - quiet blazars. As such, the results of this survey are appropriate for both phenomenological population studies and for tests of theoretical population models. We have measured polarization fractions as low as $0.015$ down to $R$ magnitude of 17 and as low as $0.035$ down to 18 magnitude. The hypothesis that the polarization fractions of gamma-ray - loud and gamma-ray - quiet blazars are drawn from the same distribution is rejected at the $10^{-3}$ level. We therefore conclude that gamma-ray - loud and gamma-ray - quiet sources have different optical polarization properties. This is the first time this statistical difference is demonstrated in optical wavelengths. The polarization fraction distributions of both samples are well-described by exponential distributions with averages of $\langle p \rangle =6.4 ^{+0.9}_{-0.8}\times 10^{-2}$ for gamma-ray--loud blazars, and $\langle p \rangle =3.2 ^{+2.0}_{-1.1}\times 10^{-2}$ for gamma-ray--quiet blazars. The most probable value for the difference of the means is $3.4^{+1.5}_{-2.0}\times 10^{-2}$. The distribution of polarization angles is statistically consistent with being uniform.

preprint2014arXiv

The RoboPol Pipeline and Control System

We describe the data reduction pipeline and control system for the RoboPol project. The RoboPol project is monitoring the optical $R$-band magnitude and linear polarization of a large sample of active galactic nuclei that is dominated by blazars. The pipeline calibrates and reduces each exposure frame, producing a measurement of the magnitude and linear polarization of every source in the $13&#39;\times 13&#39;$ field of view. The control system combines a dynamic scheduler, real-time data reduction, and telescope automation to allow high-efficiency unassisted observations.

preprint2014arXiv

The Seyfert 2 galaxy NGC 2110: hard X-ray emission observed by NuSTAR and variability of the iron K$α$ line

We present NuSTAR observations of the bright Seyfert 2 galaxy NGC 2110 obtained in 2012, when the source was at the highest flux level ever observed, and in 2013, when the source was at a more typical flux level. We include archival observations from other X-ray satellites, namely XMM-Newton, Suzaku, BeppoSAX, Chandra and Swift. Simultaneous NuSTAR and Swift broad band spectra (in the 3-80 keV range) indicate a cutoff energy $E_{\rm c}>210$ keV, with no detectable contribution from Compton reflection. NGC 2110 is one of the very few sources where no evidence for distant Compton thick scattering is found and, by using temporal information collected over more than a decade, we investigate variations of the iron K$α$ line on time scales of years. The Fe K$α$ line is likely the sum of two components: one constant (originating from distant Compton-thick material) and the other one variable and linearly correlated with the source flux (possibly arising from Compton-thin material much closer to the black hole).

preprint2013arXiv

First Results from NuSTAR Observations of Mkn 421

Mkn 421 is a nearby active galactic nucleus dominated at all wavelengths by a very broad non-thermal continuum thought to arise from a relativistic jet seen at a small angle to the line of sight. Its spectral energy distribution peaks in the X-ray and TeV gamma-ray bands, where the energy output is dominated by cooling of high-energy electrons in the jet. In order to study the electron distribution and its evolution, we carried out a dedicated multi-wavelength campaign, including extensive observations by the recently launched highly sensitive hard X-ray telescope NuSTAR, between December 2012 and May 2013. Here we present some initial results based on NuSTAR data from January through March 2013, as well as calibration observations conducted in 2012. Although the observations cover some of the faintest hard X-ray flux states ever observed for Mkn 421, the sensitivity is high enough to resolve intra-day spectral variability. We find that in this low state the dominant flux variations are smooth on timescales of hours, with typical intra-hour variations of less than 5%. We do not find evidence for either a cutoff in the hard X-ray spectrum, or a rise towards a high-energy component, but rather that at low flux the spectrum assumes a power law shape with a photon index of approximately 3. The spectrum is found to harden with increasing brightness.