TURBOMOLE: Today and Tomorrow

Which Options Exist for NISQ-Friendly Linear Response Formulations?

Linear response (LR) theory is a powerful tool in classic quantum chemistry crucial to understanding photoinduced processes in chemistry and biology. However, performing simulations for large systems and in the case of strong electron correlation remains challenging. Quantum computers are poised to facilitate the simulation of such systems, and recently, a quantum linear response formulation (qLR) was introduced [Kumar et al., J. Chem. Theory Comput. 2023, 19, 9136-9150]. To apply qLR to near-term quantum computers beyond a minimal basis set, we here introduce a resource-efficient qLR theory, using a truncated active-space version of the multiconfigurational self-consistent field LR ansatz. Therein, we investigate eight different near-term qLR formalisms that utilize novel operator transformations that allow the qLR equations to be performed on near-term hardware. Simulating excited state potential energy curves and absorption spectra for various test cases, we identify two promising candidates, dubbed "proj LRSD" and "all-proj LRSD".

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.

Transferable Anisotropic Mie Potential Force Field for Alkanediols

The development of force fields for polyfunctional molecules, such as alkanediols, requires a careful account of different average intramolecular conformations for gas states compared to dense liquid states, where intra- and intermolecular hydrogen bonds compete. In the present work, the transferable anisotropic Mie (TAMie) potential is extended to 1,n-alkanediols. Using the convention that intramolecular nonbonded interactions up to and including the third neighbor are excluded, all force field parameters developed previously for 1-alcohols were transferred to 1,5-pentanediol and beyond, with good agreement with experimental phase equilibrium data. To obtain trans-gauche ratios of 1,2-ethanediol and 1,3-propanediol that are consistent with experimental results, the propensities for intra- and intermolecular hydrogen bonds had to be balanced. This was achieved by parameterizing the intramolecular dihedral energy functions governing the O-C-C-O and O-C-C-C angles while intramolecular charge-charge interactions were active. All partial charges belonging to a functional group are collected in a charge group and all interactions among two charge groups are evaluated even if they are separated by less than three bonds. With this approach, it is possible to apply the nonbonded parameters from 1-alcohols to alkanediols without further refinement. The agreement with experimental phase equilibrium and shear viscosity data is of similar quality as for the 1-alcohols and the trans-gauche ratio agrees with literature results from spectroscopic measurements and ab initio calculations.

Optical study of Te8 ring clusters: comparison with density functional theory and a step towards materials design using nanoporous zeolite space

The Te8 ring molecule (cluster) is poorly investigated due to the lack of experimental data. Here, we report an experimental and theoretical study of a regular array of oriented Te8 rings formed in the ∼1.14 nm diameter cavities of zeolite LTA, which are arranged in a cubic lattice with a spacing of ∼1.2 nm. Single crystals of LTA with encapsulated tellurium (LTA-Te) were studied using Raman spectroscopy (RS) and optical absorption spectroscopy (OAS). The experimental LTA-Te spectra were found to be in agreement with those calculated using density functional theory (PBE0 hybrid functional and def2-TZVP basis sets) for the crown-shaped Te8 ring molecule with D4d symmetry. Using polarization-orientation RS, we show that the Te8 rings are oriented by their major axes along the 4-fold axes of cubic LTA. We also show that the site symmetry of Te8 in LTA-Te is lower than D4d. Te8 bond-bending modes are well described in the harmonic approximation, while bond-stretching modes are mixed due to the reduced ring symmetry and, probably, anharmonicity. Importantly, OAS data of LTA-Te display dependence on the Te8 concentration, implying the interaction of the rings from neighbouring LTA cavities with the generation of the valence and conduction electron bands of such a cluster crystal.

Interoperable workflows by exchanging grid-based data between quantum-chemical program packages

