Search

Multiwavelength studies indicate that nuclear activity and bulge properties are closely related, but the details remain unclear. To study this further, we combine $Hubble~Space~Telescope$ bulge structural and photometric properties with 1.5 GHz, $e$-MERLIN nuclear radio continuum data from the LeMMINGs survey for a large sample of 173 `active' galaxies (LINERs and Seyferts) and `inactive' galaxies (H IIs and absorption line galaxies, ALGs). Dividing our sample into active and inactive, they define distinct (radio core luminosity)$-$(bulge mass), L_R,core-M_*,bulge, relations, with a mass turnover at M_*, bulge ~ 10^(9.8 +- 0.3) M_sun (supermassive black hole mass M_BH ~ 10^(6.8 +- 0.3) M_sun), which marks the transition from AGN-dominated nuclear radio emission in more massive bulges to that mainly driven by stellar processes in low-mass bulges. None of our 10/173 bulgeless galaxies host an AGN. The AGN fraction increases with increasing M_*, bulge such that f_optical_AGN $\propto$ M_*,bulge^(0.24 +- 0.06) and f_radio_AGN $\propto$ M_*,bulge^(0.24 +- 0.05). Between M_*,bulge ~ 10^8.5 and 10^11.3 M_sun, f_optical_AGN steadily rises from 15 +- 4 to 80 +- 5 per cent. We find that at fixed bulge mass, the radio loudness, nuclear radio activity and the (optical and radio) AGN fraction exhibit no dependence on environment. Radio-loud hosts preferentially possess an early-type morphology than radio-quiet hosts, the two types are however indistinguishable in terms of bulge S\'ersic index and ellipticity, while results on the bulge inner logarithmic profile slope are inconclusive. We finally discuss the importance of bulge mass in determining the AGN triggering processes, including potential implications for the nuclear radio emission in nearby galaxies., Comment: 27 pages, 15 figures, 4 tables, accepted for publication in MNRAS

We study accretion disk-corona connection in Seyfert 1 galaxies using simultaneous UV/X-ray observations of NGC 4593 (July 14-18, 2016) and NGC 7469 (October 15-19, 2017) performed with AstroSat. We use the X-ray (0.5-7.0 keV) data acquired with the Soft X-ray Telescope (SXT) and the UV (FUV: 130-180 nm, NUV: 200-300 nm) data obtained with the Ultra-Violet Imaging Telescope (UVIT). We also use the contemporaneous Swift observations of NGC 4593 and demonstrate AstroSat's capability for X-ray/UV correlation studies. We performed UV/X-ray cross-correlation analysis using the Interpolated and the Discrete Cross-Correlation Functions and found similar results. In the case of NGC 4593, we found that the variations in the X-rays lead to those in the FUV and NUV bands by ~ 38 ks and ~ 44 ks, respectively. These UV lags favour the disk reprocessing model, they are consistent with the previous results within uncertainties. In contrast, we found an opposite trend in NGC 7469 where the soft X-ray variations lag those in the FUV and NUV bands by ~ 41 ks and ~ 49 ks, respectively. The hard lags in NGC 7469 favour the Thermal Comptonization model. Our results may provide direct observational evidence for the variable intrinsic UV emission from the accretion disk which acts as the seed for thermal Comptonization in a hot corona in a lamp-post like geometry. The non-detection of disk reverberation photons in NGC 7469, using AstroSat data, is most likely due to a high accretion rate resulting in a hot accretion disk and large intrinsic emission., 14 pages, 20 figures, Accepted for publication in MNRAS

We present fast (~200s sampling) ugriz photometry of the low mass AGN NGC 4395 with the Liverpool Telescope, followed by very fast (3s sampling) us, gs, rs, is and zs simultaneous monitoring with HiPERCAM on the 10.4m GTC. These observations provide the fastest ever AGN multiband photometry and very precise lag measurements. Unlike in all other AGN, gs lags us by a large amount, consistent with disc reprocessing but not with reprocessing in the Broad Line Region (BLR). There is very little increase in lag with wavelength at long wavelengths, indicating an outer edge (Rout) to the reprocessor. We have compared truncated disc reprocessing models to the combined HiPERCAM and previous X-ray/UV lags. For the normally accepted mass of 3.6E5 solar, we obtain reasonable agreement with zero spin, Rout ~1700 Rg, and the DONE physically-motivated temperature-dependent disc colour correction factor (fcol). A smaller mass of 4E4 solar can only be accomodated if fcol=2.4, which is probably unrealistically high. Disc self gravity is probably unimportant in this low mass AGN but an obscuring wind may provide an edge. For the small mass the dust sublimation radius is similar to Rout, so the wind could be dusty. However for the more likely large mass the sublimation radius is further out so the optically-thick base of a line-driven gaseous wind is more likely. The inner edge of the BLR is close to Rout in both cases. These observations provide the first good evidence for a truncated AGN disc and caution that truncation should be included in reverberation lag modelling., 19 pages, 27 figures

Full list of authors: Baldi, R. D.; Williams, D. R. A.; McHardy, I. M.; Beswick, R. J.; Brinks, E.; Dullo, B. T.; Knapen, J. H.; Argo, M. K.; Aalto, S.; Alberdi, A.; Baan, W. A.; Bendo, G. J.; Corbel, S.; Fenech, D. M.; Gallagher, J. S.; Green, D. A.; Kennicutt, R. C.; Klöckner, H. -R.; Körding, E.; Maccarone, T. J.; Muxlow, T. W. B.; Mundell, C. G.; Panessa, F.; Peck, A. B.; Pérez-Torres, M. A.; Romero-Cañizales, C.; Saikia, P.; Shankar, F.; Spencer, R. E.; Stevens, I. R.; Varenius, E.; Ward, M. J.; Yates, J.; Uttley, P., We present the second data release of high-resolution (≤0.2 arcsec) 1.5-GHz radio images of 177 nearby galaxies from the Palomar sample, observed with the e-MERLIN array, as part of the Legacy e-MERLIN Multi-band Imaging of Nearby Galaxies Sample (LeMMINGs) survey. Together with the 103 targets of the first LeMMINGs data release, this represents a complete sample of 280 local active (LINER and Seyfert) and inactive galaxies (H ii galaxies and absorption line galaxies, ALG). This large program is the deepest radio survey of the local Universe, ≳1017.6 W Hz-1, regardless of the host and nuclear type: we detect radio emission ≳0.25 mJy beam-1 for 125/280 galaxies (44.6 per cent) with sizes of typically ≲100 pc. Of those 125, 106 targets show a core which coincides within 1.2 arcsec with the optical nucleus. Although we observed mostly cores, around one third of the detected galaxies features jetted morphologies. The detected radio core luminosities of the sample range between ∼1034 and 1040 erg s-1. LINERs and Seyferts are the most luminous sources, whereas H ii galaxies are the least. LINERs show FR I-like core-brightened radio structures while Seyferts reveal the highest fraction of symmetric morphologies. The majority of H ii galaxies have single radio core or complex extended structures, which probably conceal a nuclear starburst and/or a weak active nucleus (seven of them show clear jets). ALGs, which are typically found in evolved ellipticals, although the least numerous, exhibit on average the most luminous radio structures, similar to LINERs. © 2020 The Author(s). Published by Oxford University Press on behalf of Royal Astronomical Society., The authors thank the anonymous referee for his/her helpful comments to improve the manuscript. AA and MAPT acknowledge support from the Spanish MCIU through grant PGC2018-098915-B-C21 and from the State Agency for Research of the Spanish MCIU through the 'Center of Excellence Severo Ochoa' award for the Instituto de Astrof ' isica de Andalucia (SEV-2017-0709). BTD acknowledges support from a Spanish postdoctoral fellowship 'Ayudas 1265 para la atraccion del talento investigador. Modalidad 2: jovenes investigadores.' funded by Comunidad de Madrid under grant number 2016-T2/TIC-2039. BTD also acknowledges support from grant 'Ayudas para la realizaci on de proyectos de I + D para jovenes doctores 2019.' funded by Comunidad de Madrid and Universidad Complutense de Madrid under grant number PR65/19-22417. JHK acknowledges financial support from the European Union's Horizon 2020 research and innovation programme under Marie Sklodowska-Curie grant agreement No. 721463 to the SUNDIAL ITN network, from the State Research Agency (AEI-MCINN) of the Spanish Ministry of Science and Innovation under the grant 'The structure and evolution of galaxies and their central regions' with reference PID2019-105602GB-I00/10.13039/501100011033, and from IAC project P/300724, financed by the Ministry of Science and Innovation, through the State Budget and by the Canary Islands Department of Economy, Knowledge and Employment, through the Regional Budget of the Autonomous Community. JSG thanks the University of Wisconsin-Madison and its Foundation for support of this research through his Rupple Bascom Professorship. FS acknowledges partial support from a Leverhulme Trust Research fellowship. CGM acknowledges support from the University of Bath and Jim and Hiroko Sherwin. e-MERLIN is a National Facility operated by the University of Manchester at Jodrell Bank Observatory on behalf of STFC, part of UK Research and Innovation.

