Virtual-crystal approximation that works: Locating a compositional phase boundary in Pb(Zr 1-x Ti x )O 3