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 function G 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, the G W ${GW}$ quasiparticle energies are needed to study light-matter interactions with the Bethe-Salpeter equation.

Polarization-dependent effects in vibrational absorption spectra of 2D finite-size adsorbate islands on dielectric substrates.

In the last few years, infrared reflection-absorption spectroscopy (IRRAS) has become a standard technique to study vibrational excitations of molecules. These investigations are strongly motivated by potential applications in monitoring chemical processes. For a better understanding of the adsorption mechanism of molecules on dielectrics, the polarization-dependence of an interaction of infrared light with adsorbates on dielectric surfaces is commonly used. Thus, the peak positions in absorption spectra could be different for s- and p-polarized light. This shift between the peak positions depends on both the molecule itself and the dielectric substrate. While the origin of this shift is well understood for infinite two-dimensional adsorbate layers, finite-size samples, which consist of 2D islands of a small number of molecules, have never been considered. Here, we present a study on polarization-dependent finite-size effects in the optical response of such islands on dielectric substrates. The study uses a multi-scale modeling approach that connects quantum chemistry calculations with Maxwell scattering simulations. We distinguish the optical response of a single molecule, a finite number of molecules, and a two-dimensional adsorbate layer. We analyze CO and CO2 molecules deposited on CeO2 and Al2O3 substrates. The evolution of the shift between the polarization-dependent absorbance peaks is first studied for a single molecule, which does not exhibit any shifting at all, and for finite molecular islands, where it increases with increasing island size, as well as for an infinite two-dimensional adsorbate layer. In the latter case, the agreement between the obtained results and the experimental IRRAS data and more traditional three/four-layer model theoretical studies supports the predictive power of the multi-scale approach.

Paramagnetic Nuclear Magnetic Resonance Shifts for Triplet Systems and Beyond with Modern Relativistic Density Functional Methods.

An efficient framework for the calculation of paramagnetic NMR (pNMR) shifts within exact two-component (X2C) theory and (current-dependent) density functional theory (DFT) up to the class of local hybrid functionals (LHFs) is presented. Generally, pNMR shifts for systems with more than one unpaired electron depend on the orbital shielding contribution and a temperature-dependent term. The latter includes zero-field splitting (ZFS), hyperfine coupling (HFC), and the g-tensor. For consistency, we calculate these three tensors at the same level of theory, i.e., using scalar-relativistic X2C augmented with spin-orbit perturbation theory. Results for pNMR chemical shifts of transition-metal complexes reveal that this X2C-DFT framework can yield good results for both the shifts and the individual tensor contributions of metallocenes and related systems, especially if the HFC constant is large. For small HFC constants, the relative error is often large, and sometimes the sign may be off. 4d and 5d complexes with more complicated structures demonstrate the limitations of a fully DFT-based approach. Additionally, a Co-based complex with a very large ZFS and pronounced multireference character is not well described. Here, a hybrid DFT-multireference framework is necessary for accurate results. Our results show that X2C is sufficient to describe relativistic effects and computationally cheaper than a fully relativistic approach. Thus, it allows use of large basis sets for converged HFCs. Overall, current-dependent meta-generalized gradient approximations and LHFs show some potential; however, the currently available functionals leave a lot to be desired, and the predictive power is limited.

Current density functional framework for spin-orbit coupling: Extension to periodic systems.

Spin-orbit coupling induces a current density in the ground state, which consequently requires a generalization for meta-generalized gradient approximations. That is, the exchange-correlation energy has to be constructed as an explicit functional of the current density, and a generalized kinetic energy density has to be formed to satisfy theoretical constraints. Herein, we generalize our previously presented formalism of spin-orbit current density functional theory [Holzer et al., J. Chem. Phys. 157, 204102 (2022)] to non-magnetic and magnetic periodic systems of arbitrary dimension. In addition to the ground-state exchange-correlation potential, analytical derivatives such as geometry gradients and stress tensors are implemented. The importance of the current density is assessed for band gaps, lattice constants, magnetic transitions, and Rashba splittings. In the latter, the impact of the current density may be larger than the deviation between different density functional approximations.

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.

Activating Electroluminescence of Charged Naphthalene Diimide Complexes Directly Adsorbed on a Metal Substrate.

