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

We present ionizing photon production efficiencies (${\xi}_{\rm ion}$) 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${\alpha}$ luminosity, and thus the intrinsic ionizing output within our sample. We find evidence for positive correlations between ${\xi}_{\rm ion}$ and redshift, equivalent width of [OIII]${\lambda}$5007, and O32=[OIII]${\lambda}$5007/[OII]${\lambda}$3726,3729, and negative correlations between ${\xi}_{\rm ion}$ and stellar attenuation, UV luminosity (L$_{\rm UV}$), 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 ${\xi}_{\rm ion}$ and galaxy property, notably inducing a flat trend between ${\xi}_{\rm ion}$ and L$_{\rm UV}$ 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 ${\xi}_{\rm ion}$ for inclusion in reionization models.