We consider a realistic bosonic N -particle model with unitary interactions relevant for Efimov physics. Using quantum Monte Carlo methods, we find that the critical temperature for BoseEinstein condensation is decreased with respect to the ideal Bose gas. We also determine the full momentum distribution of the gas, including its universal asymptotic behavior, and compare this crucial observable to recent experimental data. Similar to the experiments with different atomic species, differentiated solely by a three-body length scale, our model only depends on a single parameter. We establish a weak influence of this parameter on physical observables. In current experiments, the thermodynamic instability of our model from the atomic gas towards an Efimov liquid could be masked by the dynamical instability due to three-body losses.First predicted in 1970[1], the Efimov effect describes the behavior of three strongly interacting bosons when any two of them cannot bind. At unitarity, when the scattering length diverges, the three-body bound states are scale invariant and they form a sequence up to vanishing binding energy and infinite spatial extension. Efimov trimers had been intensely discussed in nuclear physics, but it was in an ultracold gas of caesium atoms that they were finally discovered [2]. To observe Efimov trimers, experiments in atomic physics rely on Feshbach resonances [3], which permits to instantly switch a gas between weak interactions and the unitary limit. Such a control of interactions lacks in nuclear physics or condensed matter experiments, and singular interactions can be probed there only in the presence of accidental fine tuning [4]. Beyond the original system[2], Efimov trimers have now been observed for several multi-component systems, including bosonic, fermionic and Bose-Fermi mixtures [5][6][7]. These experimental findings are interpreted in terms of the theory of few-body strongly interacting quantum systems. For three identical bosons in three dimensions, a complete universal theory is available, on and off unitarity [4]. Further theoretical work is aimed at understanding bound states for more than three bosons, mixtures, and the effects of dimensionality.Near-unitary interparticle interactions also impact the thermodynamics of the atomic gas, the description of which presents a challenge beyond the traditional theory of the Efimov effect. In addition, mean-field theory does not apply to infinite interactions [8], and the virial expansion [9] fails to describe the low-temperature state. Moreover, in atomic-physics experiments, strong interactions enhance the three-body loss rate, making the gas of bosons unstable. A characterization of the universal dynamics of these losses has been recently achieved [10][11][12]. On the other hand, a single breakthrough experiment [13] has addressed the low-temperature thermodynamics for a unitary bosonic gas, coming to the conclusion that equilibrium was approached faster than the system life-time. The importance of this system stems from its univer...