Showing 1–50 of 93 results for author: Kelly, B

A Parameter Study of the Electromagnetic Signatures of an Analytical Mini-Disk Model for Supermassive Binary Black Hole Systems

Authors: Kaitlyn Porter, Scott C. Noble, Eduardo M. Gutierrez, Joaquin Pelle, Manuela Campanelli, Jeremy Schnittman, Bernard J. Kelly

Abstract: Supermassive black holes (SMBHs) are thought to be located at the centers of most galactic nuclei. When galaxies merge they form supermassive black hole binary (SMBHB) systems and these central SMBHs will also merge at later times, producing gravitational waves (GWs). Because galaxy mergers are likely gas-rich environments, SMBHBs are also potential sources of electromagnetic (EM) radiation. The E… ▽ More Supermassive black holes (SMBHs) are thought to be located at the centers of most galactic nuclei. When galaxies merge they form supermassive black hole binary (SMBHB) systems and these central SMBHs will also merge at later times, producing gravitational waves (GWs). Because galaxy mergers are likely gas-rich environments, SMBHBs are also potential sources of electromagnetic (EM) radiation. The EM signatures depend on gas dynamics, orbital dynamics, and radiation processes. The gas dynamics are governed by general relativistic magnetohydrodynamics (MHD) in a time-dependent spacetime. Numerically solving the MHD equations for a time-dependent binary spacetime is computationally expensive. Therefore, it is challenging to conduct a full exploration of the parameter space of these systems and the resulting EM signatures. We have developed an analytical accretion disk model for the mini-disks of an SMBHB system and produced images and light curves using a general relativistic ray-tracing code and a superimposed harmonic binary black hole metric. This analytical model greatly reduces the time and computational resources needed to explore these systems, while incorporating some key information from simulations. We present a parameter space exploration of the SMBHB system in which we have studied the dependence of the EM signatures on the spins of the black holes (BHs), the mass ratio, the accretion rate, the viewing angle, and the initial binary separation. Additionally, we study how the commonly used fast-light approximation affects the EM signatures and evaluate its validity in GRMHD simulations. △ Less

Submitted 4 July, 2024; originally announced July 2024.

Waveform Modelling for the Laser Interferometer Space Antenna

Authors: LISA Consortium Waveform Working Group, Niayesh Afshordi, Sarp Akçay, Pau Amaro Seoane, Andrea Antonelli, Josu C. Aurrekoetxea, Leor Barack, Enrico Barausse, Robert Benkel, Laura Bernard, Sebastiano Bernuzzi, Emanuele Berti, Matteo Bonetti, Béatrice Bonga, Gabriele Bozzola, Richard Brito, Alessandra Buonanno, Alejandro Cárdenas-Avendaño, Marc Casals, David F. Chernoff, Alvin J. K. Chua, Katy Clough, Marta Colleoni, Mekhi Dhesi, Adrien Druart , et al. (121 additional authors not shown)

Abstract: LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmologic… ▽ More LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the beginnings of inspirals that will venture into the ground-based detectors' view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA's discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This white paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome. △ Less

Submitted 20 December, 2023; v1 submitted 2 November, 2023; originally announced November 2023.

Comments: 239 pages, 11 figures, white paper from the LISA Consortium Waveform Working Group, invited for submission to Living Reviews in Relativity, updated with comments from community

Addition of tabulated equation of state and neutrino leakage support to IllinoisGRMHD

Authors: Leonardo R. Werneck, Zachariah B. Etienne, Ariadna Murguia-Berthier, Roland Haas, Federico Cipolletta, Scott C. Noble, Lorenzo Ennoggi, Federico G. Lopez Armengol, Bruno Giacomazzo, Thiago Assumpção, Joshua Faber, Tanmayee Gupte, Bernard J. Kelly, Julian H. Krolik

Abstract: We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These new features are provided by two new, NRPy+-based codes: NRPyEOS, which performs highly efficient EOS table lookups and interpolations, and NRPyLeakage, which imp… ▽ More We have added support for realistic, microphysical, finite-temperature equations of state (EOS) and neutrino physics via a leakage scheme to IllinoisGRMHD, an open-source GRMHD code for dynamical spacetimes in the Einstein Toolkit. These new features are provided by two new, NRPy+-based codes: NRPyEOS, which performs highly efficient EOS table lookups and interpolations, and NRPyLeakage, which implements a new, AMR-capable neutrino leakage scheme in the Einstein Toolkit. We have performed a series of strenuous validation tests that demonstrate the robustness of these new codes, particularly on the Cartesian AMR grids provided by Carpet. Furthermore, we show results from fully dynamical GRMHD simulations of single unmagnetized neutron stars, and magnetized binary neutron star mergers. This new version of IllinoisGRMHD, as well as NRPyEOS and NRPyLeakage, is pedagogically documented in Jupyter notebooks and fully open source. The codes will be proposed for inclusion in an upcoming version of the Einstein Toolkit. △ Less

Submitted 14 December, 2023; v1 submitted 30 August, 2022; originally announced August 2022.

Comments: 20 pages, 9 figures. v3 erratum: fixed typo in Eq. (21)

Journal ref: Phys. Rev. D 107, 044037 (2023)

Handing off the outcome of binary neutron star mergers for accurate and long-term post-merger simulations

Authors: Federico G. Lopez Armengol, Zachariah B. Etienne, Scott C. Noble, Bernard J. Kelly, Leonardo R. Werneck, Brendan Drachler, Manuela Campanelli, Federico Cipolletta, Yosef Zlochower, Ariadna Murguia-Berthier, Lorenzo Ennoggi, Mark Avara, Riccardo Ciolfi, Joshua Faber, Grace Fiacco, Bruno Giacomazzo, Tanmayee Gupte, Trung Ha, Julian H. Krolik, Vassilios Mewes, Richard O'Shaughnessy, Jesús M. Rueda-Becerril, Jeremy Schnittman

Abstract: We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at larg… ▽ More We perform binary neutron star (BNS) merger simulations in full dynamical general relativity with IllinoisGRMHD, on a Cartesian grid with adaptive-mesh refinement. After the remnant black hole has become nearly stationary, the evolution of the surrounding accretion disk on Cartesian grids over long timescales (1s) is suboptimal, as Cartesian coordinates over-resolve the angular coordinates at large distances, and the accreting plasma flows obliquely across coordinate lines dissipating angular momentum artificially from the disk. To address this, we present the Handoff, a set of computational tools that enables the transfer of general relativistic magnetohydrodynamic (GRMHD) and spacetime data from IllinoisGRMHD to HARM3D, a GRMHD code that specializes in modeling black hole accretion disks in static spacetimes over long timescales, making use of general coordinate systems with spherical topology. We demonstrate that the Handoff allows for a smooth and reliable transition of GRMHD fields and spacetime data, enabling us to efficiently and reliably evolve BNS dynamics well beyond merger. We also discuss future plans, which involve incorporating advanced equations of state and neutrino physics into BNS simulations using the \handoff approach. △ Less

Submitted 31 October, 2022; v1 submitted 17 December, 2021; originally announced December 2021.

Comments: 27 pages, 22 figures

Journal ref: Phys. Rev. D 106, 083015 (2022)

HARM3D+NUC: A new method for simulating the post-merger phase of binary neutron star mergers with GRMHD, tabulated EOS and neutrino leakage

Authors: Ariadna Murguia-Berthier, Scott C. Noble, Luke F. Roberts, Enrico Ramirez-Ruiz, Leonardo R. Werneck, Michael Kolacki, Zachariah B. Etienne, Mark Avara, Manuela Campanelli, Riccardo Ciolfi, Federico Cipolletta, Brendan Drachler, Lorenzo Ennoggi, Joshua Faber, Grace Fiacco, Bruno Giacomazzo, Tanmayee Gupte, Trung Ha, Bernard J. Kelly, Julian H. Krolik, Federico G. Lopez Armengol, Ben Margalit, Tim Moon, Richard O'Shaughnessy, Jesús M. Rueda-Becerril , et al. (3 additional authors not shown)

Abstract: The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merge… ▽ More The first binary neutron star merger has already been detected in gravitational waves. The signal was accompanied by an electromagnetic counterpart including a kilonova component powered by the decay of radioactive nuclei, as well as a short $γ$-ray burst. In order to understand the radioactively-powered signal, it is necessary to simulate the outflows and their nucleosynthesis from the post-merger disk. Simulating the disk and predicting the composition of the outflows requires general relativistic magnetohydrodynamical (GRMHD) simulations that include a realistic, finite-temperature equation of state (EOS) and self-consistently calculating the impact of neutrinos. In this work, we detail the implementation of a finite-temperature EOS and the treatment of neutrinos in the GRMHD code HARM3D+NUC, based on HARM3D. We include formal tests of both the finite-temperature EOS and the neutrino leakage scheme. We further test the code by showing that, given conditions similar to those of published remnant disks following neutron star mergers, it reproduces both recombination of free nucleons to a neutron-rich composition and excitation of a thermal wind. △ Less

Submitted 31 January, 2022; v1 submitted 9 June, 2021; originally announced June 2021.

Comments: 19 pages, 14 figures, published in ApJ

GPU-Accelerated Hierarchical Bayesian Inference with Application to Modeling Cosmic Populations: CUDAHM

Authors: János M. Szalai-Gindl, Thomas J. Loredo, Brandon C. Kelly, István Csabai, Tamás Budavári, László Dobos

Abstract: We describe a computational framework for hierarchical Bayesian inference with simple (typically single-plate) parametric graphical models that uses graphics processing units (GPUs) to accelerate computations, enabling deployment on very large datasets. Its C++ implementation, CUDAHM (CUDA for Hierarchical Models) exploits conditional independence between instances of a plate, facilitating massive… ▽ More We describe a computational framework for hierarchical Bayesian inference with simple (typically single-plate) parametric graphical models that uses graphics processing units (GPUs) to accelerate computations, enabling deployment on very large datasets. Its C++ implementation, CUDAHM (CUDA for Hierarchical Models) exploits conditional independence between instances of a plate, facilitating massively parallel exploration of the replication parameter space using the single instruction, multiple data architecture of GPUs. It provides support for constructing Metropolis-within-Gibbs samplers that iterate between GPU-accelerated robust adaptive Metropolis sampling of plate-level parameters conditional on upper-level parameters, and Metropolis-Hastings sampling of upper-level parameters on the host processor conditional on the GPU results. CUDAHM is motivated by demographic problems in astronomy, where density estimation and linear and nonlinear regression problems must be addressed for populations of thousands to millions of objects whose features are measured with possibly complex uncertainties. We describe a thinned latent point process framework for modeling such demographic data. We demonstrate accurate GPU-accelerated parametric conditional density deconvolution for simulated populations of up to 300,000 objects in ~1 hour using a single NVIDIA Tesla K40c GPU. Supplementary material provides details about the CUDAHM API and the demonstration problem. △ Less

Submitted 17 May, 2021; originally announced May 2021.

Comments: 28 pages, 7 figures, 2 appendices

Journal ref: Abridged version published as "GPU-accelerated hierarchical Bayesian estimation of luminosity functions using flux-limited observations with photometric noise" in Astronomy & Computing (2018)

Electromagnetic Emission from a Binary Black Hole Merger Remnant in Plasma: Field Alignment and Plasma Temperature

Authors: Bernard J. Kelly, Zachariah B. Etienne, Jacob Golomb, Jeremy D. Schnittman, John G. Baker, Scott C. Noble, Geoffrey Ryan

