Exchange-correlation kernel for perturbation dependent auxiliary functions in auxiliary density perturbation theory.

CONTEXT: Analytic exchange-correlation kernel formulations are of the outermost importance for density functional theory (DFT) perturbation calculations. In this paper, the working equation for the exchange-correlation kernel of the generalized gradient approximation (GGA) for perturbation dependent auxiliary functions is derived and discussed in the framework of auxiliary density functional theory (ADFT). The presented new formulation is extended to the unrestricted approach, too. A comprehensive discussion of the implementation of the GGA ADFT kernel, using either the native exchange-correlation functional implementations in deMon2k or the ones from the LibXC library, is given. Calculations with analytic exchange-correlation kernels are compared to their finite difference counterparts. The obtained results are in quantitative agreement. Nevertheless, analytic GGA ADFT kernel implementations show substantial improvement in the computational performance. Similar results are reported for analytic second derivatives of effective core potential (ECP) and model core potential (MCP) matrix elements when compared to their finite difference counterparts in molecular frequency analyses. METHOD: All calculations are performed in the framework of ADFT as implemented in deMon2k. In the ADFT analytic frequency calculations, auxiliary density perturbation theory was used. The underlying two-center exchange-correlation kernel matrix elements are calculated by numerical integration either with analytic or finite difference kernel expressions. Validation calculations are performed with the VWN and PBE functionals employing DFT-optimized DZVP basis sets in conjunction with automatically generated GEN-A2 auxiliary density function sets. In the (Pt3Cu)n cluster benchmark calculations, the RPBE functional was used. For Pt atoms, the quasi-relativistic LANL2DZ effective core potential with the corresponding valence basis set was employed, whereas for Cu atoms, the all-electron DFT-optimized TZVP basis was applied. The auxiliary density was expanded by the automatically generated GEN-A2* auxiliary function set. We run all benchmark calculations in parallel on 24 cores.

Electron superhalogens as positronium superhalogens.

Positronium (Ps) exhibits the ability to form energetically stable complexes with atoms and molecules before annihilation occurs. In particular, F, a halogen, shows the highest reported positronium binding energy (2.95 eV) in the periodic table. Superhalogens are defined as molecules with electron affinities exceeding that of Cl (3.61 eV), the atom with the highest electron affinity. Building upon the concept of superhalogens, we can define Ps-superhalogens as molecules with Ps binding energies surpassing that of F. This study explores structural and energetic aspects of positronium and positron binding to neutral and anionic superhalogen molecules of the MXk+1 family (M = Li, Na, Be, Mg, B, Al, Si, P; X = F, Cl, Br), respectively and where k represents the highest formal valence of M. We perform multicomponent MP2 calculations for positron systems, which reveal how positron affinities vary with the type and number of halogen atoms present. The analysis of the results emphasizes the predominant role of electrostatic interactions in determining the positron affinity, with negligible effects of electronic and geometric relaxation upon positron attachment. We predict the energetic stability of 22 of the 24 PsMXk+1 complexes with respect to the chemically relevant dissociation channels: e+ emission, Ps emission and M-X bond breaking. Our findings reveal six MFk+1 systems that qualify as Ps-superhalogens, showing a positronium binding energy exceeding 2.95 eV. Of these, AlF4 stands out by setting a new record for the highest positronium binding energy among neutral molecules, reaching 4.36 eV.