Showing 1–50 of 130 results for author: Haggard, D

Multistructured accretion flow of Sgr A* II: Signatures of a Cool Accretion Disk in Hydrodynamic Simulations of Stellar Winds

Authors: Mayura Balakrishnan, Christopher M. P. Russell, Lia Corrales, Diego Calderón, Jorge Cuadra, Daryl Haggard, Sera Markoff, Joey Neilsen, Michael Nowak, Q. Daniel Wang, Fred Baganoff

Abstract: Hydrodynamic simulations of the stellar winds from Wolf-Rayet stars within the Galactic Center can provide predictions for the X-ray spectrum of supermassive black hole Sgr A*. Herein, we present results from updated smooth particle hydrodynamics simulations, building on the architecture of Cuadra et al. (2015); Russell et al. (2017), finding that a cold gas disk forms around Sgr A* with a simulat… ▽ More Hydrodynamic simulations of the stellar winds from Wolf-Rayet stars within the Galactic Center can provide predictions for the X-ray spectrum of supermassive black hole Sgr A*. Herein, we present results from updated smooth particle hydrodynamics simulations, building on the architecture of Cuadra et al. (2015); Russell et al. (2017), finding that a cold gas disk forms around Sgr A* with a simulation runtime of 3500 years. This result is consistent with previous grid-based simulations, demonstrating that a cold disk can form regardless of numerical method. We examine the plasma scenarios arising from an environment with and without this cold disk, by generating synthetic spectra for comparison to the quiescent Fe K alpha Sgr A* spectrum from Chandra HETG-S, taken through the Chandra X-ray Visionary Program. We find that current and future X-ray missions are unlikely to distinguish between the kinematic signatures in the plasma in these two scenarios. Nonetheless, the stellar wind plasma model presents a good fit to the dispersed Chandra spectra within 1.5" of Sgr A*. We compare our results to the Radiatively Inefficient Accretion Flow (RIAF) model fit to the HETG-S spectrum presented in Paper I and find that the Bayesian model evidence does not strongly favor either model. With 9" angular resolution and high spectral resolution of the X-IFU, NewAthena will offer a clearer differentiation between the RIAF plasma model and hydrodynamic simulations, but only a future X-ray mission with arcsecond resolution will significantly advance our understanding of Sgr A*'s accretion flow in X-rays. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 14 pages, 6 figures. Accepted for publication in ApJ

Multistructured accretion flow of Sgr A* I: Examination of a RIAF model

Authors: Mayura Balakrishnan, Lia Corrales, Sera Markoff, Michael Nowak, Daryl Haggard, Q. Daniel Wang, Joey Neilsen, Christopher M. P. Russell, Diego Calderón, Jorge Cuadra, Fred Baganoff

Abstract: The extreme low-luminosity supermassive black hole Sgr A* provides a unique laboratory in which to test radiatively inefficient accretion flow (RIAF) models. Previous fits to the quiescent Chandra ACIS-S spectrum found a RIAF model with an equal inflow-outflow balance works well. In this work, we apply the RIAF model to the Chandra HETG-S spectrum obtained through the Chandra X-ray Visionary Progr… ▽ More The extreme low-luminosity supermassive black hole Sgr A* provides a unique laboratory in which to test radiatively inefficient accretion flow (RIAF) models. Previous fits to the quiescent Chandra ACIS-S spectrum found a RIAF model with an equal inflow-outflow balance works well. In this work, we apply the RIAF model to the Chandra HETG-S spectrum obtained through the Chandra X-ray Visionary Program, which displays features suggestive of temperature and velocity structures within the plasma. A comprehensive forward model analysis accounting for the accretion flow geometry and HETG-S instrumental effects is required for a full interpretation of the quiescent Chandra HETG-S spectrum. We present a RIAF model that takes these effects into account. Our fits to the high-resolution gratings spectrum indicate an inflow balanced by an outflow ($s \sim 1$) alongside a temperature profile that appears shallower than what would be expected from a gravitational potential following $1/r$. The data require that the abundance of Iron relative to solar is $Z_{Fe} < 0.32 Z_\odot$ (90\% credible interval), much lower than the $2~Z_\odot$ metallicity measured in nearby late-type giants. While future missions like NewAthena will provide higher spectral resolution, source separation will continue to be a problem. Leveraging Chandra's unparalleled spatial resolution, which is not expected to be surpassed for decades, remains essential for detailed investigations of the densely populated Galactic Center in X-rays. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 11 pages, 5 figures, 1 table. Accepted for publication in ApJ

WISDOM project XX -- Strong shear tearing molecular clouds apart in NGC 524

Authors: Anan Lu, Daryl Haggard, Martin Bureau, Jindra Gensior, Sarah Jeffreson, Carmelle Robert, Thomas G. Williams, Fu-Heng Liang, Woorak Choi, Timothy A. Davis, Sara Babic, Hope Boyce, Benjamin Cheung, Laurent Drissen, Jacob S. Elford, Lijie Liu, Thomas Martin, Carter Rhea, Laurie Rousseau-Nepton, Ilaria Ruffa

Abstract: Early-type galaxies (ETGs) are known to harbour dense spheroids of stars but scarce star formation (SF). Approximately a quarter of these galaxies have rich molecular gas reservoirs yet do not form stars efficiently. We study here the ETG NGC~524, with strong shear suspected to result in a smooth molecular gas disc and low star-formation efficiency (SFE). We present new spatially-resolved observat… ▽ More Early-type galaxies (ETGs) are known to harbour dense spheroids of stars but scarce star formation (SF). Approximately a quarter of these galaxies have rich molecular gas reservoirs yet do not form stars efficiently. We study here the ETG NGC~524, with strong shear suspected to result in a smooth molecular gas disc and low star-formation efficiency (SFE). We present new spatially-resolved observations of the \textsuperscript{12}CO(2-1)-emitting cold molecular gas from the Atacama Large Millimeter/sub-millimeter Array (ALMA) and of the warm ionised-gas emission lines from SITELLE at the Canada-France-Hawaii Telescope. Although constrained by the resolution of the ALMA observations ($\approx37$~pc), we identify only $52$ GMCs with radii ranging from $30$ to $140$~pc, a low mean molecular gas mass surface density $\langleΣ_{\rm gas}\rangle\approx125$~M$_\odot$~pc$^{-2}$ and a high mean virial parameter $\langleα_{\rm obs,vir}\rangle\approx5.3$. We measure spatially-resolved molecular gas depletion times ($τ_{\rm dep}\equiv1/{\rm SFE}$) with a spatial resolution of $\approx100$~pc within a galactocentric distance of $1.5$~kpc. The global depletion time is $\approx2.0$~Gyr but $τ_{\rm dep}$ increases toward the galaxy centre, with a maximum $τ_{\rm dep,max}\approx5.2$~Gyr. However, no pure \ion{H}{II} region is identified in NGC~524 using ionised-gas emission-line ratio diagnostics, so the $τ_{\rm dep}$ inferred are in fact lower limits. Measuring the GMC properties and dynamical states, we conclude that shear is the dominant mechanism shaping the molecular gas properties and regulating SF in NGC~524. This is supported by analogous analyses of the GMCs in a simulated ETG similar to NGC~524. △ Less

Submitted 3 June, 2024; originally announced June 2024.

Comments: 17 pages, 10 figures. To be published in MNRAS, accepted on May 27

Simulation-based inference of radio millisecond pulsars in globular clusters

Authors: Joanna Berteaud, Christopher Eckner, Francesca Calore, Maïca Clavel, Daryl Haggard

Abstract: Millisecond pulsars (MSPs) are abundant in globular clusters (GCs), which offer favorable environments for their creation. While the advent of recent, powerful facilities led to a rapid increase in MSP discoveries in GCs through pulsation searches, detection biases persist. In this work, we investigate the ability of current and future detections in GCs to constrain the parameters of the MSP popul… ▽ More Millisecond pulsars (MSPs) are abundant in globular clusters (GCs), which offer favorable environments for their creation. While the advent of recent, powerful facilities led to a rapid increase in MSP discoveries in GCs through pulsation searches, detection biases persist. In this work, we investigate the ability of current and future detections in GCs to constrain the parameters of the MSP population in GCs through a careful study of their luminosity function. Parameters of interest are the number of MSPs hosted by a GC, as well as the mean and the width of their luminosity function, which are typically affected by large uncertainties. While, as we show, likelihood-based studies can lead to ill-behaved posterior on the size of the MSP population, we introduce a novel, likelihood-free analysis, based on Marginal Neural Ratio Estimation, which consistently produces well-behaved posteriors. We focus on the GC Terzan 5, which currently counts 48 detected MSPs. We find that about 158 MSPs should be hosted in this GC, but the uncertainty on this number remains large. We explore the performance of our new method on simulated Terzan 5-like datasets mimicking possible future observational outcomes. We find that significant improvement on the posteriors can be obtained by adding a reliable measurement of the diffuse radio emission of the GC to the analysis or by improving the detection threshold of current radio pulsation surveys by at least a factor two. △ Less

Submitted 24 May, 2024; originally announced May 2024.

Comments: 24 pages, 16 figures, 5 tables. Comments are welcome

Report number: LAPTH-028/24

Chandra Study of the Proper Motion of HST-1 in the Jet of M87

Authors: Rameshan Thimmappa, Joey Neilsen, Daryl Haggard, Michael A. Nowak, Sera Markoff

Abstract: The radio galaxy M87 is well known for its jet, which features a series of bright knots observable from radio to X-ray wavelengths. We analyze the X-ray image and flux variability of the knot HST-1 in the jet. Our analysis includes all 112 available Chandra ACIS-S observations from 2000-2021, with a total exposure time of $\sim$887 ks. We use de-convolved images to study the brightness profile of… ▽ More The radio galaxy M87 is well known for its jet, which features a series of bright knots observable from radio to X-ray wavelengths. We analyze the X-ray image and flux variability of the knot HST-1 in the jet. Our analysis includes all 112 available Chandra ACIS-S observations from 2000-2021, with a total exposure time of $\sim$887 ks. We use de-convolved images to study the brightness profile of the X-ray jet and measure the relative separation between the core and HST-1. From 2003-2005 (which coincides with a bright flare from HST-1), we find a correlation between the flux of HST-1 and its offset from the core. In subsequent data, we find a steady increase in this offset, which implies a bulk superluminal motion for HST-1 of 6.6$\pm$0.9 c (2.0$\pm$0.3 pc yr$^{-1}$), in keeping with prior results. We discuss models for the flux-offset correlation that feature either two or four emission regions separated by tens of parsecs. We attribute these results to moving shocks in the jet, that allow us to measure the internal structure of the jet. △ Less

Submitted 1 May, 2024; v1 submitted 30 April, 2024; originally announced April 2024.

Comments: Accepted for publication in ApJ. 9 pages, 13 figures, 1 table

Radio Scrutiny of the X-ray-Weak Tail of Low-Mass Active Galactic Nuclei: A Novel Signature of High-Eddington Accretion?

Authors: Jeremiah D. Paul, Richard M. Plotkin, W. N. Brandt, Christopher H. Ellis, Elena Gallo, Jenny E. Greene, Luis C. Ho, Amy E. Kimball, Daryl Haggard

Abstract: The supermassive black holes ($M_{\rm BH} \sim 10^{6}$$-$$10^{10}~M_\odot$) that power luminous active galactic nuclei (AGNs), i.e., quasars, generally show a correlation between thermal disk emission in the ultraviolet (UV) and coronal emission in hard X-rays. In contrast, some "massive" black holes (mBHs; $M_{\rm BH} \sim 10^{5}$$-$$10^{6}~M_\odot$) in low-mass galaxies present curious X-ray pro… ▽ More The supermassive black holes ($M_{\rm BH} \sim 10^{6}$$-$$10^{10}~M_\odot$) that power luminous active galactic nuclei (AGNs), i.e., quasars, generally show a correlation between thermal disk emission in the ultraviolet (UV) and coronal emission in hard X-rays. In contrast, some "massive" black holes (mBHs; $M_{\rm BH} \sim 10^{5}$$-$$10^{6}~M_\odot$) in low-mass galaxies present curious X-ray properties with coronal radiative output up to 100$\times$ weaker than expected. To examine this issue, we present new and archival 10 GHz Very Large Array observations of a sample of high-accretion-rate (Eddington ratios $L_{\rm bol}/L_{\rm Edd} > 0.1$), mBH-powered AGNs with Chandra X-ray coverage. Empirical correlations previously revealed in samples of radio-quiet, high-Eddington AGNs indicate that the radio$-$X-ray luminosity ratio, $L_{\rm R}/L_{\rm X}$, is approximately constant. Through multiwavelength analysis, we instead find that the X-ray-weaker mBHs in our sample tend toward larger values of $L_{\rm R}/L_{\rm X}$ even though they remain radio-quiet per their optical$-$UV properties. This trend results in a tentative but highly intriguing correlation between $L_{\rm R}/L_{\rm X}$ and X-ray weakness, which we argue is consistent with a scenario in which X-rays may be preferentially obscured from our line of sight by a "slim" accretion disk. We compare this observation to weak emission-line quasars (AGNs with exceptionally weak broad-line emission and a significant X-ray-weak fraction) and conclude by suggesting that our results may offer a new observational signature for finding high-accretion-rate AGNs. △ Less

