Segmented Contracted Error-Consistent Basis Sets of Quadruple-ζ Valence Quality for One- and Two-Component Relativistic All-Electron Calculations

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.

Chelate Complexes of 3d Transition Metal Ions─A Challenge for Electronic-Structure Methods?

Different electronic-structure methods were assessed for their ability to predict two important properties of the industrially relevant chelating agent nitrilotriacetic acid (NTA): its selectivity with respect to six different first-row transition metal ions and the spin-state energetics of its complex with Fe(III). The investigated methods encompassed density functional theory (DFT), the random phase approximation (RPA), coupled cluster (CC) theory, and the auxiliary-field quantum Monte Carlo (AFQMC) method, as well as the complete active space self-consistent field (CASSCF) method and the respective on-top methods: second-order N-electron valence state perturbation theory (NEVPT2) and multiconfiguration pair-density functional theory (MC-PDFT). Different strategies for selecting active spaces were explored, and the density matrix renormalization group (DMRG) approach was used to solve the largest active spaces. Despite somewhat ambiguous multi-reference diagnostics, most methods gave relatively good agreement with experimental data for the chemical reactions connected to the selectivity, which only involved transition-metal complexes in their high-spin state. CC methods yielded the highest accuracy followed by range-separated DFT and AFQMC. We discussed in detail that even higher accuracies can be obtained with NEVPT2, under the prerequisite that consistent active spaces along the entire chemical reaction can be selected, which was not the case for reactions involving Fe(III). A bigger challenge for electronic-structure methods was the prediction of the spin-state energetics, which additionally involved lower spin states that exhibited larger multi-reference diagnostics. Conceptually different, typically accurate methods ranging from CC theory via DMRG-NEVPT2 in combination with large active spaces to AFQMC agreed well that the high-spin state is energetically significantly favored over the other spin states. This was in contrast to most DFT functionals and RPA which yielded a smaller stabilization and some common DFT functionals and MC-PDFT even predicting the low-spin state to be energetically most favorable.

The p-block challenge: assessing quantum chemistry methods for inorganic heterocycle dimerizations

The elements of the p-block of the periodic table are of high interest in various chemical and technical applications like frustrated Lewis-pairs (FLP) or opto-electronics. However, high-quality benchmark data to assess approximate density functional theory (DFT) for their theoretical description are sparse. In this work, we present a benchmark set of 604 dimerization energies of 302 "inorganic benzenes" composed of all non-carbon p-block elements of main groups III to VI up to polonium. This so-called IHD302 test set comprises two classes of structures formed by covalent bonding and by weaker donor-acceptor (WDA) interactions, respectively. Generating reliable reference data with ab initio methods is challenging due to large electron correlation contributions, core-valence correlation effects, and especially the slow basis set convergence. To compute reference values for these dimerization reactions, after thorough testing, we applied a computational protocol using state-of-the-art explicitly correlated local coupled cluster theory termed PNO-LCCSD(T)-F12/cc-VTZ-PP-F12(corr.). It includes a basis set correction at the PNO-LMP2-F12/aug-cc-pwCVTZ level. Based on these reference data, we assess 26 DFT methods in combination with three different dispersion corrections and the def2-QZVPP basis set, five composite DFT approaches, and five semi-empirical quantum mechanical methods. For the covalent dimerizations, the r2SCAN-D4 meta-GGA, the r2SCAN0-D4 and ωB97M-V hybrids, and the revDSD-PBEP86-D4 double-hybrid functional are found to be the best-performing methods among the evaluated functionals of the respective class. However, since def2 basis sets for the 4th period are not associated to relativistic pseudo-potentials, we obtained significant errors in the covalent dimerization energies (up to 6 kcal mol-1) for molecules containing p-block elements of the 4th period. Significant improvements were achieved for systems containing 4th row elements by using ECP10MDF pseudopotentials along with re-contracted aug-cc-pVQZ-PP-KS basis sets introduced in this work with the contraction coefficients taken from atomic DFT (PBE0) calculations. Overall, the IHD302 set represents a challenge to contemporary quantum chemical methods. This is due to a large number of spatially close p-element bonds which are underrepresented in other benchmark sets, and the partial covalent bonding character for the WDA interactions. The IHD302 set may be helpful to develop more robust and transferable approximate quantum chemical methods in the future.

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.

