I've been running some MTS DD type models for an FCC material. I'm finding that our global solver either using the plan newton raphson or newton raphson with a line search method is having trouble converging. It seems as if we could really benefit from that addition for these tougher problems.
Once, I've ported that solver we should probably donate it to the MFEM project as their community would also benefit from it as well.