Quantum-chemical subsystem and embedding methods require complex workflows that may involve multiple quantum-chemical program packages. Moreover, such workflows require the exchange of voluminous data that go beyond simple quantities, such as molecular structures and energies. Here, we describe our approach for addressing this interoperability challenge by exchanging electron densities and embedding potentials as grid-based data. We describe the approach that we have implemented to this end in a dedicated code, PyEmbed, currently part of a Python scripting framework. We discuss how it has facilitated the development of quantum-chemical subsystem and embedding methods and highlight several applications that have been enabled by PyEmbed, including wave-function theory (WFT) in density-functional theory (DFT) embedding schemes mixing non-relativistic and relativistic electronic structure methods, real-time time-dependent DFT-in-DFT approaches, the density-based many-body expansion, and workflows including real-space data analysis and visualization. Our approach demonstrates, in particular, the merits of exchanging (complex) grid-based data and, in general, the potential of modular software development in quantum chemistry, which hinges upon libraries that facilitate interoperability.

Consistent Analytical Second Derivatives of the Kohn-Sham DFT Energy in the Framework of the Conductor-Like Screening Model through Gaussian Charge Distributions

The use of implicit solvation models such as the conductor-like screening model (COSMO) in quantum chemical calculations is very common, as both a rough estimate of solvation effects as well as a general tool for stabilizing ionic molecular structures. In order to generate a smooth potential energy surface as well as consistent gradients, it is necessary to apply the Gaussian charge model (GCM) for the COSMO charges. This work introduces an efficient implementation for consistent analytical second derivatives of the electronic energy with COSMO-GCM in the framework of the Kohn-Sham density functional theory. This is used to investigate the infrared spectroscopy of amino acids in aqueous solution, where the impact of pH on the molecular structure and vibrational spectra is examined. Furthermore, the structure and stability of selected all-metal aromatic cluster ions are assessed.

The effect of particle size on the optical and electronic properties of hydrogenated silicon nanoparticles

We use a combination of many-body perturbation theory and time-dependent density functional theory to study the optical and electronic properties of hydrogen terminated silicon nanoparticles. We predict that the lowest excited states of these silicon nanoparticles are excitonic in character and that the corresponding excitons are completely delocalised over the volume of the particle. The size of the excitons is predicted to increase proportionally with the particle size. Conversely, we predict that the fundamental gap, the optical gap, and the exciton binding energy increase with decreasing particle size. The exciton binding energy is predicted to counter-act the variation in the fundamental gap and hence to reduce the variation of the optical gap with particle size. The variation in the exciton binding energy itself is probably caused by a reduction in the dielectric screening with decreasing particle size. The intensity of the excited state corresponding to the optical gap and other low energy excitations are predicted to increase with decreasing particle size. We explain this increase in terms of the 'band structure' becoming smeared out in reciprocal space with decreasing particle size, increasing the 'overlap' between the occupied and unoccupied quasiparticle states and thus, the oscillator strength. Fourier transforms of the lowest excitons show that they inherit the periodicity of the frontier quasiparticle states. This, combined with the delocalisation of the exciton and the large exciton binding energy, means that the excitons in silicon nanoparticles combine aspects of Wannier-Mott, delocalisation and effect of periodicity of the underlying structure, and Frenkel, large exciton binding energy, excitons.

Improved CPS and CBS Extrapolation of PNO-CCSD(T) Energies: The MOBH35 and ISOL24 Data Sets