We present the long term, frequency-dependent, X-ray time lags of the Seyfert galaxy NGC 4051. We used Rossi X-ray Timing Explorer light curves in the 2-4, 4-7 and 7-10 keV bands and we measured the time-lags at $10^{-7}-10^{-6}$ Hz. This is the lowest frequency range that AGN time-lags have ever been measured. The variations in the higher energy bands are delayed with respect to the variations we observe at 2-4 keV, in agreement with the time-lags at high frequencies. When we combine our results with those from the model fitting of the time lags at higher frequencies we find that that the X-ray hard lags in NGC 4051 follow a power law of slope $\sim -1$ over a broad frequency range, from $\sim 10^{-7}$ to $\sim 10^{-3}$ Hz. This is consistent with the time-lags of Cyg X-1, a result which further supports the analogy between active galaxies and Galactic X-ray black hole binaries., Accepted for publication in MNRAS

We present the first intensive continuum reverberation mapping study of the high accretion rate Seyfert galaxy Mrk 110. The source was monitored almost daily for more than 200 days with the Swift X-ray and UV/optical telescopes, supported by ground-based observations from Las Cumbres Observatory, the Liverpool Telescope, and the Zowada Observatory, thus extending the wavelength coverage to 9100 \r{A}. Mrk 110 was found to be significantly variable at all wavebands. Analysis of the intraband lags reveals two different behaviours, depending on the timescale. On timescales shorter than 10 days the lags, relative to the shortest UV waveband ($\sim1928$ \r{A}), increase with increasing wavelength up to a maximum of $\sim2$days lag for the longest waveband ($\sim9100$ \r{A}), consistent with the expectation from disc reverberation. On longer timescales, however, the g-band lags the Swift BAT hard X-rays by $\sim10$ days, with the z-band lagging the g-band by a similar amount, which cannot be explained in terms of simple reprocessing from the accretion disc. We interpret this result as an interplay between the emission from the accretion disc and diffuse continuum radiation from the broad line region., Comment: 19 pages, 11 figures. Accepted for publication in MNRAS

We report velocity-delay maps for prominent broad emission lines, Ly_alpha, CIV, HeII and H_beta, in the spectrum of NGC5548. The emission-line responses inhabit the interior of a virial envelope. The velocity-delay maps reveal stratified ionization structure. The HeII response inside 5-10 light-days has a broad single-peaked velocity profile. The Ly_alpha, CIV, and H_beta responses peak inside 10 light-days, extend outside 20 light-days, and exhibit a velocity profile with two peaks separated by 5000 km/s in the 10 to 20 light-day delay range. The velocity-delay maps show that the M-shaped lag vs velocity structure found in previous cross-correlation analysis is the signature of a Keplerian disk with a well-defined outer edge at R=20 light-days. The outer wings of the M arise from the virial envelope, and the U-shaped interior of the M is the lower half of an ellipse in the velocity-delay plane. The far-side response is weaker than that from the near side, so that we see clearly the lower half, but only faintly the upper half, of the velocity--delay ellipse. The delay tau=(R/c)(1-sin(i))=5 light-days at line center is from the near edge of the inclined ring, giving the inclination i=45 deg. A black hole mass of M=7x10^7 Msun is consistent with the velocity-delay structure. A barber-pole pattern with stripes moving from red to blue across the CIV and possibly Ly_alpha line profiles suggests the presence of azimuthal structure rotating around the far side of the broad-line region and may be the signature of precession or orbital motion of structures in the inner disk. Further HST observations of NGC 5548 over a multi-year timespan but with a cadence of perhaps 10 days rather than 1 day could help to clarify the nature of this new AGN phenomenon., 19 pages, 9 figures, ApJ in press

Full list of authors: Baldi, R. D.; Williams, D. R. A.; Beswick, R. J.; McHardy, I.; Dullo, B. T.; Knapen, J. H.; Zanisi, L.; Argo, M. K.; Aalto, S.; Alberdi, A.; Baan, W. A.; Bendo, G. J.; Fenech, D. M.; Green, D. A.; Klöckner, H. -R.; Körding, E.; Maccarone, T. J.; Marcaide, J. M.; Mutie, I.; Panessa, F.; Pérez-Torres, M. A.; Romero-Cañizales, C.; Saikia, D. J.; Saikia, P.; Shankar, F.; Spencer, R. E.; Stevens, I. R.; Uttley, P.; Brinks, E.; Corbel, S.; Martí-Vidal, I.; Mundell, C. G.; Pahari, M.; Ward, M. J., What determines the nuclear radio emission in local galaxies? To address this question, we combine optical [O iii] line emission, robust black hole (BH) mass estimates, and high-resolution e-MERLIN 1.5-GHz data, from the LeMMINGs survey, of a statistically complete sample of 280 nearby optically active (LINER and Seyfert) and inactive [H ii and absorption line galaxies (ALGs)] galaxies. Using [O iii] luminosity (L[O III]) as a proxy for the accretion power, local galaxies follow distinct sequences in the optical-radio planes of BH activity, which suggest different origins of the nuclear radio emission for the optical classes. The 1.5-GHz radio luminosity of their parsec-scale cores (Lcore) is found to scale with BH mass (MBH) and [O iii] luminosity. Below MBH ∼106.5 M⊙, stellar processes from non-jetted H ii galaxies dominate with Lcore ∝ MBH0.61 ± 0.33 and Lcore ∝ L[O III]0.79 ± 0.30. Above MBH ∼106.5 M⊙, accretion-driven processes dominate with Lcore ∝ MBH1.5-1.65 and Lcore ∝ L[O III]0.99-1.31 for active galaxies: radio-quiet/loud LINERs, Seyferts, and jetted H ii galaxies always display (although low) signatures of radio-emitting BH activity, with L1.5 GHz ≳ 1019.8 W Hz-1 and MBH ≳ 107 M⊙, on a broad range of Eddington-scaled accretion rates (m). Radio-quiet and radio-loud LINERs are powered by low-m discs launching sub-relativistic and relativistic jets, respectively. Low-power slow jets and disc/corona winds from moderately high to high-m discs account for the compact and edge-brightened jets of Seyferts, respectively. Jetted H ii galaxies may host weakly active BHs. Fuel-starved BHs and recurrent activity account for ALG properties. In conclusion, specific accretion-ejection states of active BHs determine the radio production and the optical classification of local active galaxies. © 2021 The Author(s) Published by Oxford University Press on behalf of Royal Astronomical Society., This research was supported by a Newton Fund project, DARA (Development in Africa with Radio Astronomy), and awarded by the UK’s Science and Technology Facilities Council (STFC) – grant reference ST/R001103/1. This research was supported by European Commission Horizon 2020 Research and Innovation Programme under grant agreement No. 730884 (JUMPING JIVE). FS acknowledges support from a Leverhulme Trust Research Fellowship. JHK acknowledges support from the European Union’s Horizon 2020 research and innovation programme under Marie Skłodowska-Curie grant agreement no. 721463 to the SUNDIAL ITN network, from the State Research Agency (AEI) of the Spanish Ministry of Science and Innovation and the European Regional Development Fund (FEDER) under the grant with reference PID2019-105602GB-I00, and from IAC project P/300724, financed by the Ministry of Science and Innovation, through the State Budget and by the Canary Islands Department of Economy, Knowledge and Employment, through the Regional Budget of the Autonomous Community. IMV thanks the Generalitat Valenciana (funding from the GenT Project CIDEGENT/2018/021) and the MICINN (funding from the Research Project PID2019-108995GB-C22). AA and MAPT acknowledge support from the Spanish MCIU through grant PGC2018-098915-B-C21 and from the State Agency for Research of the Spanish MCIU through the ‘Center of Excellence Severo Ochoa’ award for the Instituto de Astrofìsica de Andalucìa (SEV-2017-0709). BTD acknowledges financial support from grant ‘Ayudas para la realización de proyectos de I + D para jóvenes doctores 2019’ funded by Comunidad de Madrid and Universidad Complutense de Madrid under grant no. PR65/19-22417.

We present geometric and dynamical modeling of the broad line region for the multi-wavelength reverberation mapping campaign focused on NGC 5548 in 2014. The dataset includes photometric and spectroscopic monitoring in the optical and ultraviolet, covering the H$\beta$, C IV, and Ly$\alpha$ broad emission lines. We find an extended disk-like H$\beta$ BLR with a mixture of near-circular and outflowing gas trajectories, while the C IV and Ly$\alpha$ BLRs are much less extended and resemble shell-like structures. There is clear radial structure in the BLR, with C IV and Ly$\alpha$ emission arising at smaller radii than the H$\beta$ emission. Using the three lines, we make three independent black hole mass measurements, all of which are consistent. Combining these results gives a joint inference of $\log_{10}(M_{\rm BH}/M_\odot) = 7.64^{+0.21}_{-0.18}$. We examine the effect of using the $V$ band instead of the UV continuum light curve on the results and find a size difference that is consistent with the measured UV-optical time lag, but the other structural and kinematic parameters remain unchanged, suggesting that the $V$ band is a suitable proxy for the ionizing continuum when exploring the BLR structure and kinematics. Finally, we compare the H$\beta$ results to similar models of data obtained in 2008 when the AGN was at a lower luminosity state. We find that the size of the emitting region increased during this time period, but the geometry and black hole mass remain unchanged, which confirms that the BLR kinematics suitably gauge the gravitational field of the central black hole., Comment: 26 pages, 19 figures, 1 table, accepted for publication in ApJ

