Under the auspices of the International Atomic Energy Agency, a series of well-documented benchmark experiments were performed at the Proteus facility of the Swiss Paul Scherrer Institute. Thirteen critical pebble bed reactor configurations were assembled, with ten of them deterministic with a precise location of the low-enriched fuel and moderator pebbles. Seven of these configurations were modeled with a very high spatial resolution with the Monte Carlo code MCNP with details that go from the fuel kernel (0.5 mm in diameter) to the walls surrounding the facility. The calculations of the k's of the configurations agree quite well with the experiments (within a fraction of a dollar). A sensitivity analysis is included to discuss the possibility of a small bias; also biases introduced by customary approximations of production codes were calculated. The experiments and the analysis of this paper might be very useful tools to check the calculational accuracy of procedures used in the emerging work related to pebble bed modular gas-cooled reactors.