Abstract: Comparable-mass black-hole mergers generically result in moderate to highly spinning holes, whose spacetime curvature will significantly affect nearby matter in observable ways. We investigate how the moderate spin of a post-merger Kerr black hole immersed in a plasma with initially uniform density and uniform magnetic field affects potentially observable accretion rates and energy fluxes. Varying… ▽ More Comparable-mass black-hole mergers generically result in moderate to highly spinning holes, whose spacetime curvature will significantly affect nearby matter in observable ways. We investigate how the moderate spin of a post-merger Kerr black hole immersed in a plasma with initially uniform density and uniform magnetic field affects potentially observable accretion rates and energy fluxes. Varying the initial specific internal energy of the plasma over two decades, we find very little change in steady-state mass accretion rate or Poynting luminosity, except at the lowest internal energies, where fluxes do not exhibit steady-state behavior during the simulation timescale. Fixing the internal energy and varying the initial fixed magnetic-field amplitude and orientation, we find that the steady-state Poynting luminosity depends strongly on the initial field angle with respect to the black hole spin axis, while the matter accretion rate is more stable until the field angle exceeds $\sim 45\degree$. The proto-jet formed along the black hole spin-axis conforms to a thin, elongated cylinder near the hole, while aligning with the asymptotic magnetic field at large distances. △ Less

Submitted 26 March, 2021; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 13 pages, 14 figures. Minor changes to match published version

Journal ref: Phys. Rev. D 103, 063039 (2021)

Meta-Learning for One-Class Classification with Few Examples using Order-Equivariant Network

Authors: Ademola Oladosu, Tony Xu, Philip Ekfeldt, Brian A. Kelly, Miles Cranmer, Shirley Ho, Adrian M. Price-Whelan, Gabriella Contardo

Abstract: This paper presents a meta-learning framework for few-shots One-Class Classification (OCC) at test-time, a setting where labeled examples are only available for the positive class, and no supervision is given for the negative example. We consider that we have a set of `one-class classification' objective-tasks with only a small set of positive examples available for each task, and a set of trainin… ▽ More This paper presents a meta-learning framework for few-shots One-Class Classification (OCC) at test-time, a setting where labeled examples are only available for the positive class, and no supervision is given for the negative example. We consider that we have a set of `one-class classification' objective-tasks with only a small set of positive examples available for each task, and a set of training tasks with full supervision (i.e. highly imbalanced classification). We propose an approach using order-equivariant networks to learn a 'meta' binary-classifier. The model will take as input an example to classify from a given task, as well as the corresponding supervised set of positive examples for this OCC task. Thus, the output of the model will be 'conditioned' on the available positive example of a given task, allowing to predict on new tasks and new examples without labeled negative examples. In this paper, we are motivated by an astronomy application. Our goal is to identify if stars belong to a specific stellar group (the 'one-class' for a given task), called \textit{stellar streams}, where each stellar stream is a different OCC-task. We show that our method transfers well on unseen (test) synthetic streams, and outperforms the baselines even though it is not retrained and accesses a much smaller part of the data per task to predict (only positive supervision). We see however that it doesn't transfer as well on the real stream GD-1. This could come from intrinsic differences from the synthetic and real stream, highlighting the need for consistency in the 'nature' of the task for this method. However, light fine-tuning improve performances and outperform our baselines. Our experiments show encouraging results to further explore meta-learning methods for OCC tasks. △ Less

Submitted 21 May, 2021; v1 submitted 8 July, 2020; originally announced July 2020.

Building A Field: The Future of Astronomy with Gravitational Waves, A State of The Profession Consideration for Astro2020

Authors: Kelly Holley-Bockelmann, Joey Shapiro Key, Brittany Kamai, Robert Caldwell, Warren Brown, Bill Gabella, Karan Jani, Quentin Baghi, John Baker, Jillian Bellovary, Pete Bender, Emanuele Berti, T. J. Brandt, Curt Cutler, John W. Conklin, Michael Eracleous, Elizabeth C. Ferrara, Bernard J. Kelly, Shane L. Larson, Jeff Livas, Maura McLaughlin, Sean T. McWilliams, Guido Mueller, Priyamvada Natarajan, Norman Rioux , et al. (6 additional authors not shown)

Abstract: Harnessing the sheer discovery potential of gravitational wave astronomy will require bold, deliberate, and sustained efforts to train and develop the requisite workforce. The next decade requires a strategic plan to build -- from the ground up -- a robust, open, and well-connected gravitational wave astronomy community with deep participation from traditional astronomers, physicists, data scienti… ▽ More Harnessing the sheer discovery potential of gravitational wave astronomy will require bold, deliberate, and sustained efforts to train and develop the requisite workforce. The next decade requires a strategic plan to build -- from the ground up -- a robust, open, and well-connected gravitational wave astronomy community with deep participation from traditional astronomers, physicists, data scientists, and instrumentalists. This basic infrastructure is sorely needed as an enabling foundation for research. We outline a set of recommendations for funding agencies, universities, and professional societies to help build a thriving, diverse, and inclusive new field. △ Less

Submitted 16 December, 2019; originally announced December 2019.

The Laser Interferometer Space Antenna: Unveiling the Millihertz Gravitational Wave Sky

Authors: John Baker, Jillian Bellovary, Peter L. Bender, Emanuele Berti, Robert Caldwell, Jordan Camp, John W. Conklin, Neil Cornish, Curt Cutler, Ryan DeRosa, Michael Eracleous, Elizabeth C. Ferrara, Samuel Francis, Martin Hewitson, Kelly Holley-Bockelmann, Ann Hornschemeier, Craig Hogan, Brittany Kamai, Bernard J. Kelly, Joey Shapiro Key, Shane L. Larson, Jeff Livas, Sridhar Manthripragada, Kirk McKenzie, Sean T. McWilliams , et al. (17 additional authors not shown)

Abstract: The first terrestrial gravitational wave interferometers have dramatically underscored the scientific value of observing the Universe through an entirely different window, and of folding this new channel of information with traditional astronomical data for a multimessenger view. The Laser Interferometer Space Antenna (LISA) will broaden the reach of gravitational wave astronomy by conducting the… ▽ More The first terrestrial gravitational wave interferometers have dramatically underscored the scientific value of observing the Universe through an entirely different window, and of folding this new channel of information with traditional astronomical data for a multimessenger view. The Laser Interferometer Space Antenna (LISA) will broaden the reach of gravitational wave astronomy by conducting the first survey of the millihertz gravitational wave sky, detecting tens of thousands of individual astrophysical sources ranging from white-dwarf binaries in our own galaxy to mergers of massive black holes at redshifts extending beyond the epoch of reionization. These observations will inform - and transform - our understanding of the end state of stellar evolution, massive black hole birth, and the co-evolution of galaxies and black holes through cosmic time. LISA also has the potential to detect gravitational wave emission from elusive astrophysical sources such as intermediate-mass black holes as well as exotic cosmological sources such as inflationary fields and cosmic string cusps. △ Less

Submitted 26 July, 2019; v1 submitted 15 July, 2019; originally announced July 2019.

Comments: White Paper submitted to Astro2020 (2020 Decadal Survey on Astronomy and Astrophysics). v2: fixed a reference

Space Telescope and Optical Reverberation Mapping Project. VIII. Time Variability of Emission and Absorption in NGC 5548 Based on Modeling the Ultraviolet Spectrum

Authors: G. A. Kriss, G. De Rosa, J. Ely, B. M. Peterson, J. Kaastra, M. Mehdipour, G. J. Ferland, M. Dehghanian, S. Mathur, R. Edelson, K. T. Korista, N. Arav, A. J. Barth, M. C. Bentz, W. N. Brandt, D. M. Crenshaw, E. Dalla Bontà, K. D. Denney, C. Done, M. Eracleous, M. M. Fausnaugh, E. Gardner, M. R. Goad, C. J. Grier, Keith Horne , et al. (142 additional authors not shown)

Abstract: 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-correcte… ▽ More 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$α$ 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$α$ 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. △ Less

Submitted 12 July, 2019; v1 submitted 8 July, 2019; originally announced July 2019.

Comments: 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

The Event Horizon General Relativistic Magnetohydrodynamic Code Comparison Project

Authors: Oliver Porth, Koushik Chatterjee, Ramesh Narayan, Charles F. Gammie, Yosuke Mizuno, Peter Anninos, John G. Baker, Matteo Bugli, Chi-kwan Chan, Jordy Davelaar, Luca Del Zanna, Zachariah B. Etienne, P. Chris Fragile, Bernard J. Kelly, Matthew Liska, Sera Markoff, Jonathan C. McKinney, Bhupendra Mishra, Scott C. Noble, Héctor Olivares, Ben Prather, Luciano Rezzolla, Benjamin R. Ryan, James M. Stone, Niccolò Tomei , et al. (3 additional authors not shown)

Abstract: Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic… ▽ More Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at close to the event horizon scale by the GRAVITY experiment motivate the development of numerical source models that solve the equations of general relativistic magnetohydrodynamics (GRMHD). Here we compare GRMHD solutions for the evolution of a magnetized accretion flow where turbulence is promoted by the magnetorotational instability from a set of nine GRMHD codes: Athena++, BHAC, Cosmos++, ECHO, H-AMR, iharm3D, HARM-Noble, IllinoisGRMHD and KORAL. Agreement between the codes improves as resolution increases, as measured by a consistently applied, specially developed set of code performance metrics. We conclude that the community of GRMHD codes is mature, capable, and consistent on these test problems. △ Less

Submitted 5 August, 2019; v1 submitted 9 April, 2019; originally announced April 2019.

Comments: Accepted version for publication in ApJS (May 28th, 2019)

Multi-Messenger Astrophysics with Pulsar Timing Arrays

Authors: Luke Zoltan Kelley, Maria Charisi, Sarah Burke-Spolaor, Joseph Simon, Laura Blecha, Tamara Bogdanovic, Monica Colpi, Julie Comerford, Daniel J. D'Orazio, Massimo Dotti, Michael Eracleous, Matthew Graham, Jenny E. Greene, Zoltán Haiman, Kelly Holley-Bockelmann, Erin Kara, Bernard Kelly, S. Komossa, Shane L. Larson, Xin Liu, Chung-Pei Ma, Scott Noble, Vasileios Paschalidis, Roman R. Rafikov, Vikram Ravi , et al. (7 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will… ▽ More Pulsar timing arrays (PTAs) are on the verge of detecting low-frequency gravitational waves (GWs) from supermassive black hole binaries (SMBHBs). With continued observations of a large sample of millisecond pulsars, PTAs will reach this major milestone within the next decade. Already, SMBHB candidates are being identified by electromagnetic surveys in ever-increasing numbers; upcoming surveys will enhance our ability to detect and verify candidates, and will be instrumental in identifying the host galaxies of GW sources. Multi-messenger (GW and electromagnetic) observations of SMBHBs will revolutionize our understanding of the co-evolution of SMBHs with their host galaxies, the dynamical interactions between binaries and their galactic environments, and the fundamental physics of accretion. Multi-messenger observations can also make SMBHBs 'standard sirens' for cosmological distance measurements out to $z\simeq0.5$. LIGO has already ushered in breakthrough insights in our knowledge of black holes. The multi-messenger detection of SMBHBs with PTAs will be a breakthrough in the years $2020-2030$ and beyond, and prepare us for LISA to help complete our views of black hole demographics and evolution at higher redshifts. △ Less

Submitted 18 March, 2019; originally announced March 2019.

Comments: Submitted to the Astro2020 decadal review. This is one of five core white papers written by members of the NANOGrav Collaboration, see also: J.Cordes et al., S.R.Taylor et al., X.Siemens et al., and E.Fonseca et al

Tests of General Relativity and Fundamental Physics with Space-based Gravitational Wave Detectors

Authors: Emanuele Berti, Enrico Barausse, Ilias Cholis, Juan Garcia-Bellido, Kelly Holley-Bockelmann, Scott A. Hughes, Bernard Kelly, Ely D. Kovetz, Tyson B. Littenberg, Jeffrey Livas, Guido Mueller, Priya Natarajan, David H. Shoemaker, Deirdre Shoemaker, Jeremy D. Schnittman, Michele Vallisneri, Nicolas Yunes