In 2014 the NGC 5548 Space Telescope and Optical Reverberation Mapping campaign discovered a two-month anomaly when variations in the absorption and emission lines decorrelated from continuum variations. During this time the soft X-ray part of the intrinsic spectrum had been strongly absorbed by a line-of-sight (LOS) obscurer, which was interpreted as the upper part of a disk wind. Our first paper showed that changes in the LOS obscurer produce the decorrelation between the absorption lines and the continuum. A second study showed that the base of the wind shields the BLR, leading to the emission-line decorrelation. In that study, we proposed the wind is normally transparent with no effect on the spectrum. Changes in the wind properties alter its shielding and affect the SED striking the BLR, producing the observed decorrelations. In this work, we investigate the impact of a translucent wind on the emission lines. We simulate the obscuration using XMM-Newton, NuSTAR, and HST observations to determine the physical characteristics of the wind. We find that a translucent wind can contribute a part of the He II and Fe K? emission. It has a modest optical depth to electron scattering, which explains the fainter far-side emission in the observed velocity delay maps. The wind produces the very broad base seen in the UV emission lines and may also be present in the Fe K? line. Our results highlight the importance of accounting for the effects of such winds in the analysis of the physics of the central engine., Comment: Accepted for publication in ApJ

Submillimetre-selected galaxies (SMGs) at high redshift ($z$ $\sim$ 2) are potential host galaxies of active galactic nuclei (AGN). If the local Universe is a good guide, $\sim$ 50$\%$ of the obscured AGN amongst the SMG population could be missed even in the deepest X-ray surveys. Radio observations are insensitive to obscuration; therefore, very long baseline interferometry (VLBI) can be used as a tool to identify AGN in obscured systems. A well-established upper limit to the brightness temperature of 10$^5$ K exists in star-forming systems, thus VLBI observations can distinguish AGN from star-forming systems via brightness temperature measurements. We present 1.6 GHz European VLBI Network (EVN) observations of four SMGs (with measured redshifts) to search for evidence of compact radio components associated with AGN cores. For two of the sources, e-MERLIN images are also presented. Out of the four SMGs observed, we detect one source, J123555.14, that has an integrated EVN flux density of 201 $\pm$ 15.2 $\mu$Jy, corresponding to a brightness temperature of 5.2 $\pm$ 0.7 $\times$ 10$^5$ K. We therefore identify that the radio emission from J123555.14 is associated with an AGN. We do not detect compact radio emission from a possible AGN in the remaining sources (J123600.10, J131225.73, and J163650.43). In the case of J131225.73, this is particularly surprising, and the data suggest that this may be an extended, jet-dominated AGN that is resolved by VLBI. Since the morphology of the faint radio source population is still largely unknown at these scales, it is possible that with a $\sim$ 10 mas resolution, VLBI misses (or resolves) many radio AGN extended on kiloparsec scales., Comment: 11 pages, 5 figures

