Similar to the natural orbitals, which are used to truncate the size of the orbital basis set, one can use the natural auxiliary function (NAF) approach of Kallay and co-workers to reduce the dimension of the auxiliary basis set. In the NAF approach, a singles value decomposition of the three-center two-electron integral matrix is performed.