Submitted 2 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ. 26 pages, 8 figures, 5 tables

An Extreme Ultra-Compact X-ray Binary in a Globular Cluster: Multi-Wavelength Observations of RZ2109 Explored in a Triple System Framework

Authors: Kristen C. Dage, Arash Bahramian, Smadar Naoz, Alexey Bobrick, Wasundara Athukoralalage, McKinley C. Brumback, Daryl Haggard, Arunav Kundu, Stephen E. Zepf

Abstract: The globular cluster ultraluminous X-ray source, RZ2109, is a complex and unique system which has been detected at X-ray, ultra-violet, and optical wavelengths. Based on almost 20 years of Chandra and XMM-Newton observations, the X-ray luminosity exhibits order-of-magnitude variability, with the peak flux lasting on the order of a few hours. We perform robust time series analysis on the archival X… ▽ More The globular cluster ultraluminous X-ray source, RZ2109, is a complex and unique system which has been detected at X-ray, ultra-violet, and optical wavelengths. Based on almost 20 years of Chandra and XMM-Newton observations, the X-ray luminosity exhibits order-of-magnitude variability, with the peak flux lasting on the order of a few hours. We perform robust time series analysis on the archival X-ray observations and find that this variability is periodic on a timescale of 1.3 $\pm 0.04$ days. The source also demonstrates broad [OIII] 5007 Angstrom emission, which has been observed since 2004, suggesting a white dwarf donor and therefore an ultra-compact X-ray binary. We present new spectra from 2020 and 2022, marking eighteen years of observed [OIII] emission from this source. Meanwhile, we find that the globular cluster counterpart is unusually bright in the NUV/UVW2 band. Finally, we discuss RZ2109 in the context of the eccentric Kozai Lidov mechanism and show that the observed 1.3 day periodicity can be used to place constraints on the tertiary configuration, ranging from 20 minutes (for a 0.1 ${\rm M}_\odot$ companion) to approximately 95 minutes (for a 1 ${\rm M}_\odot$ companion), if the eccentric Kozai Lidov mechanism is at the origin of the periodic variability. △ Less

Submitted 14 March, 2024; originally announced March 2024.

Comments: Accepted to MNRAS

Journal ref: 2024MNRAS.tmp..625D

Ordered magnetic fields around the 3C 84 central black hole

Authors: G. F. Paraschos, J. -Y. Kim, M. Wielgus, J. Röder, T. P. Krichbaum, E. Ros, I. Agudo, I. Myserlis, M. Moscibrodzka, E. Traianou, J. A. Zensus, L. Blackburn, C. -K. Chan, S. Issaoun, M. Janssen, M. D. Johnson, V. L. Fish, K. Akiyama, A. Alberdi, W. Alef, J. C. Algaba, R. Anantua, K. Asada, R. Azulay, U. Bach , et al. (258 additional authors not shown)

Abstract: 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures a… ▽ More 3C84 is a nearby radio source with a complex total intensity structure, showing linear polarisation and spectral patterns. A detailed investigation of the central engine region necessitates the use of VLBI above the hitherto available maximum frequency of 86GHz. Using ultrahigh resolution VLBI observations at the highest available frequency of 228GHz, we aim to directly detect compact structures and understand the physical conditions in the compact region of 3C84. We used EHT 228GHz observations and, given the limited (u,v)-coverage, applied geometric model fitting to the data. We also employed quasi-simultaneously observed, multi-frequency VLBI data for the source in order to carry out a comprehensive analysis of the core structure. We report the detection of a highly ordered, strong magnetic field around the central, SMBH of 3C84. The brightness temperature analysis suggests that the system is in equipartition. We determined a turnover frequency of $ν_m=(113\pm4)$GHz, a corresponding synchrotron self-absorbed magnetic field of $B_{SSA}=(2.9\pm1.6)$G, and an equipartition magnetic field of $B_{eq}=(5.2\pm0.6)$G. Three components are resolved with the highest fractional polarisation detected for this object ($m_\textrm{net}=(17.0\pm3.9)$%). The positions of the components are compatible with those seen in low-frequency VLBI observations since 2017-2018. We report a steeply negative slope of the spectrum at 228GHz. We used these findings to test models of jet formation, propagation, and Faraday rotation in 3C84. The findings of our investigation into different flow geometries and black hole spins support an advection-dominated accretion flow in a magnetically arrested state around a rapidly rotating supermassive black hole as a model of the jet-launching system in the core of 3C84. However, systematic uncertainties due to the limited (u,v)-coverage, however, cannot be ignored. △ Less

Submitted 1 February, 2024; originally announced February 2024.

Comments: 15 pages, 6 figures, published in A&A

Journal ref: Issue: A&A Volume 682, February 2024; Article number: L3; Number of pages: 15

Chasing the break: Tracing the full evolution of a black hole X-ray binary jet with multi-wavelength spectral modeling

Authors: Constanza Echiburú-Trujillo, Alexandra J. Tetarenko, Daryl Haggard, Thomas D. Russell, Karri I. I. Koljonen, Arash Bahramian, Jingyi Wang, Michael Bremer, Joe Bright, Piergiorgio Casella, David M. Russell, Diego Altamirano, M. Cristina Baglio, Tomaso Belloni, Chiara Ceccobello, Stephane Corbel, Maria Diaz Trigo, Dipankar Maitra, Aldrin Gabuya, Elena Gallo, Sebastian Heinz, Jeroen Homan, Erin Kara, Elmar Körding, Fraser Lewis , et al. (13 additional authors not shown)

Abstract: Black hole X-ray binaries (BH XRBs) are ideal targets to study the connection between accretion inflow and jet outflow. Here we present quasi-simultaneous, multi-wavelength observations of the Galactic black hole system MAXI J1820+070, throughout its 2018-2019 outburst. Our data set includes coverage from the radio through X-ray bands from 17 different instruments/telescopes, and encompasses 19 ep… ▽ More Black hole X-ray binaries (BH XRBs) are ideal targets to study the connection between accretion inflow and jet outflow. Here we present quasi-simultaneous, multi-wavelength observations of the Galactic black hole system MAXI J1820+070, throughout its 2018-2019 outburst. Our data set includes coverage from the radio through X-ray bands from 17 different instruments/telescopes, and encompasses 19 epochs over a 7 month time period, resulting in one of the most well-sampled multi-wavelength data sets of a BH XRB outburst to date. With our data, we compile and model the broad-band spectra of this source using a phenomenological model that includes emission from the jet, companion star, and accretion flow. This modeling allows us to track the evolution of the spectral break in the jet spectrum, a key observable that samples the jet launching region. We find that the spectral break location changes over at least $\approx3$ orders of magnitude in electromagnetic frequency over this period. Using these spectral break measurements, we link the full cycle of jet behavior, including the rising, quenching, and re-ignition, to the changing accretion flow properties as the source evolves through its different accretion states. Our analyses show a consistent jet behavior with other sources in similar phases of their outbursts, reinforcing that the jet quenching and recovery may be a global feature of BH XRB systems in outburst. Our results also provide valuable evidence supporting a close connection between the geometry of the inner accretion flow and the base of the jet. △ Less

Submitted 30 January, 2024; v1 submitted 19 November, 2023; originally announced November 2023.

Comments: Accepted for publication in ApJ

Journal ref: 2024 ApJ 962 116

Prospects for Time-Domain and Multi-Messenger Science with AXIS

Authors: The AXIS Time-Domain, Multi-Messenger Science Working Group, :, Riccardo Arcodia, Franz E. Bauer, S. Bradley Cenko, Kristen C. Dage, Daryl Haggard, Wynn C. G. Ho, Erin Kara, Michael Koss, Tingting Liu, Labani Mallick, Michela Negro, Pragati Pradhan, J. Quirola-Vasquez, Mark T. Reynolds, Claudio Ricci, Richard E. Rothschild, Navin Sridhar, Eleonora Troja, Yuhan Yao

Abstract: The Advanced X-ray Imaging Satellite (AXIS) promises revolutionary science in the X-ray and multi-messenger time domain. AXIS will leverage excellent spatial resolution (<1.5 arcsec), sensitivity (80x that of Swift), and a large collecting area (5-10x that of Chandra) across a 24-arcmin diameter field of view to discover and characterize a wide range of X-ray transients from supernova-shock breako… ▽ More The Advanced X-ray Imaging Satellite (AXIS) promises revolutionary science in the X-ray and multi-messenger time domain. AXIS will leverage excellent spatial resolution (<1.5 arcsec), sensitivity (80x that of Swift), and a large collecting area (5-10x that of Chandra) across a 24-arcmin diameter field of view to discover and characterize a wide range of X-ray transients from supernova-shock breakouts to tidal disruption events to highly variable supermassive black holes. The observatory's ability to localize and monitor faint X-ray sources opens up new opportunities to hunt for counterparts to distant binary neutron star mergers, fast radio bursts, and exotic phenomena like fast X-ray transients. AXIS will offer a response time of <2 hours to community alerts, enabling studies of gravitational wave sources, high-energy neutrino emitters, X-ray binaries, magnetars, and other targets of opportunity. This white paper highlights some of the discovery science that will be driven by AXIS in this burgeoning field of time domain and multi-messenger astrophysics. △ Less

Submitted 13 November, 2023; originally announced November 2023.

Comments: This White Paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS White Papers can be found at http://axis.astro.umd.edu

Overview of the Advanced X-ray Imaging Satellite (AXIS)

Authors: Christopher S. Reynolds, Erin A. Kara, Richard F. Mushotzky, Andrew Ptak, Michael J. Koss, Brian J. Williams, Steven W. Allen, Franz E. Bauer, Marshall Bautz, Arash Bodaghee, Kevin B. Burdge, Nico Cappelluti, Brad Cenko, George Chartas, Kai-Wing Chan, Lía Corrales, Tansu Daylan, Abraham D. Falcone, Adi Foord, Catherine E. Grant, Mélanie Habouzit, Daryl Haggard, Sven Herrmann, Edmund Hodges-Kluck, Oleg Kargaltsev , et al. (18 additional authors not shown)

Abstract: The Advanced X-ray Imaging Satellite (AXIS) is a Probe-class concept that will build on the legacy of the Chandra X-ray Observatory by providing low-background, arcsecond-resolution imaging in the 0.3-10 keV band across a 450 arcminute$^2$ field of view, with an order of magnitude improvement in sensitivity. AXIS utilizes breakthroughs in the construction of lightweight segmented X-ray optics usin… ▽ More The Advanced X-ray Imaging Satellite (AXIS) is a Probe-class concept that will build on the legacy of the Chandra X-ray Observatory by providing low-background, arcsecond-resolution imaging in the 0.3-10 keV band across a 450 arcminute$^2$ field of view, with an order of magnitude improvement in sensitivity. AXIS utilizes breakthroughs in the construction of lightweight segmented X-ray optics using single-crystal silicon, and developments in the fabrication of large-format, small-pixel, high readout rate CCD detectors with good spectral resolution, allowing a robust and cost-effective design. Further, AXIS will be responsive to target-of-opportunity alerts and, with onboard transient detection, will be a powerful facility for studying the time-varying X-ray universe, following on from the legacy of the Neil Gehrels (Swift) X-ray observatory that revolutionized studies of the transient X-ray Universe. In this paper, we present an overview of AXIS, highlighting the prime science objectives driving the AXIS concept and how the observatory design will achieve these objectives. △ Less

Submitted 1 November, 2023; originally announced November 2023.

Comments: Published in Proceedings of SPIE Optics & Photonics 2023, San Diego

New Black Hole Spin Values for Sagittarius A* Obtained with the Outflow Method

Authors: Ruth A. Daly, Megan Donahue, Christopher P. O'Dea, Biny Sebastian, Daryl Haggard, Anan Lu

Abstract: Six archival Chandra observations are matched with eight sets of radio data and studied in the context of the outflow method to measure and study the spin properties of $\rm{Sgr ~A^*}$. Three radio and X-ray data sets obtained simultaneously, or partially simultaneously, are identified as preferred for the purpose of measuring the spin properties of $\rm{Sgr ~A^*}$. Similar results are obtained wi… ▽ More Six archival Chandra observations are matched with eight sets of radio data and studied in the context of the outflow method to measure and study the spin properties of $\rm{Sgr ~A^*}$. Three radio and X-ray data sets obtained simultaneously, or partially simultaneously, are identified as preferred for the purpose of measuring the spin properties of $\rm{Sgr ~A^*}$. Similar results are obtained with other data sets. Results obtained with the preferred data sets are combined and indicate a weighted mean value of the spin function of $\rm{F} = 0.62 \pm 0.10$ and dimensionless spin angular momentum of $\rm{a_*} = 0.90 \pm 0.06$. The spin function translates into measurements of the black hole rotational mass, $\rm{M_{rot}}$, irreducible mass, $\rm{M_{irr}}$, and spin mass-energy available for extraction, $\rm{M_{spin}}$, relative to the total black hole dynamical mass, $\rm{M_{dyn}}$. Weighted mean values of $\rm{(M_{rot}/M_{dyn}) = (0.53 \pm 0.06)}$, $\rm{({M_{irr}/M_{dyn})} = (0.85 \pm 0.04)}$, $\rm{({M_{spin}/M_{dyn})} = (0.15 \pm 0.04)}$, $\rm{M_{rot} = (2.2 \pm 0.3) \times 10^6 ~M_{\odot}}$, $\rm{M_{irr} = (3.5 \pm 0.2) \times 10^6 ~M_{\odot}}$, and $\rm{M_{spin} = (6.2 \pm 1.6) \times 10^5 ~M_{\odot}}$ are obtained; of course $\rm{{(M_{rot}/M_{irr})} = (0.62 \pm 0.10)}$ since $\rm{{(M_{rot}/M_{irr})} = F}$. Values obtained for $\rm{Sgr ~A^*}$ are compared with those obtained for M87 based on the published spin function which indicate that M87 carries substantially more rotational energy and spin mass-energy relative to the total (i.e., dynamical) black hole mass, the irreducible black hole mass, and in absolute terms. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS on October 16, 2023