Full list of authors: Muxlow, T. W. B.; Thomson, A. P.; Radcliffe, J. F.; Wrigley, N. H.; Beswick, R. J.; Smail, Ian; McHardy, I. M.; Garrington, S. T.; Ivison, R. J.; Jarvis, M. J.; Prandoni, I.; Bondi, M.; Guidetti, D.; Argo, M. K.; Bacon, David; Best, P. N.; Biggs, A. D.; Chapman, S. C.; Coppin, K.; Chen, H.; Garratt, T. K.; Garrett, M. A.; Ibar, E.; Kneib, Jean-Paul; Knudsen, Kirsten K.; Koopmans, L. V. E.; Morabito, L. K.; Murphy, E. J.; Njeri, A.; Pearson, Chris; Pérez-Torres, M. A.; Richards, A. M. S.; Röttgering, H. J. A.; Sargent, M. T.; Serjeant, Stephen; Simpson, C.; Simpson, J. M.; Swinbank, A. M.; Varenius, E.; Venturi, T., We present an overview and description of the e-MERGE Survey (e-MERLIN Galaxy Evolution Survey) Data Release 1 (DR1), a large program of high-resolution 1.5-GHz radio observations of the GOODS-N field comprising similar to 140 h of observations with enhanced-Multi-Element Remotely Linked Interferometer Network (e-MERLIN) and similar to 40 h with the Very Large Array (VLA). We combine the long baselines of e-MERLIN (providing high angular resolution) with the relatively closely packed antennas of the VLA (providing excellent surface brightness sensitivity) to produce a deep 1.5-GHz radio survey with the sensitivity (similar to 1.5 mu Jy beam(-1)), angular resolution (0.2-0.7 arcsec) and field-of-view (similar to 15x15 arcmin(2)) to detect and spatially resolve star-forming galaxies and active galactic nucleus (AGN) at z greater than or similar to 1. The goal of e-MERGE is to provide new constraints on the deep, sub-arcsecond radio sky which will be surveyed by SKA1-mid. In this initial publication, we discuss our data analysis techniques, including steps taken to model in-beam source variability over an similar to 20-yr baseline and the development of newpoint spread function/primary beam models to seamlessly merge e-MERLIN and VLA data in the uv plane. We present early science results, including measurements of the luminosities and/or linear sizes of similar to 500 galaxies selected at 1.5 GHz. In combination with deep Hubble Space Telescope observations, we measure a mean radio-to-optical size ratio of r(e-MERGE)/r(HST) similar to 1.02 +/- 0.03, suggesting that in most high-redshift galaxies, the similar to GHz continuum emission traces the stellar light seen in optical imaging. This is the first in a series of papers that will explore the similar to kpc-scale radio properties of star-forming galaxies and AGN in the GOODS-N field observed by e-MERGE DR1. © 2020 The Author(s)., APT, TWBM, RJB, and MAG acknowledge support from STFC (ST/P000649/1). APT further acknowledges and thanks Robert Dickson, Peter Draper, David Hempston, Anthony Holloway, and Alan Lotts for their extensive support managing the highperformance computing infrastructure, which has been essential for the successful delivery of e-MERGE DR1. JFR acknowledges the financial assistance of the South African Radio Astronomy Observatory (SARAO) towards this research (http://sarao.ac.za).IP acknowledges support from INAF under the PRIN SKA/CTA project 'FORECaST' and the PRIN MAIN STREAM project 'SAuROS'. EI acknowledges partial support from FONDECYT through grant no 1171710. We thank the anonymous reviewer for their detailed and constructive referee report, which has undoubtedly improved the quality of the final manuscript. e-MERLIN is a National Facility operated by the University of Manchester at Jodrell Bank Observatory on behalf of STFC. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. The authors acknowledge the use of the IRIDIS High Performance Computing Facility, and associated support services at the University of Southampton, in the completion of this work. We thank IRIS (www.iris.ac.uk) for provision of highperformance computing facilities. STFC IRIS is investing in the UK's Radio and mm/sub-mm Interferometry Services in order to improve data quality and throughput. The European VLBI Network is a joint facility of independent European, African, Asian, and North American radio astronomy institutes. Scientific results from data presented in this publication are derived from the following EVN project code(s): EG078B. The research leading to these results has received funding from the European Commission Horizon 2020 Research and Innovation Programme under grant agreement no. 730562 (RadioNet)., With funding from the Spanish government through the ‘Severo Ochoa Centre of Excellence’ accreditation SEV-2017-0709.

We present the results of X-ray spectral-timing analysis of a 90ks XMM-Newton observation of the nearby, broad line, early-type AGN NGC5273. The X-ray spectrum revealed the clear presence of a reflection component at high energies, with a clear signature of a narrow iron line at 6.4 keV, consistent with distant reflection. Applying a relativistic reflection model, we found only marginal evidence for a broader relativistic line component. However, cross-spectral analysis revealed that, between 4 and 6 $\times 10^{-4}$ Hz, the 5-8 keV band lagged the 2-3 keV band, implying reflection of the iron line from material close to the black hole. From the analysis of the lag-energy spectrum, we found a broad, but skewed line with a peak of $\approx$ 1000s at 7.5 keV relative to the continuum, which we interpret as the iron line in the reverberation spectrum from an illuminated accretion disc. From the asymmetry in the shape of lag-energy spectrum, we also found that the source is consistent with having an inclination $\geq 45^\circ$., Accepted for publication in MNRAS. 7 pages, 6 figures

We model the ultraviolet spectra of the Seyfert 1 galaxy NGC~5548 obtained with the Hubble Space Telescope during the 6-month reverberation-mapping campaign in 2014. Our model of the emission from NGC 5548 corrects for overlying absorption and deblends the individual emission lines. Using the modeled spectra, we measure the response to continuum variations for the deblended and absorption-corrected individual broad emission lines, the velocity-dependent profiles of Ly$\alpha$ and C IV, and the narrow and broad intrinsic absorption features. We find that the time lags for the corrected emission lines are comparable to those for the original data. The velocity-binned lag profiles of Ly$\alpha$ and C IV have a double-peaked structure indicative of a truncated Keplerian disk. The narrow absorption lines show delayed response to continuum variations corresponding to recombination in gas with a density of $\sim 10^5~\rm cm^{-3}$. The high-ionization narrow absorption lines decorrelate from continuum variations during the same period as the broad emission lines. Analyzing the response of these absorption lines during this period shows that the ionizing flux is diminished in strength relative to the far-ultraviolet continuum. The broad absorption lines associated with the X-ray obscurer decrease in strength during this same time interval. The appearance of X-ray obscuration in $\sim\,2012$ corresponds with an increase in the luminosity of NGC 5548 following an extended low state. We suggest that the obscurer is a disk wind triggered by the brightening of NGC 5548 following the decrease in size of the broad-line region during the preceding low-luminosity state., Comment: 50 pages, 30 figures, uses aastex62.cls. Accepted for publication in ApJ, 07/06/2019. High-level products page in MAST will go live after 7/15/2019. Replaced Figure 4 on 7/12/2019 to be more red/green color-blind friendly

We broadly explore the effects of systematic errors on reverberation mapping lag uncertainty estimates from {\tt JAVELIN} and the interpolated cross-correlation function (ICCF) method. We focus on simulated lightcurves from random realizations of the lightcurves of five intensively monitored AGNs. Both methods generally work well even in the presence of systematic errors, although {\tt JAVELIN} generally provides better error estimates. Poorly estimated lightcurve uncertainties have less effect on the ICCF method because, unlike {\tt JAVELIN}, it does not explicitly assume Gaussian statistics. Neither method is sensitive to changes in the stochastic process driving the continuum or the transfer function relating the line lightcurve to the continuum. The only systematic error we considered that causes significant problems is if the line lightcurve is not a smoothed and shifted version of the continuum lightcurve but instead contains some additional sources of variability., Comment: 22 pages, 17 figures

We present high-sensitivity eMERLIN radio images of the Seyfert galaxy NGC 4151 at 1.5 GHz. We compare the new eMERLIN images to those from archival MERLIN observations in 1993 to determine the change in jet morphology in the 22 years between observations. We report an increase by almost a factor of 2 in the peak flux density of the central core component, C4, thought to host the black hole, but a probable decrease in some other components, possibly due to adiabatic expansion. The core flux increase indicates an AGN which is currently active and feeding the jet. We detect no significant motion in 22 years between C4 and the component C3, which is unresolved in the eMERLIN image. We present a spectral index image made within the 512 MHz band of the 1.5 GHz observations. The spectrum of the core, C4, is flatter than that of other components further out in the jet. We use HST emission line images (H$\alpha$, [O III] and [O II]) to study the connection between the jet and the emission line region. Based on the changing emission line ratios away from the core and comparison with the eMERLIN radio jet, we conclude that photoionisation from the central AGN is responsible for the observed emission line properties further than 4" (360 pc) from the core, C4. Within this region, several evidences (radio-line co-spatiality, low [O III]/H$\alpha$ and estimated fast shocks) suggest additional ionisation from the jet., Comment: 13 pages, 8 figures, 2 tables. Accepted to MNRAS: MN-17-2603-MJ.R1

We studied the nearby, interacting galaxy NGC 5195 (M51b) in the radio, optical and X-ray bands. We mapped the extended, low-surface-brightness features of its radio-continuum emission; determined the energy content of its complex structure of shock-ionized gas; constrained the current activity level of its supermassive nuclear black hole. In particular, we combined data from the European Very Long Baseline Interferometry Network (~1-pc scale), from our new e-MERLIN observations (~10-pc scale), and from the Very Large Array (~100-1000-pc scale), to obtain a global picture of energy injection in this galaxy. We put an upper limit to the luminosity of the (undetected) flat-spectrum radio core. We find steep-spectrum, extended emission within 10 pc of the nuclear position, consistent with optically-thin synchrotron emission from nuclear star formation or from an outflow powered by an active galactic nucleus (AGN). A linear spur of radio emission juts out of the nuclear source towards the kpc-scale arcs (detected in radio, Halpha and X-ray bands). From the size, shock velocity, and Balmer line luminosity of the kpc-scale bubble, we estimate that it was inflated by a long-term-average mechanical power ~3-6 x 10^{41} erg/s over the last 3-6 Myr. This is an order of magnitude more power than can be provided by the current level of star formation, and by the current accretion power of the supermassive black hole. We argue that a jet-inflated bubble scenario associated with previous episodes of AGN activity is the most likely explanation for the kpc-scale structures., Comment: 16 pages, 7 figures. Accepted by MNRAS on 2017 Feb 12

We perform the NuSTAR and Swift/XRT joint energy spectral fitting of simultaneous observations from the broad-line Seyfert 1.5 galaxy NGC 5273. When fitted with the combination of an exponential cut-off power-law and a reflection model, a high energy cut-off is detected at 143$^{+96}_{-40}$ keV with 2-sigma significance. Existence of such cut-off is also consistent with the observed Comptonizing electron temperature when fitted with a Comptonization model independently. We observe a moderate hard X-ray variability of the source over the time-scale of ~12 years using INTEGRAL/ISGRI observations in the energy range of 20-100 keV. When the hard band count rate (6-20 keV) is plotted against the soft band count rate (3-6 keV), a hard offset is observed. Our results indicate that the cut-off energy may not correlate with the coronal X-ray luminosity in a simple manner. Similarities in parameters that describe coronal properties indicate that the coronal structure of NGC 5273 may be similar to that of the broad-line radio galaxy 3C 390.3 and another galaxy MCG-5-23-16 where the coronal plasma is dominated by electrons, rather than electron-positron pairs. Therefore, the coronal cooling is equally efficient to the heating mechanism keeping the cut-off energy at low even at the low accretion rate., Comment: 11 pages, 6 figures, 1 table, accepted for publication in MNRAS

We present high-resolution e-MERLIN radio continuum maps of the Dwarf Irregular galaxy IC10 at 1.5GHz and 5GHz. We detect 11 compact sources at 1.5GHz, 5 of which have complementary detections at 5GHz. We classify 3 extended sources as compact H{\sc ii} regions within IC10, 5 sources as contaminating background galaxies and identify 3 sources which require additional observations to classify. We do not expect that any of these 3 sources are Supernova Remnants as they will likely be resolved out at the assumed distance of IC10 (0.7Mpc). We correct integrated flux densities of IC10 from the literature for contamination by unrelated background sources and obtain updated flux density measurements of $354\pm11$\,mJy at 1.5GHz and $199\pm9$\,mJy at 4.85GHz. The background contamination does not contribute significantly to the overall radio emission from IC10, so previous analysis concerning its integrated radio properties remain valid., Comment: 15 pages, 7 Figures, 5 Tables, Accepted for publication in MNRAS

We report the results of intensive X-ray, UV and optical monitoring of the Seyfert 1 galaxy NGC 4593 with Swift. There is no intrinsic flux-related spectral change in the the variable components in any band with small apparent variations due only to contamination by a second constant component, possibly a (hard) reflection component in the X-rays and the (red) host galaxy in the UV/optical bands. Relative to the shortest wavelength band, UVW2, the lags of the other UV and optical bands are mostly in agreement with the predictions of reprocessing of high energy emission from an accretion disc. The U-band lag is, however, far larger than expected, almost certainly because of reprocessed Balmer continuum emission from the more distant broad line region gas. The UVW2 band is well correlated with the X-rays but lags by ~6x more than expected if the UVW2 results from reprocessing of X-rays on the accretion disc. However, if the lightcurves are filtered to remove variations on timescales >5d, the lag approaches the expectation from disc reprocessing. MEMEcho analysis shows that direct X-rays can be the driver of most of the variations in the UV/optical bands as long as the response functions for those bands all have long tails (up to 10d) in addition to a strong peak (from disc reprocessing) at short lag (, 18 pages, 17 figures. Submitted to Monthly Notices of the Royal Astronomical Society

We present new observations of the GOODS-N field obtained at 5.5 GHz with the Karl G. Jansky Very Large Array (VLA). The central region of the field was imaged to a median r.m.s. of 3 microJy/beam with a resolution of 0.5 arcsec. From a 14-arcmin diameter region we extracted a sample of 94 radio sources with signal-to-noise ratio greater than 5. Near-IR identifications are available for about 88 percent of the radio sources. We used different multi-band diagnostics to separate active galactic nuclei (AGN), both radiatively efficient and inefficient, from star-forming galaxies. From our analysis, we find that about 80 percent of our radio-selected sample is AGN-dominated, with the fraction raising to 92 percent when considering only the radio sources with redshift >1.5. This large fraction of AGN-dominated radio sources at very low flux densities (the median flux density at 5.5 GHz is 42 microJy), where star-forming galaxies are expected to dominate, is somewhat surprising and at odds with other results. Our interpretation is that both the frequency and angular resolution of our radio observations strongly select against radio sources whose brightness distribution is diffuse on scale of several kpc. Indeed, we find that the median angular sizes of the AGN-dominated sources is around 0.2-0.3 arcsec against 0.8 arcsec for star-forming galaxies. This highlights the key role that high frequency radio observations can play in pinpointing AGN-driven radio emission at microJy levels. This work is part of the eMERGE legacy project., 18 pages, 11 figures. Online material (table 1, table 4 and appendix) available in the source file. Accepted for publication in MNRAS

We present the results of an optical spectroscopic monitoring program targeting NGC 5548 as part of a larger multi-wavelength reverberation mapping campaign. The campaign spanned six months and achieved an almost daily cadence with observations from five ground-based telescopes. The H$\beta$ and He II $\lambda$4686 broad emission-line light curves lag that of the 5100 $\AA$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$\beta$ lag relative to the 1158 $\AA$ ultraviolet continuum light curve measured by the Hubble Space Telescope is roughly $\sim$50% longer than that measured against the optical continuum, and the lag difference is consistent with the observed lag between the optical and ultraviolet continua. This suggests that the characteristic radius of the broad-line region is $\sim$50% larger than the value inferred from optical data alone. We also measured velocity-resolved emission-line lags for H$\beta$ and found a complex velocity-lag structure with shorter lags in the line wings, indicative of a broad-line region dominated by Keplerian motion. The responses of both the H$\beta$ and He II $\lambda$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $\alpha$, He II(+O III]), and Si IV(+O IV]) during the same monitoring period. Finally, given the optical luminosity of NGC 5548 during our campaign, the measured H$\beta$ lag is a factor of five shorter than the expected value implied by the $R_\mathrm{BLR} - L_\mathrm{AGN}$ relation based on the past behavior of NGC 5548., Comment: 22 pages, 13 figures, accepted to ApJ

