Search

The interaction between supermassive black hole (SMBH) feedback and the circumgalactic medium (CGM) continues to be an open question in galaxy evolution. In our study, we use smoothed particle hydrodynamics simulations to explore the impact of SMBH feedback on galactic metal retention and the motion of metals and gas into and through the CGM of L* galaxies. We examine 140 galaxies from the 25 Mpc cosmological volume Romulus25, with stellar masses between log(M */ M ⊙) = 9.5–11.5. We measure the fraction of metals remaining in the interstellar medium (ISM) and CGM of each galaxy and calculate the expected mass of each SMBH based on the M BH– σ relation (Kormendy & Ho 2013). The deviation of each SMBH from its expected mass, Δ M BH, is compared to the potential of its host via σ. We find that SMBHs with accreted mass above M BH– σ are more effective at removing metals from the ISM than undermassive SMBHs in star-forming galaxies. Overall, overmassive SMBHs suppress the total star formation of their host galaxies and more effectively move metals from the ISM into the CGM. However, we see little to no evacuation of gas from the CGM out of their halos, in contrast with other simulations. Finally, we predict that C iv column densities in the CGM of L* galaxies are unlikely to depend on host galaxy SMBH mass. Our results show that the scatter in the low-mass end of the M BH– σ relation may indicate how effective an SMBH is in the local redistribution of mass in its host galaxy. [ABSTRACT FROM AUTHOR]

We examine the quenching characteristics of 328 isolated dwarf galaxies 10 8 < M star / M ⊙ < 10 10 within the Romulus25 cosmological hydrodynamic simulation. Using mock-observation methods, we identify isolated dwarf galaxies with quenched star formation and make direct comparisons to the quenched fraction in the NASA Sloan Atlas (NSA). Similar to other cosmological simulations, we find a population of quenched, isolated dwarf galaxies below M star < 109 M ⊙ not detected within the NSA. We find that the presence of massive black holes (MBHs) in Romulus25 is largely responsible for the quenched, isolated dwarfs, while isolated dwarfs without an MBH are consistent with quiescent fractions observed in the field. Quenching occurs between z = 0.5–1, during which the available supply of star-forming gas is heated or evacuated by MBH feedback. Mergers or interactions seem to play little to no role in triggering the MBH feedback. At low stellar masses, M star ≲ 109.3 M ⊙, quenching proceeds across several Gyr as the MBH slowly heats up gas in the central regions. At higher stellar masses, M star ≳ 109.3 M ⊙, quenching occurs rapidly within 1 Gyr, with the MBH evacuating gas from the central few kpc of the galaxy and driving it to the outskirts of the halo. Our results indicate the possibility of substantial star formation suppression via MBH feedback within dwarf galaxies in the field. On the other hand, the apparent overquenching of dwarf galaxies due to MBH suggests that higher-resolution and/or better modeling is required for MBHs in dwarfs, and quenched fractions offer the opportunity to constrain current models. [ABSTRACT FROM AUTHOR]

We study satellite counts and quenched fractions for satellites of Milky Way analogs in Romulus25, a large-volume cosmological hydrodynamic simulation. Depending on the definition of a Milky Way analog, we have between 66 and 97 Milky Way analogs in Romulus25, a 25 Mpc per-side uniform volume simulation. We use these analogs to quantify the effect of environment and host properties on satellite populations. We find that the number of satellites hosted by a Milky Way analog increases predominantly with host stellar mass, while environment, as measured by the distance to a Milky Way–mass or larger halo, may have a notable impact in high isolation. Similarly, we find that the satellite quenched fraction for our analogs also increases with host stellar mass, and potentially in higher-density environments. These results are robust for analogs within 3 Mpc of another Milky Way–mass or larger halo, the environmental parameter space where the bulk of our sample resides. We place these results in the context of observations through comparisons to the Exploration of Local VolumE Satellites and Satellites Around Galactic Analogs surveys. Our results are robust to changes in Milky Way analog selection criteria, including those that mimic observations. Finally, as our samples naturally include Milky Way–Andromeda pairs, we examine quenched fractions in pairs versus isolated systems. We find potential evidence, though not conclusive, that pairs, defined as being within 1 Mpc of another Milky Way–mass or larger halo, may have higher satellite quenched fractions. [ABSTRACT FROM AUTHOR]

