Spin-orbit effects on the (119)Sn magnetic-shielding tensor in solids: a ZORA/DFT investigation

Understanding Cu(i) local environments in MOFs via63/65Cu NMR spectroscopy

The field of metal-organic frameworks (MOFs) includes a vast number of hybrid organic and inorganic porous materials with wide-ranging applications. In particular, the Cu(i) ion exhibits rich coordination chemistry in MOFs and can exist in two-, three-, and four-coordinate environments, which gives rise to many structural motifs and potential applications. Direct characterization of the structurally and chemically important Cu(i) local environments is essential for understanding the sources of specific MOF properties. For the first time, 63/65Cu solid-state NMR has been used to investigate a variety of Cu(i) sites and local coordination geometries in Cu MOFs. This approach is a sensitive probe of the local Cu environment, particularly when combined with density functional theory calculations. A wide range of structurally-dependent 63/65Cu NMR parameters have been observed, including 65Cu quadrupolar coupling constants ranging from 18.8 to 74.8 MHz. Using the data from this and prior studies, a correlation between Cu quadrupolar coupling constants, Cu coordination number, and local Cu coordination geometry has been established. Links between DFT-calculated and experimental Cu NMR parameters are also presented. Several case studies illustrate the feasibility of 63/65Cu NMR for investigating and resolving inequivalent Cu sites, monitoring MOF phase changes, interrogating the Cu oxidation number, and characterizing the product of a MOF chemical reaction involving Cu(ii) reduction to Cu(i). A convenient avenue to acquire accurate 65Cu NMR spectra and NMR parameters from Cu(i) MOFs at a widely accessible magnetic field of 9.4 T is described, with a demonstrated practical application for tracking Cu(i) coordination evolution during MOF anion exchange. This work showcases the power of 63/65Cu solid-state NMR spectroscopy and DFT calculations for molecular-level characterization of Cu(i) centers in MOFs, along with the potential of this protocol for investigating a wide variety of MOF structural changes and processes important for practical applications. This approach has broad applications for examining Cu(i) centers in other weight-dilute systems.

Structure and bonding in rhodium coordination compounds: a 103Rh solid-state NMR and relativistic DFT study

This study demonstrates the application of 103Rh solid-state NMR (SSNMR) spectroscopy to inorganic and organometallic coordination compounds, in combination with relativistic density functional theory (DFT) calculations of 103Rh chemical shift tensors and their analysis with natural bond orbital (NBO) and natural localized molecular orbital (NLMO) protocols, to develop correlations between 103Rh chemical shift tensors, molecular structure, and Rh-ligand bonding. 103Rh is one of the least receptive NMR nuclides, and consequently, there are very few reports in the literature. We introduce robust 103Rh SSNMR protocols for stationary samples, which use the broadband adiabatic inversion-cross polarization (BRAIN-CP) pulse sequence and wideband uniform-rate smooth-truncation (WURST) pulses for excitation, refocusing, and polarization transfer, and demonstrate the acquisition of 103Rh SSNMR spectra of unprecedented signal-to-noise and uniformity. The 103Rh chemical shift tensors determined from these spectra are complemented by NBO/NLMO analyses of contributions of individual orbitals to the 103Rh magnetic shielding tensors to understand their relationship to structure and bonding. Finally, we discuss the potential for these experimental and theoretical protocols for investigating a wide range of materials containing the platinum group elements.

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.

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.

Reversible Activation and Transfer of White Phosphorus by Silyl-Stannylene

