-
Theoretical strong line metallicity diagnostics for the JWST era
Authors:
Prerak Garg,
Desika Narayanan,
Ryan L. Sanders,
Romeel Davè,
Gergö Popping,
Alice E. Shapley,
Daniel P. Stark,
Jonathan R. Trump
Abstract:
The ratios of strong rest-frame optical emission lines are the dominant indicator of metallicities in high-redshift galaxies. Since typical strong-line based metallicity indicators are calibrated on auroral lines at $z=0$, their applicability for galaxies in the distant Universe is unclear. In this paper, we make use of mock emission line data from cosmological simulations to investigate the calib…
▽ More
The ratios of strong rest-frame optical emission lines are the dominant indicator of metallicities in high-redshift galaxies. Since typical strong-line based metallicity indicators are calibrated on auroral lines at $z=0$, their applicability for galaxies in the distant Universe is unclear. In this paper, we make use of mock emission line data from cosmological simulations to investigate the calibration of rest-frame optical emission lines as metallicity indicators at high redshift. Our model, which couples the SIMBA cosmological galaxy formation simulation with cloudy photoionization calculations, includes contributions from HII regions, post-AGB stars and Diffuse Ionized Gas (DIG). We find mild redshift evolution in the 12 indicators that we study, which implies that the dominant physical properties that evolve in our simulations do have a discernible impact on the metallicity calibrations at high redshifts. When comparing our calibrations with high redshift auroral line observations from James Webb Space Telescope we find a slight offset between our model results and the observations and find that a higher ionization parameter at high redshifts can be one of the possible explanations. We explore the physics that drives the shapes of strong-line metallicity relationships and propose calibrations for hitherto unexplored low-metallicity regimes. Finally, we study the contribution of DIG to total line fluxes. We find that the contribution of DIG increases with metallicity at z $\sim$ 0 for singly ionized oxygen and sulfur lines and can be as high as 70% making it crucial to include their contribution when modeling nebular emission.
△ Less
Submitted 8 November, 2023; v1 submitted 12 October, 2023;
originally announced October 2023.
-
The BPT Diagram in Cosmological Galaxy Formation Simulations: Understanding the Physics Driving Offsets at High-Redshift
Authors:
Prerak Garg,
Desika Narayanan,
Nell Byler,
Ryan L. Sanders,
Alice E. Shapley,
Allison L. Strom,
Romeel Davé,
Michaela Hirschmann,
Christopher C. Lovell,
Justin Otter,
Gergö Popping,
George C. Privon
Abstract:
The Baldwin, Philips, & Terlevich diagram of [O III]/H$β$ vs. [N II]/H$α$ (hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies at $z\sim2$ reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understa…
▽ More
The Baldwin, Philips, & Terlevich diagram of [O III]/H$β$ vs. [N II]/H$α$ (hereafter N2-BPT) has long been used as a tool for classifying galaxies based on the dominant source of ionizing radiation. Recent observations have demonstrated that galaxies at $z\sim2$ reside offset from local galaxies in the N2-BPT space. In this paper, we conduct a series of controlled numerical experiments to understand the potential physical processes driving this offset. We model nebular line emission in a large sample of galaxies, taken from the SIMBA cosmological hydrodynamic galaxy formation simulation, using the CLOUDY photoionization code to compute the nebular line luminosities from H II regions. We find that the observed shift toward higher [O III]/H$β$ and [N II]/H$α$ values at high redshift arises from sample selection: when we consider only the most massive galaxies $M_* \sim 10^{10-11} M_\odot$, the offset naturally appears, due to their high metallicities. We predict that deeper observations that probe lower-mass galaxies will reveal galaxies that lie on a locus comparable to $z\sim 0$ observations. Even when accounting for sample selection effects, we find that there is a subtle mismatch between simulations and observations. To resolve this discrepancy, we investigate the impact of varying ionization parameters, H II region densities, gas-phase abundance patterns, and increasing radiation field hardness on N2-BPT diagrams. We find that either decreasing the ionization parameter or increasing the N/O ratio of galaxies at fixed O/H can move galaxies along a self-similar arc in N2-BPT space that is occupied by high-redshift galaxies.
△ Less
Submitted 10 January, 2022;
originally announced January 2022.
-
Powderday: Dust Radiative Transfer for Galaxy Simulations
Authors:
Desika Narayanan,
Matthew J. Turk,
Thomas Robitaille,
Ashley J. Kelly,
B. Connor McClellan,
Ray S. Sharma,
Prerak Garg,
Matthew Abruzzo,
Ena Choi,
Charlie Conroy,
Benjamin D. Johnson,
Benjamin Kimock,
Qi Li,
Christopher C. Lovell,
Sidney Lower,
George C. Privon,
Jonathan Roberts,
Snigdaa Sethuram,
Gregory F. Snyder,
Robert Thompson,
John H. Wise
Abstract:
We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time fle…
▽ More
We present Powderday, a flexible, fast, open-source dust radiative transfer package designed to interface with galaxy formation simulations. Powderday builds on FSPS population synthesis models, Hyperion dust radiative transfer, and employs yt to interface between different software packages. We include our stellar population synthesis modeling on the fly, which allows for significant run-time flexibility in the assumed stellar physics. We include a model for nebular line emission that can employ either precomputed Cloudy lookup tables (for efficiency), or direct photoionization calculations for all young stars (for flexibility). The dust content follows either observationally-motivated prescriptions, direct modeling from galaxy formation simulations, or a novel approach that includes the dust content via learning-based algorithms from the SIMBA cosmological galaxy formation simulation. AGN can additionally be included via a range of prescriptions. The output of these models are broadband SEDs, as well as filter-convolved images. Powderday is designed to eliminate last-mile efforts by researchers that employ different hydrodynamic galaxy formation models, and seamlessly interfaces with GIZMO, AREPO, GASOLINE, CHANGA, and ENZO. We demonstrate the capabilities of the code via three applications: a model for the star formation rate (SFR) - infrared luminosity relation in galaxies (including the impact of AGN); the impact of circumstellar dust around AGB stars on the mid-infrared emission from galaxy SEDs; and the impact of galaxy inclination angle on dust attenuation laws.
△ Less
Submitted 18 June, 2020;
originally announced June 2020.
-
A Hard X-ray Test of HCN Enhancements as a Tracer of Embedded Black Hole Growth
Authors:
G. C. Privon,
C. Ricci,
S. Aalto,
S. Viti,
L. Armus,
T. Díaz-Santos,
E. González-Alfonso,
K. Iwasawa,
D. L. Jeff,
E. Treister,
F. Bauer,
A. S. Evans,
P. Garg,
R. Herrero-Illana,
J. M. Mazzarella,
K. Larson,
L. Blecha,
L. Barcos-Muñoz,
V. Charmandaris,
S. Stierwalt,
M. A. Pérez-Torres
Abstract:
Enhanced emission from the dense gas tracer HCN (relative to HCO$^+ $) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified galaxies with HCN/HCO$ ^+ $ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( $\lesssim15\% $) contribution of an…
▽ More
Enhanced emission from the dense gas tracer HCN (relative to HCO$^+ $) has been proposed as a signature of active galactic nuclei (AGN). In a previous single-dish millimeter line survey we identified galaxies with HCN/HCO$ ^+ $ (1-0) intensity ratios consistent with those of many AGN but whose mid-infrared spectral diagnostics are consistent with little to no ( $\lesssim15\% $) contribution of an AGN to the bolometric luminosity. To search for putative heavily obscured AGN, we present and analyze \nustar hard X-ray (3-79 keV) observations of four such galaxies from the Great Observatories All-sky LIRG Survey. We find no X-ray evidence for AGN in three of the systems and place strong upper limits on the energetic contribution of any heavily obscured ($N_{\rm H}>10^{24}$ cm$^{-2}$) AGN to their bolometric luminosity. The X-ray flux upper limits are presently an order of magnitude below what XDR-driven chemistry model predict are necessary to drive HCN enhancements. In a fourth system we find a hard X-ray excess consistent with the presence of an AGN, but contributing only $\sim3\%$ of the bolometric luminosity. It is also unclear if the AGN is spatially associated with the HCN enhancement. We further explore the relationship between HCN/HCO$^+$ (for several $\mathrm{J}_\mathrm{upper}$ levels) and $L_\mathrm{AGN}/L_\mathrm{IR}$ for a larger sample of systems in the literature. We find no evidence for correlations between the line ratios and the AGN fraction derived from X-rays, indicating that HCN/HCO$^+$ intensity ratios are not driven by the energetic dominance of AGN, nor are they reliable indicators of whether SMBH accretion is ongoing.
△ Less
Submitted 24 April, 2020; v1 submitted 6 April, 2020;
originally announced April 2020.
-
Origin of low surface brightness galaxies: A dynamical study
Authors:
Prerak Garg,
Arunima Banerjee
Abstract:
Low Surface Brighness Galaxies (LSBs), inspite of being gas rich, have low star formation rates and are, therefore, low surface brightness in nature. We calculate Q$_{\rm{RW}}$, the 2-component disc stability parameter as proposed by Romeo \& Wiegert (2011), as a function of galactocentric radius $R$ for a sample of five LSBs, for which mass models, as obtained from HI 21cm radio-sythesis observat…
▽ More
Low Surface Brighness Galaxies (LSBs), inspite of being gas rich, have low star formation rates and are, therefore, low surface brightness in nature. We calculate Q$_{\rm{RW}}$, the 2-component disc stability parameter as proposed by Romeo \& Wiegert (2011), as a function of galactocentric radius $R$ for a sample of five LSBs, for which mass models, as obtained from HI 21cm radio-sythesis observations and R-band photometry, were available in the literature. We find that the median value of Q$_{\rm{RW}}^{\rm{min}}$, the minimum of Q$_{\rm{RW}}$ over $R$, lies between 2.6 and 3.1 for our sample LSBs, which is higher than the median value of 1.8 $\pm$ 0.3 for Q$_{\rm{RW}}^{\rm{min}}$ for a sampleof high surface brightness galaxies (HSBs) as obtained in earlier studies. This clearly shows that LSBs have more stable discs than HSBs, which could explain their low star formation rates and, possibly, their low surface brightness nature. Interestingly, the calculated values of Q$_{\rm{RW}}$ decrease only slightly (median Q$_{\rm{RW}}^{\rm{min}}$ $\sim$ 2.3 - 3) if the discs were taken to respond to the gravitational potential of the dark matter halo only, but reduce by $\sim$ a factor of 2-3 (median Q$_{\rm{RW}}^{\rm{min}}$ $\sim$ 0.7 - 1.5) if they respond to their self-gravity alone. This implies that the dark matter halo is crucial in regulating disc stability in LSBs, which may have important implications for models of galaxy formation and evolution.
△ Less
Submitted 25 July, 2017;
originally announced July 2017.
-
AstroStat - A VO Tool for Statistical Analysis
Authors:
Ajit K. Kembhavi,
Ashish A. Mahabal,
Tejas Kale,
Santosh Jagade,
Ajay Vibhute,
Prerak Garg,
Kaustubh Vaghmare,
Sharmad Navelkar,
Tushar Agrawal,
Deoyani Nandrekar,
Mohasin Shaikh
Abstract:
AstroStat is an easy-to-use tool for performing statistical analysis on data. It has been designed to be compatible with Virtual Observatory (VO) standards thus enabling it to become an integral part of the currently available collection of VO tools. A user can load data in a variety of formats into AstroStat and perform various statistical tests using a menu driven interface. Behind the scenes, a…
▽ More
AstroStat is an easy-to-use tool for performing statistical analysis on data. It has been designed to be compatible with Virtual Observatory (VO) standards thus enabling it to become an integral part of the currently available collection of VO tools. A user can load data in a variety of formats into AstroStat and perform various statistical tests using a menu driven interface. Behind the scenes, all analysis is done using the public domain statistical software - R and the output returned is presented in a neatly formatted form to the user. The analyses performable include exploratory tests, visualizations, distribution fitting, correlation & causation, hypothesis testing, multivariate analysis and clustering. The tool is available in two versions with identical interface and features - as a web service that can be run using any standard browser and as an offline application. AstroStat will provide an easy-to-use interface which can allow for both fetching data and performing power statistical analysis on them.
△ Less
Submitted 10 March, 2015;
originally announced March 2015.