Scalable Spin Squeezing in Two-Dimensional Arrays of Dipolar Large-<mml:math xmlns:mml="http://www.w3.org/1998/Math/MathML" display="inline"><mml:mi>S</mml:mi></mml:math> Spins

We theoretically show that the spin-spin interactions realized in two-dimensional Mott insulators of large-spin magnetic atoms (such as Cr, Er, or Dy) lead to scalable spin squeezing along the nonequilibrium unitary evolution initialized in a coherent spin state. An experimentally relevant perturbation to the collective squeezing dynamics is offered by a quadratic Zeeman shift, which leads instead to squeezing of individual spins. Making use of a truncated cumulant expansion for the quantum fluctuations of the spin array, we show that, for sufficiently small quadratic shifts, the spin squeezing dynamics is akin to that produced by the paradigmatic one-axis-twisting model-as expected from an effective separation between collective-spin and spin-wave variables. Scalable spin squeezing is shown to be protected by the robustness of long-range ferromagnetic order to quadratic shifts in the equilibrium phase diagram of the system that we reconstruct via quantum Monte Carlo and mean-field theory.