Search papers, labs, and topics across Lattice.
This paper details the implementation of the machine-learned Skala exchange-correlation functional within the CP2K simulation framework via the GauXC library, aiming to enhance Kohn-Sham density-functional theory without the computational burden of explicitly correlated methods. By establishing a robust interface that accommodates both all-electron and valence-only density matrices, the authors ensure accurate evaluations of XC energy and atomic-orbital potentials. The implementation achieves a mean absolute deviation of 1.255 kcal/mol, closely aligning with the Skala reference value, thereby validating its effectiveness for molecular simulations.
Achieving a mean absolute deviation of just 0.020 kcal/mol from the Skala reference value, this implementation sets a new standard for accuracy in density-functional theory simulations.
Machine-learned exchange--correlation (XC) functionals offer a route to improve Kohn--Sham density-functional theory without incurring the cost of explicitly correlated electronic-structure methods. Their use in production simulation codes, however, requires a well-defined mapping between the learned model and the host-code density representation. We formulate and implement a Skala-1.1 interface in CP2K through the external GauXC library. CP2K supplies the geometry, Gaussian basis, spin-resolved atomic-orbital density matrix, and communicator, while GauXC evaluates the XC energy, atomic-orbital potential matrix, and available nuclear derivatives. The interface accepts both all-electron and valence-only density matrices. The latter may arise from separable dual-space pseudopotentials or molecular effective-core potentials. Implementation errors are isolated from functional differences by comparing the Perdew--Burke--Ernzerhof (PBE) functional evaluated through GauXC with native CP2K PBE. The resulting interface gives consistent energies, forces validated against finite-difference total-energy checks, and force-based molecular-virial diagnostics for representative molecular cases. The dietGMTKN55 benchmark suite is evaluated with an all-electron Gaussian augmented plane-wave treatment for elements up to bromine and def2 effective-core potentials for the heavier elements. The resulting aggregate mean absolute deviation of 1.255 kcal/mol is within 0.020 kcal/mol of the corresponding Skala reference value of 1.235 kcal/mol. This work establishes a validated molecular implementation of Skala in CP2K through GauXC.
