Further savings can be achieved by truncating the virtual space as well. This requires more advanced treatments since not only the virtual space is much harder to localize, one needs to build a robust selection of orbitals a and b for each ij pair. For this purpose, one makes use of pair natural orbitals, which are obtained by diagonalizing (approximate) external MP2 pair density matrices