Track the frontier, from condensed matter to cosmology.
Search papers, labs, and topics across Lattice.
The authors natively implement the machine-learned Skala exchange-correlation functional into CP2K using unified Gaussian and plane-wave/augmented-plane-wave (GPW/GAPW) frameworks with joint one-centre reconstruction. Extending neural-network functionals to condensed-phase systems bridges a critical gap between molecular machine-learned DFT and periodic solids while avoiding the severe plane-wave cutoffs that typically hinder augmented grid schemes. Across benchmark solids, the implementation matches diffusion Monte Carlo lattice energies for ice polymorphs and molecular crystals within ~1 kJ mol$^{-1}$ and yields near-hybrid band gap accuracy (0.43 eV MAE), despite exhibiting systematic equilibrium lattice overbinding.
Machine-learned exchange-correlation functionals can now scale directly to condensed-phase matter in CP2K, matching diffusion Monte Carlo lattice energy accuracy across 13 ice phases while rivaling hybrid functionals in band-gap prediction.
We implement the machine-learned Skala exchange-correlation (XC) functional natively in CP2K using its Gaussian and plane-wave (GPW) and Gaussian and augmented-plane-wave (GAPW) methods. A joint one-centre reconstruction of density, density gradients, and kinetic-energy density before functional evaluation preserves mixed gradient terms and nonlocal couplings between smooth and atom-local contributions in all-electron and pseudopotential calculations. An XC-specific GAPW representation resolves rapidly varying local contributions on atom-centred grids, reducing plane-wave requirements for pseudopotential calculations. The implementation includes Brillouin-zone sampling and point-group symmetry reduction. With sufficiently flexible orbital bases, it retains the molecular benchmark accuracy of our earlier GauXC formulation. All-electron Skala-D3(BJ) calculations yield a mean absolute error (MAE) of 1.54 kJ mol$^{-1}$ against diffusion Monte Carlo for crystalline CO$_2$, NH$_3$, and urea, spanning dispersion, quadrupolar electrostatics, and hydrogen bonding. For all thirteen DMC-ICE13 phases, the MAEs are 1.16 kJ mol$^{-1}$ for absolute lattice energies and 0.82 kJ mol$^{-1}$ for the twelve relative energies to ice Ih. At fixed experimental geometries, all-electron/mixed-core calculations give a band-gap MAE of 0.43 eV on a 15-material set, substantially below commonly used semilocal functionals and comparable to widely used hybrid functionals. The LC10 benchmark nevertheless reveals systematically underestimated equilibrium lattice constants, indicating structural overbinding. This framework connects molecular Skala to condensed-phase electronic structure and enables future development using periodic many-body reference data.