Search

We present a velocity-resolved reverberation mapping analysis of the hypervariable quasar RM160 (SDSS J141041.25+531849.0) at z = 0.359 with 153 spectroscopic epochs of data representing a ten-year baseline (2013-2023). We split the baseline into two regimes based on the 3x flux increase in the light curve: a 'low state' phase during the years 2013-2019 and a 'high state' phase during the years 2022-2023. The velocity-resolved lag profiles (VRLP) indicate that gas with different kinematics dominates the line emission in different states. The H\b{eta} VRLP begins with a signature of inflow onto the BLR in the 'low state', while in the 'high state' it is flatter with less signature of inflow. The H{\alpha} VRLP begins consistent with a virialized BLR in the 'low state', while in the 'high state' shows a signature of inflow. The differences in the kinematics between the Balmer lines and between the 'low state' and the 'high state' suggests complex BLR dynamics. We find that the BLR radius and velocity (both FWHM and {\sigma}) do not obey a constant virial product throughout the monitoring period. We find that BLR lags and continuum luminosity are correlated, consistent with rapid response of the BLR gas to the illuminating continuum. The BLR kinematic profile changes in unpredictable ways that are not related to continuum changes and reverberation lag. Our observations indicate that non-virial kinematics can significantly contribute to observed line profiles, suggesting caution for black-hole mass estimation in luminous and highly varying quasars like RM160., Comment: 23 pages, 17 figures

Dual and off-nucleus active supermassive black holes are expected to be common in the hierarchical structure formation paradigm, but their identification at parsec scales remains a challenge due to strict angular resolution requirements. We conduct a systematic study using the Very Long Baseline Array (VLBA) to examine 23 radio-bright candidate dual and off-nucleus quasars. The targets are selected by a novel astrometric technique ("varstrometry") from Gaia, aiming to identify dual or off-nucleus quasars at (sub)kilo-parsec scales. Among these quasars, 8 exhibit either multiple radio components or significant (>3$\sigma$) positional offsets between the VLBA and Gaia positions. The radio emission from the three candidates which exhibit multiple radio components is likely to originate from small-scale jets based on their morphology. Among the remaining five candidates with significant VLBA-Gaia offsets, three are identified as potential dual quasars at parsec scales, one is likely attributed to small-scale jets, and the origin of the last candidate remains unclear. We explore alternative explanations for the observed VLBA-Gaia offsets. We find no evidence for optical jets at kilo-parsec scales, nor any contamination to Gaia astrometric noise from the host galaxy; misaligned coordinate systems are unlikely to account for our offsets. Our study highlights the promise of the varstrometry technique in discovering candidate dual or off-nucleus quasars and emphasizes the need for further confirmation and investigation to validate and understand these intriguing candidates., Comment: 17 pages, 8 figures, 4 tables, accepted to ApJ

Unexpectedly strong X-ray emission from extragalactic radio jets on kiloparsec scales has been one of the major discoveries of Chandra, the only X-ray observatory capable of sub-arcsecond-scale imaging. The origin of this X-ray emission, which appears as a second spectral component from that of the radio emission, has been debated for over two decades. The most commonly assumed mechanism is inverse Compton upscattering of the Cosmic Microwave Background (IC-CMB) by very low-energy electrons in a still highly relativistic jet. Under this mechanism, no variability in the X-ray emission is expected. Here we report the detection of X-ray variability in the large-scale jet population, using a novel statistical analysis of 53 jets with multiple Chandra observations. Taken as a population, we find that the distribution of p-values from a Poisson model is strongly inconsistent with steady emission, with a global p-value of 1.96e-4 under a Kolmogorov-Smirnov test against the expected Uniform (0,1) distribution. These results strongly imply that the dominant mechanism of X-ray production in kpc-scale jets is synchrotron emission by a second population of electrons reaching multi-TeV energies. X-ray variability on the time-scale of months to a few years implies extremely small emitting volumes much smaller than the cross-section of the jet., Comment: Published in Nature Astronomy 29 May 2023; Supplemental Information and Excel File included in 1st submission; Published article in most recent

