The $^1S_0$ pairing in neutron matter has been investigated in presence of realistic two-- and three--nucleon interactions. We have adopted the Argonne $v_{8^\prime}$ NN and the Urbana IX 3N potentials. Quantum Monte Carlo theory, specifically the Auxiliary Field Diffusion Monte Carlo method, and Correlated Basis Function theory are employed in order to get quantitative and reliable estimates of the gap. They both fully take into account the medium modifications due to the interaction induced correlations. The two methods are in good agreement up to the maximum gap density and both point to a slight reduction with respect to the standard BCS value. In fact, the maximum gap is about $2.5 \text{MeV}$ at $k_F \sim 0.8 \text{fm}^{-1}$ in BCS an...