Use of a silyl supported stannylene (Mes TerSn(Sit Bu3 ) [Mes Ter=2,6-(2,4,6-Me3 C6 H2 )2 C6 H3 ] enables activation of white phosphorus under mild conditions, which is reversible under UV light. The reaction of a silylene chloride with the activated P4 complex results in facile P-atom transfer. The computational analysis rationalizes the electronic features and high reactivity of the heteroleptic silyl-substituted stannylene in contrast to the previously reported bis(aryl)stannylene.

A GIPAW versus GIAO-ZORA-SO study of 13C and 15N CPMAS NMR chemical shifts of aromatic and heterocyclic bromo derivatives

Theoretical simulation of NMR parameters in compounds bearing heavy atoms generally requires the application of relativistic corrections. We report herein the theoretical characterization of ¹³C and ¹⁵N CPMAS NMR of known bromo-derivative crystals by using both the GIPAW and the combined GIAO-ZORA-SO approximation methods. Several statistical analyses were performed to compare both approaches, with non-relativistic GIPAW method being more useful to predict the ¹³C and ¹⁵N chemical shifts. The problem of applying GIPAW to crystal structures showing static or dynamic crystalline disorder of the special class resulting in half-protons will be discussed in detail.

A DFT/ZORA Study of Cadmium Magnetic Shielding Tensors: Analysis of Relativistic Effects and Electronic-State Approximations

Theoretical considerations are discussed for the accurate prediction of cadmium magnetic shielding tensors using relativistic density functional theory (DFT). Comparison is made between calculations that model the extended lattice of the cadmium-containing solids using periodic boundary conditions and pseudopotentials with calculations that use clusters of atoms. The all-electron cluster-based calculations afford an opportunity to examine the importance of (i) relativistic effects on cadmium magnetic shielding tensors, as introduced through the ZORA Hamiltonian at either the scalar (SC) or spin-orbit (SO) levels and (ii) variation in the class of the DFT approximation. Twenty-three combinations of pseudopotentials or all-electron methods, DFT functionals, and relativistic treatments are assessed for the prediction of the principal components of the magnetic shielding tensors of 30 cadmium sites. We find that the inclusion of SO coupling can increase the cadmium magnetic shielding by as much as ca. 1100 ppm for a certain principal values; these effects are most pronounced for cadmium sites featuring bonds to other heavy atoms such as cadmium, iodine, or selenium. The best agreement with experimental values is found at the ZORA SO level in combination with a hybrid DFT method featuring a large admixture of Hartree-Fock exchange such as BH&HLYP. Finally, a theoretical examination is presented of the magnetic shielding tensor of the Cd(I) site in Cd₂(AlCl₄)₂.

Spin-orbit effects on the 125Te magnetic-shielding tensor: A cluster-based ZORA/DFT investigation

Cluster-based calculations of ¹²⁵Te magnetic-shielding tensors demonstrate that inclusion of spin-orbit effects is necessary to obtain the best agreement of theoretical predictions with experiment. The spin-orbit contribution to shielding depends on the oxidation state and stereochemistry of the ¹²⁵Te site. Comparison of the performance of various density functionals indicates that GGA functionals behave similarly to each other in predicting NMR magnetic shielding. The use of hybrid functionals improves the predictive ability on average for a large set of ¹²⁵Te-containing materials. The amount of Hartree-Fock exchange affects the predicted parameters. Inclusion of larger Hartree-Fock exchange contributions in hybrid functionals results in larger slopes of the correlation between calculated magnetic-shielding and experimental chemical-shift principal components, by 10-15% from the ideal value.

Tin chemical shift anisotropy in tin dioxide: On ambiguity of CSA asymmetry derived from MAS spectra

Two different axial symmetries of the ¹¹⁹Sn chemical shift anisotropy (CSA) in tin dioxide with the asymmetry parameter (η) of 0 and 0.27 were reported previously based on the analysis of MAS NMR spectra. By analyzing the static powder pattern, we show that the ¹¹⁹Sn CSA is axially symmetric. A nearly axial symmetry and the principal axis system of the ¹¹⁹Sn chemical shift tensor in SnO₂ were deduced from periodic scalar-relativistic density functional theory (DFT) calculations of NMR parameters. The implications of fast small-angle motions on CSA parameters were also considered, which could potentially lead to a CSA symmetry in disagreement with a crystal symmetry. Our analysis of experimental spectra using spectral simulations and iterative fittings showed that MAS spectra recorded at relatively high frequencies do not show sufficiently distinct features in order to distinguish CSAs with η ≈ 0 and η ≈ 0.4. The example of SnO₂ shows that both the MAS lineshape and spinning sideband analyses may overestimate the η value by as much as ∼0.3 and ∼0.4, respectively. The results confirm that a static powder pattern must be analysed in order to improve the accuracy of the CSA asymmetry measurements. The measurements on SnO₂ nanoparticles showed that the asymmetry parameter of the ¹¹⁹Sn CSA increases for nm-sized particles with a larger surface area compared to μm-sized particles. The increase of the η value for tin atoms near the surface in SnO₂ was also confirmed by DFT calculations.