Abstract: Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as… ▽ More Low-frequency gravitational-wave astronomy can perform precision tests of general relativity and probe fundamental physics in a regime previously inaccessible. A space-based detector will be a formidable tool to explore gravity's role in the cosmos, potentially telling us if and where Einstein's theory fails and providing clues about some of the greatest mysteries in physics and astronomy, such as dark matter and the origin of the Universe. △ Less

Submitted 7 March, 2019; originally announced March 2019.

Comments: White Paper submitted on March 7, 2019 to Astro2020 (2020 Decadal Survey on Astronomy and Astrophysics)

Prompt Electromagnetic Transients from Binary Black Hole Mergers

Authors: Bernard J. Kelly, John G. Baker, Zachariah B. Etienne, Bruno Giacomazzo, Jeremy Schnittman

Abstract: Binary black hole (BBH) mergers provide a prime source for current and future interferometric GW observatories. Massive BBH mergers may often take place in plasma-rich environments, leading to the exciting possibility of a concurrent electromagnetic (EM) signal observable by traditional astronomical facilities. However, many critical questions about the generation of such counterparts remain unans… ▽ More Binary black hole (BBH) mergers provide a prime source for current and future interferometric GW observatories. Massive BBH mergers may often take place in plasma-rich environments, leading to the exciting possibility of a concurrent electromagnetic (EM) signal observable by traditional astronomical facilities. However, many critical questions about the generation of such counterparts remain unanswered. We explore mechanisms that may drive EM counterparts with magnetohydrodynamic simulations treating a range of scenarios involving equal-mass black-hole binaries immersed in an initially homogeneous fluid with uniform, orbitally aligned magnetic fields. We find that the time development of Poynting luminosity, which may drive jet-like emissions, is relatively insensitive to aspects of the initial configuration. In particular, over a significant range of initial values, the central magnetic field strength is effectively regulated by the gas flow to yield a Poynting luminosity of $10^{45}-10^{46} ρ_{-13} M_8^2 \, {\rm erg}\,{\rm s}^{-1}$, with BBH mass scaled to $M_8 \equiv M/(10^8 M_{\odot})$ and ambient density $ρ_{-13} \equiv ρ/(10^{-13} \, {\rm g} \, {\rm cm}^{-3})$. We also calculate the direct plasma synchrotron emissions processed through geodesic ray-tracing. Despite lensing effects and dynamics, we find the observed synchrotron flux varies little leading up to merger. △ Less

Submitted 16 January, 2018; v1 submitted 5 October, 2017; originally announced October 2017.

Comments: 22 pages, 21 figures; additional reference + clarifying text added to match published version

Journal ref: Phys. Rev. D 96, 123003 (2017)

Electromagnetic Chirps from Neutron Star-Black Hole Mergers

Authors: Jeremy D. Schnittman, Tito Dal Canton, Jordan Camp, David Tsang, Bernard J. Kelly

Abstract: We calculate the electromagnetic signal of a gamma-ray flare coming from the surface of a neutron star shortly before merger with a black hole companion. Using a new version of the Monte Carlo radiation transport code Pandurata that incorporates dynamic spacetimes, we integrate photon geodesics from the neutron star surface until they reach a distant observer or are captured by the black hole. The… ▽ More We calculate the electromagnetic signal of a gamma-ray flare coming from the surface of a neutron star shortly before merger with a black hole companion. Using a new version of the Monte Carlo radiation transport code Pandurata that incorporates dynamic spacetimes, we integrate photon geodesics from the neutron star surface until they reach a distant observer or are captured by the black hole. The gamma-ray light curve is modulated by a number of relativistic effects, including Doppler beaming and gravitational lensing. Because the photons originate from the inspiraling neutron star, the light curve closely resembles the corresponding gravitational waveform: a chirp signal characterized by a steadily increasing frequency and amplitude. We propose to search for these electromagnetic chirps using matched filtering algorithms similar to those used in LIGO data analysis. △ Less

Submitted 25 April, 2017; originally announced April 2017.

Comments: 13 pages, 5 figures, submitted to ApJ

Space Telescope and Optical Reverberation Mapping Project. VII. Understanding the UV anomaly in NGC 5548 with X-Ray Spectroscopy

Authors: S. Mathur, A. Gupta, K. Page, R. W. Pogge, Y. Krongold, M. R. Goad, S. M. Adams, M. D. Anderson, P. Arevalo, A. J. Barth, C. Bazhaw, T. G. Beatty, M. C. Bentz, A. Bigley, S. Bisogni, G. A. Borman, T. A. Boroson, M. C. Bottorff, W. N. Brandt, A. A. Breeveld, J. E. Brown, J. S. Brown, E. M. Cackett, G. Canalizo, M. T. Carini , et al. (125 additional authors not shown)

Abstract: During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuu… ▽ More During the Space Telescope and Optical Reverberation Mapping Project (STORM) observations of NGC 5548, the continuum and emission-line variability became de-correlated during the second half of the 6-month long observing campaign. Here we present Swift and Chandra X-ray spectra of NGC 5548 obtained as a part of the campaign. The Swift spectra show that excess flux (relative to a power-law continuum) in the soft X-ray band appears before the start of the anomalous emission-line behavior, peaks during the period of the anomaly, and then declines. This is a model-independent result suggesting that the soft excess is related to the anomaly. We divide the Swift data into on- and off-anomaly spectra to characterize the soft excess via spectral fitting. The cause of the spectral differences is likely due to a change in the intrinsic spectrum rather than being due to variable obscuration or partial covering. The Chandra spectra have lower signal-to-noise ratios, but are consistent with Swift data. Our preferred model of the soft excess is emission from an optically thick, warm Comptonizing corona, the effective optical depth of which increases during the anomaly. This model simultaneously explains all the three observations: the UV emission line flux decrease, the soft-excess increase, and the emission line anomaly. △ Less

Submitted 1 August, 2017; v1 submitted 20 April, 2017; originally announced April 2017.

Comments: ApJ in press. Replaced with the accepted version

Space Telescope and Optical Reverberation Mapping Project. V. Optical Spectroscopic Campaign and Emission-Line Analysis for NGC 5548

Authors: L. Pei, M. M. Fausnaugh, A. J. Barth, B. M. Peterson, M. C. Bentz, G. De Rosa, K. D. Denney, M. R. Goad, C. S. Kochanek, K. T. Korista, G. A. Kriss, R. W. Pogge, V. N. Bennert, M. Brotherton, K. I. Clubb, E. Dalla Bontà, A. V. Filippenko, J. E. Greene, C. J. Grier, M. Vestergaard, W. Zheng, Scott M. Adams, Thomas G. Beatty, A. Bigley, Jacob E. Brown , et al. (131 additional authors not shown)

Abstract: 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$β$ and He II $λ$4686 broad emission-line light curves lag that of the 5100 $Å$ optical continuum by… ▽ More 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$β$ and He II $λ$4686 broad emission-line light curves lag that of the 5100 $Å$ optical continuum by $4.17^{+0.36}_{-0.36}$ days and $0.79^{+0.35}_{-0.34}$ days, respectively. The H$β$ lag relative to the 1158 $Å$ 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$β$ 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$β$ and He II $λ$4686 emission lines to the driving continuum changed significantly halfway through the campaign, a phenomenon also observed for C IV, Ly $α$, 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$β$ 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. △ Less

Submitted 3 February, 2017; originally announced February 2017.

Comments: 22 pages, 13 figures, accepted to ApJ

Space Telescope and Optical Reverberation Mapping Project VI: reverberating Disk Models for NGC 5548

Authors: D. Starkey, Keith Horne, M. M. Fausnaugh, B. M. Peterson, M. C. Bentz, C. S. Kochanek, K. D. Denney, R. Edelson, M. R. Goad, G. De Rosa, M. D. Anderson, P. Arevalo, A. J. Barth, C. Bazhaw, G. A. Borman, T. A. Boroson, M. C. Bottorff, W. N. Brandt, A. A. Breeveld, E. M. Cackett, M. T. Carini, K. V. Croxall, D. M. Crenshaw, E. Dalla Bonta, A. De Lorenzo-Caceres , et al. (68 additional authors not shown)

Abstract: 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 respo… ▽ More 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 △ Less

Submitted 24 November, 2016; v1 submitted 18 November, 2016; originally announced November 2016.

Comments: V2: Oops wrong title! V1: Accepted for publication in ApJ, 20 Pages, 11 Figures

Space Telescope and Optical Reverberation Mapping Project. IV. Anomalous behavior of the broad ultraviolet emission lines in NGC 5548

Authors: M. R. Goad, K. T. Korista, G. De Rosa, G. A. Kriss, R. Edelson, A. J. Barth, G. J. Ferland, C. S. Kochanek, H. Netzer, B. M. Peterson, M. C. Bentz, S. Bisogni, D. M. Crenshaw, K. D. Denney, J. Ely, M. M. Fausnaugh, C. J. Grier, A. Gupta, K. D. Horne, J. Kaastra, A. Pancoast, L. Pei, R. W. Pogge, A. Skielboe, D. Starkey , et al. (77 additional authors not shown)

Abstract: During an intensive Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) UV monitoring campaign of the Seyfert~1 galaxy NGC 5548 performed from 2014 February to July, the normally highly correlated far-UV continuum and broad emission-line variations decorrelated for ~60 to 70 days, starting ~75 days after the first HST/COS observation. Following this anomalous state, the flux and variabi… ▽ More During an intensive Hubble Space Telescope (HST) Cosmic Origins Spectrograph (COS) UV monitoring campaign of the Seyfert~1 galaxy NGC 5548 performed from 2014 February to July, the normally highly correlated far-UV continuum and broad emission-line variations decorrelated for ~60 to 70 days, starting ~75 days after the first HST/COS observation. Following this anomalous state, the flux and variability of the broad emission lines returned to a more normal state. This transient behavior, characterised by significant deficits in flux and equivalent width of the strong broad UV emission lines, is the first of its kind to be unambiguously identified in an active galactic nucleus reverberation mapping campaign. The largest corresponding emission-line flux deficits occurred for the high-ionization collisionally excited lines, C IV and Si IV(+O IV]), and also He II(+O III]), while the anomaly in Ly-alpha was substantially smaller. This pattern of behavior indicates a depletion in the flux of photons with E_{\rm ph} > 54 eV, relative to those near 13.6 eV. We suggest two plausible mechanisms for the observed behavior: (i) temporary obscuration of the ionizing continuum incident upon BLR clouds by a moving veil of material lying between the inner accretion disk and inner BLR, perhaps resulting from an episodic ejection of material from the disk, or (ii) a temporary change in the intrinsic ionizing continuum spectral energy distribution resulting in a deficit of ionizing photons with energies > 54 eV, possibly due to a transient restructuring of the Comptonizing atmosphere above the disk. Current evidence appears to favor the latter explanation. △ Less

Submitted 29 March, 2016; originally announced March 2016.

Comments: 12 pages, 2 figures, accepted for publication in ApJ, Friday 25th March 2016. A movie of the anomalous emission-line behavior can be found in the ancilliary documents

Space Telescope and Optical Reverberation Mapping Project. III. Optical Continuum Emission and Broad-Band Time Delays in NGC 5548

Authors: M. M. Fausnaugh, K. D. Denney, A. J. Barth, M. C. Bentz, M. C. Bottorff, M. T. Carini, K. V. Croxall, G. De Rosa, M. R. Goad, Keith Horne, M. D. Joner, S. Kaspi, M. Kim, S. A. Klimanov, C. S. Kochanek, D. C. Leonard, H. Netzer, B. M. Peterson, K. Schnulle, S. G. Sergeev, M. Vestergaard, W. -K. Zheng, Y. Zu, M. D. Anderson, P. Arevalo , et al. (72 additional authors not shown)

