Realistic synthetic observations of theoretical source models are essential for our understanding of real observational data. In using synthetic data, one can verify the extent to which source parameters can be recovered and evaluate how various data corruption effects can be calibrated. These studies are important when proposing observations of new sources, in the characterization of the capabilities of new or upgraded instruments, and when verifying model-based theoretical predictions in a comparison with observational data. We present the SYnthetic Measurement creator for long Baseline Arrays (SYMBA), a novel synthetic data generation pipeline for Very Long Baseline Interferometry (VLBI) observations. SYMBA takes into account several realistic atmospheric, instrumental, and calibration effects. We used SYMBA to create synthetic observations for the Event Horizon Telescope (EHT), a mm VLBI array, which has recently captured the first image of a black hole shadow. After testing SYMBA with simple source and corruption models, we study the importance of including all corruption and calibration effects. Based on two example general relativistic magnetohydrodynamics (GRMHD) model images of M87, we performed case studies to assess the attainable image quality with the current and future EHT array for different weather conditions. The results show that the effects of atmospheric and instrumental corruptions on the measured visibilities are significant. Despite these effects, we demonstrate how the overall structure of the input models can be recovered robustly after performing calibration steps. With the planned addition of new stations to the EHT array, images could be reconstructed with higher angular resolution and dynamic range. In our case study, these improvements allowed for a distinction between a thermal and a non-thermal GRMHD model based on salient features in reconstructed images.

We present results from modeling the broadband X-ray spectrum of the Type 1 AGN 2MASX J193013.80+341049.5 using NuSTAR, Swift and archival XMM-Newton observations. We find this source to be highly X-ray obscured, with column densities exceeding 10$^{23}$ cm$^{-2}$ across all epochs of X-ray observations, spanning an 8 year period. However, the source exhibits prominent broad optical emission lines, consistent with an unobscured Type 1 AGN classification. We fit the X-ray spectra with both phenomenological reflection models and physically-motivated torus models to model the X-ray absorption. We examine the spectral energy distribution of this source and investigate some possible scenarios to explain the mismatch between X-ray and optical classifications. We compare the ratio of reddening to X-ray absorbing column density ($E_{B-V}/N_{H}$) and find that 2MASX J193013.80+341049.5 likely has a much lower dust-to-gas ratio relative to the Galactic ISM, suggesting that the Broad Line Region (BLR) itself could provide the source of extra X-ray obscuration, being composed of low-ionization, dust-free gas.

We present a multi-frequency study of PKS J1222$+$0413 (4C$+$04.42), currently the highest redshift $\gamma$-ray emitting narrow-line Seyfert 1 ($\gamma$-NLS1). We assemble a broad spectral energy distribution (SED) including previously unpublished datasets: X-ray data obtained with the NuSTAR and Neil Gehrels Swift observatories; near-infrared, optical and UV spectroscopy obtained with VLT X-shooter; and multiband radio data from the Effelsberg telescope. These new observations are supplemented by archival data from the literature. We apply physical models to the broadband SED, parameterising the accretion flow and jet emission to investigate the disc-jet connection. PKS J1222$+$0413 has a much greater black hole mass than most other NLS1s, $M_\mathrm{BH}\approx2\times10^{8}$ M$_\odot$, similar to those found in flat spectrum radio quasars (FSRQs). Therefore this source provides insight into how the jets of $\gamma$-NLS1s relate to those of FSRQs.

We present an analysis of a sample of 69 local obscured Swift/Burst Alert Telescope active galactic nuclei (AGNs) with X-ray spectra from NuSTAR and infrared (IR) spectral energy distributions from Herschel and WISE. We combine this X-ray and IR phenomenological modeling and find a significant correlation between reflected hard X-ray emission and IR AGN emission, with suggestive indications that this correlation may be stronger than the one between intrinsic hard X-ray and IR emissions. This relation between the IR and reflected X-ray emission suggests that both are the result of the processing of intrinsic emission from the corona and accretion disk by the same structure. We explore the resulting implications on the underlying distribution of covering fraction for all AGNs, by generating mock observables for the reflection parameter and IR luminosity ratio using empirical relations found for the covering fraction with each quantity. We find that the observed distributions of the reflection parameter and IR-to-X-ray ratio are reproduced with broad distributions centered around covering fractions of at least ~40%-50%, whereas narrower distributions match our observations only when centered around covering fractions of ~70%-80%. Our results are consistent with both independent estimates of the covering fractions of individual objects and the typical covering fraction obtained on the basis of obscured fractions for samples of AGNs. These results suggest that the level of reprocessing in AGNs, including X-ray reflection, is related in a relatively straightforward way to the geometry of the obscuring material.

