-
The SDSS-V Black Hole Mapper Reverberation Mapping Project: Multi-Line Dynamical Modeling of a Highly Variable Active Galactic Nucleus with Decade-long Light Curves
Authors:
Zachary Stone,
Yue Shen,
Scott F. Anderson,
Franz Bauer,
W. N. Brandt,
Priyanka Chakraborty,
Megan C. Davis,
Logan B. Fries,
Catherine J. Grier,
P. B. Hall,
Anton M. Koekemoer,
Mary Loli Martínez-Aldama,
Knox Long,
Sean Morrison,
Claudio Ricci,
Donald P. Schneider,
Matthew J. Temple,
Jonathan R. Trump
Abstract:
We present dynamical modeling of the broad-line region (BLR) for the highly variable AGN SDSS J141041.25+531849.0 ($z = 0.359$) using photometric and spectroscopic monitoring data from the Sloan Digital Sky Survey Reverberation Mapping project and the SDSS-V Black Hole Mapper program, spanning from early 2013 to early 2023. We model the geometry and kinematics of the BLR in the H$β$, H$α$, and MgI…
▽ More
We present dynamical modeling of the broad-line region (BLR) for the highly variable AGN SDSS J141041.25+531849.0 ($z = 0.359$) using photometric and spectroscopic monitoring data from the Sloan Digital Sky Survey Reverberation Mapping project and the SDSS-V Black Hole Mapper program, spanning from early 2013 to early 2023. We model the geometry and kinematics of the BLR in the H$β$, H$α$, and MgII, emission lines for three different time periods to measure the potential change of structure within the BLR across time and line species. We consistently find a moderately edge-on $(i_{\rm full-state} = 53.29^{\circ} \,{}^{+7.29}_{-6.55})$ thick-disk $(θ_{\rm opn, \; full-state} = 54.86^{\circ} \,{}^{+5.83}_{-4.74})$ geometry for all BLRs, with a joint estimate for the mass of the supermassive black hole (SMBH) for each of three time periods, yielding $\log_{10}(M_{\rm BH} / M_{\odot}) = 7.66^{+0.12}_{-0.13}$ when using the full dataset. The inferred individual virial factor $f$ $\sim 1$ is significantly smaller than the average factor for a local sample of dynamically modeled AGNs. There is strong evidence for non-virial motion, with over $80\%$ of clouds on inflowing/outflowing orbits. We analyze the change in model parameters across emission lines, finding the radii of BLRs for the emission lines are consistent with the following relative sizes $R_{\rm Hβ} \lesssim R_{\rm MgII } \lesssim R_{\rm Hα}$. Comparing results across time, we find $R_{\rm low-state} \lesssim R_{\rm high-state}$, with the change in BLR size for H$β$, being more significant than for the other two lines. The data also reveal complex, time-evolving, and potentially transient dynamics of the BLR gas over decade-long timescales, encouraging for future dynamical modeling of fine-scale BLR kinematics.
△ Less
Submitted 8 August, 2024;
originally announced August 2024.
-
A Comprehensive Overview of the Lebesgue Differentiation Theorem in Coq
Authors:
Reynald Affeldt,
Zachary Stone
Abstract:
Formalization of real analysis offers a chance to rebuild traditional proofs of important theorems as unambiguous theories that can be interactively explored. This paper provides a comprehensive overview of the Lebesgue Differentiation Theorem formalized in the Coq proof assistant, from which the first Fundamental Theorem of Calculus (FTC) for the Lebesgue integral is obtained as a corollary. Prov…
▽ More
Formalization of real analysis offers a chance to rebuild traditional proofs of important theorems as unambiguous theories that can be interactively explored. This paper provides a comprehensive overview of the Lebesgue Differentiation Theorem formalized in the Coq proof assistant, from which the first Fundamental Theorem of Calculus (FTC) for the Lebesgue integral is obtained as a corollary. Proving the first FTC in this way has the advantage of decomposing into loosely-coupled theories of moderate size and of independent interest that lend themselves well to incremental and collaborative development. We explain how we formalize all the topological constructs and all the standard lemmas needed to eventually relate the definitions of derivability and of Lebesgue integration of MathComp-Analysis, a formalization of analysis developed on top of the Mathematical Components library. In the course of this experiment, we substantially enrich MathComp-Analysis and even devise a new proof for Urysohn's lemma.
△ Less
Submitted 30 June, 2024; v1 submitted 26 March, 2024;
originally announced March 2024.
-
High-quality Extragalactic Legacy-field Monitoring (HELM) with DECam
Authors:
Ming-Yang Zhuang,
Qian Yang,
Yue Shen,
Monika Adamow,
Douglas N. Friedel,
R. A. Gruendl,
Xin Liu,
Paul Martini,
Timothy M. C. Abbott,
Scott F. Anderson,
Roberto J. Assef,
Franz E. Bauer,
Rich Bielby,
W. N. Brandt,
Colin J. Burke,
Jorge Casares,
Yu-Ching Chen,
Gisella De Rosa,
Alex Drlica-Wagner,
Tom Dwelly,
Alice Eltvedt,
Gloria Fonseca Alvarez,
Jianyang Fu,
Cesar Fuentes,
Melissa L. Graham
, et al. (23 additional authors not shown)
Abstract:
High-quality Extragalactic Legacy-field Monitoring (HELM) is a long-term observing program that photometrically monitors several well-studied extragalactic legacy fields with the Dark Energy Camera (DECam) imager on the CTIO 4m Blanco telescope. Since Feb 2019, HELM has been monitoring regions within COSMOS, XMM-LSS, CDF-S, S-CVZ, ELAIS-S1, and SDSS Stripe 82 with few-day cadences in the…
▽ More
High-quality Extragalactic Legacy-field Monitoring (HELM) is a long-term observing program that photometrically monitors several well-studied extragalactic legacy fields with the Dark Energy Camera (DECam) imager on the CTIO 4m Blanco telescope. Since Feb 2019, HELM has been monitoring regions within COSMOS, XMM-LSS, CDF-S, S-CVZ, ELAIS-S1, and SDSS Stripe 82 with few-day cadences in the $(u)gri(z)$ bands, over a collective sky area of $\sim 38$ deg${\rm ^2}$. The main science goal of HELM is to provide high-quality optical light curves for a large sample of active galactic nuclei (AGNs), and to build decades-long time baselines when combining past and future optical light curves in these legacy fields. These optical images and light curves will facilitate the measurements of AGN reverberation mapping lags, as well as studies of AGN variability and its dependences on accretion properties. In addition, the time-resolved and coadded DECam photometry will enable a broad range of science applications from galaxy evolution to time-domain science. We describe the design and implementation of the program and present the first data release that includes source catalogs and the first $\sim 3.5$ years of light curves during 2019A--2022A.
△ Less
Submitted 8 February, 2024;
originally announced February 2024.
-
TESS photometry of the nova eruption in V606 Vul: asymmetric photosphere and multiple ejections?
Authors:
Kirill V. Sokolovsky,
Elias Aydi,
Konstantin Malanchev,
Colin J. Burke,
Koji Mukai,
J. L. Sokoloski,
Brian D. Metzger,
Kirill E. Atapin,
Aleksandre A. Belinski,
Yu-Ching Chen,
Laura Chomiuk,
Pavol A. Dubovsky,
Claude-Andre Faucher-Giguere,
Rebekah A. Hounsell,
Natalia P. Ikonnikova,
Vsevolod Yu. Lander,
Junyao Li,
Justin D. Linford,
Amy J. Mioduszewski,
Isabella Molina,
Ulisse Munari,
Sergey A. Potanin,
Robert M. Quimby,
Michael P. Rupen,
Simone Scaringi
, et al. (46 additional authors not shown)
Abstract:
Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The ligh…
▽ More
Lightcurves of many classical novae deviate from the canonical "fast rise - smooth decline" pattern and display complex variability behavior. We present the first TESS-space-photometry-based investigation of this phenomenon. We use TESS Sector 41 full-frame images to extract a lightcurve of the slow Galactic nova V606 Vul that erupted nine days prior to the start of the TESS observations. The lightcurve covers the first of two major peaks of V606 Vul that was reached 19 days after the start of the eruption. The nova reached its brightest visual magnitude V=9.9 in its second peak 64 days after the eruption onset, following the completion of Sector 41 observations. To increase the confidence level of the extracted lightcurve, we performed the analysis using four different codes implementing the aperture photometry (Lightkurve, VaST) and image subtraction (TESSreduce, tequila_shots) and find good agreement between them. We performed ground-based photometric and spectroscopic monitoring to complement the TESS data. The TESS lightcurve reveals two features: periodic variations (0.12771d, 0.01mag average peak-to-peak amplitude) that disappeared when the source was within 1mag of peak optical brightness and a series of isolated mini-flares (with peak-to-peak amplitudes of up to 0.5mag) appearing at seemingly random times. We interpret the periodic variations as the result of azimuthal asymmetry of the photosphere engulfing the nova-hosting binary that was distorted by and rotating with the binary. Whereas we use spectra to associate the two major peaks in the nova lightcurve with distinct episodes of mass ejection, the origin of mini-flares remains elusive.
△ Less
Submitted 8 November, 2023;
originally announced November 2023.
-
The Sloan Digital Sky Survey Reverberation Mapping Project: Key Results
Authors:
Yue Shen,
Catherine J. Grier,
Keith Horne,
Zachary Stone,
Jennifer I. Li,
Qian Yang,
Yasaman Homayouni,
Jonathan R. Trump,
Scott F. Anderson,
W. N. Brandt,
Patrick B. Hall,
Luis C. Ho,
Linhua Jiang,
Patrick Petitjean,
Donald P. Schneider,
Charling Tao,
Fergus. R. Donnan,
Yusra AlSayyad,
Matthew A. Bershady,
Michael R. Blanton,
Dmitry Bizyaev,
Kevin Bundy,
Yuguang Chen,
Megan C. Davis,
Kyle Dawson
, et al. (22 additional authors not shown)
Abstract:
We present the final data from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a precursor to the SDSS-V Black Hole Mapper Reverberation Mapping program. This data set includes 11-year photometric and 7-year spectroscopic light curves for 849 broad-line quasars over a redshift range of 0.1<z<4.5 and a luminosity range of Lbol=1E44-47.5 erg/s, along with spectral and variabili…
▽ More
We present the final data from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a precursor to the SDSS-V Black Hole Mapper Reverberation Mapping program. This data set includes 11-year photometric and 7-year spectroscopic light curves for 849 broad-line quasars over a redshift range of 0.1<z<4.5 and a luminosity range of Lbol=1E44-47.5 erg/s, along with spectral and variability measurements. We report 23, 81, 125, and 110 reverberation mapping lags (relative to optical continuum variability) for broad Halpha, Hbeta, MgII and CIV using the SDSS-RM sample, spanning much of the luminosity and redshift ranges of the sample. Using 30 low-redshift RM AGNs with dynamical-modeling black hole masses, we derive a new estimate of the average virial factor of <log f>=0.62+-0.07 for the line dispersion measured from the RMS spectrum. The intrinsic scatter of individual virial factors is 0.31+-0.07 dex, indicating a factor of two systematic uncertainty in RM black hole masses. Our lag measurements reveal significant R-L relations for Hbeta and MgII at high redshift, consistent with the latest measurements based on heterogeneous samples. While we are unable to robustly constrain the slope of the R-L relation for CIV given the limited dynamical range in luminosity, we found substantially larger scatter in CIV lags at fixed L1350. Using the SDSS-RM lag sample, we derive improved single-epoch (SE) mass recipes for Hbeta, MgII and CIV, which are consistent with their respective RM masses as well as between the SE recipes from two different lines, over the luminosity range probed by our sample. The new Hbeta and MgII recipes are approximately unbiased estimators at given RM masses, but there are systematic biases in the CIV recipe. The intrinsic scatter of SE masses around RM masses is ~0.45 dex for Hbeta and MgII, increasing to ~0.58 dex for CIV.
△ Less
Submitted 1 April, 2024; v1 submitted 1 May, 2023;
originally announced May 2023.
-
Temperature Fluctuations in Quasar Accretion Discs from Spectroscopic Monitoring Data
Authors:
Zachary Stone,
Yue Shen
Abstract:
Neustadt & Kochanek (2022, hereafter NK22) proposed a new method to reconstruct the temperature perturbation map (as functions of time and disc radius) of AGN accretion discs using multi-wavelength photometric light curves. We apply their technique to 100 quasars at $z=0.5-2$ from the Sloan Digital Sky Survey Reverberation Mapping project, using multi-epoch spectroscopy that covers rest-frame UV-o…
▽ More
Neustadt & Kochanek (2022, hereafter NK22) proposed a new method to reconstruct the temperature perturbation map (as functions of time and disc radius) of AGN accretion discs using multi-wavelength photometric light curves. We apply their technique to 100 quasars at $z=0.5-2$ from the Sloan Digital Sky Survey Reverberation Mapping project, using multi-epoch spectroscopy that covers rest-frame UV-optical continuum emission from the quasar and probes days to months timescales. Consistent with NK22 for low-redshift AGNs, we find that the dominant pattern of disc temperature perturbations is either slow inward/outward moving waves with typical amplitudes $δT/T_0\sim 10\%$ traveling at $\sim 0.01-0.1c$, with a typical radial frequency of $\sim$ 0.5 dex in $\log R$, or incoherent perturbations. In nearly none of the cases do we find clear evidence for coherent, fast outgoing temperature perturbations at the speed of light, reminiscent of the lamppost model; but such lamppost signals may be present in some quasars for limited periods of the monitoring data. Using simulated data, we demonstrate that high-fidelity temperature perturbation maps can be recovered with high-quality monitoring spectroscopy, with limited impact from seasonal gaps in the data. On the other hand, reasonable temperature perturbation maps can be reconstructed with high-cadence photometric light curves from the Vera C. Rubin Observatory Legacy Survey of Space and Time. Our findings, together with NK22, suggest that internal disc processes are the main driver for temperature fluctuations in AGN accretion discs over days to months timescales.
△ Less
Submitted 20 September, 2023; v1 submitted 13 October, 2022;
originally announced October 2022.
-
Optical Variability of Quasars with 20-Year Photometric Light Curves
Authors:
Zachary Stone,
Yue Shen,
Colin J. Burke,
Yu-Ching Chen,
Qian Yang,
Xin Liu,
R. A. Gruendl,
M. Adamów,
F. Andrade-Oliveira,
J. Annis,
D. Bacon,
E. Bertin,
S. Bocquet,
D. Brooks,
D. L. Burke,
A. Carnero Rosell,
M. Carrasco Kind,
J. Carretero,
L. N. da Costa,
M. E. S. Pereira,
J. De Vicente,
S. Desai,
H. T. Diehl,
P. Doel,
I. Ferrero
, et al. (30 additional authors not shown)
Abstract:
We study the optical $gri$ photometric variability of a sample of 190 quasars within the SDSS Stripe 82 region that have long-term photometric coverage during $\sim 1998-2020$ with SDSS, PanSTARRS-1, the Dark Energy Survey, and dedicated follow-up monitoring with Blanco 4m/DECam. With on average $\sim 200$ nightly epochs per quasar per filter band, we improve the parameter constraints from a Dampe…
▽ More
We study the optical $gri$ photometric variability of a sample of 190 quasars within the SDSS Stripe 82 region that have long-term photometric coverage during $\sim 1998-2020$ with SDSS, PanSTARRS-1, the Dark Energy Survey, and dedicated follow-up monitoring with Blanco 4m/DECam. With on average $\sim 200$ nightly epochs per quasar per filter band, we improve the parameter constraints from a Damped Random Walk (DRW) model fit to the light curves over previous studies with 10-15 yr baselines and $\lesssim 100$ epochs. We find that the average damping timescale $τ_{\rm DRW}$ continues to rise with increased baseline, reaching a median value of $\sim 750$ days ($g$ band) in the rest-frame of these quasars using the 20-yr light curves. Some quasars may have gradual, long-term trends in their light curves, suggesting that either the DRW fit requires very long baselines to converge, or that the underlying variability is more complex than a single DRW process for these quasars. Using a subset of quasars with better-constrained $τ_{\rm DRW}$ (less than 20\% of the baseline), we confirm a weak wavelength dependence of $τ_{\rm DRW}\propto λ^{0.51\pm0.20}$. We further quantify optical variability of these quasars over days to decades timescales using structure function (SF) and power spectrum density (PSD) analyses. The SF and PSD measurements qualitatively confirm the measured (hundreds of days) damping timescales from the DRW fits. However, the ensemble PSD is steeper than that of a DRW on timescales less than $\sim$ a month for these luminous quasars, and this second break point correlates with the longer DRW damping timescale.
△ Less
Submitted 8 June, 2022; v1 submitted 7 January, 2022;
originally announced January 2022.