Abstract: 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 conti… ▽ More 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\,Å to the $z$ band ($\sim\!9160$\,Å). 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}\,$λ1640$ and $λ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 $τ\propto λ^{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. △ Less

Submitted 29 February, 2016; v1 submitted 19 October, 2015; originally announced October 2015.

Comments: 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

The Time Domain Spectroscopic Survey: Variable Object Selection and Anticipated Results

Authors: Eric Morganson, Paul J. Green, Scott F. Anderson, John J. Ruan, Adam D. Myers, Michael Eracleous, Brandon Kelly, Carlos Badenes, Eduardo Banados, Michael R. Blanton, Matthew A. Bershady, Jura Borissova, William Nielsen Brandt, William S. Burgett, Kenneth Chambers, Peter W. Draper, James R. A. Davenport, Heather Flewelling, Peter Garnavich, Suzanne L. Hawley, Klaus W. Hodapp, Jedidah C. Isler, Nick Kaiser, Karen Kinemuchi, Rolf P. Kudritzki , et al. (15 additional authors not shown)

Abstract: We present the selection algorithm and anticipated results for the Time Domain Spectroscopic Survey (TDSS). TDSS is an SDSS-IV eBOSS subproject that will provide initial identification spectra of approximately 220,000 luminosity-variable objects (variable stars and AGN) across 7,500 square degrees selected from a combination of SDSS and multi-epoch Pan-STARRS1 photometry. TDSS will be the largest… ▽ More We present the selection algorithm and anticipated results for the Time Domain Spectroscopic Survey (TDSS). TDSS is an SDSS-IV eBOSS subproject that will provide initial identification spectra of approximately 220,000 luminosity-variable objects (variable stars and AGN) across 7,500 square degrees selected from a combination of SDSS and multi-epoch Pan-STARRS1 photometry. TDSS will be the largest spectroscopic survey to explicitly target variable objects, avoiding pre-selection on the basis of colors or detailed modeling of specific variability characteristics. Kernel Density Estimate (KDE) analysis of our target population performed on SDSS Stripe 82 data suggests our target sample will be 95% pure (meaning 95% of objects we select have genuine luminosity variability of a few magnitudes or more). Our final spectroscopic sample will contain roughly 135,000 quasars and 85,000 stellar variables, approximately 4,000 of which will be RR Lyrae stars which may be used as outer Milky Way probes. The variability-selected quasar population has a smoother redshift distribution than a color-selected sample, and variability measurements similar to those we develop here may be used to make more uniform quasar samples in large surveys. The stellar variable targets are distributed fairly uniformly across color space, indicating that TDSS will obtain spectra for a wide variety of stellar variables including pulsating variables, stars with significant chromospheric activity, cataclysmic variables and eclipsing binaries. TDSS will serve as a pathfinder mission to identify and characterize the multitude of variable objects that will be detected photometrically in even larger variability surveys such as LSST. △ Less

Submitted 4 May, 2015; originally announced May 2015.

Comments: 22 pages, 18 figures, accepted by ApJ

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

Authors: G. De Rosa, B. M. Peterson, J. Ely, G. A. Kriss, D. M. Crenshaw, Keith Horne, K. T. Korista, H. Netzer, R. W. Pogge, P. Arevalo, A. J. Barth, M. C. Bentz, W. N. Brandt, A. A. Breeveld, B. J. Brewer, E. Dalla Bonta, A. De Lorenzo-Caceres, K. D. Denney, M. Dietrich, R. Edelson, P. A. Evans, M. M. Fausnaugh, N. Gehrels, J. M. Gelbord, M. R. Goad , et al. (25 additional authors not shown)

Abstract: 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 f… ▽ More 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. △ Less

Submitted 27 April, 2015; v1 submitted 23 January, 2015; originally announced January 2015.

Comments: 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

Journal ref: ApJ 806 (2015) 128

Space Telescope and Optical Reverberation Mapping Project. II. Swift and HST Reverberation Mapping of the Accretion Disk of NGC 5548

Authors: R. Edelson, J. M. Gelbord, K. Horne, I. M. McHardy, B. M. Peterson, P. Arevalo, A. A. Breeveld, G. De Rosa, P. A. Evans, M. R. Goad, G. A. Kriss, W. N. Brandt, N. Gehrels, D. Grupe, J. A. Kennea, C. S. Kochanek, J. A. Nousek, I. Papadakis, M. Siegel, D. Starkey, P. Uttley, S. Vaughan, S. Young, A. J. Barth, M. C. Bentz , et al. (25 additional authors not shown)

Abstract: 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 <0.5 day. Approximately daily HST UV sampling was also obtained. The UV/optical light curves show strong correlations (r_max = 0.57 - 0.90) and t… ▽ More 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 <0.5 day. Approximately daily HST UV sampling was also obtained. The UV/optical light curves show strong correlations (r_max = 0.57 - 0.90) and the clearest measurement to date of interband lags. These lags are well-fit by a τpropto λ^4/3 wavelength dependence, with a normalization that indicates an unexpectedly large disk radius of 0.35 +/- 0.05 lt-day at 1367 A, assuming a simple face-on model. The U-band shows a marginally larger lag than expected from the fit and surrounding bands, which could be due to Balmer continuum emission from the broad-line region as suggested by Korista and Goad. The UV/X-ray correlation is weaker (r_max < 0.45) and less consistent over time. This indicates that while Swift is beginning to measure UV/optical lags in general agreement with accretion disk theory (although the derived size is larger than predicted), the relationship with X-ray variability is less well understood. Combining this accretion disk size estimate with those from quasar microlensing studies suggests that AGN disk sizes scale approximately linearly with central black hole mass over a wide range of masses. △ Less

Submitted 7 April, 2015; v1 submitted 23 January, 2015; originally announced January 2015.

Comments: 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

Journal ref: ApJ 806 (2015) 129

Data Mining for Gravitationally Lensed Quasars

Authors: Adriano Agnello, Brandon C. Kelly, Tommaso Treu, Philip J. Marshall

Abstract: Gravitationally lensed (GL) quasars are brighter than their unlensed counterparts and produce images with distinctive morphological signatures. Past searches and target selection algorithms, in particular the Sloan Quasar Lens Search (SQLS), have relied on basic morphological criteria, which were applied to samples of bright, spectroscopically confirmed quasars. The SQLS techniques are not suffici… ▽ More Gravitationally lensed (GL) quasars are brighter than their unlensed counterparts and produce images with distinctive morphological signatures. Past searches and target selection algorithms, in particular the Sloan Quasar Lens Search (SQLS), have relied on basic morphological criteria, which were applied to samples of bright, spectroscopically confirmed quasars. The SQLS techniques are not sufficient for searching into new surveys (e.g. DES, PS1, LSST), because spectroscopic information is not readily available and the large data volume requires higher purity in target/candidate selection. We carry out a systematic exploration of machine learning techniques and demonstrate that a two step strategy can be highly effective. In the first step we use catalog-level information ($griz$+WISE magnitudes, second moments) to preselect targets, using artificial neural networks. The accepted targets are then inspected with pixel-by-pixel pattern recognition algorithms (Gradient-Boosted Trees), to form a final set of candidates. The results from this procedure can be used to further refine the simpler SQLS algorithms, with a twofold (or threefold) gain in purity and the same (or $80\%$) completeness at target-selection stage, or a purity of $70\%$ and a completeness of $60\%$ after the candidate-selection step. Simpler photometric searches in $griz$+WISE based on colour cuts would provide samples with $7\%$ purity or less. Our technique is extremely fast, as a list of candidates can be obtained from a stage III experiment (e.g. DES catalog/database) in {a few} CPU hours. The techniques are easily extendable to Stage IV experiments like LSST with the addition of time domain information. △ Less

Submitted 16 October, 2014; originally announced October 2014.

Comments: MNRAS submitted, positive report received; u-band selection to be added in new version; 19 pages, 10 figures

Discovery of a ~5 day characteristic timescale in the Kepler power spectrum of Zw 229-15

Authors: Rick Edelson, Simon Vaughan, Matt Malkan, Brandon Kelly, Krista Smith, Padi Boyd, Richard Mushotzky

Abstract: We present time series analyses of the full Kepler dataset of Zw 229-15. This Kepler light curve --- with a baseline greater than three years, composed of virtually continuous, evenly sampled 30-minute measurements --- is unprecedented in its quality and precision. We utilize two methods of power spectral analysis to investigate the optical variability and search for evidence of a bend frequency a… ▽ More We present time series analyses of the full Kepler dataset of Zw 229-15. This Kepler light curve --- with a baseline greater than three years, composed of virtually continuous, evenly sampled 30-minute measurements --- is unprecedented in its quality and precision. We utilize two methods of power spectral analysis to investigate the optical variability and search for evidence of a bend frequency associated with a characteristic optical variability timescale. Each method yields similar results. The first interpolates across data gaps to use the standard Fourier periodogram. The second, using the CARMA-based time-domain modeling technique of Kelly et al. (2014), does not need evenly-sampled data. Both methods find excess power at high frequencies that may be due to Kepler instrumental effects. More importantly both also show strong bends (Δα ~ 2) at timescales of ~5 days, a feature similar to those seen in the X-ray PSDs of AGN but never before in the optical. This observed ~5 day timescale may be associated with one of several physical processes potentially responsible for the variability. A plausible association could be made with light-crossing, dynamical or thermal timescales, depending on the assumed value of the accretion disk size and on unobserved disk parameters such as α and H/R. This timescale is not consistent with the viscous timescale, which would be years in a ~10^7 Solar mass AGN such as Zw 229-15. However there must be a second bend on long (>~1 year) timescales, and that feature could be associated with the viscous timescale. △ Less

Submitted 4 September, 2014; originally announced September 2014.

Comments: 10 pages, 5 figures, 1 table. To appear in the Astrophysical Journal, Part 1

Journal ref: ApJ, 795, 2 (2014)

The Sloan Digital Sky Survey Reverberation Mapping Project: Technical Overview

Authors: Yue Shen, W. N. Brandt, Kyle S. Dawson, Patrick B. Hall, Ian D. McGreer, Scott F. Anderson, Yuguang Chen, Kelly D. Denney, Sarah Eftekharzadeh, Xiaohui Fan, Yang Gao, Paul J. Green, Jenny E. Greene, Luis C. Ho, Keith Horne, Linhua Jiang, Brandon C. Kelly, Karen Kinemuchi, Christopher S. Kochanek, Isabelle Pâris, Christina M. Peters, Bradley M. Peterson, Patrick Petitjean, Kara Ponder, Gordon T. Richards , et al. (14 additional authors not shown)

Abstract: The Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM) is a dedicated multi-object RM experiment that has spectroscopically monitored a sample of 849 broad-line quasars in a single 7 deg$^2$ field with the SDSS-III BOSS spectrograph. The RM quasar sample is flux-limited to i_psf=21.7 mag, and covers a redshift range of 0.1<z<4.5. Optical spectroscopy was performed during 2014 Jan-Jul… ▽ More The Sloan Digital Sky Survey Reverberation Mapping project (SDSS-RM) is a dedicated multi-object RM experiment that has spectroscopically monitored a sample of 849 broad-line quasars in a single 7 deg$^2$ field with the SDSS-III BOSS spectrograph. The RM quasar sample is flux-limited to i_psf=21.7 mag, and covers a redshift range of 0.1<z<4.5. Optical spectroscopy was performed during 2014 Jan-Jul dark/grey time, with an average cadence of ~4 days, totaling more than 30 epochs. Supporting photometric monitoring in the g and i bands was conducted at multiple facilities including the CFHT and the Steward Observatory Bok telescopes in 2014, with a cadence of ~2 days and covering all lunar phases. The RM field (RA, DEC=14:14:49.00, +53:05:00.0) lies within the CFHT-LS W3 field, and coincides with the Pan-STARRS 1 (PS1) Medium Deep Field MD07, with three prior years of multi-band PS1 light curves. The SDSS-RM 6-month baseline program aims to detect time lags between the quasar continuum and broad line region (BLR) variability on timescales of up to several months (in the observed frame) for ~10% of the sample, and to anchor the time baseline for continued monitoring in the future to detect lags on longer timescales and at higher redshift. SDSS-RM is the first major program to systematically explore the potential of RM for broad-line quasars at z>0.3, and will investigate the prospects of RM with all major broad lines covered in optical spectroscopy. SDSS-RM will provide guidance on future multi-object RM campaigns on larger scales, and is aiming to deliver more than tens of BLR lag detections for a homogeneous sample of quasars. We describe the motivation, design and implementation of this program, and outline the science impact expected from the resulting data for RM and general quasar science. △ Less