Prediction of 19F NMR chemical shifts for organic compounds with ORCA

Accurate assignment of 19F NMR has long been a challenge, and quantum chemical methods are possible solutions. Herein we reported a scaling method for the prediction of 19F NMR chemical shift with freely available ORCA program package. Performance of 31 DFT functionals coupled with 11 basis sets were evaluated and influence of geometry optimization was also studied with five functionals coupled with three basis sets. The significance of geometry was further examined through the execution of relaxed surface scans of seven flexible compounds, and averaged shieldings of obtained conformers yielded notable improvement of the correlation between calculated isotropic shielidings and experimental chemical shifts. Utilization of the best scaling factor obtained successfully assigned of fluorine atoms in multifluorinated molecules with different conformations. The method reported here was computationally inexpensive, easily available with acceptable accuracy.

Novel Computational Chemistry Infrastructure for Simulating Astatide in Water: From Basis Sets to Force Fields Using Particle Swarm Optimization

Using the example of astatine, the heaviest naturally occurring halogen whose isotope At-211 has promising medical applications, we propose a new infrastructure for large-scale computational models of heavy elements with strong relativistic effects. In particular, we focus on developing an accurate force field for At- in water based on reliable relativistic density functional theory (DFT) calculations. To ensure the reliability of such calculations, we design novel basis sets for relativistic DFT, via the particle swarm optimization algorithm to optimize the coefficients of the new basis sets and the polarization-consistent basis set idea's extension to heavy elements to eliminate the basis set error from DFT calculations. The resulting basis sets enable the well-grounded evaluation of relativistic DFT against "gold-standard" CCSD(T) results. Accounting for strong relativistic effects, including spin-orbit interaction, via our redesigned infrastructure, we elucidate a noticeable dissimilarity between At- and I- in halide-water force field parameters, radial distribution functions, diffusion coefficients, and hydration energies. This work establishes the framework for the systematic development of polarization-consistent basis sets for relativistic DFT and accurate force fields for molecular dynamics simulations to be used in large-scale models of complex molecular systems with elements from the bottom of the periodic table, including actinides and even superheavy elements.

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.

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.

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.

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.

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

The applicability of reduced scaling algorithms based on auxiliary subspace methods for the correlation energy from the random phase approximation (RPA) as well as the correlation part of the self-energy obtained from the GW method is outlined for time-reversal symmetry breaking Kohn-Sham (KS) references. The updated algorithms allow for an efficient evaluation of RPA energies and GW quasiparticle energies for molecular systems with KS references that break time-reversal symmetry. The latter occur, for example, in magnetic fields. Furthermore, KS references for relativistic open-shell molecules also break time-reversal symmetry due to the single determinant ansatz used. Errors of the updated reduced-scaling algorithms are shown to be negligible compared to reference implementations, while the overall computational scaling is reduced by 2 orders of magnitude. Ionization energies obtained from the GW approximation are shown to be robust even for the electronically complicated group of trivalent lanthanoid ions. Starting from GW quasiparticle energies, it is subsequently shown that light-matter interactions of these systems can be calculated using the Bethe-Salpeter equation (BSE). Using the combined GW-BSE method, the absorption and emission spectra of a molecular europium(III) complex can be obtained including spin-orbit coupling.

Reducing Exact Two-Component Theory for NMR Couplings to a One-Component Approach: Efficiency and Accuracy

The self-consistent and complex spin-orbit exact two-component (X2C) formalism for NMR spin-spin coupling constants [ J. Chem. Theory Comput. 17, 2021, 3874-3994] is reduced to a scalar one-component ansatz. This way, the first-order response term can be partitioned into the Fermi-contact (FC) and spin-dipole (SD) interactions as well as the paramagnetic spin-orbit (PSO) contribution. The FC+SD terms are real and symmetric, while the PSO term is purely imaginary and antisymmetric. The relativistic one-component approach is combined with a modern density functional treatment up to local hybrid functionals including the response of the current density. Computational demands are reduced by factors of 8-24 as shown for a large tin compound consisting of 137 atoms. Limitations of the current ansatz are critically assessed for Sn, Pb, Pd, and Pt compounds, i.e. the one-component treatment is not sufficient for tin compounds featuring a few heavy halogen atoms.

