Résumé
The bispectrum of galaxy number counts is a key probe of large-scale structure (LSS), offering insights into the initial conditions of the Universe, the nature of gravity, and cosmological parameters. In this work, we derive the theoretical angular bispectrum of number counts for the first time without relying on the Limber approximation, while incorporating redshift binning. Notably, our analysis includes all Newtonian effects, leading relativistic projection effects, and general relativistic contributions, including radiation dynamics, up to second order in perturbation theory. For simplicity, however, we neglect any biasing effects. We have implemented these expressions in an open access code to evaluate the bispectrum for two redshift bins, $z=2 \pm 0.25$ and $z=0.6 \pm 0.05$, and compare our analytical results with simulations. For the contributions that already appear in a Newtonian treatment we find an interesting cancellation between the quadratic terms. At $z=2$, the projection effects and the dynamical effects have similar amplitude on large scales as we approach $k \sim \mathcal{H}$. In the squeezed limit, radiation effects are found to be the leading relativistic effects, by one order of magnitude. At $z=0.6$, we find the correct weak-field hierarchy between the terms, controlled by the ratio $\mathcal{H}/k$, but we still find that dynamical effects (from nonlinear evolution) are only a factor $2-3$ smaller than the projection effects. We compare our results with simulation measurements and find good agreement for the total bispectrum.