We perform a variational Monte Carlo study of the nematic state in iron-pnictide superconductors within a realistic five-orbital model. Our numerical results show that the nematic state, formed by introducing an anisotropic hopping order into the projected wave function, is not stable unless the off-site Coulomb interaction V exceeds a critical value. This demonstrates that V plays a key role in forming the nematic state in iron-pnictide superconductors. In the nematic state,the orbital order and the anisotropic spin correlations are consistent with the experimental observations. We argue that the experimentally observed anisotropic magnetic couplings and structural transition are associated with the nematic state and can be understood in a unified framework.
