The two-dimensional Gross-Neveu model is anticipated to undergo a crystalline phase transition at high baryon charge densities. This conclusion is drawn from the mean-field approximation, which closely resembles models of Peierls instability. We demonstrate that this transition indeed occurs when both the rank of the symmetry group and the dimension of the particle representation contributing to the baryon density are large (the large N limit). We derive this result through the exact solution of the model, developing the large N limit of the Bethe ansatz. Our analytical construction of the large-N solution of the Bethe ansatz equations aligns perfectly with the periodic (finite-gap) solution of the Korteweg-de Vries (KdV) of the mean-field analysis.