We present results from a detailed spectral-timing analysis of a long ~486 ks XMM-Newton observation of the bare Seyfert 1 galaxy Ark 120 which showed alternating diminution and increment in the 0.3-10 keV X-ray flux over four consecutive orbits in 2014. We study the energy-dependent variability of Ark 120 through broad-band X-ray spectroscopy, fractional root-mean-squared (rms) spectral modelling, hardness-intensity diagram and flux-flux analysis. The X-ray (0.3-10 keV) spectra are well fitted by a thermally Comptonized primary continuum with two (blurred and distant) reflection components and an optically thick, warm Comptonization component for the soft X-ray excess emission below ~2 keV. During the first and third observations, the fractional X-ray variability amplitude decreases with energy while for second and fourth observations, X-ray variability spectra are found to be inverted-crescent and crescent shaped respectively. The rms variability spectra are well modelled by two constant reflection components, a soft excess component with variable luminosity and a variable intrinsic continuum with the normalization and spectral slope being correlated. The spectral softening of the source with both the soft excess and UV luminosities favour Comptonization models where the soft excess and primary X-ray emission are produced through Compton up-scattering of the UV and UV/soft X-ray seed photons in the putative warm and hot coronae, respectively. Our analyses imply that the observed energy-dependent variability of Ark 120 is most likely due to variations in the spectral shape and luminosity of the hot corona and to variations in the luminosity of the warm corona, both of which are driven by variations in the seed photon flux., Comment: 17 pages, 19 figures, 6 tables, accepted for publication in MNRAS

We present the first systematic physical modelling of the time-lag spectra between the soft (0.3-1 keV) and the hard (1.5-4 keV) X-ray energy bands, as a function of Fourier frequency, in a sample of 12 active galactic nuclei which have been observed by XMM-Newton. We concentrate particularly on the negative X-ray time-lags (typically seen above $10^{-4}$ Hz) i.e. soft band variations lag the hard band variations, and we assume that they are produced by reprocessing and reflection by the accretion disc within a lamp-post X-ray source geometry. We also assume that the response of the accretion disc, in the soft X-ray bands, is adequately described by the response in the neutral iron line (Fe k$\alpha$) at 6.4 keV for which we use fully general relativistic ray-tracing simulations to determine its time evolution. These response functions, and thus the corresponding time-lag spectra, yield much more realistic results than the commonly-used, but erroneous, top-hat models. Additionally we parametrize the positive part of the time-lag spectra (typically seen below $10^{-4}$ Hz) by a power-law. We find that the best-fitting BH masses, M, agree quite well with those derived by other methods, thus providing us with a new tool for BH mass determination. We find no evidence for any correlation between M and the BH spin parameter, $\alpha$, the viewing angle, $\theta$, or the height of the X-ray source above the disc, $h$. Also on average, the X-ray source lies only around 3.7 gravitational radii above the accretion disc and the viewing angles are distributed uniformly between 20 and 60 degrees. Finally, there is a tentative indication that the distribution of spin parameters may be bimodal above and below 0.62., Comment: Accepted for publication in MNRAS. The paper is 22 pages long and contains 19 figures and 2 tables

The production of artificial light curves with known statistical and variability properties is of great importance in astrophysics. Consolidating the confidence levels during cross-correlation studies, understanding the artefacts induced by sampling irregularities, establishing detection limits for future observatories are just some of the applications of simulated data sets. Currently, the widely used methodology of amplitude and phase randomisation is able to produce artificial light curves which have a given underlying power spectral density (PSD) but which are strictly Gaussian distributed. This restriction is a significant limitation, since the majority of the light curves e.g. active galactic nuclei, X-ray binaries, gamma-ray bursts show strong deviations from Gaussianity exhibiting `burst-like' events in their light curves yielding long-tailed probability distribution functions (PDFs). In this study we propose a simple method which is able to precisely reproduce light curves which match both the PSD and the PDF of either an observed light curve or a theoretical model. The PDF can be representative of either the parent distribution or the actual distribution of the observed data, depending on the study to be conducted for a given source. The final artificial light curves contain all of the statistical and variability properties of the observed source or theoretical model i.e. same PDF and PSD, respectively. Within the framework of Reproducible Research, the code, together with the illustrative example used in this manuscript, are both made publicly available in the form of an interactive Mathematica notebook., Accepted for publication in MNRAS. The paper is 23 pages long and contains 21 figures and 2 tables. The Mathematica notebook can be found in the web as part of this paper (Online Material) or at http://www.astro.soton.ac.uk/~de1e08/ArtificialLightCurves/

We present a detailed X-ray variability study of the low mass Active Galactic Nuclei (AGN) NGC 7314 using the two newly obtained XMM-Newton observations ($140$ and $130$ ks), together with two archival data sets of shorter duration ($45$ and $84$ ks). The relationship between the X-ray variability characteristics and other physical source properties (such as the black hole mass) are still relatively poorly defined, especially for low-mass AGN. We perform a new, fully analytical, power spectral density (PSD) model analysis method, which will be described in detail in a forthcoming paper, that takes into consideration the spectral distortions, caused by red-noise leak. We find that the PSD in the $0.5-10$ keV energy range, can be represented by a bending power-law with a bend around $6.7\times10^{-5}$ Hz, having a slope of $0.51$ and $1.99$ below and above the bend, respectively. Adding our bend time-scale estimate, to an already published ensemble of estimates from several AGN, supports the idea that the bend time-scale depends linearly only on the black hole mass and not on the bolometric luminosity. Moreover, we find that as the energy range increases, the PSD normalization increases and there is a hint that simultaneously the high frequency slope becomes steeper. Finally, the X-ray time-lag spectrum of NGC 7314 shows some very weak signatures of relativistic reflection, and the energy resolved time-lag spectrum, for frequencies around $3\times10^{-4}$ Hz, shows no signatures of X-ray reverberation. We show that the previous claim about ks time-delays in this source, is simply an artefact induced by the minuscule number of points entering during the time-lag estimation in the low frequency part of the time-lag spectrum (i.e. below $10^{-4}$ Hz)., Accepted for publication in MNRAS. The paper is 21 pages long and contains 15 figures and 3 tables

We present ground-based optical photometric monitoring data for NGC 5548, part of an extended multi-wavelength reverberation mapping campaign. The light curves have nearly daily cadence from 2014 January to July in nine filters (\emph{BVRI} and \emph{ugriz}). Combined with ultraviolet data from the \emph{Hubble Space Telescope} and \emph{Swift}, we confirm significant time delays between the continuum bands as a function of wavelength, extending the wavelength coverage from 1158\,\AA\ to the $z$ band ($\sim\!9160$\,\AA). We find that the lags at wavelengths longer than the {\it V} band are equal to or greater than the lags of high-ionization-state emission lines (such as He\,{\sc ii}\,$\lambda 1640$ and $\lambda 4686$), suggesting that the continuum-emitting source is of a physical size comparable to the inner broad-line region (BLR). The trend of lag with wavelength is broadly consistent with the prediction for continuum reprocessing by an accretion disk with $\tau \propto \lambda^{4/3}$. However, the lags also imply a disk radius that is 3 times larger than the prediction from standard thin-disk theory, assuming that the bolometric luminosity is 10\% of the Eddington luminosity ($L = 0.1L_{\rm Edd}$). Using optical spectra from the Large Binocular Telescope, we estimate the bias of the interband continuum lags due to BLR emission observed in the filters. We find that the bias for filters with high levels of BLR contamination ($\sim\! 20\%$) can be important for the shortest continuum lags, and likely has a significant impact on the {\it u} and {\it U} bands owing to Balmer continuum emission., Comment: 26 pages, 10 figures, accepted to ApJ. For a brief video describing the main results of this paper, please see: https://www.youtube.com/watch?v=yaYtcDvIoP0&feature=youtu.be

We present the results of a detailed study of the X-ray power spectra density (PSD) functions of twelve X-ray bright AGN, using almost all the archival XMM-Newton data. The total net exposure of the EPIC-pn light curves is larger than 350 ks in all cases (and exceeds 1 Ms in the case of 1H 0707-497). In a physical scenario in which X-ray reflection occurs in the inner part of the accretion disc of AGN, the X-ray reflection component should be a filtered echo of the X-ray continuum signal and should be equal to the convolution of the primary emission with the response function of the disc. Our primary objective is to search for these reflection features in the 5-7 keV (iron line) and 0.5-1 keV (soft) bands, where the X-ray reflection fraction is expected to be dominant. We fit to the observed periodograms two models: a simple bending power law model (BPL) and a BPL model convolved with the transfer function of the accretion disc assuming the lamp-post geometry and X-ray reflection from a homogeneous disc. We do not find any significant features in the best-fitting BPL model residuals either in individual PSDs in the iron band, soft and full band (0.3-10 keV) or in the average PSD residuals of the brightest and more variable sources (with similar black hole mass estimates). The typical amplitude of the soft and full-band residuals is around 3-5 per cent. It is possible that the expected general relativistic effects are not detected because they are intrinsically lower than the uncertainty of the current PSDs, even in the strong relativistic case in which X-ray reflection occurs on a disc around a fast rotating black hole having an X-ray source very close above it. However, we could place strong constrains to the X-ray reflection geometry with the current data sets if we knew in advance the intrinsic shape of the X-ray PSDs, particularly its high frequency slope., Comment: Accepted for publication in MNRAS. The paper is 15 pages long and contains 7 figures and 6 tables

We present X-ray spectral variability of 24 local active galactic nuclei (AGN) from the Palomar sample of nearby galaxies, as observed mainly by Swift. From hardness ratio measurements, we find that 18 AGN with low accretion rates show hardening with increasing count rate, converse to the softer-when-brighter behaviour normally observed in AGN with higher accretion rates. Two AGN show softening with increasing count rate, two show more complex behaviour, and two do not show any simple relationship. Sufficient data were available for the spectra of 13 AGN to be summed in flux-bins. In 9 of these sources, correlated luminosity-dependent changes in the photon index ($\Gamma$) of a power-law component are found to be the main cause of hardness variability. For 6 objects, with a low accretion rate as a fraction of the Eddington rate (\.m$\mathrm{_{Edd}}$), $\Gamma$ is anticorrelated with \.m$\mathrm{_{Edd}}$, i.e. `harder-when-brighter' behaviour is observed. The 3 higher-\.m$\mathrm{_{Edd}}$-rate objects show a positive correlation between $\Gamma$ and \.m$\mathrm{_{Edd}}$. This transition from harder-when-brighter at low \.m$\mathrm{_{Edd}}$s to softer-when-brighter at high \.m$\mathrm{_{Edd}}$s can be explained by a change in the dominant source of seed-photons for X-ray emission from cyclo-synchrotron emission from the Comptonising corona itself to thermal seed-photons from the accretion disc. This transition is also seen in the `hard state' of black hole X-ray binaries (BHXRBs). The results support the idea that LINERs are analogues of BHXRBs in the hard state and that Seyferts are analogues of BHXRBs in either the high-accretion-rate end of the hard state or in the hard-intermediate state., Comment: 25 pages, 10 figures, accepted for publication in the Monthly Notices of the Royal Astronomical Society

