Entangled atomic states, such as spin squeezed states, represent a promising resource for a new generation of quantum sensors and atomic clocks. We demonstrate that optimal control techniques can be used to substantially enhance the degree of spin squeezing in strongly interacting many-body systems, even in the presence of noise and imperfections. Specifically, we present a protocol that is robust to noise which outperforms conventional methods. Potential experimental implementations are discussed.PACS numbers: 42.50. Dv, 02.30.Yy, 32.80.Qk Spin squeezed states are among the most interesting examples of entangled states. In quantum metrology they allow for measurements with an improved precision, ultimately limited only by the Heisenberg limit. Since the early theoretical proposals to realize them with non linear interactions [2,3], spin squeezed states have been implemented in several experiments. Specific examples include generation of spin squeezed states in cavity QED [4][5][6], in trapped ions through shared motional modes [7,8] or using a Bose-Einstein condensate [9,10].In this Letter we demonstrate that optimal control can be effectively employed to produce highly squeezed spin states in many-body quantum systems, drastically reducing the impact of relaxation and decoherence. Other approaches applied control techniques creating spin squeezing as a succession of unitary pulses of a constant Hamiltonian [11][12][13]. We employ the Chopped Random Basis (CRAB) technique [14,15] to optimally control the evolution of a collection of N two-level systems mutually coupled through a time-dependent non linear (i.e. quadratic) interaction. linear (i.e. quadratic) interaction. We calculate optimized evolutions occurring on time scales several orders of magnitude shorter than the corresponding adiabatic evolutions, with a speed-up increasing with the system size. Such a speed-up translates directly into an enhanced robustness of the squeezing in the presence of noise, as schematically depicted in Fig. 1. We illustrate this enhanced robustness by modelling two practical experimental implementations of squeezed state preparations: cavity QED and trapped ions [6,8].We will focus on two methods realizing spin squeezed states, both with advantages in different situations. The first is based on the so called one-axis twisting protocol, consisting in letting a collection of two-level systems evolve under the effect of a collective non linear interaction [3], described by a Hamiltonian of the formWhere ω is the precession frequency and χ is the strength of the nonlinear interaction and J is a collective spin operator (defined below). The relative simplicity of the one-axis twisting scheme has been at the basis of its ubiquitous presence in squeezing experiments; however such a scheme is known to be non optimal [3], the spherical nature of the angular momentum phase space limiting the maximal squeezing achievable. Such a bound is intrinsic for the one-axis twisting protocol with fixed χ. It nevertheless allows to achieve spi...