Endofullerenes and Dispersion-Corrected Density Functional Approximations: A Cautionary Tale

Abstract: A recent study by Panchagnula et al. [J. Chem. Phys. 161, 054308 (2024)] illustrated the non-concordance of a variety of electronic structure methods at describing the symmetric double-well potential expected along the anisotropic direction of the endofullerene Ne@C${70}$. In this article we delve deeper into the difficulties of accurately capturing the dispersion interaction within this system, scrutinising a variety of state-of-the-art density-functional approximations (DFAs) and dispersion corrections (DCs). We identify rigorous criteria for the double-well potential and compare the shapes, barrier heights, and minima positions obtained with the DFAs and DCs to the correlated wavefunction data in the previous study, alongside new coupled-cluster calculations. We show that many of the DFAs are extremely sensitive to the numerical integration grid used, and note that the choice of DC is not independent of the DFA. Functionals with many empirical parameters tuned for main-group thermochemistry do not necessarily result in a reasonable PES, while improved performance can be obtained using nearly dispersionless DFAs with very few empirical parameters and allowing the DC to compensate. We pose the Ne@C${70}$ system as a challenge to functional developers and as a diagnostic system for testing dispersion corrections, and reiterate the need for more experimental data for comparison.