We present an algorithm for computing explicitly correlated second- and third-order Møller–Plesset energies near the basis set limit for large molecules with a cost that scales formally as 풩4 with system size 풩. This is achieved through a hybrid approach where locality is exploited first through orbital specific virtuals (OSVs) and subsequently through pair natural orbitals (PNOs) and integrals are approximated using density fitting. Our method combines the low orbital transformation costs of the OSVs with the compactness of the PNO representation of the doubles amplitude vector. The 풩4 scaling does not rely upon the a priori definition of domains, enforced truncation of pair lists, or even screening and the energies converge smoothly to the canonical values with decreasing occupation number thresholds, used in the selection of the PNO basis. For MP2.5 intermolecular interaction energies, we find that 99% of benchmark basis set limit correlation energy contributions are recovered using an aug-cc-pVTZ basis and that on average only 50 PNOs are required to correlate the significant orbital pairs.