Efficient and low-scaling linear-response time-dependent density functional theory implementation for core-level spectroscopy of large and periodic systems

A new parameterization of the DFT/CIS method with applications to core-level spectroscopy

Time-dependent density functional theory (TD-DFT) within a restricted excitation space is an efficient means to compute core-level excitation energies using only a small subset of the occupied orbitals. However, core-to-valence excitation energies are significantly underestimated when standard exchange-correlation functionals are used, which is partly traceable to systemic issues with TD-DFT's description of Rydberg and charge-transfer excited states. To mitigate this, we have implemented an empirically modified combination of configuration interaction with single substitutions (CIS) based on Kohn-Sham orbitals, which is known as "DFT/CIS." This semi-empirical approach is well-suited for simulating x-ray near-edge spectra, as it contains sufficient exact exchange to model charge-transfer excitations yet retains DFT's low-cost description of dynamical electron correlation. Empirical corrections to the matrix elements enable semi-quantitative simulation of near-edge x-ray spectra without the need for significant a posteriori shifts; this should be useful in complex molecules and materials with multiple overlapping x-ray edges. Parameter optimization for use with a specific range-separated hybrid functional makes this a black-box method intended for both core and valence spectroscopy. Results herein demonstrate that realistic K-edge absorption and emission spectra can be obtained for second- and third-row elements and 3d transition metals, with promising results for L-edge spectra as well. DFT/CIS calculations require absolute shifts that are considerably smaller than what is typical in TD-DFT.

TDDFT and the x-ray absorption spectrum of liquid water: Finding the "best" functional

We investigate the performance of time-dependent density functional theory (TDDFT) for reproducing high-level reference x-ray absorption spectra of liquid water and water clusters. For this, we apply the integrated absolute difference (IAD) metric, previously used for x-ray emission spectra of liquid water [T. Fransson and L. G. M. Pettersson, J. Chem. Theory Comput. 19, 7333-7342 (2023)], in order to investigate which exchange-correlation (xc) functionals yield TDDFT spectra in best agreement to reference, as well as to investigate the suitability of IAD for x-ray absorption spectroscopy spectrum calculations. Considering highly asymmetric and symmetric six-molecule clusters, it is seen that long-range corrected xc-functionals are required to yield good agreement with the reference coupled cluster (CC) and algebraic-diagrammatic construction spectra, with 100% asymptotic Hartree-Fock exchange resulting in the lowest IADs. The xc-functionals with best agreement to reference have been adopted for larger water clusters, yielding results in line with recently published CC theory, but which still show some discrepancies in the relative intensity of the features compared to experiment.

Understanding X-ray absorption in liquid water using triple excitations in multilevel coupled cluster theory

X-ray absorption (XA) spectroscopy is an essential experimental tool to investigate the local structure of liquid water. Interpretation of the experiment poses a significant challenge and requires a quantitative theoretical description. High-quality theoretical XA spectra require reliable molecular dynamics simulations and accurate electronic structure calculations. Here, we present the first successful application of coupled cluster theory to model the XA spectrum of liquid water. We overcome the computational limitations on system size by employing a multilevel coupled cluster framework for large molecular systems. Excellent agreement with the experimental spectrum is achieved by including triple excitations in the wave function and using molecular structures from state-of-the-art path-integral molecular dynamics. We demonstrate that an accurate description of the electronic structure within the first solvation shell is sufficient to successfully model the XA spectrum of liquid water within the multilevel framework. Furthermore, we present a rigorous charge transfer analysis of the XA spectrum, which is reliable due to the accuracy and robustness of the electronic structure methodology. This analysis aligns with previous studies regarding the character of the prominent features of the XA spectrum of liquid water.

A Perspective on the Simulation of Electronic Circular Dichroism and Circularly Polarized Luminescence Spectra in Chiral Solid Materials

Chiral materials have shown tremendous potential for many technological applications, such as optoelectronics, sensing, magnetism, information technology, and imaging. Characterization of these materials is mostly based on chiroptical spectroscopies, such as electronic circular dichroism (ECD) and circularly polarized luminescence (CPL). These experimental measurements would greatly benefit from theoretical simulations for interpretation of the spectra as well as predictions on new materials. While ECD and CPL simulations are well established for molecular systems, they are not for materials. In this Perspective, we describe the theoretical quantities necessary to simulate ECD and CPL spectra in oriented systems. Then, we discuss the approximate strategies currently used to perform these calculations, what computational machinery is already available to develop more general approaches, and some of the open challenges for the simulation of ECD and CPL spectra in solid materials. When methods that are as reliable and computationally efficient as those for molecules are developed, these simulations will provide invaluable insight and guidance for the rational design of optically active materials.

Time-Dependent Density Functional Theory for X-ray Absorption Spectra: Comparing the Real-Time Approach to Linear Response

