The polarization of a nearest-neighbor (nn) Pb-O vacancy pair [(VPb-VO)nn] in PbTiO3 is calculated, using the modern theory of polarization, implemented in the density-functional-theory ultrasoft pseudopotential formalism. ...