Practical Post-Kohn-Sham Methods for Time-Reversal Symmetry Breaking References

Beyond Electrons: Correlation and Self-Energy in Multicomponent Density Functional Theory

Post-Kohn-Sham methods are used to evaluate the ground-state correlation energy and the orbital self-energy of systems consisting of multiple flavors of different fermions. Starting from multicomponent density functional theory, suitable ways to arrive at the corresponding multicomponent random-phase approximation and the multicomponent Green's functionG W ${GW}$ approximation, including relativistic effects, are outlined. Given the importance of both of this methods in the development of modern Kohn-Sham density functional approximations, this work will provide a foundation to design advanced multicomponent density functional approximations. Additionally, theG W ${GW}$ quasiparticle energies are needed to study light-matter interactions with the Bethe-Salpeter equation.

Non-linear light-matter interactions from the Bethe-Salpeter equation

A route to assess non-linear light-matter interactions from the increasingly popular GW-Bethe-Salpeter equation (GW-BSE) method is outlined. In the present work, the necessary analytic expressions within the static-screened exchange approximation of the BSE are derived. This enables a straightforward implementation of the computation of the first hyperpolarizability as well as two-photon absorption processes for molecular systems. Benchmark calculations on small molecular systems reveal that the GW-BSE method is intriguingly accurate for predicting both first hyperpolarizabilities and two-photon absorption strengths. Using state-of-the-art Kohn-Sham references as a starting point, the accuracy of the GW-BSE method rivals that of the coupled-cluster singles-and-doubles method, outperforming both second-order coupled-cluster and time-dependent density-functional theory.

Insight on Gaussian Basis Set Truncation Errors in Weak to Intermediate Magnetic Fields with an Approximate Hamiltonian

Strong magnetic fields such as those found on white dwarfs have significant effects on the electronic structures of atoms and molecules. However, the vast majority of molecular studies in the literature in such fields are carried out with Gaussian basis sets designed for zero field, leading to large basis set truncation errors [Lehtola et al., Mol. Phys. 2020, 118, e1597989]. In this work, we aim to identify the failures of the Gaussian basis sets in atomic calculations to guide the design of new basis sets for strong magnetic fields. We achieve this by performing fully numerical electronic structure calculations at the complete basis set (CBS) limit for the ground state and low lying excited states of the atoms 1 ≤ Z ≤ 18 in weak to intermediate magnetic fields. We also carry out finite-field calculations for a variety of Gaussian basis sets, introducing a real-orbital approximation for the magnetic-field Hamiltonian. Our primary focus is on the aug-cc-pVTZ basis set, which has been used in many works in the literature. A study of the differences in total energies of the fully numerical CBS limit calculations and the approximate Gaussian basis calculations is carried out to provide insight into basis set truncation errors. Examining a variety of states over the range of magnetic field strengths from B = 0 to B = 0.6B0, we observe significant differences for the aug-cc-pVTZ basis set, while much smaller errors are afforded by the benchmark-quality AHGBSP3-9 basis set [Lehtola, J. Chem. Phys. 2020, 152, 134108]. This suggests that there is considerable room to improve Gaussian basis sets for calculations at finite magnetic fields.

TURBOMOLE: Today and Tomorrow

TURBOMOLE is a highly optimized software suite for large-scale quantum-chemical and materials science simulations of molecules, clusters, extended systems, and periodic solids. TURBOMOLE uses Gaussian basis sets and has been designed with robust and fast quantum-chemical applications in mind, ranging from homogeneous and heterogeneous catalysis to inorganic and organic chemistry and various types of spectroscopy, light-matter interactions, and biochemistry. This Perspective briefly surveys TURBOMOLE's functionality and highlights recent developments that have taken place between 2020 and 2023, comprising new electronic structure methods for molecules and solids, previously unavailable molecular properties, embedding, and molecular dynamics approaches. Select features under development are reviewed to illustrate the continuous growth of the program suite, including nuclear electronic orbital methods, Hartree-Fock-based adiabatic connection models, simplified time-dependent density functional theory, relativistic effects and magnetic properties, and multiscale modeling of optical properties.

Two-Component GW Calculations: Cubic Scaling Implementation and Comparison of Vertex-Corrected and Partially Self-Consistent GW Variants

We report an all-electron, atomic orbital (AO)-based, two-component (2C) implementation of the GW approximation (GWA) for closed-shell molecules. Our algorithm is based on the space-time formulation of the GWA and uses analytical continuation (AC) of the self-energy, and pair-atomic density fitting (PADF) to switch between AO and auxiliary basis. By calculating the dynamical contribution to the GW self-energy at a quasi-one-component level, our 2C-GW algorithm is only about a factor of 2-3 slower than in the scalar relativistic case. Additionally, we present a 2C implementation of the simplest vertex correction to the self-energy, the statically screened G3W2 correction. Comparison of first ionization potentials (IPs) of a set of 67 molecules with heavy elements (a subset of the SOC81 set) calculated with our implementation against results from the WEST code reveals mean absolute deviations (MAD) of around 70 meV for G0W0@PBE and G0W0@PBE0. We check the accuracy of our AC treatment by comparison to full-frequency GW calculations, which shows that in the absence of multisolution cases, the errors due to AC are only minor. This implies that the main sources of the observed deviations between both implementations are the different single-particle bases and the pseudopotential approximation in the WEST code. Finally, we assess the performance of some (partially self-consistent) variants of the GWA for the calculation of first IPs by comparison to vertical experimental reference values. G0W0@PBE0 (25% exact exchange) and G0W0@BHLYP (50% exact exchange) perform best with mean absolute deviations (MAD) of about 200 meV. Explicit treatment of spin-orbit effects at the 2C level is crucial for systematic agreement with experiment. On the other hand, eigenvalue-only self-consistent GW (evGW) and quasi-particle self-consistent GW (qsGW) significantly overestimate the IPs. Perturbative G3W2 corrections increase the IPs and therefore improve the agreement with experiment in cases where G0W0 alone underestimates the IPs. With a MAD of only 140 meV, 2C-G0W0@PBE0 + G3W2 is in best agreement with the experimental reference values.

Robust relativistic many-body Green's function based approaches for assessing core ionized and excited states

A two-component contour deformation (CD) based GW method that employs frequency sampling to drastically reduce the computational effort when assessing quasiparticle states far away from the Fermi level is outlined. Compared to the canonical CD-GW method, computational scaling is reduced by an order of magnitude without sacrificing accuracy. This allows for an efficient calculation of core ionization energies. The improved computational efficiency is used to provide benchmarks for core ionized states, comparing the performance of 15 density functional approximations as Kohn-Sham starting points for GW calculations on a set of 65 core ionization energies of 32 small molecules. Contrary to valence states, GW calculations on core states prefer functionals with only a moderate amount of Hartree-Fock exchange. Moreover, modern ab initio local hybrid functionals are also shown to provide excellent generalized Kohn-Sham references for core GW calculations. Furthermore, the core-valence separated Bethe-Salpeter equation (CVS-BSE) is outlined. CVS-BSE is a convenient tool to probe core excited states. The latter is tested on a set of 40 core excitations of eight small inorganic molecules. Results from the CVS-BSE method for excitation energies and the corresponding absorption cross sections are found to be in excellent agreement with those of reference damped response BSE calculations.