Search

We present IR photometry and HST imaging and spectroscopy of Sab galaxy NGC 4826. Schwarzschild dynamical modeling is used to measure its central black hole mass $M$. Photometric decomposition is used to enable a comparison of $M$ to published scaling relations between black hole masses and properties of host bulges. This decomposition implies that NGC 4826 contains classical and pseudo bulges of approximately equal mass. The classical bulge has best-fit S\'ersic index $n=3.27$. The pseudobulge is made up of three parts, an inner lens ($n=0.18$ at $r\lesssim4^{\prime\prime}$), an outer lens ($n=0.17$ at $r \lesssim 45^{\prime\prime}$), and a $n=0.58$ component required to match the surface brightness between the lens components. The total $V$-band luminosity of the galaxy is $M_{VT}=-21.07$, the ratio of classical bulge to total light is $B/T\simeq0.12$, and the ratio of pseudobulge to total light is $PB/T\simeq0.13$. The outer disk is exponential ($n=1.07$) and makes up $D/T=0.75$ of the light of the galaxy. Our best-fit Schwarzschild model has a black hole mass with $1\sigma$ uncertainties of $M=8.4^{+1.7}_{-0.6}\times10^6\ M_\odot$ and a stellar $K$-band mass-to-light ratio of $\Upsilon_K=0.46\pm0.03\ M_{\odot}\ \mathrm{L}_{\odot}^{-1}$ at the assumed distance of 7.27 Mpc. Our modeling is marginally consistent with $M=0$ at the $3\sigma$ limit. These best-fit parameters were calculated assuming the black hole is located where the velocity dispersion is largest; this is offset from the maximum surface brightness, probably because of dust absorption. The black hole mass -- one of the smallest measured by modeling stellar dynamics -- satisfies the well known correlations of $M$ with the $K$-band luminosity, stellar mass, and velocity dispersion of the classical bulge only in contrast to total (classical plus pseudo) bulge luminosity., Comment: Accepted by ApJ. 27 or so pages, 18 figures I think. Online data is available at https://doi.org/10.7302/kr8z-fj98

We present 0$.\!\!^{\prime\prime}{22}$-resolution CO(2$-$1) observations of the circumnuclear gas disk in the local compact galaxy NGC 384 with the Atacama Large Millimeter/submillimeter Array (ALMA). While the majority of the disk displays regular rotation with projected velocities rising to $370$ km s$^{-1}$, the inner $\sim$0\farcs{5} exhibits a kinematic twist. We develop warped disk gas-dynamical models to account for this twist, fit those models to the ALMA data cube, and find a stellar mass-to-light ratio in the $H$-band of \mlabstract\ and a supermassive black hole (BH) mass ($M_{\mathrm{BH}}$) of $M_{\mathrm{BH}}$ $= (7.26^{+0.43}_{-0.48}$ [$1\sigma$ statistical] $^{+0.55}_{-1.00}$ [systematic])$\times 10^8$ $M_\odot$. In contrast to most previous dynamical $M_{\mathrm{BH}}$ measurements in local compact galaxies, which typically found over-massive BHs compared to the local BH mass$-$bulge luminosity and BH mass$-$bulge mass relations, NGC 384 lies within the scatter of those scaling relations. NGC 384 and other local compact galaxies are likely relics of $z\sim2$ red nuggets, and over-massive BHs in these relics indicate BH growth may conclude before the host galaxy stars have finished assembly. Our NGC 384 results may challenge this evolutionary picture, suggesting there may be increased scatter in the scaling relations than previously thought. However, this scatter could be inflated by systematic differences between stellar- and gas-dynamical measurement methods, motivating direct comparisons between the methods for NGC 384 and the other compact galaxies in the sample., Comment: 19 pages, 8 figures, accepted for publication in ApJ