While supermassive black holes are ubiquitous features of galactic nuclei, only a small minority are observed during episodes of luminous accretion. The physical mechanism(s) driving the onset of fueling and ignition in these active galactic nuclei (AGN) are still largely unknown for many galaxies and AGN-selection criteria. Attention has focused on AGN triggering by means of major galaxy mergers gravitationally funneling gas towards the galactic center, with evidence both for and against this scenario. However, several recent studies have found that radio-loud AGN overwhelmingly reside in ongoing or recent major galaxy mergers. In this study, we test the hypothesis that major galaxy mergers are important triggers for radio-loud AGN activity in powerful quasars during cosmic noon (1 < z < 2). To this end, we compare Hubble Space Telescope WFC3/IR observations of the z > 1 3CR radio-loud broad-lined quasars to three matched radio-quiet quasar control samples. We find strong evidence for major-merger activity in nearly all radio-loud AGN, in contrast to the much lower merger fraction in the radio-quiet AGN. These results suggest major galaxy mergers are key ingredients to launching powerful radio jets. Given many of our radio-loud quasars are blue, our results present a possible challenge to the "blow-out" paradigm of galaxy evolution models in which blue quasars are the quiescent end result following a period of red quasar feedback initiated by a galaxy merger. Finally, we find a tight correlation between black hole mass and host galaxy luminosity for these different high-redshift AGN samples inconsistent with those observed for local elliptical galaxies., Comment: Published by ApJ

Over 150 resolved, kpc-scale X-ray jets hosted by active galactic nuclei have been discovered with the Chandra X-ray Observatory. A significant fraction of these jets have an X-ray spectrum either too high in flux or too hard to be consistent with the high-energy extension of the radio-to-optical synchrotron spectrum, a subtype we identify as Multiple Spectral Component (MSC) X-ray jets. A leading hypothesis for the origin of the X-rays is the inverse-Compton scattering of the cosmic microwave background by the same electron population producing the radio-to-optical synchrotron spectrum (known as the IC/CMB model). In this work, we test the IC/CMB model in 45 extragalactic X-ray jets using observations from the Fermi Large Area Telescope to look for the expected high level of gamma-ray emission, utilizing observations from the Atacama Large Millimeter/submillimeter Array (ALMA) and the Hubble Space Telescope (HST) when possible to best constrain the predicted gamma-ray flux. Including this and previous works, we now find the IC/CMB model to be ruled out in a total of 24/45 MSC X-ray jets due to its overprediction for the observed MeV-to-GeV gamma-ray flux. We present additional evidence against the IC/CMB model, including the relative X-ray-to-radio relativistic beaming in these sources, and the general mismatch between radio and X-ray spectral indexes. Finally, we present upper limits on the large scale bulk-flow Lorentz factors for all jets based on the Fermi upper limits, which suggest that these jets are at most mildly relativistic., Comment: Accepted for publication by MNRAS (October 21, 2022)

Previous Ku-band (15 GHz) imaging with data obtained from the Very Long Baseline Array (VLBA) had shown two compact, sub-pc components at the location of a presumed kpc-scale radio core in the Seyfert galaxy NGC 7674. It was then presumed that these two unresolved and compact components were dual radio cores corresponding to two supermassive black holes (SMBHs) accreting surrounding gas and launching radio-bright relativistic jets. However, utilizing the original VLBA dataset used to claim the detection of a binary SMBH, in addition to later multi-epoch/multi-frequency datatsets obtained from both the VLBA and the European VLBI Network, we find no evidence to support the presence of a binary SMBH. We place stringent upper limits to the flux densities of any sub-pc-scale radio cores which are at least an order of magnitude lower than the original VLBI radio-core detections, directly challenging the original binary SMBH detection claim. With this in mind, we discuss the possible reasons for the non-detection of any VLBI radio cores in our imaging, the possibility of a binary SMBH still residing in NGC 7674, and the prospect of future observations shedding further light on the true nature of this active galactic nucleus., Comment: Accepted by ApJ

The broad high-energy spectral component in blazars is usually attributed to various inverse Compton scattering processes in the relativistic jet, but has not been clearly identified in most cases due to degeneracies in physical models. AP Librae, a low-synchrotron-peaking BL Lac object (LBL) detected in 2015 by H.E.S.S. at very high energies (VHE; $>$ 0.5 TeV), has an extremely broad high-energy spectrum, covering $\sim$ 9 decades in energy. Standard synchrotron self-Compton models generally fail to reproduce the VHE emission, which has led to the suggestion that it might arise not from the blazar core, but on kiloparsec scales from inverse Compton scattering of cosmic microwave background (CMB) photons by a still-relativistic jet (IC/CMB). IC/CMB models for the TeV emission of AP Librae in prior works have implied a high level of infrared emission from the kpc-scale jet. With newly obtained Hubble Space Telescope imaging, we obtain a deep upper limit on the kpc-scale jet emission at 1.6 $\mu$m, well below the expected level. High-resolution ALMA imaging in bands 3-9 reveals a residual dust disk signature after core subtraction, with a clearly thermal spectrum, and an extent ($\sim$500 pc) which matches with a non-jet residual emission seen after PSF subtraction in our 1.6 $\mu$m HST imaging. We find that the unusually broad GeV and VHE emission in AP Librae can be reproduced through the combined IC scattering of photons from the CMB and the dust disk, respectively, by electrons in both the blazar core and sub-kpc jet., Comment: 25 pages, 14 figures. Minor revisions. Accepted for publication in the Astrophysical Journal

