Showing 1–50 of 92 results for author: Gültekin, K

A Stellar Dynamical Mass Measurement of the Supermassive Black Hole in NGC 3258

Authors: Thomas K. Waters, Kayhan Gültekin, Karl Gebhardt, Neil Nagar, Vanessa Ávila

Abstract: We present a stellar dynamical mass measurement of the supermassive black hole in the elliptical (E1) galaxy NGC 3258. Our findings are based on Integral Field Unit spectroscopy from the Multi Unit Spectroscopic Explorer observations in Narrow Field Mode with adaptive optics and in the MUSE Wide Field Mode, from which we extract kinematic information by fitting the Ca II and Mg $b$ triplets, respe… ▽ More We present a stellar dynamical mass measurement of the supermassive black hole in the elliptical (E1) galaxy NGC 3258. Our findings are based on Integral Field Unit spectroscopy from the Multi Unit Spectroscopic Explorer observations in Narrow Field Mode with adaptive optics and in the MUSE Wide Field Mode, from which we extract kinematic information by fitting the Ca II and Mg $b$ triplets, respectively. Using axisymmetric, three-integral Schwarzschild orbit-library models, we fit the observed line-of-sight velocity distributions to infer the supermassive black hole mass, the $H$-band mass-to-light ratio, the asymptotic circular velocity, and the dark matter halo scale radius of the galaxy. We report a black hole mass of $(2.2 \pm 0.2)\times10^9 \ \rm M_{\scriptscriptstyle\odot}$ at an assumed distance of $31.9 \ \rm Mpc$. This value is in close agreement with a previous measurement from ALMA CO observations. The consistency between these two measurements provides strong support for both the gas dynamical and stellar dynamical methods. △ Less

Submitted 20 June, 2024; originally announced June 2024.

Comments: 18 pages, 12 figures, 1 table

The NANOGrav 15 yr Data Set: Looking for Signs of Discreteness in the Gravitational-wave Background

Authors: Gabriella Agazie, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Lucas Brown, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Heling Deng, Timothy Dolch, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman, Nate Garver-Daniels , et al. (58 additional authors not shown)

Abstract: The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analyt… ▽ More The cosmic merger history of supermassive black hole binaries (SMBHBs) is expected to produce a low-frequency gravitational wave background (GWB). Here we investigate how signs of the discrete nature of this GWB can manifest in pulsar timing arrays through excursions from, and breaks in, the expected $f_{\mathrm{GW}}^{-2/3}$ power-law of the GWB strain spectrum. To do this, we create a semi-analytic SMBHB population model, fit to NANOGrav's 15 yr GWB amplitude, and with 1,000 realizations we study the populations' characteristic strain and residual spectra. Comparing our models to the NANOGrav 15 yr spectrum, we find two interesting excursions from the power-law. The first, at $2 \; \mathrm{nHz}$, is below our GWB realizations with $p$-value significance $p = 0.05$ to $0.06$ ($\approx 1.8 σ- 1.9 σ$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 σ- 2.1 σ$). We explore the properties of a loud SMBHB which could cause such an excursion. Our simulations also show that the expected number of SMBHBs decreases by three orders of magnitude, from $\sim 10^6$ to $\sim 10^3$, between $2\; \mathrm{nHz}$ and $20 \; \mathrm{nHz}$. This causes a break in the strain spectrum as the stochasticity of the background breaks down at $26^{+28}_{-19} \; \mathrm{nHz}$, consistent with predictions pre-dating GWB measurements. The diminished GWB signal from SMBHBs at frequencies above the $26$~nHz break opens a window for PTAs to detect continuous GWs from individual SMBHBs or GWs from the early universe. △ Less

Submitted 10 April, 2024; originally announced April 2024.

Comments: 10 pages, 8 figures, 1 appendix, submitted to ApJ

Searching for the Highest-z Dual AGN in the Deepest Chandra Surveys

Authors: Brandon Sandoval, Adi Foord, Steven W. Allen, Marta Volonteri, Aaron Stemo, Nianyi Chen, Tiziana Di Matteo, Kayhan Gultekin, Melanie Habouzit, Clara Puerto-Sanchez, Edmund Hodges-Kluck, Yohan Dubois

Abstract: We present an analysis searching for dual AGN among 62 high-redshift ($2.5 < z < 3.5$) X-ray sources selected from publicly available deep Chandra fields. We aim to quantify the frequency of dual AGN in the high-redshift Universe, which holds implications for black hole merger timescales and low-frequency gravitational wave detection rates. We analyze each X-ray source using BAYMAX, an analysis to… ▽ More We present an analysis searching for dual AGN among 62 high-redshift ($2.5 < z < 3.5$) X-ray sources selected from publicly available deep Chandra fields. We aim to quantify the frequency of dual AGN in the high-redshift Universe, which holds implications for black hole merger timescales and low-frequency gravitational wave detection rates. We analyze each X-ray source using BAYMAX, an analysis tool that calculates the Bayes factor for whether a given archival Chandra AGN is more likely a single or dual point source. We find no strong evidence for dual AGN in any individual source in our sample. We then increase our sensitivity to search for dual AGN across the sample by comparing our measured distribution of Bayes factors to that expected from a sample composed entirely of single point sources, and again find no evidence for dual AGN in the observed sample distribution. Although our analysis utilizes one of the largest Chandra catalogs of high-$z$ X-ray point sources available to study, the findings remain limited by the modest number of sources observed at the highest spatial resolution with Chandra and the typical count rates of the detected sources. Our non-detection allows us to place an upper-limit on the X-ray dual AGN fraction between $2.5<z<3.5$ of 4.8\%. Expanding substantially on these results at X-ray wavelengths will require future surveys spanning larger sky areas and extending to fainter fluxes than has been possible with Chandra. We illustrate the potential of the AXIS mission concept in this regard. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: 19 pages, 5 figures, 2 tales. Submitted to ApJ

Connecting Core Galaxy Properties to the Massive Black Hole Binary Population

Authors: CJ Harris, Kayhan Gultekin

Abstract: We investigate how the properties of massive black hole binaries influence the observed properties of core galaxies. We compare the observed trend in stellar mass deficit as a function of total stellar mass in the core galaxy with predicted trends in IllustrisTNG. We calculate mass deficits in simulated galaxies by applying sub-grid, post-processing physics based on the results of literature N-bod… ▽ More We investigate how the properties of massive black hole binaries influence the observed properties of core galaxies. We compare the observed trend in stellar mass deficit as a function of total stellar mass in the core galaxy with predicted trends in IllustrisTNG. We calculate mass deficits in simulated galaxies by applying sub-grid, post-processing physics based on the results of literature N-body experiments. We find the median value of the posterior distribution for the minimum binary mass ratio capable of creating a core is 0.7. For the gas mass fraction above which a core is erased we find a median value of 0.6. Thus low mass ratio binaries do not contribute to core formation and gas-rich mergers must lead to star formation within the nucleus, ultimately erasing the core. Such constraints have important implications for the overall massive black hole binary population, black hole-galaxy co-evolution, and the origin of the gravitational wave background. △ Less

Submitted 8 November, 2023; originally announced November 2023.

Comments: 10 pages, 12 figures, Accepted by MNRAS

The NANOGrav 15-year data set: Search for Transverse Polarization Modes in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Jeremy Baier, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest , et al. (74 additional authors not shown)

Abstract: Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correl… ▽ More Recently we found compelling evidence for a gravitational wave background with Hellings and Downs (HD) correlations in our 15-year data set. These correlations describe gravitational waves as predicted by general relativity, which has two transverse polarization modes. However, more general metric theories of gravity can have additional polarization modes which produce different interpulsar correlations. In this work we search the NANOGrav 15-year data set for evidence of a gravitational wave background with quadrupolar Hellings and Downs (HD) and Scalar Transverse (ST) correlations. We find that HD correlations are the best fit to the data, and no significant evidence in favor of ST correlations. While Bayes factors show strong evidence for a correlated signal, the data does not strongly prefer either correlation signature, with Bayes factors $\sim 2$ when comparing HD to ST correlations, and $\sim 1$ for HD plus ST correlations to HD correlations alone. However, when modeled alongside HD correlations, the amplitude and spectral index posteriors for ST correlations are uninformative, with the HD process accounting for the vast majority of the total signal. Using the optimal statistic, a frequentist technique that focuses on the pulsar-pair cross-correlations, we find median signal-to-noise-ratios of 5.0 for HD and 4.6 for ST correlations when fit for separately, and median signal-to-noise-ratios of 3.5 for HD and 3.0 for ST correlations when fit for simultaneously. While the signal-to-noise-ratios for each of the correlations are comparable, the estimated amplitude and spectral index for HD are a significantly better fit to the total signal, in agreement with our Bayesian analysis. △ Less

Submitted 18 October, 2023; originally announced October 2023.

Comments: 11 pages, 5 figures

ALMA gas-dynamical mass measurement of the supermassive black hole in the red nugget relic galaxy PGC 11179

Authors: Jonathan H. Cohn, Maeve Curliss, Jonelle L. Walsh, Kyle M. Kabasares, Benjamin D. Boizelle, Aaron J. Barth, Karl Gebhardt, Kayhan Gültekin, Akın Yıldırım, David Buote, Jeremy Darling, Andrew J. Baker, Luis Ho

Abstract: We present 0$.\!\!^{\prime\prime}22$-resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO(2$-$1) emission from the circumnuclear gas disk in the red nugget relic galaxy PGC 11179. The disk shows regular rotation, with projected velocities near the center of 400 km s$^{-1}$. We assume the CO emission originates from a dynamically cold, thin disk and fit gas-dynamical mo… ▽ More We present 0$.\!\!^{\prime\prime}22$-resolution Atacama Large Millimeter/submillimeter Array (ALMA) observations of CO(2$-$1) emission from the circumnuclear gas disk in the red nugget relic galaxy PGC 11179. The disk shows regular rotation, with projected velocities near the center of 400 km s$^{-1}$. We assume the CO emission originates from a dynamically cold, thin disk and fit gas-dynamical models directly to the ALMA data. In addition, we explore systematic uncertainties by testing the impacts of various model assumptions on our results. The supermassive black hole (BH) mass ($M_\mathrm{BH}$) is measured to be $M_\mathrm{BH} = (1.91\pm0.04$ [$1σ$ statistical] $^{+0.11}_{-0.51}$ [systematic])$\times 10^9$ $M_\odot$, and the $H$-band stellar mass-to-light ratio $M/L_H=1.620\pm0.004$ [$1σ$ statistical] $^{+0.211}_{-0.107}$ [systematic] $M_\odot/L_\odot$. This $M_\mathrm{BH}$ is consistent with the BH mass$-$stellar velocity dispersion relation but over-massive compared to the BH mass$-$bulge luminosity relation by a factor of 3.7. PGC 11179 is part of a sample of local compact early-type galaxies that are plausible relics of $z\sim2$ red nuggets, and its behavior relative to the scaling relations echoes that of three relic galaxy BHs previously measured with stellar dynamics. These over-massive BHs could suggest BHs gain most of their mass before their host galaxies do. However, our results could also be explained by greater intrinsic scatter at the high-mass end of the scaling relations, or by systematic differences in gas- and stellar-dynamical methods. Additional $M_\mathrm{BH}$ measurements in the sample, including independent cross-checks between molecular gas- and stellar-dynamical methods, will advance our understanding of the co-evolution of BHs and their host galaxies. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: 18 pages, 8 figures, 2 tables. Accepted in ApJ

Shock-Driven Periodic Variability in a Low-Mass-Ratio Supermassive Black Hole Binary

Authors: K. Whitley, A. Kuznetsova, K. Gültekin, M. Ruszkowski

