When do machine-learned exchange-correlation improvements inherit into density-functional tight binding?
Authors: Can Polat, Mustafa Kurban, Erchin Serpedin, Hasan Kurban
Organizations: Department of Computer and Electrical Engineering, Texas A&M University, College Station, TX, USA · Department of Electrical and Computer Engineering, Texas A&M University at Qatar, Doha, Qatar · Department of Prosthetics and Orthotics, Ankara University, Ankara, Turkey · College of Science and Engineering, Hamad Bin Khalifa University, Doha, Qatar
Machine-learned exchange-correlation functionals correct band gaps at near-semilocal cost, while density-functional tight binding reaches the 103-106-atom regime; combining them assumes that a better parent yields a better parameterization, but we show it does not. Current-generation functionals are orbital-dependent generalized Kohn-Sham operators, whereas the parameterization channel is built on a multiplicative potential, preventing exact representation. Using the transfer ratio, the surviving fraction of a parent-level change, we find anti-transfer: coherently negative ratios across four covalent semiconductors move the gap in the wrong direction, consistent with a molecular proxy and an r2SCAN control. The minimal-basis overgap is dominated by the on-site convention rather than basis incompleteness; correcting the on-site block removes most of it, while one d-polarization shell closes a further 16-40, depending on the placement of the empty d level, which no free-atom eigenvalue uniquely fixes. Occupied-manifold enhancements, ionic and closed-shell repulsive potentials, and rocksalt-oxide gaps inherit, whereas elemental and III-V covalent networks inherit neither gaps nor repulsive potentials and oxide networks inherit only the latter. We screen 23 elements and release the parameter sets, showing that the transfer ratio provides a cheap pre-test before any parameterization campaign.
Figures & tables
Figure 1: The functional-agnostic parameterization bridge, and the three checks that certify it. a , The pipeline, with what each stage produces. The functional enters at exactly one point, the Kohn–Sham potential inversion, which is also where the compression becomes lossy: the two-center integrator consumes a multiplicative v(r) and nothing else, and an orbital-dependent functional admits no such v . b , Fidelity of the inversion against a stored reference potential. The route the bridge uses and a full potential reconstruction both recover it to within a few meV; the density-gradient-only approximation fails by two orders of magnitude and was rejected. c , The confinement sweep for the carbon reference. The maximum overlap deviation has a sharp minimum at the certified radius, which is what fixes the convention; the Hamiltonian root-mean-square difference does not, falling monotonically through it. The convention is certified against S and against S alone. d , Coverage. Elements passing the atomic feasibility gates across three element sets; the periodic machine-learned self-consistent field is the binding constraint.
Parent DFT
DFTB (free-atom)
Gap transfer ratio
System
ΔEg (eV)
ΔEg (eV) 1
Free-atom
Confined 2
Covalent semiconductors
Diamond
+2.459
−1.896
−0.771
+0.128
Silicon
+2.211
−2.323
−1.051
−0.325
3C-SiC
+3.595
−3.900
−1.085
−0.058
BP
+2.082
−2.029
−0.974
−0.024
Table 1: The band-gap correction does not merely fail to survive, it moves the tight-binding gap the wrong way. All free-atom entries are evaluated on the fundamental band-path gap at both levels. The magnesium-oxide gap change is of order 0.01 eV, so the ionic case is functional-insensitive rather than positively transferring. A ratio of one is full inheritance and 0.5 the useful-transfer threshold; the confined ratios, which the on-site convention error drove toward zero, are given for comparison in the last column.
Figure 2: The band-gap correction reverses under the compression, and what survives instead. a , Input against output: the parent gap opens (grey; dashed outlines mark full inheritance) while the compressed DFTB gap closes (red). The upper segment carries the scale this is read against, the residual overgap of the same tables under the same convention: +2.49 , +0.21 and +1.31 eV for diamond, silicon and 3C-SiC. b , The same result as a ratio, against the 0.5 useful-transfer threshold. Gap ratios for the covalent solids, the near-zero ionic case (MgO), and the bandwidth ratios (a width quantity, convention-independent) reported for the three solids that carry them; none clears the gate. c , The molecular proxy on N 2 , CO and CO 2 , with a parent correction five times larger. Parent and tight-binding responses are drawn as ranges. d , What survives. The free-atom on-site s – p opening persists as an occupied-width enhancement into the dimer and solid while the solid gap ratio barely moves, the partition the theory predicts.
Figure 3: Why the parameterization channel cannot carry an orbital-dependent functional, and the two measurements that decide whether it does. κ≡f′(0)/Zf(0) is the dimensionless number the functional fixes at the nucleus. a , The l -resolved local representative with the cusp condition imposed, at κ=0 . Each angular channel sees a different potential and the difference diverges as 1/r , so no multiplicative v(r) lies in the shaded band; only the s divergence cs=0 is parameter-free. b , The cusp coefficients against κ . Both are linear and vanish at different κ ( 1 and 2/3 ), so they never vanish together; the ratio cd/cp=4/3 at κ=0 only. c , Where the discarded gap-shift term acts. Sketched for a homonuclear pair; the argument does not run for ionic solids. d , Two measurements on a common scale. At the production inversion’s matching radius the channel splitting ∣vs−vp∣ is 5.92 eV, three orders above that inversion’s 6.4 meV noise floor; under PBE the theory predicts exactly zero and the negative control confirms it at the floor.
Figure 4: The overgap of the unfitted tables, what a single d shell removes of it, and the two things that would undo the claim. a , Γ -point gaps; red rules mark the parent (identical settings). The minimal sp tables overgap the parent by ∼10 eV; adding a d -polarized spd set (confinement radii and on-site eigenvalues held fixed) reduces the overgap several-fold (per-solid values in Table 2 ). The published-parameter baseline is given in Fig. 6 ; under the free-atom convention the excess is shared with fitted sets, not a property of our convention alone. b , The closure and its residual, drawn at the confined on-site convention, which isolates the basis-completeness effect; under the free-atom convention adopted in this work the closure is smaller and depends on the placement of the empty d level, as Fig. 5 shows. The residual grows with polarity and is intrinsic to the two-center on-site approximation, not removed by further basis enlargement (per-solid closures and residuals in Table 2 ). c , The on-site convention, tested and corrected. The standard scheme puts free -atom eigenvalues on the diagonal Hμμ0 ; the confined convention places confined eigenvalues there instead (every on-site εp>0 ). Rebuilding with free-atom eigenvalues collapses the sp overgap, confirming the convention was the dominant contributor to panel a .
Figure 5: Under the convention the overgap closure tracks the placement of the empty d level, which no free-atom eigenvalue fixes. a , Mean and per-solid closure as the carbon Ed is swept from the confined value ( +40.96 eV) toward the valence shell, with silicon shifted in proportion. The dashed rule marks the closure obtained when the confined convention is used on both the sp and spd sides. b , The Mulliken d population across the same sweep, against the eight valence electrons of the cell; placements further down, at the free-atom p level, take the d manifold to a valence-shell population and are excluded on that basis. All gaps are fundamental band-path gaps; n=3 solids, confinement radii and free-atom s and p on-site energies held fixed throughout.
Parent
Minimal sp (eV)
d -polarized spd (eV)
Closure
System
DFT (eV)
Gap
Overgap 1
Gap
Residual 2
× par.
(%) 3
Diamond
5.20
14.97
+9.78
6.83
+1.63
1.31
83
Silicon
2.28
11.81
+9.53
4.71
+2.43
2.07
75
3C-SiC
5.54
16.89
+11.34
10.46
+4.92
1.89
57
Mean
–
–
–
–
–
–
71
Table 2: One d shell removes most of the overgap at fixed on-site energies, and what it leaves behind grows with polarity. Gaps are reported at the confined-atom on-site convention, held fixed between the sp and spd sets so that this table isolates the basis-completeness effect alone. Both columns are Γ -point gaps for parent and model alike, so this is an overgap decomposition and not a transfer ratio. The free-atom closure and its dependence on the d -channel placement are given in Fig. 5 .
Figure 6: Band structures, densities of states, and the fitted-table baseline, in the free-atom on-site convention. a–c , Band structures along L– Γ –X–W–K– Γ with the total density of states, for diamond, silicon, and 3C-SiC under the semilocal reference (PBE, sp and spd ) and the machine-learned functional (CIDER23X, sp ); energies referenced to the valence-band maximum. d–f , Fundamental gaps of each functional-consistent table and of the fitted published sets, run in the same DFTB+ binary, against the experimental gap (dashed).
Figure 7: The bonding-class structure of the transfer failure. a , Gap inheritance per compound, against the 0.30 eV gate. Only the rocksalt oxide and chloride approach it; every covalent network diverges, and two ionic compounds do so as well, by distinct minimal-basis mechanisms. b , Held-out repulsive-potential transfer, scored as the mean-centered correlation between a potential fitted on one configuration set and tested on another, separating cleanly by bonding class. Mg–O is drawn in its own class colour as the one ionic pair that fails; the text diagnoses it as a near-flat residual rather than a fitting instability. c , The class separation and the swept criterion, which are different things. The shaded bands are the two class ranges; the gap between them is where a class-separating bar must fall. The swept criterion lies inside the covalent band and does not separate the classes, moving only the sulfur dimer. d , The per-battery means against the gate.
Figure 8: The application boundary: what the compressed model can and cannot be asked for. a , Three regimes on a common scale of deviation from the parent: same-manifold spin flips track it closely, gap-crossing promotions carry the ground-state overgap with them, and multi-determinant multiplets return an identically zero splitting under single-reference Δ SCF, which is a method boundary rather than a parameterization one. b , The silicon divacancy gap against supercell size, the inter-functional difference changing sign as the cell grows.
Quantity (Ha)
Carbon
Silicon
Permitted by the residual identity 1
l -resolved ∣vs−vp∣ , median
0.099
0.218
–
spσ at bond, machine-learned 2
+0.283
−6.21
10−2 to 10−1
spσ at bond, semilocal
−0.249
−0.13
–
Table S1 : l -resolved wall signature for carbon against silicon. Median l -resolved on-site potential difference from separate s - and p -channel inversions of the machine-learned confined atom, and the spσ two-center integral at the bond length under the inverted machine-learned and the semilocal potentials.
System
ΔWDFT (eV)
ΔWDFTB (eV)
TRbw
Diamond
+5.952
+2.734
+0.459
Silicon
+2.633
−1.328
−0.504
3C-SiC
+4.164
+0.294
+0.071
Table S2 : Occupied-bandwidth changes under the change of functional (PBE → CIDER23X), and the transfer ratio they form. Parent-DFT values are GPAW occupied bandwidths; DFTB values are the change in the tight-binding occupied bandwidth at the same geometry.
CIDER23X
CIDER24Xe
PBE control
Atom
p10
p50
p90
p10
p50
p90
p10
p50
p90
C
0.319
1.404
4.582
0.439
1.415
7.087
0.0074
0.0383
0.3041
N
0.295
2.040
7.663
0.294
2.041
11.878
0.0016
0.0064
0.0196
O
0.346
1.458
4.625
0.354
1.380
4.635
0.0113
0.0564
0.4586
Si
0.485
3.180
8.325
0.824
3.235
7.387
0.0112
0.0610
0.4496
Table S3 : Percentiles of the node-masked channel-potential asymmetry ∣vs−vp∣ (Ha), for two machine-learned functionals and the PBE control.
Prediction
What a failure would mean
Outcome
P1
∣vs−vp∣=0 exactly under PBE, where f≡0
the inversion is ill-conditioned, not the channel
holds: median 3.85 and 2.72 meV, at or below the 6.4 meV noise floor
P2
cs=0 for every κ ; cp and cd linear in κ with the functional’s own κ
the near-nucleus law is not the mechanism
holds: κ≈6.0 – 6.3 , far from the degenerate 2/3 ; naive l(l+1) excluded
Table S4 : The two predictions the theory makes that this study tests, and what a failure of each would mean. Both are computable from quantities the parameterization already produces, namely the functional evaluated on its own converged density and the confined-atom radial solutions, so neither requires new data. P1 is parameter-free: it decides whether the inversion is orbital-faithful at all, and everything downstream rests on it passing.
Parent gap (eV)
Parent
Transfer ratio
System
PBE
r 2 SCAN
ΔEg (eV)
Free-atom
Confined
Diamond
4.130
4.339
+0.210
−1.049
−1.731
Silicon
0.571
0.765
+0.194
−2.674
−1.051
3C-SiC
1.352
1.746
+0.394
−0.442
−0.343
Table S5 : Independent meta-GGA control. r 2 SCAN, carried through the same pipeline, anti-transfers like the machine-learned functionals: the parent gap opens while the compressed DFTB gap closes, giving a negative transfer ratio for every solid under the physical free-atom convention. Parent gaps are GPAW band-path fundamental gaps; the transfer ratio is the DFTB gap change divided by the parent gap change, r2SCAN−PBE .
Set
Elements
Files
Expected 1
Parent functional
PBE
C, N, O, Si
12
16
PBE (semilocal reference)
CIDER23X
C, Si
4
4
nonlocal-feature meta-GGA, energy-trained
CIDER24Xe
C, N, O
9
9
density-matrix features, gap-trained (molecular)
Total
25
29
Table S6 : Released Slater–Koster parameter sets. Files counts the extended-format .skf files (homonuclear plus ordered heteronuclear pairs). Provenance is the parent functional bridged through the confined-atom inversion.
Density Functional Theory (DFT) is widely used for first-principles simulations in chemistry and materials science, but its computational cost remains a key limitation for large systems. Motivated by recent advances in ML-based exchange-correlation (XC) functionals, this paper introduces a differentiable framework that integrates machine learning models into density functional theory (DFT) for solids and other periodic systems. The framework defines a clean API for neural network models that can act as drop in replacements for conventional exchange-correlation (XC) functionals and enables gradients to flow through the full self-consistent DFT workflow. The framework is implemented in Python using a PyTorch backend, making it fully differentiable and easy to use with standard deep learning tools. We integrate the implementation with the DeepChem library to promote the reuse of established models and to lower the barrier for experimentation. In initial benchmarks against established electronic structure packages (GPAW and PySCF), our models achieve relative errors on the order of 5-10%.
Rakshit Kumar Singh, Aryan Amit Barsainyan, Bharath Ramsundar
We introduce STRUCTURES-P, a machine-learned all-electron orbital-free density functional for periodic CHNO systems. Electrostatics are evaluated analytically in an atom-centered Gaussian basis, while an equivariant network learns the kinetic and exchange-correlation energy from energies and gradients of non-ground-state densities. Minimizing the functional from a superposition of atomic densities recovers the Kohn-Sham ground state: on held-out crystals, the density error falls from 6.7% to 0.4% of the electron count and the energy error reaches 0.25 mHa per atom, half that of a same-architecture interatomic potential. A single functional can be trained for both isolated molecules and periodic systems, and the model generalizes to larger unit cells. It largely preserves relative-energy rankings and reproduces convex-hull vertices on the tested H-N composition slice. Together, these results demonstrate the viability of machine-learned all-electron orbital-free DFT for variational density optimization and energy evaluation in fixed-geometry solids.
Marc K. Ickler, Andreas Hermann, Fred A. Hamprecht
Interdisciplinary Center for Scientific Computing (IWR), Heidelberg University, Heidelberg, 69120, Germany. · Centre for Science at Extreme Conditions and SUPA, School of Physics and Astronomy, University of Edinburgh, Edinburgh, EH9 3FD, United Kingdom.
Machine-learned (ML) exchange-correlation (XC) functionals aim to replace human-designed density functional approximations by learning directly from reference data, but they still do not consistently outperform traditional O(N4)-scaling hybrid functionals. We study a hybrid-distillation setting in which O(N3)-scaling ML-XC functionals are trained to reproduce B3LYP/def2-SVP targets. We introduce Derivative Informed XC-Loss (DI-Loss), a loss that incorporates additional information from the reference hybrid functional by supervising first and second derivatives of the energy on the Grassmannian of admissible density matrices. Rather than only matching the self-consistent fixed point, DI-Loss aligns the local first- and second-order response of the learned functional with that of the target functional. Across four evaluated architectures, DI-Loss consistently improves the main energy metrics. Averaged uniformly across architectures, the total-energy MAE decreases by 66% relative to energy and density supervision alone. The density-sensitive mean-field energy metric Eρ improves from 1.2 to 0.8 mEh on average, while dipole and L2 density errors do not improve uniformly. We further show that densities from the distilled functionals reduce hybrid-functional SCF iterations by up to 50%. In downstream TDDFT calculations, Hessian supervision improves excited-state predictions, with XCdiff reducing the mean excitation-energy MAE by 19 - 35%.
Eike S. Eberhard, Luca A. Thiede, Abdul Aldossary +5
Technical University of Munich, Munich, Germany · Munich Data Science Institute, Munich, Germany · Munich Center for Machine Learning, Munich, Germany +4