PSIXAS: A Psi4 plugin for efficient simulations of X-ray absorption spectra based on the transition-potential and Δ-Kohn–Sham method

Near edge X-ray absorption fine structure (NEXAFS) spectra and their pump-probe extension (PP-NEXAFS) offer insights into valence- and core-excited states. We present PSIXAS, a recent implementation for simulating NEXAFS and PP-NEXAFS spectra by means of the transition-potential and the Δ-Kohn–Sham method. The approach is implemented in form of a software plugin for the Psi4 code, which provides access to a wide selection of basis sets as well as density functionals. We briefly outline the theoretical foundation and the key aspects of the plugin. Then, we use the plugin to simulate PP-NEXAFS spectra of thymine, a system already investigated by others and us. It is found that larger, extended basis sets are needed to obtain more accurate absolute resonance positions. We further demonstrate that, in contrast to ordinary NEXAFS simulations, where the choice of the density functional plays a minor role for the shape of the spectrum, for PP-NEXAFS simulations the choice of the density functional is important. Especially hybrid functionals (which could not be used straightforwardly before to simulate PP-NEXAFS spectra) and their amount of “Hartree-Fock like” exact exchange affects relative resonance positions in the spectrum.     

Bringing about matrix sparsity in linear-scaling electronic structure calculations

The performance of linear-scaling electronic structure calculations depends critically on matrix sparsity. This article gives an overview of different strategies for removal of small matrix elements, with emphasis on schemes that allow for rigorous control of errors. In particular, a novel scheme is proposed that has significantly smaller computational overhead compared with the Euclidean norm-based truncation scheme of Rubensson et al. (J Comput Chem 2009, 30, 974) while still achieving the desired asymptotic behavior required for linear scaling. Small matrix elements are removed while ensuring that the Euclidean norm of the error matrix stays below a desired value, so that the resulting error in the occupied subspace can be controlled. The efficiency of the new scheme is investigated in benchmark calculations for water clusters including up to 6523 water molecules. Furthermore, the foundation of matrix sparsity is investigated. This includes a study of the decay of matrix element magnitude with distance between basis function centers for different molecular systems and different methods. The studied methods include Hartree–Fock and density functional theory using both pure and hybrid functionals. The relation between band gap and decay properties of the density matrix is also discussed.