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.

Benchmark Study on the Calculation of 207Pb NMR Chemical Shifts

A benchmark set for the computation of 207Pb nuclear magnetic resonance (NMR) chemical shifts is presented. The PbS50 set includes conformer ensembles of 50 lead-containing molecular compounds and their experimentally measured 207Pb NMR chemical shifts. Various bonding motifs at the Pb center with up to seven bonding partners are included. Six different solvents were used in the measurements. The respective shifts lie in the range between +10745 and -5030 ppm. Several calculation settings are assessed by evaluating computed 207Pb NMR shifts for the use with different density functional approximations (DFAs), relativistic approaches, treatment of the conformational space, and levels for geometry optimization. Relativistic effects were included explicitly with the zeroth order regular approximation (ZORA), for which only the spin-orbit variant was able to yield reliable results. In total, seven GGAs and three hybrid DFAs were tested. Hybrid DFAs significantly outperform GGAs. The most accurate DFAs are mPW1PW with a mean absolute deviation (MAD) of 429 ppm and PBE0 with an MAD of 446 ppm. Conformational influences are small as most compounds are rigid, but more flexible structures still benefit from Boltzmann averaging. Including explicit relativistic treatments such as SO-ZORA in the geometry optimization does not show any significant improvement over the use of effective core potentials (ECPs).

Confined Lewis Pairs: Investigation of the X- →Si20 Interaction in Halogen-Encapsulating Silafulleranes

A joint theoretical and experimental study on 32 endohedral silafullerane derivatives [X@Si20 Y20 ]- (X=F-I; Y=F-I, H, Me, Et) andT h ${T_h }$ -[Cl@Si20 H12 Y8 ]- (Y=F-I) is presented. First, we evaluated the structure-determining template effect of Cl- in a systematic series of concave silapolyquinane model systems. Second, we investigated the X- →Si20 interaction energy (E int ${E_{{\rm{int}}} }$ ) as a function of X- and Y and found the largestE int ${E_{{\rm{int}}} }$ values for electron-withdrawing exohedral substituents Y. Given that X- ions can be considered as Lewis bases and empty Si20 Y20 clusters as Lewis acids, we classify our inseparable host-guest complexes [X@Si20 Y20 ]- as "confined Lewis pairs". Third, 35 Cl NMR spectroscopy proved to be highly diagnostic for an experimental assessment of the Cl- →Si20 interaction as the paramagnetic shielding and, in turn, δ ${\delta }$ (35 Cl) of the endohedral Cl- ion correlate inversely withE int ${E_{{\rm{int}}} }$ . Finally, we disclose the synthesis of [PPN][Cl@Si20 Y20 ] (Y=Me, Et, Br) and provide a thorough characterization of these new silafulleranes.

Brominated [20]silafulleranes: pushing the limits of steric loading

Starting from the perhydrogenated silafullerane [nBu₄N][Cl@Si₂₀(SiH₃)₁₂H₈], treatment with BBr₃ leads to partially and exhaustively brominated clusters, [nBu₄N][Cl@Si₂₀(SiBr₂H)₁₂Br₈] (120 eq. BBr₃, room temperature, 30 min) and [nBu₄N][Cl@Si₂₀(SiBr₃)₁₂Br₈] (300 eq. BBr₃, 130 °C, 3 d). Perbromination is accompanied by a massively increased steric strain on the cluster surface, which explains why our approach achieves regioselective derivatization of the Si₃₂ framework when mild conditions are maintained. Partial Br/H exchange on [nBu₄N][Cl@Si₂₀(SiBr₂H)₁₂Br₈] (30 eq. iBu₂AlH, room temperature, 16 h) affords [nBu₄N][Cl@Si₂₀(SiH₃)₁₂Br₈].

Regioselective Derivatization of Silylated [20]Silafulleranes

Silafulleranes with endohedral Cl^- ions are a unique, scarcely explored class of structurally well-defined silicon clusters and host-guest complexes. Herein, we report regioselective derivatization reactions on the siladodecahedrane [nBu₄N][Cl@Si₂₀(SiCl₃)₁₂Cl₈] ([nBu₄N][1]), which has its cluster surface decorated with 12 SiCl₃ and 8 Cl substituents in perfect T_h symmetry. The room-temperature reaction of [nBu₄N][1] with excess iBu₂AlH in ortho-difluorobenzene (oDFB) furnishes perhydrogenated [nBu₄N][Cl@Si₂₀(SiH₃)₁₂H₈] ([nBu₄N][2]) in 50% yield; the non-pyrophoric [2]^- is the largest structurally authenticated (by X-ray diffraction) hydridosilane known to date. A simple switch from pure oDFB to an oDFB/Et₂O solvent mixture suppresses core hydrogenation and results in the formation of [nBu₄N][Cl@Si₂₀(SiH₃)₁₂Cl₈] ([nBu₄N][3]). In addition to the exhaustive Cl/H exchange at all 44 Si-Cl bonds of [1]^- and the regioselective 36-fold silyl group hydrogenation, we achieved the simultaneous introduction of Me substituents at all 8 SiCl vertices along with the conversion of all 12 SiCl₃ to SiH₃ groups by treating [nBu₄N][1] with Me₂AlH/Me₃Al in oDFB ([nBu₄N][Cl@Si₂₀(SiH₃)₁₂Me₈], [nBu₄N][4]; 73%). Quantum-chemical free-energy calculations find an S_N2-Si-type hydrogenation of the exohedral SiCl₃ moieties in [1]^- (trigonal-bipyramidal intermediate) slightly preferred over metathesis-like S_Ni-Si substitutions (four-membered transition state). Cage hydrogenation likely occurs via S_Ni-Si processes. The experimentally demonstrated influence of an Et₂O co-solvent, which drastically increases the respective reaction barriers, is attributed to the increased stability of the resulting iBu₂AlH-OEt₂ adduct and its higher steric bulk compared to free iBu₂AlH.

