Relativistic pseudopotentials: their development and scope of applications

Structure and Excitation Spectra of Third-Row Transition Metal Hexafluorides Based on Multi-Reference Exact Two-Component Theory

The structures and some vertical excitation energies of third-row transition metal hexafluorides (MF6, M = Re, Os, Ir, Pt, Au, Hg) were calculated using the generalized-active-space configuration interaction (GASCI) theory based on the exact two-component (X2C) Hamiltonian. The spin-orbit coupling (SOC) was included at the Hartree-Fock level, enabling us to analyze the SOC at the orbital level (spinor-representation). The excitation spectra were assigned based on the double group, a relativistic group theory applicable to states with the SOC. This study provides a fundamental understanding of the ligand field splitting, including the SOC, that is useful for the photochemistry and spin chemistry involving heavy elements.

Lanthanide-Dependent Photochemical and Photophysical Properties of Lanthanide-Anthracene Complexes: Experimental and Theoretical Approaches

The structural, photophysical, and photochemical properties of Ln(depma)(hmpa)2(NO3)3 (Ln = La, Ce, Nd, Sm, Eu, Tb, Ho, Er, and Yb) complexes 1-Ln were investigated with a multidisciplinary approach involving synthesis, photocycloaddition-based crystal engineering, spectroscopic analytical techniques and quantum chemical ab initio calculations. Depending on the Ln3+ ion the isostructural 1-Ln complexes exhibit quite different behavior upon excitation at 350-400 nm. Some 1-Ln complexes (Ln = La, Ce, Sm, Tb, Yb) emit a broad and strong band near 533 nm arising from paired anthracene moieties, whereas others (Ln = Nd, Eu, Ho, Er) do not. 1-Eu is not emissive at all, whereas 1-Nd, 1-Ho, and 1-Er exhibit a Ln3+ based luminescence. Upon irradiation with 365 nm ultraviolet (UV) light 1-Ln (Ln = La, Ce, Sm, Tb, Yb) dimerize by means of a photochemically induced [4 + 4] cycloaddition of the anthracene moieties, whereas 1-Ln (Ln = Nd, Eu, Ho, Er) remain monomers. We propose three models, based on the matching of the energy levels between the Ln3+ ion and the paired or dimerized anthracene units in the energy-resonance crossing region, as well as on internal conversion-driven and intersystem crossing-driven energy transfer, which explain the Ln3+ ion regulated photophysics and photochemistry of the 1-Ln complexes.

Actinide endohedral inter-metalloid clusters of the group 15 elements

Inter-metalloid clusters in Zintl chemistry have been extensively studied due to their unique electronic structures and potential applications. In this work, we explored a series of actinide endohedral inter-metalloid clusters of the group 15 elements [An@Bi12]4- (An = Th-U) and [An@Sb12]4- using density functional theory (DFT). [Th@Bi12]4- and [U@Bi12]4- exhibit Cs symmetry, while [Pa@Bi12]4- and [An@Sb12]4- (An = Th-U) have C1 structures. Bonding analyses such as bond order, molecular orbitals (MO) and quantum theory of atoms in molecules (QTAIM) show covalent An-Bi/An-Sb bonding in the clusters. All these clusters are highly stable according to the studied formation reactions and may be accessible experimentally. Compared with [An@Bi12]4-, [An@Sb12]4- possesses stronger bonding interactions, mainly arising from the higher electrostatic interaction energy. For clusters with the same group 15 elements, the bonding interactions increase gradually from Th to U, which is mainly determined by the covalent interactions of An-Bi/An-Sb bonding. This work is expected to provide potential avenues for the construction of robust inter-metalloid clusters of the group 15 elements.

Unraveling the enigma of Craig-type Möbius-aromatic osmium compounds

Nuclear magnetic resonance (NMR) chemical shifts and the magnetically induced current density (MICD) susceptibility of four osmium containing molecules have been calculated at the density functional theory (DFT) level using three relativistic levels of theory. The calculations were performed at the quasi-relativistic level using an effective core potential (ECP) for Os, at the all-electron scalar exact two-component (X2C) relativistic level, and at the relativistic X2C level including spin-orbit coupling (SO-X2C). In earlier studies, the osmapentalene (1) and the osmapentalynes (2 and 3) were considered Craig-type Möbius aromatic and it was suggested that the analogous osmium compound (4) is Craig-type Möbius antiaromatic. Here, the ring-current strengths were obtained with the gauge including magnetically induced currents (GIMIC) method by integrating the MICD susceptibility passing through planes that intersect chemical bonds and by line integration of the induced magnetic field using Ampère-Maxwell's law. The ring-current calculations suggest that 1, 2 and 3 are weakly aromatic and that 4 is nonaromatic. The accuracy of the MICD susceptibility was assessed by comparing calculated NMR chemical shifts to available experimental data. Visualization of the MICD susceptibility shows that the ring current does not pass from one side of the molecular plane to the other, which means that the MICD susceptibility of the studied molecules does not exhibit any Möbius topology as one would expect for Craig-type Möbius aromatic and for Craig-type Möbius antiaromatic molecules. Thus, molecules 1-3 are not Craig-type Möbius aromatic and molecule 4 is not Craig-type Möbius antiaromatic as previously suggested. Calculations of the 1H NMR and 13C NMR chemical shifts of atoms near the Os atom show the importance of including spin-orbit effects. Overall, our study revisits the understanding of the aromaticity of organometallic molecules containing transition metals.

The aromatic nature of auracycles and diauracycles based on calculated ring-current strengths

We have calculated the magnetically induced current density susceptibility for gold-containing organometallic molecular rings using the gauge-including magnetically induced currents (GIMIC) method. The aromatic nature has been determined by calculating the strength of the magnetically induced ring current susceptibility, which is often called ring current. To our knowledge, we show here for the first time that gold-containing organometallic rings may be aromatic or antiaromatic sustaining ring currents in the presence of an external magnetic field. The calculated aromatic character of the rings agrees with the aromatic nature one expects when using Hückel's aromaticity rules. The studied auracycles and diauracycles with 4n electrons in the conjugated orbitals generally sustain a weak paratropic ring current, whereas those having 4n + 2 electrons in the conjugated orbitals sustain a diatropic ring current that is almost as strong as that of benzene. The number of electrons are obtained by assuming that each C, N and Au atom of the ring contribute one electron, and a H atom connected to a N atom in the ring increases the number of electrons by one. An electron-attracting ligand at Au removes one electron from the ring. Formation of a short Au-Au bonding diauracycles reduces the number of electrons in the ring by two.

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.

