  We revisit the relation between the spherical model of Berlin-Kac and the
spin $O(N)$ model in the limit $N \to \infty$ and explain how they are related
via the discrete Gaussian free field (GFF).
  More precisely, using probabilistic limit theorems and concentration of
measure we first prove that the infinite volume limit of the spherical model on
a $d$-dimensional torus is a massive GFF in the high-temperature regime, a
standard GFF at the critical temperature, and a standard GFF plus a Rademacher
random constant in the low-temperature regime. The proof in the case of the
critical temperature appears to be new and requires a fine understanding of the
zero-average Green's function on the torus.
  For the spin $O(N)$ model, we study the model in the double limit of the
spin-dimensionality and the torus size. We take the limit as first the
spin-dimension $N$ goes to infinity, and then the size of the torus, and obtain
that the different spin coordinates become i.i.d. fields, whose distribution in
the high-temperature regime is a massive GFF, a standard GFF at the critical
temperature, and a standard GFF plus a Gaussian random constant in the
low-temperature regime.
  In particular, this means that although the limiting free energies per site
of the two models agree at all temperatures, their actual finite-dimensional
laws still differ in terms of their zero modes in the low-temperature regime.