Computation of heats of reaction of large molecules is now feasible using the domain-based pair natural orbital (PNO)-coupled-cluster singles, doubles, and perturbative triples [CCSD(T)] theory. However, to obtain agreement within 1 kcal/mol of experiment, it is necessary to eliminate basis set incompleteness error, which comprises both the AO basis set error and the PNO truncation error. Our investigation into the convergence to the canonical limit of PNO-CCSD(T) energies with the PNO truncation threshold T shows that errors follow the model E ( T ) = E + A T 1 / 2 . Therefore, PNO truncation errors can be eliminated using a simple two-point CPS extrapolation to the canonical limit so that subsequent CBS extrapolation is not limited by the residual PNO truncation error. Using the ISOL24 and MOBH35 data sets, we find that PNO truncation errors are larger for molecules with significant static correlation and that it is necessary to use very tight thresholds of T = 10 - 8 to ensure that errors do not exceed 1 kcal/mol. We present a lower-cost extrapolation scheme that uses information from small basis sets to estimate the PNO truncation errors for larger basis sets. In this way, the canonical limit of CCSD(T) calculations on sizable molecules with large basis sets can be reliably estimated in a practical way. Using this approach, we report near complete basis set (CBS)-CCSD(T) reaction energies for the full ISOL24 and MOBH35 data sets.

Role of Singles Amplitudes in ADC(2) and CC2 for Low-Lying Electronically Excited States

The closely related second-order methods CC2 and ADC(2) usually perform very similarly for single excitations of organic molecules. However, as rationalized in this work, significant deviations between these two methods can arise if the ground state and a low-lying singly excited state arise from a strong coupling between their leading configurations. Such a configuration mixing is partially accounted for in CC2 through the ground-state singles amplitudes but is omitted in ADC(2). This can cause unusual deviations between the results obtained with these methods. In this work, we study how severe this effect can become at the example of two solvatochromic dyes: the negatively solvatochromic betaine dye N1-tBu and the positively solvatochromic bithiophene P1. These two dyes allow one to study the limits of both small and somewhat larger excitation energies and configuration mixing by tuning the S0 → S1 transition energy through the polarity of the environment. Higher-level calculations at the CC3 level provide information on the accuracy of ADC(2) and CC2 in these cases. The most extreme deviation between ADC(2) and CC2 is found for N1-tBu in vacuum, where the ADC(2) result is 0.45 eV below that of CC2. In this case, the methodical error of CC2 with respect to CC3 is only 0.05 eV. With increasing excitation energy in polar solvents, the CC2-ADC(2) deviation decreases and reaches a value of only 0.15 eV. For P1, which has larger excitation energies, these effects are reversed due to the opposite solvatochromism but also smaller in magnitude: the deviation increases from 0.08 eV in vacuum to 0.16 eV in the so-called conductor limit of the continuum solvation model. Although for these two dyes larger deviations are observed for smaller excitation energies, the extent of configuration mixing does not generally correlate with only the size of excitation energy. For example, s-triazine (0.15 eV), formamide (0.19 eV), and formaldehyde (0.23 eV) also show large deviations between CC2 and ADC(2) despite their much higher excitation energies compared to those of N1-tBu and P1.

Implementation and First Evaluation of Strong-Correlation-Corrected Local Hybrid Functionals for the Calculation of NMR Shieldings and Shifts

Local hybrid functionals containing strong-correlation factors (scLHs) and range-separated local hybrids (RSLHs) have been integrated into an efficient coupled-perturbed Kohn-Sham implementation for the calculation of nuclear shielding constants. Several scLHs and the ωLH22t RSLH have then been evaluated for the first time for the extended NS372 benchmark set of main-group shieldings and shifts and the TM70 benchmark of 3d transition-metal shifts. The effects of the strong-correlation corrections have been analyzed with respect to the spatial distribution of the sc-factors, which locally diminish exact-exchange admixture at certain regions in a molecule. The scLH22t, scLH23t-mBR, and scLH23t-mBR-P functionals, which contain a "damped" strong-correlation factor to retain the excellent performance of the underlying LH20t functional for weakly correlated situations, tend to make smaller corrections to shieldings and shifts than the "undamped" scLH22ta functional. While the latter functional can also deteriorate agreement with the reference data in certain weakly correlated cases, it provides overall better performance, in particular for systems where static correlation is appreciable. This pertains only to a minority of systems in the NS372 main-group test set but to many more systems in the TM70 transition-metal test set, in particular for high-oxidation-state complexes, e.g., Cr(+VI) complexes and other systems with stretched bonds. Another undamped scLH, the simpler LDA-based scLH21ct-SVWN-m, also tends to provide significant improvements in many cases. The differences between the functionals and species can be rationalized on the basis of one-dimensional plots of the strong-correlation factors, augmented by isosurface plots of the fractional orbital density (FOD). Position-dependent exact-exchange admixture is thus shown to provide substantial flexibility in treating response properties like NMR shifts for both weakly and strongly correlated systems.