Abstract: We investigate the time-varying electromagnetic emission of a low-mass-ratio supermassive black hole binary (SMBHB) embedded in a circumprimary disk, with a particular interest in variability of shocks driven by the binary. We perform a 2D, locally isothermal hydrodynamics simulation of a SMBHB with mass ratio $q=0.01$ and separation $a=100\;R_g$, using a physically self-consistent steady disk mod… ▽ More We investigate the time-varying electromagnetic emission of a low-mass-ratio supermassive black hole binary (SMBHB) embedded in a circumprimary disk, with a particular interest in variability of shocks driven by the binary. We perform a 2D, locally isothermal hydrodynamics simulation of a SMBHB with mass ratio $q=0.01$ and separation $a=100\;R_g$, using a physically self-consistent steady disk model. We estimate the electromagnetic variability from the system by monitoring accretion onto the secondary and using an artificial viscosity scheme to capture shocks and monitor the energy dissipated. The SMBHB produces a wide, eccentric gap in the disk, previously only observed for larger mass ratios, which we attribute to our disk model being much thinner ($H/R\approx0.01$ near the secondary) than is typical of previous works. The eccentric gap drives periodic accretion onto the secondary SMBH on a timescale matching the orbital period of the binary, $t_{\rm{bin}}\approx0.1\;\rm{yr}$, implying that the variable accretion regime of the SMBHB parameter space extends to lower mass ratios than previously established. Shocks driven by the binary are periodic, with a period matching the orbital period, and the shocks are correlated with the accretion rate, with peaks in the shock luminosity lagging peaks in the accretion rate by $0.43\;t_{\rm{bin}}$. We propose that the correlation of these quantities represents a useful identifier of SMBHB candidates, via observations of correlated variability in X-ray and UV monitoring of AGN, rather than single-waveband periodicity alone. △ Less

Submitted 18 October, 2023; v1 submitted 10 October, 2023; originally announced October 2023.

Comments: 12 pages, 8 figures, accepted by MNRAS

The NANOGrav 12.5-year data set: A computationally efficient eccentric binary search pipeline and constraints on an eccentric supermassive binary candidate in 3C 66B

Authors: Gabriella Agazie, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Belinda D. Cheeseboro, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Lankeswar Dey, Timothy Dolch, Justin A. Ellis, Robert D. Ferdman, Elizabeth C. Ferrara , et al. (63 additional authors not shown)

Abstract: The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for contin… ▽ More The radio galaxy 3C 66B has been hypothesized to host a supermassive black hole binary (SMBHB) at its center based on electromagnetic observations. Its apparent 1.05-year period and low redshift ($\sim0.02$) make it an interesting testbed to search for low-frequency gravitational waves (GWs) using Pulsar Timing Array (PTA) experiments. This source has been subjected to multiple searches for continuous GWs from a circular SMBHB, resulting in progressively more stringent constraints on its GW amplitude and chirp mass. In this paper, we develop a pipeline for performing Bayesian targeted searches for eccentric SMBHBs in PTA data sets, and test its efficacy by applying it on simulated data sets with varying injected signal strengths. We also search for a realistic eccentric SMBHB source in 3C 66B using the NANOGrav 12.5-year data set employing PTA signal models containing Earth term-only as well as Earth+Pulsar term contributions using this pipeline. Due to limitations in our PTA signal model, we get meaningful results only when the initial eccentricity $e_0<0.5$ and the symmetric mass ratio $η>0.1$. We find no evidence for an eccentric SMBHB signal in our data, and therefore place 95% upper limits on the PTA signal amplitude of $88.1\pm3.7$ ns for the Earth term-only and $81.74\pm0.86$ ns for the Earth+Pulsar term searches for $e_0<0.5$ and $η>0.1$. Similar 95% upper limits on the chirp mass are $(1.98 \pm 0.05) \times 10^9\,M_{\odot}$ and $(1.81 \pm 0.01) \times 10^9\,M_{\odot}$. These upper limits, while less stringent than those calculated from a circular binary search in the NANOGrav 12.5-year data set, are consistent with the SMBHB model of 3C 66B developed from electromagnetic observations. △ Less

Submitted 15 January, 2024; v1 submitted 29 September, 2023; originally announced September 2023.

Comments: 27 Pages, 10 Figures, 1 Table, Accepted for publication in ApJ

How to Detect an Astrophysical Nanohertz Gravitational-Wave Background

Authors: Bence Bécsy, Neil J. Cornish, Patrick M. Meyers, Luke Zoltan Kelley, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Tyler Cohen, James M. Cordes, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch , et al. (71 additional authors not shown)

Abstract: Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, nam… ▽ More Analysis of pulsar timing data have provided evidence for a stochastic gravitational wave background in the nHz frequency band. The most plausible source of such a background is the superposition of signals from millions of supermassive black hole binaries. The standard statistical techniques used to search for such a background and assess its significance make several simplifying assumptions, namely: i) Gaussianity; ii) isotropy; and most often iii) a power-law spectrum. However, a stochastic background from a finite collection of binaries does not exactly satisfy any of these assumptions. To understand the effect of these assumptions, we test standard analysis techniques on a large collection of realistic simulated datasets. The dataset length, observing schedule, and noise levels were chosen to emulate the NANOGrav 15-year dataset. Simulated signals from millions of binaries drawn from models based on the Illustris cosmological hydrodynamical simulation were added to the data. We find that the standard statistical methods perform remarkably well on these simulated datasets, despite their fundamental assumptions not being strictly met. They are able to achieve a confident detection of the background. However, even for a fixed set of astrophysical parameters, different realizations of the universe result in a large variance in the significance and recovered parameters of the background. We also find that the presence of loud individual binaries can bias the spectral recovery of the background if we do not account for them. △ Less

Submitted 1 December, 2023; v1 submitted 8 September, 2023; originally announced September 2023.

Comments: 14 pages, 8 figures, version matching published paper

Journal ref: ApJ 959 9 (2023)

Comparing recent PTA results on the nanohertz stochastic gravitational wave background

Authors: The International Pulsar Timing Array Collaboration, G. Agazie, J. Antoniadis, A. Anumarlapudi, A. M. Archibald, P. Arumugam, S. Arumugam, Z. Arzoumanian, J. Askew, S. Babak, M. Bagchi, M. Bailes, A. -S. Bak Nielsen, P. T. Baker, C. G. Bassa, A. Bathula, B. Bécsy, A. Berthereau, N. D. R. Bhat, L. Blecha, M. Bonetti, E. Bortolas, A. Brazier, P. R. Brook, M. Burgay , et al. (220 additional authors not shown)

Abstract: The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTA… ▽ More The Australian, Chinese, European, Indian, and North American pulsar timing array (PTA) collaborations recently reported, at varying levels, evidence for the presence of a nanohertz gravitational wave background (GWB). Given that each PTA made different choices in modeling their data, we perform a comparison of the GWB and individual pulsar noise parameters across the results reported from the PTAs that constitute the International Pulsar Timing Array (IPTA). We show that despite making different modeling choices, there is no significant difference in the GWB parameters that are measured by the different PTAs, agreeing within $1σ$. The pulsar noise parameters are also consistent between different PTAs for the majority of the pulsars included in these analyses. We bridge the differences in modeling choices by adopting a standardized noise model for all pulsars and PTAs, finding that under this model there is a reduction in the tension in the pulsar noise parameters. As part of this reanalysis, we "extended" each PTA's data set by adding extra pulsars that were not timed by that PTA. Under these extensions, we find better constraints on the GWB amplitude and a higher signal-to-noise ratio for the Hellings and Downs correlations. These extensions serve as a prelude to the benefits offered by a full combination of data across all pulsars in the IPTA, i.e., the IPTA's Data Release 3, which will involve not just adding in additional pulsars, but also including data from all three PTAs where any given pulsar is timed by more than as single PTA. △ Less

Submitted 1 September, 2023; originally announced September 2023.

Comments: 21 pages, 9 figures, submitted to ApJ

Bose-Einstein Condensation and Black Holes in Dark Matter and Dark Energy

Authors: Kemal Gültekin

Abstract: The main aim of this study is to reveal curved space and particle physics effects on the formation of Bose-Einstein condensate (BEC) scalar fields in cosmology and around a black hole. Cosmological scalar fields for dark energy and dark matter may be considered as a result of Bose-Einstein condensation. In this regard, our main attention will be devoted to BECs in curved space. By considering the… ▽ More The main aim of this study is to reveal curved space and particle physics effects on the formation of Bose-Einstein condensate (BEC) scalar fields in cosmology and around a black hole. Cosmological scalar fields for dark energy and dark matter may be considered as a result of Bose-Einstein condensation. In this regard, our main attention will be devoted to BECs in curved space. By considering the dynamics of a BEC scalar field at a microscopic level, we first study the initial phase of the formation of condensation in cosmology. To this end, we initially introduce an effective Minkowski space formulation that enables considering only the effect of particle physics processes, excluding the effect of gravitational particle production and enabling us to see cosmological evolution more easily. Then, by using this formulation, we study a model with a trilinear coupling that induces the processes. After considering the phase evolution of the produced particles, we find that they evolve towards the formation of a BEC if some specific conditions are satisfied. In principle, the effective Minkowski space formulation introduced in this study can be applied to particle physics processes in any sufficiently smooth spacetime. In this regard, we also analyse if a BEC scalar field is realized in the spacetime around a Reissner - Nordström black hole. We find that the produced particles of particle physics processes are localized in a region around the black hole and have a tendency toward condensation if the emerged particles are much heavier than ingoing particles. We also find that such a configuration is phenomenologically viable only if the scalars and the black hole have dark electric charges. Finally, we consider gravitational collapse around Schwarzschild black holes and form a first step towards a study in future about the effects of gravitational collapse on Bose-Einstein condensation. △ Less

Submitted 10 August, 2023; v1 submitted 30 July, 2023; originally announced July 2023.

Comments: Ph.D thesis; 86 pages; 8 figures

The NANOGrav 12.5-year Data Set: Search for Gravitational Wave Memory

Authors: Gabriella Agazie, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Dallas DeGan, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis , et al. (65 additional authors not shown)

Abstract: We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5-yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set (Bayes factor = 2.8). As such, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a sign… ▽ More We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5-yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set (Bayes factor = 2.8). As such, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a signal model that assumes the existence of a common, spatially uncorrelated red noise in addition to a GW memory signal. The median strain upper limit as a function of sky position is approximately $3.3 \times 10^{-14}$. We also find that there are some differences in the upper limits as a function of sky position centered around PSR J0613$-$0200. This suggests that this pulsar has some excess noise which can be confounded with GW memory. Finally, the upper limits as a function of burst epoch continue to improve at later epochs. This improvement is attributable to the continued growth of the pulsar timing array. △ Less

Submitted 25 July, 2023; originally announced July 2023.

Comments: 29 pages, 5 figures

The Impact of Black Hole Scaling Relation Assumptions on the Mass Density of Black Holes

Authors: Cayenne Matt, Kayhan Gültekin, Joseph Simon

Abstract: We examine the effect of supermassive black hole (SMBH) mass scaling relation choice on the inferred SMBH mass population since redshift $z \sim 3$. To make robust predictions for the gravitational wave background (GWB) we must have a solid understanding of the underlying SMBH demographics. Using the SDSS and 3D HST+CANDELS surveys for $0 < z < 3$ we evaluate the inferred SMBH masses from two SMBH… ▽ More We examine the effect of supermassive black hole (SMBH) mass scaling relation choice on the inferred SMBH mass population since redshift $z \sim 3$. To make robust predictions for the gravitational wave background (GWB) we must have a solid understanding of the underlying SMBH demographics. Using the SDSS and 3D HST+CANDELS surveys for $0 < z < 3$ we evaluate the inferred SMBH masses from two SMBH-galaxy scaling relations: $\mathrm{M_{BH}}$-$\mathrm{M_{bulge}}$ and $\mathrm{M_{BH}}$-$σ$. Our SMBH mass functions come directly from stellar mass measurements for $\mathrm{M_{BH}}$-$\mathrm{M_{bulge}}$, and indirectly from stellar mass and galaxy radius measurements along with the galaxy mass fundamental plane for $\mathrm{M_{BH}}$-$σ$. We find that there is a substantial difference in predictions especially for $z > 1$, and this difference increases out to $z = 3$. In particular we find that using velocity dispersion predicts a greater number of SMBHs with masses greater than $10^9 \mathrm{M}_\odot$. The GWB that pulsar timing arrays find evidence for is higher in amplitude than expected from GWB predictions which rely on high redshift extrapolations of local SMBH mass-galaxy scaling relations. The difference in SMBH demographics resulting from different scaling relations may be the origin for the mismatch between the signal amplitude and predictions. Generally, our results suggest that a deeper understanding of the potential redshift evolution of these relations is needed if we are to draw significant insight from their predictions at $z > 1$ △ Less

Submitted 10 July, 2023; originally announced July 2023.

Comments: Accepted by MNRAS