Electroluminescence from single molecules adsorbed on a conducting surface imposes conflicting demands for the molecule-electrode coupling. To conduct electrons, the molecular orbitals need to be hybridized with the electrodes. To emit light, they need to be decoupled from the electrodes to prevent fluorescence quenching. Here, we show that fully quenched 2,6-core-substituted naphthalene diimide derivative in a self-assembled monolayer directly deposited on a Au(111) surface can be activated with the tip of a scanning tunneling microscope to decouple the relevant frontier orbitals from the metallic substrate. In this way, individual molecules can be driven from a strongly hybridized state with quenched luminescence to a light-emitting state. The emission performance compares in terms of quantum efficiency, stability, and reproducibility to that of single molecules deposited on thin insulating layers. Quantum chemical calculations suggest that the emitted light originates from the singly charged cationic pair of the molecules.

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.

Natural virtual orbitals for the GW method in the random-phase approximation and beyond.

The increasingly popular GW method is becoming a convenient tool to determine vertical ionization energies in molecular systems. However, depending on the formalism used and the range of orbitals investigated, it may be hampered by a steep computational scaling. To alleviate this issue, correlated natural virtual orbitals (NVOs) based on second-order Møller-Plesset (MP2) and direct MP2 correlation energies are implemented, and the resulting correlated NVOs are tested on GW quasiparticle energies. Test cases include the popular GW variants G0W0 and evGW0 as well as more elaborate vertex corrections. We find that for increasingly larger molecular systems and basis sets, NVOs considerably improve efficiency. Furthermore, we test the performance of the truncated (frozen) NVO ansatz on the GW100 test set. For the latter, it is demonstrated that, using a carefully chosen truncation threshold, NVOs lead to a negligible loss in accuracy while providing speedups of one order of magnitude. Furthermore, we compare the resulting quasiparticle energies to very accurate vertical ionization energies obtained from coupled-cluster theory with singles, doubles, and noniterative triples [CCSD(T)], confirming that the loss in accuracy introduced by truncating the NVOs is negligible compared to the methodical errors in the GW approximation. It is also demonstrated that the choice of basis set impacts results far more than using a suitably truncated NVO space. Therefore, at the same computational expense, more accurate results can be obtained using NVOs. Finally, we provide improved reference CCSD(T) values for the GW100 test set, which have been obtained using the def2-QZVPP basis set.

Exact two-component theory becoming an efficient tool for NMR shieldings and shifts with spin-orbit coupling.

We present a gauge-origin invariant exact two-component (X2C) approach within a modern density functional framework, supporting meta-generalized gradient approximations such as TPSS and range-separated hybrid functionals such as CAM-B3LYP. The complete exchange-correlation kernel is applied, including the direct contribution of the field-dependent basis functions and the reorthonormalization contribution from the perturbed overlap matrix. Additionally, the finite nucleus model is available for the electron-nucleus potential and the vector potential throughout. Efficiency is ensured by the diagonal local approximation to the unitary decoupling transformation in X2C as well as the (multipole-accelerated) resolution of the identity approximation for the Coulomb term (MARI-J, RI-J) and the seminumerical exchange approximation. Errors introduced by these approximations are assessed and found to be clearly negligible. The applicability of our implementation to large-scale calculations is demonstrated for a tin pincer-type system as well as low-valent tin and lead complexes. Here, the calculation of the Sn nuclear magnetic resonance shifts for the pincer-type ligand with about 2400 basis functions requires less than 1 h for hybrid density functionals. Further, the impact of spin-orbit coupling on the nucleus-independent chemical shifts and the corresponding ring currents of all-metal aromatic systems is studied.

Zero-field splitting parameters within exact two-component theory and modern density functional theory using seminumerical integration.