The corona is a key component of most luminous accreting black holes, carrying 5 - 30 % of the power and in non-jetted Active Galactic Nuclei (AGN), creating all the X-ray emission above $\simeq 1-2$ keV. Its emission illuminates the inner accretion disc, creating the atomic line-rich reflection spectrum used to diagnose and map the accretion flow and measure black hole spin. The corona is likely powered magnetically by the strong differential rotation of the accretion disc and it may be intimately related to relativistic jets. Recent work shows that many black hole coronae may be dominated by electron-positron pairs produced by photon-photon collisions in the compact coronal environment. Despite the corona being an integral component of AGN and black hole binary systems, much is still unknown about the nature of the corona, such as its geometry, location, and the physical mechanisms powering the emission. In this white paper we explore our current understanding of coronal properties, such as its temperature, obtained from measurements with existing hard X-ray telescopes such as NuSTAR, and discuss important questions to be addressed in the coming decade surrounding the nature of the corona. Hard X-ray observations will continue to dispel the mystery of coronae and open up this part of the quasar engine to full understanding.

Accreting supermassive black holes (SMBHs) can exhibit variable emission across the electromagnetic spectrum and over a broad range of timescales. The variability of active galactic nuclei (AGNs) in the ultraviolet and optical is usually at the few tens of per cent level over timescales of hours to weeks. Recently, rare, more dramatic changes to the emission from accreting SMBHs have been observed, including tidal disruption events, 'changing look' AGNs and other extreme variability objects. The physics behind the 're-ignition', enhancement and 'shut-down' of accretion onto SMBHs is not entirely understood. Here we present a rapid increase in ultraviolet-optical emission in the centre of a nearby galaxy, marking the onset of sudden increased accretion onto a SMBH. The optical spectrum of this flare, dubbed AT 2017bgt, exhibits a mix of emission features. Some are typical of luminous, unobscured AGNs, but others are likely driven by Bowen fluorescence - robustly linked here with high-velocity gas in the vicinity of the accreting SMBH. The spectral features and increased ultraviolet flux show little evolution over a period of at least 14 months. This disfavours the tidal disruption of a star as their origin, and instead suggests a longer-term event of intensified accretion. Together with two other recently reported events with similar properties, we define a new class of SMBH-related flares. This has important implications for the classification of different types of enhanced accretion onto SMBHs.

We present a broadband X-ray spectral analysis of the M51 system, including the dual active galactic nuclei (AGN) and several off-nuclear point sources. Using a deep observation by NuSTAR, new high-resolution coverage of M51b by Chandra, and the latest X-ray torus models, we measure the intrinsic X-ray luminosities of the AGN in these galaxies. The AGN of M51a is found to be Compton thick, and both AGN have very low accretion rates ($\lambda_{\rm Edd} <10^{-4}$). The latter is surprising considering that the galaxies of M51 are in the process of merging, which is generally predicted to enhance nuclear activity. We find that the covering factor of the obscuring material in M51a is $0.26 \pm 0.03$, consistent with the local AGN obscured fraction at $L_{\rm X}\sim 10^{40}$ erg s$^{-1}$. The substantial obscuring column does not support theories that the torus, presumed responsible for the obscuration, disappears at these low accretion luminosities. However, the obscuration may have resulted from the gas infall driven by the merger rather than the accretion process. We report on several extra-nuclear sources with $L_{\rm X}>10^{39}$ erg s$^{-1}$ and find that a spectral turnover is present below 10 keV in most such sources, in line with recent results on ultraluminous X-ray sources.

The NuSTAR observatory, with its high sensitivity in hard X-rays, has enabled detailed broadband modeling of the X-ray spectra of Active Galactic Nuclei (AGN), thereby allowing constraints to be placed on the high-energy cutoff of the X-ray coronal continuum. We investigate the spectral properties of a sample of 46 NuSTAR-observed Seyfert 1 AGN selected from the Swift/BAT 70-month hard X-ray survey. Our measurements of the high-energy cutoff of the continuum from modeling the NuSTAR X-ray spectra are used to map out the temperature - compactness ($\theta-l$) plane for AGN coronae. We find that most of the coronae lie clustered near the boundary for runaway pair production, suggesting that annihilation and pair production act to regulate the temperature of the corona. We discuss the implications of coronae whose high-energy cutoff may indicate a low coronal temperature on the heating and thermalization mechanisms in the corona.

