skip to main content


Search for: All records

Creators/Authors contains: "Bellovary, Jillian"

Note: When clicking on a Digital Object Identifier (DOI) number, you will be taken to an external site maintained by the publisher. Some full text articles may not yet be available without a charge during the embargo (administrative interval).
What is a DOI Number?

Some links on this page may take you to non-federal websites. Their policies may differ from this site.

  1. ABSTRACT

    Supermassive black holes (SMBHs) merging in dwarf galaxies will be detectable by the Laser Interferometer Space Antenna (LISA) in the mid-2030s. Previous cosmological hydrodynamic simulations have shown the prediction of massive BHs merging in dwarf galaxies, but these simulations are limited by their resolution and cannot follow BH pairs all the way to coalescence. We calculate the delay time between BH pairing and merger based on the properties of the BHs and their host galaxies, and use these properties to calculate gravitational wave strains for eleven different binary BHs that merge inside dwarf galaxies from eight cosmological simulations. This delay time calculation accounts for dynamical friction due to gas and stars, loss-cone scattering, and hardening of the binary due to gravitational radiation. Out of the eleven BH mergers in the simulations, five BH pairs will merge within 0.8–8 Gyr of forming a close pair and could be observed by LISA, and the remaining six are unresolved due to resolution limitations of the simulation. As all five of the resolved close pairs merge within a Hubble time, we make the broad estimate that close SMBH pairs in dwarf galaxies will merge and be detectable by LISA, but this estimate depends on either the presence of gas during orbital decay or a solution to the dynamical buoyancy problem in cored potentials.

     
    more » « less
  2. ABSTRACT

    Stars and stellar remnants orbiting a supermassive black hole (SMBH) can interact with an active galactic nucleus (AGN) disc. Over time, prograde orbiters (inclination i < 90°) decrease inclination, as well as semimajor axis (a) and eccentricity (e) until orbital alignment with the gas disc (‘disc capture’). Captured stellar-origin black holes (sBH) add to the embedded AGN population that drives sBH–sBH mergers detectable in gravitational waves using LIGO–Virgo–KAGRA or sBH–SMBH mergers detectable with Laser Interferometer Space Antenna. Captured stars can be tidally disrupted by sBH or the SMBH or rapidly grow into massive ‘immortal’ stars. Here, we investigate the behaviour of polar and retrograde orbiters (i ≥ 90°) interacting with the disc. We show that retrograde stars are captured faster than prograde stars, flip to prograde orientation (i < 90°) during capture, and decrease a dramatically towards the SMBH. For sBH, we find a critical angle iret ∼ 113°, below which retrograde sBH decay towards embedded prograde orbits (i → 0°), while for io > iret sBH decay towards embedded retrograde orbits (i → 180°). sBH near polar orbits (i ∼ 90°) and stars on nearly embedded retrograde orbits (i ∼ 180°) show the greatest decreases in a. Whether a star is captured by the disc within an AGN lifetime depends primarily on disc density, and secondarily on stellar type and initial a. For sBH, disc capture time is longest for polar orbits, low-mass sBH, and lower density discs. Larger mass sBH should typically spend more time in AGN discs, with implications for the spin distribution of embedded sBH.

     
    more » « less
  3. ABSTRACT

    We present a phenomenological forward Monte Carlo model for forecasting the population of active galactic nuclei (AGNs) in dwarf galaxies observable via their optical variability. Our model accounts for expected changes in the spectral energy distribution of AGNs in the intermediate-mass black hole (IMBH) mass range and uses observational constraints on optical variability as a function of black hole (BH) mass to generate mock light curves. Adopting several different models for the BH occupation function, including one for off-nuclear IMBHs, we quantify differences in the predicted local AGN mass and luminosity functions in dwarf galaxies. As a result, we are able to model the fraction of variable AGNs as a function of important galaxy host properties, such as host galaxy stellar mass, in the presence of selection effects. We find that our adopted occupation fractions for the ‘heavy’ and ‘light’ initial BH seeding scenarios can be distinguished with variability at the 2–3σ level for galaxy host stellar masses below ∼108M⊙ with data from the upcoming Vera C. Rubin Observatory. We also demonstrate the prevalence of a selection bias whereby recovered IMBH masses fall, on average, above the predicted value from the local host galaxy–BH mass scaling relation with the strength of this bias dependent on the survey sensitivity. Our methodology can be used more broadly to calibrate AGN demographic studies in synoptic surveys. Finally, we show that a targeted ∼ hourly cadence program over a few nights with the Rubin Observatory can provide strong constraints on IMBH masses given their expected rapid variability time-scales.

     
    more » « less
  4. Abstract Current observational facilities have yet to conclusively detect 10 3 –10 4 M ⊙ intermediate-mass black holes (IMBHs) that fill in the evolutionary gap between seed black holes in the early universe and z ∼ 0 supermassive black holes. Dwarf galaxies present an opportunity to reveal active IMBHs amidst persistent star formation. We introduce photoionization simulations tailored to address key physical uncertainties: coincident versus noncoincident mixing of IMBH and starlight excitation, open versus closed geometries of surrounding gas clouds, and different shapes of the spectral energy distribution of active galactic nuclei (AGN). We examine possible AGN emission line diagnostics in the optical and mid-IR, and find that the diagnostics are often degenerate with respect to the investigated physical uncertainties. In spite of these setbacks, and in contrast to recent work, we are able to show that [O iii ]/H β typically remains bright for dwarf AGN powered by IMBHs down to 10 3 M ⊙ . Dwarf AGN are predicted to have inconsistent star-forming and Seyfert/LINER classifications using the most common optical diagnostics. In the mid-IR, [O iv ] 25.9 μ m and [Ar ii ] 6.98 μ m are less sensitive to physical uncertainties than are optical diagnostics. Based on these emission lines, we provide several diagrams of mid-IR emission line diagnostic diagrams with demarcations for separating starbursts and AGN with varying levels of activity. The diagrams are valid over a wide range of ionization parameters and metallicities out to z ∼ 0.1, so will prove useful for future JWST observations of local dwarf AGN in the search for IMBHs. We make our photoionization simulation suite freely available. 
    more » « less
  5. Abstract

    Existing star-forming vs. active galactic nucleus (AGN) classification schemes using optical emission-line diagnostics mostly fail for low-metallicity and/or highly star-forming galaxies, missing AGN in typicalz∼ 0 dwarfs. To recover AGN in dwarfs with strong emission lines (SELs), we present a classification scheme optimizing the use of existing optical diagnostics. We use Sloan Digital Sky Survey emission-line catalogs overlapping the volume- and mass-limited REsolved Spectroscopy Of a Local VolumE (RESOLVE) and Environmental COntex (ECO) surveys to determine the AGN percentage in SEL dwarfs. Our photoionization grids show that the [Oiii]/Hβversus [Sii]/Hαdiagram (Siiplot) and [Oiii]/Hβversus [Oi]/Hαdiagram (Oiplot) are less metallicity sensitive and more successful in identifying dwarf AGN than the popular [Oiii]/Hβversus [Nii]/Hαdiagnostic (Niiplot or “BPT diagram”). We identify a new category of “star-forming AGN” (SF-AGN) classified as star-forming by the Niiplot but as AGN by the Siiand/or Oiplots. Including SF-AGN, we find thez∼ 0 AGN percentage in dwarfs with SELs to be ∼3%–16%, far exceeding most previous optical estimates (∼1%). The large range in our dwarf AGN percentage reflects differences in spectral fitting methodologies between catalogs. The highly complete nature of RESOLVE and ECO allows us to normalize strong emission-line galaxy statistics to the full galaxy population, reducing the dwarf AGN percentage to ∼0.6%–3.0%. The newly identified SF-AGN are mostly gas-rich dwarfs with halo mass <1011.5M, where highly efficient cosmic gas accretion is expected. Almost all SF-AGN also have low metallicities (Z≲ 0.4Z), demonstrating the advantage of our method.

     
    more » « less
  6. null (Ed.)
    ABSTRACT Massive black holes often exist within dwarf galaxies, and both simulations and observations have shown that a substantial fraction of these may be off-centre with respect to their hosts. We trace the evolution of off-centre massive black holes (MBHs) in dwarf galaxies using cosmological hydrodynamical simulations, and show that the reason for off-centre locations is mainly due to galaxy–galaxy mergers. We calculate dynamical time-scales and show that off-centre MBHs are unlikely to sink to their galaxys’ centres within a Hubble time, due to the shape of the hosts’ potential wells and low stellar densities. These wandering MBHs are unlikely to be detected electromagnetically, nor is there a measurable dynamical effect on the galaxy’s stellar population. We conclude that off-centre MBHs may be common in dwarfs, especially if the mass of the MBH is small or the stellar mass of the host galaxy is large. However, detecting them is extremely challenging, because their accretion luminosities are very low and they do not measurably alter the dynamics of their host galaxies. 
    more » « less
  7. null (Ed.)
    ABSTRACT Active galactic nuclei (AGN) are powered by the accretion of discs of gas on to supermassive black holes (SMBHs). Stars and stellar remnants orbiting the SMBH in the nuclear star cluster (NSC) will interact with the AGN disc. Orbiters plunging through the disc experience a drag force and, through repeated passage, can have their orbits captured by the disc. A population of embedded objects in AGN discs may be a significant source of binary black hole mergers, supernovae, tidal disruption events, and embedded gamma-ray bursts. For two representative AGN disc models, we use geometric drag and Bondi–Hoyle–Littleton drag to determine the time to capture for stars and stellar remnants. We assume a range of initial inclination angles and semimajor axes for circular Keplerian prograde orbiters. Capture time strongly depends on the density and aspect ratio of the chosen disc model, the relative velocity of the stellar object with respect to the disc, and the AGN lifetime. We expect that for an AGN disc density $\rho \gtrsim 10^{-11}{\rm g\, cm^{-3}}$ and disc lifetime ≥1 Myr, there is a significant population of embedded stellar objects, which can fuel mergers detectable in gravitational waves with LIGO-Virgo and LISA. 
    more » « less
  8. null (Ed.)
  9. Abstract The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA’s first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultra-compact stellar-mass binaries, massive black hole binaries, and extreme or interme-diate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe. 
    more » « less
    Free, publicly-accessible full text available December 1, 2024