Distinguishing the Geometric and Electronic Structures of Actinide Carbides AnxC8 (An = Th, U; x = 2, 3) through Exchange Interactions

Global-minimum optimizations combined with relativistic quantum chemistry calculations have been performed to characterize the ground-state stable structures of four titled compounds and to analyze the bonding properties. Th2C8 was identified as being a ThC4-Th(C2)2 structure, U2C8 has been found to favor the U-U(C8) structure, and both Th3C8 and U3C8 adopt the (AnC3)2-(AnC2) structure. Then, the wave function analyses reveal that the interactions between the Th 7s-based orbital and the σg molecular orbital of the C2 unit compensate for the excitation energy of 7s16d1 → 6d2 and lead to the stabilization of two Th(IV)s in the ThC4-Th(C2)2 structure. It also reveals that the U species exhibit magnetic exchange coupling behavior in UxC8, for instance, as seen in the direct interaction of U2C8 and the superexchange pathway of U3C8, which effectively stabilizes their low-spin states. This interpretation indicates that the geometric and electronic structures of AnxC8 species are largely influenced by the local magnetic moment and spin correlation.

Level-Shifted Embedded Cluster Method for Modeling the Chemistry of Metal Oxides

The embedded cluster method has been used extensively in the study of the chemical and physical properties of metal oxides. This method has been a popular tool due to its relatively high accuracy and low computational cost. An even more promising option may entail integrating the embedded cluster method with the combined quantum mechanical and molecular mechanical (QM/MM) approach, thereby enabling further consideration of interactions within the entire system for superior results. We aim to accurately model the chemistry of metal oxides using this combined scheme. Here, using the prototypical MgO(100) surface as a test system, with Mg9O14 as the cluster in the quantum mechanical region, we show that the embedded cluster with untailored boundary effective core potentials (ECPs) can have frontier orbital energy levels that substantially deviate from the quantum mechanical reference results. This occurs even when Mg9O9, which retains the stoichiometry of MgO, is used as the cluster in the quantum mechanical region. As a result, the chemical properties of the embedded cluster models differ from those of the quantum mechanical reference model. To address this issue, we propose a new variant of the embedded cluster method called the level-shifted embedded cluster (LSEC) method, which allows the energy levels to be shifted to match the reference levels by tuning the boundary ECPs. Our validation calculations on the adsorption of various adsorbates with different properties on the MgO(100) surface show that the overall performance of QM/MM with the LSEC method is excellent for the adsorption energies, geometries, and charge properties. The excellent performance holds for both the nonstoichiometric and stoichiometric clusters (i.e., Mg9O14 and Mg9O9, respectively), demonstrating the robustness of the LSEC method. We expect that the LSEC method can be combined with QM/MM or used separately for future chemical studies of metal oxides and other ionically bonded systems.

Pd8 Cluster: Too Small to Melt? A BOMD Study

The question of whether a solid-liquid phase transition occurs in small clusters poses a fundamental challenge. In this study, we attempt to elucidate this phenomenon through a thorough examination of the thermal behavior and structural stability of Pd8 clusters employing ab initio simulations. Initially, a systematic global search is carried out to identify the various isomers of the Pd8 cluster. This is accomplished by employing an ab initio basin-hopping algorithm and using the PBE/SDD scheme integrated in the Gaussian code. The resulting isomers are further refined through reoptimization using the deMon2k package. To ensure the structural firmness of the lowest-energy isomer, we calculated normal modes. The structural stability as a function of temperature is analyzed through the Born-Oppenheimer molecular dynamics (BOMD) approach. Multiple BOMD trajectories at distinct simulated temperatures are examined with data clustering analysis to determine cluster isomers. This analysis establishes a connection between the potential energy landscape and the simulated temperature. To address the question of cluster melting, canonical parallel-tempering BOMD runs are performed and analyzed with the multiple-histogram method. A broad maximum in the heat capacity curve indicates a melting transition between 500 and 600 K. To further examine this transition, the mean-squared displacement and the pair-distance distribution function are calculated. The results of these calculations confirm the existence of a solid-liquid phase transition, as indicated by the heat capacity curve.

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.

Error of relativistic effective core potentials for closed-shell diatomic molecules of p-block heavy and superheavy elements in DFT and TDDFT calculations

Pseudopotentials (PP) are extensively used in electronic structure calculations, particularly for molecules containing heavy elements. Parameters in PPs are mainly determined from ab initio results, and errors of such PPs in density functional theory (DFT) calculations have been studied previously. However, PP errors on results with spin-orbit coupling and those in time-dependent DFT (TDDFT) calculations have not been reported previously. In this work, we investigate the error of the small-core energy-consistent Stuttgart/Koln pseudopotentials in DFT and TDDFT calculations with and without spin-orbit coupling. Ground state bond lengths, harmonic frequencies, dissociation energies, and vertical excitation energies for a series of closed-shell diatomic heavy and superheavy p-block molecules are calculated using several popular exchange-correlation functionals. PP errors are estimated by comparing with results using the all-electron Dirac-Coulomb (-Gaunt) Hamiltonian. Our results show that the difference between ground state properties and most excitation energies in scalar-relativistic calculations with the PP and those of all-electron calculations is quite small. This difference becomes somewhat larger when spin-orbit coupling (SOC) is present, especially for properties that are affected by SOC to some extent. In addition, the errors of the PPs are insensitive to the employed exchange-correlation functionals in most cases. Our results indicate that reasonable DFT and TDDFT results can be obtained using the small-core energy-consistent Stuttgart/Koln pseudopotentials for heavy and super-heavy p-block molecules.

Quantum computations in heavy noble-gas hydride cations: Reference energies and new spectroscopic data

