Published October 2021 | Version Accepted Version + Published
Journal Article Open

Which AGN jets quench star formation in massive galaxies?

  • 1. ROR icon Columbia University
  • 2. ROR icon California Institute of Technology
  • 3. ROR icon University of Connecticut
  • 4. ROR icon Northwestern University
  • 5. ROR icon Tel Aviv University
  • 6. ROR icon Harvard-Smithsonian Center for Astrophysics
  • 7. ROR icon Durham University
  • 8. ROR icon University of California, San Diego
  • 9. ROR icon Rutgers, The State University of New Jersey
  • 10. ROR icon University of Zurich

Abstract

Without additional heating, radiative cooling of the halo gas of massive galaxies (Milky Way-mass and above) produces cold gas or stars exceeding that observed. Heating from active galactic nucleus (AGN) jets is likely required, but the jet properties remain unclear. This is particularly challenging for galaxy simulations, where the resolution is orders-of-magnitude insufficient to resolve jet formation and evolution. On such scales, the uncertain parameters include the jet energy form [kinetic, thermal, cosmic ray (CR)]; energy, momentum, and mass flux; magnetic fields; opening angle; precession; and duty cycle. We investigate these parameters in a 10¹⁴M⊙ halo using high-resolution non-cosmological magnetohydrodynamic simulations with the FIRE-2 (Feedback In Realistic Environments) stellar feedback model, conduction, and viscosity. We explore which scenarios qualitatively meet observational constraints on the halo gas and show that CR-dominated jets most efficiently quench the galaxy by providing CR pressure support and modifying the thermal instability. Mildly relativistic (∼MeV or ∼10¹⁰K) thermal plasma jets work but require ∼10 times larger energy input. For fixed energy flux, jets with higher specific energy (longer cooling times) quench more effectively. For this halo mass, kinetic jets are inefficient at quenching unless they have wide opening or precession angles. Magnetic fields also matter less except when the magnetic energy flux reaches ≳ 10⁴⁴ erg s⁻¹ in a kinetic jet model, which significantly widens the jet cocoon. The criteria for a successful jet model are an optimal energy flux and a sufficiently wide jet cocoon with a long enough cooling time at the cooling radius.

Additional Information

© 2021 The Author(s). Published by Oxford University Press on behalf of Royal Astronomical Society. This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model). Accepted 2021 July 10. Received 2021 July 9; in original form 2021 February 3. Published: 19 July 2021. We thank Eliot Quataert for useful discussion. KS acknowledges financial support from the Simons Foundation. We thank Lucy Reading-Ikkanda for preparing the summarizing cartoon Figure (Fig. 7). Support for PFH and co-authors was provided by an Alfred P. Sloan Research Fellowship, NSF Collaborative Research Grant #1715847, and CAREER grant #1455342, and NASA grants NNX15AT06G, JPL 1589742, 17-ATP17-0214. GLB acknowledges financial support from the NSF (grant AST-1615955, OAC-1835509, AST-2006176) and computing support from NSF XSEDE. RSS, CCH, and DAA were supported by the Simons Foundation through the Flatiron Institute. DAA was supported by NSF grant AST-2009687. CAFG was supported by NSF through grants AST-1715216 and CAREER award AST-1652522; by NASA through grant 17-ATP17-0067; and by a Cottrell Scholar Award and a Scialog Award from the Research Corporation for Science Advancement. RF acknowledges financial support from the Swiss National Science Foundation (grant no 194814). DK was supported by NSF grant AST-1715101 and the Cottrell Scholar Award from the Research Corporation for Science Advancement. TKC is supported by Science and Technology Facilities Council (STFC) astronomy consolidated grant ST/T000244/1. Numerical calculations were run on the Flatiron Institute cluster 'popeye' and 'rusty', Caltech cluster 'Wheeler', allocations from XSEDE TG-AST120025, TG-AST130039, and PRAC NSF.1713353 supported by the NSF, and NASA HEC SMD-16-7592. This work was carried out as part of the FIRE project and in collaboration with the SMAUG collaboration. SMAUG gratefully acknowledges support from the Center for Computational Astrophysics at the Flatiron Institute, which is supported by the Simons Foundation. Data Availability Statement: The data supporting the plots within this article are available on reasonable request to the corresponding author. A public version of the gizmo code is available at http://www.tapir.caltech.edu/~phopkins/Site/GIZMO.html.

Attached Files

Published - stab2021.pdf

Accepted Version - 2102.02206.pdf

Files

2102.02206.pdf

Files (20.8 MB)

Name Size
md5:477f6a3207d3232906022db75f88c1e3
8.7 MB Preview Download
md5:ed196322887edbcf78932ff16f5dbfb8
12.1 MB Preview Download

Additional details

Identifiers

Eprint ID
111556
Resolver ID
CaltechAUTHORS:20211020-174301603

Funding

Simons Foundation
Alfred P. Sloan Foundation
NSF
AST-1715847
NSF
AST-1455342
NASA
NNX15AT06G
JPL
1589742
JPL
17-ATP17-0214
NSF
AST-1615955
NSF
OAC-1835509
NSF
AST-2006176
Flatiron Institute
NSF
AST-2009687
NSF
AST-1715216
NSF
AST-1652522
NASA
17-ATP17-0067
Cottrell Scholar of Research Corporation
Swiss National Science Foundation (SNSF)
194814
NSF
AST-1715101
Science and Technology Facilities Council (STFC)
ST/T000244/1
NSF
TG-AST120025
NSF
TG-AST130039
NSF
OAC-1713353
NASA
SMD-16-7592

Dates

Created
2021-10-20
Created from EPrint's datestamp field
Updated
2021-10-21
Created from EPrint's last_modified field

Caltech Custom Metadata

Caltech groups
Astronomy Department , TAPIR