Fast Baryonic Field Painting for Sunyaev-Zel'dovich Analyses: Transfer Function vs. Hybrid Effective Field Theory

We present two approaches for "painting" baryonic properties relevant to the Sunyaev-Zel'dovich (SZ) effect - optical depth and Compton-$y$ - onto 3-dimensional $N$-body simulations, using the MillenniumTNG suite as a benchmark. The goal of these methods is to produce fast and accurate reconstruction methods to aid future analyses of baryonic feedback using the SZ effect. The first approach employs a Gaussian Process emulator to model the SZ quantities via a transfer function, while the second utilizes Hybrid Effective Field Theory (HEFT) to reproduce these quantities within the simulation. Our analysis involves comparing both methods to the true MillenniumTNG optical depth and Compton-$y$ fields using several metrics, including the cross-correlation coefficient, power spectrum, and power spectrum error. Additionally, we assess how well the reconstructed fields correlate with dark matter haloes across various mass thresholds. The results indicate that the transfer function method yields more accurate reconstructions for fields with initially high correlations ($r \approx 1$), such as between the optical depth and dark matter fields. Conversely, the HEFT-based approach proves more effective in enhancing correlations for fields with weaker initial correlations ($r \sim 0.5$), such as between the Compton-$y$ and dark matter fields. Lastly, we discuss extensions of our methods to improve the reconstruction performance at the field level.