In this letter we investigate the connection between the cosmic web and the halo distribution through the gravitational potential. We combine three fields of research, cosmic web classification, perturbation theory expansions of the halo bias, and halo (galaxy) mock catalog making methods. In particular, we use the invariants of the tidal field tensor as generating functions (dubbed I-web), to reproduce the halo number counts of a reference catalog from full gravity calculations, populating the dark matter field on a mesh well into the non-linear regime ($\sim5$ Mpc scales). Our results show an unprecedented agreement with the reference power spectrum within 0.5% up to $k=0.72\,h$ Mpc$^{-1}$. By analysing the three point statistics on large scales (configurations of up to $k=0.2\,h$ Mpc$^{-1}$), we find evidence for non-local bias at the 4.8 $\sigma$ confidence level (an information gain of $\sim$ 3.4 $\sigma$ over the commonly used T-web), being fully compatible with the reference catalog. In particular, we find that a detailed description of tidal anisotropic clustering on large scales is crucial to achieve this accuracy. We conclude that the I-web can potentially be useful to study the cosmic web, to improve the generation of mock galaxy catalogs, to improve on halo mass reconstructions, to study primordial non-Gaussianities, to develop new effective Eulerian galaxy bias models at the field level, and to investigate galaxy evolution improving on environmental studies.

Providing accurate predictions for the spatial distribution of matter and luminous tracers in the presence of massive neutrinos is an important task, given the imminent arrival of highly accurate large-scale structure observations. In this work, we address this challenge by extending cosmology-rescaling algorithms to massive neutrino cosmologies. In this way, a $\Lambda$CDM simulation can be modified to provide nonlinear structure formation predictions in the presence a hot component of arbitrary mass, and, if desired, to include non-gravitational modifications to the clustering of matter on large scales. We test the accuracy of the method by comparing its predictions to a suite of simulations carried out explicitly including a neutrino component in its evolution equations. We find that, for neutrino masses in the range $M_\nu \in [0.06, 0.3] ~ \mathrm{eV}$ the matter power spectrum is recovered to better than $1\%$ on all scales $k<2~h~\mathrm{Mpc}^{-1}$. Similarly, the halo mass function is predicted at a few percent level over the range $M_{\rm halo} \in [10^{12}, 10^{15}] ~ h^{-1} ~ \mathrm{M}_{\odot}$, and so do also the multipoles of the galaxy 2-point correlation function in redshift space over $r \in [0.1, 200] ~ h^{-1} ~ \mathrm{Mpc}$. We provide parametric forms for the necessary transformations, as a function of $\Omega_{\rm m}$ and $\Omega_{\nu}$ for various target redshifts.

In this paper we demonstrate that the information encoded in \emph{one} single (sufficiently large) $N$-body simulation can be used to reproduce arbitrary numbers of halo catalogues, using approximated realisations of dark matter density fields with different initial conditions. To this end we use as a reference one realisation (from an ensemble of $300$) of the Minerva $N$-body simulations and the recently published Bias Assignment Method to extract the local and non-local bias linking the halo to the dark matter distribution. We use an approximate (and fast) gravity solver to generate $300$ dark matter density fields from the down-sampled initial conditions of the reference simulation and sample each of these fields using the halo-bias and a kernel, both calibrated from the arbitrarily chosen realisation of the reference simulation. We show that the power spectrum, its variance and the three-point statistics are reproduced within $\sim 2\%$ (up to $k\sim1.0\,h\,{\rm Mpc}^{-1}$), $\sim 5-10\%$ and $\sim 10\%$, respectively. Using a model for the real space power spectrum (with three free bias parameters), we show that the covariance matrices obtained from our procedure lead to parameter uncertainties that are compatible within $\sim 10\%$ with respect to those derived from the reference covariance matrix, and motivate approaches that can help to reduce these differences to $\sim 1\%$. Our method has the potential to learn from one simulation with moderate volumes and high-mass resolution and extrapolate the information of the bias and the kernel to larger volumes, making it ideal for the construction of mock catalogues for present and forthcoming observational campaigns such as Euclid or DESI.