We use high resolution N-body simulations to investigate the Lagrangian bias of cold dark matter haloes within the LCDM cosmology. Our analysis focuses on "proto-haloes", which we identify in the simulation initial conditions with the subsets of particles belonging to individual redshift-zero haloes. We then calculate the number-density and velocity-divergence fields of proto-haloes and estimate their auto spectral densities. We also measure the corresponding cross spectral densities with the linear matter distribution. We use our results to test a Lagrangian-bias model presented by Desjacques and Sheth which is based on the assumption that haloes form out of local density maxima of a specific height. Our comparison validates the predicted functional form for the scale-dependence of the bias for both the density and velocity fields. We also show that the bias coefficients are accurately predicted for the velocity divergence. On the contrary, the theoretical values for the density bias parameters do not accurately match the numerical results as a function of halo mass. This is likely due to the simplistic assumptions that relate virialized haloes to density peaks of a given height in the model. We also detect appreciable stochasticity for the Lagrangian density bias, even on very large scales. These are not included in the model at leading order but correspond to higher order corrections.

We address the issue of the cosmological bias between matter and galaxy distributions, looking at dark-matter haloes as a first step to characterize galaxy clustering. Starting from the linear density field at high redshift, we follow the centre of mass trajectory of the material that will form each halo at late times (proto-halo). We adopt a fluid-like description for the evolution of perturbations in the proto-halo distribution, which is coupled to the matter density field via gravity. We present analytical solutions for the density and velocity fields, in the context of renormalized perturbation theory. We start from the linear solution, then compute one-loop corrections for the propagator and the power spectrum. Finally we analytically resum the propagator and we use a suitable extension of the time-renormalization-group method (Pietroni 2008) to resum the power spectrum. For halo masses M<10^{14} Msol/h our results at z=0 are in good agreement with N-body simulations. Our model is able to predict the halo-matter cross spectrum with an accuracy of 5 per cent up to k = 0.1 h/Mpc approaching the requirements of future galaxy redshift surveys.