skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


This content will become publicly available on October 16, 2026

Title: Emission line tracers of galactic outflows driven by stellar feedback in simulations of isolated disc galaxies
ABSTRACT Galactic outflows strongly influence galactic evolution and have been detected in a range of observations. Hydrodynamic simulations can help interpret these by connecting direct observables to the physical conditions of the outflowing gas. Here we use simulations of isolated disc galaxies ranging from dwarf mass ($$M_{200} = 10^{10}\, \mathrm{M}_{\odot }$$) to Milky Way mass ($$M_{200} = 10^{12}\, \mathrm{M}_{\odot }$$), based on the FIRE-2 subgrid models to investigate multiphase galactic outflows. We use the chimes non-equilibrium chemistry module to create synthetic spectra of common outflow tracers ([C ii]$$_{158\, \mu\rm m}$$, $$\mathrm{CO}_{J(1-0)}$$, H$$\alpha$$ and $$[\mathrm{O}{\small III}]_{5007\, \rm{\mathring{\rm A}}}$$). Using our synthetic spectra we measure the mass outflow rate, kinetic power and momentum flux using observational techniques. In [C ii]$$_{158\, \mu\rm m}$$ we measure outflow rates of $$10^{-4}$$ to 1 $$\mathrm{\, {\rm M}_{\odot }\, \rm yr^{-1}}$$ across an SFR range of $$10^{-3}$$ to 1 $$\text{M}_{\odot }\text{yr}^{-1}$$, which is in reasonable agreement with observations. The significant discrepancy is in $$\mathrm{CO}_{J(1-0)}$$, with the simulations lying $$\approx 1$$ dex below the observational sample. We test observational assumptions used to derive outflow properties from synthetic spectra. We find the greatest uncertainty lies in measurements of electron density, as estimates using the SII doublet can overestimate the actual electron density by up to 2 dex, which changes mass outflow rates by up to 4 dex. We also find that molecular outflows are especially sensitive to the conversion factor between CO luminosity and H2 mass, with outflow rates changing by up to 4 dex in our least massive galaxy. Comparing the outflow properties derived from the synthetic spectra to those derived directly from the simulation, we find that [C ii]$$_{158\, \mu\rm m}$$ probes outflows at greater distances from the disc, whilst we find that molecular gas does not survive at large distances within outflows within our modestly star-forming disc galaxies simulated in this work.  more » « less
Award ID(s):
2108230 2307327
PAR ID:
10655042
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ;
Publisher / Repository:
MNRAS
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
543
Issue:
4
ISSN:
0035-8711
Page Range / eLocation ID:
3428 to 3446
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; (, Monthly Notices of the Royal Astronomical Society)
    null (Ed.)
    ABSTRACT Observations of emission lines in active galactic nuclei (AGNs) often find fast (∼1000 km s−1) outflows extending to kiloparsec scales, seen in ionized, neutral atomic and molecular gas. In this work we present radiative transfer calculations of emission lines in hydrodynamic simulations of AGN outflows driven by a hot wind bubble, including non-equilibrium chemistry, to explore how these lines trace the physical properties of the multiphase outflow. We find that the hot bubble compresses the line-emitting gas, resulting in higher pressures than in the ambient interstellar medium or that would be produced by the AGN radiation pressure. This implies that observed emission line ratios such as [O iv]$$_{25 \, \rm {\mu m}}$$ / [Ne ii]$$_{12 \, \rm {\mu m}}$$, [Ne v]$$_{14 \, \rm {\mu m}}$$ / [Ne ii]$$_{12 \, \rm {\mu m}}$$, and [N iii]$$_{57 \, \rm {\mu m}}$$ / [N ii]$$_{122 \, \rm {\mu m}}$$ constrain the presence of the bubble and hence the outflow driving mechanism. However, the line-emitting gas is under-pressurized compared to the hot bubble itself, and much of the line emission arises from gas that is out of pressure, thermal and/or chemical equilibrium. Our results thus suggest that assuming equilibrium conditions, as commonly done in AGN line emission models, is not justified if a hot wind bubble is present. We also find that ≳50 per cent of the mass outflow rate, momentum flux, and kinetic energy flux of the outflow are traced by lines such as [N ii]$$_{122 \, \rm {\mu m}}$$ and [Ne iii]$$_{15 \, \rm {\mu m}}$$ (produced in the 10$$^{4} \, \rm {K}$$ phase) and [C ii]$$_{158 \, \rm {\mu m}}$$ (produced in the transition from 10$$^{4} \, \rm {K}$$ to 100 K). 
    more » « less
  2. ; ; ; ; ; ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT We present band 6 ALMA observations of a heavily obscured radio-loud (L1.4 GHz = 1025.4 W Hz−1) active galactic nucleus (AGN) candidate at zphot = 6.83 ± 0.06 found in the 1.5 deg2 COSMOS field. The ALMA data reveal detections of exceptionally strong [C ii]158 $$\mu$$m (z[C ii] = 6.8532) and underlying dust continuum emission from this object (COS-87259), where the [C ii] line luminosity, line width, and 158 $$\mu$$m continuum luminosity are comparable to those seen from z ∼ 7 sub-mm galaxies and quasar hosts. The 158 $$\mu$$m continuum detection suggests a total infrared luminosity of $$9\times 10^{12}\, \mathrm{ L}_\odot$$ with corresponding very large obscured star formation rate (1300 M⊙ yr−1) and dust mass ($$2\times 10^9\, \mathrm{ M}_\odot$$). The strong break seen between the VIRCam and IRAC photometry perhaps suggests that COS-87259 is an extremely massive reionization-era galaxy with $$M_\ast \approx 1.7\times 10^{11}\, \mathrm{ M}_\odot$$. Moreover, the MIPS, PACS, and SPIRE detections imply that this object harbours an AGN that is heavily obscured ($$\tau _{_{\mathrm{9.7\,\mu m}}}=2.3$$) with a bolometric luminosity of approximately $$5\times 10^{13}\, \mathrm{ L}_\odot$$. Such a very high AGN luminosity suggests that this object is powered by an ≈1.6 × 10$$^9\, \mathrm{ M}_\odot$$ black hole if accreting near the Eddington limit, and is effectively a highly obscured version of an extremely ultraviolet (UV)-luminous (M1450 ≈ −27.3) z ∼ 7 quasar. Notably, these z ∼ 7 quasars are an exceedingly rare population (∼0.001 deg−2), while COS-87259 was identified over a relatively small field. Future very wide area surveys with e.g. Roman and Euclid have the potential to identify many more extremely red yet UV-bright z ≳ 7 objects similar to COS-87259, providing richer insight into the occurrence of intense obscured star formation and supermassive black hole growth among this population. 
    more » « less
  3. ; ; ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT In recent years, cosmological hydrodynamical simulations have proven their utility as key interpretative tools in the study of galaxy formation and evolution. In this work, we present a comparative analysis of the baryon cycle in three publicly available, leading cosmological simulation suites: EAGLE, IllustrisTNG, and SIMBA. While these simulations broadly agree in terms of their predictions for the stellar mass content and star formation rates of galaxies at $$z\approx 0$$, they achieve this result for markedly different reasons. In EAGLE and SIMBA, we demonstrate that at low halo masses ($$M_{\rm 200c}\lesssim 10^{11.5}\, \mathrm{M}_{\odot }$$), stellar feedback (SF)-driven outflows can reach far beyond the scale of the halo, extending up to $$2\!-\!3\times R_{\rm 200c}$$. In contrast, in TNG, SF-driven outflows, while stronger at the scale of the interstellar medium, recycle within the circumgalactic medium (within $$R_{\rm 200c}$$). We find that active galactic nucleus (AGN)-driven outflows in SIMBA are notably potent, reaching several times $$R_{\rm 200c}$$ even at halo masses up to $$M_{\rm 200c}\approx 10^{13.5}\, \mathrm{M}_{\odot }$$. In both TNG and EAGLE, AGN feedback can eject gas beyond $$R_{\rm 200c}$$ at this mass scale, but seldom beyond $$2\!-\!3\times R_{\rm 200c}$$. We find that the scale of feedback-driven outflows can be directly linked with the prevention of cosmological inflow, as well as the total baryon fraction of haloes within $$R_{\rm 200c}$$. This work lays the foundation to develop targeted observational tests that can discriminate between feedback scenarios, and inform subgrid feedback models in the next generation of simulations. 
    more » « less
  4. ; ; ; ; ; ; ; ; ; ; et al (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT We present the spatially resolved measurements of a cool galactic outflow in the gravitationally lensed galaxy RCS0327 at z ≈ 1.703 using VLT/MUSE IFU observations. We probe the cool outflowing gas, traced by blueshifted Mg ii and Fe ii absorption lines, in 15 distinct regions of the same galaxy in its image-plane. Different physical regions, 5 – 7 kpc apart within the galaxy, drive the outflows at different velocities (Vout ∼ −161 to −240 km s−1), and mass outflow rates ($$\dot{M}_{out} \sim 183$$ – 527 $${\rm M}_{\odot }\, \mathrm{yr}^{-1}$$). The outflow velocities from different regions of the same galaxy vary by 80 km s−1, which is comparable to the variation seen in a large sample of star-burst galaxies in the local universe. Using multiply lensed images of RCS0327, we probe the same star-forming region at different spatial scales (0.5–25 kpc2), we find that outflow velocities vary between ∼ −120 and −242 km s−1, and the mass outflow rates vary between ∼37 and 254 $${\rm M}_{\odot }\, \mathrm{yr}^{-1}$$. The outflow momentum flux in this galaxy is ≥ 100% of the momentum flux provided by star formation in individual regions, and outflow energy flux is ≈ 10% of the total energy flux provided by star formation. These estimates suggest that the outflow in RCS0327 is energy driven. This work shows the importance of small scale variations of outflow properties due to the variations of local stellar properties of the host galaxy in the context of galaxy evolution. 
    more » « less
  5. ; ; ; ; ; ; ; ; (, Monthly Notices of the Royal Astronomical Society)
    ABSTRACT While the first “seeds” of supermassive black holes (BH) can range from $$\sim 10^2-10^6 \rm ~{\rm M}_{\odot }$$, the lowest mass seeds ($$\lesssim 10^3~\rm {\rm M}_{\odot }$$) are inaccessible to most cosmological simulations due to resolution limitations. We present our new BRAHMA simulations that use a novel flexible seeding approach to predict the $$z\ge 7$$ BH populations for low-mass seeds. We ran two types of boxes that model $$\sim 10^3~\rm {\rm M}_{\odot }$$ seeds using two distinct but mutually consistent seeding prescriptions at different simulation resolutions. First, we have the highest resolution $$[9~\mathrm{Mpc}]^3$$ (BRAHMA-9-D3) boxes that directly resolve $$\sim 10^3~\rm {\rm M}_{\odot }$$ seeds and place them within haloes with dense, metal-poor gas. Second, we have lower resolution, larger volume $$[18~\mathrm{Mpc}]^3$$ (BRAHMA-18-E4), and $$\sim [36~\mathrm{Mpc}]^3$$ (BRAHMA-36-E5) boxes that seed their smallest resolvable $$\sim 10^4~\&~10^5~\mathrm{{\rm M}_{\odot }}$$ BH descendants using new stochastic seeding prescriptions calibrated using BRAHMA-9-D3. The three boxes together probe key BH observables between $$\sim 10^3\,\mathrm{ and}\,10^7~\rm {\rm M}_{\odot }$$. The active galactic nuclei (AGN) luminosity function variations are small (factors of $$\sim 2-3$$) at the anticipated detection limits of potential future X-ray facilities ($$\sim 10^{43}~ \mathrm{ergs~s^{-1}}$$ at $$z\sim 7$$). Our simulations predict BHs $$\sim 10-100$$ times heavier than the local $$M_*$$ versus $$M_{\mathrm{ bh}}$$ relations, consistent with several JWST-detected AGN. For different seed models, our simulations merge binaries at $$\sim 1-15~\mathrm{kpc}$$, with rates of $$\sim 200-2000$$ yr−1 for $$\gtrsim 10^3~\rm {\rm M}_{\odot }$$ BHs, $$\sim 6-60$$ yr−1 for $$\gtrsim 10^4~\rm {\rm M}_{\odot }$$ BHs, and up to $$\sim 10$$ yr−1 amongst $$\gtrsim 10^5~\rm {\rm M}_{\odot }$$ BHs. These results suggest that Laser Interferometer Space Antenna mission has promising prospects for constraining seed models. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email