We investigate the properties of a dilute gas of impurities embedded in an ultracold gas of bosons that forms a Bose–Einstein condensate (BEC). This work focuses mainly on the equation of state (EoS) of the impurity gas at zero temperature and the induced interaction between impurities mediated by the host bath. We use perturbative field-theory approaches, such as Hugenholtz– Pines formalism, in the weakly interacting regime. In turn, for strong interactions, we aim at non-perturbative techniques such as quantum–Monte Carlo (QMC) methods. Our findings agree with experimental observations for an ultra dilute gas of impurities, modeled in the framework of the single impurity problem; however, as the density of impurities increases, systematic deviations are displayed with respect to the one-body Bose polaron problem. © 2022 by the authors. Licensee MDPI, Basel, Switzerland.