Accurate Spin-Orbit Coupling by Relativistic Mixed-Reference Spin-Flip-TDDFT

Relativistic mixed-reference spin-flip (MRSF)-TDDFT is developed considering the spin-orbit coupling (SOC) within the mean-field approximation. The resulting SOC-MRSF faithfully reproduces the experiments with very high accuracy, which is also consistent with the values by four-component (4c) relativistic CASSCF and 4c-CASPT2 in the spin-orbit-energy splitting calculations of the C, Si, and Ge atoms. Even for the fifth-row element Sn, the SOC-MRSF yielded accurate splittings (∼ 3 % error). In the SOC calculations of the molecular 4-thiothymine with a third-row element, SOC-MRSF values are in excellent agreement with those of the SO-GMC-QDPT2 level, regardless of geometries and exchange-correlation functionals. The same SOC-MRSF predicted the anticipated chance of S₁ (nπ*) → T₁ (ππ*) intersystem crossing, even in thymine with only second-row elements. With its accuracy and practicality, thus, SOC-MRSF is a promising electronic structure protocol in challenging situations such as nonadiabatic molecular dynamics (NAMD) incorporating both internal conversions and intersystem crossings in large systems.

Investigating the Heaviest Halogen: Lessons Learned from Modeling the Electronic Structure of Astatine's Small Molecules

We present a systematic study of electron-correlation and relativistic effects in diatomic molecular species of the heaviest halogen astatine (At) within relativistic single- and multireference coupled-cluster approaches and relativistic density functional theory. We establish revised reference ab initio data for the ground states of At₂, HAt, AtAu, and AtO⁺ using a highly accurate relativistic effective core potential model and in-house basis sets developed for accurate modeling of molecules with large spin-orbit effects. Spin-dependent relativistic effects on chemical bonding in the ground state are comparable to the binding energy or even exceed it in At₂. Electron-correlation effects near the equilibrium internuclear separation are mostly dynamical and can be adequately captured using single-reference CCSD(T). However, bond elongation in At₂ and, especially, AtO⁺ results in rapid manifestation of its multireference character. While useful for evaluating the spin-orbit effects on the ground-state bonding and properties, the two-component density functional theory lacks predictive power, especially in combination with popular empirically adjusted exchange-correlation functionals. This drawback supports the necessity to develop new functionals for reliable quantum-chemical models of heavy-element compounds with strong relativistic effects.

The CH-3Σ+ Anion: Inelastic Rate Coefficients from Collisions with He at Interstellar Conditions

We present accurate ab initio calculations on several properties of a gas-phase system of interest in the interstellar medium (ISM), where the title molecular anion has been often surmised but not yet confirmed by observations. The CH^-3Σ⁺ constitutes the smallest term in the series of longer anionic polyynes which have been observed in the ISM (e.g., C₄H^- and several others). Hence, its dynamical behavior in collision with He atoms, one of the most abundant atoms in that environment, can provide quantitative indicators on the changes which can occur in the rotational state population of the title anion when driven by this collision dynamics. We therefore report an accurate evaluation of the full potential energy surface (PES) which acts between the molecular anion in its ground vibrational state and the He atom. The relevant inelastic scattering cross sections and the corresponding inelastic rate coefficients are then computed within a quantum treatment of the collisions. We find that the fairly small values of the final inelastic rate coefficients indicate state-changing processes by collisions to be inefficient paths for modifying the rotational state populations of this anion and therefore to aid its possible observation from direct radiative emission in the microwave region.

Electron Correlation or Basis Set Quality: How to Obtain Converged and Accurate NMR Shieldings for the Third-Row Elements?

