Bayesian analysis of quasar lightcurves with a running optimal average : new time delay measurements of COSMOGRAIL gravitationally lensed quasars

Funding: KH and JVHS acknowledge support from UK Science and Technology Facilities Council grant ST/R00824/1.We present a new method of modelling time-series data based on the running optimal average (ROA). By identifying the effective number of parameters for the ROA model, in terms of the shape and width of its window function and the times and accuracies of the data, we enable a Bayesian analysis, optimising the ROA width, along with other model parameters, by minimising the Bayesian Information Criterion (BIC) and sampling joint posterior parameter distributions using MCMC methods. For analysis of quasar lightcurves, our implementation of ROA modelling can measure time delays among lightcurves at different wavelengths or from different images of a lensed quasar and, in future work, be used to inter-calibrate lightcurve data from different telescopes and estimate the shape and thus the power-density spectrum of the lightcurve. Our noise model implements a robust treatment of outliers and error-bar adjustments to account for additional variance or poorly-quantified uncertainties. Tests with simulated data validate the parameter uncertainty estimates. We compare ROA delay measurements with results from cross-correlation and from JAVELIN, which models lightcurves with a prior on the power-density spectrum. We analyse published COSMOGRAIL lightcurves of multi-lensed quasar lightcurves and present the resulting measurements of the inter-image time delays and detection of microlensing effects.PostprintPeer reviewe

Optical continuum reverberation in the dwarf Seyfert nucleus of NGC 4395

Funding information: LCO telescope time was granted by NOIRLab through the Mid-Scale Innovations Program (MSIP). MSIP is funded by NSF. This Letter is based on observations made with the MuSCAT3 instrument, developed by the Astrobiology Center and under financial supports by JSPS KAKENHI (JP18H05439) and JST PRESTO (JPMJPR1775), at Faulkes Telescope North on Maui, HI, operated by the Las Cumbres Observatory. Research at UC Irvine was supported by NSF grant AST-1907290. D.H.G.-B. acknowledges CONACYT support #319800. J.V.H.S. acknowledges support from STFC grant ST/R000824/1.In order to constrain the size of the optical continuum emission region in the dwarf Seyfert 1 galaxy NGC 4395 through reverberation mapping, we carried out high-cadence photometric monitoring in the griz filter bands on two consecutive nights in 2022 April using the four-channel MuSCAT3 camera on the Faulkes Telescope North at Haleakalā Observatory. Correlated variability across the griz bands is clearly detected, and the r-, i-, and z-band light curves show lags of - + 7.72 +1.01 -1.09, 14.16 +1.22-1.25, and 20.78  +1.99-2.09 minutes with respect to the g band when measured using the full-duration light curves. When lags are measured for each night separately, the Night 2 data exhibit lower cross-correlation amplitudes and shorter lags than the Night 1 light curves. Using the full-duration lags, we find that the lag–wavelength relationship is consistent with the τ ∝ λ4/3 dependence found for more luminous active galactic nuclei. Combining our results with continuum lags measured for other objects, the lag between g and z band scales with optical continuum luminosity as τgz ∝ L0.56±0.05, similar to the scaling of broad-line region size with luminosity, reinforcing recent evidence that diffuse continuum emission from the broad-line region may contribute substantially to optical continuum variability and reverberation lags.Publisher PDFPeer reviewe

Intensive swift and LCO monitoring of PG 1302–102 : active galactic nucleus disk reverberation mapping of a supermassive black hole binary candidate

This work is supported by the NANOGrav National Science Foundation Physics Frontiers Center award No. 2020265 and NASA grant 80NSSC24K0251. Research at UC Irvine was supported by NSF grant AST-1907290. E.M.C. gratefully acknowledges support from the NSF through grant No. AST-1909199.We present an intensive multiwavelength monitoring campaign of the quasar PG 1302−102 with Swift and the Las Cumbres Observatory network telescopes. At z ∼ 0.3, it tests the limits of the reverberation mapping (RM) technique in probing the accretion disk around a supermassive black hole (SMBH) and extends the parameter space to high masses and high accretion rates. This is also the first time the RM technique has been applied to test disk structures predicted in the SMBH binary model that has been suggested for this source. PG 1302−102 was observed at a ∼daily cadence for ∼9 months in 14 bands spanning from X-ray to UV and optical wavelengths, and it shows moderate to significant levels of variability correlated between wavelengths. We measure the interband time lags, which are consistent with a τ ∝ λ 4/3 relation as expected from standard disk reprocessing, albeit with large uncertainties. The disk size implied by the lag spectrum is consistent with the expected disk size for its black hole mass within uncertainties. While the source resembles other reverberation-mapped active galactic nuclei in many respects, and we do not find evidence supporting the prevalent hypothesis that it hosts an SMBH binary, we demonstrate the feasibility of studying SMBH binaries from this novel angle and suggest possibilities for the LSST Deep Drilling Fields.Peer reviewe

A complex dust morphology in the high-luminosity AGN Mrk 876

Recent models for the inner structure of active galactic nuclei (AGN) advocate the presence of a radiatively accelerated, dusty outflow launched from the outer regions of the accretion disk. Here we present the first near-infrared (near-IR) variable (rms) spectrum for the high-luminosity, nearby AGN Mrk 876. We find that it tracks the accretion disk spectrum out to longer wavelengths than the mean spectrum due to a reduced dust emission. The implied outer accretion disk radius is consistent with the infrared results predicted by a contemporaneous optical accretion disk reverberation mapping campaign and much larger than the self-gravity radius. The reduced flux variability of the hot dust could be either due to the presence of a secondary, constant dust component in the mean spectrum or introduced by the destructive superposition of the dust and accretion disk variability signals or some combination of both. Assuming thermal equilibrium for optically thin dust, we derive the luminosity-based dust radius for different grain properties using our measurement of the temperature. We find that in all cases considered the values are significantly larger than the dust response time measured by IR photometric monitoring campaigns, with the least discrepancy present relative to the result for a wavelength-independent dust emissivity law, i.e. a blackbody, which is appropriate for large grain sizes. This result can be well explained by assuming a flared, disk-like structure for the hot dust.Comment: 18 pages, 7 figures; accepted to Ap