We explore the characteristics of actively accreting massive black holes (MBHs) within dwarf galaxies in the R omulus25 cosmological hydrodynamic simulation. We examine the MBH occupation fraction, X-ray active fractions, and active galactic nucleus (AGN) scaling relations within dwarf galaxies of stellar mass 108 M ⊙ < M star < 1010 M ⊙ out to redshift z = 2. In the local universe, the MBH occupation fraction is consistent with observed constraints, dropping below unity at M star < 3 × 1010 M ⊙, M 200 < 3 × 1011 M ⊙. Local dwarf AGN in R omulus25 follow observed scaling relations between AGN X-ray luminosity, stellar mass, and star formation rate, though they exhibit slightly higher active fractions and number densities than comparable X-ray observations. Since z = 2, the MBH occupation fraction has decreased, the population of dwarf AGN has become overall less luminous, and as a result the overall number density of dwarf AGN has diminished. We predict the existence of a large population of MBHs in the local universe with low X-ray luminosities and high contamination from X-ray binaries and the hot interstellar medium that are undetectable by current X-ray surveys. These hidden MBHs make up 76% of all MBHs in local dwarf galaxies and include many MBHs that are undermassive relative to their host galaxy’s stellar mass. Their detection relies on not only greater instrument sensitivity but also better modeling of X-ray contaminants or multiwavelength surveys. Our results indicate that dwarf AGN were substantially more active in the past, despite having low luminosity today, and that future deep X-ray surveys may uncover many hidden MBHs in dwarf galaxies out to at least z = 2. [ABSTRACT FROM AUTHOR]

We use hydrodynamical simulations of two Milky Wayâ€"mass galaxies to demonstrate the impact of cosmic-ray pressure on the kinematics of cool and warm circumgalactic gas. Consistent with previous studies, we find that cosmic-ray pressure can dominate over thermal pressure in the inner 50 kpc of the circumgalactic medium (CGM), creating an overall cooler CGM than that of similar galaxy simulations run without cosmic rays. We generate synthetic sight lines of the simulated galaxies’ CGM and use Voigt profile-fitting methods to extract ion column densities, Doppler- b parameters, and velocity centroids of individual absorbers. We directly compare these synthetic spectral line fits with HST/COS CGM absorption-line data analyses, which tend to show that metallic species with a wide range of ionization potential energies are often kinematically aligned. Compared to the Milky Way simulation run without cosmic rays, the presence of cosmic-ray pressure in the inner CGM creates narrower O vi absorption features and broader Si iii absorption features, a quality that is more consistent with observational data. Additionally, because the cool gas is buoyant due to nonthermal cosmic-ray pressure support, the velocity centroids of both cool and warm gas tend to align in the simulated Milky Way with feedback from cosmic rays. Our study demonstrates that detailed, direct comparisons between simulations and observations, focused on gas kinematics, have the potential to reveal the dominant physical mechanisms that shape the CGM. [ABSTRACT FROM AUTHOR]

Using the N-body+Smoothed particle hydrodynamics code, ChaNGa, we identify two merger-driven processes—disk disruption and supermassive black hole (SMBH) feedback—which work together to quench L* galaxies for over 7 Gyr. Specifically, we examine the cessation of star formation in a simulated Milky Way (MW) analog, driven by an interaction with two minor satellites. Both interactions occur within ∼100 Myr of each other, and the satellites both have masses 5–20 times smaller than that of their MW-like host galaxy. Using the genetic modification process of Roth et al., we generate a set of four zoom-in, MW-mass galaxies all of which exhibit unique star formation histories due to small changes to their assembly histories. In two of these four cases, the galaxy is quenched by z = 1. Because these are controlled modifications, we are able to isolate the effects of two closely spaced minor merger events, the relative timing of which determines whether the MW-mass main galaxy quenches. This one–two punch works to: (1) fuel the SMBH at its peak accretion rate and (2) disrupt the cold, gaseous disk of the host galaxy. The end result is that feedback from the SMBH thoroughly and abruptly ends the star formation of the galaxy by z ≈ 1. We search for and find a similar quenching event in R omulus 25, a hydrodynamical (25 Mpc)3 volume simulation, demonstrating that the mechanism is common enough to occur even in a small sample of MW-mass quenched galaxies at z = 0. [ABSTRACT FROM AUTHOR]