Computational quantum chemistry has become a powerful tool with a wide range of possibilities to solve chemical-physical problems. As a result of this, the interest in the applications of computational quantum chemistry has expanded considerably, and has opened up novel research opportunities. In particular, those related to the characterization of heavy-atoms complexes, as most electronic structure calculations for such systems struggle with the problem posed by the large number of electrons present in them, and consequently, the introduction of relativistic effects. The present study performed an exhaustive assess to characterized the uncommon NgH+ (Ng = Kr, Xe, and Rn) hydride cations in order to provide accurate rovibrational data of their isotopes to assist in the laboratory characterization or even their astronomical detection. Scalar relativistic effects were included, and the ground and first electronically exited states potential curves were obtained from benchmark ab initio CCSD(T)/CBS and MRCI+Q electronic structure calculations. Next, such interaction potentials, correctly extended to long-range asymptotic regions, were employed in quantum bound state calculations and molecular spectroscopic constants were determined for the most abundance 84Kr, 132Xe, and 222Rn isotopes. Our results were discussed in comparison with available experimental and previous theoretical estimates, aiming to treat accuracy issues. The new sets provide reference data that could serve for spectroscopic characterization of such low abundance and high radioactive species.

CRYSTAL23: A Program for Computational Solid State Physics and Chemistry

The Crystal program for quantum-mechanical simulations of materials has been bridging the realm of molecular quantum chemistry to the realm of solid state physics for many years, since its first public version released back in 1988. This peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals (LCAO) approach and from the corresponding efficiency in the evaluation of the exact Fock exchange series. In particular, this has led to the implementation of a rich variety of hybrid density functional approximations since 1998. Nowadays, it is acknowledged by a broad community of solid state chemists and physicists that the inclusion of a fraction of Fock exchange in the exchange-correlation potential of the density functional theory is key to a better description of many properties of materials (electronic, magnetic, mechanical, spintronic, lattice-dynamical, etc.). Here, the main developments made to the program in the last five years (i.e., since the previous release, Crystal17) are presented and some of their most noteworthy applications reviewed.

Gd(III) and Yb(III) Complexes Derived from a New Water-Soluble Dioxopolyazacyclohexane Macrocycle

A new macrocyclic ligand was synthesized by a reaction between diethylenetriaminepentaacetic (DTPA) dianhydride and trans-1,4-diaminocyclohexane, and the Gd(III) and Yb(III) complexes were prepared. The compounds were characterized by spectroscopic methods. Structural calculation by DFT shows that the amide linkages are arranged in such a way that a conformational strain is minimized in the macrocyclic frame. The coordination modes of the ligand and water in the metal complexes were also determined by DFT. The longitudinal relaxation time T1 was measured for aqueous solutions of the Gd(III) complex. The T1 relaxivity arises from the structural feature that a water molecule coordinated to the paramagnetic metal is surrounded by a large open space, through which the exchange of water occurs readily to shorten the relaxation time of water in the entire region, as a result of the chelate conformation defined strictly by the amide groups and the cyclohexane ring.

Excited Electronic States of Sr2: Ab Initio Predictions and Experimental Observation of the 21Σu+ State

Despite its apparently simple nature with four valence electrons, the strontium dimer constitutes a challenge for modern electronic structure theory. Here we focus on excited electronic states of Sr₂, which we investigate theoretically up to 25000 cm^-1 above the ground state, to guide and explain new spectroscopic measurements. In particular, we focus on potential energy curves for the 1¹Σ_u⁺, 2¹Σ_u⁺, 1¹Π_u, 2¹Π_u, and 1¹Δ_u states computed using several variants of ab initio coupled-cluster and configuration-interaction methods to benchmark them. In addition, a new experimental study of the excited 2¹Σ_u⁺ state using polarization labeling spectroscopy is presented, which extends knowledge of this state to high vibrational levels, where perturbation by higher electronic states is observed. The available experimental observations are compared with the theoretical predictions and help to assess the accuracy and limitations of employed theoretical models. The present results pave the way for future more accurate theoretical and experimental spectroscopic studies.

The Good, the Bad, and the Ugly: Pseudopotential Inconsistency Errors in Molecular Applications of Density Functional Theory

The pseudopotential (PP) approximation is one of the most common techniques in computational chemistry. Despite its long history, the development of custom PPs has not tracked with the explosion of different density functional approximations (DFAs). As a result, the use of PPs with exchange/correlation models for which they were not developed is widespread, although this practice is known to be theoretically unsound. The extent of PP inconsistency errors (PPIEs) associated with this practice has not been systematically explored across the types of energy differences commonly evaluated in chemical applications. We evaluate PPIEs for a number of PPs and DFAs across 196 chemically relevant systems of both transition-metal and main-group elements, as represented by the W4-11, TMC34, and S22 data sets. Near the complete basis set limit, these PPs are found to cleanly approach all-electron (AE) results for noncovalent interactions but introduce root-mean-squared errors (RMSEs) upwards of 15 kcal mol^-1 into predictions of covalent bond energies for a number of popular DFAs. We achieve significant improvements through the use of empirical atom- and DFA-specific PP corrections, indicating considerable systematicity of the PPIEs. The results of this work have implications for chemical modeling in both molecular contexts and for DFA design, which we discuss.

Synergetic argentophilic and through space electronic interactions in a single-crystal-to-single-crystal photocycloaddition reaction: a mechanistic study

Ag(I) is able to mediate single-crystal-to-single-crystal transformation through [2+2] photocycloaddition to prepare high-conductivity materials. However, the intrinsic mechanism of Ag(I) mediation, the detailed photophysical and photochemical processes as well as the origin of the enhanced conductivity of nanocrystals are still unclear. In this work, the comprehensive kinetic scheme and regulation mechanism are established by the accurate QM/MM calculations at the CASPT2//CASSCF/AMBER level of theory with consideration of the crystal environment. We find that the argentophilic interaction and through space electronic interaction are the key factors that promote Ag(I)-mediated [2+2] PCA reactions and may account for the enhancement of conductivity. These mechanistic insights into the Ag(I)-regulated photo-dimerization in the crystal surrounding are beneficial for the design of the structurally and electrically favorable skeletons of a metal-organic coordination polymer.