We conduct a multiwavelength continuum variability study of the Seyfert 1 galaxy NGC 5548 to investigate the temperature structure of its accretion disk. The 19 overlapping continuum light curves (1158 to 9157 angstroms) combine simultaneous HST , Swift , and ground-based observations over a 180 day period from 2014 January to July. Light-curve variability is interpreted as the reverberation response of the accretion disk to irradiation by a central time-varying point source. Our model yields the disk inclination, i, temperature T1 at 1 light day from the black hole, and a temperature-radius slope, alpha. We also infer the driving light curve and find that it correlates poorly with both the hard and soft X-ray light curves, suggesting that the X-rays alone may not drive the ultraviolet and optical variability over the observing period. We also decompose the light curves into bright, faint, and mean accretion-disk spectra. These spectra lie below that expected for a standard blackbody accretion disk accreting at L/LEdd = 0.1, Comment: V2: Oops wrong title! V1: Accepted for publication in ApJ, 20 Pages, 11 Figures

We present the spectral results from a 130 ks observation, obtained from the X-ray Multi-Mirror Mission-Newton (XMM-Newton) observatory, of the type I Seyfert galaxy Fairall 9. An X-ray hardness-ratio analysis of the light-curves, reveals a `softer-when-brighter' behaviour which is typical for radio-quiet type I Seyfert galaxies. Moreover, we analyse the high spectral-resolution data of the reflection grating spectrometer and we did not find any significant evidence supporting the presence of warm-absorber in the low X-ray energy part of the source's spectrum. This means that the central nucleus of Fairall 9 is `clean' and thus its X-ray spectral properties probe directly the physical conditions of the central engine. The overall X-ray spectrum in the 0.5-10 keV energy-range, derived from the EPIC data, can be modelled by a relativistically blurred disc-reflection model. This spectral model yields for Fairall 9 an intermediate black-hole best-fit spin parameter of $\alpha=0.39^{+0.48}_{-0.30}$.

We present the results of extensive multi-waveband monitoring of the blazar 3C~279 between 1996 and 2007 at X-ray energies (2-10 keV), optical R band, and 14.5 GHz, as well as imaging with the Very Long Baseline Array (VLBA) at 43 GHz. In all bands the power spectral density corresponds to "red noise" that can be fit by a single power law over the sampled time scales. Variations in flux at all three wavebands are significantly correlated. The time delay between high and low frequency bands changes substantially on time scales of years. A major multi-frequency flare in 2001 coincided with a swing of the jet toward a more southerly direction, and in general the X-ray flux is modulated by changes in the position angle of the jet near the core. The flux density in the core at 43 GHz--increases in which indicate the appearance of new superluminal knots--is significantly correlated with the X-ray flux. We decompose the X-ray and optical light curves into individual flares, finding that X-ray leads optical variations (XO) in 6 flares, the reverse occurs in 3 flares (OX), and there is essentially zero lag in 4 flares. Upon comparing theoretical expectations with the data, we conclude that (1) XO flares can be explained by gradual acceleration of radiating electrons to the highest energies; (2) OX flares can result from either light-travel delays of the seed photons (synchrotron self-Compton scattering) or gradients in maximum electron energy behind shock fronts; and (3) events with similar X-ray and optical radiative energy output originate well upstream of the 43 GHz core, while those in which the optical radiative output dominates occur at or downstream of the core., Accepted for publication in ApJ. 50 pages (including 23 figures and 8 tables)

Recently, it has been shown that soft-state black hole X-ray binaries and active galactic nuclei populate a plane in the space defined by the black hole mass, accretion rate and characteristic frequency. We show that this plane can be extended to hard-state objects if one allows a constant offset for the frequencies in the soft and the hard state. During a state transition the frequencies rapidly move from one scaling to the other depending on an additional parameter, possibly the disk-fraction. The relationship between frequency, mass and accretion rate can be further extended by including weakly accreting neutron stars. We explore if the lower kHz QPOs of neutron stars and the dwarf nova oscillations of white dwarfs can be included as well and discuss the physical implications of the found correlation.

We present the results of a deep 610-MHz survey of the 1^HXMM–Newton/Chandra survey area with the Giant Metre-wave Radio Telescope. The resulting maps have a resolution of ~7 arcsec and an rms noise limit of 60 μJy. To a 5σ detection limit of 300 μJy, we detect 223 sources within a survey area of 64 arcmin in diameter. We compute the 610-MHz source counts and compare them to those measured at other radio wavelengths. The well-known flattening of the Euclidean-normalized 1.4-GHz source counts below ~2 mJy, usually explained by a population of starburst galaxies undergoing luminosity evolution, is seen at 610 MHz. The 610-MHz source counts can be modelled by the same populations that explain the 1.4-GHz source counts, assuming a spectral index of −0.7 for the starburst galaxies and the steep spectrum active galactic nucleus (AGN) population. We find a similar dependence of luminosity evolution on redshift for the starburst galaxies at 610 MHz as is found at 1.4 GHz (i.e. 'Q'= 2.45^(+0.3)_(−0.4)).