Optimization of the r2SCAN-3c Composite Electronic-Structure Method for Use with Slater-Type Orbital Basis Sets

The "Swiss army knife" composite density functional electronic-structure method r²SCAN-3c (J. Chem. Phys. 2021, 154, 064103) is extended and optimized for the use with Slater-type orbital basis sets. The meta generalized-gradient approximation (meta-GGA) functional r²SCAN by Furness et al. is combined with a tailor-made polarized triple-ζ Slater-type atomic orbital (STO) basis set (mTZ2P), the semiclassical London dispersion correction (D4), and a geometrical counterpoise (gCP) correction. Relativistic effects are treated explicitly with the scalar-relativistic zeroth-order regular approximation (SR-ZORA). The performance of the new implementation is assessed on eight geometry and 74 energy benchmark sets, including the extensive GMTKN55 database as well as recent sets such as ROST61 and IONPI19. In geometry optimizations, the STO-based r²SCAN-3c is either on par with or more accurate than the hybrid density functional approximation M06-2X-D3(0)/TZP. In energy calculations, the overall accuracy is similar to the original implementation of r²SCAN-3c with Gaussian-type atomic orbitals (GTO), but basic properties, intermolecular noncovalent interactions, and barrier heights are better described with the STO approach, resulting in a lower weighted mean absolute deviation (WTMAD-2(STO) = 7.15 vs 7.50 kcal mol^-1 with the original method) for the GMTKN55 database. The STO-optimized r²SCAN-3c outperforms many conventional hybrid/QZ approaches in most common applications at a fraction of their cost. The reliable, robust, and accurate r²SCAN-3c implementation with STOs is a promising alternative to the original implementation with GTOs and can be generally used for a broad field of quantum chemical problems.

Benchmark Study on the Calculation of 119Sn NMR Chemical Shifts

A new benchmark set termed SnS51 for assessing quantum chemical methods for the computation of ¹¹⁹Sn NMR chemical shifts is presented. It covers 51 unique ¹¹⁹Sn NMR chemical shifts for a selection of 50 tin compounds with diverse bonding motifs and ligands. The experimental reference data are in the spectral range of ±2500 ppm measured in seven different solvents. Fifteen common density functional approximations, two scalar- and one spin-orbit relativistic approach are assessed based on conformer ensembles generated using the CREST/CENSO scheme and state-of-the-art semiempirical (GFN2-xTB), force field (GFN-FF), and composite DFT methods (r²SCAN-3c). Based on the results of this study, the spin-orbit relativistic method combinations of SO-ZORA with PBE0 or revPBE functionals are generally recommended. Both yield mean absolute deviations from experimental data below 100 ppm and excellent linear regression determination coefficients of ≤0.99. If spin-orbit calculations are not affordable, the use of SR-ZORA with B3LYP or X2C with ωB97X or M06 may be considered to obtain qualitative predictions if no severe spin-orbit effects, for example, due to heavy nuclei containing ligands, are expected. An empirical linear scaling correction is demonstrated to be applicable for further improvement, and respective empirical parameters are given. Conformational effects on chemical shifts are studied in detail but are mostly found to be small. However, in specific cases when the ligand sphere differs substantially between conformers, chemical shifts can change by up to several hundred ppm.

Comprehensive Benchmark Study on the Calculation of 29Si NMR Chemical Shifts

A comprehensive and diverse benchmark set for the calculation of ²⁹Si NMR chemical shifts is presented. The SiS146 set includes 100 silicon containing compounds with 146 experimentally determined reference ²⁹Si NMR chemical shifts measured in nine different solvents in a range from -400 to +828 ppm. Silicon atoms bound to main group elements as well as transition metals with coordination numbers of 2-6 in various bonding patterns including multiple bonds and coordinative and aromatic bonding are represented. The performance of various common and specialized density functional approximations including (meta-)GGA, hybrid, and double-hybrid functionals in combination with different AO basis sets and for differently optimized geometries is evaluated. The role of scalar-relativistic effects is further investigated by inclusion of the zeroth order regular approximation (ZORA) method into the calculations. GGA density functional approximations (DFAs) are found to outperform hybrid DFAs with B97-D3 performing best with an MAD of 7.2 ppm for the subset including only light atoms (Z < 18), while TPSSh is the best tested hybrid functional with an MAD of 10.3 ppm. For ²⁹Si cores in the vicinity of heavier atoms, the application of ZORA proved indispensable. Inclusion of spin-orbit effects into the ²⁹Si NMR chemical shift calculation decreases the mean absolute deviations by up to 74% compared to calculations applying effective core potentials.