Francesco Marino (U. Milan) Ab initio-based nuclear energy functionals: Constraints from the nuclear matter response Ab initio methods [1] hold the promise of allowing to determine all the structural properties of nuclei and infinite nucleonic matter, starting from the individual interactions between protons and neutrons. Despite recent remarkable progress, full-scale studies of heavy nuclei are still out of reach. By contrast, Density Functional Theory (DFT) [2] can be readily applied to ground-state and excited-state properties across the whole nuclear chart. However, the Energy Density Functionals (EDFs), the key quantities in DFT, are phenomenological and biased towards stable nuclei close to magicity. Shortcomings of the empirical EDFs manifest themselves far from the stability valley, e.g. in neutron-rich nuclei or pure neutron matter (PNM). To accompany the growing experimental research efforts devoted to unstable nuclei at the limits of the nuclear chart, there is a strong need for improving the theoretical tools at our disposal. In this contribution, we wish to discuss our approach, inspired by the "Jacob's ladder" of condensed matter DFT [3], that aims at constructing ab initio-constrained EDFs. We shall first present the first rung [4], called local density approximation (LDA), that exploits the equation of state of nucleonic matter as sole input. Then, we describe our current work, in which the static response [5] of both PNM and symmetric nuclear matter to an external potential is computed by means of ab initio techniques, namely the Quantum Monte Carlo and the Self-consistent Green's functions approach, and exploited to constrain the gradient terms of the nuclear EDF. [1] H. Hergert, Front. Phys. 8, 00379 (2020) [2] G. Colò, Advances in Physics: X 5, 1740061 (2020) [3] J. P. Perdew and K. Schmidt, AIP Conf. Proc. 577, 1 (2001) [4] F. Marino et al., Phys. Rev. C 104, 024315 (2021). [5] M. Buraczynski et al., Physics Letters B 818, 136347 (2021)`