We simulate X-ray absorption spectra at elemental K-edges using time-dependent density functional theory (TDDFT) in both its conventional linear-response implementation and its explicitly time-dependent or "real-time" formulation. Real-time TDDFT simulations enable broadband spectra calculations without the need to invoke frozen occupied orbitals ("core/valence separation"), but we find that these spectra are often contaminated by transitions to the continuum that originate from lower-energy core and semicore orbitals. This problem becomes acute in triple-ζ basis sets, although it is sometimes sidestepped in double-ζ basis sets. Transitions to the continuum acquire surprisingly large dipole oscillator strengths, leading to spectra that are difficult to interpret. Meaningful spectra can be recovered by means of a filtering technique that decomposes the spectrum into contributions from individual occupied orbitals, and the same procedure can be used to separate L- and K-edge spectra arising from different elements within a given molecule. In contrast, conventional linear-response TDDFT requires core/valence separation but is free of these artifacts. It is also significantly more efficient than the real-time approach, even when hundreds of individual states are needed to reproduce near-edge absorption features and even when Padé approximants are used to reduce the real-time simulations to just 2-4 fs of time propagation. Despite the cost, the real-time approach may be useful to examine the validity of the core/valence separation approximation.

The Structure Properties of Carbon Materials Formed in 2,4,6-Triamino-1,3,5-Trinitrobenzene Detonation: A Theoretical Insight for Nucleation of Diamond-like Carbon

The structure and properties of nano-carbon materials formed in explosives detonation are always a challenge, not only for the designing and manufacturing of these materials but also for clearly understanding the detonation performance of explosives. Herein, we study the dynamic evolution process of condensed-phase carbon involved in 2,4,6-Triamino-1,3,5-trinitrobenzene (TATB) detonation using the quantum-based molecular dynamics method. Various carbon structures such as, graphene-like, diamond-like, and "diaphite", are obtained under different pressures. The transition from a C sp2- to a sp3-hybrid, driven by the conversion of a hexatomic to a non-hexatomic ring, is detected under high pressure. A tightly bound nucleation mechanism for diamond-like carbon dominated by a graphene-like carbon layer is uncovered. The graphene-like layer is readily constructed at the early stage, which would connect with surrounding carbon atoms or fragments to form the tetrahedral structure, with a high fraction of sp3-hybridized carbon. After that, the deformed carbon layers further coalesce with each other by bonding between carbon atoms within the five-member ring, to form the diamond-like nucleus. The complex "diaphite" configuration is detected during the diamond-like carbon nucleation, which illustrates that the nucleation and growth of detonation nano-diamond would accompany the intergrowth of graphene-like layers.

Electron-Affinity Time-Dependent Density Functional Theory: Formalism and Applications to Core-Excited States

The lack of particle-hole attraction and orbital relaxation within time-dependent density functional theory (TDDFT) lead to extreme errors in the prediction of K-edge X-ray absorption spectra (XAS). We derive a linear-response formalism that uses optimized orbitals of the n - 1-electron system as the reference, building orbital relaxation and a proper hole into the initial density. Our approach is an exact generalization of the static-exchange approximation that ameliorates the particle-hole interaction error associated with the adiabatic approximation and reduces errors in TDDFT XAS by orders of magnitude. With a statistical performance of just 0.5 eV root-mean-square error and the same computational scaling as TDDFT under the core-valence separation approximation, we anticipate that this approach will be of great utility in XAS calculations of large systems.

Computational approaches for XANES, VtC-XES, and RIXS using linear-response time-dependent density functional theory based methods

The emergence of state-of-the-art X-ray light sources has paved the way for novel spectroscopies that take advantage of their atomic specificity to shed light on fundamental physical, chemical, and biological processes both in the static and time domains. The success of these experiments hinges on the ability to interpret and predict core-level spectra, which has opened avenues for theory to play a key role. Over the last two decades, linear-response time-dependent density functional theory (LR-TDDFT), despite various theoretical challenges, has become a computationally attractive and versatile framework to study excited-state spectra including X-ray spectroscopies. In this context, we focus our discussion on LR-TDDFT approaches for the computation of X-ray Near-Edge Structure (XANES), Valence-to-Core X-ray Emission (VtC-XES), and Resonant Inelastic X-ray Scattering (RIXS) spectroscopies in molecular systems with an emphasis on Gaussian basis set implementations. We illustrate these approaches with applications and provide a brief outlook of possible new directions.

Spin-Orbit Couplings for Nonadiabatic Molecular Dynamics at the ΔSCF Level

A procedure for the calculation of spin-orbit coupling (SOC) at the delta self-consistent field (ΔSCF) level of theory is presented. Singlet and triplet excited electronic states obtained with the ΔSCF method are expanded into a linear combination of singly excited Slater determinants composed of ground electronic state Kohn-Sham orbitals. This alleviates the nonorthogonality between excited and ground electronic states and introduces a framework, similar to the auxiliary wave function at the time-dependent density functional theory (TD-DFT) level, for the calculation of observables. The ΔSCF observables of the formaldehyde system were compared to reference TD-DFT values. Our procedure gives all components (energies, gradients, nonadiabatic couplings, and SOC terms) at the ΔSCF level of theory for conducting efficient, full-atomistic nonadiabatic molecular dynamics with intersystem crossing, particularly in condensed phase systems.