Trends in the Electronic Structure and Chemical Bonding of a Series of Porphyrinoid-Uranyl Complexes

In this paper, we have explored the relativistic density functional theory study on a series of deprotonated porphyrinoid (Lⁿ) complexes of uranyl to investigate the geometrical structures and chemical bonding. The ligands bound with uranyl in the 1:1 complexes [UO₂(Lⁿ)]^x (n = 4, 5, 6; x = 0, -1, -2), showing more thermodynamic stability for "in-cavity" structures of L⁵ and L⁶ than that of the "side-on" structure of L⁴ and an increase in stability with the increase of negative charges, L^2- < L^3- < L^4-. Among the six ligands, the cyclo[6]pyrrole presents the best selectivity toward uranyl. Based on chemical bonding analyses, the U-N_L bond in the in-cavity complexes adopts a typical dative N_L → U bond with mainly ionic bonding and significant covalency, which comes from the significant orbital interaction of U 5f_ϕ6d_δ7s hybrid AOs and N_L 2p-based MOs. This work provides a systematic understanding of the coordination chemistry in uranyl pyrrole-containing macrocycle complexes and the nature of chemical bonding in such systems, which may provide inspirations for the future design of synthetic targets that could be relevant to actinide separations or in the remediation of spent nuclear fuel.

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.

Perturbation Theory Treatment of Spin-Orbit Coupling. III: Coupled Perturbed Method for Solids

A previously proposed noncanonical coupled-perturbed Kohn-Sham density functional theory (KS-DFT)/Hartree-Fock (HF) treatment for spin-orbit coupling is here generalized to infinite periodic systems. The scalar-relativistic periodic KS-DFT/HF solution, obtained with a relativistic effective core potential, is taken as the zeroth-order approximation. Explicit expressions are given for the total energy through third-order, which satisfy the 2N + 1 rule (i.e., requiring only the first-order perturbed wave function for determining the energy through third-order). Expressions for additional second-order corrections to the perturbed wave function (as well as related one-electron properties) are worked out at the uncoupled-perturbed level of theory. The approach is implemented in the Crystal program and validated with calculations of the total energy, electronic band structure, and density variables of spin-current DFT on the tungsten dichalcogenide hexagonal bilayer series (i.e., WSe₂, WTe₂, WPo₂, WLv₂), including 6p and 7p elements as a stress test. The computed properties through second- or third-order match well with those from reference two-component self-consistent field (2c-SCF) calculations. For total energies, E⁽³⁾ was found to consistently improve the agreement against the 2c-SCF reference values. For electronic band structures, visible differences w.r.t. 2c-SCF remained through second-order in only the single-most difficult case of WLv₂. As for density variables of spin-current DFT, the perturbed electron density, being vanishing in first-order, is the most challenging for the perturbation theory approach. The visible differences in the electron densities are, however, largest close to the core region of atoms and smaller in the valence region. Perturbed spin-current densities, on the other hand, are well reproduced in all tested cases.

The Importance of Tight f Basis Functions for Heavy p-Block Oxides and Halides: A Parallel With Tight d functions in the Second Row

It is well-known that both wave function ab initio and DFT calculations on second-row compounds exhibit anomalously slow basis set convergence unless the basis sets are augmented with additional "tight" (high-exponent) d functions, as in the cc-pV(n+d)Z and aug-cc-pV(n+d)Z basis sets. This has been rationalized as being necessary for a better description of the low-lying 3d orbital, which as the oxidation state increases sinks low enough to act as a back-donation acceptor from chalcogen and halogen lone pairs. This prompts the question whether a similar phenomenon exists for the isovalent compounds of the heavy p-block. We show that for the fourth and fifth row, this is the case, but this time for tight f functions enhancing the description of the low-lying 4f and 5f Rydberg orbitals, respectively. In the third-row heavy p block, the 4f orbitals are too far up, while the 4d orbitals are adequately covered by the basis functions already present to describe the 3d subvalence orbitals.

Mechanistic insights into aerobic oxidative cleavage of glycol catalyzed by an Anderson-type polyoxometalate [IMo6O24]5

A computational investigation of the aerobic oxidative C-C bond cleavage reaction of glycol catalyzed by an Anderson-type heteropolyanion HPA [IMo₆O₂₄]^5- in the presence of acetonitrile as solvent has been performed at the WB97XD/6-31G(d,p)/lanl2dz level. Two reaction pathways have been identified. The catalytic cycle of each pathway consists of three steps: oxidation cleavage of a glycol molecule by the HPA, oxidation of the HPA by one dioxygen molecule, and, finally, oxidation of a second glycol and regeneration of the catalyst. These reaction pathways have been thoroughly investigated in terms of energetic, natural bond orbital (NBO), natural charges, and geometrical parameters. It is found that (i) even though the top oxygen atoms of the Anderson heteropolyanion are not the most negatively charged ones, they are more likely to react with the diol hydroxyl groups, (ii) a direct relationship between the presence of the iodine ion I(VII) and the studied oxidation reaction could not be identified, and (iii) in terms of energy, the transfer of the two hydrogen atoms is the most energetic step.

Norm-Conserving 4f-in-Core Pseudopotentials and Basis Sets Optimized for Trivalent Lanthanides (Ln = Ce-Lu)