An efficient implementation of zero-field splitting parameters based on the work of Schmitt et al. [J. Chem. Phys. 134, 194113 (2011)] is presented. Seminumerical integration techniques are used for the two-electron spin-dipole contribution and the response equations of the spin-orbit perturbation. The original formulation is further generalized. First, it is extended to meta-generalized gradient approximations and local hybrid functionals. For these functional classes, the response of the paramagnetic current density is considered in the coupled-perturbed Kohn-Sham equations for the spin-orbit perturbation term. Second, the spin-orbit perturbation is formulated within relativistic exact two-component theory and the screened nuclear spin-orbit (SNSO) approximation. The accuracy of the implementation is demonstrated for transition-metal and diatomic main-group compounds. The efficiency is assessed for Mn and Mo complexes. Here, it is found that coarse integration grids for the seminumerical schemes lead to drastic speedups while introducing clearly negligible errors. In addition, the SNSO approximation substantially reduces the computational demands and leads to very similar results as the spin-orbit mean field Ansatz.

A local hybrid exchange functional approximation from first principles.

Local hybrid functionals are a more flexible class of density functional approximations, allowing for a position-dependent admixture of exact exchange. This additional flexibility, however, comes with a more involved mathematical form and a more complicated design. A common denominator for previously constructed local hybrid functionals is the usage of thermochemical benchmark data to construct these functionals. Herein, we design a local hybrid functional without relying on benchmark data. Instead, we construct it in a more ab initio manner, following the principles of modern meta-generalized gradient approximations and considering theoretical constraints. To achieve this, we make use of the density matrix expansion and a local mixing function based on an approximate correlation length. The accuracy of the developed density functional approximation is assessed for thermochemistry, excitation energies, polarizabilities, magnetizabilities, nuclear magnetic resonance (NMR) spin-spin coupling constants, NMR shieldings, and shifts, as well as EPR g-tensors and hyperfine coupling constants. Here, the new exchange functional shows a robust performance and is especially well suited for atomization energies, barrier heights, excitation energies, NMR coupling constants, and EPR properties, whereas it loses some ground for the NMR shifts. Therefore, the designed functional is a major step forward for functionals that have been designed from first principles.

Impact of the current density on paramagnetic NMR properties.

Meta-generalized gradient approximations (meta-GGAs) and local hybrid functionals generally depend on the kinetic energy density τ. For magnetic properties, this necessitates generalizations to ensure gauge invariance. In most implementations, τ is generalized by incorporating the external magnetic field. However, this introduces artifacts in the response of the density matrix and does not satisfy the iso-orbital constraint. Here, we extend previous approaches based on the current density to paramagnetic nuclear magnetic resonance (NMR) shieldings and electron paramagnetic resonance (EPR) g-tensors. The impact is assessed for main-group compounds and transition-metal complexes considering 25 density functional approximations. It is shown that the current density leads to substantial improvements-especially for the popular Minnesota and strongly constrained and appropriately normed (SCAN) functional families. Thus, we strongly recommend to use the current density generalized τ in paramagnetic NMR and EPR calculations with meta-GGAs.

Current density functional framework for spin-orbit coupling.

Relativistic two-component density functional calculations are carried out in a non-collinear formalism to describe spin-orbit interactions, where the exchange-correlation functional is constructed as a generalization of the non-relativistic density functional approximation. Contrary to non-relativistic density functional theory (DFT), spin-orbit coupling, however, leads to a non-vanishing paramagnetic current density. Density functionals depending on the kinetic energy density, such as meta-generalized gradient approximations, should therefore be constructed in the framework of current DFT (CDFT). The latter has previously exclusively been used in the regime of strong magnetic fields. Herein, we present a consistent CDFT approach for relativistic DFT, including spin-orbit coupling. Furthermore, we assess the importance of the current density terms for ground-state energies, excitation energies, nuclear magnetic resonance shielding, and spin-spin coupling constants, as well as hyperfine coupling constants, Δg-shifts, and the nuclear quadrupole interaction tensor in electron paramagnetic resonance (EPR) spectroscopy. The most notable changes are found for EPR properties. The impact of the current-dependent terms rises with the number of unpaired electrons, and consequently, the EPR properties are more sensitive toward CDFT. Considerable changes are observed for the strongly constrained and appropriately normed functionals, as well as the B97M family and TASK. The current density terms are less important when exact exchange is incorporated. At the same time, the current-dependent kernel ensures the stability of response calculations in all cases. We, therefore, strongly recommend to use the framework of CDFT for self-consistent spin-orbit calculations.

Modeling and measuring plasmonic excitations in hollow spherical gold nanoparticles.