Submitted 25 August, 2014; originally announced August 2014.

Comments: 25 pages, submitted to ApJS; project website at http://www.sdssrm.org

Improved Moving Puncture Gauge Conditions for Compact Binary Evolutions

Authors: Zachariah B. Etienne, John G. Baker, Vasileios Paschalidis, Bernard J. Kelly, Stuart L. Shapiro

Abstract: Robust gauge conditions are critically important to the stability and accuracy of numerical relativity (NR) simulations involving compact objects. Most of the NR community use the highly robust---though decade-old---moving-puncture (MP) gauge conditions for such simulations. It has been argued that in binary black hole (BBH) evolutions adopting this gauge, noise generated near adaptive-mesh-refine… ▽ More Robust gauge conditions are critically important to the stability and accuracy of numerical relativity (NR) simulations involving compact objects. Most of the NR community use the highly robust---though decade-old---moving-puncture (MP) gauge conditions for such simulations. It has been argued that in binary black hole (BBH) evolutions adopting this gauge, noise generated near adaptive-mesh-refinement (AMR) boundaries does not converge away cleanly with increasing resolution, severely limiting gravitational waveform accuracy at computationally feasible resolutions. We link this noise to a sharp (short-wavelength), initial outgoing gauge wave crossing into progressively lower resolution AMR grids, and present improvements to the standard MP gauge conditions that focus on stretching, smoothing, and more rapidly settling this outgoing wave. Our best gauge choice greatly reduces gravitational waveform noise during inspiral, yielding less fluctuation in convergence order and $\sim 40%$ lower waveform phase and amplitude errors at typical resolutions. Noise in other physical quantities of interest is also reduced, and constraint violations drop by more than an order of magnitude. We expect these improvements will carry over to simulations of all types of compact binary systems, as well as other $N$+1 formulations of gravity for which MP-like gauge conditions can be chosen. △ Less

Submitted 5 October, 2014; v1 submitted 25 April, 2014; originally announced April 2014.

Comments: 25 pages, 16 figures, 2 tables. Matches published version

Journal ref: Phys. Rev. D 90, 064032 (2014)

Stochastic Modeling of the Fermi/LAT Gamma-ray Blazar Variability

Authors: Malgorzata A. Sobolewska, Aneta Siemiginowska, Brandon C. Kelly, Krzysztof Nalewajko