We present here a set of scalar-relativistic norm-conserving 4f-in-core pseudopotentials, together with complementary valence-shell Gaussian basis sets, for the lanthanide (Ln) series (Ce-Lu). The Goedecker, Teter, and Hutter (GTH) formalism is adopted with the generalized gradient approximation (GGA) for the exchange-correlation Perdew-Burke-Ernzerhof (PBE) functional. The 4f-in-core pseudopotentials are built through attributing 4f-subconfiguration 4fⁿ (n = 1-14) for Ln (Ln = Ce-Lu) into the atomic core region, making it possible to circumvent the difficulty of the description of the open 4fⁿ valence shell. A wide variety of computational benchmarks and tests have been carried out on lanthanide systems including Ln³⁺-containing molecular complexes, aqueous solutions, and bulk solids to validate the accuracy, reliability, and efficiency of the optimized 4f-in-core GTH pseudopotentials and basis sets. The 4f-in-core GTH pseudopotentials successfully replicate the main features of lanthanide structural chemistry and reaction energetics, particularly for nonredox reactions. The chemical bonding features and solvation shells, hydrolysis energetics, acidity constants, and solid-state properties of selected lanthanide systems are also discussed in detail by utilizing these new 4f-in-core GTH pseudopotentials. This work bridges the idea of keeping highly localized 4f electrons in the atomic core and efficient pseudopotential formalism of GTH, thus providing a highly efficient approach for studying lanthanide chemistry in multi-scale modeling of constituent-wise and structurally complicated systems, including electronic structures of the condensed phase and first-principles molecular dynamics simulations.

Computational insight into a mechanistic overview of water exchange kinetics and thermodynamic stabilities of bis and tris-aquated complexes of lanthanides

A thorough investigation of Ln³⁺ complexes with more than one inner-sphere water molecule is crucial for designing high relaxivity contrast agents (CAs) used in magnetic resonance imaging (MRI). This study accomplished a comparative stability analysis of two hexadentate (H₃cbda and H₃dpaa) and two heptadentate (H₄peada and H₃tpaa) ligands with Ln³⁺ ions. The higher stability of the hexadentate H₃cbda and heptadentate H₄peada ligands has been confirmed by the binding affinity and Gibbs free energy analysis in aqueous solution. In addition, energy decomposition analysis (EDA) reveals the higher binding affinity of the peada^4- ligand than the cbda^3- ligand towards Ln³⁺ ions due to the higher charge density of the peada^4- ligand. Moreover, a mechanistic overview of water exchange kinetics has been carried out based on the strength of the metal-water bond. The strength of the metal-water bond follows the trend Gd-O47 (w) > Gd-O39 (w) > Gd-O36 (w) in the case of the tris-aquated [Gd(cbda)(H₂O)₃] and Gd-O43 (w) > Gd-O40 (w) for the bis-aquated [Gd(peada)(H₂O)₂]^- complex, which was confirmed by bond length, electron density (ρ), and electron localization function (ELF) at the corresponding bond critical points. Our analysis also predicts that the activation energy barrier decreases with the decrease in bond strength; hence k _ex increases. The ¹⁷O and ¹H hyperfine coupling constant values of all the coordinated water molecules were different, calculated by using the second-order Douglas-Kroll-Hess (DKH2) approach. Furthermore, the ionic nature of the bonding in the metal-ligand (M-L) bond was confirmed by the Quantum Theory of Atoms-In-Molecules (QTAIM) and ELF along with energy decomposition analysis (EDA). We hope that the results can be used as a basis for the design of highly efficient Gd(iii)-based high relaxivity MRI contrast agents for medical applications.

Theoretical Exploration of Energy Transfer and Single Electron Transfer Mechanisms to Understand the Generation of Triplet Nitrene and the C(sp3)-H Amidation with Photocatalysts

Mechanistic explorations and kinetic evaluations were performed based on electronic structure calculations at the CASPT2//CASSCF level of theory, the Fermi's golden rule combined with the Dexter model, and the Marcus theory to unveil the key factors regulating the processes of photocatalytic C(sp³)-H amidation starting from the newly emerged nitrene precursor of hydroxamates. The highly reactive nitrene was found to be generated efficiently via a triplet-triplet energy transfer process and to be benefited from the advantages of hydroxamates with long-range charge-transfer (CT) excitation from the N-centered lone pair to the 3,5-bis(trifluoromethyl)benzoyl group. The properties of the metal-to-ligand charge-transfer (MLCT) state of photocatalysts, the functionalization of chemical moieties for substrates involved in the charge-transfer (CT) excitation, such as the electron-withdrawing trifluoromethyl group, and the energetic levels of singlet and triplet reaction pathways may regulate the reaction yield of C(sp³)-H amidation. Kinetic evaluations show that the triplet-triplet energy transfer is the main driving force of the reaction rather than the single electron transfer process. The effects of electronic coupling, molecular rigidity, and excitation energies on the energy transfer efficiency were further discussed. Finally, we investigated the inverted behavior of single-electron transfer, which is correlated unfavorably to the catalytic efficiency and amidation reaction. All theoretical explorations allow us to better understand the generation of nitrene with visible-light photocatalysts, to expand highly efficient substrate sources, and to broaden our scope of available photosensitizers for various cross-coupling reactions and the construction of N-heterocycles.

Correlation Consistent Basis Sets and Core Polarization Potentials for Al-Ar with ccECP Pseudopotentials

New correlation consistent basis sets for the second-row atoms (Al–Ar) to be used with the neon-core correlation consistent effective core potentials (ccECPs) have been developed. The basis sets, denoted cc-pV(n+d)Z-ccECP (n = D, T, Q), include the “tight”-d functions that are known to be important for second-row elements. Sets augmented with additional diffuse functions are also reported. Effective core polarization potentials (CPPs) to account for the effect of core–valence correlation have been adjusted for the same elements, and two different forms of the CPP cutoff function have been analyzed. The accuracy of both the basis sets and the CPPs is assessed through benchmark calculations at the coupled-cluster level of theory for atomic and molecular properties. Agreement with all-electron results is much improved relative to the basis sets that originally accompanied the ccECPs; moreover, the combination of cc-pV(n+d)Z-ccECP and CPPs is found to be a computationally efficient and accurate alternative to including core electrons in the correlation treatment.

A new generation of effective core potentials from correlated and spin-orbit calculations: selected heavy elements

