ABSTRACTWe investigate the exponential f(Q) symmetric teleparallel gravitation, namely $f(Q)=Q+\alpha Q_0(1-e^{-\beta \sqrt{Q/Q_0}})$ using ME-GADGET code to probe the structure formation with box sizes LBox = 10/100 Mpc h–1 and middle resolution $N_p^{1/3}=512$. To reproduce viable cosmology within the aforementioned modified gravity theory, we first perform Markov Chain Monte Carlo (MCMC) sampling on OHD/BAO/Pantheon data sets and constrain a parameter space. Furthermore, we also derive theoretical values for deceleration parameter q(z), statefinder pair {r, s} and effective gravitational constant Geff, perform Om(z) diagnostics. While carrying out N-body+SPH simulations, we derive CDM+baryons over density/temperature/mean molecular weight fields, matter power spectrum (both 2/3D, with/without redshift space distortions), bispectrum, two-point correlation function and halo mass function. Results for small and big simulation box sizes are therefore properly compared, halo mass function is related to the Seth–Tormen theoretical prediction and matter power spectrum to the standard CAMB output.