We investigate molecular plasmonic excitations sustained in hollow spherical gold nanoparticles using time-dependent density functional theory (TD-DFT). Specifically, we consider Au60 spherical, hollow molecules as a toy model for single-shell plasmonic molecules. To quantify the plasmonic character of the excitations obtained from TD-DFT, the energy-based plasmonicity index is generalized to the framework of DFT, validated on simple systems such as the sodium Na20 chain and the silver Ag20 compound, and subsequently successfully applied to more complex molecules. We also compare the quantum mechanical TD-DFT simulations to those obtained from a classical Mie theory that relies on macroscopic electrodynamics to model the light-matter interaction. This comparison allows us to distinguish those features that can be explained classically from those that require a quantum-mechanical treatment. Finally, a double-shell system obtained by placing a C60 buckyball inside the hollow spherical gold particle is further considered. It is found that the double-shell, while increasing the overall plasmonic character of the excitations, leads to significantly lowered absorption cross sections.

18-Crown-6 Coordinated Metal Halides with Bright Luminescence and Nonlinear Optical Effects.

The crown-ether coordination compounds ZnX2(18-crown-6), EuX2(18-crown-6) (X: Cl, Br, I), MnI2(18-crown-6), Mn3Cl6(18-crown-6)2, Mn3I6(18-crown-6)2, and Mn2I4(18-crown-6) are obtained by ionic-liquid-based synthesis. Whereas MX2(18-crown-6) (M: Zn, Eu) show conventional structural motives, Mn3Cl6(18-crown-6)2, Mn3I6(18-crown-6)2, and Mn2I4(18-crown-6) exhibit unusual single MnX4 tetrahedra coordinated to the crown-ether complex. Surprisingly, some compounds show outstanding photoluminescence. Thus, rare Zn2+-based luminescence is observed and unexpectedly efficient for ZnI2(18-crown-6) with a quantum yield of 54%. Unprecedented quantum yields are also observed for Mn3I6(18-crown-6)2, EuBr2(18-crown-6), and EuI2(18-crown-6) with values of 98, 72, and 82%, respectively, which can be rationalized based on the specific structural features. Most remarkable, however, is Mn2I4(18-crown-6). Its specific structural features with finite sensitizer-activator couples result in an extremely strong emission with an outstanding quantum yield of 100%. Consistent with its structural features, moreover, anisotropic angle-dependent emission under polarized light and nonlinear optical (NLO) effects occur, including second-harmonic generation (SHG). The title compounds and their optical properties are characterized by single-crystal structure analysis, X-ray powder diffraction, chemical analysis, density functional theory (DFT) calculations, and advanced spectroscopic methods.

Boosting Light Emission from Single Hydrogen Phthalocyanine Molecules by Charging.

Interest in electroluminescence of single molecules is stimulated by the prospect of possible applications in novel light emitting devices. Recent studies provide valuable insights into the mechanisms leading to single molecule electroluminescence. Concrete information on how to boost the intensity of the emitted light, however, is rare. By combining scanning tunnelling microscopy (STM) and quantum chemical calculations, we show that the light emission efficiencies of an individual hydrogen-phthalocyanine molecule can be increased by a factor of ≈19 upon charging. This boost in intensity can be explained by the development of a vertical dipole moment normal to the substrate facilitating out-coupling of the local excitation to the far field. As this effect is not related to the specific nature of hydrogen-phthalocyanine, it opens up a general way to increase light emission from molecular junctions.

Synthesis of New Donor-Substituted Biphenyls: Pre-ligands for Highly Luminescent (C

We herein report on new synthetic strategies for the preparation of pyridine and imidazole substituted 2,2'-dihalo biphenyls. These structures are pre-ligands suitable for the preparation of respective stannoles. The latter can successfully be transmetalated to K[AuCl4 ] forming non-palindromic [(C^C^D)AuIII ] pincer complexes featuring a lateral pyridine (D=N) or N-heterocyclic carbene (NHC, D=C') donor. The latter is the first report on a pincer complex with two formally anionic sp2 and one carbenic carbon donor. The [(C^C^D)AuIII ] complexes show intense phosphorescence in solution at room temperature. We discuss the developed multistep strategy and touch upon synthetic challenges. The prepared complexes have been fully characterized including X-ray diffraction analysis. The gold(III) complexes' photophysical properties have been investigated by absorption and emission spectroscopy as well as quantum chemical calculations on the quasi-relativistic two-component TD-DFT and GW/Bethe-Salpeter level including spin-orbit coupling. Thus, we shed light on the electronic influence of the non-palindromic pincer ligand and reveal non-radiative relaxation pathways of the different ligands employed.