We introduce new correlation consistent effective core potentials (ccECPs) for the elements I, Te, Bi, Ag, Au, Pd, Ir, Mo, and W with $4d$, $5d$, $6s$ and $6p$ valence spaces. These ccECPs are given as a sum of spin-orbit averaged relativistic effective potential (AREP) and effective spin-orbit (SO) terms. The construction involves several steps with increasing refinements from more simple to fully correlated methods. The optimizations are carried out with objective functions that include weighted many-body atomic spectra, norm-conservation criteria, and spin-orbit splittings. Transferability tests involve molecular binding curves of corresponding hydride and oxide dimers. The constructed ccECPs are systematically better and in a few cases on par with previous effective core potential (ECP) tables on all tested criteria and provide a significant increase in accuracy for valence-only calculations with these elements. Our study confirms the importance of the AREP part in determining the overall quality of the ECP even in the presence of sizable spin-orbit effects. The subsequent quantum Monte Carlo (QMC) calculations point out the importance of accurate trial wave functions which in some cases (mid series transition elements) require treatment well beyond single-reference.

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.

Excited-State Dynamics of Crossing-Controlled Energy Transfer in Europium Complexes

Photosensitized energy transfer (EnT) phenomena occur frequently in a variety of photophysical and photochemical processes and have traditionally been treated with the donor-acceptor distance-dependent Förster and Dexter models. However, incorrect arguments and formulae were employed by ignoring energy resonance conditions and the selection rules of the state-to-state transition in special cases, especially for the sensitive intramolecular EnT of lanthanide complexes. Herein, we proposed an innovative model of energy-degeneracy-crossing-controlled EnT, which can be experimentally confirmed by time-resolved two-dimensional photoluminescence measurements. The computationally determined energy resonance region provides the most effective channel to achieve metal-to-ligand EnT beyond the distance-dependent model and sensitively bifurcates into symmetry-allowed or -forbidden channels for some representative europium antenna complexes. The outcomes of the multidisciplinary treatment contribute to a complementary EnT model that can be tuned by introducing a phosphorescence modulator and altering the antenna-related parameters of the ligand-centered energy level of the ³ππ* state and its spin-orbit coupling for the ³ππ* → S₀ ^* transition through mechanism-guided crystal engineering and should motivate further development of mechanistic models and applications.

Small-Basis Set Density-Functional Theory Methods Corrected with Atom-Centered Potentials

Density functional theory (DFT) is currently the most popular method for modeling noncovalent interactions and thermochemistry. The accurate calculation of noncovalent interaction energies, reaction energies, and barrier heights requires choosing an appropriate functional and, typically, a relatively large basis set. Deficiencies of the density-functional approximation and the use of a limited basis set are the leading sources of error in the calculation of noncovalent and thermochemical properties in molecular systems. In this article, we present three new DFT methods based on the BLYP, M06-2X, and CAM-B3LYP functionals in combination with the 6-31G* basis set and corrected with atom-centered potentials (ACPs). ACPs are one-electron potentials that have the same form as effective-core potentials, except they do not replace any electrons. The ACPs developed in this work are used to generate energy corrections to the underlying DFT/basis-set method such that the errors in predicted chemical properties are minimized while maintaining the low computational cost of the parent methods. ACPs were developed for the elements H, B, C, N, O, F, Si, P, S, and Cl. The ACP parameters were determined using an extensive training set of 118655 data points, mostly of complete basis set coupled-cluster level quality. The target molecular properties for the ACP-corrected methods include noncovalent interaction energies, molecular conformational energies, reaction energies, barrier heights, and bond separation energies. The ACPs were tested first on the training set and then on a validation set of 42567 additional data points. We show that the ACP-corrected methods can predict the target molecular properties with accuracy close to complete basis set wavefunction theory methods, but at a computational cost of double-ζ DFT methods. This makes the new BLYP/6-31G*-ACP, M06-2X/6-31G*-ACP, and CAM-B3LYP/6-31G*-ACP methods uniquely suited to the calculation of noncovalent, thermochemical, and kinetic properties in large molecular systems.

Intersystem Crossings in Late-Row Elements: A Perspective

Intersystem crossing (ISC), a vital component of the electronic and nuclear transitions that compose photophysics, has been successfully simulated in light elements and transition metal complexes. Derived from the Z-dependent spin-orbit coupling (SOC), ISC is expected to be of greater importance in heavier elements, but few attempts have been made at the simulation of ISC in lanthanides or actinides. In this work, we explore several of the challenges that will need to be overcome in order to treat ISC in late-row elements, including the loss of spin as a good quantum number, the need to include SOC variationally via two- or four-component electronic structure, and the high density of states present in late-row complexes. Density functional theory (DFT) calculations are used to illustrate several of these effects, while a model Hamiltonian is used to illustrate the importance of momentum rescaling in surface hopping simulations of strongly coupled states.

Intermediate Hamiltonian Fock-space coupled-cluster theory for excitation energies, double ionization potentials, and double electron attachments with spin–orbit coupling

The intermediate Hamiltonian Fock-space coupled-cluster methods at the singles and doubles level (IHFSCCSD) for excitation energies in the (1p, 1h) sector, double ionization potentials in the (0p, 2h) sector, and double electron attachments in the (2p, 0h) sector of the Fock space are implemented based on the CCSD method with spin–orbit coupling (SOC) included in the post-Hartree–Fock treatment using a closed-shell reference in this work. The active space is chosen to contain those orbitals that have the largest contribution to principal ionized or electron-attached states obtained from the equation-of-motion coupled-cluster calculations. Both time-reversal symmetry and spatial symmetry are exploited in the implementation. Our results show that the accuracy of IHFSCCSD results is closely related to the active space, and the sufficiency of the active space can be assessed from the percentage of transitions within the active space. In addition, unreasonable results may be encountered when the ionized or electron-attached states with a somewhat larger contribution from double excitations are included to determine the active space and cluster operators in the (0p, 1h) or (1p, 0h) sector of the Fock space. A larger active space may be required to describe SO splitting reliably than that in the scalar-relativistic calculations in some cases. The IHFSCCSD method with SOC developed in this work can provide reliable results for heavy-element systems when a sufficient active space built upon the principal ionization potential/electron affinity states is adopted.

Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials

There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.

Employing Pseudopotentials to Tackle Excited-State Electron Spill-Out in Frozen Density Embedding Calculations

