Abstract We study the behavior of Lorentzian Engle-Pereira-Rovelli-Livine spinfoam amplitude with homogeneous boundary data, under a graph refinement going from five to twenty tetrahedra. This can be interpreted as wave function Universe, for which we compute geometrical operators, correlation functions, and entanglement entropy. The numerical calculation is made possible by adapting Metropolis...