With the aim to make multicomponent wave function based correlation methods accessible for large systems, further approximations have to be applied to lower the computational scaling. This can for example be achieved by localizing the delocalized canonical orbitals within pair-natural orbitals. This localization lowers the quintic scaling of the canonical DF-NEO-MP2 to an asymptotic linear scaling of the local NEO-PNO-LMP2 method. In the context of this work, we make use of the PNO-LMP2 approach as proposed by Werner and coworkers and apply it to the electron-electron correlation. We shortly summarize the methodology, further details and associated works can be found in the original reference. The MP2 correlation energy can be constructed as the sum of orbital pair energies, with