The AURORA Survey: Ionizing Photon Production Efficiency with Minimal Nebular Dust Attenuation Systematics

We present ionizing photon production efficiencies (ξᵢₒₙ) for 63 z=1.5-6.9 star-forming galaxies using precise nebular dust attenuation corrections from the JWST/AURORA survey. A subset of objects within AURORA have individually-determined nebular dust attenuation curves, which vary significantly in shape and normalization, resulting in reduced systematic uncertainty when constraining the total attenuation of H$α$ luminosity, and thus the intrinsic ionizing output within our sample. We find evidence for positive correlations between ξᵢₒₙ and redshift, equivalent width of [OIII]$λ$5007, and O32=[OIII]$λ$5007/[OII]$λ$3726,3729, and negative correlations between ξᵢₒₙ and stellar attenuation, UV luminosity (LUV), stellar mass, and direct-method metallicity. We test alternate dust prescriptions within this sample, and find that the total attenuation is lower when using the commonly-assumed Galactic extinction curve or when assuming that stellar attenuation is equal to nebular attenuation. We also find that assuming either of these alternate dust prescriptions can change the slope of relationships between ξᵢₒₙ and galaxy property, notably inducing a flat trend between ξᵢₒₙ and LUV within AURORA. While the novel nebular dust curves derived from AURORA spectroscopy reveal obscured ionizing photon production within star-forming galaxies at these redshifts, a more complete understanding of stellar attenuation is required to fully reduce dust systematics on ξᵢₒₙ for inclusion in reionization models.