Abstract: We study the gamma-ray variability of 13 blazars observed with the Fermi Large Area Telescope (LAT). These blazars have the most complete light curves collected during the first 4 years of the Fermi sky survey. We model them with the Ornstein-Uhlenbeck (OU) process or a mixture of the OU processes. The OU process has power spectral density (PSD) proportional to 1/f^alpha with alpha changing at a c… ▽ More We study the gamma-ray variability of 13 blazars observed with the Fermi Large Area Telescope (LAT). These blazars have the most complete light curves collected during the first 4 years of the Fermi sky survey. We model them with the Ornstein-Uhlenbeck (OU) process or a mixture of the OU processes. The OU process has power spectral density (PSD) proportional to 1/f^alpha with alpha changing at a characteristic time scale, tau_0, from 0 (tau>>tau_0) to 2 (tau<<tau_0). The PSD of the mixed OU process has two characteristic time scales and an additional intermediate region with 0<alpha<2. We show that the OU model provides a good description of the Fermi/LAT light curves of three blazars in our sample. For the first time we constrain a characteristic gamma-ray time scale of variability in two BL Lac sources, 3C 66A and PKS 2155-304 (tau_0=25 day and tau_0=43 day, respectively, in the observer's frame), which are longer than the soft X-ray time scales detected in blazars and Seyfert galaxies. We find that the mixed OU process approximates the light curves of the remaining 10 blazars better than the OU process. We derive limits on their long and short characteristic time scales, and infer that their Fermi/LAT PSDs resemble a power-law function. We constrain the PSD slopes for all but one source in the sample. We find hints for sub-hour Fermi/LAT variability in four flat spectrum radio quasars. We discuss the implications of our results for theoretical models of blazar variability. △ Less

Submitted 20 March, 2014; originally announced March 2014.

Comments: Accepted for publication in ApJ, 16 pages, 11 figures. High resolution figures at http://thedata.harvard.edu/dvn/dv/blazars

X-ray constraints on the local supermassive black hole occupation fraction

Authors: Brendan P. Miller, Elena Gallo, Jenny E. Greene, Brandon C. Kelly, Tommaso Treu, Jong-Hak Woo, Vivienne Baldassare

Abstract: Distinct seed formation mechanisms are imprinted upon the fraction of dwarf galaxies currently containing a central supermassive black hole. Seeding by Pop III remnants is expected to produce a higher occupation fraction than is generated with direct gas collapse precursors. Chandra observations of nearby early-type galaxies can directly detect even low-level supermassive black hole activity, and… ▽ More Distinct seed formation mechanisms are imprinted upon the fraction of dwarf galaxies currently containing a central supermassive black hole. Seeding by Pop III remnants is expected to produce a higher occupation fraction than is generated with direct gas collapse precursors. Chandra observations of nearby early-type galaxies can directly detect even low-level supermassive black hole activity, and the active fraction immediately provides a firm lower limit to the occupation fraction. Here, we use the volume-limited AMUSE surveys of ~200 optically-selected early-type galaxies to characterize simultaneously, for the first time, the occupation fraction and the scaling of nuclear X-ray luminosity with stellar mass, accounting for intrinsic scatter, measurement uncertainties, and X-ray limits. For early-type galaxies with log(M_star/M_sun)<10, we obtain a lower limit to the occupation fraction of >20% (at 95% confidence), but full occupation cannot be excluded. The preferred dependence of log(L_X) upon log(M_star) has a slope of about 0.7-0.8, consistent with the "downsizing" trend previously identified from the AMUSE dataset, and a uniform Eddington efficiency is disfavored at ~2 sigma. We provide guidelines for the future precision with which these parameters may be refined with larger or more sensitive samples. △ Less

Submitted 23 October, 2014; v1 submitted 17 March, 2014; originally announced March 2014.

Comments: Changes from v1: slightly improved galaxy distances (analysis rerun with nearly identical results), added table with sample properties, expanded discussion of methodology, updated references. Resubmitted to ApJ

Flexible and Scalable Methods for Quantifying Stochastic Variability in the Era of Massive Time-Domain Astronomical Data Sets

Authors: Brandon C. Kelly, Andrew C. Becker, Malgosia Sobolewska, Aneta Siemiginowska, Phil Uttley

Abstract: We present the use of continuous-time autoregressive moving average (CARMA) models as a method for estimating the variability features of a light curve, and in particular its power spectral density (PSD). CARMA models fully account for irregular sampling and measurement errors, making them valuable for quantifying variability, forecasting and interpolating light curves, and for variability-based c… ▽ More We present the use of continuous-time autoregressive moving average (CARMA) models as a method for estimating the variability features of a light curve, and in particular its power spectral density (PSD). CARMA models fully account for irregular sampling and measurement errors, making them valuable for quantifying variability, forecasting and interpolating light curves, and for variability-based classification. We show that the PSD of a CARMA model can be expressed as a sum of Lorentzian functions, which makes them extremely flexible and able to model a broad range of PSDs. We present the likelihood function for light curves sampled from CARMA processes, placing them on a statistically rigorous foundation, and we present a Bayesian method to infer the probability distribution of the PSD given the measured lightcurve. Because calculation of the likelihood function scales linearly with the number of data points, CARMA modeling scales to current and future massive time-domain data sets. We conclude by applying our CARMA modeling approach to light curves for an X-ray binary, two AGN, a long-period variable star, and an RR-Lyrae star, in order to illustrate their use, applicability, and interpretation. △ Less

Submitted 7 May, 2014; v1 submitted 24 February, 2014; originally announced February 2014.

Comments: 18 pages, 16 figures, in press at ApJ. Revised to match accepted version. Code for utilizing the CARMA models is at https://github.com/bckelly80/carma_pack

The Luminosity Function at z~8 from 97 Y-band dropouts: Inferences About Reionization

Authors: Kasper B. Schmidt, Tommaso Treu, Michele Trenti, Larry D. Bradley, Brandon C. Kelly, Pascal A. Oesch, Benne W. Holwerda, J. Michael Shull, Massimo Stiavelli

Abstract: [Abbreviated] We present the largest search to date for $z\sim8$ Lyman break galaxies (LBGs) based on 350 arcmin$^2$ of HST observations in the V-, Y-, J- and H-bands from the Brightest of Reionizing Galaxies (BoRG) survey. The BoRG dataset includes $\sim$50 arcmin$^2$ of new data and deeper observations of two previous BoRG pointings, from which we present 9 new $z\sim8$ LBG candidates, bringing… ▽ More [Abbreviated] We present the largest search to date for $z\sim8$ Lyman break galaxies (LBGs) based on 350 arcmin$^2$ of HST observations in the V-, Y-, J- and H-bands from the Brightest of Reionizing Galaxies (BoRG) survey. The BoRG dataset includes $\sim$50 arcmin$^2$ of new data and deeper observations of two previous BoRG pointings, from which we present 9 new $z\sim8$ LBG candidates, bringing the total number of BoRG LBGs to 38 with $25.5\leqslant m_{J} \leqslant 27.6$ (AB system). We introduce a new Bayesian formalism for estimating the galaxy luminosity function (LF), which does not require binning (and thus smearing) of the data and includes a likelihood based on the formally correct binomial distribution as opposed to the often used approximate Poisson distribution. We demonstrate the utility of the new method on a sample of $97$ LBGs that combines the bright BoRG galaxies with the fainter sources published in Bouwens et al. (2012) from the HUDF and ERS programs. We show that the $z\sim8$ LF is well described by a Schechter function with a characteristic magnitude $M^\star = -20.15^{+0.29}_{-0.38}$, a faint-end slope of $α= -1.87^{+0.26}_{-0.26}$, and a number density of $\log_{10} φ^\star [\textrm{Mpc}^{-3}] = -3.24^{+0.25}_{-0.24}$. Integrated down to $M=-17.7$ this LF yields a luminosity density, $\log_{10} ε[\textrm{erg}/\textrm{s/Hz/Mpc}^{3}] = 25.52^{+0.05}_{-0.05}$. Our LF analysis is consistent with previously published determinations within 1$σ$. We discuss the implication of our study for the physics of reionization. By assuming theoretically motivated priors on the clumping factor and the photon escape fraction we show that the UV LF from galaxy samples down to $M=-17.7$ can ionize only 10-50% of the neutral hydrogen at $z\sim8$. Full reionization would require extending the LF down to $M=-15$. △ Less

Submitted 24 March, 2014; v1 submitted 17 February, 2014; originally announced February 2014.

Comments: Accepted for publication in ApJ, 22 pages, 15 figures

Jet Emission in Young Radio Sources: a Fermi-LAT Gamma-ray View

Authors: Giulia Migliori, Aneta Siemiginowska, Brandon C. Kelly, Lukasz Stawarz, Annalisa Celotti, Mitchell C. Begelman

Abstract: We investigate the contribution of the beamed jet component to the high energy emission in young and compact extragalactic radio sources, focusing for the first time on the gamma-ray band. We derive predictions on the gamma-ray luminosities associated with the relativistic jet assuming a leptonic radiative model. The high energy emission is produced via Compton scattering by the relativistic elect… ▽ More We investigate the contribution of the beamed jet component to the high energy emission in young and compact extragalactic radio sources, focusing for the first time on the gamma-ray band. We derive predictions on the gamma-ray luminosities associated with the relativistic jet assuming a leptonic radiative model. The high energy emission is produced via Compton scattering by the relativistic electrons in a spherical region at the considered scales ($\lesssim$10 kpc). Simulations show a wide range of gamma-ray luminosities, with intensities up to $\sim10^{46}-10^{48}$ erg s$^{-1}$ depending on the assumed jet parameters. We find a highly linear relation between the simulated X-ray and gamma-ray luminosities that can be used to select candidates for a gamma-ray detection. We compare the simulated luminosity distributions in the radio, X-ray and gamma-ray regimes with observations for the largest sample of X-ray detected young radio quasars. Our analysis of $\sim$4-year Fermi Large Area Telescope (LAT) data does not give any statistically significant detection. However, the majority of the model-predicted gamma-ray fluxes for the sample are near or below the current Fermi-LAT flux threshold and compatible with the derived upper limits. Our study gives constraints on the minimum jet power ($L_{jet,kin}/L_{disk}>0.01$), on a potential jet contribution to the X-ray emission in the most compact sources ($\lesssim1$ kpc) and on the particles to magnetic field energy density ratio in broad agreement with equipartition assumption. △ Less

Submitted 29 November, 2013; originally announced November 2013.

Comments: 21 pages, 8 pages, accepted for publication in ApJ

An 8 hour characteristic time-scale in submillimetre light curves of Sagittarius A*

Authors: Jason Dexter, Brandon Kelly, Geoffrey C. Bower, Daniel P. Marrone, Jordan Stone, Richard Plambeck, Sheperd S. Doeleman

Abstract: We compile and analyse long term (~10 year) submillimetre (1.3, 0.87, 0.43 mm, submm) wavelength light curves of the Galactic centre black hole, Sagittarius A*. The 0.87 and 0.43 mm data are taken from the literature, while the majority of the 1.3 mm light curve is from previously unpublished SMA and CARMA data. We use Monte Carlo simulations to show that on minute to few hour time-scales the vari… ▽ More We compile and analyse long term (~10 year) submillimetre (1.3, 0.87, 0.43 mm, submm) wavelength light curves of the Galactic centre black hole, Sagittarius A*. The 0.87 and 0.43 mm data are taken from the literature, while the majority of the 1.3 mm light curve is from previously unpublished SMA and CARMA data. We use Monte Carlo simulations to show that on minute to few hour time-scales the variability is consistent with a red noise process with a 230 GHz power spectrum slope of 2.3+0.8-0.6 at 95% confidence. The light curve is de-correlated (white noise) on very long (month to year) times. In order to identify the transition time between red and white noise, we model the light curves as a stochastic damped random walk process. The models allow a quantitative estimate of this physical characteristic time-scale of 8-4+3 hours at 230 GHz at 95% confidence, with consistent results at 345 and 690 GHz. This corresponds to ~10 orbital times or ~1 inflow (viscous) time at R = 3 Rs, a typical radius producing the 230 GHz emission as measured by very long baseline interferometry and found in theoretical accretion flow and jet models. This time-scale is significantly shorter (longer) than those measured at radio (near-infrared, NIR) wavelengths, and is marginally inconsistent with the same variability mechanism operating in the submm and NIR for the expected t ~ R^3/2 scaling. It is crudely consistent with the analogous time-scale inferred in studies of quasar optical light curves after accounting for the difference in emission radius. We find evidence that the submm variability persists at least down to the ISCO, if not the event horizon. These results can be compared quantitatively with similar analyses at different wavebands to test for connections between the variability mechanisms, and with light curves from theoretical models of accreting black holes. △ Less

Submitted 27 August, 2013; originally announced August 2013.

Comments: 12 pages, 9 figures, submitted to MNRAS 6/2013

Quasar Black Hole Mass Estimates in the Era of Time Domain Astronomy

Authors: Brandon C. Kelly, Tommaso Treu, Matthew Malkan, Anna Pancoast, Jong-Hak Woo

Abstract: We investigate the dependence of the normalization of the high-frequency part of the X-ray and optical power spectral densities (PSD) on black hole mass for a sample of 39 active galactic nuclei (AGN) with black hole masses estimated from reverberation mapping or dynamical modeling. We obtained new Swift observations of PG 1426+015, which has the largest estimated black hole mass of the AGN in our… ▽ More We investigate the dependence of the normalization of the high-frequency part of the X-ray and optical power spectral densities (PSD) on black hole mass for a sample of 39 active galactic nuclei (AGN) with black hole masses estimated from reverberation mapping or dynamical modeling. We obtained new Swift observations of PG 1426+015, which has the largest estimated black hole mass of the AGN in our sample. We develop a novel statistical method to estimate the PSD from a lightcurve of photon counts with arbitrary sampling, eliminating the need to bin a lightcurve to achieve Gaussian statistics, and we use this technique to estimate the X-ray variability parameters for the faint AGN in our sample. We find that the normalization of the high-frequency X-ray PSD is inversely proportional to black hole mass. We discuss how to use this scaling relationship to obtain black hole mass estimates from the short time-scale X-ray variability amplitude with precision ~ 0.38 dex. The amplitude of optical variability on time scales of days is also anti-correlated with black hole mass, but with larger scatter. Instead, the optical variability amplitude exhibits the strongest anti-correlation with luminosity. We conclude with a discussion of the implications of our results for estimating black hole mass from the amplitude of AGN variability. △ Less

Submitted 19 July, 2013; originally announced July 2013.

Comments: 19 pages, 10 figures, emulateapj format, submitted to ApJ

Indications of a sub-linear and non-universal Kennicutt-Schmidt relationship

Authors: Rahul Shetty, Brandon C. Kelly, Nurur Rahman, Frank Bigiel, Alberto D. Bolatto, Paul C. Clark, Ralf S. Klessen, Lukas K. Konstandin

Abstract: We estimate the parameters of the Kennicutt-Schmidt (KS) relationship, linking the star formation rate (Sigma_SFR) to the molecular gas surface density (Sigma_mol), in the STING sample of nearby disk galaxies using a hierarchical Bayesian method. This method rigorously treats measurement uncertainties, and provides accurate parameter estimates for both individual galaxies and the entire population… ▽ More We estimate the parameters of the Kennicutt-Schmidt (KS) relationship, linking the star formation rate (Sigma_SFR) to the molecular gas surface density (Sigma_mol), in the STING sample of nearby disk galaxies using a hierarchical Bayesian method. This method rigorously treats measurement uncertainties, and provides accurate parameter estimates for both individual galaxies and the entire population. Assuming standard conversion factors to estimate Sigma_SFR and Sigma_mol from the observations, we find that the KS parameters vary between galaxies, indicating that no universal relationship holds for all galaxies. The KS slope of the whole population is 0.76, with the 2sigma range extending from 0.58 to 0.94. These results imply that the molecular gas depletion time is not constant, but varies from galaxy to galaxy, and increases with the molecular gas surface density. Therefore, other galactic properties besides just Sigma_mol affect Sigma_SFR, such as the gas fraction or stellar mass. The non-universality of the KS relationship indicates that a comprehensive theory of star formation must take into account additional physical processes that may vary from galaxy to galaxy. △ Less

Submitted 10 October, 2013; v1 submitted 12 June, 2013; originally announced June 2013.

Comments: 7 pages, 2 figures, 1 table. Updated to match MNRAS accepted version

Empirical Links between XRB and AGN accretion using the complete z<0.4 spectroscopic CSC/SDSS Catalog

Authors: Markos Trichas, Paul Green, Anca Constantin, Tom Aldcroft, Malgosia Sobolewska, Ashley K. Hyde, Hongyan Zhou, Dong-Woo Kim, Daryl Haggard, Brandon Kelly, Eleni Kalfountzou

Abstract: Striking similarities have been seen between accretion signatures of Galactic X-ray binary (XRB) systems and active galactic nuclei (AGN). XRB spectral states show a V-shaped correlation between X-ray spectral hardness and Eddington ratio as they vary, and some AGN samples reveal a similar trend, implying analogous processes at vastly larger masses and timescales. To further investigate the analog… ▽ More Striking similarities have been seen between accretion signatures of Galactic X-ray binary (XRB) systems and active galactic nuclei (AGN). XRB spectral states show a V-shaped correlation between X-ray spectral hardness and Eddington ratio as they vary, and some AGN samples reveal a similar trend, implying analogous processes at vastly larger masses and timescales. To further investigate the analogies, we have matched 617 sources from the Chandra Source Catalog to SDSS spectroscopy, and uniformly measured both X-ray and optical spectral characteristics across a broad range of AGN and galaxy types. We provide useful tabulations of X-ray spectral slope for broad and narrow line AGN, star-forming and passive galaxies and composite systems, also updating relationships between optical (Ha and [OIII]) line emission and X-ray luminosity. We further fit broadband spectral energy distributions with a variety of templates to estimate bolometric luminosity. Our results confirm a significant trend in AGN between X-ray spectral hardness and Eddington ratio expressed in X-ray luminosity, albeit with significant dispersion. The trend is not significant when expressed in the full bolometric or template-estimated AGN luminosity. We also confirm a relationship between the X-ray/optical spectral slope aox, and Eddington ratio, but it may not follow the trend predicted by analogy with XRB accretion states. △ Less

Submitted 3 June, 2013; originally announced June 2013.

Comments: 12 pages, 10 figures, submitted to ApJ. Full online table will be provided via ApJ

Stacked reverberation lags at high redshift

Authors: S. Fine, T. Shanks, P. Green, B. C. Kelly, S. M. Croom, R. L. Webster, E. Berger, R. Chornock, W. S. Burgett, K. C. Chambers, N. Kaise, P. A. Price

Abstract: Over the past 20years reverberation mapping has proved one of the most successful techniques for studying the local (<1pc) environment of super-massive black holes that drive active galactic nuclei. Key successes of reverberation mapping have been direct black hole mass estimates, the radius-luminosity relation for the Hbeta line and the calibration of single-epoch mass estimators commonly employe… ▽ More Over the past 20years reverberation mapping has proved one of the most successful techniques for studying the local (<1pc) environment of super-massive black holes that drive active galactic nuclei. Key successes of reverberation mapping have been direct black hole mass estimates, the radius-luminosity relation for the Hbeta line and the calibration of single-epoch mass estimators commonly employed up to z~7. However, observing constraints mean that few studies have been successful at z>0.1, or for the more-luminous quasars that make up the majority of current spectroscopic samples, or for the rest-frame ultra-violet emission lines available in optical spectra of z>0.5 objects. Previously we described a technique for stacking cross correlations to obtain reverberation mapping results at high z. Here we present the first results from a campaign designed for this purpose. We construct stacked cross-correlation functions for the CIV and MgII lines and find a clear peak in both. We find the peak in the MgII correlation is at longer lags than CIV consistent with previous results at low redshift. For the CIV sample we are able to bin by luminosity and find evidence for increasing lags for more-luminous objects. This CIV radius-luminosity relation is consistent with previous studies but with a fraction of the observational cost. △ Less

Submitted 8 May, 2013; originally announced May 2013.

Comments: 5 pages, 5 figures, accepted for publication in MNRAS letters

SCUBA-2: The 10000 pixel bolometer camera on the James Clerk Maxwell Telescope

Authors: W. S. Holland, D. Bintley, E. L. Chapin, A. Chrysostomou, G. R. Davis, J. T. Dempsey, W. D. Duncan, M. Fich, P. Friberg, M. Halpern, K. D. Irwin, T. Jenness, B. D. Kelly, M. J. MacIntosh, E. I. Robson, D. Scott, P. A. R. Ade, E. Atad-Ettedgui, D. S. Berry, S. C. Craig, X. Gao, A. G. Gibb, G. C. Hilton, M. I. Hollister, J. B. Kycia , et al. (24 additional authors not shown)

Abstract: SCUBA-2 is an innovative 10000 pixel bolometer camera operating at submillimetre wavelengths on the James Clerk Maxwell Telescope (JCMT). The camera has the capability to carry out wide-field surveys to unprecedented depths, addressing key questions relating to the origins of galaxies, stars and planets. With two imaging arrays working simultaneously in the atmospheric windows at 450 and 850 micro… ▽ More SCUBA-2 is an innovative 10000 pixel bolometer camera operating at submillimetre wavelengths on the James Clerk Maxwell Telescope (JCMT). The camera has the capability to carry out wide-field surveys to unprecedented depths, addressing key questions relating to the origins of galaxies, stars and planets. With two imaging arrays working simultaneously in the atmospheric windows at 450 and 850 microns, the vast increase in pixel count means that SCUBA-2 maps the sky 100-150 times faster than the previous SCUBA instrument. In this paper we present an overview of the instrument, discuss the physical characteristics of the superconducting detector arrays, outline the observing modes and data acquisition, and present the early performance figures on the telescope. We also showcase the capabilities of the instrument via some early examples of the science SCUBA-2 has already undertaken. In February 2012, SCUBA-2 began a series of unique legacy surveys for the JCMT community. These surveys will take 2.5 years and the results are already providing complementary data to the shorter wavelength, shallower, larger-area surveys from Herschel. The SCUBA-2 surveys will also provide a wealth of information for further study with new facilities such as ALMA, and future telescopes such as CCAT and SPICA. △ Less

Submitted 16 January, 2013; originally announced January 2013.

Comments: 23 pages, 20 figures, 3 tables. Accepted by MNRAS

A Quasar-Galaxy Mixing Diagram: Quasar Spectral Energy Distribution Shapes in the Optical to Near-Infrared

Authors: Heng Hao, Martin Elvis, Angela Bongiorno, Gianni Zamorani, Andrea Merloni, Brandon C. Kelly, Francesca Civano, Annalisa Celotti, Luis C. Ho, Knud Jahnke, Andrea Comastri, Jonathan R. Trump, Vincenzo Mainieri, Mara Salvato, Marcella Brusa, Chris D. Impey, Anton M. Koekemoer, Giorgio Lanzuisi, Cristian Vignali, John D. Silverman, C. Megan Urry, Kevin Schawinski

Abstract: We define a quasar-galaxy mixing diagram using the slopes of their spectral energy distributions (SEDs) from 1μm to 3000Å and from 1μm to 3μm in the rest frame. The mixing diagram can easily distinguish among quasar-dominated, galaxy-dominated and reddening-dominated SED shapes. By studying the position of the 413 XMM selected Type 1 AGN in the wide-field "Cosmic Evolution Survey" (COSMOS) in the… ▽ More We define a quasar-galaxy mixing diagram using the slopes of their spectral energy distributions (SEDs) from 1μm to 3000Å and from 1μm to 3μm in the rest frame. The mixing diagram can easily distinguish among quasar-dominated, galaxy-dominated and reddening-dominated SED shapes. By studying the position of the 413 XMM selected Type 1 AGN in the wide-field "Cosmic Evolution Survey" (COSMOS) in the mixing diagram, we find that a combination of the Elvis et al. (1994, hereafter E94) quasar SED with various contributions from galaxy emission and some dust reddening is remarkably effective in describing the SED shape from 0.3-3μm for large ranges of redshift, luminosity, black hole mass and Eddington ratio of type 1 AGN. In particular, the location in the mixing diagram of the highest luminosity AGN is very close (within 1σ) to that of the E94 SED. The mixing diagram can also be used to estimate the host galaxy fraction and reddening in quasar. We also show examples of some outliers which might be AGN in different evolutionary stages compared to the majority of AGN in the quasar-host galaxy co-evolution cycle. △ Less

Submitted 3 July, 2013; v1 submitted 10 October, 2012; originally announced October 2012.

Comments: 19 pages, 13 figures, 2 tables, accepted by MNRAS

Evidence for a non-universal Kennicutt-Schmidt relationship using hierarchical Bayesian linear regression

Authors: Rahul Shetty, Brandon C. Kelly, Frank Bigiel

Abstract: For investigating the relationship between the star formation rate and gas surface density, we develop a Bayesian linear regression method that rigorously treats measurement uncertainties and accounts for hierarchical data structure. The hierarchical method simultaneously estimates the intercept, slope, and scatter about the regression line of each individual subject and the population. Using synt… ▽ More For investigating the relationship between the star formation rate and gas surface density, we develop a Bayesian linear regression method that rigorously treats measurement uncertainties and accounts for hierarchical data structure. The hierarchical method simultaneously estimates the intercept, slope, and scatter about the regression line of each individual subject and the population. Using synthetic datasets, we demonstrate that the method recovers the underlying parameters of both the individuals and the population, especially when compared to commonly employed ordinary least squares techniques, such as the bisector fit. We apply the hierarchical method to estimate the Kennicutt-Schmidt (KS) parameters of a sample of spiral galaxies compiled by Bigiel et al. (2008). We find significant variation in the KS parameters, indicating that no single relationship holds for all galaxies. This suggests that the relationship between molecular gas and star formation differs between galaxies, possibly due to the influence of other properties within a given galaxy, such as metallicity, molecular gas fraction, or stellar mass. In four galaxies the slope estimates are sub-linear (at the 95% level). We estimate the mean index of the KS relationship for the population to be 0.84, with 95% range [0.63, 1.0]. For the galaxies with sub-linear KS relationships, a possible interpretation is that CO emission is tracing some molecular gas that is not directly associated with star formation. Equivalently, a sub-linear KS relationship may be indicative of an increasing gas depletion time at higher surface densities, as traced by CO emission. The Bayesian method can account for all sources of uncertainties, including variations in the conversion of observed intensities to star formation rates and gas surface densities, and is therefore well suited for a thorough statistical analysis of the KS law. △ Less

Submitted 12 December, 2012; v1 submitted 3 October, 2012; originally announced October 2012.

Comments: 17 pages, 3 figures, 8 tables, Accepted to MNRAS - updated to match accepted version. (Abstract abridged)

Recalibration of the virial factor and M-sigma relation for local active galaxies

Authors: Daeseong Park, Brandon C. Kelly, Jong-Hak Woo, Tommaso Treu

Abstract: Determining the virial factor of the broad-line region (BLR) gas is crucial for calibrating AGN black hole mass estimators, since the measured line-of-sight velocity needs to be converted into the intrinsic virial velocity. The average virial factor has been empirically calibrated based on the M-sigma relation of quiescent galaxies, but the claimed values differ by a factor of two in recent studie… ▽ More Determining the virial factor of the broad-line region (BLR) gas is crucial for calibrating AGN black hole mass estimators, since the measured line-of-sight velocity needs to be converted into the intrinsic virial velocity. The average virial factor has been empirically calibrated based on the M-sigma relation of quiescent galaxies, but the claimed values differ by a factor of two in recent studies. We investigate the origin of the difference by measuring the M-sigma relation using an updated galaxy sample from the literature, and explore the dependence of the virial factor on various fitting methods. We find that the discrepancy is primarily caused by the sample selection, while the difference stemming from the various regression methods is marginal. However, we generally prefer the FITEXY and Bayesian estimators based on Monte Carlo simulations for the M-sigma relation. In addition, the choice of independent variable in the regression leads to ~0.2 dex variation in the virial factor inferred from the calibration process. Based on the determined virial factor, we present the updated M-sigma relation of local active galaxies. △ Less

Submitted 17 September, 2012; originally announced September 2012.

Comments: Accepted for publication in ApJS. 13 pages, 9 figures

Disclosing the Radio Loudness Distribution Dichotomy in Quasars: An Unbiased Monte Carlo Approach Applied to the SDSS-FIRST Quasar Sample

Authors: Mislav Balokovic, Vernesa Smolcic, Zeljko Ivezic, Gianni Zamorani, Eva Schinnerer, Brandon C. Kelly

Abstract: We investigate the dichotomy in the radio loudness distribution of quasars by modelling their radio emission and various selection effects using a Monte Carlo approach. The existence of two physically distinct quasar populations, the radio-loud and radio-quiet quasars, is controversial and over the last decade a bimodal distribution of radio loudness of quasars has been both affirmed and disputed.… ▽ More We investigate the dichotomy in the radio loudness distribution of quasars by modelling their radio emission and various selection effects using a Monte Carlo approach. The existence of two physically distinct quasar populations, the radio-loud and radio-quiet quasars, is controversial and over the last decade a bimodal distribution of radio loudness of quasars has been both affirmed and disputed. We model the quasar radio luminosity distribution with simple unimodal and bimodal distribution functions. The resulting simulated samples are compared to a fiducial sample of 8,300 quasars drawn from the SDSS DR7 Quasar Catalog and combined with radio observations from the FIRST survey. Our results indicate that the SDSS-FIRST sample is best described by a radio loudness distribution which consists of two components, with 12+/-1 % of sources in the radio-loud component. On the other hand, the evidence for a local minimum in the loudness distribution (bimodality) is not strong and we find that previous claims for its existence were probably affected by the incompleteness of the FIRST survey close to its faint limit. We also investigate the redshift and luminosity dependence of the radio loudness distribution and find tentative evidence that at high redshift radio-loud quasars were rarer, on average "louder", and exhibited a smaller range in radio loudness. In agreement with other recent work, we conclude that the SDSS-FIRST sample strongly suggests that the radio loudness distribution of quasars is not a universal function, and that more complex models than presented here are needed to fully explain available observations. △ Less

Submitted 5 September, 2012; originally announced September 2012.

Comments: Accepted for publication in ApJ; 13 pages, 10 figures

Journal ref: ApJ, 759, 30 (2012)

The Demographics of Broad Line Quasars in the Mass-Luminosity Plane II. Black Hole Mass and Eddington Ratio Functions

Authors: Brandon C. Kelly, Yue Shen

Abstract: We employ a flexible Bayesian technique to estimate the black hole mass and Eddington ratio functions for Type 1 (i.e., broad line) quasars from a uniformly-selected data set of ~58,000 quasars from the SDSS DR7. We find that the SDSS becomes significantly incomplete at M_{BH} < 3 x 10^8 M_{Sun} or L / L_{Edd} < 0.07, and that the number densities of Type 1 quasars continue to increase down to the… ▽ More We employ a flexible Bayesian technique to estimate the black hole mass and Eddington ratio functions for Type 1 (i.e., broad line) quasars from a uniformly-selected data set of ~58,000 quasars from the SDSS DR7. We find that the SDSS becomes significantly incomplete at M_{BH} < 3 x 10^8 M_{Sun} or L / L_{Edd} < 0.07, and that the number densities of Type 1 quasars continue to increase down to these limits. Both the mass and Eddington ratio functions show evidence of downsizing, with the most massive and highest Eddington ratio black holes experiencing Type 1 quasar phases first, although the Eddington ratio number densities are flat at z < 2. We estimate the maximum Eddington ratio of Type 1 quasars in the observable Universe to be L / L_{Edd} ~ 3. Consistent with our results in Paper I, we do not find statistical evidence for a so-called "sub-Eddington boundary" in the mass-luminosity plane of broad line quasars, and demonstrate that such an apparent boundary in the observed distribution can be caused by selection effect and errors in virial BH mass estimates. Based on the typical Eddington ratio in a given mass bin, we estimate typical growth times for the black holes in Type 1 quasars and find that they are typically comparable to or longer than the age of the universe, implying an earlier phase of accelerated (i.e., with higher Eddington ratios) and possibly obscured growth. The large masses probed by our sample imply that most of our black holes reside in what are locally early type galaxies, and we interpret our results within the context of models of self-regulated black hole growth. △ Less

Submitted 18 December, 2012; v1 submitted 3 September, 2012; originally announced September 2012.

Comments: Submitted to ApJ, 25 pages (emulateapj), 15 figures; revised to match accepted version with primary changes to the introduction and discussion, replaced Fig 12

The linewidth-size relationship in the dense ISM of the Central Molecular Zone

Authors: Rahul Shetty, Christopher N. Beaumont, Michael G. Burton, Brandon C. Kelly, Ralf S. Klessen

Abstract: The linewidth (sigma) - size (R) relationship has been extensively measured and analysed, in both the local ISM and in nearby normal galaxies. Generally, a power-law describes the relationship well with an index ranging from 0.2-0.6, now referred to as one of "Larson's Relationships." The nature of turbulence and star formation is considered to be intimately related to these relationships, so eval… ▽ More The linewidth (sigma) - size (R) relationship has been extensively measured and analysed, in both the local ISM and in nearby normal galaxies. Generally, a power-law describes the relationship well with an index ranging from 0.2-0.6, now referred to as one of "Larson's Relationships." The nature of turbulence and star formation is considered to be intimately related to these relationships, so evaluating the sigma-R correlations in various environments is important for developing a comprehensive understanding of the ISM. We measure the sigma-R relationship in the Central Molecular Zone (CMZ) of the Galactic Centre using spectral line observations of the high density tracers N2H+, HCN, H13CN, and HCO+. We use dendrograms, which map the hierarchical nature of the position-position-velocity (PPV) data, to compute sigma and R of contiguous structures. The dispersions range from ~2-30 km/s in structures spanning sizes 2-40 pc, respectively. By performing Bayesian inference, we show that a power-law with exponent 0.3-1.1 can reasonably describe the sigma-R trend. We demonstrate that the derived sigma-R relationship is independent of the locations in the PPV dataset where sigma and R are measured. The uniformity in the sigma-R relationship suggests turbulence in the CMZ is driven on the large scales beyond >30 pc. We compare the CMZ sigma-R relationship to that measured in the Galactic molecular cloud Perseus. The exponents between the two systems are similar, suggestive of a connection between the turbulent properties within a cloud to its ambient medium. Yet, the velocity dispersion in the CMZ is systematically higher, resulting in a coefficient that is nearly five times larger. The systematic enhancement of turbulent velocities may be due to the combined effects of increased star formation activity, larger densities, and higher pressures relative to the local ISM. △ Less

Submitted 25 June, 2012; originally announced June 2012.

Comments: 11 pages, 8 figures, Accepted for publication in MNRAS

Composite reverberation mapping

Authors: S. Fine, T. Shanks, S. M. Croom, P. Green, B. C. Kelly, E. Berge, R. Chornock, W. S. Burgett, E. A. Magnier, P. A. Price

Abstract: Reverberation mapping offers one of the best techniques for studying the inner regions of QSOs. It is based on cross-correlating continuum and emission-line light curves. New time-resolved optical surveys will produce well sampled light curves for many thousands of QSOs. We explore the potential of stacking samples to produce composite cross-correlations for groups of objects that have well sample… ▽ More Reverberation mapping offers one of the best techniques for studying the inner regions of QSOs. It is based on cross-correlating continuum and emission-line light curves. New time-resolved optical surveys will produce well sampled light curves for many thousands of QSOs. We explore the potential of stacking samples to produce composite cross-correlations for groups of objects that have well sampled continuum light curves, but only a few (~2) emission-line measurements. This technique exploits current and future wide-field optical monitoring surveys (e.g. Pan-STARRS, LSST) and the multiplexing capability of multi-object spectrographs (e.g. 2dF, Hectospec) to significantly reduce the observational expense of reverberation mapping, in particular at high redshift (0.5 to 2.5). We demonstrate the technique using simulated QSO light curves and explore the biases involved when stacking cross-correlations in some simplified situations. We show that stacked cross correlations have smaller amplitude peaks compared to well sampled correlation functions as the mean flux of the emission light curve is poorly constrained. However, the position of the peak remains intact. We find there can be `kinks' in stacked correlation functions due to different measurements contributing to different parts of the correlation function. Using the Pan-STARRS Medium-Deep Survey (MDS) as a template we show that cross-correlation lags should be measurable in a sample size of 500 QSOs that have weekly photometric monitoring and two spectroscopic observations. Finally we apply the technique to a small sample (42) of QSOs that have light curves from the MDS. We find no indication of a peak in the stacked cross-correlation. A larger spectroscopic sample is required to produce robust reverberation lags. △ Less

Submitted 7 May, 2012; originally announced May 2012.

Comments: 11 pages, 11 figures, accepted for publication in MNRAS

Dust SEDs in the era of Herschel and Planck: a Hierarchical Bayesian fitting technique

Authors: Brandon C. Kelly, Rahul Shetty, Amelia M. Stutz, Jens Kauffmann, Alyssa A. Goodman, Ralf Launhardt

Abstract: We present a hierarchical Bayesian method for fitting infrared spectral energy distributions (SEDs) of dust emission to observed fluxes. Under the standard assumption of optically thin single temperature (T) sources the dust SED as represented by a power--law modified black body is subject to a strong degeneracy between T and the spectral index beta. The traditional non-hierarchical approaches, ty… ▽ More We present a hierarchical Bayesian method for fitting infrared spectral energy distributions (SEDs) of dust emission to observed fluxes. Under the standard assumption of optically thin single temperature (T) sources the dust SED as represented by a power--law modified black body is subject to a strong degeneracy between T and the spectral index beta. The traditional non-hierarchical approaches, typically based on chi-square minimization, are severely limited by this degeneracy, as it produces an artificial anti-correlation between T and beta even with modest levels of observational noise. The hierarchical Bayesian method rigorously and self-consistently treats measurement uncertainties, including calibration and noise, resulting in more precise SED fits. As a result, the Bayesian fits do not produce any spurious anti-correlations between the SED parameters due to measurement uncertainty. We demonstrate that the Bayesian method is substantially more accurate than the chi-square fit in recovering the SED parameters, as well as the correlations between them. As an illustration, we apply our method to Herschel and sub millimeter ground-based observations of the star-forming Bok globule CB244. This source is a small, nearby molecular cloud containing a single low-mass protostar and a starless core. We find that T and beta are weakly positively correlated -- in contradiction with the chi-square fits, which indicate a T-beta anti-correlation from the same data-set. Additionally, in comparison to the chi-square fits the Bayesian SED parameter estimates exhibit a reduced range in values. △ Less

Submitted 17 April, 2012; v1 submitted 29 February, 2012; originally announced March 2012.

Comments: 20 pages, 9 figures, ApJ format, revised version matches ApJ-accepted version

Measurement Error Models in Astronomy

Authors: Brandon C. Kelly

Abstract: I discuss the effects of measurement error on regression and density estimation. I review the statistical methods that have been developed to correct for measurement error that are most popular in astronomical data analysis, discussing their advantages and disadvantages. I describe functional models for accounting for measurement error in regression, with emphasis on the methods of moments approac… ▽ More I discuss the effects of measurement error on regression and density estimation. I review the statistical methods that have been developed to correct for measurement error that are most popular in astronomical data analysis, discussing their advantages and disadvantages. I describe functional models for accounting for measurement error in regression, with emphasis on the methods of moments approach and the modified loss function approach. I then describe structural models for accounting for measurement error in regression and density estimation, with emphasis on maximum-likelihood and Bayesian methods. As an example of a Bayesian application, I analyze an astronomical data set subject to large measurement errors and a non-linear dependence between the response and covariate. I conclude with some directions for future research. △ Less

Submitted 7 December, 2011; originally announced December 2011.

Comments: 15 pages, 2 figures, review based on invited talk at "Statistical Challenges in Modern Astronomy V", held at Penn State in June 2011, to appear in the proceedings

Mass Functions of Supermassive Black Holes Across Cosmic Time

Authors: Brandon C. Kelly, Andrea Merloni

Abstract: The black hole mass function of supermassive black holes describes the evolution of the distribution of black hole mass. It is one of the primary empirical tools available for mapping the growth of supermassive black holes and for constraining theoretical models of their evolution. In this review we discuss methods for estimating the black hole mass function, including their advantages and disadva… ▽ More The black hole mass function of supermassive black holes describes the evolution of the distribution of black hole mass. It is one of the primary empirical tools available for mapping the growth of supermassive black holes and for constraining theoretical models of their evolution. In this review we discuss methods for estimating the black hole mass function, including their advantages and disadvantages. We also review the results of using these methods for estimating the mass function of both active and inactive black holes. In addition, we review current theoretical models for the growth of supermassive black holes that predict the black hole mass function. We conclude with a discussion of directions for future research which will lead to improvement in both empirical and theoretical determinations of the mass function of supermassive black holes. △ Less

Submitted 6 December, 2011; originally announced December 2011.

Comments: 40 pages, 7 figures, review paper accepted for the Advances in Astronomy Special Issue "Seeking for the Leading Actor on the Cosmic Stage: Galaxies versus Supermassive Black Holes"

A Description of Quasar Variability Measured Using Repeated SDSS and POSS Imaging

Authors: Chelsea L. MacLeod, Zeljko Ivezic, Branimir Sesar, Wim de Vries, Christopher S. Kochanek, Brandon C. Kelly, Andrew C. Becker, Robert H. Lupton, Patrick B. Hall, Gordon T. Richards, Scott F. Anderson, Donald P. Schneider

Abstract: We provide a quantitative description and statistical interpretation of the optical continuum variability of quasars. The Sloan Digital Sky Survey (SDSS) has obtained repeated imaging in five UV-to-IR photometric bands for 33,881 spectroscopically confirmed quasars. About 10,000 quasars have an average of 60 observations in each band obtained over a decade along Stripe 82 (S82), whereas the remain… ▽ More We provide a quantitative description and statistical interpretation of the optical continuum variability of quasars. The Sloan Digital Sky Survey (SDSS) has obtained repeated imaging in five UV-to-IR photometric bands for 33,881 spectroscopically confirmed quasars. About 10,000 quasars have an average of 60 observations in each band obtained over a decade along Stripe 82 (S82), whereas the remaining ~25,000 have 2-3 observations due to scan overlaps. The observed time lags span the range from a day to almost 10 years, and constrain quasar variability at rest-frame time lags of up to 4 years, and at rest-frame wavelengths from 1000A to 6000A. We publicly release a user-friendly catalog of quasars from the SDSS Data Release 7 that have been observed at least twice in SDSS or once in both SDSS and the Palomar Observatory Sky Survey, and we use it to analyze the ensemble properties of quasar variability. Based on a damped random walk (DRW) model defined by a characteristic time scale and an asymptotic variability amplitude that scale with the luminosity, black hole mass, and rest wavelength for individual quasars calibrated in S82, we can fully explain the ensemble variability statistics of the non-S82 quasars such as the exponential distribution of large magnitude changes. All available data are consistent with the DRW model as a viable description of the optical continuum variability of quasars on time scales of ~5-2000 days in the rest frame. We use these models to predict the incidence of quasar contamination in transient surveys such as those from PTF and LSST. △ Less

Submitted 4 May, 2012; v1 submitted 3 December, 2011; originally announced December 2011.

Comments: 33 pages, 19 figures, replaced with accepted version. Catalog is available at http://www.astro.washington.edu/users/ivezic/macleod/qso_dr7/

Journal ref: 2012, ApJ, 753, 106