We present results from the "Mint" resolution DC Justice League suite of Milky Way–like zoom-in cosmological simulations, which extend our study of nearby galaxies down into the ultrafaint dwarf (UFD) regime for the first time. The mass resolution of these simulations is the highest ever published for cosmological Milky Way zoom-in simulations run to z = 0, with initial star (dark matter) particle masses of 994 (17900) M⊙, and a force resolution of 87 pc. We study the surrounding dwarfs and UFDs, and find that the simulations match the observed dynamical properties of galaxies with −3 > MV > −19, and reproduce the scatter seen in the size–luminosity plane for rh ≳ 200 pc. We predict the vast majority of nearby galaxies will be observable by the Vera Rubin Observatory's coadded Legacy Survey of Space and Time. We additionally show that faint dwarfs with velocity dispersions ≲5 km s−1 result from severe tidal stripping of the host halo. We investigate the quenching of UFDs in a hydrodynamical Milky Way context and find that the majority of UFDs are quenched prior to interactions with the Milky Way, though some of the quenched UFDs retain their gas until infall. Additionally, these simulations yield some unique dwarfs that are the first of their kind to be simulated, e.g., an H i -rich field UFD, a late-forming UFD that has structural properties similar to Crater 2, as well as a compact dwarf satellite that has no dark matter at z = 0. [ABSTRACT FROM AUTHOR]

We revisit the predictions for the merger rate of massive black hole binaries detectable by the Laser Interferometer Space Antenna (LISA) and their background signal for pulsar-timing arrays. We focus on the effect of the delays between the merger of galaxies and the final coalescence of black hole binaries, and on supernova feedback on the black hole growth. By utilizing a semianalytic galaxy formation model, not only do we account for the driving of the evolution of binaries at separations ≲1 pc (gas-driven migration, stellar hardening, and triple/quadruple massive black hole systems), but we also improve on previous studies by accounting for the time spent by black hole pairs from kiloparsec down to parsec separation. We also include the effect of supernova feedback, which may eject gas from the nuclear region of low-mass galaxies, thus hampering the growth of black holes via accretion and suppressing their orbital migration in circumbinary disks. Despite including these novel physical effects, we predict that the LISA detection rate should still be , irrespective of the model for the black hole seeds at high redshifts. Scenarios where black holes form from seeds are more significantly impacted by supernova feedback. We also find that for detectable events, the merging black holes typically have mass ratios between ∼0.1 and 1. Predictions for the stochastic background in the band of pulsar-timing array experiments are instead rather robust and show only a mild dependence on the model. [ABSTRACT FROM AUTHOR]

We investigate the effects of massive black hole growth on the structural evolution of dwarf galaxies within the R omulus 25 cosmological hydrodynamical simulation. We study a sample of 205 central, isolated dwarf galaxies with stellar masses and a central BH. We find that the local MBH–Mstar relation exhibits a high degree of scatter below Mstar < 1010M⊙, which we use to classify BHs as overmassive or undermassive relative to their host Mstar. Within isolated dwarf galaxies, only 8% of undermassive BHs ever undergo a BH merger, while 95% of overmassive BHs grow through a mixture of BH mergers and accretion. We find that isolated dwarf galaxies that host overmassive BHs also follow different evolutionary tracks relative to their undermassive BH counterparts, building up their stars and dark matter earlier and experiencing star formation suppression starting around z = 2. By z = 0.05, overmassive BH hosts above Mstar > 109M⊙ are more likely to exhibit lower central stellar mass density, lower H i gas content, and lower star formation rates than their undermassive BH counterparts. Our results suggest that overmassive BHs in isolated galaxies above Mstar > 109M⊙ are capable of driving feedback, in many cases suppressing and even quenching star formation by late times. [ABSTRACT FROM AUTHOR]