In frozen density embedding (FDE), the properties of a target molecule are computed in the presence of an effective embedding potential, which accounts for the attractive and repulsive contributions of the environment. The formally exact embedding potential, however, is in practice calculated using explicit kinetic-energy functionals for which the resulting potentials are in many cases not repulsive enough to account fully for Pauli repulsion by the electrons of the environment and to compensate thereby the strong electron-nuclear attraction. For the excited states on the target molecule, this leads to charge spill-out when diffuse basis functions are included, which allow that valence electrons are excited to those regions of the environment where the strong nuclear attraction is not sufficiently compensated by repulsive contributions. To reduce this insufficiency, we propose in the present work the inclusion of atomic all-electron pseudopotentials for all environment atoms on top of the conventional embedding potential. In the current work, the pseudopotentials are applied for computing vertical excitation energies of local excited states in complex systems employing the second-order algebraic diagrammatic construction (ADC(2)) scheme. The proposed approach leads to significantly reduced charge spill-out and an improved agreement of FDE and supermolecular calculations in the frozen solvent approximation. In particular, when diffuse functions are employed, the mean absolute deviation (MAD) is reduced from 0.27 to 0.05 eV for the investigated cases.

The MOBH35 Metal–Organic Barrier Heights Reconsidered: Performance of Local-Orbital Coupled Cluster Approaches in Different Static Correlation Regimes

We have revisited the MOBH35 (Metal–Organic Barrier Heights, 35 reactions) benchmark [Iron, Janes, J. Phys. Chem. A, 2019, 123 ( (17), ), 3761−378130973722; ibid. 2019, 123, 6379–6380] for realistic organometallic catalytic reactions, using both canonical CCSD(T) and localized orbital approximations to it. For low levels of static correlation, all of DLPNO-CCSD(T), PNO-LCCSD(T), and LNO-CCSD(T) perform well; for moderately strong levels of static correlation, DLPNO-CCSD(T) and (T₁) may break down catastrophically, and PNO-LCCSD(T) is vulnerable as well. In contrast, LNO-CCSD(T) converges smoothly to the canonical CCSD(T) answer with increasingly tight convergence settings. The only two reactions for which our revised MOBH35 reference values differ substantially from the original ones are reaction 9 and to a lesser extent 8, both involving iron. For the purpose of evaluating density functional theory (DFT) methods for MOBH35, it would be best to remove reaction 9 entirely as its severe level of static correlation makes it just too demanding for a test. The magnitude of the difference between DLPNO-CCSD(T) and DLPNO-CCSD(T₁) is a reasonably good predictor for errors in DLPNO-CCSD(T₁) compared to canonical CCSD(T); otherwise, monitoring all of T₁, D₁, max|t_i^A|, and 1/(ε_LUMO – ε_HOMO) should provide adequate warning for potential problems. Our conclusions are not specific to the def2-SVP basis set but are largely conserved for the larger def2-TZVPP, as they are for the smaller def2-SV(P): the latter may be an economical choice for calibrating against canonical CCSD(T). Finally, diagnostics for static correlation are statistically clustered into groups corresponding to (1) importance of single excitations in the wavefunction; (2a) the small band gap, weakly separated from (2b) correlation entropy; and (3) thermochemical importance of correlation energy, as well as the slope of the DFT reaction energy with respect to the percentage of HF exchange. Finally, a variable reduction analysis reveals that much information on the multireference character is provided by T₁, I_ND/I_tot, and the exchange-based diagnostic A₁₀₀[TPSS].

Mechanistic Insights into Visible-Light-Driven Dearomative Fluoroalkylation Mediated by an Electron Donor-Acceptor Complex

Electron donor-acceptor (EDA) complex photochemistry has become a burgeoning topic in the synthetic radical chemistry mediated by visible light; however, the theoretical insights into the reaction mechanisms are limited. Herein, accurate electronic structure calculations at the CASPT2//CASSCF/PCM level of theory were performed to investigate the paradigm example of EDA complex-enabled photoreaction for visible-light-driven dearomative perfluoroalkylation of β-naphthol. The excitation energy levels of the EDA complex are controlled by noncovalent interactions because the photoinduced intermolecular charge is enhanced when the noncovalent interaction becomes weaker, leading to the broad spectra ranging from UVA (<380 nm) to visible light (>500 nm). The competitiveness of the radical-radical coupling over the radical chain pathway is also regulated due to the tunable radical concentrations varying the excitation wavelength.

Ab initio calculations in atoms, molecules, and solids, treating spin-orbit coupling and electron interaction on an equal footing

We incorporate explicit, non-perturbative treatment of spin-orbit coupling into ab initio auxiliary-field quantum Monte Carlo (AFQMC) calculations. The approach allows a general computational framework for molecular and bulk systems in which material specificity, electron correlation, and spin-orbit coupling effects can be captured accurately and on an equal footing, with favorable computational scaling vs system size. We adopt relativistic effective-core potentials that have been obtained by fitting to fully relativistic data and that have demonstrated a high degree of reliability and transferability in molecular systems. This results in a two-component spin-coupled Hamiltonian, which is then treated by generalizing the ab initio AFQMC approach. We demonstrate the method by computing the electron affinity in Pb, the bond dissociation energy in Br₂ and I₂, and solid Bi.

Development and Application of a Complete Active Space Spin-Orbit Configuration Interaction Program Designed for Single Molecule Magnets

We present a spin-orbit configuration interaction program which has been tailored for the description of the magnetic properties of polynuclear metal complexes with partially filled d- and f-shells. The spin-orbit operators are directly included in the configuration interaction program based on Slater-determinants. The lowest states are obtained by a Block-Davidson-type diagonalisation. The usage of localised active orbitals enables the construction of start vectors from tensor products of single-center wave functions that already include spin-orbit interaction. This allows for an analysis of the role and the interplay of the different metal centres. Furthermore, in case of weak coupling of the metal centres these tensor products are already close to the final wave functions ensuring fast convergence. In combination with a two-layer hybrid parallelisation, this makes the program highly efficient. Based on the spin-orbit coupled wave functions, magnetic D-tensors, g-tensors and temperature-dependent susceptibilities can be calculated. The applicability and performance of the program is shown exemplarily on a trinuclear transition metal (CoII VII CoII ) complex.