Galaxy mergers occur frequently in the early universe and bring multiple supermassive black holes (SMBHs) into the nucleus, where they may eventually coalesce. Identifying post-merger-scale (i.e., <~a few kpc) dual SMBHs is a critical pathway to understanding their dynamical evolution and successive mergers. While serendipitously discovering kpc-scale dual SMBHs at z<1 is possible, such systems are elusive at z>2, but critical to constraining the progenitors of SMBH mergers. The redshift z~2 also marks the epoch of peak activity of luminous quasars, hence probing this spatial regime at high redshift is of particular significance in understanding the evolution of quasars. However, given stringent resolution requirements, there is currently no confirmed <10 kpc physical SMBH pair at z>2. Here we report two sub-arcsec double quasars at z>2 discovered from a targeted search with a novel astrometric technique, demonstrating a high success rate (~50%) in this systematic approach. These high-redshift double quasars could be the long-sought kpc-scale dual SMBHs, or sub-arcsec gravitationally-lensed quasar images. One of these double quasars (at z=2.95) was spatially resolved with optical spectroscopy, and slightly favors the scenario of a physical quasar pair with a projected separation of 3.5 kpc (0.46"). Follow-up observations of double quasars discovered by this targeted approach will be able to provide the first observational constraints on kpc-scale dual SMBHs at z>2., Comment: Published in Nature Astronomy (April 2021)

In several previous studies, quasars exhibiting broad emission lines with >1000 km/s velocity offsets with respect to the host galaxy rest frame have been discovered. One leading hypothesis for the origin of these velocity-offset broad lines is the dynamics of a binary supermassive black hole (SMBH). We present high-resolution radio imaging of 34 quasars showing these velocity-offset broad lines with the Very Long Baseline Array (VLBA), aimed at finding evidence for the putative binary SMBHs (such as dual radio cores), and testing the competing physical models. We detect exactly half of the target sample from our VLBA imaging, after implementing a 5 detection limit. While we do not resolve double radio sources in any of the targets, we obtain limits on the instantaneous projected separations of a radio-emitting binary for all of the detected sources under the assumption that a binary still exists within our VLBA angular resolution limits. We also assess the likelihood that a radio-emitting companion SMBH exists outside of our angular resolution limits, but its radio luminosity is too weak to produce a detectable signal in the VLBA data. Additionally, we compare the precise sky positions afforded by these data to optical positions from both the SDSS and Gaia DR2 source catalogs. We find projected radio/optical separations on the order of 10 pc for three quasars. Finally, we explore how future multi-wavelength campaigns with optical, radio, and X-ray observatories can help discriminate further between the competing physical models., Comment: Accepted for publishing by ApJ

Over the past two decades, the most commonly adopted explanation for high and hard X-ray emission in resolved quasar jets has been inverse Compton upscattering of the Cosmic Microwave Background (IC/CMB), which requires jets which remain highly relativistic on 10-1000 kpc scales. In more recent years various lines of observational evidence, including gamma-ray upper limits, have disfavored this explanation in favor of a synchrotron origin. While the IC/CMB model generally predicts a high level of gamma-ray emission, it has never been detected. Here we report the detection of a low-state Fermi/LAT gamma-ray spectrum associated with two jetted AGN which is consistent with the predictions of the IC/CMB model for their X-ray emission. We have used archival multiwavelength observations to make precise predictions for the expected minimum flux in the GeV band, assuming that the X-ray emission from the kpc-scale jet is entirely due to the IC/CMB process. In both sources -- OJ 287 and PKS 1510-089 -- the minimum-detected gamma-ray flux level agrees with predictions. Both sources exhibit extreme superluminal proper motions relative to their jet power, which argues for the well-aligned jets required by the IC/CMB model. In the case of PKS~1510-089, it cannot be ruled out that the minimum gamma-ray flux level is due to a low state of the variable core which only matches the IC/CMB prediction by chance. Continued long-term monitoring with the Fermi/LAT could settle this issue by detecting a plateau signature in the recombined light-curve which would clearly signal the presence of a non-variable emission component., Comment: Accepted for publication in ApJ Letters

We report an optical/UV jet and counterjet in M84, previously unreported in archival HST imaging. With archival VLA, ALMA, and Chandra imaging, we examine the first well-sampled spectral energy distribution of the inner jet of M84, where we find that multiple co-spatial spectral components are required. In particular, the ALMA data reveal that the radio spectrum of all four knots in the jet turns over at approximately 100 GHz, which requires a second component for the bright optical/UV emission. Further, the optical/UV has a soft spectrum and is inconsistent with the relatively flat X-ray spectrum, which indicates a third component at higher energies. Using archival VLA imaging, we have measured the proper motion of the innermost knots at 0.9+/-0.6 and 1.1+/-0.4 c, which when combined with the low jet-to-counterjet flux ratio yields an orientation angle for the system of 74 (+9,-18) degrees. In the radio, we find high fractional polarization of the inner jet of up to 30% while in the optical no polarization is detected (< 8%). We investigate different scenarios for explaining the particular multi-component SED of the knots. Inverse Compton models are ruled out due to the extreme departure from equipartition and the unrealistically high total jet power required. The multi-component SED can be naturally explained within a leptohadronic scenario, but at the cost of very high power in relativistic protons. A two-component synchrotron model remains a viable explanation, but more theoretical work is needed to explain the origin and properties of the electron populations., Comment: 20 pages, 12 figures; accepted for publication in ApJ

Recent multi-wavelength work led by the Boston University blazar group (e.g., Marscher et al. 2010) strongly suggests that a fraction of the blazar flares seen by the Fermi Large Area Telescope (LAT) take place a few to several pc away from the central engine. However, at such distances from the central engine, there is no adequate external photon field to provide the seed photons required for producing the observed GeV emission under leptonic inverse Compton (IC) models. A possible solution is a spine-sheath geometry for the emitting region (MacDonald et al. 2015, but see Nalewajko et al. 2014). Here we use the current view of the molecular torus (e.g., Elitzur, 2012; Netzer 2015) in which the torus extends a few pc beyond the dust sublimation radius with dust clouds distributed with a declining density for decreasing polar angle. We show that for a spine-sheath blazar jet embedded in the torus, the wide beaming pattern of the synchrotron radiation of the relatively slow sheath will heat molecular clouds whose subsequent IR radiation will be seen highly boosted in the spine comoving frame, and that under reasonable conditions this photon field can dominate over the sheath photons directly entering the spine. If the sheath is sufficiently luminous it will sublimate the dust, and if the sheath synchrotron radiation extends to optical-UV energies (as may happen during flares), this will illuminate the sublimated dust clouds to produce emission lines that will vary in unison with the optical-UV continuum, as has been very recently reported for blazar CTA 102 (Jorstad et al. 2017)., Comment: Accepted for Publishing in the Astrophysical Journal

Since its launch in 1999, the Chandra X-ray observatory has discovered several dozen X-ray jets associated with powerful quasars. In many cases the X-ray spectrum is hard and appears to come from a second spectral component. The most popular explanation for the kpc-scale X-ray emission in these cases has been inverse-Compton (IC) scattering of Cosmic Microwave Background (CMB) photons by relativistic electrons in the jet (the IC/CMB model). Requiring the IC/CMB emission to reproduce the observed X-ray flux density inevitably predicts a high level of gamma-ray emission which should be detectable with the Fermi Large Area Telescope (LAT). In previous work, we found that gamma-ray upper limits from the large scale jets of 3C 273 and PKS 0637-752 violate the predictions of the IC/CMB model. Here we present Fermi/LAT flux density upper limits for the X-ray jets of four additional sources: PKS 1136-135, PKS 1229-021, PKS 1354+195, and PKS 2209+080, and show that these limits violate the IC/CMB predictions at a very high significance level. We also present new Hubble Space Telescope (HST) observations of the quasar PKS 2209+080 showing a newly detected optical jet, and Atacama Large Millimeter/submillimeter Array (ALMA) band 3 and 6 observations of all four sources, which provide key constraints on the spectral shape that enable us to rule out the IC/CMB model.

The Chandra X-ray observatory has discovered several dozen anomalously X-ray-bright jets associated with powerful quasars. A popular explanation for the X-ray flux from the knots in these jets is that relativistic synchrotron-emitting electrons inverse-Compton scatter Cosmic Microwave Background (CMB) photons to X-ray energies (the IC/CMB model). This model predicts a high gamma-ray flux which should be detectable by the Fermi Large Area Telescope (LAT) for many sources. GeV-band upper limits from Fermi/LAT for the well-known anomalous X-ray jet in PKS 0637-752 were previously shown in Meyer et al., (2015) to violate the predictions of the IC/CMB model. Previously, measurements of the jet synchrotron spectrum, important for accurately predicting the gamma-ray flux level, were lacking between radio and infrared wavelengths. Here we present new Atacama Large Millimeter/submillimeter Array (ALMA) observations of the large-scale jet at 100, 233, and 319 GHz which further constrain the synchrotron spectrum, supporting the previously published empirical model. We also present updated limits from the Fermi/LAT using the new `Pass 8' calibration and approximately 30% more time on source. With these deeper limits we rule out the IC/CMB model at the 8.7 sigma level. Finally, we demonstrate that complete knowledge of the synchrotron SED is critical in evaluating the IC/CMB model., Comment: Accepted for publication in ApJ Letters (v2 corrects typo in author list)