We report a new changing-look quasar, WISE~J105203.55+151929.5 at $z=0.303$, found by identifying highly mid-IR variable quasars in the WISE/NEOWISE data stream. Compared to multi-epoch mid-IR photometry of a large sample of SDSS-confirmed quasars, WISE J1052+1519 is an extreme photometric outlier, fading by more than a factor of two at $3.4$ and $4.6 \mu$m since 2009. Swift target-of-opportunity observations in 2017 show even stronger fading in the soft X-rays compared to the ROSAT detection of this source in 1995, with at least a factor of fifteen decrease. We obtained second-epoch spectroscopy with the Palomar telescope in 2017 which, when compared with the 2006 archival SDSS spectrum, reveals that the broad H$\beta$ emission has vanished and that the quasar has become significantly redder. The two most likely interpretations for this dramatic change are source fading or obscuration, where the latter is strongly disfavored by the mid-IR data. We discuss various physical scenarios that could cause such changes in the quasar luminosity over this timescale, and favor changes in the innermost regions of the accretion disk that occur on the thermal and heating/cooling front timescales. We discuss possible physical triggers that could cause these changes, and predict the multiwavelength signatures that could distinguish these physical scenarios.

The basic unified model of active galactic nuclei (AGN) invokes an anisotropic obscuring structure, usually referred to as a torus, to explain AGN obscuration as an angle-dependent effect. We present a new grid of X-ray spectral templates based on radiative transfer calculations in neutral gas in an approximately toroidal geometry, appropriate for CCD-resolution X-ray spectra (FWHM > 130 eV). Fitting the templates to broadband X-ray spectra of AGN provides constraints on two important geometrical parameters of the gas distribution around the supermassive black hole: the average column density and the covering factor. Compared to the currently available spectral templates, our model is more flexible, and capable of providing constraints on the main torus parameters in a wider range of AGN. We demonstrate the application of this model using hard X-ray spectra from NuSTAR (3-79 keV) for four AGN covering a variety of classifications: 3C 390.3, NGC 2110, IC 5063, and NGC 7582. This small set of examples was chosen to illustrate the range of possible torus configurations, from disk-like to sphere-like geometries with column densities below, as well as above, the Compton-thick threshold. This diversity of torus properties challenges the simple assumption of a standard geometrically and optically thick toroidal structure commonly invoked in the basic form of the unified model of AGN. Finding broad consistency between the our constraints and those from infrared modeling, we discuss how the approach from the X-ray band complements similar measurements of AGN structures at other wavelengths.

We use results of our 3 year polarimetric monitoring program to investigate the previously suggested connection between rotations of the polarization plane in the optical emission of blazars and their gamma-ray flares in the GeV band. The homogeneous set of 40 rotation events in 24 sources detected by {\em RoboPol} is analysed together with the gamma-ray data provided by {\em Fermi}-LAT. We confirm that polarization plane rotations are indeed related to the closest gamma-ray flares in blazars and the time lags between these events are consistent with zero. Amplitudes of the rotations are anticorrelated with amplitudes of the gamma-ray flares. This is presumably caused by higher relativistic boosting (higher Doppler factors) in blazars that exhibit smaller amplitude polarization plane rotations. Moreover, the time scales of rotations and flares are marginally correlated.

We present X-ray bolometric correction factors, $\kappa_{Bol}$ ($\equiv L_{Bol}/L_X$), for Compton-thick (CT) active galactic nuclei (AGN) with the aim of testing AGN torus models, probing orientation effects, and estimating the bolometric output of the most obscured AGN. We adopt bolometric luminosities, $L_{Bol}$, from literature infrared (IR) torus modeling and compile published intrinsic 2--10 keV X-ray luminosities, $L_{X}$, from X-ray torus modeling of NuSTAR data. Our sample consists of 10 local CT AGN where both of these estimates are available. We test for systematic differences in $\kappa_{Bol}$ values produced when using two widely used IR torus models and two widely used X-ray torus models, finding consistency within the uncertainties. We find that the mean $\kappa_{Bol}$ of our sample in the range $L_{Bol}\approx10^{42}-10^{45}$ erg/s is log$_{10}\kappa_{Bol}=1.44\pm0.12$ with an intrinsic scatter of $\sim0.2$ dex, and that our derived $\kappa_{Bol}$ values are consistent with previously established relationships between $\kappa_{Bol}$ and $L_{Bol}$ and $\kappa_{Bol}$ and Eddington ratio. We investigate if $\kappa_{Bol}$ is dependent on $N_H$ by comparing our results on CT AGN to published results on less-obscured AGN, finding no significant dependence. Since many of our sample are megamaser AGN, known to be viewed edge-on, and furthermore under the assumptions of AGN unification whereby unobscured AGN are viewed face-on, our result implies that the X-ray emitting corona is not strongly anisotropic. Finally, we present $\kappa_{Bol}$ values for CT AGN identified in X-ray surveys as a function of their observed $L_X$, where an estimate of their intrinsic $L_{X}$ is not available, and redshift, useful for estimating the bolometric output of the most obscured AGN across cosmic time.