We present long term monitoring of MCG-6-30-15 in X-rays, optical and near-IR wavelengths, collected over five years of monitoring. We determine the power spectrum density of all the observed bands and show that after taking into account the host contamination similar power is observed in the optical and near-IR bands. There is evidence for a correlation between the light curves of the X-ray photon flux and the optical B-band, but it is not possible to determine a lag with certainty, with the most likely value being around zero days. Strong correlation is seen between the optical and near-IR bands. Cross correlation analysis shows some complex probability distributions and lags that range from 10 to 20 days, with the near-IR following the optical variations. Filtering the light curves in frequency space shows that the strongest correlations are those corresponding to the shortest time-scales. We discuss the nature of the X-ray variability and conclude that this is intrinsic and cannot be accounted for by absorption episodes due to material intervening in the line of sight. It is also found that the lags agree with the relation tau ~ lambda^(4/3), as expected for an optically thick geometrically thin accretion disc, although for a larger disc than that predicted by the estimated black hole mass and accretion rate in MCG-6-30-15. The cross correlation analysis suggests that the torus is located at ~20 light-days from the central source and at most at ~50 light-days from the central region. This implies an AGN bolometric luminosity of ~3x10^(43) ergs/s/cm-2., 13 pages (double column), 11 figures, 1 table, accepted for publication (MNRAS)

Recent intensive Swift monitoring of the Seyfert 1 galaxy NGC 5548 yielded 282 usable epochs over 125 days across six UV/optical bands and the X-rays. This is the densest extended AGN UV/optical continuum sampling ever obtained, with a mean sampling rate, Comment: Accepted for publication in the Astrophysical Journal. Seventeen pages, 10 figures, 6 tables. See also STORM Paper I: "Space Telescope and Optical Reverberation Mapping Project. I. Ultraviolet Observations of the Seyfert 1 Galaxy NGC 5548 with the Cosmic Origins Spectrograph on Hubble Space Telescope" by G. De Rosa et al., http://arxiv.org/abs/1501.05954

We describe the first results from a six-month long reverberation-mapping experiment in the ultraviolet based on 170 observations of the Seyfert 1 galaxy NGC 5548 with the Cosmic Origins Spectrograph on the Hubble Space Telescope. Significant correlated variability is found in the continuum and broad emission lines, with amplitudes ranging from ~30% to a factor of two in the emission lines and a factor of three in the continuum. The variations of all the strong emission lines lag behind those of the continuum, with He II 1640 lagging behind the continuum by ~2.5 days and Lyman alpha 1215, C IV 1550, and Si IV 1400 lagging by ~5-6 days. The relationship between the continuum and emission lines is complex. In particular, during the second half of the campaign, all emission-line lags increased by a factor of 1.3-2 and differences appear in the detailed structure of the continuum and emission-line light curves. Velocity-resolved cross-correlation analysis shows coherent structure in lag versus line-of-sight velocity for the emission lines; the high-velocity wings of C IV respond to continuum variations more rapidly than the line core, probably indicating higher velocity BLR clouds at smaller distances from the central engine. The velocity-dependent response of Lyman alpha, however, is more complex and will require further analysis., Accepted for publication in The Astrophysical Journal. 20 pages, 6 figures, 4 tables. See also STORM Paper II "Space telescope and optical reverberation mapping project. II. Reverberation mapping of the accretion disk with SWIFT and HST" by R. Edelson et al

We present an analysis of the submillimetre/X-ray properties of 19 X-ray absorbed, Compton-thin quasi-stellar objects (QSOs) selected to have luminosities and redshifts that represent the peak of cosmic QSO activity, i.e. ∼ L* objects at 1 < z < 3. Of these, we present new data for 11 objects not previously observed at submillimetre wavelengths and additional data for a further three. The detection rate is 42 per cent, much higher than typically reported for samples of QSOs. Detection statistics show (at the 3-4σ level) that this sample of absorbed QSOs has a higher submillimetre output than a matched sample of unabsorbed QSOs. We argue that the farinfrared luminosity is produced by massive star formation. In this case, the correlation found between far-infrared luminosity and redshift can be interpreted as cosmological evolution of the star formation rate in the QSO host galaxies. Because the submillimetre luminous phase is confined to z > 1.5, the high star formation rates are consistent with a scenario in which the QSOs evolve to become local luminous elliptical galaxies. Combining these results with previously published data for X-ray unabsorbed QSOs and submillimetre-selected galaxies, we propose the following evolutionary sequence: (i) the forming galaxy is initially far-infrared luminous but X-ray weak similar to the sources discovered by the Submillimetre Common-User Bolometer Array (SCUBA); (ii) as the black hole and spheroid grow with time, a point is reached when the central QSO becomes powerful enough to terminate the star formation and eject the bulk of the fuel supply (the Compton-thin absorbed QSO phase); (iii) this transition is followed by a period of unobscured QSO activity, which subsequently declines to leave a quiescent spheroidal galaxy., FJC acknowledges support from the Spanish Ministerio de Cienca y Technología, under project ESP2003-00812. IS acknowledges support from the Royal Society.

We present a detailed study of the long-timescale X-ray variability of the Narrow Line Seyfert 1 (NLS1) Galaxy MCG-6-30-15, based on 8 years of frequent observations with the Rossi X-ray Timing Explorer. Combined with short timescale XMM-Newton observations, we derive the powerspectral density (PSD) covering 6 decades from ~10^{-8} to ~10^{-2} Hz. As with NGC4051, another NLS1, the PSD of MCG-6-30-15 is a close analogue of the PSD of a Galactic Black Hole X-ray binary system (GBH) in a `high' rather than a `low' state. Like high state GBHs, its PSD is better fitted by a smoothly bending rather than sharply breaking powerlaw model and the break frequency is 7.6^{+10}_{-3} x 10^{-5} Hz. Assuming linear scaling of break frequency with mass, we estimate the black hole (BH) mass in MCG-6-30-15 to be ~2.9^{+1.8}_{-1.6} X 10^{6} solar masses. We also derive the BH mass using a variety of new optical observations and find a value between 3.6 and 6 x 10^{6} solar masses, consistent with the mass derived from the PSD. Combining these data with revised reverberation masses for other AGN we update the BH mass/break timescale diagram. The observations are consistent with NLS1s having shorter break timescales, for a given mass, than broad line AGN, probably reflecting a higher accretion rate. However the data are also consistent with all of the X-ray bright AGN being high state objects. This result may simply be a selection effect, based on their selection from X-ray all sky catalogues, and their consequent typically high X-ray/radio ratios, which indicate high state systems.

In order to determine the relationship between the faint X-ray and faint radio source populations, and hence to help understand the X-ray and radio emission mechanisms in those faint source populations, we have made a deep 1.4GHz Very Large Array radio survey of the 13h +38deg XMM/ROSAT X-ray Survey Area. From a combined data set of 10hours B configuration data and 14hours A configuration data, maps with 3.35'' resolution and a noise limit of 7.5uJy were constructed. A complete sample of 449 sources was detected within a 30' diameter region above a 4sigma detection limit of 30uJy, at the map centre, making this one of the deepest radio surveys at this frequency. The differential source count shows a significant upturn at sub-milliJansky flux densities, similar to that seen in other deep surveys at 1.4GHz (eg the Phoenix survey), but larger than that seen in the HDF which may have been selected to be under-dense. This upturn is well modelled by the emergence of a population of medium redshift starforming galaxies which dominate at faint flux densities. The brighter source counts are well modelled by AGNs.

We present XMM-Newton Reflection Grating Spectrometer observations of the soft X-ray spectrum of NGC 4051 and explore their implications for the inner accretion disk and ionized outflow in the active galactic nucleus. We fit the soft X-ray excess with a relativistically broadened O VIII recombination spectrum, including the entire line series and recombination continuum. This plus an underlying power-law continuum provides a much better fit to the soft X-ray spectrum than a single temperature or disk blackbody plus power law. The emission-line profiles, computed for a Kerr metric around a maximally rotating black hole, reveal a sharply peaked disk emissivity law and inner radius smaller than 1.7RG. The spectrum also includes narrow absorption and emission lines from C, N, O, Ne, and Fe in an ionized outflow. Outflow column densities are relatively low and do not create significant edges in the spectrum. The small amount of absorption bolsters confidence in the detection of relativistic emission-line features. The narrow-line emitter has a large (76%) global covering fraction, leading to strong forbidden lines and filling in of the resonance absorption lines. We also find broad C VI Lyα and very broad O VII emission from the broad-line region. The narrow- and broad-line regions span large ranges in ionization parameter and may arise in a disk outflow. The ionized absorber has a large ionization range, which is inconsistent with pressure equilibrium in a multiphase medium. The mass outflow rate exceeds the accretion rate by a factor of 1000.

