Besides the periodic boundary conditions, another primary difficulty lies in the scale of the models required for a relevant description of such a long-range effect. To accurately capture these oscillations, especially in metallic systems, models must be exceptionally large, often exceeding 1000 metal atoms or requiring extremely large unit cells. Such extensive models are necessary to ensure that the oscillations are not artificially influenced by the boundaries of the simulation cell, which can lead to erroneous results. The limitations of current computational resources further exacerbate these challenges. Simulating large models with high precision demands substantial computational power and memory. While advances in supercomputing and novel algorithms optimization schemes are continuously advancing, the stage where such extensive DFT calculations can be performed routinely and efficiently is not yet reached. Despite these difficulties, ongoing research in the field holds promise.