The quality of theoretical NMR shieldings calculated at the quantum-chemical level depends on various theoretical aspects, of which the basis set type and size are among the most important factors. Nevertheless, not much information is available on the basis set effect on theoretical shieldings of the NMR-active nuclei of the third row. Here, we report on the importance of proper basis set selection to obtain accurate and reliable NMR shielding parameters for nuclei from the third row of the periodic table. All calculations were performed on a set of eleven compounds containing the elements Na, Mg, Al, Si, P, S, or Cl. NMR shielding tensors were calculated using the SCF-HF, DFT-B3LYP, and CCSD(T) methods, combined with the Dunning valence aug-cc-pVXZ, core-valence aug-cc-pCVXZ, Jensen polarized-convergent aug-pcSseg-n and Karlsruhe x2c-Def2 basis set families. We also estimated the complete basis set limit (CBS) values of the NMR parameters. Widely scattered nuclear shieldings were observed for the Dunning polarized-valence basis set, which provides irregular convergence. We show that the use of Dunning core-valence or Jensen basis sets effectively reduces the scatter of theoretical NMR results and leads to their exponential-like convergence to CBS. We also assessed the effect of vibrational, temperature, and relativistic corrections on the predicted shieldings. For systems with single bonds, all corrections are relatively small, amounting to less than 4% of the CCSD(T)/CBS value. Vibrational and temperature corrections were less reliable for H₃PO and HSiCH due to the high anharmonicity of the molecules. An abnormally high relativistic correction was observed for phosphorus in PN, reaching ~20% of the CCSD(T)/CBS value, while the correction was less than 7% for other tested molecules.

Paramagnetic NMR Shielding Tensors Based on Scalar Exact Two-Component and Spin-Orbit Perturbation Theory

The temperature-dependent Fermi-contact and pseudocontact terms are important contributions to the paramagnetic NMR shielding tensor. Herein, we augment the scalar-relativistic (local) exact two-component (X2C) framework with spin-orbit perturbation theory including the screened nuclear spin-orbit correction for the EPR hyperfine coupling and g tensor to compute these temperature-dependent terms. The accuracy of this perturbative ansatz is assessed with the self-consistent spin-orbit two-component and four-component treatments serving as reference. This shows that the Fermi-contact and pseudocontact interaction is sufficiently described for paramagnetic NMR shifts; however, larger deviations are found for the EPR spectra and the principle components of the EPR properties of heavy elements. The impact of the perturbative treatment is further compared to that of the density functional approximation and the basis set. Large-scale calculations are routinely possible with the multipole-accelerated resolution of the identity approximation and the seminumerical exchange approximation, as shown for [CeTi₆O₃(OiPr)₉(salicylate)₆].

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 funtionals is 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 constrains. 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, NMR spin-spincoupling 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 looses some ground for the NMR shifts.Therefore, the designed functional is a major step forwards for functionals that have been designed from first principles.

Communication: 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 NMR shieldings and 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 SCAN functional families. Thus, we strongly recommend to use the current density generalized τ in paramagnetic NMR and EPR calculations with meta-GGAs.

Quasi-Relativistic Calculation of EPR g Tensors with Derivatives of the Decoupling Transformation, Gauge-Including Atomic Orbitals, and Magnetic Balance

We present an exact two-component (X2C) ansatz for the EPR g tensor using gauge-including atomic orbitals (GIAOs) and a magnetically balanced basis set expansion. In contrast to previous X2C and four-component relativistic ansätze for the g tensor, this implementation results in a gauge-origin-invariant formalism. Furthermore, the derivatives of the relativistic decoupling matrix are incorporated to form the complete analytical derivative of the X2C Hamiltonian. To reduce the associated computational costs, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The quasi-relativistic X2C and DLU-X2C Hamiltonians accurately reproduce the results of the parent four-component relativistic theory when accounting for two-electron picture-change effects with the modified screened nuclear spin-orbit approximation in the respective one-electron integrals and integral derivatives. According to our benchmark studies, the uncontracted Dyall and segmented-contracted Karlsruhe x2c-type basis sets perform well when compared to large even-tempered basis sets. Moreover, (range-separated) hybrid density functional approximations such as LC-ωPBE and ωB97X-D are needed to match the experimental findings. The impact of the GIAOs depends on the distribution of the spin density, and their use may change the Δg shifts by 10-50% as shown for [(C₅Me₅)₂Y(μ-S)₂Mo(μ-S)₂Y(C₅Me₅)₂]^-. Routine calculations of large molecules are possible with widely available and comparably low-cost hardware as demonstrated for [Pt(C₆Cl₅)₄]^- with 3003 basis functions and three spin-(1/2) La(II) and Lu(II) compounds, for which we observe good agreement with the experimental findings.

