Galaxy Assembly and Evolution in the P-Millennium simulation: galaxy clustering

[abridged] We present results from the latest version of the GAEA theoretical model of galaxy formation coupled with merger trees extracted from the Planck Millennium Simulation (PMS). With respect to the Millennium Simulation, the PMS provides a better mass resolution, a larger volume and assumes cosmological parameters consistent with latest results from the Planck mission. The model includes a treatment for the partition of cold gas into atomic and molecular (H$_2$) components; a better treatment for environmental processes; an updated modelling of cold gas accretion on Super-Massive Black Holes. We compare GAEA predictions based on the PMS, with model realizations based on other simulations in the Millennium Suite, showing that the new model provides a remarkable consistency for most statistical properties of galaxy populations. We interpret this as due to the interplay between AGN feedback and H$_2$-based SFR, as model versions considering only one of the two mechanisms do not show the same level of consistency. We then compare model predictions with available data for the galaxy 2-point correlation function (2pCF) in the redshift range 0<z$\lesssim$3. We show that GAEA runs correctly recover the main dependencies of the 2pCF as a function of stellar mass, star formation activity, HI-content and redshift for galaxies more massive than 10$^{9}$ M$_\odot$. These results suggest that our model correctly captures both the distribution of galaxy populations in the Large Scale Structure and the interplay between the main physical processes regulating their baryonic content, both for central and satellite galaxies. The model predicts a small redshift evolution of the clustering amplitude, that results in an overprediction of z$\sim$3 clustering strength with respect to the available estimates, but is still consistent with data within 1-$\sigma$ uncertainties.