  The design of better exchange-correlation (XC) functionals for Density
Functional Theory (DFT) is a central challenge of modern electronic structure
theory. However, current developments are limited by the mathematical form of
the functional, with efficient semi-local functionals being inaccurate for many
technologically important systems and the more accurate hybrid functionals
being too expensive for large solid-state systems due to the use of the exact
exchange operator. In this work, we use physics-informed machine learning to
design an exchange functional that is both orbital-dependent and nonlocal, but
which can be evaluated at roughly the cost of semi-local functionals and is
significantly faster than hybrid DFT in plane-wave codes. By training
functionals with several different feature sets, we elucidate the roles of
orbital-dependent and nonlocal features in learning the exchange energy and
determine that both types of features provide vital and independently important
information to the model. Having trained our new exchange functional with an
expressive, nonlocal feature set, we substitute it into existing hybrid
functionals to achieve hybrid-DFT accuracy on thermochemical benchmark sets and
improve the accuracy of band gap predictions over semi-local DFT. To
demonstrate the scalability of our approach as well as the practical benefits
of improved band gap prediction, we compute charged defect transition levels in
silicon using large supercells. Due to its transferability and computational
efficiency for both molecular and extended systems, our model overcomes the
cost-accuracy trade-off between semi-local and hybrid DFT, and our general
approach provides a feasible path toward a universal exchange-correlation
functional with post-hybrid DFT accuracy and semi-local DFT cost.
