When it comes to improving the accuracy of DFT calculations, there is a well known hierarchy:
- starting from LDA using $E_{xc}[n]$,
- proceeding with GGA and its $E_{xc}[ n, \nabla n]$ dependency
- and finally meta-GGAs with $E_{xc}[ n, \nabla n, \nabla^2n~\text{or}~\nabla^2\phi_i]$
I always wondered why higher order derivatives are not used but never really looked into it. Recently I heard from an expert in the field that it can be shown that $E_{xc}[n]$ can not converge in the sense of a Taylor series.
How? Why?
I just couldn't find more on this matter.