The JWST Galactic Center Survey -- A White Paper

Authors: Rainer Schoedel, Steve Longmore, Jonny Henshaw, Adam Ginsburg, John Bally, Anja Feldmeier, Matt Hosek, Francisco Nogueras Lara, Anna Ciurlo, Mélanie Chevance, J. M. Diederik Kruijssen, Ralf Klessen, Gabriele Ponti, Pau Amaro-Seoane, Konstantina Anastasopoulou, Jay Anderson, Maria Arias, Ashley T. Barnes, Cara Battersby, Giuseppe Bono, Lucía Bravo Ferres, Aaron Bryant, Miguel Cano Gonzáalez, Santi Cassisi, Leonardo Chaves-Velasquez , et al. (85 additional authors not shown)

Abstract: The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of… ▽ More The inner hundred parsecs of the Milky Way hosts the nearest supermassive black hole, largest reservoir of dense gas, greatest stellar density, hundreds of massive main and post main sequence stars, and the highest volume density of supernovae in the Galaxy. As the nearest environment in which it is possible to simultaneously observe many of the extreme processes shaping the Universe, it is one of the most well-studied regions in astrophysics. Due to its proximity, we can study the center of our Galaxy on scales down to a few hundred AU, a hundred times better than in similar Local Group galaxies and thousands of times better than in the nearest active galaxies. The Galactic Center (GC) is therefore of outstanding astrophysical interest. However, in spite of intense observational work over the past decades, there are still fundamental things unknown about the GC. JWST has the unique capability to provide us with the necessary, game-changing data. In this White Paper, we advocate for a JWST NIRCam survey that aims at solving central questions, that we have identified as a community: i) the 3D structure and kinematics of gas and stars; ii) ancient star formation and its relation with the overall history of the Milky Way, as well as recent star formation and its implications for the overall energetics of our galaxy's nucleus; and iii) the (non-)universality of star formation and the stellar initial mass function. We advocate for a large-area, multi-epoch, multi-wavelength NIRCam survey of the inner 100\,pc of the Galaxy in the form of a Treasury GO JWST Large Program that is open to the community. We describe how this survey will derive the physical and kinematic properties of ~10,000,000 stars, how this will solve the key unknowns and provide a valuable resource for the community with long-lasting legacy value. △ Less

Submitted 14 March, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: This White Paper will be updated when required (e.g. new authors joining, editing of content). Most recent update: 24 Oct 2023

Spectroscopic r-Process Abundance Retrieval for Kilonovae II: Lanthanides in the Inferred Abundance Patterns of Multi-Component Ejecta from the GW170817 Kilonova

Authors: Nicholas Vieira, John J. Ruan, Daryl Haggard, Nicole M. Ford, Maria R. Drout, Rodrigo Fernández

Abstract: In kilonovae, freshly-synthesized $r$-process elements imprint features on optical spectra, as observed in AT2017gfo, the counterpart to the GW170817 binary neutron star merger. However, measuring the $r$-process compositions of the merger ejecta is computationally challenging. Vieira et al. (2023) introduced Spectroscopic $r$-Process Abundance Retrieval for Kilonovae (SPARK), a software tool to i… ▽ More In kilonovae, freshly-synthesized $r$-process elements imprint features on optical spectra, as observed in AT2017gfo, the counterpart to the GW170817 binary neutron star merger. However, measuring the $r$-process compositions of the merger ejecta is computationally challenging. Vieira et al. (2023) introduced Spectroscopic $r$-Process Abundance Retrieval for Kilonovae (SPARK), a software tool to infer elemental abundance patterns of the ejecta, and associate spectral features with particular species. Previously, we applied SPARK to the 1.4 day spectrum of AT2017gfo and inferred its abundance pattern for the first time, characterized by electron fraction $Y_e=0.31$, a substantial abundance of strontium, and a dearth of lanthanides and heavier elements. This ejecta is consistent with wind from a remnant hypermassive neutron star and/or accretion disk. We now extend our inference to spectra at 2.4 and 3.4 days, and test the need for multicomponent ejecta, where we stratify the ejecta in composition. The ejecta at 1.4 and 2.4 days is described by the same single blue component. At 3.4 days, a new redder component with lower $Y_e=0.16$ and a significant abundance of lanthanides emerges. This new redder component is consistent with dynamical ejecta and/or neutron-rich ejecta from a magnetized accretion disk. As expected from photometric modelling, this component emerges as the ejecta expands, the photosphere recedes, and the earlier bluer component dims. At 3.4 days, we find an ensemble of lanthanides, with the presence of cerium most concrete. This presence of lanthanides has important implications for the contribution of kilonovae to the $r$-process abundances observed in the Universe. △ Less

Submitted 22 December, 2023; v1 submitted 31 August, 2023; originally announced August 2023.

Comments: Accepted to ApJ; 19 pages, 8 figures in body; 10 pages, 7 figures in appendix

KilonovAE: Exploring Kilonova Spectral Features with Autoencoders

Authors: N. M. Ford, Nicholas Vieira, John J. Ruan, Daryl Haggard

Abstract: Kilonovae are likely a key site of heavy r-process element production in the Universe, and their optical/infrared spectra contain insights into both the properties of the ejecta and the conditions of the r-process. However, the event GW170817/AT2017gfo is the only kilonova so far with well-observed spectra. To understand the diversity of absorption features that might be observed in future kilonov… ▽ More Kilonovae are likely a key site of heavy r-process element production in the Universe, and their optical/infrared spectra contain insights into both the properties of the ejecta and the conditions of the r-process. However, the event GW170817/AT2017gfo is the only kilonova so far with well-observed spectra. To understand the diversity of absorption features that might be observed in future kilonovae spectra, we use the TARDIS Monte Carlo radiative transfer code to simulate a suite of optical spectra spanning a wide range of kilonova ejecta properties and r-process abundance patterns. To identify the most common and prominent absorption lines, we perform dimensionality reduction using an autoencoder, and we find spectra clusters in the latent space representation using a Bayesian Gaussian Mixture model. Our synthetic kilonovae spectra commonly display strong absorption by strontium Sr II, yttrium Y II, and zirconium Zr I - II, with strong lanthanide contributions at low electron fractions (Ye < 0.25). When a new kilonova is observed, our machine learning framework will provide context on the dominant absorption lines and key ejecta properties, helping to determine where this event falls within the larger 'zoo' of kilonovae spectra. △ Less

Submitted 18 January, 2024; v1 submitted 29 August, 2023; originally announced August 2023.

Comments: 17 pages, 13 figures, published by ApJ. Replaced to match the accepted version

Journal ref: ApJ 961 119 (2024)

A search for pulsars around Sgr A* in the first Event Horizon Telescope dataset

