2023MNRAS.523..453W

We present the first hydrodynamical cosmological simulations in the νHDM framework based on Milgromian dynamics (MOND) with light (11 eV) sterile neutrinos. νHDM can explain the expansion history, CMB anisotropies, and galaxy cluster dynamics similarly to standard cosmology while preserving MOND's successes on galaxy scales, making this the most conservative Milgromian framework. We generate initial conditions including sterile neutrinos using CAMB and MUSIC and modify the publicly available code PHANTOM OF RAMSES to run νHDM models. The simulations start at redshift $z_e$ = 199, when the gravitational fields are stronger than $a_{_0}$ provided this does not vary. We analyse the growth of structure and investigate the impact of resolution and box size, which is at most 600 comoving Mpc. Large density contrasts arise at late times, which may explain the KBC void and Hubble tension. We quantify the mass function of formed structures at different redshifts. We show that the sterile neutrino mass fraction in these structures is similar to the cosmic fraction at high masses (consistent with MOND dynamical analyses) but approaches zero at lower masses, as expected for galaxies. We also identify structures with a low peculiar velocity comparable to the Local Group, but these are rare. The onset of group/cluster-scale structure formation at $z_e$ ≈ 4 appears to be in tension with observations of high redshift galaxies, which we discuss in comparison to prior analytical work in a MONDian framework. The formation of a cosmic web of filaments and voids demonstrates that this is not unique to standard Einstein/Newton-based cosmology.