Computational NMR of the iron pyrazolylborate complexes [Tp2Fe]+ and Tp2Fe including solvation and spin-crossover effects.

Transition metal complexes have important roles in many biological processes as well as applications in fields such as pharmacy, chemistry and materials science. Paramagnetic nuclear magnetic resonance (pNMR) is a valuable tool in understanding such molecules, and theoretical computations are often advantageous or even necessary in the assignment of experimental pNMR signals. We have employed density functional theory (DFT) and the domain-based local pair natural orbital coupled-cluster method with single and double excitations (DLPNO-CCSD), as well as a number of model improvements, to determine the critical hyperfine part of the chemical shifts of the iron pyrazolylborate complexes [Tp2Fe]+ and Tp2Fe using a modern version of the Kurland-McGarvey theory, which is based on parameterising the hyperfine, electronic Zeeman and zero-field splitting interactions via the parameters of the electron paramagnetic resonance Hamiltonian. In the doublet [Tp2Fe]+ system, the calculations suggest a re-assignment of the 13C signal shifts. Consideration of solvent via the conductor-like polarisable continuum model (C-PCM) versus explicit solvent molecules reveals C-PCM alone to be insufficient in capturing the most important solvation effects. Tp2Fe exhibits a spin-crossover effect between a high-spin quintet (S = 2) and a low-spin singlet (S = 0) state, and its recorded temperature dependence can only be reproduced theoretically by accounting for the thermal Boltzmann distribution of the open-shell excited state and the closed-shell ground-state occupations. In these two cases, DLPNO-CCSD is found, in calculating the hyperfine couplings, to be a viable alternative to DFT, the demonstrated shortcomings of which have been a significant issue in the development of computational pNMR.

Paramagnetic Pyrazolylborate Complexes Tp2M and Tp*2M: 1H, 13C, 11B, and 14N NMR Spectra and First-Principles Studies of Chemical Shifts.

The paramagnetic pyrazolylborates Tp2M and Tp*2M (M = Cu, Ni, Co, Fe, Mn, Cr, V) as well as [Tp2M]+ and [Tp*2M]+ (M = Fe, Cr, V) have been synthesized and their NMR spectra recorded. The 1H signal shift ranges vary from ∼30 ppm (Cu(II) and V(III)) to ∼220 ppm (Co(II)), and the 13C signal shift ranges from ∼180 ppm (Fe(III)) to ∼1150 ppm (Cr(II)). The 11B and 14N shifts are ∼360 and ∼730 ppm, respectively. Both negative and positive shifts have been observed for all nuclei. The narrow NMR signals of the Co(II), Fe(II), Fe(III), and V(III) derivatives provide resolved 13C,1H couplings. All chemical shifts have been calculated from first-principles on a modern version of Kurland-McGarvey theory which includes optimized structures, zero-field splitting, and g tensors, as well as signal shift contributions. Temperature dependence in the Fe(II) spin-crossover complex results from the equilibrium of the ground singlet and the excited quintet. We illustrate both the assignment and analysis capabilities, as well as the shortcomings of the current computational methodology.

Remarkable reversal of 13C-NMR assignment in d1, d2 compared to d8, d9 acetylacetonate complexes: analysis and explanation based on solid-state MAS NMR and computations.

13C solid-state MAS NMR spectra of a series of paramagnetic metal acetylacetonate complexes; [VO(acac)2] (d1, S = ½), [V(acac)3] (d2, S = 1), [Ni(acac)2(H2O)2] (d8, S = 1), and [Cu(acac)2] (d9, S = ½), were assigned using modern NMR shielding calculations. This provided a reliable assignment of the chemical shifts and a qualitative insight into the hyperfine couplings. Our results show a reversal of the isotropic 13C shifts, δiso(13C), for CH3 and CO between the d1 and d2versus the d8 and d9 acetylacetonate complexes. The CH3 shifts change from about -150 ppm (d1,2) to roughly 1000 ppm (d8,9), whereas the CO shifts decrease from 800 ppm to about 150 ppm for d1,2 and d8,9, respectively. This was rationalized by comparison of total spin-density plots and computed contact couplings to those corresponding to singly occupied molecular orbitals (SOMOs). This revealed the interplay between spin delocalization of the SOMOs and spin polarization of the lower-energy MOs, influenced by both the molecular symmetry and the d-electron configuration. A large positive chemical shift results from spin delocalization and spin polarization acting in the same direction, whereas their cancellation corresponds to a small shift. The SOMO(s) for the d8 and d9 complexes are σ-like, implying spin-delocalization on the CH3 and CO groups of the acac ligand, cancelled only for CO by spin polarization. In contrast, the SOMOs of the d1 and d2 systems are π-like and a large CO-shift results from spin polarization, which accounts for the reversed assignment of δiso(13C) for CH3 and CO.