> > Unfortunately, forces are not correct in LDA+U (nor with spin-orbit > coupling). These terms have extra contributions to E-tot and thus also to > dE/dR = -F this means that optimization of atomic coordinates can be done within LDA+U only using an energy surface minimum criterion and not a zero-force criterion?