NMR Coupling Constants Based on the Bethe-Salpeter Equation in the GW Approximation

We present the first steps to extend the Green's function GW method and the Bethe-Salpeter equation (BSE) to molecular response properties such as nuclear magnetic resonance (NMR) indirect spin-spin coupling constants. We discuss both a nonrelativistic one-component and a quasi-relativistic two-component formalism. The latter describes scalar-relativistic and spin-orbit effects and allows us to study heavy-element systems with reasonable accuracy. Efficiency is maintained by the application of the resolution of the identity approximation throughout. The performance is demonstrated using conventional central processing units (CPUs) and modern graphics processing units (GPUs) for molecules involving several thousand basis functions. Our results show that a large amount of Hartree-Fock exchange is vital to provide a sufficient Kohn-Sham starting point to compute the GW quasi-particle energies. As the GW-BSE approach is generally less accurate for triplet excitations or related properties such as the Fermi-contact interaction, the admixture of the Kohn-Sham correlation kernel through the contracted BSE (cBSE) method improves the results for NMR coupling constants. This leads to remarkable results when combined with the eigenvalue-only self-consistent variant (evGW) and Becke's half and half functional (BH&HLYP) or the CAM-QTP family. The developed methodology is used to calculate the Karplus curve of tin molecules, illustrating its applicability to extended chemically relevant molecules. Here, the GW-cBSE method improves upon the chosen BH&HLYP Kohn-Sham starting points.

Hyperfine Coupling Constants in Local Exact Two-Component Theory

We present a highly efficient implementation of the electron-nucleus hyperfine coupling matrix within the one-electron exact two-component (X2C) theory. The complete derivative of the X2C Hamiltonian is formed, that is, the derivatives of the unitary decoupling transformation are considered. This requires the solution of the response and Sylvester equations, consequently increasing the computational costs. Therefore, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The finite nucleus model is employed for both the scalar potential and the vector potential. Two-electron picture-change effects are modeled with the (modified) screened nuclear spin-orbit approach. Our implementation is fully integral direct and OpenMP-parallelized. An extensive benchmark study regarding the Hamiltonian, the basis set, and the density functional approximation is carried out for a set of 12-17 transition-metal compounds. The error introduced by DLU is negligible, and the DLU-X2C Hamiltonian accurately reproduces its four-component "fully" relativistic parent results. Functionals with a large amount of Hartree-Fock exchange such as CAM-QTP-02 and ωB97X-D are generally favorable. The pure density functional r²SCAN performs remarkably and even outperforms the common hybrid functionals TPSSh and CAM-B3LYP. Fully uncontracted basis sets or contracted quadruple-ζ bases are required for accurate results. The capability of our implementation is demonstrated for [Pt(C₆Cl₅)₄]^- with more than 4700 primitive basis functions and four rare-earth single-molecule magnets: [La(OAr*)₃]^-, [Lu(NR₂)₃]^-, [Lu(OAr*)₃]^-, and [TbPc₂]^-. Here, the results with the spin-orbit DLU-X2C Hamiltonian are in an excellent agreement with the experimental findings of all Pt, La, Lu, and Tb molecules.

Paramagnetic NMR Shielding Tensors and Ring Currents: Efficient Implementation and Application to Heavy Element Compounds

We present an efficient implementation of paramagnetic NMR shielding tensors and shifts in a nonrelativistic and scalar-relativistic density functional theory framework. For the latter, we make use of the scalar exact two-component Hamiltonian in its local approximation, and generally we apply the well established (multipole-accelerated) resolution of the identity approximation and the seminumerical exchange approximation. The perturbed density matrix of a paramagnetic NMR shielding calculation is further used to study the magnetically induced current density and ring currents of open-shell systems as illustrated for [U@Bi₁₂]^3-. [U@Bi₁₂]^3- features delocalized highest occupied molecular orbitals and sustains a net diatropic ring current of ca. 18 nA/T through the Bi₁₂ torus similar to the all-metal aromatic heavy-element cluster [Th@Bi₁₂]^4-.

Property-optimized Gaussian basis sets for lanthanides

