Searching for observational evidence of primordial non-Gaussianty (PNG) is one of the highest priority objectives of contemporary large-scale structure surveys. If detected, PNG would provide crucial evidence for inflation and possibly rule out a wide range of early universe models [1]}, [2]}. The simulations that our emulator was trained on have purely Gaussian initial conditions. Non-Gaussianity can be added to the initial conditions by first drawing the modes of a random Gaussian scalar potential \(\Phi _{\mathrm {g}}(\mathbf {k})\) , and then constructing the non-Gaussian potential [3]},
\(\Phi _{\mathrm {ng}}(\mathbf {k}) = \Phi _{\mathrm {g}}(\mathbf {k}) + f_{\mathrm {NL}} \!\int \! \frac{\mathrm {d}^3k^{\prime }}{(2\pi )^3} K(\mathbf {k}^{\prime }, \mathbf {k} - \mathbf {k}^{\prime }) \Phi _{\mathrm {g}}(\mathbf {k}^{\prime }) \Phi _{\mathrm {g}}(\mathbf {k} - \mathbf {k}^{\prime }) \, .\) 