ReveaLLAGN 0: First Look at JWST MIRI data of Sombrero and NGC 1052

Authors: K. Goold, A. Seth, M. Molina, D. Ohlson, J. C. Runnoe, T. Boeker, T. A. Davis, A. Dumont, M. Eracleous, J. A. Fernández-Ontiveros, E. Gallo, A. D. Goulding, J. E. Greene, L. C. Ho, S. B. Markoff, N. Neumayer, R. Plotkin, A. Prieto, S. Satyapal, G. Van De Ven, J. L. Walsh, F. Yuan, A. Feldmeier-Krause, K. Gültekin, S. Hoenig , et al. (6 additional authors not shown)

Abstract: We present the first results from the Revealing Low-Luminosity Active Galactic Nuclei (ReveaLLAGN) survey, a JWST survey of seven nearby LLAGN. We focus on two observations with the Mid-Infrared Instrument's (MIRI) Medium Resolution Spectrograph (MRS) of the nuclei of NGC 1052 and Sombrero (NGC 4594 / M104). We also compare these data to public JWST data of a higher-luminosity AGN, NGC 7319 and NG… ▽ More We present the first results from the Revealing Low-Luminosity Active Galactic Nuclei (ReveaLLAGN) survey, a JWST survey of seven nearby LLAGN. We focus on two observations with the Mid-Infrared Instrument's (MIRI) Medium Resolution Spectrograph (MRS) of the nuclei of NGC 1052 and Sombrero (NGC 4594 / M104). We also compare these data to public JWST data of a higher-luminosity AGN, NGC 7319 and NGC 7469. JWST clearly separates the AGN spectrum from the galaxy light even in Sombrero, the faintest target in our survey; the AGN components have very red spectra. We find that the emission-line widths in both NGC 1052 and Sombrero increase with increasing ionization potential, with FWHM > 1000 km/s for lines with ionization potential > 50 eV. These lines are also significantly blue-shifted in both LLAGN. The high ionization potential lines in NGC 7319 show neither broad widths or significant blue shifts. Many of the lower ionization potential emission lines in Sombrero show significant blue wings extending > 1000 km/s. These features and the emission-line maps in both galaxies are consistent with outflows along the jet direction. Sombrero has the lowest luminosity high-ionization potential lines ([Ne V] and [O IV]) ever measured in the mid-IR, but the relative strengths of these lines are consistent with higher luminosity AGN. On the other hand, the [Ne V] emission is much weaker relative to the [Ne III] and [Ne II] lines of higher-luminosity AGN. These initial results show the great promise that JWST holds for identifying and studying the physical nature of LLAGN. △ Less

Submitted 1 March, 2024; v1 submitted 3 July, 2023; originally announced July 2023.

Comments: Accepted to ApJ Feb 28, 2024

The NANOGrav 15-year Gravitational-Wave Background Analysis Pipeline

Authors: Aaron D. Johnson, Patrick M. Meyers, Paul T. Baker, Neil J. Cornish, Jeffrey S. Hazboun, Tyson B. Littenberg, Joseph D. Romano, Stephen R. Taylor, Michele Vallisneri, Sarah J. Vigeland, Ken D. Olum, Xavier Siemens, Justin A. Ellis, Rutger van Haasteren, Sophie Hourihane, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, J. Andrew Casey-Clyde , et al. (71 additional authors not shown)

Abstract: This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{P… ▽ More This paper presents rigorous tests of pulsar timing array methods and software, examining their consistency across a wide range of injected parameters and signal strength. We discuss updates to the 15-year isotropic gravitational-wave background analyses and their corresponding code representations. Descriptions of the internal structure of the flagship algorithms \texttt{Enterprise} and \texttt{PTMCMCSampler} are given to facilitate understanding of the PTA likelihood structure, how models are built, and what methods are currently used in sampling the high-dimensional PTA parameter space. We introduce a novel version of the PTA likelihood that uses a two-step marginalization procedure that performs much faster when the white noise parameters remain fixed. We perform stringent tests of consistency and correctness of the Bayesian and frequentist analysis software. For the Bayesian analysis, we test prior recovery, injection recovery, and Bayes factors. For the frequentist analysis, we test that the cross-correlation-based optimal statistic, when modified to account for a non-negligible gravitational-wave background, accurately recovers the amplitude of the background. We also summarize recent advances and tests performed on the optimal statistic in the literature from both GWB detection and parameter estimation perspectives. The tests presented here validate current and future analyses of PTA data. △ Less

Submitted 7 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 10 figures, 1 table; Companion paper to "The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background"; For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan DeCesar, Paul B. Demorest, Matthew C. Digman, Timothy Dolch, Brendan Drachler , et al. (74 additional authors not shown)

Abstract: Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav… ▽ More Evidence for a low-frequency stochastic gravitational wave background has recently been reported based on analyses of pulsar timing array data. The most likely source of such a background is a population of supermassive black hole binaries, the loudest of which may be individually detected in these datasets. Here we present the search for individual supermassive black hole binaries in the NANOGrav 15-year dataset. We introduce several new techniques, which enhance the efficiency and modeling accuracy of the analysis. The search uncovered weak evidence for two candidate signals, one with a gravitational-wave frequency of $\sim$4 nHz, and another at $\sim$170 nHz. The significance of the low-frequency candidate was greatly diminished when Hellings-Downs correlations were included in the background model. The high-frequency candidate was discounted due to the lack of a plausible host galaxy, the unlikely astrophysical prior odds of finding such a source, and since most of its support comes from a single pulsar with a commensurate binary period. Finding no compelling evidence for signals from individual binary systems, we place upper limits on the strain amplitude of gravitational waves emitted by such systems. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 23 pages, 13 figures, 2 tables. Accepted for publication in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Search for Anisotropy in the Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore , et al. (68 additional authors not shown)

Abstract: The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these… ▽ More The North American Nanohertz Observatory for Gravitational Waves (NANOGrav) has reported evidence for the presence of an isotropic nanohertz gravitational wave background (GWB) in its 15 yr dataset. However, if the GWB is produced by a population of inspiraling supermassive black hole binary (SMBHB) systems, then the background is predicted to be anisotropic, depending on the distribution of these systems in the local Universe and the statistical properties of the SMBHB population. In this work, we search for anisotropy in the GWB using multiple methods and bases to describe the distribution of the GWB power on the sky. We do not find significant evidence of anisotropy, and place a Bayesian $95\%$ upper limit on the level of broadband anisotropy such that $(C_{l>0} / C_{l=0}) < 20\%$. We also derive conservative estimates on the anisotropy expected from a random distribution of SMBHB systems using astrophysical simulations conditioned on the isotropic GWB inferred in the 15-yr dataset, and show that this dataset has sufficient sensitivity to probe a large fraction of the predicted level of anisotropy. We end by highlighting the opportunities and challenges in searching for anisotropy in pulsar timing array data. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 19 pages, 11 figures; submitted to Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Constraints on Supermassive Black Hole Binaries from the Gravitational Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Paul T. Baker, Bence Bécsy, Laura Blecha, Alexander Bonilla, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler , et al. (89 additional authors not shown)

Abstract: The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary popul… ▽ More The NANOGrav 15-year data set shows evidence for the presence of a low-frequency gravitational-wave background (GWB). While many physical processes can source such low-frequency gravitational waves, here we analyze the signal as coming from a population of supermassive black hole (SMBH) binaries distributed throughout the Universe. We show that astrophysically motivated models of SMBH binary populations are able to reproduce both the amplitude and shape of the observed low-frequency gravitational-wave spectrum. While multiple model variations are able to reproduce the GWB spectrum at our current measurement precision, our results highlight the importance of accurately modeling binary evolution for producing realistic GWB spectra. Additionally, while reasonable parameters are able to reproduce the 15-year observations, the implied GWB amplitude necessitates either a large number of parameters to be at the edges of expected values, or a small number of parameters to be notably different from standard expectations. While we are not yet able to definitively establish the origin of the inferred GWB signal, the consistency of the signal with astrophysical expectations offers a tantalizing prospect for confirming that SMBH binaries are able to form, reach sub-parsec separations, and eventually coalesce. As the significance grows over time, higher-order features of the GWB spectrum will definitively determine the nature of the GWB and allow for novel constraints on SMBH populations. △ Less

Submitted 18 July, 2023; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: Accepted by Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org. Edited to fix two equation typos (Eq.13 & 21), and minor text typos

The NANOGrav 15-year Data Set: Search for Signals from New Physics

Authors: Adeela Afzal, Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Jose Juan Blanco-Pillado, Laura Blecha, Kimberly K. Boddy, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie , et al. (98 additional authors not shown)

Abstract: The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic string… ▽ More The 15-year pulsar timing data set collected by the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) shows positive evidence for the presence of a low-frequency gravitational-wave (GW) background. In this paper, we investigate potential cosmological interpretations of this signal, specifically cosmic inflation, scalar-induced GWs, first-order phase transitions, cosmic strings, and domain walls. We find that, with the exception of stable cosmic strings of field theory origin, all these models can reproduce the observed signal. When compared to the standard interpretation in terms of inspiraling supermassive black hole binaries (SMBHBs), many cosmological models seem to provide a better fit resulting in Bayes factors in the range from 10 to 100. However, these results strongly depend on modeling assumptions about the cosmic SMBHB population and, at this stage, should not be regarded as evidence for new physics. Furthermore, we identify excluded parameter regions where the predicted GW signal from cosmological sources significantly exceeds the NANOGrav signal. These parameter constraints are independent of the origin of the NANOGrav signal and illustrate how pulsar timing data provide a new way to constrain the parameter space of these models. Finally, we search for deterministic signals produced by models of ultralight dark matter (ULDM) and dark matter substructures in the Milky Way. We find no evidence for either of these signals and thus report updated constraints on these models. In the case of ULDM, these constraints outperform torsion balance and atomic clock constraints for ULDM coupled to electrons, muons, or gluons. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 74 pages, 31 figures, 4 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-Year Data Set: Detector Characterization and Noise Budget

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Bécsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. Decesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Elizabeth C. Ferrara, William Fiore, Emmanuel Fonseca , et al. (66 additional authors not shown)

Abstract: Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated… ▽ More Pulsar timing arrays (PTAs) are galactic-scale gravitational wave detectors. Each individual arm, composed of a millisecond pulsar, a radio telescope, and a kiloparsecs-long path, differs in its properties but, in aggregate, can be used to extract low-frequency gravitational wave (GW) signals. We present a noise and sensitivity analysis to accompany the NANOGrav 15-year data release and associated papers, along with an in-depth introduction to PTA noise models. As a first step in our analysis, we characterize each individual pulsar data set with three types of white noise parameters and two red noise parameters. These parameters, along with the timing model and, particularly, a piecewise-constant model for the time-variable dispersion measure, determine the sensitivity curve over the low-frequency GW band we are searching. We tabulate information for all of the pulsars in this data release and present some representative sensitivity curves. We then combine the individual pulsar sensitivities using a signal-to-noise-ratio statistic to calculate the global sensitivity of the PTA to a stochastic background of GWs, obtaining a minimum noise characteristic strain of $7\times 10^{-15}$ at 5 nHz. A power law-integrated analysis shows rough agreement with the amplitudes recovered in NANOGrav's 15-year GW background analysis. While our phenomenological noise model does not model all known physical effects explicitly, it provides an accurate characterization of the noise in the data while preserving sensitivity to multiple classes of GW signals. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 67 pages, 73 figures, 3 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Observations and Timing of 68 Millisecond Pulsars

Authors: Gabriella Agazie, Md Faisal Alam, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Victoria Bonidie, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, Christopher Chapman, Maria Charisi, Shami Chatterjee, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler , et al. (75 additional authors not shown)

Abstract: We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival… ▽ More We present observations and timing analyses of 68 millisecond pulsars (MSPs) comprising the 15-year data set of the North American Nanohertz Observatory for Gravitational Waves (NANOGrav). NANOGrav is a pulsar timing array (PTA) experiment that is sensitive to low-frequency gravitational waves. This is NANOGrav's fifth public data release, including both "narrowband" and "wideband" time-of-arrival (TOA) measurements and corresponding pulsar timing models. We have added 21 MSPs and extended our timing baselines by three years, now spanning nearly 16 years for some of our sources. The data were collected using the Arecibo Observatory, the Green Bank Telescope, and the Very Large Array between frequencies of 327 MHz and 3 GHz, with most sources observed approximately monthly. A number of notable methodological and procedural changes were made compared to our previous data sets. These improve the overall quality of the TOA data set and are part of the transition to new pulsar timing and PTA analysis software packages. For the first time, our data products are accompanied by a full suite of software to reproduce data reduction, analysis, and results. Our timing models include a variety of newly detected astrometric and binary pulsar parameters, including several significant improvements to pulsar mass constraints. We find that the time series of 23 pulsars contain detectable levels of red noise, 10 of which are new measurements. In this data set, we find evidence for a stochastic gravitational-wave background. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 90 pages, 74 figures, 6 tables; published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