Property-optimized Gaussian basis sets of split-valence, triple-zeta valence, and quadruple-zeta valence quality are developed for the lanthanides Ce-Lu for use with small-core relativistic effective core potentials. They are constructed in a systematic fashion by augmenting def2 orbital basis sets with diffuse basis functions and minimizing negative static isotropic polarizabilities of lanthanide atoms with respect to basis set exponents within the unrestricted Hartree-Fock method. The basis set quality is assessed using a test set of 70 molecules containing the lanthanides in their common oxidation states and f electron occupations. 5d orbital occupation turns out to be the determining factor for the basis set convergence of polarizabilities in lanthanide atoms and the molecular test set. Therefore, two series of property-optimized basis sets are defined. The augmented def2-SVPD, def2-TZVPPD, and def2-QZVPPD basis sets balance the accuracy of polarizabilities across lanthanide oxidation states. The relative errors in atomic and molecular polarizability calculations are ≤8% for augmented split-valence basis sets, ≤ 2.5% for augmented triple-zeta valence basis sets, and ≤1% for augmented quadruple-zeta valence basis sets. In addition, extended def2-TZVPPDD and def2-QZVPPDD are provided for accurate calculations of lanthanide atoms and neutral clusters. The property-optimized basis sets developed in this work are shown to accurately reproduce electronic absorption spectra of a series of LnCp₃ ^'- complexes (Cp' = C₅H₄SiMe₃, Ln = Ce-Nd, Sm) with time-dependent density functional theory.

NMR Indirect Spin-Spin Coupling Constants in a Modern Quasi-Relativistic Density Functional Framework

A quasi-relativistic implementation of NMR indirect spin-spin coupling constants is presented. The exact two-component (X2C) Hamiltonian and its diagonal local approximation to the unitary decoupling transformation (DLU) are utilized together with the (modified) screened nuclear spin-orbit approach. In a restricted kinetic balance, the finite nucleus model is available for both the scalar and vector potentials. The implementation supports density functionals up to the fourth rung of Jacob's ladder, i.e., (range-separated) hybrid and local hybrid functionals based on a seminumerical ansatz. We assess the quality of our quasi-relativistic X2C approach by comparison with "fully" relativistic four-component results for small main-group molecules and alkynyl compounds. The mean absolute error introduced by the DLU scheme is less than 0.05 × 10¹⁹ T J^-2 of the reduced coupling constant for the small main-group molecules and 0.5 Hz for the alkynyl compounds. Thus, the error is significantly smaller than finite nucleus size effects for heavy elements. The basis set convergence and the impact of different density functional approximations are further studied. We propose a simple scheme to develop segmented-contracted relativistic all-electron basis sets for NMR spin-spin couplings. Our implementation allows us to perform calculations of extended molecules with reasonable computational effort, which is illustrated for the ¹J(¹¹⁹Sn, ³¹P) coupling constant of a low-valent tin phosphinidenide complex. The corresponding results are in good agreement with the experimental findings.

Assessing the Accuracy of Local Hybrid Density Functional Approximations for Molecular Response Properties

A comprehensive overview of the performance of local hybrid functionals for molecular properties like excited states, ionization potentials within the GW framework, polarizabilities, magnetizabilities, NMR chemical shifts, and NMR spin-spin coupling constants is presented. We apply the generalization of the kinetic energy, τ, with the paramagnetic current density to all magnetic properties and the excitation energies from time-dependent density functional theory. This restores gauge invariance for these properties. Different ansätze for local mixing functions such as the iso-orbital indicator, the correlation length, the Görling-Levy second-order limit, and the spin polarization are compared. For the latter, we propose a modified version of the corresponding hyper-generalized gradient approximation functional of Perdew, Staroverov, Tao, and Scuseria (PSTS) [Phys. Rev. A 2008, 78, 052513] to allow for a numerically stable evaluation of the exchange-correlation kernel and hyperkernel. The PSTS functional leads to a very consistent improvement compared to the related TPSSh functional. It is further shown that the "best" choice of the local mixing function depends on the studied property and molecular class. While functionals based on the iso-orbital indicator lead to rather accurate excitation energies and ionization energies, the results are less impressive for NMR properties, for which a considerable dependence on the considered molecular test set and nuclei is observed. Johnson's local hybrid functional based on the correlation length yields remarkable results for NMR shifts of compounds featuring heavy elements and also for the excitation energies of organic compounds.