Insights into localization, energy ordering, and substituent effect in excited states of azobenzenes from coupled cluster calculations of nuclear spin-induced circular dichroism

Nuclear spin-induced circular dichroism (NSCD) is a molecular effect of differential absorption of left- and right-circularly polarized light due to nuclear spins in the molecule. In this work, new tools for its calculation are presented. Specifically, analytic expressions for the computation of the K term of NSCD have been derived and implemented for the second-order coupled cluster singles and doubles (CC2) model. NSCD results obtained thereby for three derivatives of azobenzenes have been compared with results from time-dependent density functional theory (TD-DFT). The complementary information that could be obtained from NSCD measurements compared to NMR for these three species is discussed.

Spin-Symmetry Breaking and Hyperfine Couplings in Transition-Metal Complexes Revisited Using Density Functionals Based on the Exact-Exchange Energy Density

A small set of mononuclear manganese complexes evaluated previously for their Mn hyperfine couplings (HFCs) has been analyzed using density functionals based on the exact-exchange energy density─in particular, the spin symmetry breaking (SSB) found previously when using hybrid functionals. Employing various strong-correlation corrected local hybrids (scLHs) and strong-correlation corrected range-separated local hybrids (scRSLHs) with or without additional corrections to their local mixing functions (LMFs) to mitigate delocalization errors (DE), the SSB and the associated dipolar HFCs of [Mn(CN)4]2-, MnO3, [Mn(CN)4N]-, and [Mn(CN)5NO]2- (the latter with cluster embedding) have been examined. Both strong-correlation (sc)-correction and DE-correction terms help to diminish SSB and correct the dipolar HFCs. The DE corrections are more effective, and the effects of the sc corrections depend on their damping factors. Interestingly, the DE-corrections reduce valence-shell spin polarization (VSSP) and thus SSB by locally enhancing exact-exchange (EXX) admixture near the metal center and thereby diminishing spin-density delocalization onto the ligand atoms. In contrast, sc corrections diminish EXX admixture locally, mostly on specific ligand atoms. This then reduces VSSP and SSB as well. The performance of scLHs and scRSLHs for the isotropic Mn HFCs has also been analyzed, with particular attention to core-shell spin-polarization contributions. Further sc-corrected functionals, such as the KP16/B13 construction and the DM21 deep-neural-network functional, have been examined.

Understanding Aromaticity in [5]Helicene-Bridged Cyclophanes: A Comprehensive Study

This study explores the aromaticity of doubly [5]helicene-bridged (1,4)cyclophane and triply [5]helicene-bridged (1,3,5)cyclophane via calculations of the magnetic response and of electronic aromaticity indices. The primary objective is to assess the π-electron delocalization to determine whether they sustain global ring currents associated with π aromaticity. The molecules show local ring currents in the presence of an external magnetic field. The ring currents flow diatropically in the stacked six-membered rings and in the helicene arms. However, these π currents are not interconnected due to the discontinuity of the π delocalization at the C-C single bonds connecting the central six-membered rings to the helicene arms. Electronic indices suggest that the helicene-arm systems have significantly smaller electron delocalization than benzene. The reduction in the delocalization does not compromise their ability to exhibit ring currents in the presence of an external magnetic field. The analysis provides further evidence that the magnetic criteria yield a different degree of aromaticity for the helicene arms than obtained in the calculation of the electronic aromaticity indices. However, both approaches confirm that the studied molecules are not globally aromatic.

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.