We use new X-ray data obtained with the Nuclear Spectroscopic Telescope Array (NuSTAR), near-infrared (NIR) fluxes, and mid-infrared (MIR) spectra of a sample of 24 unobscured type 1 active galactic nuclei (AGN) to study the correlation between various hard X-ray bands between 3 and 80 keV and the infrared (IR) emission. The IR to X-ray correlation spectrum (IRXCS) shows a maximum at ~15-20 micron, coincident with the peak of the AGN contribution to the MIR spectra of the majority of the sample. There is also a NIR correlation peak at ~2 micron, which we associate with the NIR bump observed in some type 1 AGN at ~1-5 micron and is likely produced by nuclear hot dust emission. The IRXCS shows practically the same behaviour in all the X-ray bands considered, indicating a common origin for all of them. We finally evaluated correlations between the X-ray luminosities and various MIR emission lines. All the lines show a good correlation with the hard X-rays (rho>0.7), but we do not find the expected correlation between their ionization potentials and the strength of the IRXCS.

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.

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.

Aims. We study the long-term variability of the well-known Seyfert 2 galaxy Mrk 1210 (a.k.a. UGC 4203, or the Phoenix galaxy). Methods. The source was observed by many X-ray facilities in the last 20 years. Here we present a NuSTAR observation and put the results in context of previously published observations. Results. NuSTAR observed Mrk 1210 in 2012 for 15.4 ks. The source showed Compton-thin obscuration similar to that observed by Chandra, Suzaku, BeppoSAX and XMM-Newton over the past two decades, but different from the first observation by ASCA in 1995, in which the active nucleus was caught in a low flux state - or obscured by Compton-thick matter, with a reflection-dominated spectrum. Thanks to the high-quality hard X-ray spectrum obtained with NuSTAR and exploiting the long-term spectral coverage spanning 16.9 years, we can precisely disentangle the transmission and reflection components and put constraints on both the intrinsic long-term variability and hidden nucleus scenarios. In the former case, the distance between the reflector and the source must be at least ~ 2 pc, while in the latter one the eclipsing cloud may be identified with a water maser-emitting clump.

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.

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.

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, {\lambda}$_{Edd}$, and the X-ray power-law index, {\Gamma}, 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 {\Gamma} 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 {\Gamma} values and X-ray luminosities, L$_X$. With a simple bolometric correction to L$_X$ to calculate {\lambda}$_{Edd}$, we find a statistically significant correlation between {\Gamma} and {\lambda}$_{Edd}$ (p = 0.007). A linear fit to the data yields {\Gamma} = (0.41$\pm$0.18)log$_{10}${\lambda}$_{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 {\Gamma}-{\lambda}$_{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 {\Gamma} as a growth-rate indicator for accreting black holes, even for Compton-thick AGN.

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.

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.

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.

We present coordinated multiwavelength observations of the bright, nearby BL Lac object Mrk 421 taken in 2013 January-March, involving GASP-WEBT, Swift, NuSTAR, Fermi-LAT, MAGIC, VERITAS, and other collaborations and instruments, providing data from radio to very-high-energy (VHE) gamma-ray bands. NuSTAR yielded previously unattainable sensitivity in the 3-79 keV range, revealing that the spectrum softens when the source is dimmer until the X-ray spectral shape saturates into a steep power law with a photon index of approximately 3, with no evidence for an exponential cutoff or additional hard components up to about 80 keV. For the first time, we observed both the synchrotron and the inverse-Compton peaks of the spectral energy distribution (SED) simultaneously shifted to frequencies below the typical quiescent state by an order of magnitude. The fractional variability as a function of photon energy shows a double-bump structure which relates to the two bumps of the broadband SED. In each bump, the variability increases with energy which, in the framework of the synchrotron self-Compton model, implies that the electrons with higher energies are more variable. The measured multi-band variability, the significant X-ray-to-VHE correlation down to some of the lowest fluxes ever observed in both bands, the lack of correlation between optical/UV and X-ray flux, the low degree of polarization and its significant (random) variations, the short estimated electron cooling time, and the significantly longer variability timescale observed in the NuSTAR light curves point toward in-situ electron acceleration, and suggest that there are multiple compact regions contributing to the broadband emission of Mrk 421 during low-activity states.

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.

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

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

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.

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 $\gamma$-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 $\gamma$-ray spectrum during one of the flares was very hard, with an index of $\Gamma_\gamma = 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_\gamma/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 $\Delta\Gamma_{\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 $\gamma$-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 $\gamma$-ray and optical spectra cannot satisfactorily explain the production of X-ray emission.

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.

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 $\tau_e=1.23\pm0.08$ assuming a slab (disk-like) geometry, and $kT_e=25\pm2$ keV and $\tau_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.

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$\alpha$ line on time scales of years. The Fe K$\alpha$ 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).

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.

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 'gold standard' 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.

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.

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.

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'\times 13'$ field of view. The control system combines a dynamic scheduler, real-time data reduction, and telescope automation to allow high-efficiency unassisted observations.

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.