Active galactic nuclei (AGNs) that show strong rest-frame optical/UV variability in their blue continuum and broad line emission are classified as changing-look AGN, or at higher luminosities, ...changing-look quasars (CLQs). These surprisingly large and sometimes rapid transitions challenge accepted models of quasar physics and duty cycles, offer several new avenues for study of quasar host galaxies, and open a wider interpretation of the cause of differences between broad and narrow-line AGN. To better characterize extreme quasar variability, we present follow-up spectroscopy as part of a comprehensive search for CLQs across the full Sloan Digital Sky Survey (SDSS) footprint using spectroscopically confirmed quasars from the SDSS DR7 catalog. Our primary selection requires large-amplitude ( mag, mag) variability over any of the available time baselines probed by the SDSS and Pan-STARRS 1 surveys. We employ photometry from the Catalina Sky Survey to verify variability behavior in CLQ candidates where available, and confirm CLQs using optical spectroscopy from the William Herschel, MMT, Magellan, and Palomar telescopes. For our adopted signal-to-noise ratio threshold on variability of broad Hβ emission, we find 17 new CLQs, yielding a confirmation rate of 20%. These candidates are at lower Eddington ratio relative to the overall quasar population, which supports a disk-wind model for the broad line region. Based on our sample, the CLQ fraction increases from 10% to roughly half as the continuum flux ratio between repeat spectra at 3420 increases from 1.5 to 6. We release a catalog of more than 200 highly variable candidates to facilitate future CLQ searches.
The origin of the X-ray emission from neutron star coalescence GW170817/GRB 170817A is a key diagnostic of the unsettled post-merger narrative, and different scenarios predict distinct evolution in ...its X-ray light curve. Due to its sky proximity to the Sun, sensitive X-ray monitoring of GW170817/GRB 170817A has not been possible since a previous detection at 16 days post-burst. We present new, deep Chandra observations of GW170817/GRB 170817A at 109 days post-burst, immediately after Sun constraints were lifted. The X-ray emission has brightened from a 0.3-8.0 keV flux of erg s−1 cm−2 at 16 days to erg s−1 cm−2 at 109 days, at a rate similar to the radio observations. This confirms that the X-ray and radio emission have a common origin. We show that the X-ray light curve is consistent with models of outflow afterglows, in which the outflow can be a cocoon shocked by the jet, dynamical ejecta from the merger, or an off-axis structured jet. Further deep X-ray monitoring can place powerful constraints on the physical parameters of these models, by both timing the passing of a synchrotron cooling break through the X-ray band and detecting the associated steepening of the X-ray photon index. Finally, the X-ray brightening strengthens the argument that simple off-axis top-hat jet models are not consistent with the latest observations of GW170817/GRB 170817A.
We report Chandra observations of GW170817, the first neutron star-neutron star merger discovered by the joint LIGO-Virgo Collaboration, and the first direct detection of gravitational radiation ...associated with an electromagnetic counterpart, Fermi short γ-ray burst GRB 170817A. The event occurred on 2017 August 17 and subsequent observations identified an optical counterpart, SSS17a, coincident with NGC 4993 (∼10″ separation). Early Chandra ( days) and Swift ( days) observations yielded non-detections at the optical position, but ∼9 days post-trigger Chandra monitoring revealed an X-ray point source coincident with SSS17a. We present two deep Chandra observations totaling ∼95 ks, collected on 2017 September 01-02 ( days). We detect X-ray emission from SSS17a with erg s−1, and a power law spectrum of . We find that the X-ray light curve from a binary NS coalescence associated with this source is consistent with the afterglow from an off-axis short γ-ray burst, with a jet angled 23° from the line of sight. This event marks both the first electromagnetic counterpart to a LIGO-Virgo gravitational-wave source and the first identification of an off-axis short GRB. We also confirm extended X-ray emission from NGC 4993 ( erg s−1) consistent with its E/S0 galaxy classification, and report two new Chandra point sources in this field, CXOU J130948 and CXOU J130946.
Despite their factor of ∼108 difference in black hole mass, several lines of evidence suggest possible similarities between black hole accretion flows in active galactic nuclei (AGN) and Galactic ...X-ray binaries. However, it is still unclear whether the geometry of the disk-corona system in X-ray binaries directly scales up to AGN and whether this analogy still holds in different accretion states. We test this AGN/X-ray binary analogy by comparing the observed correlations between the UV-to-X-ray spectral index ( OX) and Eddington ratio in AGN to those predicted from observations of X-ray binary outbursts. This approach probes the geometry of their disk-corona systems as they transition between different accretion states. We use new Chandra X-ray and ground-based rest-UV observations of faded "changing-look" quasars to extend this comparison to lower Eddington ratios of <10−2, where observations of X-ray binaries predict a softening of OX in AGN. We find that the observed correlations between the OX and Eddington ratio of AGN displays a remarkable similarity to accretion state transitions in prototypical X-ray binary outbursts, including an inversion of this correlation at a critical Eddington ratio of ∼10−2. Our results suggest that the structures of black hole accretion flows directly scale across a factor of ∼108 in black hole mass and across different accretion states, enabling us to apply theoretical models of X-ray binaries to explain AGN phenomenology.
The multi-wavelength electromagnetic afterglow from the binary neutron star merger GW170817/GRB 170817A has displayed long-term power-law brightening, and has presented challenges to post-merger ...models of the non-thermal emission. The most recent radio observations up to 200 days post-merger suggest that the afterglow has finally peaked and may now be fading, but fading has not been confirmed in the X-rays. We present new, deep Chandra observations of GW170817/GRB 170817A at 260 days post-merger that reveal an X-ray flux of erg s−1 cm−2, and confirm that the X-ray light curve is now also fading. Through rigorous comparisons with previous Chandra observations of GW170817/GRB 170817A, X-ray fading is detected between 160 and 260 days post-merger at a 4.4 significance on the basis of the X-ray data alone. We further constrain the X-ray photon index to steepen by <0.5 at 3.1 significance during this period, which disfavors the passing of the synchrotron cooling frequency through the X-ray band as the cause of the observed fading. These observations remain consistent with optically thin synchrotron afterglow emission. If this afterglow emission arises from a quasi-spherical mildly relativistic outflow, the X-ray fading suggests that the outflow is now decelerating. Alternatively, if this afterglow arises from a successful off-axis structured jet, the X-ray fading suggests that emission from the jet core has already entered the line of sight.
Abstract
Kilonovae are likely a key site of heavy
r
-process element production in the Universe, and their optical/infrared spectra contain insights into both the properties of the ejecta and the ...conditions of the
r
-process. However, the event GW170817/AT2017gfo is the only kilonova so far with well-observed spectra. To understand the diversity of absorption features that might be observed in future kilonovae spectra, we use the TARDIS Monte Carlo radiative transfer code to simulate a suite of optical spectra spanning a wide range of kilonova ejecta properties and
r
-process abundance patterns. To identify the most common and prominent absorption lines, we perform dimensionality reduction using an autoencoder, and we find spectra clusters in the latent space representation using a Bayesian Gaussian Mixture model. Our synthetic kilonovae spectra commonly display strong absorption by strontium
38
Sr
ii
, yttrium
38
Y
ii
, and zirconium
40
Zr
i–ii
, with strong lanthanide contributions at low electron fractions (
Y
e
≲ 0.25). When a new kilonova is observed, our machine-learning framework will provide context on the dominant absorption lines and key ejecta properties, helping to determine where this event falls within the larger “zoo” of kilonovae spectra.
ABSTRACT
The recent association of two flaring blazars (TXS 0506+056 and 3HSP J095507.9+355101) with IceCube neutrinos strongly suggests that blazars are at least one source of high-energy neutrinos. ...A peculiar characteristic of these two blazars is that they both appear to lie above the well-known blazar sequence, and they move further above the blazar sequence during the neutrino emission. Although the reason for this peculiarity is unclear, it may indicate a link between blazar SED properties and neutrino emission. Here, we test whether blazars lying above the blazar sequence are overrepresented in IceCube neutrino error regions. We spatially cross-match IceCube error regions with a sample of γ-ray blazars from the Fermi-LAT 4LAC catalogue, as a function of the blazars’ distance in synchrotron peak frequency above the blazar sequence. We do not find a statistically significant excess of counterparts, and thus conclude that the current data cannot confirm that the sources above the blazar sequence are a significant source of IceCube neutrinos. However, we show that the inclusion of more recently detected IceCube events increases the likelihood of an excess of blazars above the blazar sequence, and suggest that this excess may reach statistical significance in the future as the sample of IceCube neutrino detections increases.
Abstract
Low-frequency gravitational-wave experiments such as the Laser Interferometer Space Antenna and pulsar timing arrays are expected to detect individual massive black hole (MBH) binaries and ...mergers. However, secure methods of identifying the exact host galaxy of each MBH merger among the large number of galaxies in the gravitational-wave localization region are currently lacking. We investigate the distinct morphological signatures of MBH merger host galaxies, using the Romulus25 cosmological simulation. We produce mock telescope images of 201 simulated galaxies in Romulus25 hosting recent MBH mergers through stellar population synthesis and dust radiative transfer. Based on comparisons to mass- and redshift-matched control samples, we show that combining multiple morphological statistics via a linear discriminant analysis enables identification of the host galaxies of MBH mergers, with accuracies that increase with chirp mass and mass ratio. For mergers with high chirp masses (≳10
8.2
M
⊙
) and high mass ratios (≳0.5), the accuracy of this approach reaches ≳80%, and does not decline for at least ∼1 Gyr after numerical merger. We argue that these trends arise because the most distinctive morphological characteristics of MBH merger and binary host galaxies are prominent classical bulges, rather than relatively short-lived morphological disturbances from their preceding galaxy mergers. Since these bulges are formed though major mergers of massive galaxies, they lead to (and become permanent signposts for) MBH binaries and mergers that have high chirp masses and mass ratios. Our results suggest that galaxy morphology can aid in identifying the host galaxies of future MBH binaries and mergers.
Abstract
Freshly synthesized
r
-process elements in kilonovae ejecta imprint absorption features on optical spectra, as observed in the GW170817 binary neutron star merger. These spectral features ...encode insights into the physical conditions of the
r
-process and the origins of the ejected material, but associating features with particular elements and inferring the resultant abundance pattern is computationally challenging. We introduce Spectroscopic
r
-Process Abundance Retrieval for Kilonovae (
SPARK
), a modular framework to perform Bayesian inference on kilonova spectra with the goals of inferring elemental abundance patterns and identifying absorption features at early times.
SPARK
inputs an atomic line list and abundance patterns from reaction network calculations into the
TARDIS
radiative transfer code. It then performs fast Bayesian inference on observed kilonova spectra by training a Gaussian process surrogate for the approximate posteriors of kilonova ejecta parameters, via active learning. We use the spectrum of GW170817 at 1.4 days to perform the first inference on a kilonova spectrum, and recover a complete abundance pattern. Our inference shows that this ejecta was generated by an
r
-process with either (1) high electron fraction
Y
e
∼ 0.35 and high entropy
s
/
k
B
∼ 25, or, (2) a more moderate
Y
e
∼ 0.30 and
s
/
k
B
∼ 14. These parameters are consistent with a shocked, polar dynamical component, and a viscously driven outflow from a remnant accretion disk, respectively. We also recover previous identifications of strontium absorption at ∼8000 Å, and tentatively identify yttrium and/or zirconium at ≲4500 Å. Our approach will enable computationally tractable inference on the spectra of future kilonovae discovered through multimessenger observations.
We present a wide-field optical imaging search for electromagnetic counterparts to the likely neutron star-black hole (NS-BH) merger GW190814/S190814bv. This compact binary merger was detected ...through gravitational waves by the LIGO/Virgo interferometers, with masses suggestive of an NS-BH merger. We imaged the LIGO/Virgo localization region using the MegaCam instrument on the Canada-France-Hawaii Telescope (CFHT). We describe our hybrid observing strategy of both tiling and galaxy-targeted observations, as well as our image differencing and transient detection pipeline. Our observing campaign produced some of the deepest multiband images of the region between 1.7 and 8.7 days post-merger, reaching a 5 depth of g > 22.8 (AB mag) at 1.7 days and i > 23.1 and i > 23.9 at 3.7 and 8.7 days, respectively. These observations cover a mean total integrated probability of 67.0% of the localization region. We find no compelling candidate transient counterparts to this merger in our images, which suggests that the lighter object was tidally disrupted inside of the BH's innermost stable circular orbit, the transient lies outside of the observed sky footprint, or the lighter object is a low-mass BH. We use 5 source detection upper limits from our images in the NS-BH interpretation of this merger to constrain the mass of the kilonova ejecta to be Mej 0. 015M for a "blue" ( ) kilonova and Mej 0. 04M for a "red" ( ) kilonova. Our observations emphasize the key role of large-aperture telescopes and wide-field imagers such as CFHT MegaCam in enabling deep searches for electromagnetic counterparts to gravitational-wave events.