Beating the Egg-box Effect in Plane-Wave DFT Simulations

TY - JOUR

T1 - Beating the Egg-box Effect in Plane-Wave DFT Simulations

AU - Durham, Ben

AU - Probert, Matt

AU - Hasnip, Philip James

N1 - © 2025 The Author(s).

PY - 2025/5/7

Y1 - 2025/5/7

N2 - The "egg-box effect" is a known challenge in Density Functional Theory (DFT) calculations which arises from the discretization of continuous quantities, e.g. the electron density. This effect is observed when the system is moved relative to the underlying computational grid and causes an unphysical change in the system's total energy, violating translational invariance. The pattern of energy change with translation is reminiscent of an egg-box. This effect can cause unphysical results such as geometry relaxations finding incorrect crystal symmetries or imaginary phonon modes in vibrational calculations. The egg-box effect can be mitigated by using finer grids for all continuous quantities, but this greatly increases computational cost. For plane-wave DFT the effect's origin is the evaluation of the Exchange and Correlation (XC) energy (EXC ). We present a novel technique for estimating the violation of translational invariance in EXC , using a Fourier interpolation scheme, providing an estimate of the uncertainty. Our results also show that the numerical behaviour of the XC approximations is strongly linked to the magnitude of the violation of the translation invariance. The more numerically ill behaved XC functionals, such as the more advanced meta-GGAs functionals, exhibit changes in EXC (due to violation of translation invariance) that are orders of magnitude worse than less advanced functionals. Performing this analysis at an early stage of a workflow can inform the user about the expected accuracy of subsequent calculations. Further, our results demonstrate that by selectively computing EXC and its corresponding potential (VXC ) on a finer grid, the egg-box effect can be significantly reduced. Coupled with our uncertainty quantification method, egg-box related inaccuracies can be avoided more conveniently and efficiently than just increasing grid resolution until inaccuracies appear suppressed. This work offers a promising pathway towards mitigating the egg-box effect in a diverse range of materials modelling applications.

AB - The "egg-box effect" is a known challenge in Density Functional Theory (DFT) calculations which arises from the discretization of continuous quantities, e.g. the electron density. This effect is observed when the system is moved relative to the underlying computational grid and causes an unphysical change in the system's total energy, violating translational invariance. The pattern of energy change with translation is reminiscent of an egg-box. This effect can cause unphysical results such as geometry relaxations finding incorrect crystal symmetries or imaginary phonon modes in vibrational calculations. The egg-box effect can be mitigated by using finer grids for all continuous quantities, but this greatly increases computational cost. For plane-wave DFT the effect's origin is the evaluation of the Exchange and Correlation (XC) energy (EXC ). We present a novel technique for estimating the violation of translational invariance in EXC , using a Fourier interpolation scheme, providing an estimate of the uncertainty. Our results also show that the numerical behaviour of the XC approximations is strongly linked to the magnitude of the violation of the translation invariance. The more numerically ill behaved XC functionals, such as the more advanced meta-GGAs functionals, exhibit changes in EXC (due to violation of translation invariance) that are orders of magnitude worse than less advanced functionals. Performing this analysis at an early stage of a workflow can inform the user about the expected accuracy of subsequent calculations. Further, our results demonstrate that by selectively computing EXC and its corresponding potential (VXC ) on a finer grid, the egg-box effect can be significantly reduced. Coupled with our uncertainty quantification method, egg-box related inaccuracies can be avoided more conveniently and efficiently than just increasing grid resolution until inaccuracies appear suppressed. This work offers a promising pathway towards mitigating the egg-box effect in a diverse range of materials modelling applications.

U2 - 10.1088/2516-1075/adc056

DO - 10.1088/2516-1075/adc056

M3 - Article

SN - 2516-1075

VL - 7

JO - Electronic Structure

JF - Electronic Structure

M1 - 025004

ER -