NuSTAR Spectroscopy of Multi-component X-Ray Reflection from NGC 1068

We report on high-energy X-ray observations of the Compton-thick Seyfert 2 galaxy NGC 1068 with NuSTAR, which provide the best constraints to date on its >10 keV spectral shape. The NuSTAR data are consistent with those from past and current instruments to within cross-calibration uncertainties, and we find no strong continuum or line variability over the past two decades, which is in line with its X-ray classification as a reflection-dominated Compton-thick active galactic nucleus. The combined NuSTAR, Chandra, XMM-Newton, and Swift BAT spectral data set offers new insights into the complex secondary emission seen instead of the completely obscured transmitted nuclear continuum. The critical combination of the high signal-to-noise NuSTAR data and the decomposition of the nuclear and extranuclear emission with Chandra allow us to break several model degeneracies and greatly aid physical interpretation. When modeled as a monolithic (i.e., a single N_H) reflector, none of the common Compton reflection models are able to match the neutral fluorescence lines and broad spectral shape of the Compton reflection hump without requiring unrealistic physical parameters (e.g., large Fe overabundances, inconsistent viewing angles, or poor fits to the spatially resolved spectra). A multi-component reflector with three distinct column densities (e.g., with best-fit values of N_H of 1.4 × 10²³, 5.0 × 10²⁴, and 10²⁵ cm^-2) provides a more reasonable fit to the spectral lines and Compton hump, with near-solar Fe abundances. In this model, the higher N_H component provides the bulk of the flux to the Compton hump, while the lower N_H component produces much of the line emission, effectively decoupling two key features of Compton reflection. We find that ≈30% of the neutral Fe Kα line flux arises from >2″ (≈140 pc) and is clearly extended, implying that a significant fraction (and perhaps most) of the <10 keV reflected component arises from regions well outside a parsec-scale torus. These results likely have ramifications for the interpretation of Compton-thick spectra from observations with poorer signal-to-noise and/or more distant objects.