Highly stable actinide(III) complexes supported by doubly aromatic ligands

Owing to the electron-deficient nature of boron atom, the structures and properties of boron clusters can be enriched by doping various metal atoms, including lanthanide metal atoms. Nevertheless, the viability...

Interaction-determined extraction capacity between rare earth ions and extractants: taking lanthanum and lutetium as models through theoretical calculations

Extractant plays an important role in the separation and purification of rare earth elements (REEs), whereas, extraction performance is the most effective tactic to evaluate whether an extractant is complete...

Electronic structure and bonding in endohedral Zintl clusters

Endohedral Zintl clusters-multi-metallic anionic molecules in which a d-block or f-block metal atom is enclosed by p-block (semi)metal atoms-are very topical in contemporary inorganic chemistry. Not only do they provide insight into the embryonic states of intermetallic compounds and show promise in catalytic applications, they also shed light on the nature of chemical bonding between metal atoms. Over the past two decades, a plethora of endohedral Zintl clusters have been synthesized, revealing a fascinating diversity of molecular architectures. Many different perspectives on the bonding in them have emerged in the literature, sometimes complementary and sometimes conflicting, and there has been no concerted effort to classify the entire family based on a small number of unifying principles. A closer look, however, reveals distinct patterns in structure and bonding that reflect the extent to which valence electrons are shared between the endohedral atom and the cluster shell. We show that there is a much more uniform relationship between the total valence electron count and the structure and bonding patterns of these clusters than previously anticipated. All of the p-block (semi)metal shells can be placed on a ladder of total valence electron count that ranges between 4n+2 (closo deltahedra), 5n (closed, three-bonded polyhedra) and 6n (crown-like structures). Although some structural isomerism can occur for a given electron count, the presence of a central metal cation imposes a preference for rather regular and approximately spherical structures which maximise electrostatic interactions between the metal and the shell. In cases where the endohedral metal has relatively accessible valence electrons (from the d or f shells), it can also contribute its valence electrons to the total electron count of the cluster shell, raising the effective electron count and often altering the structural preferences. The electronic situation in any given cluster is considered from different perspectives, some more physical and some more chemical, in a way that highlights the important point that, in the end, they explain the same situation. This article provides a unifying perspective of bonding that captures the structural diversity across this diverse family of multimetallic clusters.

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-.

Material-Specific Optimization of Gaussian Basis Sets against Plane Wave Data

Since in periodic systems, a given element may be present in different spatial arrangements displaying vastly different physical and chemical properties, an elemental basis set that is independent of physical properties of materials may lead to significant simulation inaccuracies. To avoid such a lack of material specificity within a given basis set, we present a material-specific Gaussian basis optimization scheme for solids, which simultaneously minimizes the total energy of the system and optimizes the band energies when compared to the reference plane wave calculation while taking care of the overlap matrix condition number. To assess this basis set optimization scheme, we compare the quality of the Gaussian basis sets generated for diamond, graphite, and silicon via our method against the existing basis sets. The optimization scheme of this work has also been tested on the existing Gaussian basis sets for periodic systems such as MoS₂ and NiO, yielding improved results.

Perturbation Theory Treatment of Spin-Orbit Coupling, Part I: Double Perturbation Theory Based on a Single-Reference Initial Approximation

We develop a perturbation theory for solving the many-body Dirac equation within a given relativistic effective-core potential approximation. Starting from a scalar-relativistic unrestricted Hartree-Fock (SR UHF) solution, we carry out a double perturbation expansion in terms of spin-orbit coupling (SOC) and the electron fluctuation potential. Computationally convenient energy expressions are derived through fourth order in SOC, second order in the electron fluctuation potential, and a total of third order in the coupling between the two. Illustrative calculations on the halogen series of neutral and singly positive diatomic molecules show that the perturbation expansion is well-converged by taking into account only the leading (nonvanishing) term at each order of the electron fluctuation potential. Our perturbation theory approach provides a computationally attractive alternative to a two-component self-consistent field treatment of SOC. In addition, it includes coupling with the fluctuation potential through third order and can be extended (in principle) to multireference calculations, when necessary for both closed- and open-shell cases, using quasi-degenerate perturbation theory.

Computational modelling of Pd-catalysed alkoxycarbonylation of alkenes and alkynes

This perspective highlights the computational modelling of alkene and alkyne alkoxycarbonylation at palladium catalysts. We cover studies on Pd-catalysed alkoxycarbonylation of alkenes with bidentate diphosphine ligands, which reveal a hydride pathway is operating with an intermolecular alcoholysis step, where explicit solvation is mandatory to estimate the overall barriers correctly and model alcoholysis/copolymerisation selectivities. Subsequently, we discuss Pd-catalysed alkyne alkoxycarbonylation with P,N-chelating ligands, where an in situ base mechanism is operating involving ketene-type intermediates. We also discuss catalyst poisoning due to allene and designing a potential new catalyst tolerant towards allene poisoning.

Perturbation Theory Treatment of Spin-Orbit Coupling II: A Coupled Perturbed Kohn-Sham Method

A noncanonical coupled perturbed Kohn-Sham density functional theory (KS-DFT)/Hartree-Fock (HF) treatment of spin-orbit coupling (SOC) is provided. We take the scalar-relativistic KS-DFT/HF solution, obtained with a relativistic effective core potential, as the zeroth-order approximation. Explicit expressions are given for the total energy through the 4th order, which satisfy the 2n + 1 rule. Second-order expressions are provided for orbital energies and density variables of spin-current DFT. Test calculations are carried out on the halogen homonuclear diatomic and hydride molecules, including 6p and 7p elements, as well as open-shell negative ions. The computed properties through second or third order match well with those from reference two-component self-consistent field calculations for total and orbital energies as well as spin-current densities. In only one case (At₂^-) did a significant deviation occur for the remaining density variables. Our coupled perturbation theory approach provides an efficient way of adding the effect of SOC to a scalar-relativistic single-reference KS-DFT/HF treatment, in particular because it does not require diagonalization in the two-component spinor basis, leading to saving factors on the number of required floating-point operations that may exceed one order of magnitude.