2026ApJ...996...37P

We present ionizing photon production efficiencies (ξ_ion) for 63 z = 1.5-6.9 star-forming galaxies using precise nebular dust attenuation corrections from the JWST Assembly of Ultradeep Rest-optical Observations Revealing Astrophysics (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 ξ_ion and redshift, equivalent width of [O III] λ5007, and O32 = [O III] λ5007/[O II] λλ3726, 3729, and negative correlations between ξ_ion and stellar attenuation, UV luminosity (L_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 the relationships between ξ_ion and galaxy property, notably inducing a flat trend between ξ_ion and L_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 the dust systematics on ξ_ion for inclusion in reionization models.