Track synthesis, materials and catalysis as it is published.
Search papers, labs, and topics across Lattice.
To determine whether neural network corrections to density functional approximations learn transferable physics or merely fit baseline-specific artifacts, the authors evaluated an identical semilocal machine-learning correction trained across four parent functionals: PBE, B3LYP, SCAN, and r$^2$SCAN. They discovered that the learned correction universally exhibits a strictly linear size-dependent drift in atomization energies ($|r| \ge 0.9999$ on $n$-alkanes), which fortuitously offsets baseline errors for PBE and B3LYP but actively degrades the already well-calibrated extensivity of SCAN and r$^2$SCAN. Spatial error decompositions further reveal that ML corrections target fundamentally divergent physical domains depending on the parent functional—core regions for PBE versus bonding regions for the others—confirming that ML augmentations act as bespoke error-cancellation masks rather than universal representations of missing correlation.
Machine learning corrections to density functionals do not learn transferable physics; instead, they function as parent-dependent error sinks that can catastrophically undermine size-extensivity when paired with modern meta-GGAs.
Machine-learning (ML) corrections to density functional approximations (DFAs) offer a route to improving electronic-structure predictions while retaining the efficiency of semilocal functionals. An important question is whether such corrections learn transferable improvements or instead compensate for errors specific to the parent functional. Here, we address this question by applying the semilocal ML correction of Wang \textit{et al.} [J.~Chem.~Phys.~\textbf{158}, 154107 (2023)] to PBE, B3LYP, SCAN, and r$^2$SCAN, while keeping the network architecture, loss function, training set, and optimization protocol unchanged. We find that the learned correction develops a systematic size-dependent contribution to atomization energies for all four parents, with $|r|\geq0.9999$ along the $n$-alkane series, where r is the Pearson coefficient of the linear fit. Its magnitude and sign depend systematically on the parent: the correction partially compensates the size-dependent errors of PBE and B3LYP, but introduces substantial size-dependent contributions for SCAN and r$^2$SCAN, whose parent errors show little size dependence. In contrast, ionization potentials, electron affinities, and isomerization energies are largely unaffected. Spatial decomposition further reveals distinct origins of this extensive contribution, with bonding regions dominating for B3LYP, SCAN, and r$^2$SCAN and core regions dominating for PBE. These results show that the global-energy ML correction is strongly dependent on the error structure of its parent DFA, highlighting size-dependent error accumulation as a critical consideration in the assessment and development of ML-enhanced density functionals.