"Natural" atomic orbitals are the optimal (maximum occupancy) effective AOs in the molecular environment. They are obtained as local eigenvectors of the 1-particle reduced density operator, and are hence intrinsic to the wavefunction itself.