The NANOGrav 15-year Data Set: Evidence for a Gravitational-Wave Background

Authors: Gabriella Agazie, Akash Anumarlapudi, Anne M. Archibald, Zaven Arzoumanian, Paul T. Baker, Bence Becsy, Laura Blecha, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Rand Burnette, Robin Case, Maria Charisi, Shami Chatterjee, Katerina Chatziioannou, Belinda D. Cheeseboro, Siyuan Chen, Tyler Cohen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Kathryn Crowter, Curt J. Cutler, Megan E. DeCesar , et al. (89 additional authors not shown)

Abstract: We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectr… ▽ More We report multiple lines of evidence for a stochastic signal that is correlated among 67 pulsars from the 15-year pulsar-timing data set collected by the North American Nanohertz Observatory for Gravitational Waves. The correlations follow the Hellings-Downs pattern expected for a stochastic gravitational-wave background. The presence of such a gravitational-wave background with a power-law-spectrum is favored over a model with only independent pulsar noises with a Bayes factor in excess of $10^{14}$, and this same model is favored over an uncorrelated common power-law-spectrum model with Bayes factors of 200-1000, depending on spectral modeling choices. We have built a statistical background distribution for these latter Bayes factors using a method that removes inter-pulsar correlations from our data set, finding $p = 10^{-3}$ (approx. $3σ$) for the observed Bayes factors in the null no-correlation scenario. A frequentist test statistic built directly as a weighted sum of inter-pulsar correlations yields $p = 5 \times 10^{-5} - 1.9 \times 10^{-4}$ (approx. $3.5 - 4σ$). Assuming a fiducial $f^{-2/3}$ characteristic-strain spectrum, as appropriate for an ensemble of binary supermassive black-hole inspirals, the strain amplitude is $2.4^{+0.7}_{-0.6} \times 10^{-15}$ (median + 90% credible interval) at a reference frequency of 1/(1 yr). The inferred gravitational-wave background amplitude and spectrum are consistent with astrophysical expectations for a signal from a population of supermassive black-hole binaries, although more exotic cosmological and astrophysical sources cannot be excluded. The observation of Hellings-Downs correlations points to the gravitational-wave origin of this signal. △ Less

Submitted 28 June, 2023; originally announced June 2023.

Comments: 30 pages, 18 figures. Published in Astrophysical Journal Letters as part of Focus on NANOGrav's 15-year Data Set and the Gravitational Wave Background. For questions or comments, please email comments@nanograv.org

Revisiting the black hole mass of M87* using VLT/MUSE Adaptive Optics Integral Field Unit data I: Ionized gas kinematics

Authors: J. Osorno, N. Nagar, T. Richtler, P. Humire, K. Gebhardt, K. Gultekin

Abstract: The stellar dynamic-based black hole mass measurements of M87 are twice that determined via ionized gas kinematics; the former is closer to the estimation from the diameter of the gravitationally-lensed ring around the black hole. Using deeper and more comprehensive ionized gas kinematic data, we aim to better constrain the morphology and kinematics of the nuclear ionized gas, thus gaining insight… ▽ More The stellar dynamic-based black hole mass measurements of M87 are twice that determined via ionized gas kinematics; the former is closer to the estimation from the diameter of the gravitationally-lensed ring around the black hole. Using deeper and more comprehensive ionized gas kinematic data, we aim to better constrain the morphology and kinematics of the nuclear ionized gas, thus gaining insights into the reasons behind the disagreement of the measurements. We use both Narrow and Wide Field Mode integral field spectroscopic data from the Multi Unit Spectroscopic Explorer instrument, to model the morphology and kinematics of multiple ionized gas emission lines in the nucleus of M87. The new deep dataset reveals complexities in the nuclear ionized gas kinematics. Several ionized gas filaments can be traced down into the projected sphere of influence. We also found evidence of a partially-filled biconical outflow. The velocity isophotes of the ionized gas disk are twisted and the position angle of the innermost gas disk tends toward a value perpendicular to the radio jet axis. The complexity of the nuclear morphology and kinematics precludes the measurement of an accurate black hole mass. The results support a 6.0 $\times 10^{9}\rm M_{\odot}$ black hole in a 25°disk, rather than a 3.5 $\times 10^{9}\rm M_{\odot}$ black hole in a 42°disk. The specific RIAF model earlier proposed to reconcile the mass measurement discrepancy was also tested. In general, Keplerian disk models perform better than the RIAF model when fitting the sub-arcsec ionized gas disk. A disk inclination close to 25°for the nuclear gas disk, and the warp in the sub-arcsec ionized gas disk, help to reconcile the contradictory nature of the mass discrepancy between stellar and ionized gas black hole masses, and the mis-orientation between the axes of the ionized gas disk and the jet. △ Less

Submitted 31 August, 2023; v1 submitted 21 April, 2023; originally announced April 2023.

Comments: 21 pages, 22 figures (5 of them in the appendix). Accepted in Astronomy & Astrophysics

Journal ref: A&A 679, A37 (2023)

Massive black holes in galactic nuclei: Observations

Authors: Marianne Vestergaard, Kayhan Gültekin

Abstract: Since a black hole does not emit light from its interior, nor does it have a surface on which light from nearby sources can be reflected, observational study of black hole physics requires observing the gravitational impact of the black hole on its surroundings. A massive black hole leaves a dynamical imprint on stars and gas close by. Gas in the immediate vicinity of an accreting massive black ho… ▽ More Since a black hole does not emit light from its interior, nor does it have a surface on which light from nearby sources can be reflected, observational study of black hole physics requires observing the gravitational impact of the black hole on its surroundings. A massive black hole leaves a dynamical imprint on stars and gas close by. Gas in the immediate vicinity of an accreting massive black hole can, due to the presence of the black hole, shine so brightly that it outshines the light of the billions of stars in its host galaxy and be detected across the Universe. By observing the emission from stars and gas and determining their kinematics scientists can extract vital information not only on the fundamental properties of the black holes themselves but also the impact they have on their surroundings. As it turns out, supermassive black holes appear to play a vital role in shaping the Universe as we know it, as they can profoundly impact the star formation history in galaxies. As a consequence, these black holes indirectly impact the cosmic build up of chemical elements heavier than Helium and thus affect when and where life can form. For these reasons alone, observations of massive black holes constitute a very active research area of modern astrophysics. In this chapter we aim to provide a general overview -- fit for a non-expert -- of what scientists have learned, and hope to learn, from analyzing observations of massive black holes and the material around them. △ Less

Submitted 20 April, 2023; originally announced April 2023.

Comments: This chapter is the pre-print of the version currently in production. Please cite this chapter as the following: M. Vestergaard and K. Gültekin. "Massive black holes in galactic nuclei: observations,'' in The Encyclopedia of Cosmology (Set 2): Black Holes, edited by Z. Haiman (World Scientific, New Jersey, 2023). 77 pages, 21 figures

Searching for continuous Gravitational Waves in the second data release of the International Pulsar Timing Array

Authors: M. Falxa, S. Babak, P. T. Baker, B. Bécsy, A. Chalumeau, S. Chen, Z. Chen, N. J. Cornish, L. Guillemot, J. S. Hazboun, C. M. F. Mingarelli, A. Parthasarathy, A. Petiteau, N. S. Pol, A. Sesana, S. B. Spolaor, S. R. Taylor, G. Theureau, M. Vallisneri, S. J. Vigeland, C. A. Witt, X. Zhu, J. Antoniadis, Z. Arzoumanian, M. Bailes , et al. (102 additional authors not shown)

Abstract: The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evi… ▽ More The International Pulsar Timing Array 2nd data release is the combination of datasets from worldwide collaborations. In this study, we search for continuous waves: gravitational wave signals produced by individual supermassive black hole binaries in the local universe. We consider binaries on circular orbits and neglect the evolution of orbital frequency over the observational span. We find no evidence for such signals and set sky averaged 95% upper limits on their amplitude h 95 . The most sensitive frequency is 10nHz with h 95 = 9.1 10-15 . We achieved the best upper limit to date at low and high frequencies of the PTA band thanks to improved effective cadence of observations. In our analysis, we have taken into account the recently discovered common red noise process, which has an impact at low frequencies. We also find that the peculiar noise features present in some pulsars data must be taken into account to reduce the false alarm. We show that using custom noise models is essential in searching for continuous gravitational wave signals and setting the upper limit. △ Less

Submitted 19 March, 2023; originally announced March 2023.

The NANOGrav 12.5-year Data Set: Bayesian Limits on Gravitational Waves from Individual Supermassive Black Hole Binaries

Authors: Zaven Arzoumanian, Paul T. Baker, Laura Blecha, Harsha Blumer, Adam Brazier, Paul R. Brook, Sarah Burke-Spolaor, Bence Bécsy, J. Andrew Casey-Clyde, Maria Charisi, Shami Chatterjee, Siyuan Chen, James M. Cordes, Neil J. Cornish, Fronefield Crawford, H. Thankful Cromartie, Megan E. DeCesar, Paul B. Demorest, Timothy Dolch, Brendan Drachler, Justin A. Ellis, E. C. Ferrara, William Fiore, Emmanuel Fonseca, Gabriel E. Freedman , et al. (53 additional authors not shown)

Abstract: Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data s… ▽ More Pulsar timing array collaborations, such as the North American Nanohertz Observatory for Gravitational Waves (NANOGrav), are seeking to detect nanohertz gravitational waves emitted by supermassive black hole binaries formed in the aftermath of galaxy mergers. We have searched for continuous waves from individual circular supermassive black hole binaries using the NANOGrav's recent 12.5-year data set. We created new methods to accurately model the uncertainties on pulsar distances in our analysis, and we implemented new techniques to account for a common red noise process in pulsar timing array data sets while searching for deterministic gravitational wave signals, including continuous waves. As we found no evidence for continuous waves in our data, we placed 95\% upper limits on the strain amplitude of continuous waves emitted by these sources. At our most sensitive frequency of 7.65 nanohertz, we placed a sky-averaged limit of $h_0 < $ $(6.82 \pm 0.35) \times 10^{-15}$, and $h_0 <$ $(2.66 \pm 0.15) \times 10^{-15}$ in our most sensitive sky location. Finally, we placed a multi-messenger limit of $\mathcal{M} <$ $(1.41 \pm 0.02) \times 10^9 M_\odot$ on the chirp mass of the supermassive black hole binary candidate 3C~66B. △ Less

Submitted 6 June, 2023; v1 submitted 9 January, 2023; originally announced January 2023.

Comments: 22 pages, 12 figures. Accepted by ApJL

Intermediate-mass black holes and the fundamental plane of black hole accretion

Authors: Kayhan Gültekin, Kristina Nyland, Nichole Gray, Greg Fehmer, Tianchi Huang, Matthew Sparkman, Amy E. Reines, Jenny E. Greene, Edward M. Cackett, Vivienne Baldassare