Synthesis and Properties of Azahomocorannulenyl Cations and Radicals

Cycloheptatrienyl (tropyl) molecules are representative non-alternant hydrocarbons that offer interesting chemistry because of their unique structures and properties. However, there have been a limited number of polycyclic aromatic tropyl cations and radicals reported in the literature. Herein, we report the synthesis of a series of azahomocorannulene derivatives, where the key reactions are a 1,3-dipolar cycloaddition of polycyclic aromatic azomethine ylides with dibenzotropone and a subsequent palladium-catalyzed cyclization. X-ray diffraction analysis revealed that the obtained azahomocorannulenyl cation and radical adopt planar structures and exhibit unique packing structures. Their electronic and optical properties were investigated experimentally and theoretically to reveal their aromatic character.

Continuous-wave cavity ringdown for high-sensitivity polarimetry and magnetometry measurements

We report the development of a novel variant of cavity ringdown polarimetry using a continuous-wave laser operating at 532 nm for highly precise chiroptical activity and magnetometry measurements. The key methodology of the apparatus relies upon the external modulation of the laser frequency at the frequency splitting between non-degenerate left- and right-circularly polarized cavity modes. The method is demonstrated by the evaluation of the Verdet constants of crystalline CeF3 and fused silica, in addition to the observation of gas- and solution-phase optical rotations of selected chiral molecules. Specifically, optical rotations of (i) vapors of α-pinene and R-(+)-limonene, (ii) mutarotating D-glucose in water, and (iii) acidified L-histidine solutions are determined. The detection sensitivities for the gas- and solution-phase chiral activity measurements are ∼30 and ∼120μdeg over a 30 s detection period per cavity round trip pass, respectively. Furthermore, the measured optical rotations for R-(+)-limonene are compared with computations performed using the TURBOMOLE quantum chemistry package. The experimentally observed optically rotatory dispersion of this cyclic monoterpene was thus rationalized via a consideration of its room temperature conformer distribution as determined by the aforementioned single-point energy calculations.

Synthesis of Crystallographically Characterizable Bis(cyclopentadienyl) Sc(II) Complexes: (C5H2tBu3)2Sc and {[C5H3(SiMe3)2]2ScI}1