Authors: Pablo Torne, Kuo Liu, Ralph P. Eatough, Jompoj Wongphechauxsorn, James M. Cordes, Gregory Desvignes, Mariafelicia De Laurentis, Michael Kramer, Scott M. Ransom, Shami Chatterjee, Robert Wharton, Ramesh Karuppusamy, Lindy Blackburn, Michael Janssen, Chi-kwan Chan, Geoffrey B. Crew, Lynn D. Matthews, Ciriaco Goddi, Helge Rottmann, Jan Wagner, Salvador Sanchez, Ignacio Ruiz, Federico Abbate, Geoffrey C. Bower, Juan J. Salamanca , et al. (261 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission… ▽ More The Event Horizon Telescope (EHT) observed in 2017 the supermassive black hole at the center of the Milky Way, Sagittarius A* (Sgr A*), at a frequency of 228.1 GHz ($λ$=1.3 mm). The fundamental physics tests that even a single pulsar orbiting Sgr A* would enable motivate searching for pulsars in EHT datasets. The high observing frequency means that pulsars - which typically exhibit steep emission spectra - are expected to be very faint. However, it also negates pulse scattering, an effect that could hinder pulsar detections in the Galactic Center. Additionally, magnetars or a secondary inverse Compton emission could be stronger at millimeter wavelengths than at lower frequencies. We present a search for pulsars close to Sgr A* using the data from the three most-sensitive stations in the EHT 2017 campaign: the Atacama Large Millimeter/submillimeter Array, the Large Millimeter Telescope and the IRAM 30 m Telescope. We apply three detection methods based on Fourier-domain analysis, the Fast-Folding-Algorithm and single pulse search targeting both pulsars and burst-like transient emission; using the simultaneity of the observations to confirm potential candidates. No new pulsars or significant bursts were found. Being the first pulsar search ever carried out at such high radio frequencies, we detail our analysis methods and give a detailed estimation of the sensitivity of the search. We conclude that the EHT 2017 observations are only sensitive to a small fraction ($\lesssim$2.2%) of the pulsars that may exist close to Sgr A*, motivating further searches for fainter pulsars in the region. △ Less

Submitted 29 August, 2023; originally announced August 2023.

Comments: 33 pages, 7 figures, 6 Tables. Accepted for publication in ApJ

Explaining the GWSkyNet-Multi machine learning classifier predictions for gravitational-wave events

Authors: Nayyer Raza, Man Leong Chan, Daryl Haggard, Ashish Mahabal, Jess McIver, Thomas C. Abbott, Eitan Buffaz, Nicholas Vieira

Abstract: GWSkyNet-Multi is a machine learning model developed for classification of candidate gravitational-wave events detected by the LIGO and Virgo observatories. The model uses limited information released in the low-latency Open Public Alerts to produce prediction scores indicating whether an event is a merger of two black holes, a merger involving a neutron star, or a non-astrophysical glitch. This f… ▽ More GWSkyNet-Multi is a machine learning model developed for classification of candidate gravitational-wave events detected by the LIGO and Virgo observatories. The model uses limited information released in the low-latency Open Public Alerts to produce prediction scores indicating whether an event is a merger of two black holes, a merger involving a neutron star, or a non-astrophysical glitch. This facilitates time sensitive decisions about whether to perform electromagnetic follow-up of candidate events during LIGO-Virgo-KAGRA (LVK) observing runs. However, it is not well understood how the model is leveraging the limited information available to make its predictions. As a deep learning neural network, the inner workings of the model can be difficult to interpret, impacting our trust in its validity and robustness. We tackle this issue by systematically perturbing the model and its inputs to explain what underlying features and correlations it has learned for distinguishing the sources. We show that the localization area of the 2D sky maps and the computed coherence versus incoherence Bayes factors are used as strong predictors for distinguishing between real events and glitches. The estimated distance to the source is further used to discriminate between binary black hole mergers and mergers involving neutron stars. We leverage these findings to show that events misclassified by GWSkyNet-Multi in LVK's third observing run have distinct sky area, coherence factor, and distance values that influence the predictions and explain these misclassifications. The results help identify the model's limitations and inform potential avenues for further optimization. △ Less

Submitted 23 August, 2023; originally announced August 2023.

Comments: 22 pages, 11 figures, submitted to ApJ

Signatures of Massive Black Hole Merger Host Galaxies from Cosmological Simulations I: Unique Galaxy Morphologies in Imaging

Authors: Jaeden Bardati, John J. Ruan, Daryl Haggard, Michael Tremmel

Abstract: Low-frequency gravitational wave experiments such as the Laser Interferometer Space Antenna and pulsar timing arrays are expected to detect individual massive black hole (MBH) binaries and mergers. However, secure methods of identifying the exact host galaxy of each MBH merger amongst the large number of galaxies in the gravitational wave localization region are currently lacking. We investigate t… ▽ More Low-frequency gravitational wave experiments such as the Laser Interferometer Space Antenna and pulsar timing arrays are expected to detect individual massive black hole (MBH) binaries and mergers. However, secure methods of identifying the exact host galaxy of each MBH merger amongst the large number of galaxies in the gravitational wave localization region are currently lacking. We investigate the distinct morphological signatures of MBH merger host galaxies, using the Romulus25 cosmological simulation. We produce mock telescope images of 201 simulated galaxies in Romulus25 hosting recent MBH mergers, through stellar population synthesis and dust radiative transfer. Based on comparisons to mass- and redshift-matched control samples, we show that combining multiple morphological statistics via a linear discriminant analysis enables identification of the host galaxies of MBH mergers, with accuracies that increase with chirp mass and mass ratio. For mergers with high chirp masses (>10^8.2 Msun) and high mass ratios (>0.5), the accuracy of this approach reaches >80%, and does not decline for at least >1 Gyr after numerical merger. We argue that these trends arise because the most distinctive morphological characteristics of MBH merger and binary host galaxies are prominent classical bulges, rather than relatively short-lived morphological disturbances from their preceding galaxy mergers. Since these bulges are formed though major mergers of massive galaxies, they lead to (and become permanent signposts for) MBH binaries and mergers that have high chirp masses and mass ratios. Our results suggest that galaxy morphology can aid in identifying the host galaxies of future MBH binaries and mergers. △ Less

Submitted 13 October, 2023; v1 submitted 7 August, 2023; originally announced August 2023.

Comments: 19 pages, 10 figures. ApJ accepted version

A fundamental plane of black hole accretion at millimetre wavelengths

Authors: Ilaria Ruffa, Timothy A. Davis, Jacob S. Elford, Martin Bureau, Michele Cappellari, Jindra Gensior, Daryl Haggard, Satoru Iguchi, Federico Lelli, Fu-Heng Liang, Lijie Liu, Marc Sarzi, Thomas G. Williams, Hengyue Zhang

Abstract: We report the discovery of the ``mm fundamental plane of black-hole accretion'', which is a tight correlation between the nuclear 1 mm luminosity ($L_{\rm ν, mm}$), the intrinsic $2$ -- $10$~keV X-ray luminosity ($L_{\rm X,2-10}$) and the supermassive black hole (SMBH) mass ($M_{\rm BH}$) with an intrinsic scatter ($σ_{\rm int}$) of $0.40$ dex. The plane is found for a sample of 48 nearby galaxies… ▽ More We report the discovery of the ``mm fundamental plane of black-hole accretion'', which is a tight correlation between the nuclear 1 mm luminosity ($L_{\rm ν, mm}$), the intrinsic $2$ -- $10$~keV X-ray luminosity ($L_{\rm X,2-10}$) and the supermassive black hole (SMBH) mass ($M_{\rm BH}$) with an intrinsic scatter ($σ_{\rm int}$) of $0.40$ dex. The plane is found for a sample of 48 nearby galaxies, most of which are low-luminosity active galactic nuclei (LLAGN). Combining these sources with a sample of high-luminosity (quasar-like) nearby AGN, we find that the plane still holds. We also find that $M_{\rm BH}$ correlates with $L_{\rm ν, mm}$ at a highly significant level, although such correlation is less tight than the mm fundamental plane ($σ_{\rm int}=0.51$ dex). Crucially, we show that spectral energy distribution (SED) models for both advection-dominated accretion flows (ADAFs) and compact jets can explain the existence of these relations, which are not reproduced by the standard torus-thin accretion disc models usually associated to quasar-like AGN. The ADAF models reproduces the observed relations somewhat better than those for compact jets, although neither provides a perfect prediction. Our findings thus suggest that radiatively-inefficient accretion processes such as those in ADAFs or compact (and thus possibly young) jets may play a key role in both low- and high-luminosity AGN. This mm fundamental plane also offers a new, rapid method to (indirectly) estimate SMBH masses. △ Less

Submitted 6 November, 2023; v1 submitted 25 July, 2023; originally announced July 2023.

Comments: 9 pages, 5 figures

Massive Black Hole Binaries as LISA Precursors in the Roman High Latitude Time Domain Survey

Authors: Zoltán Haiman, Chengcheng Xin, Tamara Bogdanović, Pau Amaro Seoane, Matteo Bonetti, J. Andrew Casey-Clyde, Maria Charisi, Monica Colpi, Jordy Davelaar, Alessandra De Rosa, Daniel J. D'Orazio, Kate Futrowsky, Poshak Gandhi, Alister W. Graham, Jenny E. Greene, Melanie Habouzit, Daryl Haggard, Kelly Holley-Bockelmann, Xin Liu, Alberto Mangiagli, Alessandra Mastrobuono-Battisti, Sean McGee, Chiara M. F. Mingarelli, Rodrigo Nemmen, Antonella Palmese , et al. (5 additional authors not shown)

Abstract: With its capacity to observe $\sim 10^{5-6}$ faint active galactic nuclei (AGN) out to redshift $z\approx 6$, Roman is poised to reveal a population of $10^{4-6}\, {\rm M_\odot}$ black holes during an epoch of vigorous galaxy assembly. By measuring the light curves of a subset of these AGN and looking for periodicity, Roman can identify several hundred massive black hole binaries (MBHBs) with 5-12… ▽ More With its capacity to observe $\sim 10^{5-6}$ faint active galactic nuclei (AGN) out to redshift $z\approx 6$, Roman is poised to reveal a population of $10^{4-6}\, {\rm M_\odot}$ black holes during an epoch of vigorous galaxy assembly. By measuring the light curves of a subset of these AGN and looking for periodicity, Roman can identify several hundred massive black hole binaries (MBHBs) with 5-12 day orbital periods, which emit copious gravitational radiation and will inevitably merge on timescales of $10^{3-5}$ years. During the last few months of their merger, such binaries are observable with the Laser Interferometer Space Antenna (LISA), a joint ESA/NASA gravitational wave mission set to launch in the mid-2030s. Roman can thus find LISA precursors, provide uniquely robust constraints on the LISA source population, help identify the host galaxies of LISA mergers, and unlock the potential of multi-messenger astrophysics with massive black hole binaries. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: White Paper for the Nancy Grace Roman Space Telescope's Core Community Surveys (https://roman.gsfc.nasa.gov/science/ccs_white_papers.html)

The X-ray Point Source Population Hosted by Globular Clusters in the Elliptical Galaxy NGC 4261

Authors: Sneha Nair, Kristen Dage, Daryl Haggard, Arunav Kundu, Richard M. Plotkin, Katherine L. Rhode, Stephen E. Zepf

Abstract: Utilising archival Chandra X-ray Observatory data and Hubble Space Telescope globular cluster catalogues, we probe the time-domain properties of the low mass X-ray binary population in the elliptical galaxy NGC 4261. Of the 98 unique X-ray sources identified in this study, 62 sources are within the optical field of view and, of those, 33% are aligned with an optical cluster counterpart. We find tw… ▽ More Utilising archival Chandra X-ray Observatory data and Hubble Space Telescope globular cluster catalogues, we probe the time-domain properties of the low mass X-ray binary population in the elliptical galaxy NGC 4261. Of the 98 unique X-ray sources identified in this study, 62 sources are within the optical field of view and, of those, 33% are aligned with an optical cluster counterpart. We find twenty X-ray sources coincident with globular clusters; two are previously discovered ultra-luminous X-ray sources (ULXs) and eighteen are low mass X-ray binaries (GCLMXBs) with $L_X < 10^{39}$ erg s$^{-1}$. ULXs are a heterogeneous class of extremely bright X-ray binaries ($L_X > 10^{39}$ erg s$^{-1}$) and ULXs located in globular clusters (GCULXs) and may be indicators of black holes. Identifying these unusually X-ray bright sources and measuring their optical properties can provide valuable constraints on the progenitors of gravitational wave sources. We compare observations of these sources to the twenty previously-studied GCULXs from five other early-type galaxies, and find that GCULXs in NGC 4261 are of similar colour and luminosity and do not significantly deviate from the rest of the sample in terms of distance from the galaxy centre or X-ray luminosity. Both the GCULX and low mass X-ray binary (GCLMXB) populations of NGC 4261 show long term variability; the former may have implications for fast radio bursts originating in globular clusters and the latter will likely introduce additional scatter into the low mass end of GCLMXB X-ray luminosity functions. △ Less

Submitted 20 June, 2023; originally announced June 2023.

Comments: Accepted to MNRAS

Key Science Goals for the Next-Generation Event Horizon Telescope

Authors: Michael D. Johnson, Kazunori Akiyama, Lindy Blackburn, Katherine L. Bouman, Avery E. Broderick, Vitor Cardoso, R. P. Fender, Christian M. Fromm, Peter Galison, José L. Gómez, Daryl Haggard, Matthew L. Lister, Andrei P. Lobanov, Sera Markoff, Ramesh Narayan, Priyamvada Natarajan, Tiffany Nichols, Dominic W. Pesce, Ziri Younsi, Andrew Chael, Koushik Chatterjee, Ryan Chaves, Juliusz Doboszewski, Richard Dodson, Sheperd S. Doeleman , et al. (20 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) has led to the first images of a supermassive black hole, revealing the central compact objects in the elliptical galaxy M87 and the Milky Way. Proposed upgrades to this array through the next-generation EHT (ngEHT) program would sharply improve the angular resolution, dynamic range, and temporal coverage of the existing EHT observations. These improvements will u… ▽ More The Event Horizon Telescope (EHT) has led to the first images of a supermassive black hole, revealing the central compact objects in the elliptical galaxy M87 and the Milky Way. Proposed upgrades to this array through the next-generation EHT (ngEHT) program would sharply improve the angular resolution, dynamic range, and temporal coverage of the existing EHT observations. These improvements will uniquely enable a wealth of transformative new discoveries related to black hole science, extending from event-horizon-scale studies of strong gravity to studies of explosive transients to the cosmological growth and influence of supermassive black holes. Here, we present the key science goals for the ngEHT and their associated instrument requirements, both of which have been formulated through a multi-year international effort involving hundreds of scientists worldwide. △ Less

Submitted 21 April, 2023; originally announced April 2023.

Comments: 32 pages, 11 figures, accepted for publication in a special issue of Galaxies on the ngEHT (https://www.mdpi.com/journal/galaxies/special_issues/ngEHT_blackholes)

Journal ref: Galaxies 2023, 11(3), 61

Comparison of Polarized Radiative Transfer Codes used by the EHT Collaboration

Authors: Ben S. Prather, Jason Dexter, Monika Moscibrodzka, Hung-Yi Pu, Thomas Bronzwaer, Jordy Davelaar, Ziri Younsi, Charles F. Gammie, Roman Gold, George N. Wong, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Uwe Bach, Anne-Kathrin Baczko, David Ball, Mislav Baloković, John Barrett, Michi Bauböck, Bradford A. Benson, Dan Bintley , et al. (248 additional authors not shown)

Abstract: Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curve… ▽ More Interpretation of resolved polarized images of black holes by the Event Horizon Telescope (EHT) requires predictions of the polarized emission observable by an Earth-based instrument for a particular model of the black hole accretion system. Such predictions are generated by general relativistic radiative transfer (GRRT) codes, which integrate the equations of polarized radiative transfer in curved spacetime. A selection of ray-tracing GRRT codes used within the EHT collaboration is evaluated for accuracy and consistency in producing a selection of test images, demonstrating that the various methods and implementations of radiative transfer calculations are highly consistent. When imaging an analytic accretion model, we find that all codes produce images similar within a pixel-wise normalized mean squared error (NMSE) of 0.012 in the worst case. When imaging a snapshot from a cell-based magnetohydrodynamic simulation, we find all test images to be similar within NMSEs of 0.02, 0.04, 0.04, and 0.12 in Stokes I, Q, U , and V respectively. We additionally find the values of several image metrics relevant to published EHT results to be in agreement to much better precision than measurement uncertainties. △ Less

Submitted 21 March, 2023; originally announced March 2023.

Comments: Accepted for publication in ApJ

Monitoring observations of SMC X-1's excursions (MOOSE)-II: A new excursion accompanies spin-up acceleration

Authors: Chin-Ping Hu, Kristen C. Dage, William I. Clarkson, McKinley Brumback, Philip A. Charles, Daryl Haggard, Ryan C. Hickox, Tatehiro Mihara, Arash Bahramian, Rawan Karam, Wasundara Athukoralalage, Diego Altamirano, Joey Neilsen, Jamie Kennea

Abstract: SMC X-1 is a high-mass X-ray binary showing superorbital modulation with an unstable period. Previous monitoring shows three excursion events in 1996--1998, 2005--2007, and 2014--2016. The superorbital period drifts from >60 days to <40 days and then evolves back during an excursion. Here we report a new excursion event of SMC X-1 in 2020--2021, indicating that the superorbital modulation has an u… ▽ More SMC X-1 is a high-mass X-ray binary showing superorbital modulation with an unstable period. Previous monitoring shows three excursion events in 1996--1998, 2005--2007, and 2014--2016. The superorbital period drifts from >60 days to <40 days and then evolves back during an excursion. Here we report a new excursion event of SMC X-1 in 2020--2021, indicating that the superorbital modulation has an unpredictable, chaotic nature. We trace the spin-period evolution and find that the spin-up rate accelerated one year before the onset of this new excursion, which suggests a possible inside-out process connecting the spin-up acceleration and the superorbital excursion. This results in a deviation of the spin period residual, similar to the behaviour of the first excursion in 1996--1998. In further analysis of the pulse profile evolution, we find that the pulsed fraction shows a long-term evolution and may be connected to the superorbital excursion. These discoveries deepen the mystery of SMC X-1 because they cannot be solely interpreted by the warped disc model. Upcoming pointed observations and theoretical studies may improve our understanding of the detailed accretion mechanisms taking place. △ Less

Submitted 1 February, 2023; originally announced February 2023.

Comments: 7 pages, 3 figures, Accepted for publication in MNRAS

Debiasing Standard Siren Inference of the Hubble Constant with Marginal Neural Ratio Estimation

Authors: Samuel Gagnon-Hartman, John Ruan, Daryl Haggard

Abstract: Gravitational wave (GW) standard sirens may resolve the Hubble tension, provided that standard siren inference of $H_0$ is free from systematic biases. However, standard sirens from binary neutron star (BNS) mergers suffer from two sources of systematic bias, one arising from the anisotropy of GW emission, and the other from the anisotropy of electromagnetic (EM) emission from the kilonova. For an… ▽ More Gravitational wave (GW) standard sirens may resolve the Hubble tension, provided that standard siren inference of $H_0$ is free from systematic biases. However, standard sirens from binary neutron star (BNS) mergers suffer from two sources of systematic bias, one arising from the anisotropy of GW emission, and the other from the anisotropy of electromagnetic (EM) emission from the kilonova. For an observed sample of BNS mergers, the traditional Bayesian approach to debiasing involves the direct computation of the detection likelihood. This is infeasible for large samples of detected BNS merger due to the high dimensionality of the parameter space governing merger detection. In this study, we bypass this computation by fitting the Hubble constant to forward simulations of the observed GW and EM data under a simulation-based inference (SBI) framework using marginal neural ratio estimation. A key innovation of our method is the inclusion of BNS mergers which were only detected in GW, which allows for estimation of the bias introduced by EM anisotropy. Our method corrects for $\sim$90$\%$ of the bias in the inferred value of $H_0$ when telescope follow-up observations of BNS mergers have extensive tiling of the merger localization region, using known telescope sensitivities and assuming a model of kilonova emission. Our SBI-based method thus enables a debiased inference of the Hubble constant of BNS mergers, including both mergers with detected EM counterparts and those without. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: Accepted MNRAS, 14 pages, 10 figures

Globular Cluster Ultraluminous X-ray Sources in the Furthest Early-Type Galaxies

Authors: E. Thygesen, Y. Sun, J. Huang, K. C. Dage, S. E. Zepf, A. Kundu, D. Haggard, T. J. Maccarone

Abstract: Ultraluminous X-ray Sources (ULXs) in globular clusters are low mass X-ray binaries that achieve high X-ray luminosities through a currently uncertain accretion mechanism. Using archival Chandra and Hubble Space Telescope observations, we perform a volume-limited search ($\lesssim$ 70 Mpc) of 21 of the most massive ($>10^{11.5} M_\odot$) early-type galaxies to identify ULXs hosted by globular clus… ▽ More Ultraluminous X-ray Sources (ULXs) in globular clusters are low mass X-ray binaries that achieve high X-ray luminosities through a currently uncertain accretion mechanism. Using archival Chandra and Hubble Space Telescope observations, we perform a volume-limited search ($\lesssim$ 70 Mpc) of 21 of the most massive ($>10^{11.5} M_\odot$) early-type galaxies to identify ULXs hosted by globular cluster (GC) candidates. We find a total of 34 ULX candidates above the expected background within 5 times the effective radius of each galaxy, with 10 of these ($\sim29.4\%$) potentially hosted by a GC. A comparison of the spatial and luminosity distributions of these new candidate GC ULXs with previously identified GC ULXs shows that they are similar: both samples peak at $L_X \sim$ a few $\times 10^{39}$ erg/s and are typically located within a few effective radii of their host galaxies. △ Less

Submitted 14 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS. 10 pages, 4 figures, 5 tables

Far Ultra-Violet Insights Into NGC 1399's Globular Cluster Population

Authors: Kristen C. Dage, Yifan Sun, Arunav Kundu, Stephen E. Zepf, Daryl Haggard

Abstract: We investigate archival Hubble Space Telescope ACS/SBC F140LP observations of NGC~1399 to search for evidence of multiple stellar populations in extragalactic globular clusters. Enhanced FUV populations are thought to be indicators of He-enhanced second generation populations in globular clusters, specifically extreme/blue horizontal branch stars. Out of 149 globular clusters in the field of view,… ▽ More We investigate archival Hubble Space Telescope ACS/SBC F140LP observations of NGC~1399 to search for evidence of multiple stellar populations in extragalactic globular clusters. Enhanced FUV populations are thought to be indicators of He-enhanced second generation populations in globular clusters, specifically extreme/blue horizontal branch stars. Out of 149 globular clusters in the field of view, 58 have far ultraviolet (FUV) counterparts with magnitudes brighter than 28.5. Six of these FUV-deteced globular clusters are also detected in X-rays, including one ultraluminous X-ray source ($L_X > 10^{39}$ erg/s). While optically bright clusters corresponded to brighter FUV counterparts, we observe FUV emission from both metal-rich and metal-poor clusters, which implies that the FUV excess is not dependent on optical colour. We also find no evidence that the cluster size influences the FUV emission. The clusters with X-ray emission are not unusually FUV bright, which suggests that even the ultraluminous X-ray source does not provide significant FUV contributions. NGC 1399 is only the fourth galaxy to have its globular cluster system probed for evidence of FUV-enhanced populations, and we compare these clusters to previous studies of the Milky Way, M31, M87, and the brightest cluster in M81. These sources indicate that many globular clusters likely host extreme HB stars and/or second generation stars, and highlight the need for more complete FUV observations of extragalactic globular cluster systems. △ Less

Submitted 26 October, 2022; originally announced October 2022.

Comments: accepted to MNRAS

Spectroscopic r-Process Abundance Retrieval for Kilonovae I: The Inferred Abundance Pattern of Early Emission from GW170817

Authors: Nicholas Vieira, John J. Ruan, Daryl Haggard, Nicole Ford, Maria R. Drout, Rodrigo Fernández, N. R. Badnell

Abstract: Freshly-synthesized r-process elements in kilonovae ejecta imprint absorption features on optical spectra, as observed in the GW170817 binary neutron star merger. These spectral features encode insights into the physical conditions of the r-process and the origins of the ejected material, but associating features with particular elements and inferring the resultant abundance pattern is computation… ▽ More Freshly-synthesized r-process elements in kilonovae ejecta imprint absorption features on optical spectra, as observed in the GW170817 binary neutron star merger. These spectral features encode insights into the physical conditions of the r-process and the origins of the ejected material, but associating features with particular elements and inferring the resultant abundance pattern is computationally challenging. We introduce Spectroscopic r-Process Abundance Retrieval for Kilonovae (SPARK), a modular framework to perform Bayesian inference on kilonova spectra with the goals of inferring elemental abundance patterns and identifying absorption features at early times. SPARK inputs an atomic line list and abundance patterns from reaction network calculations into the TARDIS radiative transfer code. It then performs fast Bayesian inference on observed kilonova spectra by training a Gaussian process surrogate for the approximate posteriors of kilonova ejecta parameters, via active learning. We use the spectrum of GW170817 at 1.4 days to perform the first inference on a kilonova spectrum, and recover a complete abundance pattern. Our inference shows that this ejecta was generated by an r-process with either (1) high electron fraction Y_e ~ 0.35 and high entropy s/k_B ~ 25, or, (2) a more moderate Y_e ~ 0.30 and s/k_B ~ 14. These parameters are consistent with a shocked, polar dynamical component, and a viscously-driven outflow from a remnant accretion disk, respectively. We also recover previous identifications of strontium absorption at ~8000 AA, and tentatively identify yttrium and/or zirconium at < 4500 AA. Our approach will enable computationally-tractable inference on the spectra of future kilonovae discovered through multi-messenger observations. △ Less

Submitted 22 December, 2022; v1 submitted 14 September, 2022; originally announced September 2022.

Comments: Accepted to ApJ; 34 pages, 13 figures

Multiwavelength astrophysics of the blazar OJ 287 and the project MOMO

Authors: S. Komossa, A. Kraus, D. Grupe, M. L. Parker, A. Gonzalez, L. C. Gallo, M. A. Gurwell, S. Laine, S. Yao, S. Chandra, L. Dey, J. L. Gomez, A. Gopakumar. K. Hada, D. Haggard, A. R. Hollett, H. Jermak, S. Jorstad, T. P. Krichbaum, S. Markoff, C. McCall, J. Neilsen, M. Nowak

Abstract: We are carrying out the densest and longest multiyear, multiwavelength monitoring project of OJ 287 ever done. The project MOMO (Multiwavelength Observations and Modelling of OJ 287) covers wavelengths from the radio to the high-energy regime. A few selected observations are simultaneous with those of the Event Horizon Telescope (EHT). MOMO aims at understanding disk-jet physics and at testing pre… ▽ More We are carrying out the densest and longest multiyear, multiwavelength monitoring project of OJ 287 ever done. The project MOMO (Multiwavelength Observations and Modelling of OJ 287) covers wavelengths from the radio to the high-energy regime. A few selected observations are simultaneous with those of the Event Horizon Telescope (EHT). MOMO aims at understanding disk-jet physics and at testing predictions of the binary black hole scenario of OJ 287. Here, we present a discussion of extreme outburst and minima states in context, and then focus on the recent flux and spectral evolution between 2021 and May 2022, including an ongoing bright radio flare. Further, we show that there is no evidence for precursor flare activity in our optical-UV-X-ray light curves that would be associated with any secondary supermassive black hole (SMBH) disk impact and that was predicted to start as thermal flare on 2021 December 23. △ Less

Submitted 22 July, 2022; originally announced July 2022.

Comments: Submitted to AN. Contribution to the XMM-Newton Workshop 2022 (Madrid, Spain). 5 figures

Millimeter light curves of Sagittarius A* observed during the 2017 Event Horizon Telescope campaign

Authors: Maciek Wielgus, Nicola Marchili, Ivan Marti-Vidal, Garrett K. Keating, Venkatessh Ramakrishnan, Paul Tiede, Ed Fomalont, Sara Issaoun, Joey Neilsen, Michael A. Nowak, Lindy Blackburn, Charles F. Gammie, Ciriaco Goddi, Daryl Haggard, Daeyoung Lee, Monika Moscibrodzka, Alexandra J. Tetarenko, Geoffrey C. Bower, Chi-Kwan Chan, Koushik Chatterjee, Paul M. Chesler, Jason Dexter, Sheperd S. Doeleman, Boris Georgiev, Mark Gurwell , et al. (6 additional authors not shown)

Abstract: The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data s… ▽ More The Event Horizon Telescope (EHT) observed the compact radio source, Sagittarius A* (Sgr A*), in the Galactic Center on 2017 April 5-11 in the 1.3 millimeter wavelength band. At the same time, interferometric array data from the Atacama Large Millimeter/submillimeter Array and the Submillimeter Array were collected, providing Sgr A* light curves simultaneous with the EHT observations. These data sets, complementing the EHT very-long-baseline interferometry, are characterized by a cadence and signal-to-noise ratio previously unattainable for Sgr A* at millimeter wavelengths, and they allow for the investigation of source variability on timescales as short as a minute. While most of the light curves correspond to a low variability state of Sgr A*, the April 11 observations follow an X-ray flare, and exhibit strongly enhanced variability. All of the light curves are consistent with a red noise process, with a power spectral density (PSD) slope measured to be between -2 and -3 on timescales between 1 min and several hours. Our results indicate a steepening of the PSD slope for timescales shorter than 0.3 h. The spectral energy distribution is flat at 220 GHz and there are no time-lags between the 213 and 229 GHz frequency bands, suggesting low optical depth for the event horizon scale source. We characterize Sgr A*'s variability, highlighting the different behavior observed just after the X-ray flare, and use Gaussian process modeling to extract a decorrelation timescale and a PSD slope. We also investigate the systematic calibration uncertainties by analyzing data from independent data reduction pipelines. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Journal ref: The Astrophysical Journal Letters, 930:L19 (2022)

Monitoring Observations of SMC X-1's Excursions (MOOSE) I: Program Description and Initial High-State Spectral Results

Authors: Kristen C. Dage, McKinley Brumback, Joey Neilsen, Chin-Ping Hu, Diego Altamirano, Arash Bahramian, Philip A. Charles, William I. Clarkson, Daryl Haggard, Ryan C. Hickox, Jamie Kennea

Abstract: SMC X-1 has exhibited three super-orbital period excursions since the onset of X-ray monitoring beginning with RXTE's launch in 1995. NICER has recently probed a fourth observed excursion beginning in 2021 with our program Monitoring Observations of SMC X-1's Excursions (MOOSE). These sensitive new MOOSE data probe different super-orbital periods and phases within them. Spectral fits to the high-s… ▽ More SMC X-1 has exhibited three super-orbital period excursions since the onset of X-ray monitoring beginning with RXTE's launch in 1995. NICER has recently probed a fourth observed excursion beginning in 2021 with our program Monitoring Observations of SMC X-1's Excursions (MOOSE). These sensitive new MOOSE data probe different super-orbital periods and phases within them. Spectral fits to the high-state continuum during April 2021 to January 2022 show that the intrinsic spectral shapes are characterised by a soft (kT~0.19 keV) disc component and a hard (Gamma~0.7) power-law tail. When the 2021-2022 NICER observations, taken during an excursion, are compared to 2016 XMM-Newton observations (outside of an excursion), we find little evidence for intrinsic spectral variability across the high-states, but find evidence for a >3 sigma change in the absorption, although we caution that there may be calibration differences between the two instruments. Thus, over different lengths of super-orbital periods, we see little evidence for intrinsic spectral changes in the high-state. Upcoming studies of the pulse profiles may shed light on the mechanism behind the excursions. △ Less

Submitted 13 June, 2022; originally announced June 2022.

Comments: accepted to MNRAS

WISDOM project -- XI. Star Formation Efficiency in the Bulge of the AGN-host Galaxy NGC 3169 with SITELLE and ALMA

Authors: Anan Lu, Hope Boyce, Daryl Haggard, Martin Bureau, Fu-Heng Liang, Lijie Liu, Woorak Choi, Michele Cappellari, Laurent Chemin, Mélanie Chevance, Timothy A. Davis, Laurent Drissen, Jacob S. Elford, Jindra Gensior, J. M. Diederik Kruijssen, Thomas Martin, Etienne Massé, Carmelle Robert, Ilaria Ruffa, Laurie Rousseau-Nepton, Marc Sarzi, Gabriel Savard Thomas G. Williams

Abstract: The star formation efficiency (SFE) has been shown to vary across different environments, particularly within galactic starbursts and deep within the bulges of galaxies. Various quenching mechanisms may be responsible, ranging from galactic dynamics to feedback from active galactic nuclei (AGN). Here, we use spatially-resolved observations of warm ionised gas emission lines from the imaging Fourie… ▽ More The star formation efficiency (SFE) has been shown to vary across different environments, particularly within galactic starbursts and deep within the bulges of galaxies. Various quenching mechanisms may be responsible, ranging from galactic dynamics to feedback from active galactic nuclei (AGN). Here, we use spatially-resolved observations of warm ionised gas emission lines from the imaging Fourier transform spectrograph SITELLE at the Canada-France-Hawaii Telescope (CFHT) and cold molecular gas (CO(2-1)) from the Atacama Large Millimeter/sub-millimeter Array (ALMA) to study the SFE in the bulge of the AGN-host galaxy NGC 3169. After distinguishing star-forming regions from AGN-ionised regions using emission-line ratio diagnostics, we measure spatially-resolved molecular gas depletion times (τ_dep = 1/SFE) with a spatial resolution of \approx 100 pc within a galactocentric radius of 1.8 kpc. We identify a star-forming ring located at radii 1.25 \pm 0.6 kpc with an average τ_dep of 0.3 Gyr. At radii < 0.9 kpc, however, the molecular gas surface densities and depletion times increase with decreasing radius, the latter reaching approximately 2.3 Gyr at a radius \approx 500 pc. Based on analyses of the gas kinematics and comparisons with simulations, we identify AGN feedback, bulge morphology and dynamics as the possible causes of the radial profile of SFE observed in the central region of NGC 3169. △ Less

Submitted 7 June, 2022; originally announced June 2022.

Comments: 25 pages, 15 figures, accepted for publication in MNRAS

Using Machine Learning to Identify Extragalactic Globular Cluster Candidates from Ground-Based Photometric Surveys of M87

Authors: Emilia Barbisan, Jeff Huang, Kristen C. Dage, Daryl Haggard, Robin Arnason, Arash Bahramian, William I. Clarkson, Arunav Kundu, Stephen E. Zepf

Abstract: Globular clusters (GCs) have been at the heart of many longstanding questions in many sub-fields of astronomy and, as such, systematic identification of GCs in external galaxies has immense impacts. In this study, we take advantage of M87's well-studied GC system to implement supervised machine learning (ML) classification algorithms - specifically random forest and neural networks - to identify G… ▽ More Globular clusters (GCs) have been at the heart of many longstanding questions in many sub-fields of astronomy and, as such, systematic identification of GCs in external galaxies has immense impacts. In this study, we take advantage of M87's well-studied GC system to implement supervised machine learning (ML) classification algorithms - specifically random forest and neural networks - to identify GCs from foreground stars and background galaxies using ground-based photometry from the Canada-France-Hawai'i Telescope (CFHT). We compare these two ML classification methods to studies of "human-selected" GCs and find that the best performing random forest model can reselect 61.2% $\pm$ 8.0% of GCs selected from HST data (ACSVCS) and the best performing neural network model reselects 95.0% $\pm$ 3.4%. When compared to human-classified GCs and contaminants selected from CFHT data - independent of our training data - the best performing random forest model can correctly classify 91.0% $\pm$ 1.2% and the best performing neural network model can correctly classify 57.3% $\pm$ 1.1%. ML methods in astronomy have been receiving much interest as Vera C. Rubin Observatory prepares for first light. The observables in this study are selected to be directly comparable to early Rubin Observatory data and the prospects for running ML algorithms on the upcoming dataset yields promising results. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: 14 pages, 9 figures, accepted to MNRAS

Multi-wavelength Variability of Sagittarius A* in July 2019

Authors: H. Boyce, D. Haggard, G. Witzel, S. von Fellenberg, S. P. Willner, E. E. Becklin, T. Do, A. Eckart, G. G. Fazio, M. A. Gurwell, J. L. Hora, S. Markoff, M. R. Morris, J. Neilsen, M. Nowak, H. A. Smith, S. Zhang

Abstract: We report timing analysis of near-infrared (NIR), X-ray, and sub-millimeter (submm) data during a three-day coordinated campaign observing Sagittarius A*. Data were collected at 4.5 micron with the Spitzer Space Telescope, 2-8 keV with the Chandra X-ray Observatory, 3-70 keV with NuSTAR, 340 GHz with ALMA, and at 2.2 micron with the GRAVITY instrument on the Very Large Telescope Interferometer. Tw… ▽ More We report timing analysis of near-infrared (NIR), X-ray, and sub-millimeter (submm) data during a three-day coordinated campaign observing Sagittarius A*. Data were collected at 4.5 micron with the Spitzer Space Telescope, 2-8 keV with the Chandra X-ray Observatory, 3-70 keV with NuSTAR, 340 GHz with ALMA, and at 2.2 micron with the GRAVITY instrument on the Very Large Telescope Interferometer. Two dates show moderate variability with no significant lags between the submm and the infrared at 99% confidence. July 18 captured a moderately bright NIR flare (F_K ~ 15 mJy) simultaneous with an X-ray flare (F ~ 0.1 cts/s) that most likely preceded bright submm flux (F ~ 5.5 Jy) by about +34 (+14 -33) minutes at 99% confidence. The uncertainty in this lag is dominated by the fact that we did not observe the peak of the submm emission. A synchrotron source cooled through adiabatic expansion can describe a rise in the submm once the synchrotron-self-Compton NIR and X-ray peaks have faded. This model predicts high GHz and THz fluxes at the time of the NIR/X-ray peak and electron densities well above those implied from average accretion rates for Sgr A*. However, the higher electron density postulated in this scenario would be in agreement with the idea that 2019 was an extraordinary epoch with a heightened accretion rate. Since the NIR and X-ray peaks can also be fit by a non-thermal synchrotron source with lower electron densities, we cannot rule out an unrelated chance coincidence of this bright submm flare with the NIR/X-ray emission. △ Less

Submitted 24 March, 2022; originally announced March 2022.

Comments: Accepted for publication in The Astrophysical Journal

Astrophysics with the Laser Interferometer Space Antenna

Authors: Pau Amaro Seoane, Jeff Andrews, Manuel Arca Sedda, Abbas Askar, Quentin Baghi, Razvan Balasov, Imre Bartos, Simone S. Bavera, Jillian Bellovary, Christopher P. L. Berry, Emanuele Berti, Stefano Bianchi, Laura Blecha, Stephane Blondin, Tamara Bogdanović, Samuel Boissier, Matteo Bonetti, Silvia Bonoli, Elisa Bortolas, Katelyn Breivik, Pedro R. Capelo, Laurentiu Caramete, Federico Cattorini, Maria Charisi, Sylvain Chaty , et al. (134 additional authors not shown)

Abstract: The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery… ▽ More The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. △ Less

Submitted 25 May, 2023; v1 submitted 11 March, 2022; originally announced March 2022.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

Extended radio emission in the galaxy cluster MS 0735.6+7421 detected with the Karl G. Jansky Very Large Array

Authors: T. Bégin, J. Hlavacek-Larrondo, C. L. Rhea, M. Gendron-Marsolais, B. McNamara, R. J. van Weeren, A. Richard-Laferrière, L. Guité, M. Prasow-Émond, D. Haggard

Abstract: MS 0735.6+7421 ($z = 0.216$) is a massive cool core galaxy cluster hosting one of the most powerful active galactic nuclei (AGN) outbursts known. The radio jets of the AGN have carved out an unusually large pair of X-ray cavities, each reaching a diameter of $200$ kpc. This makes MS 0735.6+7421 a unique case to investigate active galactic nuclei feedback processes, as well as other cluster astroph… ▽ More MS 0735.6+7421 ($z = 0.216$) is a massive cool core galaxy cluster hosting one of the most powerful active galactic nuclei (AGN) outbursts known. The radio jets of the AGN have carved out an unusually large pair of X-ray cavities, each reaching a diameter of $200$ kpc. This makes MS 0735.6+7421 a unique case to investigate active galactic nuclei feedback processes, as well as other cluster astrophysics at radio wavelengths. We present new low-radio-frequency observations of MS 0735.6+7421 taken with the Karl G. Jansky Very Large Array (VLA): 5 hours of P-band ($224-480$ MHz) and 5 hours of L-band ($1-2$ GHz) observations, both in C configuration. Our VLA P-band ($224-480$ MHz) observations reveal the presence of a new diffuse radio component reaching a scale of $\sim$ $900$ kpc in the direction of the jets and of $\sim$ $500$ kpc in the direction perpendicular to the jets. This component is centered on the cluster core and has a radio power scaled at $1.4$ GHz of $P_{1.4\text{ GHz}} = (4\pm2)\times 10^{24}$ WHz$^{-1}$. Its properties are consistent with those expected from a radio mini-halo as seen in other massive cool core clusters, although it may also be associated with radio plasma that has diffused out of the X-ray cavities. Observations at higher spatial resolution are needed to fully characterize the properties and nature of this component. We also suggest that if radio mini-halos originate from jetted activity, we may be witnessing the early stages of this process. △ Less

Submitted 19 January, 2023; v1 submitted 2 February, 2022; originally announced February 2022.

Comments: 11 pages, 7 figures, submitted to MNRAS

Journal ref: Volume 519, February 2023, Pages 767-777

A Swift study of long-term changes in the X-ray flaring properties of Sagittarius A*

Authors: A. Andrés, J. van den Eijnden, N. Degenaar, P. A. Evans, K. Chatterjee, M. Reynolds, J. M. Miller, J. Kennea, R. Wijnands, S. Markoff, D. Altamirano, C. O. Heinke, A. Bahramian, G. Ponti, D. Haggard

Abstract: The radiative counterpart of the supermassive black hole at the Galactic Centre, Sagittarius A*, displays flaring emission in the X-ray band atop a steady, quiescent level. Flares are also observed in the near-infrared band. The physical process producing the flares is not fully understood and it is unclear if the flaring rate varies, although some recent works suggest it has reached unprecedented… ▽ More The radiative counterpart of the supermassive black hole at the Galactic Centre, Sagittarius A*, displays flaring emission in the X-ray band atop a steady, quiescent level. Flares are also observed in the near-infrared band. The physical process producing the flares is not fully understood and it is unclear if the flaring rate varies, although some recent works suggest it has reached unprecedented variability in recent years. Using over a decade of regular X-ray monitoring of Neil Gehrels Swift Observatory, we studied the variations in count rate of Sgr A* on time scales of years. We decomposed the X-ray emission into quiescent and flaring emission, modelled as a constant and power law process, respectively. We found that the complete, multi-year dataset cannot be described by a stationary distribution of flare fluxes, while individual years follow this model better. In three of the ten studied years, the data is consistent with a purely Poissonian quiescent distribution, while for five years only an upper limit of the flare flux distribution parameter could be determined. We find that these possible changes cannot be explained fully by the different number of observations per year. Combined, these results are instead consistent with a changing flaring rate of Sgr A*, appearing more active between 2006-2007 and 2017-2019, than between 2008-2012. Finally, we discuss this result in the context of flare models and the passing of gaseous objects, and discuss the extra statistical steps taken, for instance to deal with the background in the Swift observations. △ Less

Submitted 19 November, 2021; originally announced November 2021.

Comments: 14 pages, 10 figures, Accepted for publication in MNRAS

GWSkyNet-Multi: A Machine Learning Multi-Class Classifier for LIGO-Virgo Public Alerts

Authors: Thomas C. Abbott, Eitan Buffaz, Nicholas Vieira, Miriam Cabero, Daryl Haggard, Ashish Mahabal, Jess McIver

Abstract: Compact object mergers which produce both detectable gravitational waves and electromagnetic emission can provide valuable insights into the neutron star equation of state, the tension in the Hubble constant, and the origin of the r-process elements. However, electromagnetic follow-up of gravitational wave sources is complicated by false positive detections, sources which do not emit light, and th… ▽ More Compact object mergers which produce both detectable gravitational waves and electromagnetic emission can provide valuable insights into the neutron star equation of state, the tension in the Hubble constant, and the origin of the r-process elements. However, electromagnetic follow-up of gravitational wave sources is complicated by false positive detections, sources which do not emit light, and the transient nature of the associated electromagnetic emission. GWSkyNet-Multi is a machine learning model that attempts to resolve these issues by providing real-time predictions of the source of a gravitational wave detection. The model uses information from Open Public Alerts (OPAs) released by LIGO-Virgo within minutes of a gravitational wave detection. GWSkyNet was first introduced in Cabero et al. (2020) as a binary classifier and uses the OPA skymaps to classify sources as either astrophysical or as glitches. In this paper, we introduce GWSkyNet-Multi, a new version of GWSkyNet which further distinguishes sources as binary black hole mergers, mergers involving a neutron star, or non-astrophysical glitches in the detectors. GWSkyNet-Multi is a sequence of three one-versus-all classifiers trained using a class-balanced and physically-motivated source mass distribution. Training on this data set, we obtain test set accuracies of 94% for BBH-versus-all, 94% for NS-versus-all, and 95% for glitch-versus-all. We obtain an overall accuracy of 93% using a hierarchical classification scheme. Furthermore, we correctly identify 36 of the 40 gravitational wave detections from the first half of LIGO-Virgo's third observing run (O3a) and present predictions for O3b sources. We also briefly discuss the interpretability of our models. As gravitational wave detections increase in number and frequency, GWSkyNet-Multi will be a powerful tool for prioritizing successful electromagnetic follow-up. △ Less

Submitted 7 November, 2021; originally announced November 2021.

Comments: 12 pages, 7 figures, submitted to ApJ

Event Horizon Telescope observations of the jet launching and collimation in Centaurus A

Authors: Michael Janssen, Heino Falcke, Matthias Kadler, Eduardo Ros, Maciek Wielgus, Kazunori Akiyama, Mislav Baloković, Lindy Blackburn, Katherine L. Bouman, Andrew Chael, Chi-kwan Chan, Koushik Chatterjee, Jordy Davelaar, Philip G. Edwards, Christian M. Fromm, José L. Gómez, Ciriaco Goddi, Sara Issaoun, Michael D. Johnson, Junhan Kim, Jun Yi Koay, Thomas P. Krichbaum, Jun Liu, Elisabetta Liuzzo, Sera Markoff , et al. (215 additional authors not shown)

Abstract: Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supe… ▽ More Very-long-baseline interferometry (VLBI) observations of active galactic nuclei at millimeter wavelengths have the power to reveal the launching and initial collimation region of extragalactic radio jets, down to $10-100$ gravitational radii ($r_g=GM/c^2$) scales in nearby sources. Centaurus A is the closest radio-loud source to Earth. It bridges the gap in mass and accretion rate between the supermassive black holes (SMBHs) in Messier 87 and our galactic center. A large southern declination of $-43^{\circ}$ has however prevented VLBI imaging of Centaurus A below $λ1$cm thus far. Here, we show the millimeter VLBI image of the source, which we obtained with the Event Horizon Telescope at $228$GHz. Compared to previous observations, we image Centaurus A's jet at a tenfold higher frequency and sixteen times sharper resolution and thereby probe sub-lightday structures. We reveal a highly-collimated, asymmetrically edge-brightened jet as well as the fainter counterjet. We find that Centaurus A's source structure resembles the jet in Messier 87 on ${\sim}500r_g$ scales remarkably well. Furthermore, we identify the location of Centaurus A's SMBH with respect to its resolved jet core at $λ1.3$mm and conclude that the source's event horizon shadow should be visible at THz frequencies. This location further supports the universal scale invariance of black holes over a wide range of masses. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 27 pages, 9 figures. This is a post-peer-review, pre-copyedit version of an article published in Nature Astronomy. The final authenticated version is available online at: https://doi.org/10.1038/s41550-021-01417-w

Journal ref: Nature Astronomy, July 2021, Volume 5, p. 1017-1028

The Variability of the Black-Hole Image in M87 at the Dynamical Time Scale

Authors: Kaushik Satapathy, Dimitrios Psaltis, Feryal Ozel, Lia Medeiros, Sean T. Dougall, Chi-kwan Chan, Maciek Wielgus, Ben S. Prather, George N. Wong, Charles F. Gammie, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David R. Ball, Mislav Baloković, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell , et al. (213 additional authors not shown)

Abstract: The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expect… ▽ More The black-hole images obtained with the Event Horizon Telescope (EHT) are expected to be variable at the dynamical timescale near their horizons. For the black hole at the center of the M87 galaxy, this timescale (5-61 days) is comparable to the 6-day extent of the 2017 EHT observations. Closure phases along baseline triangles are robust interferometric observables that are sensitive to the expected structural changes of the images but are free of station-based atmospheric and instrumental errors. We explored the day-to-day variability in closure phase measurements on all six linearly independent non-trivial baseline triangles that can be formed from the 2017 observations. We showed that three triangles exhibit very low day-to-day variability, with a dispersion of $\sim3-5^\circ$. The only triangles that exhibit substantially higher variability ($\sim90-180^\circ$) are the ones with baselines that cross visibility amplitude minima on the $u-v$ plane, as expected from theoretical modeling. We used two sets of General Relativistic magnetohydrodynamic simulations to explore the dependence of the predicted variability on various black-hole and accretion-flow parameters. We found that changing the magnetic field configuration, electron temperature model, or black-hole spin has a marginal effect on the model consistency with the observed level of variability. On the other hand, the most discriminating image characteristic of models is the fractional width of the bright ring of emission. Models that best reproduce the observed small level of variability are characterized by thin ring-like images with structures dominated by gravitational lensing effects and thus least affected by turbulence in the accreting plasmas. △ Less

Submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for Publication in ApJ

Ultraluminous X-ray Sources in Seven Edge-On Spiral Galaxies

Authors: Kristen C. Dage, Noah Vowell, Erica Thygesen, Arash Bahramian, Daryl Haggard, Konstantinos Kovlakas, Arunav Kundu, Thomas J. Maccarone, Jay Strader, Ryan Urquhart, Stephen E. Zepf

Abstract: We investigate a sample of seven edge-on spiral galaxies using Chandra observations. Edge-on spiral galaxies allow us to clearly separate source associated with their star-forming regions versus the outer edges of the system; offering a clear advantage over other systems. We uncover a number of X-ray point sources across these galaxies, and after eliminating contaminating foreground and background… ▽ More We investigate a sample of seven edge-on spiral galaxies using Chandra observations. Edge-on spiral galaxies allow us to clearly separate source associated with their star-forming regions versus the outer edges of the system; offering a clear advantage over other systems. We uncover a number of X-ray point sources across these galaxies, and after eliminating contaminating foreground and background sources, we identify 12 candidate ultraluminous X-ray sources. All of these sources are projected onto the central regions, implying that the majority of ULXs in this sample of spiral galaxies are disk/bulge, and thus not halo sources. This also includes two transient ULXs, which may be long-duration transients and low mass X-ray binaries. This finding illustrates the need for further studies of transient ULXs. △ Less

Submitted 29 September, 2021; originally announced September 2021.

Comments: accepted to MNRAS

Constraining particle acceleration in Sgr A* with simultaneous GRAVITY, Spitzer, NuSTAR and Chandra observations

Authors: R. Abuter, A. Amorim, M. Bauböck, F. Baganoff, J. P. Berge, H. Boyce, H. Bonnet, W. Brandner, Y. Clénet, R. Davies, P. T. de Zeeuw, J. Dexter, Y. Dallilar, A. Drescher, A. Eckart, F. Eisenhauer, G. G. Fazio, N. M. Förster Schreiber, K. Foster, C. Gammie, P. Garcia, F. Gao, E. Gendron, R. Genzel, G. Ghisellini , et al. (59 additional authors not shown)

Abstract: We report the time-resolved spectral analysis of a bright near-infrared and moderate X-ray flare of Sgr A*. We obtained light curves in the $M$-, $K$-, and $H$-bands in the mid- and near-infrared and in the $2-8~\mathrm{keV}$ and $2-70~\mathrm{keV}$ bands in the X-ray. The observed spectral slope in the near-infrared band is $νL_ν\propto ν^{0.5\pm0.2}$; the spectral slope observed in the X-ray ban… ▽ More We report the time-resolved spectral analysis of a bright near-infrared and moderate X-ray flare of Sgr A*. We obtained light curves in the $M$-, $K$-, and $H$-bands in the mid- and near-infrared and in the $2-8~\mathrm{keV}$ and $2-70~\mathrm{keV}$ bands in the X-ray. The observed spectral slope in the near-infrared band is $νL_ν\propto ν^{0.5\pm0.2}$; the spectral slope observed in the X-ray band is $νL_ν\propto ν^{-0.7\pm0.5}$. We tested synchrotron and synchrotron self-Compton (SSC) scenarios. The observed near-infrared brightness and X-ray faintness, together with the observed spectral slopes, pose challenges for all models explored. We rule out a scenario in which the near-infrared emission is synchrotron emission and the X-ray emission is SSC. A one-zone model in which both the near-infrared and X-ray luminosity are produced by SSC and a model in which the luminosity stems from a cooled synchrotron spectrum can explain the flare. In order to describe the mean SED, both models require specific values of the maximum Lorentz factor $γ_{max}$, which however differ by roughly two orders of magnitude: the SSC model suggests that electrons are accelerated to $γ_{max}\sim 500$, while cooled synchrotron model requires acceleration up to $γ_{max}\sim5\times 10^{4}$. The SSC scenario requires electron densities of $10^{10}~\mathrm{cm^{-3}}$ much larger than typical ambient densities in the accretion flow, and thus require in an extraordinary accretion event. In contrast, assuming a source size of $1R_s$, the cooled synchrotron scenario can be realized with densities and magnetic fields comparable with the ambient accretion flow. For both models, the temporal evolution is regulated through the maximum acceleration factor $γ_{max}$, implying that sustained particle acceleration is required to explain at least a part of the temporal evolution of the flare. △ Less

Submitted 2 July, 2021; originally announced July 2021.

Comments: accepted for publication in Astronomy & Astrophysics; preview abstract shortened due to arXiv requirements

Journal ref: A&A 654, A22 (2021)

Constraints on black-hole charges with the 2017 EHT observations of M87*

Authors: Prashant Kocherlakota, Luciano Rezzolla, Heino Falcke, Christian M. Fromm, Michael Kramer, Yosuke Mizuno, Antonios Nathanail, Hector Olivares, Ziri Younsi, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley, Lindy Blackburn, Raymond Blundell, Wilfred Boland , et al. (212 additional authors not shown)

Abstract: Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87*… ▽ More Our understanding of strong gravity near supermassive compact objects has recently improved thanks to the measurements made by the Event Horizon Telescope (EHT). We use here the M87* shadow size to infer constraints on the physical charges of a large variety of nonrotating or rotating black holes. For example, we show that the quality of the measurements is already sufficient to rule out that M87* is a highly charged dilaton black hole. Similarly, when considering black holes with two physical and independent charges, we are able to exclude considerable regions of the space of parameters for the doubly-charged dilaton and the Sen black holes. △ Less

Submitted 19 May, 2021; originally announced May 2021.

Comments: 15 pages, 4 figures, published in PRD on May 19

The Polarized Image of a Synchrotron Emitting Ring of Gas Orbiting a Black Hole

Authors: Ramesh Narayan, Daniel C. M. Palumbo, Michael D. Johnson, Zachary Gelles, Elizabeth Himwich, Dominic O. Chang, Angelo Ricarte, Jason Dexter, Charles F. Gammie, Andrew A. Chael, The Event Horizon Telescope Collaboration, :, Kazunori Akiyama, Antxon Alberdi, Walter Alef, Juan Carlos Algaba, Richard Anantua, Keiichi Asada, Rebecca Azulay, Anne-Kathrin Baczko, David Ball, Mislav Balokovic, John Barrett, Bradford A. Benson, Dan Bintley , et al. (215 additional authors not shown)

Abstract: Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equ… ▽ More Synchrotron radiation from hot gas near a black hole results in a polarized image. The image polarization is determined by effects including the orientation of the magnetic field in the emitting region, relativistic motion of the gas, strong gravitational lensing by the black hole, and parallel transport in the curved spacetime. We explore these effects using a simple model of an axisymmetric, equatorial accretion disk around a Schwarzschild black hole. By using an approximate expression for the null geodesics derived by Beloborodov (2002) and conservation of the Walker-Penrose constant, we provide analytic estimates for the image polarization. We test this model using currently favored general relativistic magnetohydrodynamic simulations of M87*, using ring parameters given by the simulations. For a subset of these with modest Faraday effects, we show that the ring model broadly reproduces the polarimetric image morphology. Our model also predicts the polarization evolution for compact flaring regions, such as those observed from Sgr A* with GRAVITY. With suitably chosen parameters, our simple model can reproduce the EVPA pattern and relative polarized intensity in Event Horizon Telescope images of M87*. Under the physically motivated assumption that the magnetic field trails the fluid velocity, this comparison is consistent with the clockwise rotation inferred from total intensity images. △ Less

Submitted 13 May, 2021; v1 submitted 4 May, 2021; originally announced May 2021.

Comments: 29 pages, 11 figures, published in ApJ on May 3

Journal ref: ApJ 912 35 (2021)

X-ray spectral components of the blazar and binary black hole candidate OJ 287 (2005-2020)

Authors: S. Komossa, D. Grupe, M. L. Parker, J. L. Gómez, M. J. Valtonen, M. A. Nowak, S. G. Jorstad, D. Haggard, S. Chandra, S. Ciprini, L. Dey, A. Gopakumar, K. Hada, S. Markoff, J. Neilsen

Abstract: We present a comprehensive analysis of all XMM-Newton spectra of OJ 287 spanning 15 years of X-ray spectroscopy of this bright blazar. We also report the latest results from our dedicated Swift UVOT and XRT monitoring of OJ 287 which started in 2015, along with all earlier public Swift data since 2005. During this time interval, OJ 287 was caught in extreme minima and outburst states. Its X-ray sp… ▽ More We present a comprehensive analysis of all XMM-Newton spectra of OJ 287 spanning 15 years of X-ray spectroscopy of this bright blazar. We also report the latest results from our dedicated Swift UVOT and XRT monitoring of OJ 287 which started in 2015, along with all earlier public Swift data since 2005. During this time interval, OJ 287 was caught in extreme minima and outburst states. Its X-ray spectrum is highly variable and encompasses all states seen in blazars from very flat to exceptionally steep. The spectrum can be decomposed into three spectral components: Inverse Compton (IC) emission dominant at low-states, super-soft synchrotron emission which becomes increasingly dominant as OJ 287 brightens, and an intermediately-soft (Gamma_x=2.2) additional component seen at outburst. This last component extends beyond 10 keV and plausibly represents either a second synchrotron/IC component and/or a temporary disk corona of the primary supermassive black hole (SMBH). Our 2018 XMM-Newton observation, quasi-simultaneous with the Event Horizon Telescope observation of OJ 287, is well described by a two-component model with a hard IC component of Gamma_x=1.5 and a soft synchrotron component. Low-state spectra limit any long-lived accretion disk/corona contribution in X-rays to a very low value of L_x/L_Edd < 5.6 times 10^(-4) (for M_(BH, primary) = 1.8 times 10^10 M_sun). Some implications for the binary SMBH model of OJ 287 are discussed. △ Less

Submitted 4 May, 2021; originally announced May 2021.

Comments: 14 pages. MNRAS, in press

Supermassive binary black holes and the case of OJ 287

Authors: S. Komossa, S. Ciprini, L. Dey, L. C. Gallo, J. L. Gomez, A. Gonzalez, D. Grupe, A. Kraus, S. J. Laine, M. L. Parker, M. J. Valtonen, S. Chandra, A. Gopakumar, D. Haggard, M. A. Nowak

Abstract: Supermassive binary black holes (SMBBHs) are laboratories par excellence for relativistic effects, including precession effects in the Kerr metric and the emission of gravitational waves. Binaries form in the course of galaxy mergers, and are a key component in our understanding of galaxy evolution. Dedicated searches for SMBBHs in all stages of their evolution are therefore ongoing and many syste… ▽ More Supermassive binary black holes (SMBBHs) are laboratories par excellence for relativistic effects, including precession effects in the Kerr metric and the emission of gravitational waves. Binaries form in the course of galaxy mergers, and are a key component in our understanding of galaxy evolution. Dedicated searches for SMBBHs in all stages of their evolution are therefore ongoing and many systems have been discovered in recent years. Here we provide a review of the status of observations with a focus on the multiwavelength detection methods and the underlying physics. Finally, we highlight our ongoing, dedicated multiwavelength program MOMO (for Multiwavelength Observations and Modelling of OJ 287). OJ 287 is one of the best candidates to date for hosting a sub-parsec SMBBH. The MOMO program carries out a dense monitoring at >13 frequencies from radio to X-rays and especially with Swift since 2015. Results so far included: (1) The detection of two major UV-X-ray outbursts with Swift in 2016/17 and 2020; exhibiting softer-when-brighter behaviour. The non-thermal nature of the outbursts was clearly established and shown to be synchrotron radiation. (2) Swift multi-band dense coverage and XMM-Newton spectroscopy during EHT campaigns caught OJ 287 at an intermediate flux level with synchrotron and IC spectral components. (3) Discovery of a remarkable, giant soft X-ray excess with XMM and NuSTAR during the 2020 outburst. (4) Spectral evidence (at 2sigma) for a relativistically shifted iron absorption line in 2020. (5) The non-thermal 2020 outburst is consistent with an after-flare predicted by the SMBBH model of OJ 287. △ Less

Submitted 26 April, 2021; originally announced April 2021.

Comments: Invited review at the COST Special Session on "The gravitational-wave Universe" during the XIX Serbian Astronomical Conference (Oct. 2020). 14 pages incl. 7 figures

Broadband Multi-wavelength Properties of M87 during the 2017 Event Horizon Telescope Campaign

Authors: J. C. Algaba, J. Anczarski, K. Asada, M. Balokovic, S. Chandra, Y. -Z. Cui, A. D. Falcone, M. Giroletti, C. Goddi, K. Hada, D. Haggard, S. Jorstad, A. Kaur, T. Kawashima, G. Keating, J. -Y. Kim, M. Kino, S. Komossa, E. V. Kravchenko, T. P. Krichbaum, S. -S. Lee, R. -S. Lu, M. Lucchini, S. Markoff, J. Neilsen , et al. (14 additional authors not shown)

Abstract: In 2017, the Event Horizon Telescope (EHT) Collaboration succeeded in capturing the first direct image of the center of the M87 galaxy. The asymmetric ring morphology and size are consistent with theoretical expectations for a weakly accreting supermassive black hole of mass approximately 6.5 x 10^9 M_solar. The EHTC also partnered with several international facilities in space and on the ground,… ▽ More In 2017, the Event Horizon Telescope (EHT) Collaboration succeeded in capturing the first direct image of the center of the M87 galaxy. The asymmetric ring morphology and size are consistent with theoretical expectations for a weakly accreting supermassive black hole of mass approximately 6.5 x 10^9 M_solar. The EHTC also partnered with several international facilities in space and on the ground, to arrange an extensive, quasi-simultaneous multi-wavelength campaign. This Letter presents the results and analysis of this campaign, as well as the multi-wavelength data as a legacy data repository. We captured M87 in a historically low state, and the core flux dominates over HST-1 at high energies, making it possible to combine core flux constraints with the more spatially precise very long baseline interferometry data. We present the most complete simultaneous multi-wavelength spectrum of the active nucleus to date, and discuss the complexity and caveats of combining data from different spatial scales into one broadband spectrum. We apply two heuristic, isotropic leptonic single-zone models to provide insight into the basic source properties, but conclude that a structured jet is necessary to explain M87's spectrum. We can exclude that the simultaneous gamma-ray emission is produced via inverse Compton emission in the same region producing the EHT mm-band emission, and further conclude that the gamma-rays can only be produced in the inner jets (inward of HST-1) if there are strongly particle-dominated regions. Direct synchrotron emission from accelerated protons and secondaries cannot yet be excluded. △ Less

Submitted 14 April, 2021; originally announced April 2021.

Comments: 43 pages, 19 figures, 10 tables; complete author list available in manuscript; The Astrophysical Journal Letters, 2021, 911, L11; publication doi: 3847/2041-8213/abef71, data doi: 10.25739/mhh2-cw46

Three Ultraluminous X-ray Sources Hosted by Globular Clusters in NGC 1316

Authors: Kristen C. Dage, Arunav Kundu, Erica Thygesen, Arash Bahramian, Daryl Haggard, Jimmy A. Irwin, Thomas J. Maccarone, Sneha Nair, Mark B. Peacock, Jay Strader, Stephen E. Zepf

Abstract: We have identified three ultraluminous X-ray sources (ULXs) hosted by globular clusters (GCs) within NGC 1316's stellar system. These discoveries bring the total number of known ULXs in GCs up to 20. We find that the X-ray spectra of the three new sources do not deviate from the established pattern of spectral behaviour of the other known GC ULXs. The consistency of the X-ray spectral behaviour fo… ▽ More We have identified three ultraluminous X-ray sources (ULXs) hosted by globular clusters (GCs) within NGC 1316's stellar system. These discoveries bring the total number of known ULXs in GCs up to 20. We find that the X-ray spectra of the three new sources do not deviate from the established pattern of spectral behaviour of the other known GC ULXs. The consistency of the X-ray spectral behaviour for these sources points to multiple paths of formation and evolution mechanisms for these rare and unique sources. Using the now larger sample of GC ULXs, we compare the optical properties of the entire known population of GC ULXs to other GCs across five galaxies and find that the properties of clusters that host ULXs are quite different from the typical clusters. Lastly, any trend of GC ULXs being preferentially hosted by metal-rich clusters is not strongly significant in this sample. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Comments: accepted to MNRAS

Probing the disk-corona systems and broad line regions of changing-look quasars with X-ray and optical observations

Authors: Xiangyu Jin, John J. Ruan, Daryl Haggard, Marie-Joëlle Gingras, Joseph Hountalas, Chelsea L. MacLeod, Scott F. Anderson, Anh Doan, Michael Eracleous, Paul J. Green, Jessie C. Runnoe

Abstract: "Changing-look" quasars are a new class of highly variable active galactic nuclei that have changed their spectral type over surprisingly short timescales of just a few years. The origin of this phenomenon is debated, but is likely to reflect some change in the accretion flow. To investigate the disk-corona systems in these objects, we measure optical/UV-X-ray spectral indices ($α_{\rm OX}$) and E… ▽ More "Changing-look" quasars are a new class of highly variable active galactic nuclei that have changed their spectral type over surprisingly short timescales of just a few years. The origin of this phenomenon is debated, but is likely to reflect some change in the accretion flow. To investigate the disk-corona systems in these objects, we measure optical/UV-X-ray spectral indices ($α_{\rm OX}$) and Eddington ratios ($λ_{\rm Edd}$) of ten previously-discovered changing-look quasars at two or more epochs. By comparing these data with simulated results based on the behavior of X-ray binaries, we find possible similarities in spectral indices below 1% Eddington ratio. We further investigate the Eddington ratios of changing-look quasars before and after their spectral type changes, and find that changing-look quasars cross the 1% Eddington ratio boundary when their broad emission lines disappear/emerge. This is consistent with the disk-wind model as the origin of broad emission lines. △ Less

Submitted 3 March, 2021; originally announced March 2021.

Comments: 31 pages, 9 figures, accepted for publication in ApJ

Rapid Variability of Sgr A* across the Electromagnetic Spectrum

Authors: G. Witzel, G. Martinez, S. P. Willner, E. E. Becklin, 4 H. Boyce, T. Do, A. Eckart, G. G. Fazio, A. Ghez, M. A. Gurwell, D. Haggard, R. Herrero-Illana, J. L. Hora, Z. Li, J. Liu, N. Marchili, Mark R. Morris, Howard A. Smith, M. Subroweit, J. A. Zensus

Abstract: Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We have analyzed a comprehensive submillimeter (including new observations simultaneous with NIR monitoring), NIR, and 2-8 keV dataset. Submillimeter variations tend to lag those in the NIR by $\sim$30 minutes. An approximate Bayesian computation (ABC)… ▽ More Sagittarius A* (Sgr A*) is the variable radio, near-infrared (NIR), and X-ray source associated with accretion onto the Galactic center black hole. We have analyzed a comprehensive submillimeter (including new observations simultaneous with NIR monitoring), NIR, and 2-8 keV dataset. Submillimeter variations tend to lag those in the NIR by $\sim$30 minutes. An approximate Bayesian computation (ABC) fit to the X-ray first-order structure function shows significantly less power at short timescales in the X-rays than in the NIR. Less X-ray variability at short timescales combined with the observed NIR-X-ray correlations means the variability can be described as the result of two strictly correlated stochastic processes, the X-ray process being the low-pass-filtered version of the NIR process. The NIR--X-ray linkage suggests a simple radiative model: a compact, self-absorbed synchrotron sphere with high-frequency cutoff close to NIR frequencies plus a synchrotron self-Compton scattering component at higher frequencies. This model, with parameters fit to the submillimeter, NIR, and X-ray structure functions, reproduces the observed flux densities at all wavelengths, the statistical properties of all light curves, and the time lags between bands. The fit also gives reasonable values for physical parameters such as magnetic flux density $B\approx13$ G, source size $L \approx2.2R_{S}$, and high-energy electron density $n_{e}\approx4\times10^{7}$ cm$^{-3}$. An animation illustrates typical light curves, and we make public the parameter chain of our Bayesian analysis, the model implementation, and the visualization code. △ Less

Submitted 7 June, 2021; v1 submitted 18 November, 2020; originally announced November 2020.

Comments: Accepted by ApJS, June 3, 2021. The journal version of figure 21 is animated. The animation can also be found here: https://doi.org/10.17617/1.kctx3s25. This version (version 2) has been revised according to the referee's suggestions and includes optimized figures and some changes and corrections of the text; no scientific conclusions have changed