Friedel applied this theorem for two cases: (1) the onedimensional square box and (2) the thin spherical well -for both, he obtained the oscillatory interference patterned decaying with the distance from the surface/well. Moreover, he obtained a similar oscillatory pattern caused by any spherical potential using partial-wave analysis from the scattering theory. His method, though, requires analyzing all of the possible electron wavefunctions, therefore effectively making it computationally expensive and impractical in real applications. Here, we should also note that the decay of the electron density disturbance manifested by Friedel oscillations depends on the dimensionality of the system studied. The latter becomes a limiting factor in the numerical studies, as the theoretical model should be sufficiently extended in each direction. Even with modern available DFT codes, such scale is prohibitively expensive. Green's function approach allows for an increase in the efficacy of such calculations since it erases the need to calculate multiple wavefunctions. However, a readyto-use computational method that uses it in the context of FOs is yet to be developed. The absence of fast and easy-to-use methods to account for FOs (and their long-range character) has since limited their computational studies.