Shuren Zhou, Zhao Chen, Yu Yu

Galaxy clustering is an important probe in the upcoming China Space Station Telescope (CSST) survey to understand the structure growth and reveal the nature of the dark sector. However, it is a long-term challenge to model this biased tracer and connect the observable to the underlying physics. In this work, we present a hybrid Lagrangian bias expansion emulator, combining the Lagrangian bias expansion and the accurate dynamical evolution from $N$-body simulation, to predict the power spectrum of the biased tracer in real space. We employ the Kun simulation suite to construct the emulator, emulating across the space of 8 cosmological parameters including dynamic dark energy $w_0$, $w_a$, and total neutrino mass $\sum m_{\nu}$. The sample variance due to the finite simulation box is further reduced using the Zel'dovich variance control, and it enables the precise measurement of the Lagrangian basis spectra up to quadratic order. The emulation of basis spectra realizes 1% level accuracy, covering wavelength $ k \leq 1 \,{\rm Mpc}^{-1}h$ and redshift $0\leq z\leq 3$ up to quadratic order field. To validate the emulator, we perform the joint fitting of the halo auto power spectrum and the halo-matter cross power spectrum from 46 independent simulations. Depending on the choice of counterpart, the joint fitting is unbiased up to $k_{\rm max}\simeq 0.7\,{\rm Mpc}^{-1}h$ with $1\sim 2$ percent accuracy, for all the redshift and halo mass samples. As one of the CSST cosmological emulator series, our emulator is expected to provide accurate theoretical predictions of the galaxy power spectrum for the upcoming CSST survey.