Neural network based emulation of galaxy power spectrum covariances: A reanalysis of BOSS DR12 data

We train neural networks to quickly generate redshift-space galaxy power spectrum covariances from a given parameter set (cosmology and galaxy bias). This covariance emulator utilizes a combination of traditional fully connected network layers and transformer architecture to accurately predict covariance matrices for the high redshift, north galactic cap sample of the BOSS DR12 galaxy catalog. We run simulated likelihood analyses with emulated and brute-force computed covariances, and we quantify the network's performance via two different metrics: (1) difference in ${\ensuremath{\chi}}^{2}$ and (2) likelihood contours for simulated BOSS DR 12 analyses. We find that the emulator returns excellent results over a large parameter range. We then use our emulator to perform a reanalysis of the BOSS HighZ NGC galaxy power spectrum, and find that varying covariance with cosmology along with the model vector produces ${\mathrm{\ensuremath{\Omega}}}_{m}=0.27{6}_{\ensuremath{-}0.015}^{+0.013}$, ${H}_{0}=70.2\ifmmode\pm\else\textpm\fi{}1.9\text{ }\text{ }\mathrm{km}/\mathrm{s}/\mathrm{Mpc}$, and ${\ensuremath{\sigma}}_{8}=0.67{4}_{\ensuremath{-}0.077}^{+0.058}$. These constraints represent an average $0.46\ensuremath{\sigma}$ shift in best-fit values and a 5% increase in constraining power compared to fixing the covariance matrix (${\mathrm{\ensuremath{\Omega}}}_{m}=0.293\ifmmode\pm\else\textpm\fi{}0.017$, ${H}_{0}=70.3\ifmmode\pm\else\textpm\fi{}2.0\text{ }\text{ }\mathrm{km}/\mathrm{s}/\mathrm{Mpc}$, ${\ensuremath{\sigma}}_{8}=0.70{2}_{\ensuremath{-}0.075}^{+0.063}$). This work demonstrates that emulators for more complex cosmological quantities than second-order statistics can be trained over a wide parameter range at sufficiently high accuracy to be implemented in realistic likelihood analyses.