  Cosmological observables in perturbation theory turn out to be plagued with
infrared divergences, which represents both a conceptual and computational
challenge. In this paper we present a proof of concept for a systematic
procedure to remove these divergences in a large class of scalar cosmological
integrals and consistently define an infrared safe computable in perturbation
theory. We provide diagrammatic rules which are based on the nestohedra
underlying the asymptotic structure of such integrals.