The synthesis of previously unknown bis(cyclopentadienyl) complexes of the first transition metal, i.e., Sc(II) scandocene complexes, has been investigated using C5H2(tBu)3 (Cpttt), C5Me5 (Cp*), and C5H3(SiMe3)2 (Cp″) ligands. Cpttt2ScI, 1, formed from ScI3 and KCpttt, can be reduced with potassium graphite (KC8) in hexanes to generate dark-red crystals of the first crystallographically characterizable bis(cyclopentadienyl) scandium(II) complex, Cpttt2Sc, 2. Complex 2 has a 170.6° (ring centroid)-Sc-(ring centroid) angle and exhibits an eight-line EPR spectrum characteristic of Sc(II) with Aiso = 82.6 MHz (29.6 G). It sublimes at 200 °C at 10-4 Torr and has a melting point of 268-271 °C. Reductions of Cp*2ScI and Cp″2ScI under analogous conditions in hexanes did not provide new Sc(II) complexes, and reduction of Cp*2ScI in benzene formed the Sc(III) phenyl complex, Cp*2Sc(C6H5), 3, by C-H bond activation. However, in Et2O and toluene, reduction of Cp*2ScI at -78 °C gives a dark-red solution, 4, which displays an eight-line EPR pattern like that of 1, but it did not provide thermally stable crystals. Reduction of Cp″2ScI, in THF or Et2O at -35 °C in the presence of 2.2.2-cryptand, yields the green Sc(II) metallocene iodide complex, [K(crypt)][Cp″2ScI], 5, which was identified by X-ray crystallography and EPR spectroscopy and is thermally unstable. The analogous reaction of Cp*2ScI with KC8 and 18-crown-6 in Et2O gave the ligand redistribution product, [Cp*2Sc(18-crown-6-κ2O,O')][Cp*2ScI2], 6, as the only crystalline product. Density functional theory calculations on the electronic structure of these compounds are reported in addition to a steric analysis using the Guzei method.

A Multi-Scale Approach to Simulate the Nonlinear Optical Response of Molecular Nanomaterials

Nonlinear optics is essential for many recent photonic technologies. Here, a novel multi-scale approach is introduced to simulate the nonlinear optical response of molecular nanomaterials combining ab initio quantum-chemical and classical Maxwell-scattering computations. In this approach, the first hyperpolarizability tensor is computed with time-dependent density-functional theory and incorporated into a multi-scattering formalism that considers the optical interaction between neighboring molecules. Such incorporation is achieved by a novel object: the Hyper-Transition(T)-matrix. With this object at hand, the nonlinear optical response from single molecules and also from entire photonic devices can be computed, including the full tensorial and dispersive nature of the optical response of the molecules, as well as the optical interaction between different molecules as, for example, in the lattice of a molecular crystal. To demonstrate the applicability of the novel approach, the generation of a second-harmonic signal from a thin film of an Urea molecular crystal is computed and compared to more traditional simulations. Furthermore, an optical cavity is designed, which enhances the second-harmonic response of the molecular film up to more than two orders of magnitude. This approach is highly versatile and accurate and can be the working horse for the future exploration of nonlinear photonic molecular materials in structured photonic environments.

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.

A DFT perspective on organometallic lanthanide chemistry

Computational studies of the coordination chemistry and bonding of lanthanides have grown in recent decades as the need for understanding the distinct physical, optical, and magnetic properties of these compounds increased. Density functional theory (DFT) methods offer a favorable balance of computational cost and accuracy in lanthanide chemistry and have helped to advance the discovery of novel oxidation states and electronic configurations. This Frontier article examines the scope and limitations of DFT in interpreting structural and spectroscopic data of low-valent lanthanide complexes, elucidating periodic trends, and predicting their properties and reactivity, presented through selected examples.

Adiabatic connection interaction strength interpolation method made accurate for the uniform electron gas

The adiabatic connection interaction strength interpolation (ISI)-like method provides a high-level expression for the correlation energy, being, in principle, exact not only in the weak-interaction limit, where it recovers the second-order Görling-Levy perturbation term, but also in the strong-interaction limit that is described by the strictly correlated electron approach. In this work, we construct a genISI functional made accurate for the uniform electron gas, a solid-state physics paradigm that is a very difficult test for ISI-like correlation functionals. We assess the genISI functional for various jellium spheres with the number of electrons Z ≤ 912 and for the non-relativistic noble atoms with Z ≤ 290. For the jellium clusters, the genISI is remarkably accurate, while for the noble atoms, it shows a good performance, similar to other ISI-like methods. Then, the genISI functional can open the path using the ISI-like method in solid-state calculations.

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.

Koopmans' theorem for acidic protons

The famous Brønsted acidity, which is relevant in many areas of experimental and synthetic chemistry, but also in biochemistry and other areas, is investigated from a new perspective. Nuclear electronic orbital methods, which explicitly account for the quantum character of selected protons, are applied. The resulting orbital energies of the proton wavefunction are interpreted and related to enthalpies of deprotonation and acid strength in analogy to the Koopmans' theorem for electrons. For a set of organic acids, we observe a correlation which indicates the validity of such a NEO-Koopmans' approach and opens up new opportunities for the computational investigation of more complex acidic systems.

Range-Separated Local Hybrid Functionals with Small Fractional-Charge and Fractional-Spin Errors: Escaping the Zero-Sum Game of DFT Functionals

Extending recent developments on strong-correlation (sc) corrections to local hybrid functionals to the recent accurate ωLH22t range-separated local hybrid, a series of highly flexible strong-correlation-corrected range-separated local hybrids (scRSLHs) has been constructed and evaluated. This has required the position-dependent reduction of both short- and long-range exact-exchange admixtures in regions of space characterized by strong static correlations. Using damping procedures provides scRSLHs that retain largely the excellent performance of ωLH22t for weakly correlated situations and, in particular, for accurate quasiparticle energies of a wide variety of systems while reducing dramatically static-correlation errors, e.g., in stretched-bond situations. An additional correction to the local mixing function to reduce delocalization errors in abnormal open-shell situations provides further improvements in thermochemical and kinetic parameters, making scRSLH functionals such as ωLH23tdE or ωLH23tdP promising tools for complex molecular or condensed-phase systems, where low fractional-charge and fractional-spin errors are simultaneously important. The proposed rung 4 functionals thereby largely escape the usual zero-sum game between these two quantities and are expected to open new areas of accurate computations by Kohn-Sham DFT. At the same time, they require essentially no extra computational effort over the underlying ωLH22t functional, which means that their use is only moderately more demanding than that of global, local, or range-separated hybrid functionals.

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.

Learning from the 4-(dimethylamino)benzonitrile twist: Two-parameter range-separated local hybrid functional with high accuracy for triplet and charge-transfer excitations

The recent ωLH22t range-separated local hybrid (RSLH) is shown to provide outstanding accuracy for the notorious benchmark problem of the two lowest excited-state potential energy curves for the amino group twist in 4-(dimethylamino)benzonitrile (DMABN). However, the design of ωLH22t as a general-purpose functional resulted in less convincing performance for triplet excitations, which is an important advantage of previous LHs. Furthermore, ωLH22t uses 8 empirical parameters to achieve broad accuracy. In this work, the RSLH ωLH23ct-sir is constructed with minimal empiricism by optimizing its local mixing function prefactor and range-separation parameter for only 8 excitation energies. ωLH23ct-sir maintains the excellent performance of ωLH22t for the DMABN twist and charge-transfer benchmarks but significantly improves the errors for triplet excitation energies (0.17 vs 0.24 eV). Additional test calculations for the AE6BH6 thermochemistry test set and large dipole moment and static polarizability test sets confirm that the focus on excitation energies in the optimization of ωLH23ct-sir has not caused any dramatic errors for ground-state properties. Although ωLH23ct-sir cannot replace ωLH22t as a general-purpose functional, it is preferable for problems requiring a universally good description of localized and charge-transfer excitations of both singlet and triplet multiplicity. Current limitations on the application of ωLH23ct-sir and other RSLHs to the study of singlet-triplet gaps of emitters for thermally activated delayed fluorescence are discussed. This work also includes the first systematic analysis of the influence of the local mixing function prefactor and the range-separation parameter in an RSLH on different types of excitations.

Efficient Protocol for Computing MCD Spectra in a Broad Frequency Range Combining Resonant and Damped CC2 Quadratic Response Theory

Coupled cluster response theory offers a path to high-accuracy calculations of spectroscopic properties, such as magnetic circular dichroism (MCD). However, divergence or slow convergence issues are often encountered for electronic transitions in high-energy regions with a high density of states. This is here addressed for MCD by an implementation of damped quadratic response theory for resolution-of-identity coupled cluster singles-and-approximate-doubles (RI-CC2), along with an implementation of the MCD A term from resonant response theory. Combined, damped and resonant response theory calculations provide an efficient strategy to obtain MCD spectra over a broad frequency range and for systems that include highly symmetric molecules with degenerate excited states. The protocol is illustrated by application to zinc tetrabenzoporphyrin in the energy region of 2-8 eV and comparison to experimental data. Timings are reported for the resonant and damped approaches, showing that a greater part of the calculation time is consumed by the construction of the building blocks for the final MCD ellipticity. A recommendation on how to use the procedure is outlined.

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.