An improved seminumerical Coulomb and exchange algorithm for properties and excited states in modern density functional theory.

A seminumerical algorithm capable of performing large-scale (time-dependent) density functional theory (TD-DFT) calculations to extract excitation energies and other ground-state and excited-state properties is outlined. The algorithm uses seminumerical integral techniques for evaluating Coulomb and exchange parts for a set of density matrices as occurring in standard TD-DFT or similar methods for the evaluation of vibrational frequencies. A suitable optimized de-aliasing procedure is introduced. The latter does not depend on further auxiliary quantities and retains the symmetry of a given density matrix. The algorithm is self-contained and applicable to any orbital basis set available without the need for further auxiliary basis sets or optimized de-aliasing grids. Relativistic two-component excited-state TD-DFT calculations are reported for the first time using the developed seminumerical algorithm for standard and local hybrid density functional approximations. Errors are compared with the widely used "resolution of the identity" (RI) approximations for Coulomb (RI-J) and exchange integrals (RI-K). The fully seminumerical algorithm does not exhibit an enlarged error for standard DFT functionals compared to the RI approximation. For the more involved local hybrid functionals and within strong external fields, accuracy is even considerably improved.

The first microsolvation step for furans: New experiments and benchmarking strategies.

The site-specific first microsolvation step of furan and some of its derivatives with methanol is explored to benchmark the ability of quantum-chemical methods to describe the structure, energetics, and vibrational spectrum at low temperature. Infrared and microwave spectra in supersonic jet expansions are used to quantify the docking preference and some relevant quantum states of the model complexes. Microwave spectroscopy strictly rules out in-plane docking of methanol as opposed to the top coordination of the aromatic ring. Contrasting comparison strategies, which emphasize either the experimental or the theoretical input, are explored. Within the harmonic approximation, only a few composite computational approaches are able to achieve a satisfactory performance. Deuteration experiments suggest that the harmonic treatment itself is largely justified for the zero-point energy, likely and by design due to the systematic cancellation of important anharmonic contributions between the docking variants. Therefore, discrepancies between experiment and theory for the isomer abundance are tentatively assigned to electronic structure deficiencies, but uncertainties remain on the nuclear dynamics side. Attempts to include anharmonic contributions indicate that for systems of this size, a uniform treatment of anharmonicity with systematically improved performance is not yet in sight.

TURBOMOLE: Modular program suite for ab initio quantum-chemical and condensed-matter simulations.

TURBOMOLE is a collaborative, multi-national software development project aiming to provide highly efficient and stable computational tools for quantum chemical simulations of molecules, clusters, periodic systems, and solutions. The TURBOMOLE software suite is optimized for widely available, inexpensive, and resource-efficient hardware such as multi-core workstations and small computer clusters. TURBOMOLE specializes in electronic structure methods with outstanding accuracy-cost ratio, such as density functional theory including local hybrids and the random phase approximation (RPA), GW-Bethe-Salpeter methods, second-order Møller-Plesset theory, and explicitly correlated coupled-cluster methods. TURBOMOLE is based on Gaussian basis sets and has been pivotal for the development of many fast and low-scaling algorithms in the past three decades, such as integral-direct methods, fast multipole methods, the resolution-of-the-identity approximation, imaginary frequency integration, Laplace transform, and pair natural orbital methods. This review focuses on recent additions to TURBOMOLE's functionality, including excited-state methods, RPA and Green's function methods, relativistic approaches, high-order molecular properties, solvation effects, and periodic systems. A variety of illustrative applications along with accuracy and timing data are discussed. Moreover, available interfaces to users as well as other software are summarized. TURBOMOLE's current licensing, distribution, and support model are discussed, and an overview of TURBOMOLE's development workflow is provided. Challenges such as communication and outreach, software infrastructure, and funding are highlighted.