Revisiting the Performance of Time-Dependent Density Functional Theory for Electronic Excitations: Assessment of 43 Popular and Recently Developed Functionals from Rungs One to Four

In this paper, the performance of more than 40 popular or recently developed density functionals is assessed for the calculation of 463 vertical excitation energies against the large and accurate QuestDB benchmark set. For this purpose, the Tamm-Dancoff approximation offers a good balance between computational efficiency and accuracy. The functionals ωB97X-D and BMK are found to offer the best performance overall with a root-mean square error (RMSE) of around 0.27 eV, better than the computationally more demanding CIS(D) wave function method with a RMSE of 0.36 eV. The results also suggest that Jacob's ladder still holds for time-dependent density functional theory excitation energies, though hybrid meta generalized-gradient approximations (meta-GGAs) are not generally better than hybrid GGAs. Effects of basis set convergence, gauge invariance correction to meta-GGAs, and nonlocal correlation (VV10) are also studied, and practical basis set recommendations are provided.

Relativistic Orbital-Optimized Density Functional Theory for Accurate Core-Level Spectroscopy

Core-level spectra of 1s electrons of elements heavier than Ne show significant relativistic effects. We combine advances in orbital-optimized density functional theory (OO-DFT) with the spin-free exact two-component (X2C) model for scalar relativistic effects to study K-edge spectra of third period elements. OO-DFT/X2C is found to be quite accurate at predicting energies, yielding a ∼0.5 eV root-mean-square error versus experiment with the modern SCAN (and related) functionals. This marks a significant improvement over the >50 eV deviations that are typical for the popular time-dependent DFT (TDDFT) approach. Consequently, experimental spectra are quite well reproduced by OO-DFT/X2C, sans empirical shifts for alignment. OO-DFT/X2C combines high accuracy with ground state DFT cost and is thus a promising route for computing core-level spectra of third period elements. We also explored K and L edges of 3d transition metals to identify limitations of the OO-DFT/X2C approach in modeling the spectra of heavier atoms.

All-Electron BSE@GW Method for K-Edge Core Electron Excitation Energies

We present an accurate computational approach to calculate absolute K-edge core electron excitation energies as measured by X-ray absorption spectroscopy. Our approach employs an all-electron Bethe-Salpeter equation (BSE) formalism based on GW quasiparticle energies (BSE@GW) using numeric atom-centered orbitals (NAOs). The BSE@GW method has become an increasingly popular method for the computation of neutral valence excitation energies of molecules. However, it was so far not applied to molecular K-edge excitation energies. We discuss the influence of different numerical approximations on the BSE@GW calculation and employ in our final setup (i) exact numeric algorithms for the frequency integration of the GW self-energy, (ii) G₀W₀ and BSE starting points with ∼50% of exact exchange, (iii) the Tamm-Dancoff approximation and (iv) relativistic corrections. We study the basis set dependence and convergence with common Gaussian-type orbital and NAO basis sets. We identify the importance of additional spatially confined basis functions as well as of diffuse augmenting basis functions. The accuracy of our BSE@GW method is assessed for a benchmark set of small organic molecules, previously used for benchmarking the equation-of-motion coupled cluster method [Peng et al., J. Chem. Theory Comput., 2015, 11, 4146], as well as the medium-sized dibenzothiophene (DBT) molecule. Our BSE@GW results for absolute excitation energies are in excellent agreement with the experiment, with a mean average error of only 0.63 eV for the benchmark set and with errors <1 eV for the DBT molecule.

On the choice of reference orbitals for linear-response calculations of solution-phase K-edge X-ray absorption spectra

X-ray absorption spectra of liquids calculated with linear-response theories like TDDFT and CIS are dramatically improved with core-ion reference orbitals.

First-principles correction scheme for linear-response time-dependent density functional theory calculations of core electronic states

Linear-response time-dependent density functional theory (LR-TDDFT) for core level spectroscopy using standard local functionals suffers from self-interaction error and a lack of orbital relaxation upon creation of the core hole. As a result, LR-TDDFT calculated x-ray absorption near edge structure spectra needed to be shifted along the energy axis to match experimental data. We propose a correction scheme based on many-body perturbation theory to calculate the shift from first-principles. The ionization potential of the core donor state is first computed and then substituted for the corresponding Kohn-Sham orbital energy, thus emulating Koopmans's condition. Both self-interaction error and orbital relaxation are taken into account. The method exploits the localized nature of core states for efficiency and integrates seamlessly in our previous implementation of core level LR-TDDFT, yielding corrected spectra in a single calculation. We benchmark the correction scheme on molecules at the K- and L-edges as well as for core binding energies and report accuracies comparable to higher order methods. We also demonstrate applicability in large and extended systems and discuss efficient approximations.