Abstract: We present new 5 GHz VLA observations of a sample of 8 active intermediate-mass black holes with masses $10^{4.9} < M < 10^{6.1}\ M_{\odot}$ found in galaxies with stellar masses $M_{*} < 3 \times 10^{9}\ M_{\odot}$. We detected 5 of the 8 sources at high significance. Of the detections, 4 were consistent with a point source, and one (SDSS J095418.15+471725.1, with black hole mass… ▽ More We present new 5 GHz VLA observations of a sample of 8 active intermediate-mass black holes with masses $10^{4.9} < M < 10^{6.1}\ M_{\odot}$ found in galaxies with stellar masses $M_{*} < 3 \times 10^{9}\ M_{\odot}$. We detected 5 of the 8 sources at high significance. Of the detections, 4 were consistent with a point source, and one (SDSS J095418.15+471725.1, with black hole mass $M < 10^{5}\ M_{\odot}$) clearly shows extended emission that has a jet morphology. Combining our new radio data with the black hole masses and literature X-ray measurements, we put the sources on the fundamental plane of black hole accretion. We find that the extent to which the sources agree with the fundamental plane depends on their star-forming/composite/AGN classification based on optical narrow emission line ratios. he single star-forming source is inconsistent with the fundamental plane. The three composite sources are consistent, and three of the four AGN sources are inconsistent with the fundamental plane. We argue that this inconsistency is genuine and not a result of misattributing star-formation to black hole activity. Instead, we identify the sources in our sample that have AGN-like optical emission line ratios as not following the fundamental plane and thus caution the use of the fundamental plane to estimate masses without additional constraints, such as radio spectral index, radiative efficiency, or the Eddington fraction. △ Less

Submitted 20 September, 2022; originally announced September 2022.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 9 pages, 2 figures. Images can be accessed in fits format from https://doi.org/10.7302/3100-6e62

Gas inflows in the polar ring of NGC 4111: the birth of an AGN

Authors: Gabriel R. Hauschild-Roier, Thaisa Storchi-Bergmann, Richard M. McDermid, Jonelle L. Walsh, Joanne Tan, Jonathan Cohn, Davor Krajnović, Jenny Greene, Monica Valluri, Kayhan Gültekin, Sabine Thater, Glenn van de Ven, Karl Gebhardt, Nora Lützgendorf, Benjamin D. Boizelle, Chung-Pei Ma, Aaron J. Barth

Abstract: We have used Hubble Space Telescope (HST) images, SAURON Integral Field Spectroscopy (IFS) and adaptative optics assisted Gemini NIFS near-infrared K-band IFS to map the stellar and gas distribution, excitation and kinematics of the inner few kpc of the nearby edge-on S0 galaxy NGC 4111. The HST images map its $\approx$ 450 pc diameter dusty polar ring, with an estimated gas mass $\ge10^7$ M… ▽ More We have used Hubble Space Telescope (HST) images, SAURON Integral Field Spectroscopy (IFS) and adaptative optics assisted Gemini NIFS near-infrared K-band IFS to map the stellar and gas distribution, excitation and kinematics of the inner few kpc of the nearby edge-on S0 galaxy NGC 4111. The HST images map its $\approx$ 450 pc diameter dusty polar ring, with an estimated gas mass $\ge10^7$ M$_\odot$. The NIFS datacube maps the inner 110 pc radius at $\approx$ 7 pc spatial resolution revealing a $\approx$ 220 pc diameter polar ring in hot ($2267\pm166$ K) molecular H$_2$ 1-0 S(1) gas embedded in the polar ring. The stellar velocity field shows disk-dominated kinematics along the galaxy plane both in the SAURON large-scale and in the NIFS nuclear-scale data. The large-scale [O III] $\lambda5007$ Åvelocity field shows a superposition of two disk kinematics: one similar to that of the stars and another along the polar ring, showing non-circular motions that seem to connect with the velocity field of the nuclear H$_2$ ring, whose kinematics indicate accelerated inflow to the nucleus. The estimated mass inflow rate is enough not only to feed an Active Galactic Nucleus (AGN) but also to trigger circumnuclear star formation in the near future. We propose a scenario in which gas from the polar ring, which probably originated from the capture of a dwarf galaxy, is moving inwards and triggering an AGN, as supported by the local X-ray emission, which seems to be the source of the H$_2$ 1-0 S(1) excitation. The fact that we see neither near-UV nor Br$γ$ emission suggests that the nascent AGN is still deeply buried under the optically thick dust of the polar ring. △ Less

Submitted 18 July, 2022; v1 submitted 4 March, 2022; originally announced March 2022.

Comments: 18 pages, 21 figures

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 512, Issue 2, May 2022, Pages 2556-2572

Investigating the Accretion Nature of Binary Supermassive Black Hole Candidate SDSS J025214.67-002813.7

Authors: Adi Foord, Xin Liu, Kayhan Gültekin, Kevin Whitley, Fangzheng Shi, Yu-Ching Chen

Abstract: We present results on a multi-wavelength analysis of SDSS J025214.67-002813.7, a system which has been previously classified as a binary AGN candidate based on periodic signals detected in the optical light curves. We use available radio-X-ray observations of the system to investigate the true accretion nature. Analyzing new observations from XMM-Newton and NuSTAR, we characterize the X-ray emissi… ▽ More We present results on a multi-wavelength analysis of SDSS J025214.67-002813.7, a system which has been previously classified as a binary AGN candidate based on periodic signals detected in the optical light curves. We use available radio-X-ray observations of the system to investigate the true accretion nature. Analyzing new observations from XMM-Newton and NuSTAR, we characterize the X-ray emission and search for evidence of circumbinary accretion. Although the 0.5-10 keV spectrum shows evidence of an additional soft emission component, possibly due to extended emission from hot nuclear gas, we find the spectral shape consistent with a single AGN. Compiling a full multi-wavelength SED, we also search for signs of circumbinary accretion, such as a "notch" in the continuum due to the presence of minidisks. We find that the radio-optical emission agrees with the SED of a standard, radio-quiet, AGN, however there is a large deficit in emission blueward of ~1400 A. Although this deficit in emission can plausibly be attributed to a binary AGN system, we find that the SED of SDSS J0252-0028 is better explained by emission from a reddened, single AGN. However, future studies on the expected hard X-ray emission associated with binary AGN (especially in the unequal-mass regime), will allow for more rigorous analyses of the binary AGN hypothesis. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Comments: Submitted to ApJ. 11 pages, 3 figures

H$α$ Reverberation Mapping of the Intermediate-Mass Active Galactic Nucleus in NGC 4395

Authors: Hojin Cho, Jong-Hak Woo, Tommaso Treu, Peter R. Williams, Stephen F. Armen, Aaron J. Barth, Vardha N. Bennert, Wanjin Cho, Alexei V. Filippenko, Elena Gallo, Jaehyuk Geum, Diego González-Buitrago, Kayhan Gültekin, Edmund Hodges-Kluck, John C. Horst, Seong Hyeon Hwang, Wonseok Kang, Minjin Kim, Taewoo Kim, Douglas C. Leonard, Matthew A. Malkan, Raymond P. Remigio, David J. Sand, Jaejin Shin, Donghoon Son , et al. (2 additional authors not shown)

Abstract: We present the results of a high-cadence spectroscopic and imaging monitoring campaign of the active galactic nucleus (AGN) of NGC 4395. High signal-to-noise-ratio spectra were obtained at the Gemini-N 8 m telescope using the GMOS integral field spectrograph (IFS) on 2019 March 7, and at the Keck-I 10 m telescope using the Low-Resolution Imaging Spectrometer (LRIS) with slitmasks on 2019 March 3 a… ▽ More We present the results of a high-cadence spectroscopic and imaging monitoring campaign of the active galactic nucleus (AGN) of NGC 4395. High signal-to-noise-ratio spectra were obtained at the Gemini-N 8 m telescope using the GMOS integral field spectrograph (IFS) on 2019 March 7, and at the Keck-I 10 m telescope using the Low-Resolution Imaging Spectrometer (LRIS) with slitmasks on 2019 March 3 and April 2. Photometric data were obtained with a number of 1 m-class telescopes during the same nights. The narrow-line region (NLR) is spatially resolved; therefore, its variable contributions to the slit spectra make the standard procedure of relative flux calibration impractical. We demonstrate that spatially-resolved data from the IFS can be effectively used to correct the slit-mask spectral light curves. While we obtained no reliable lag owing to the lack of strong variability pattern in the light curves, we constrain the broad line time lag to be less than 3 hr, consistent with the photometric lag of $\sim80$ min reported by Woo et al. (2019). By exploiting the high-quality spectra, we measure the second moment of the broad component of the H$α$ emission line to be $586\pm19$ km s$^{-1}$, superseding the lower value reported by Woo et al. (2019). Combining the revised line dispersion and the photometric time lag, we update the black hole mass as $(1.7\pm 0.3)\times10^4$ M$_{\odot}$. △ Less

Submitted 17 August, 2021; originally announced August 2021.

Comments: Accepted for publication in ApJ (Aug. 16th, 2021). 13 pages, 11 figures

Properties of cold molecular gas in four type-1 active galaxies hosting outflows

Authors: Jessie C. Runnoe, Kayhan Gültekin, David Rupke, Ana López-Sepulcre

Abstract: Feedback from active galactic nuclei (AGN) has proven to be a critical ingredient in the current picture of galaxy assembly and growth. However, observational constraints on AGN-driven outflows face technical challenges and as a result, the cold molecular gas outflow properties of type-1 AGN are not well known. We present new IRAM NOrthern Extended Milimeter Array (NOEMA) observations of CO (1-0)… ▽ More Feedback from active galactic nuclei (AGN) has proven to be a critical ingredient in the current picture of galaxy assembly and growth. However, observational constraints on AGN-driven outflows face technical challenges and as a result, the cold molecular gas outflow properties of type-1 AGN are not well known. We present new IRAM NOrthern Extended Milimeter Array (NOEMA) observations of CO (1-0) in F07599+6508, Z11598$-$0112, F13342+3932, and PG1440+356, all nearby type-1 AGN and ultraluminous infrared galaxies (ULIRGs). We achieve spatial resolution of 1-3 arcsec corresponding to physical scales of 2--8~kpc and spectral resolution of 15-60 km s$^{-1}$ which enables updated CO (1-0) redshifts and a detailed morphological view of the cold molecular gas in these sources. The CO (1-0) luminosities, $L_{CO}^{\prime}$, are in the range 2-12$\times 10^9$ K km s$^{-1}$ pc$^{2}$ and inferred molecular gas masses, $M$(H$_2$), are in the range 2-9$\times 10^9$ M$_{\odot}$. The velocity fields and gas distributions do not unambiguously identify any of these sources as having outflows. However, Z11598$-$0112 has signs of infalling material and after the subtraction of a rotating disk model PG 1440+356 shows complex kinematics in the residuals that may indicate an outflow or warped disk. △ Less

Submitted 27 May, 2021; originally announced May 2021.

Comments: 20 pages, 10 figures, accepted for publication in MNRAS

An ALMA Gas-dynamical Mass Measurement of the Supermassive Black Hole in the Local Compact Galaxy UGC 2698

Authors: Jonathan H. Cohn, Jonelle L. Walsh, Benjamin D. Boizelle, Aaron J. Barth, Karl Gebhardt, Kayhan Gültekin, Akın Yıldırım, David A. Buote, Jeremy Darling, Andrew J. Baker, Luis C. Ho, Kyle M. Kabasares

Abstract: We present 0\farcs{14}-resolution Atacama Large Millimeter/submillimeter Array (ALMA) CO(2$-$1) observations of the circumnuclear gas disk in UGC 2698, a local compact galaxy. The disk exhibits regular rotation with projected velocities rising to 450 km s$^{-1}$ near the galaxy center. We fit gas-dynamical models to the ALMA data cube, assuming the CO emission originates from a dynamically cold, t… ▽ More We present 0\farcs{14}-resolution Atacama Large Millimeter/submillimeter Array (ALMA) CO(2$-$1) observations of the circumnuclear gas disk in UGC 2698, a local compact galaxy. The disk exhibits regular rotation with projected velocities rising to 450 km s$^{-1}$ near the galaxy center. We fit gas-dynamical models to the ALMA data cube, assuming the CO emission originates from a dynamically cold, thin disk, and measured the mass of the supermassive black hole (BH) in UGC 2698 to be $M_{\mathrm{BH}} = (2.46 \pm{0.07}$ [$1σ$ stat] $^{+0.70}_{-0.78}$ [sys])$\times 10^9$ $M_\odot$. UGC 2698 is part of a sample of nearby early-type galaxies that are plausible $z\sim2$ red nugget relics. Previous stellar-dynamical modeling for three galaxies in the sample found BH masses consistent with the BH mass$-$stellar velocity dispersion ($M_{\mathrm{BH}}-σ_\star$) relation but over-massive relative to the BH mass$-$bulge luminosity ($M_{\mathrm{BH}}-L_{\mathrm{bul}}$) correlation, suggesting that BHs may gain the majority of their mass before their host galaxies. However, UGC 2698 is consistent with both $M_{\mathrm{BH}}-σ_\star$ and $M_{\mathrm{BH}}-L_{\mathrm{bul}}$. As UGC 2698 has the largest stellar mass and effective radius in the local compact galaxy sample, it may have undergone more recent mergers that brought it in line with the BH scaling relations. Alternatively, given that the three previously-measured compact galaxies are outliers from $M_{\mathrm{BH}}-L_{\mathrm{bul}}$, while UGC 2698 is not, there may be significant scatter at the poorly sampled high-mass end of the relation. Additional gas-dynamical $M_{\mathrm{BH}}$ measurements for the compact galaxy sample will improve our understanding of BH$-$galaxy co-evolution. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: 18 pages, 7 figures. Submitted

Getting Ready for LISA: The Data, Support and Preparation Needed to Maximize US Participation in Space-Based Gravitational Wave Science

Authors: Kelly Holley-Bockelmann, :, Jillian Bellovary, Peter Bender, Emanuele Berti, Warren Brown, Robert Caldwell, Neil Cornish, Jeremy Darling, Matthew Digman, Mike Eracleous, Kayhan Gultekin, Zoltan Haiman, Kelly Holley-Bockelmann, Joey Key, Shane Larson, Xin Liu, Sean McWilliams, Priyamvada Natarajan, David Shoemaker, Deirdre Shoemaker, Krista Lynne Smith, Marcelle Soares-Santos, Robin Stebbins

Abstract: The NASA LISA Study Team was tasked to study how NASA might support US scientists to participate and maximize the science return from the Laser Interferometer Space Antenna (LISA) mission. LISA is gravitational wave observatory led by ESA with NASA as a junior partner, and is scheduled to launch in 2034. Among our findings: LISA science productivity is greatly enhanced by a full-featured US scienc… ▽ More The NASA LISA Study Team was tasked to study how NASA might support US scientists to participate and maximize the science return from the Laser Interferometer Space Antenna (LISA) mission. LISA is gravitational wave observatory led by ESA with NASA as a junior partner, and is scheduled to launch in 2034. Among our findings: LISA science productivity is greatly enhanced by a full-featured US science center and an open access data model. As other major missions have demonstrated, a science center acts as both a locus and an amplifier of research innovation, data analysis, user support, user training and user interaction. In its most basic function, a US Science Center could facilitate entry into LISA science by hosting a Data Processing Center and a portal for the US community to access LISA data products. However, an enhanced LISA Science Center could: support one of the parallel independent processing pipelines required for data product validation; stimulate the high level of research on data analysis that LISA demands; support users unfamiliar with a novel observatory; facilitate astrophysics and fundamental research; provide an interface into the subtleties of the instrument to validate extraordinary discoveries; train new users; and expand the research community through guest investigator, postdoc and student programs. Establishing a US LISA Science Center well before launch can have a beneficial impact on the participation of the broader astronomical community by providing training, hosting topical workshops, disseminating mock catalogs, software pipelines, and documentation. Past experience indicates that successful science centers are established several years before launch; this early adoption model may be especially relevant for a pioneering mission like LISA. △ Less

Submitted 4 December, 2020; originally announced December 2020.

Comments: 93 pages with a lovely cover page thanks to Bernard Kelly and Elizabeth Ferrara

AGN Triality of Triple Mergers: Multi-wavelength Classifications

Authors: Adi Foord, Kayhan Gultekin, Jessie C. Runnoe, Michael J. Koss

Abstract: We present results from a multi-wavelength analysis searching for multiple AGN systems in nearby (z<0.077) triple galaxy mergers. Combining archival Chandra, SDSS, WISE, and VLA observations, we quantify the rate of nearby triple AGN, as well as investigate possible connections between SMBH accretion and merger environments. Analyzing the multi-wavelength observations of 7 triple galaxy mergers, w… ▽ More We present results from a multi-wavelength analysis searching for multiple AGN systems in nearby (z<0.077) triple galaxy mergers. Combining archival Chandra, SDSS, WISE, and VLA observations, we quantify the rate of nearby triple AGN, as well as investigate possible connections between SMBH accretion and merger environments. Analyzing the multi-wavelength observations of 7 triple galaxy mergers, we find that 1 triple merger has a single AGN (NGC 3341); we discover, for the first time, 4 likely dual AGN (SDSS J1027+1749, SDSS J1631+2352, SDSS J1708+2153, and SDSS J2356-1016); we confirm one triple AGN system, SDSS J0849+1114; and 1 triple merger in our sample remains ambiguous (SDSS J0858+1822). Analyzing the WISE data, we find a trend of increasing N_H (associated with the primary AGN) as a function of increasing W1-W2 color, reflecting that the motions of gas and dust are coupled in merging environments, where large amount of both can be funneled into the active central region during mergers. Additionally, we find that the one triple AGN system in our sample has the highest levels of N_H and W1-W2 color, while the dual AGN candidates all have lower levels; these results are consistent with theoretical merger simulations that suggest higher levels of nuclear gas are more likely to activate AGN in mergers. △ Less

Submitted 1 December, 2020; originally announced December 2020.

Comments: 19 pages, 3 figures, accepted to ApJ

AGN Triality of Triple Mergers: Detection of Faint X-ray Point Sources

Authors: Adi Foord, Kayhan Gultekin, Jessie C. Runnoe, Michael J. Koss

Abstract: We present results from our X-ray analysis of the first systematic search for triple AGN in nearby (z<0.077) triple galaxy mergers. We analyze archival Chandra observations of 7 triple galaxy mergers with BAYMAX (Bayesian Analysis of Multiple AGN in X-rays), fitting each observation with single, dual, and triple X-ray point source models. In doing so, we conclude that 1 triple merger has one X-ray… ▽ More We present results from our X-ray analysis of the first systematic search for triple AGN in nearby (z<0.077) triple galaxy mergers. We analyze archival Chandra observations of 7 triple galaxy mergers with BAYMAX (Bayesian Analysis of Multiple AGN in X-rays), fitting each observation with single, dual, and triple X-ray point source models. In doing so, we conclude that 1 triple merger has one X-ray point source (SDSS J0858+1822, although it's unlikely to be an AGN); 5 triple mergers are likely composed of two X-ray point sources (NGC 3341, SDSS J1027+1749, SDSS J1631+2352, SDSS J1708+2153, and SDSS J2356$-$1016); and one system is composed of three X-ray point sources (SDSS J0849+1114). By fitting the individual X-ray spectra of each point source, we analyze the 2-7 keV luminosities as well as the levels of obscuration associated with each potential AGN. We find that 4/5 dual X-ray point source systems have primary and secondary point sources with bright X-ray luminosities (L_2-7 kev >10^40 erg s^-1), possibly associated with 4 new undetected dual AGN. The dual and triple point source systems are found to have physical separations between 3-9 kpc and flux ratios between 2x10^-3 - 0.84. A multi-wavelength analysis to determine the origin of the X-ray point sources discovered in this work is presented in our companion paper (Foord et al. 2020c). △ Less

Submitted 1 December, 2020; originally announced December 2020.

Comments: 27 pages, 8 figures, accepted to ApJ

Chandra Observations of Abell 2261 Brightest Cluster Galaxy, a Candidate Host to a Recoiling Black Hole

Authors: Kayhan Gultekin, Sarah Burke-Spolaor, Tod R. Lauer, T. Joseph W. Lazio, Leonidas A. Moustakas, Patrick Ogle, Marc Postman

Abstract: We use Chandra X-ray observations to look for evidence of a recoiling black hole from the brightest cluster galaxy in Abell 2261 (A2261-BCG). A2261-BCG is a strong candidate for a recoiling black hole because of its large, flat stellar core, revealed by Hubble Space Telescope imaging observations. We took 100-ksec observations with Chandra and combined it with 35 ksec of archival observations to l… ▽ More We use Chandra X-ray observations to look for evidence of a recoiling black hole from the brightest cluster galaxy in Abell 2261 (A2261-BCG). A2261-BCG is a strong candidate for a recoiling black hole because of its large, flat stellar core, revealed by Hubble Space Telescope imaging observations. We took 100-ksec observations with Chandra and combined it with 35 ksec of archival observations to look for low-level accretion onto a black hole of expected mass $M\sim10^{10}\ M_{\scriptscriptstyle \odot}$ that could possibly be located in one of four off-center stellar knots near the galaxy's center or else in the optical center of the galaxy or in the location of radio emission. We found no X-ray emission arising from a point source in excess of the cluster gas and can place limits on the accretion of any black hole in the central region to a 2-7 keV flux below $4.3 \times 10^{-16}\ \mathrm{erg\ s^{-1}\ cm^{-2}}$, corresponding to a bolometric Eddington fraction of about $10^{-6}$. Thus there is either no $10^{10}\ M_{\scriptscriptstyle \odot}$ black hole in the core of A2261-BCG, or it is accreting at a low level. We also discuss the morphology of the X-ray emitting gas in the cluster and how its asymmetry is consistent with a large dynamic event. △ Less

Submitted 26 October, 2020; originally announced October 2020.

Comments: Accepted to AAS Journals. 12 pages of main text and figures, plus references and appendix tables and figures

Evidence for variability timescale dependent UV/X-ray delay in Seyfert 1 AGN NGC 7469

Authors: Mayukh Pahari, Ian M. McHardy, Federico Vincentelli, Edward M. Cackett, Bradley M. Peterson, Michael R. Goad, Kayhan Gültekin, Keith Horne

Abstract: Using a month-long X-ray lightcurve from ${\it RXTE}$/PCA and 1.5 month long UV continuum lightcurves from ${\it IUE}$ spectra in 1220$-$1970 $Å$, we performed a detailed time-lag study of the Seyfert 1 galaxy NGC 7469. Our cross-correlation analysis confirms previous results showing that the X-rays are delayed relative to the UV continuum at 1315 $Å$ by 3.49 $\pm$ 0.22 days which is possibly caus… ▽ More Using a month-long X-ray lightcurve from ${\it RXTE}$/PCA and 1.5 month long UV continuum lightcurves from ${\it IUE}$ spectra in 1220$-$1970 $Å$, we performed a detailed time-lag study of the Seyfert 1 galaxy NGC 7469. Our cross-correlation analysis confirms previous results showing that the X-rays are delayed relative to the UV continuum at 1315 $Å$ by 3.49 $\pm$ 0.22 days which is possibly caused by either propagating fluctuation or variable comptonisation. However, if variations slower than 5 days are removed from the X-ray lightcurve, the UV variations then lag behind the X-rays variations by 0.37$\pm$0.14 day, consistent with reprocessing of the X-rays by a surrounding accretion disc. A very similar reverberation delay is observed between ${\it Swift}$/XRT X-ray and ${\it Swift}$/UVOT UVW2, U lightcurves. Continuum lightcurves extracted from the ${\it Swift}$/GRISM spectra show delays with respect to X-rays consistent with reverberation. Separating the UV continuum variations faster and slower than 5 days, the slow variations at 1825 $Å$ lag those at 1315 $Å$ by $0.29\pm0.06$ day, while the fast variations are coincident ($0.04\pm0.12$ day). The UV/optical continuum reverberation lag from ${\it IUE}$, ${\it Swift}$ and other optical telescopes at different wavelengths are consistent with the relationship: $τ\propto λ^{4/3}$, predicted for the standard accretion disc theory while the best-fit X-ray delay from ${\it RXTE}$ and ${\it Swift}$/XRT shows a negative X-ray offset of $\sim$0.38 days from the standard disc delay prediction. △ Less

Submitted 16 April, 2020; originally announced April 2020.

Comments: 13 pages, 9 figures, 1 table, accepted for publication in MNRAS

A Second Look at 12 Candidate Dual AGNs using ${\tt BAYMAX}$

Authors: Adi Foord, Kayhan Gültekin, Rebecca Nevin, Julia M. Comerford, Edmund Hodges-Kluck, R. Scott Barrows, Andrew D. Goulding, Jenny E. Greene

Abstract: We present an analysis of 12 optically selected dual AGN candidates at $z < 0.34$. Each candidate was originally identified via double-peaked [O III] $λ$5007 emission lines, and have received follow-up $Chandra$ and $HST$ observations. Because the X-ray data are low-count ($<100$ counts) with small separations ($<1$"), a robust analysis is necessary for classifying each source. Pairing long-slit [… ▽ More We present an analysis of 12 optically selected dual AGN candidates at $z < 0.34$. Each candidate was originally identified via double-peaked [O III] $λ$5007 emission lines, and have received follow-up $Chandra$ and $HST$ observations. Because the X-ray data are low-count ($<100$ counts) with small separations ($<1$"), a robust analysis is necessary for classifying each source. Pairing long-slit [O III] observations with existing $Chandra$ observations, we re-analyze the X-ray observations with ${\tt BAYMAX}$ to determine whether the X-ray emission from each system is more likely a single or dual point source. We find that 4 of the 12 sources are likely dual X-ray point source systems. We examine each point source's spectra via a Monte Carlo method that probabilistically identifies the likely origin of each photon. When doing so, we find that (i) the secondary X-ray point sources in 2 of the systems have $L_{\mathrm{X}}<10^{40}$ erg s$^{-1}$, such that we cannot rule out a non-AGN origin, (ii) one source has a secondary with $L_{\mathrm{X}}>10^{40}$ erg s$^{-1}$ but a spectrum that is too soft to definitively preclude being X-ray emitting diffuse gas that was photoionized by the primary AGN, and (iii) one system (SDSS J1126+2944) is a dual AGN. Additionally, using complementary $HST$ observations, we analyze a sub-sample of systems that are visually identified as merging. Our results suggest that dual AGNs may preferentially reside in mergers with small separations, consistent with both simulations and observations. △ Less

Submitted 3 February, 2020; originally announced February 2020.

Comments: 26 pages, 9 figures, accepted to ApJ

A mechanism for formation of Bose-Einstein condensation in cosmology

Authors: Recai Erdem, Kemal Gültekin

Abstract: We introduce a toy model of scalar particles with a trilinear scalar coupling in cosmology. The trilinear coupling $φ^2χ$ causes production of non-relativistic $φ$ particles through the process $χχ\,\rightarrow\,φφ$ where, initially, only relativistic $χ$ particles are present. We consider the initial times of $χχ\,\rightarrow\,φφ$ and observe that the curved space effects promote formation of Bos… ▽ More We introduce a toy model of scalar particles with a trilinear scalar coupling in cosmology. The trilinear coupling $φ^2χ$ causes production of non-relativistic $φ$ particles through the process $χχ\,\rightarrow\,φφ$ where, initially, only relativistic $χ$ particles are present. We consider the initial times of $χχ\,\rightarrow\,φφ$ and observe that the curved space effects promote formation of Bose-Einstein condensate of $φ$ particles. △ Less

Submitted 16 October, 2019; v1 submitted 23 August, 2019; originally announced August 2019.

Comments: Some technical details in the main text are moved to two appendices. Typos are fixed, and the presentation of the manuscript is improved

A Bayesian Analysis of SDSS J0914+0853, a Low-Mass Dual AGN Candidate

Authors: Adi Foord, Kayhan Gultekin, Mark T. Reynolds, Edmund Hodges-Kluck, Edward M. Cackett, Julia M. Comerford, Ashley L. King, Jon M. Miller, Jessie C. Runnoe

Abstract: We present the first results from BAYMAX (Bayesian AnalYsis of Multiple AGN in X-rays), a tool that uses a Bayesian framework to quantitatively evaluate whether a given Chandra observation is more likely a single or dual point source. Although the most robust method of determining the presence of dual AGNs is to use X-ray observations, only sources that are widely separated relative to the instrum… ▽ More We present the first results from BAYMAX (Bayesian AnalYsis of Multiple AGN in X-rays), a tool that uses a Bayesian framework to quantitatively evaluate whether a given Chandra observation is more likely a single or dual point source. Although the most robust method of determining the presence of dual AGNs is to use X-ray observations, only sources that are widely separated relative to the instrument PSF are easy to identify. It becomes increasingly difficult to distinguish dual AGNs from single AGNs when the separation is on the order of Chandra's angular resolution (<1''). Using likelihood models for single and dual point sources, BAYMAX quantitatively evaluates the likelihood of an AGN for a given source. Specifically, we present results from BAYMAX analyzing the lowest-mass dual AGN candidate to date, SDSS J0914+0853, where archival Chandra data shows a possible secondary AGN ~0.3'' from the primary. Analyzing a new 50 ks Chandra observation, results from BAYMAX shows that SDSS J0914+0853 is most likely a single AGN with a Bayes factor of 13.5 in favor of a single point source model. Further, posterior distributions from the dual point source model are consistent with emission from a single AGN. We find the probability of SDSS J0914+0853 being a dual AGN system with a flux ratio f>0.3 and separation r>0.3'' to be very low. Overall, BAYMAX will be an important tool for correctly classifying candidate dual AGNs in the literature, and studying the dual AGN population where past spatial resolution limits have prevented systematic analyses. △ Less

Submitted 12 April, 2019; originally announced April 2019.

Comments: 16 pages, 5 figures, accepted to ApJ

Astro2020 Science White Paper: Black Holes Across Cosmic Time

Authors: Kayhan Gültekin, Aaron Barth, Karl Gebhardt, Jenny Greene, Luis Ho, Stéphanie Juneau, Chung-Pei Ma, Anil Seth, Vivian U, Monica Valluri, Jonelle Walsh

Abstract: Supermassive black holes are located at the center of most, if not all, massive galaxies. They follow close correlations with global properties of their host galaxies (scaling relations), and are thought to play a crucial role in galaxy evolution. Yet, we lack a complete understanding of fundamental aspects of their growth across cosmic time. In particular, we still do not understand: (1) whether… ▽ More Supermassive black holes are located at the center of most, if not all, massive galaxies. They follow close correlations with global properties of their host galaxies (scaling relations), and are thought to play a crucial role in galaxy evolution. Yet, we lack a complete understanding of fundamental aspects of their growth across cosmic time. In particular, we still do not understand: (1) whether black holes or their host galaxies grow faster and (2) what is the maximum mass that black holes can reach. The high angular resolution capability and sensitivity of 30-m class telescopes will revolutionize our understanding of the extreme end of the black hole and galaxy mass scale. With such facilities, we will be able to dynamically measure masses of the largest black holes and characterize galaxy properties out to redshift $z \sim 1.5$. Together with the evolution of black hole-galaxy scaling relations since $z \sim 1.5$, the maximum mass black hole will shed light on the main channels of black hole growth. △ Less

Submitted 2 April, 2019; originally announced April 2019.

Comments: Astro2020 white paper, 7 pages, 1 figure

Astro2020 Science White Paper: The Local Relics of of Supermassive Black Hole Seeds

Authors: Jenny E Greene, Aaron Barth, Andrea Bellini, Jillian Bellovary, Kelly Holley-Bockelmann, Tuan Do, Elena Gallo, Karl Gebhardt, Kayhan Gultekin, Zoltan Haiman, Matthew Hosek Jr., Dongwon Kim, Mattia Libralato, Jessica Lu, Kristina Nyland, Matthew Malkan, Amy Reines, Anil Seth, Tommaso Treu, Jonelle Walsh, Joan Wrobel

Abstract: We have compelling evidence for stellar-mass black holes (BHs) of ~5-80 M_sun that form through the death of massive stars. We also have compelling evidence for so-called supermassive BHs (10^5-10^10 M_sun) that are predominantly found in the centers of galaxies. We have very good reason to believe there must be BHs with masses in the gap between these ranges: the first ~10^9 M_sun BHs are observe… ▽ More We have compelling evidence for stellar-mass black holes (BHs) of ~5-80 M_sun that form through the death of massive stars. We also have compelling evidence for so-called supermassive BHs (10^5-10^10 M_sun) that are predominantly found in the centers of galaxies. We have very good reason to believe there must be BHs with masses in the gap between these ranges: the first ~10^9 M_sun BHs are observed only hundreds of millions of years after the Big Bang, and all theoretically viable paths to making supermassive BHs require a stage of "intermediate" mass. However, no BHs have yet been reliably detected in the 100-10}^5 M_sun mass range. Uncovering these intermediate-mass BHs of 10^3-10^5 M_sun is within reach in the coming decade. In this white paper we highlight the crucial role that 30-m class telescopes will play in dynamically detecting intermediate-mass black holes, should they exist. △ Less

Submitted 20 March, 2019; originally announced March 2019.

Comments: Science White Paper Submitted for the Astro2020 Decadal Survey on Astronomy and Astrophysics

Astro2020 Science White Paper: Black Hole Growth in Mergers and Dual AGN

Authors: Michael Koss, Vivian U, Edmund Hodges-Kluck, Ezequiel Treister, Laura Blecha, Claudio Ricci, Jeyhan Kartaltepe, Dale Kocevski, Julia M. Comerford, R. Scott Barrows, Claudia Cicone, Francisco Muller-Sanchez, Kayhan Gultekin, Adi Foord, Shobita Satyapal, Jennifer Lotz

Abstract: Hierarchical models of galaxy formation predict that galaxy mergers represent a significant transitional stage of rapid supermassive black hole (SMBH) growth. Yet, the connection between the merging process and enhanced active galactic nuclei (AGN) activity as well as the timescale of SMBH mergers remains highly uncertain. The breakthrough in reconciling the importance of galaxy mergers with black… ▽ More Hierarchical models of galaxy formation predict that galaxy mergers represent a significant transitional stage of rapid supermassive black hole (SMBH) growth. Yet, the connection between the merging process and enhanced active galactic nuclei (AGN) activity as well as the timescale of SMBH mergers remains highly uncertain. The breakthrough in reconciling the importance of galaxy mergers with black hole growth lies in a thoroughly-studied census of dual AGN across cosmic history, which will be enabled by next-generation observational capabilities, theoretical advances, and simulations. This white paper outlines the key questions in galaxy mergers, dual and offset AGN, and proposes multiwavelength solutions using future high-resolution observatories in the X-rays (AXIS, Lynx), near and mid-infrared (30 meter class telescopes, JWST), and submillimeter (ALMA). △ Less

Submitted 15 March, 2019; originally announced March 2019.

Comments: Submission to the 2020-2030 Astronomy and Astrophysics Decadal Survey (Astro2020), 10 pages, 3 figures

Probing the Jet Turnover Frequency Dependence on Black Hole Mass and Mass Accretion Rate

Authors: Erica Hammerstein, Kayhan Gultekin, Ashley King

Abstract: We have examined a sample of 13 sub-Eddington supermassive black holes hosted by galaxies spanning a variety of morphological classifications to further understand the empirical fundamental plane of black hole activity. This plane describes black holes from stellar-mass to supermassive and relates the mass of an accreting black hole and its radio and X-ray luminosities. A key factor in studying th… ▽ More We have examined a sample of 13 sub-Eddington supermassive black holes hosted by galaxies spanning a variety of morphological classifications to further understand the empirical fundamental plane of black hole activity. This plane describes black holes from stellar-mass to supermassive and relates the mass of an accreting black hole and its radio and X-ray luminosities. A key factor in studying the fundamental plane is the turnover frequency, the frequency at which the radio continuum emission becomes optically thin. We measured this turnover frequency using new VLA observations combined, when necessary, with archival Chandra observations. Radio observations are in the range of 5--40 GHz across four frequency bands in B-configuration, giving high spatial resolution to focus on the core emission. We use Markov Chain Monte Carlo methods to fit the continuum emission in order to find the turnover frequency. After testing for correlations, the turnover frequency does not display a significant dependence on either mass or mass accretion rate, indicating that more complicated physics than simple scaling and optical depth effects are at play, as has been suggested by recent theoretical work. △ Less

Submitted 10 March, 2019; originally announced March 2019.

Comments: 10 pages, 4 figures, 2 tables, accepted for publication in ApJ

The Fundamental Plane of Black Hole Accretion and its Use as a Black Hole-Mass Estimator

Authors: Kayhan Gültekin, Ashley L. King, Edward M. Cackett, Kristina Nyland, Jon M. Miller, Tiziana Di Matteo, Sera Markoff, Michael P. Rupen

Abstract: We present an analysis of the fundamental plane of black hole accretion, an empirical correlation of the mass of a black hole ($M$), its 5 GHz radio continuum luminosity ($νL_ν$), and its 2-10 keV X-ray power-law continuum luminosity ($L_X$). We compile a sample of black holes with primary, direct black hole-mass measurements that also have sensitive, high-spatial-resolution radio and X-ray data.… ▽ More We present an analysis of the fundamental plane of black hole accretion, an empirical correlation of the mass of a black hole ($M$), its 5 GHz radio continuum luminosity ($νL_ν$), and its 2-10 keV X-ray power-law continuum luminosity ($L_X$). We compile a sample of black holes with primary, direct black hole-mass measurements that also have sensitive, high-spatial-resolution radio and X-ray data. Taking into account a number of systematic sources of uncertainty and their correlations with the measurements, we use Markov chain Monte Carlo methods to fit a mass-predictor function of the form $\log(M/10^{8}\,M_{\scriptscriptstyle \odot}) = μ_0 + ξ_{μR} \log(L_R / 10^{38}\,\mathrm{erg\,s^{-1}}) + ξ_{μX} \log(L_X / 10^{40}\,\mathrm{erg\,s^{-1}})$. Our best-fit results are $μ_0 = 0.55 \pm 0.22$, $ξ_{μR} = 1.09 \pm 0.10$, and $ξ_{μX} = -0.59^{+0.16}_{-0.15}$ with the natural logarithm of the Gaussian intrinsic scatter in the log-mass direction $\lnε_μ= -0.04^{+0.14}_{-0.13}$. This result is a significant improvement over our earlier mass scaling result because of the increase in active galactic nuclei sample size (from 18 to 30), improvement in our X-ray binary sample selection, better identification of Seyferts, and improvements in our analysis that takes into account systematic uncertainties and correlated uncertainties. Because of these significant improvements, we are able to consider potential influences on our sample by including all sources with compact radio and X-ray emission but ultimately conclude that the fundamental plane can empirically describe all such sources. We end with advice for how to use this as a tool for estimating black hole masses. △ Less

Submitted 8 January, 2019; originally announced January 2019.

Comments: ApJ Accepted. Online interactive version of Figure 7 available at http://kayhan.astro.lsa.umich.edu/supplementary_material/fp/

A Multi-wavelength Analysis of Binary-AGN Candidate PSO J334.2028+01.4075

Authors: Adi Foord, Kayhan Gultekin, Mark Reynolds, Megan Ayers, Tingting Liu, Suvi Gezari, Jessie Runnoe

Abstract: We present analysis of the first Chandra observation of PSO J334.2028+01.4075 (PSO J334), targeted as a binary-AGN candidate based on periodic variations of the optical flux. With no prior targeted X-ray coverage for PSO J334, our new 40 ksec Chandra observation allows for the opportunity to differentiate between a single or binary-AGN system, and if a binary, can characterize the mode of accretio… ▽ More We present analysis of the first Chandra observation of PSO J334.2028+01.4075 (PSO J334), targeted as a binary-AGN candidate based on periodic variations of the optical flux. With no prior targeted X-ray coverage for PSO J334, our new 40 ksec Chandra observation allows for the opportunity to differentiate between a single or binary-AGN system, and if a binary, can characterize the mode of accretion. Simulations show that the two expected accretion disk morphologies for binary-AGN systems are (i) a "cavity", where the inner region of the accretion disk is mostly empty and emission is truncated blueward of the wavelength associated with the temperature of the innermost ring, or (ii) "minidisks", where there is substantial accretion from the cirumbinary disk onto one or both of the members of the binary, each with their own shock-heated thin-disk accretion system. We find the X-ray emission to be well-fit with an absorbed power-law, incompatible with the simple cavity scenario. Further, we construct an SED of PSO J334 by combining radio through X-ray observations and find that the SED agrees well with that of a normal AGN, most likely incompatible with the minidisk scenario. Other analyses, such as locating the quasar on IR color-color diagrams and analyzing the quasar mass predicted by the fundamental plane of black hole activity, further highlight the similarity of PSO J334 with respect to normal AGN. On the multi-wavelength fronts we investigated, we find no evidence supporting PSO J334 as a binary-AGN system, though our analysis remains insensitive to some binary configurations. △ Less

Submitted 15 November, 2017; originally announced November 2017.

Comments: 9 pages, 2 figures, accepted to ApJ

Does the Compact Radio Jet in PG 1700+518 Drive a Molecular Outflow?

Authors: Jessie Runnoe, Kayhan Gültekin, David Rupke

Abstract: Radio jets play an important role in quasar feedback, but direct observations showing how the jets interact with the multi-phase interstellar medium of galaxy disks are few and far between. In this work, we provide new millimeter interferometric observations of PG 1700+518 in order to investigate the effect of its radio jet on the surrounding molecular gas. PG 1700 is a radio-quiet, low-ionization… ▽ More Radio jets play an important role in quasar feedback, but direct observations showing how the jets interact with the multi-phase interstellar medium of galaxy disks are few and far between. In this work, we provide new millimeter interferometric observations of PG 1700+518 in order to investigate the effect of its radio jet on the surrounding molecular gas. PG 1700 is a radio-quiet, low-ionization broad absorption line quasar whose host galaxy has a nearby interacting companion. On sub-kiloparsec scales, the ionized gas is driven to high velocities by a compact radio jet that is identified by radio interferometry. We present observations from the NOrthern Extended Millimeter Array (NOEMA) interferometer with a 3.8 arcsec (16 kpc) synthesized beam where we detect the CO (1-0) emission line at $30σ$ significance with a total flux of $3.12\pm0.02$ Jy km s$^{-1}$ and a typical velocity dispersion of $125\pm5$ km s$^{-1}$. Despite the outflow in ionized gas, we find no concrete evidence that the CO gas is being affected by the radio jet on size scales of a kiloparsec or more. However, a $\sim\!1$ arcsec drift in the spatial centroid of the CO emission as a function of velocity across the emission line and the compact nature of the jet hint that higher spatial resolution observations may reveal a signal of interaction between the jet and molecular gas. △ Less

Submitted 6 November, 2017; v1 submitted 25 October, 2017; originally announced October 2017.

Comments: Updated version with minor changes after last round of review. Accepted to ApJ

A Radio Relic and a Search for the Central Black Hole in the Abell 2261 Brightest Cluster Galaxy

Authors: Sarah Burke-Spolaor, Kayhan Gultekin, Marc Postman, Tod R. Lauer, Joanna M. Taylor, T. Joseph W. Lazio, Leonidas A. Moustakas

Abstract: We present VLA images and HST/STIS spectra of sources within the center of the brightest cluster galaxy (BCG) in Abell 2261. These observations were obtained to test the hypothesis that its extremely large, flat core reflects the ejection of its supermassive black hole. Spectra of three of the four most luminous "knots" embedded in the core were taken to test whether one may represent stars bound… ▽ More We present VLA images and HST/STIS spectra of sources within the center of the brightest cluster galaxy (BCG) in Abell 2261. These observations were obtained to test the hypothesis that its extremely large, flat core reflects the ejection of its supermassive black hole. Spectra of three of the four most luminous "knots" embedded in the core were taken to test whether one may represent stars bound to a displaced massive black hole. The three knots have radial velocity offsets dV < ~150 km/s from the BCG. Knots 2 and 3 show kinematics, colors, and stellar masses consistent with infalling low-mass galaxies or larger stripped cluster members. Large errors in the stellar velocity dispersion of Knot 1, however, mean that we cannot rule out the hypothesis that it hosts a high-mass black hole. A2261-BCG has a compact, relic radio-source offset by 6.5 kpc (projected) from the optical core's center, but no active radio core that would pinpoint the galaxy's central black hole to a tight 10 GHz flux limit <3.6 uJy. Its spectrum and morphology are suggestive of an active galactic nucleus that switched off >48 Myr ago, with an equipartition condition magnetic field of 15 uG. These observations are still consistent with the hypothesis that the nuclear black hole has been ejected from its core, but the critical task of locating the supermassive black hole or demonstrating that A2261-BCG lacks one remains to be done. △ Less

Submitted 21 September, 2017; originally announced September 2017.

Comments: 12 pages, 9 figures, submitted to ApJ

Quasar-Mode Feedback in Nearby Type 1 Quasars: Ubiquitous Kiloparsec-Scale Outflows and Correlations with Black Hole Properties

Authors: David Rupke, Kayhan Gültekin, Sylvain Veilleux

Abstract: The prevalence and properties of kiloparsec-scale outflows in nearby Type 1 quasars have been the subject of little previous attention. This work presents Gemini integral field spectroscopy of ten Type 1 radio-quiet quasars at $z<0.3$. The excellent image quality, coupled with a new technique to remove the point spread function using spectral information, allow the fitting of the underlying host o… ▽ More The prevalence and properties of kiloparsec-scale outflows in nearby Type 1 quasars have been the subject of little previous attention. This work presents Gemini integral field spectroscopy of ten Type 1 radio-quiet quasars at $z<0.3$. The excellent image quality, coupled with a new technique to remove the point spread function using spectral information, allow the fitting of the underlying host on a spaxel-by-spaxel basis. Fits to stars, line-emitting gas, and interstellar absorption show that 100% of the sample host warm ionized and/or cool neutral outflows with spatially-averaged velocities ($\langle v_{98\%}\rangle \equiv \langle v+2σ\rangle$) of 200-1300 km/s and peak velocities (maximum $v_{98\%}$) of 500-2600 km/s. These minor-axis outflows are powered primarily by the central AGN, reach scales of 3-12 kpc, and often fill the field of view. Including molecular data and Type 2 quasar measurements, nearby quasars show a wide range in mass outflow rates ($dM/dt = 1$ to $>$1000 M$_\odot$/yr) and momentum boosts [($c$ $dp/dt$)/L$_\mathrm{AGN}$ = 0.01-20]. After extending the mass scale to Seyferts, $dM/dt$ and $dE/dt$ correlate with black hole mass ($dM/dt \sim M_\mathrm{BH}^{0.7\pm0.3}$ and $dE/dt \sim M_\mathrm{BH}^{1.3\pm0.5}$). Thus, the most massive black holes in the local universe power the most massive and energetic quasar-mode winds. △ Less

Submitted 14 September, 2017; v1 submitted 17 August, 2017; originally announced August 2017.

Comments: 33 pages of text and figures/tables + 37 pages of extended figures; new version to match the paper re-submitted to ApJ

AGN Activity in Nucleated Galaxies as Measured by Chandra

Authors: Adi Foord, Elena Gallo, Edmund Hodges-Kluck, Brendan P. Miller, Vivienne F. Baldassare, Kayhan Gultekin, Oleg Y. Gnedin

Abstract: Motivated by theoretical expectations that Nuclear Star Clusters (NSCs) in galactic centers may provide a favorable environment for super-massive black holes to form and/or efficiently grow, we set out to measure the fraction of nearby nucleated galaxies that also host an Active Galactic Nucleus (AGN). We targeted a distance-limited sample of 98 objects with the Chandra X-ray Telescope, down to a… ▽ More Motivated by theoretical expectations that Nuclear Star Clusters (NSCs) in galactic centers may provide a favorable environment for super-massive black holes to form and/or efficiently grow, we set out to measure the fraction of nearby nucleated galaxies that also host an Active Galactic Nucleus (AGN). We targeted a distance-limited sample of 98 objects with the Chandra X-ray Telescope, down to a uniform X-ray luminosity threshold of $\sim$10$^{38}$ erg s$^{-1}$. The sample is composed of 47 late-types and 51 early-types, enabling us to further investigate the active fraction as a function of galactic morphology. After correcting for contamination to the nuclear X-ray signal from bright X-ray binaries, we measure an active fraction $f$=11.2$\%^{+7.4}_{-4.9}$ (1$σ$ C.L.) across the whole sample, in agreement with previous estimates based on an heterogeneous combination of optical, X-ray and radio diagnostics, by Seth et al. (2008). After accounting for the different stellar mass distributions in our samples, we find no statistically significant difference in the active fraction of early- vs. late-type nucleated galaxies, with $f$=10.6$\%^{+11.9}_{-4.9}$ and 10.8$\%^{+11.3}_{-6.3}$, respectively. For the early-type nucleated galaxies, we are able to carry out a controlled comparison with a parent sample of non-nucleated galaxies covering the same stellar mass range, finding again no statistically significant difference in the active fraction. Taken at face value, our findings suggest that the presence of a NSC does not facilitate nor enhance accretion-powered emission from a nuclear super-massive black hole. This is true even for late-type nucleated galaxies, home to bluer NSCs and arguably larger gas reservoirs. △ Less

Submitted 12 April, 2017; originally announced April 2017.

Comments: ApJ accepted, 12 pages, 8 Figures, 4 Tables