Stellar Astrophysics and Exoplanet Science with the Maunakea Spectroscopic Explorer (MSE)

The Maunakea Spectroscopic Explorer (MSE) is a planned 11.25-m aperture facility with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. A rebirth of the 3.6m Canada-France-Hawaii Telescope on Maunakea, MSE will use 4332 fibers operating at three different resolving powers (R ~ 2500, 6000, 40000) across a wavelength range of 0.36-1.8mum, with dynamical fiber positioning that allows fibers to match the exposure times of individual objects. MSE will enable spectroscopic surveys with unprecedented scale and sensitivity by collecting millions of spectra per year down to limiting magnitudes of g ~ 20-24 mag, with a nominal velocity precision of ~100 m/s in high-resolution mode. This white paper describes science cases for stellar astrophysics and exoplanet science using MSE, including the discovery and atmospheric characterization of exoplanets and substellar objects, stellar physics with star clusters, asteroseismology of solar-like oscillators and opacity-driven pulsators, studies of stellar rotation, activity, and multiplicity, as well as the chemical characterization of AGB and extremely metal-poor stars.Comment: 31 pages, 11 figures; To appear as a chapter for the Detailed Science Case of the Maunakea Spectroscopic Explore

The Sloan Digital Sky Survey Reverberation Mapping project : photometric g and i light curves

Funding: P.H. acknowledges the support of the Natural Sciences and Engineering Research Council of Canada (NSERC), funding reference number 2017-05983. C.J.G., W.N.B., and D.P.S. acknowledge support from NSF grant AST-1517113. Y.S. acknowledges support from an Alfred P. Sloan Research Fellowship and NSF grant AST-1715579. L.C.H. acknowledges the National Key R&D Program of China (2016YFA0400702) and the National Science Foundation of China (11721303, 11991052). J.V.H.S. and K.H. acknowledge support from a STFC grant ST/R000824/1. C.S.K. is supported by NSF grants AST-1814440 and AST-1908570.The Sloan Digital Sky Survey (SDSS) Reverberation Mapping program monitors 849 active galactic nuclei (AGNs) both spectroscopically and photometrically. The photometric observations used in this work span over 4 yr and provide an excellent baseline for variability studies of these objects. We present the photometric light curves from 2014 to 2017 obtained by the Steward Observatory's Bok telescope and the Canada-France-Hawaii telescope with MegaCam. We provide details on the data acquisition and processing of the data from each telescope, the difference imaging photometry used to produce the light curves, and the calculation of a variability index to quantify each AGN's variability. We find that the Welch-Stetson J index provides a useful characterization of AGN variability and can be used to select AGNs for further study.PostprintPeer reviewe

An irradiated brown-dwarf companion to an accreting white dwarf

Interacting compact binary systems provide a natural laboratory in which to study irradiated substellar objects. As the mass-losing secondary (donor) in these systems makes a transition from the stellar to the substellar regime, it is also irradiated by the primary (compact accretor)1, 2. The internal and external energy fluxes are both expected to be comparable in these objects, providing access to an unexplored irradiation regime. The atmospheric properties of donors are largely unknown3, but could be modified by the irradiation. To constrain models of donor atmospheres, it is necessary to obtain accurate observational estimates of their physical properties (masses, radii, temperatures and albedos). Here we report the spectroscopic detection and characterization of an irradiated substellar donor in an accreting white-dwarf binary system. Our near-infrared observations allow us to determine a model-independent mass estimate for the donor of 0.055 ± 0.008 solar masses and an average spectral type of L1 ± 1, supporting both theoretical predictions and model-dependent observational constraints that suggest that the donor is a brown dwarf. Our time-resolved data also allow us to estimate the average irradiation-induced temperature difference between the dayside and nightside of the substellar donor (57 kelvin) and the maximum difference between the hottest and coolest parts of its surface (200 kelvin). The observations are well described by a simple geometric reprocessing model with a bolometric (Bond) albedo of less than 0.54 at the 2σ confidence level, consistent with high reprocessing efficiency, but poor lateral heat redistribution in the atmosphere of the brown-dwarf donor4, 5. These results add to our knowledge of binary evolution, in that the donor has survived the transition from the stellar to the substellar regime, and of substellar atmospheres, in that we have been able to test a regime in which the irradiation and the internal energy of a brown dwarf are comparable

AGN STORM 2. I. First results: A Change in the Weather of Mrk 817

We present the first results from the ongoing, intensive, multiwavelength monitoring program of the luminous Seyfert 1 galaxy Mrk 817. While this active galactic nucleus was, in part, selected for its historically unobscured nature, we discovered that the X-ray spectrum is highly absorbed, and there are new blueshifted, broad, and narrow UV absorption lines, which suggest that a dust-free, ionized obscurer located at the inner broad-line region partially covers the central source. Despite the obscuration, we measure UV and optical continuum reverberation lags consistent with a centrally illuminated Shakura–Sunyaev thin accretion disk, and measure reverberation lags associated with the optical broad-line region, as expected. However, in the first 55 days of the campaign, when the obscuration was becoming most extreme, we observe a de-coupling of the UV continuum and the UV broad emission-line variability. The correlation recovered in the next 42 days of the campaign, as Mrk 817 entered a less obscured state. The short C IV and Lyα lags suggest that the accretion disk extends beyond the UV broad-line region. Unified