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 fu...