In this work, the recently developed “second-order” self-consistent method [Liu, Y., Ponte Castaneda, P., 2004a. Second-order estimates for the effective behavior and field fluctuations in viscoplastic polycrystals. J. Mech. Phys. Solids 52 467–495] is used to simulate texture evolution in halite polycrystals. This method makes use of a suitably optimized linear comparison polycrystal and has the distinguishing property of being exact to second order in the heterogeneity contrast. The second-order model takes into consideration the effects of hardening and of the evolution of both crystallographic and morphological texture to yield reliable predictions for the macroscopic behavior of the polycrystal. Comparisons of these predictions with full-field numerical simulations [Lebensohn, R.A., Dawson, P.R., Kern, H.M., Wenk, H.R., 2003. Heterogeneous deformation and texture development in halite polycrystals: comparison of different modeling approaches and experimental data. Tectonophysics 370 287–311], as well as with predictions resulting from the earlier “variational” and “tangent” self-consistent models, included here for comparison purposes, provide insight into how the underlying assumptions of the various models affect slip in the grains, and therefore the texture predictions in highly anisotropic and nonlinear polycrystalline materials. The “second-order” self-consistent method, while giving a softer stress-strain response than the corresponding full-field results, predicts a pattern of texture evolution that is not captured by the other homogenization models and that agrees reasonably well with the full-field predictions and with the experimental measures.