Machine learning density functionals from the random-phase approximation - arXiv 2308.00665v1

Aug 4, 2026

Kohn-Sham density functional theory (DFT) is widely used in computational chemistry and materials science, yet more accurate methods like the random-phase approximation (RPA) remain difficult to apply because of their high computational cost. This work builds a DFT substitute (ML-RPA) functional for RPA by combining supervised and unsupervised machine learning, yielding a non-local extension of the standard gradient approximation. The model is trained for diamond surfaces and liquid water and demonstrates higher accuracy than standard gradient functionals. ML-RPA extends RPA’s applicability to larger systems, longer time scales, and broader chemical spaces.

LanguageEnglish

CategoryResearch & Report

FormatPDF

Pages69

Reading Time2 h 54 min

AccessFree to download

Abstract

1 Introduction

arXiv :2308 .00665v1 [ cond-mat .mtrl-sci ] 1 Aug 2023 Machine learning density functionals from the random-phase approximation Stefan Riemelmoser, ∗ ,†,‡ Carla Verdi,†,¶ Merzuk Kaltak, § and Georg Kresse†, § †University of Vienna, Faculty of Physics and Center for Computational Materials Science, Kolingasse 14-16, A-1090 Vienna, Austria ‡University of Vienna, Vienna Doctoral School in Physics, Boltzmanngasse 5, A-1090 Vienna, Austria ¶School of Physics, The University of Sydney, Sydney, NSW 2006, Australia §VASP Software GmbH, Sensengasse 8/12, A-1090 Vienna, Austria E-mail: [stefan.riemelmoser@univie.ac.at](stefan.riemelmoser@univie.ac.at) Abstract Kohn-Sham density functional theory (DFT) is the standard method for first-principles calculations in computational chemistry and materials science. More accurate theories such as the random-phase approximation (RPA) are limited in application due to their large computational cost. Here, we construct a DFT substitute functional for the RPAusing supervised and unsupervised machine learning (ML) techniques. Our ML-RPA model can be interpreted as a non-local extension to the standard gradient approximation. We train an ML-RPA functional for diamond surfaces and liquid water and show that ML-RPA can outperform the standard gradient functionals in terms of accuracy. Our work demonstrates how ML-RPA can extend the applicability of the RPA to larger system sizes, time scales and chemical spaces. 1 Introduction For over half a century, generations of researchers have been looking to find ever better approximations for the elusive exchange-correlation (xc) functional of Kohn-Sham density functional theory (DFT) . 1,2 The Hohenberg-Kohn theorems guarantee this exchange-correlation functional to be a universal functional of the electronic density n, in other words, there exists a map n →7 Exc [n] .3 However, the exact functional is of complex non-local form and in general unknown. Common approximations for the exchange-correlation energy are built on the principle of near-sightedness: 4 in the local density approximation (LDA), 5 the exchangecorrelation energy density is approximated to be locally that of a homogeneous electron gas with the same density. To improve the accuracy of the approximation, one can add more non-local information by climbing up Jacob’s ladder of DFT.6 Next to the local density, one includes its gradient |∇n| as ingredient to the exchange-correlation functional (generalized gradient approximation or GGA) . Going beyond the GGA, more non-local information is added in the form of the kinetic energy density (meta-GGA), and fractions of exact exchange (hybrids) . However, strictly speaking, most meta-GGA and hybrid functional approaches deviate from “pure”Kohn-Sham DFT, as the ingredients go beyond the electronic density alone. 7 A different route beyond GGA are non-local van der Waals (vdW) functionals, 8,9 which account for pair-wise dispersion interactions between densities n (r) and n (r′) . As the ingredients are only the electronic density and its gradient at points r and r′, the non-local vdW method stays within pure KS-DFT. Given ingredients X for the exchange-correlation functional, one still has to find the actual functional form, that is the map X →7 Exc [X] . In the analytical approach developed by Perdew and others, 10–12 those maps are found by satisfying exact constraints. The empirical approach pursued by Becke and others 13–16 optimizes a small number of adjustable parameters to experimental data or higher level theory. Broadly speaking, the analytical functionals are more universally applicable, whereas the empirical approach can achieve higher accuracy for systems similar to the ones represented in their respective training sets. 2 This comes atthe cost of worse performance for systems not represented by these training sets. For example, the widely used B3LYP functional, 17 whose parameters have been optimized for small

FAQs

Why is the random-phase approximation (RPA) difficult to apply in practice?

RPA is limited by its large computational cost, making it challenging for routine use in first-principles calculations.

How does the proposed ML-RPA functional relate to standard DFT approximations?

ML-RPA is constructed as a non-local extension to the standard gradient approximation, integrating machine learning to build the substitute for RPA.

What systems are used to train and evaluate the ML-RPA functional?

The functional is trained and assessed for diamond surfaces and liquid water.

Aug 4, 2026

Upload by Olivia Brown