The NANOGrav 15-year data provides compelling evidence for a stochastic gravitational-wave (GW) background at nanohertz frequencies. The simplest model-independent approach to characterizing the frequency spectrum of this signal consists in a simple power-law fit involving two parameters: an amplitude A and a spectral index \gamma. In this paper, we consider the next logical step beyond this minimal spectral model, allowing for a running (i.e., logarithmic frequency dependence) of the spectral index, \gamma_run(f) = \gamma + \beta \ln(f/f_ref). We fit this running-power-law (RPL) model to the NANOGrav 15-year data and perform a Bayesian model comparison with the minimal constant-power-law (CPL) model, which results in a 95% credible interval for the parameter \beta consistent with no running, \beta \in [-0.80,2.96], and an inconclusive Bayes factor, B(RPL vs. CPL) = 0.69 +- 0.01. We thus conclude that, at present, the minimal CPL model still suffices to adequately describe the NANOGrav signal; however, future data sets may well lead to a measurement of nonzero \beta. Finally, we interpret the RPL model as a description of primordial GWs generated during cosmic inflation, which allows us to combine our results with upper limits from big-bang nucleosynthesis, the cosmic microwave background, and LIGO-Virgo-KAGRA., Comment: 17 pages, 4 figures, 2 tables

Pulsar-timing-array experiments have reported evidence for a stochastic background of nanohertz gravitational waves consistent with the signal expected from a population of supermassive--black-hole binaries. Those analyses assume power-law spectra for intrinsic pulsar noise and for the background, as well as a Hellings--Downs cross-correlation pattern among the gravitational-wave--induced residuals across pulsars. These assumptions are idealizations that may not be realized in actuality. We test them in the NANOGrav 15 yr data set using Bayesian posterior predictive checks: after fitting our fiducial model to real data, we generate a population of simulated data-set replications, and use them to assess whether the optimal-statistic significance, inter-pulsar correlations, and spectral coefficients assume extreme values for the real data when compared to the replications. We confirm that the NANOGrav 15 yr data set is consistent with power-law and Hellings--Downs assumptions. We also evaluate the evidence for the stochastic background using posterior-predictive versions of the frequentist optimal statistic and of Bayesian model comparison, and find comparable significance (3.2\ $\sigma$ and 3\ $\sigma$ respectively) to what was previously reported for the standard statistics. We conclude with novel visualizations of the reconstructed gravitational waveforms that enter the residuals for each pulsar. Our analysis strengthens confidence in the identification and characterization of the gravitational-wave background as reported by NANOGrav., Comment: 20 pages, 14 Figures

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 (MUSE) observations in narrow-field mode with adaptive optics and 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 Atacama Large Millimeter/submillimeter Array CO observations. The consistency between these two measurements provides strong support for both the gas dynamical and stellar dynamical methods., Comment: 18 pages, 12 figures, 1 table

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 \sigma - 1.9 \sigma$). The second, at $16 \; \mathrm{nHz}$, is above our GWB realizations with $p = 0.04$ to $0.15$ ($\approx 1.4 \sigma - 2.1 \sigma$). 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., Comment: 10 pages, 8 figures, 1 appendix, submitted to ApJ

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., Comment: 11 pages, 5 figures

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\sigma$ 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\sigma$ 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., Comment: 18 pages, 8 figures, 2 tables. Accepted in ApJ

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 $\eta>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 $\eta>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., Comment: 27 Pages, 10 Figures, 1 Table, Accepted for publication in ApJ

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., Comment: 14 pages, 8 figures, version matching published paper

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., Comment: 29 pages, 5 figures

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}}$-$\sigma$. 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}}$-$\sigma$. 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$, Comment: Accepted by MNRAS

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., Comment: 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

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., Comment: 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 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., Comment: 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 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., Comment: 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 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., Comment: 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

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., Comment: 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

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., Comment: 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

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., Comment: This chapter is the pre-print of the version currently in production. Please cite this chapter as the following: M. Vestergaard and K. G\"ultekin. "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

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., Comment: 22 pages, 12 figures. Accepted by ApJL

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., Comment: 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

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$ \AA 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$\gamma$ emission suggests that the nascent AGN is still deeply buried under the optically thick dust of the polar ring., Comment: 18 pages, 21 figures

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., Comment: Submitted to ApJ. 11 pages, 3 figures