We present over 6.5 years of frequent X-ray observations of the NLS1 galaxy NGC4051 by RXTE together with a >100ks observation by XMM-Newton and so derive a 2-10 keV PSD covering an unprecedent frequency range of over 6.5 decades from 10^-2 Hz. This PSD is a very good match to the PSD of the galactic black hole binary system (GBH) Cyg X-1 when in a `high', rather than `low', state providing the first definite confirmation of an AGN in a `high' state. A bending, rather than sharply broken, powerlaw is a better description of the PSDs of Cyg X-1 in the high state and of NGC4051. At low frequencies the PSD of NGC4051 has a slope of -1.1 bending, at nu_B = 8 x 10^-4 Hz, to a slope of alpha_H ~ -2 (steeper at lower energies). If nu_B scales with mass, we imply a black hole mass of 3 x 10^5 Msolar in NGC4051, consistent with the 5 x 10^5 Msolar reverberation value. Hence NGC4051 is emitting at ~30% L_Edd. NGC4051 follows the same rms-flux relationship as GBHs. Phase lags and coherence between bands are also the same as in GBHs. The lag of soft by hard photons increases, and coherence decreases, as band separation increases. The lag, and coherence, are greater for variations of longer Fourier period. This behaviour suggests that higher photon energies and shorter variability timescales are associated with smaller radii. Using also data from the literature, we cannot fit all AGN to the same linear scaling of break timescale with black hole mass. However broad line AGN are consistent with a linear scaling from Cyg X-1 in its low state and NLS1 galaxies scale better with Cyg X-1 in its high state. We suggest that the relationship between black hole mass and break timescale is a function of at least one other underlying parameter, which may be accretion rate or black hole spin., Accepted for publication in MNRAS; 22 pages, 19 figures

Since 2005, the blazar 3C 454.3 has shown remarkable flaring activity at all frequencies, and during the last four years it has exhibited more than one gamma-ray flare per year, becoming the most active gamma-ray blazar in the sky. We present for the first time the multi-wavelength AGILE, SWIFT, INTEGRAL, and GASP-WEBT data collected in order to explain the extraordinary gamma-ray flare of 3C 454.3 which occurred in November 2010. On 2010 November 20 (MJD 55520), 3C 454.3 reached a peak flux (E>100 MeV) of F_gamma(p) = (6.8+-1.0)E-5 ph/cm2/s on a time scale of about 12 hours, more than a factor of 6 higher than the flux of the brightest steady gamma-ray source, the Vela pulsar, and more than a factor of 3 brighter than its previous super-flare on 2009 December 2-3. The multi-wavelength data make a thorough study of the present event possible: the comparison with the previous outbursts indicates a close similarity to the one that occurred in 2009. By comparing the broadband emission before, during, and after the gamma-ray flare, we find that the radio, optical and X-ray emission varies within a factor 2-3, whereas the gamma-ray flux by a factor of 10. This remarkable behavior is modeled by an external Compton component driven by a substantial local enhancement of soft seed photons., Accepted for publication in ApJ Letters. 18 Pages, 4 Figures, 1 Table

We report on the extreme gamma-ray activity from the FSRQ PKS 1510-089 observed by AGILE in March 2009. In the same period a radio-to-optical monitoring of the source was provided by the GASP-WEBT and REM. Moreover, several Swift ToO observations were triggered, adding important information on the source behaviour from optical/UV to hard X-rays. We paid particular attention to the calibration of the Swift/UVOT data to make it suitable to the blazars spectra. Simultaneous observations from radio to gamma rays allowed us to study in detail the correlation among the emission variability at different frequencies and to investigate the mechanisms at work. In the period 9-30 March 2009, AGILE detected an average gamma-ray flux of (311+/-21)x10^-8 ph cm^-2 s^-1 for E>100 MeV, and a peak level of (702+/-131)x10^-8 ph cm^-2 s^-1 on daily integration. The gamma-ray activity occurred during a period of increasing activity from near-IR to UV, with a flaring episode detected on 26-27 March 2009, suggesting that a single mechanism is responsible for the flux enhancement observed from near-IR to UV. By contrast, Swift/XRT observations seem to show no clear correlation of the X-ray fluxes with the optical and gamma-ray ones. However, the X-ray observations show a harder photon index (1.3-1.6) with respect to most FSRQs and a hint of harder-when-brighter behaviour, indicating the possible presence of a second emission component at soft X-ray energies. Moreover, the broad band spectrum from radio-to-UV confirmed the evidence of thermal features in the optical/UV spectrum of PKS 1510-089 also during high gamma-ray state. On the other hand, during 25-26 March 2009 a flat spectrum in the optical/UV energy band was observed, suggesting an important contribution of the synchrotron emission in this part of the spectrum during the brightest gamma-ray flare, therefore a significant shift of the synchrotron peak., 13 pages, 7 figures, 3 tables. Accepted for publication in Astronomy and Astrophysics

Competing models for broad spectral features in the soft X-ray spectrum of the Seyfert I galaxy Mrk766 are tested against data from a 130 ks XMM-Newton observation. A model including relativistically broadened Ly-alpha emission lines of OVIII, NVII and CVI is a better fit to 0.3-2 keV XMM RGS data than a dusty warm absorber. Moreover, the measured depth of neutral iron absorption lines in the spectrum is inconsistent with the magnitude of the iron edge required to produce the continuum break at 17-18Angstroms in the dusty warm absorber model. The relativistic emission line model can reproduce the broad-band (0.1-12 keV) XMM-EPIC data with the addition of a fourth line to represent emission from ionized iron at 6.7 keV and an excess due to reflection at energies above the iron line. The profile of the 6.7 keV iron line is consistent with that measured for the low energy lines. There is evidence in the RGS data at the 3sigma level for spectral features that vary with source flux. The covering fraction of warm absorber gas is estimated to be ~12%. Iron in the warm absorber is found to be overabundant with respect to CNO compared to solar values., Accepted for publication in ApJ Letters

We discuss a 1.5 day long observation of the bright variable Seyfert 1 galaxy NGC 4051 that was obtained with the X-ray and ultraviolet telescopes on XMM-Newton and extended in the X-ray band using the Rossi X-Ray Timing Explorer. The ultraviolet (2900 A) and X-ray (2-10 keV) light curves varied with peak-to-peak amplitudes of 5% and a factor of ~4, respectively. A simple interband cross-correlation function shows three broad features: a positive correlation with coefficient rM = +0.63 at a lag (Topt - TX-ray) τ ~ +0.2 days, an anticorrelation with rM = -0.65 at τ ~ +0.7 days, and a broad but slightly weaker feature with rM ≥ 0.53 over τ = 1.2-1.4 days. Because of the asymmetric sampling of the optical and X-ray data, only positive lags could be probed. Monte Carlo simulations based on synthetic uncorrelated light curves derived from the power spectrum of the real data give an ~85% confidence in features at this level. If the X-ray and UV light curves are indeed correlated, the lag is highly significant (i.e., not consistent with zero), which is suggestive of a scenario in which X-rays are being reprocessed into UV light at a distance from the X-ray source. A simple model in which the reprocessor is an inclined optically thick ring at a radius of 0.14 lt-days from the X-ray source reproduces the UV data well. A comparison of the modulated fractions in the ultraviolet and X-ray bands suggests that 5% of the 2900 A flux is produced in this way. The reprocessing site is close enough to the X-ray emitter for X-ray heating from a source with the luminosity of NGC 4051 to produce a brightening in the UV band, and the amount of UV flux emitted suggests that the covering fraction of the reprocessor is about 6%. We briefly discuss the possible physical nature of the reprocessing ring structure.

We present results from RXTE data obtained during a systematic monitoring program of four Seyfert galaxies (NGC5548, NGC5506, MCG-6-30-15 and NGC4051). We studied the variability of three hardness ratios derived from the light curves in four energy bands. All the objects show similar spectral variations in all ratios. In order to interpret the results we computed the hardness ratios corresponding to a simple spectral model of a power law plus iron line plus reflection component, considering two possibilities: a)variations with constant line equivalent width and reflection parameter R and b) variations with constant line and reflection flux. The overall, mean observed trends can be explained by spectral slope variations and a constant flux Fe line and reflection component, although the existence of a line component which is variable on short time scales cannot be excluded. Finally, we find that the data are not consistent with an increase of R with flux for individual sources, indicating that, as a single source varies, softer spectra do not correspond to larger R values.

We have developed a Monte Carlo technique to test models for the true power spectra of intermittently sampled lightcurves against the noisy, observed power spectra, and produce a reliable estimate of the goodness of fit of the given model. We apply this technique to constrain the broadband power spectra of a sample of four Seyfert galaxies monitored by the Rossi X-ray Timing Explorer (RXTE) over three years. We show that the power spectra of three of the AGN in our sample (MCG-6-30-15, NGC5506 and NGC3516) flatten significantly towards low frequencies, while the power spectrum of NGC5548 shows no evidence of flattening. We fit two models for the flattening, a `knee' model, analogous to the low-frequency break seen in the power spectra of BHXRBs in the low state (where the power-spectral slope flattens to \alpha=0) and a `high-frequency break' model (where the power-spectral slope flattens to \alpha=1), analogous to the high-frequency break seen in the high and low-state power spectra of the classic BHXRB Cyg X-1. Both models provide good fits to the power spectra of all four AGN. For both models, the characteristic frequency for flattening is significantly higher in MCG-6-30-15 than in NGC 3516 (by factor ~10) although both sources have similar X-ray luminosities, suggesting that MCG-6-30-15 has a lower black hole mass and is accreting at a higher rate than NGC 3516. Assuming linear scaling of characteristic frequencies with black hole mass, the high accretion rate implied for MCG-6-30-15 favours the high-frequency break model for this source and further suggests that MCG-6-30-15 and possibly NGC 5506, may be analogues of Cyg X-1 in the high state [ABRIDGED]., Comment: 23 pages, accepted for publication in MNRAS